Menu
September 22, 2019

Structural variants exhibit allelic heterogeneity and shape variation in complex traits

Despite extensive effort to reveal the genetic basis of complex phenotypic variation, studies typically explain only a fraction of trait heritability. It has been hypothesized that individually rare hidden structural variants (SVs) could account for a significant fraction of variation in complex traits. To investigate this hypothesis, we assembled 14 Drosophila melanogaster genomes and systematically identified more than 20,000 euchromatic SVs, of which ~40% are invisible to high specificity short read genotyping approaches. SVs are common in Drosophila genes, with almost one third of diploid individuals harboring an SV in genes larger than 5kb, and nearly a quarter harboring multiple SVs in genes larger than 10kb. We show that SV alleles are rarer than amino acid polymorphisms, implying that they are more strongly deleterious. A number of functionally important genes harbor previously hidden structural variants that likely affect complex phenotypes (e.g., Cyp6g1, Drsl5, Cyp28d1&2, InR, and Gss1&2). Furthermore, SVs are overrepresented in quantitative trait locus candidate genes from eight Drosophila Synthetic Population Resource (DSPR) mapping experiments. We conclude that SVs are pervasive in genomes, are frequently present as heterogeneous allelic series, and can act as rare alleles of large effect.


September 22, 2019

De novo assembly, delivery and expression of a 101 kb human gene in mouse cells

Design and large-scale synthesis of DNA has been applied to the functional study of viral and microbial genomes. New and expanded technology development is required to unlock the transformative potential of such bottom-up approaches to the study of larger, mammalian genomes. Two major challenges include assembling and delivering long DNA sequences. Here we describe a pipeline for de novo DNA assembly and delivery that enables functional evaluation of mammalian genes on the length scale of 100 kb. The DNA assembly step is supported by an integrated robotic workcell. We assemble the 101 kb human HPRT1 gene in yeast, deliver it to mouse cells, and show expression of the human protein from its full-length gene. This pipeline provides a framework for producing systematic, designer variants of any mammalian gene locus for functional evaluation in cells.


September 22, 2019

Parliament2: Fast structural variant calling using optimized combinations of callers

Here we present Parliament2: a structural variant caller which combines multiple best-in-class structural variant callers to create a highly accurate callset. This captures more events than the individual callers achieve independently. Parliament2 uses a call-overlap-genotype approach that is highly extensible to new methods and presents users the choice to run some or all of Breakdancer, Breakseq, CNVnator, Delly, Lumpy, and Manta to run. Parliament2 applies an additional parallelization framework to speed certain callers and executes these in parallel, taking advantage of the different resource requirements to complete structural variant calling much faster than running the programs individually. Parliament2 is available as a Docker container, which pre-installs all required dependencies. This allows users to run any caller with easy installation and execution. This Docker container can easily be deployed in cloud or local environments and is available as an app on DNAnexus.


September 22, 2019

Identification of the KPC plasmid pCT-KPC334: New insights on the evolutionary pathway of epidemic plasmids harboring fosA3-blaKPC-2 genes.

A novel, non-conjugative plasmid pKP1034 isolated from a fosfomycin-resistant, carbapenemase-producing Klebsiella pneumonia strain KP1034 was recently reported to carry fosA3, blaKPC-2, blaCTX-M-65, blaSHV-12 and rmtB genes, and was hypothesized to evolve from several recombination events of two closely related plasmids, pHN7A8 and pKPC-LK30 [1]. In this study, a plasmid pCT-KPC334 carrying fosA3, blaKPC-2, blaCTX-M-65, blaSHV-12, blaTEM-1, and rmtB genes was identified, providing evidence on the evolutionary pathway of plasmids harboring fosA3-blaKPC-2 genes.


September 22, 2019

4.5 years within-patient evolution of a colistin resistant KPC-producing Klebsiella pneumoniae ST258.

Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae (KPC-Kp) has emerged globally over the last decade as a major nosocomial pathogen that threatens patient care. These highly resistant bacteria are mostly associated with a single Kp clonal group, CG258, but the reasons for its host and hospital adaptation remain largely unknown.We analyzed the in vivo evolution of a colistin-resistant KPC-Kp CG258 strain that contaminated a patient following an endoscopy and was responsible for a fatal bacteremia 4.5 years later. Whole-genome sequencing was performed on 17 KPC-Kp isolates from this patient; single-nucleotide polymorphisms were analyzed and their implication in antimicrobial resistance and bacterial host adaptation investigated.The patient KPC-Kp strain diversified over 4.5 years at a rate of 7.5 substitutions per genome per year, resulting in broad phenotypic modifications. After 2 years of carriage, all isolates restored susceptibility to colistin. Higher expression of the fimbriae conferred the ability to produce more biofilm, and the isolate responsible for a bacteremia grew in human serum. The convergent mutations occurring in specific pathways, such as the respiratory chain and the cell envelope, revealed a complex long-term adaptation of KPC-Kp.Broad genomic and phenotypic diversification and the parallel selection of pathoadaptive mutations might contribute to long-term carriage and virulence of KPC-Kp CG258 strains and to the dissemination of this clone.


September 22, 2019

Development and validation of 58K SNP-array and high-density linkage map in Nile tilapia (O. niloticus).

Despite being the second most important aquaculture species in the world accounting for 7.4% of global production in 2015, tilapia aquaculture has lacked genomic tools like SNP-arrays and high-density linkage maps to improve selection accuracy and accelerate genetic progress. In this paper, we describe the development of a genotyping array containing more than 58,000 SNPs for Nile tilapia (Oreochromis niloticus). SNPs were identified from whole genome resequencing of 32 individuals from the commercial population of the Genomar strain, and were selected for the SNP-array based on polymorphic information content and physical distribution across the genome using the Orenil1.1 genome assembly as reference sequence. SNP-performance was evaluated by genotyping 4991 individuals, including 689 offspring belonging to 41 full-sib families, which revealed high-quality genotype data for 43,588 SNPs. A preliminary genetic linkage map was constructed using Lepmap2 which in turn was integrated with information from the O_niloticus_UMD1 genome assembly to produce an integrated physical and genetic linkage map comprising 40,186 SNPs distributed across 22 linkage groups (LGs). Around one-third of the LGs showed a different recombination rate between sexes, with the female being greater than the male map by a factor of 1.2 (1632.9 to 1359.6 cM, respectively), with most LGs displaying a sigmoid recombination profile. Finally, the sex-determining locus was mapped to position 40.53 cM on LG23, in the vicinity of the anti-Müllerian hormone (amh) gene. These new resources has the potential to greatly influence and improve the genetic gain when applying genomic selection and surpass the difficulties of efficient selection for invasively measured traits in Nile tilapia.


September 22, 2019

Amycomicin is a potent and specific antibiotic discovered with a targeted interaction screen.

The rapid emergence of antibiotic-resistant pathogenic bacteria has accelerated the search for new antibiotics. Many clinically used antibacterials were discovered through culturing a single microbial species under nutrient-rich conditions, but in the environment, bacteria constantly encounter poor nutrient conditions and interact with neighboring microbial species. In an effort to recapitulate this environment, we generated a nine-strain actinomycete community and used 16S rDNA sequencing to deconvolute the stochastic production of antimicrobial activity that was not observed from any of the axenic cultures. We subsequently simplified the community to just two strains and identified Amycolatopsis sp. AA4 as the producing strain and Streptomyces coelicolor M145 as an inducing strain. Bioassay-guided isolation identified amycomicin (AMY), a highly modified fatty acid containing an epoxide isonitrile warhead as a potent and specific inhibitor of Staphylococcus aureus Amycomicin targets an essential enzyme (FabH) in fatty acid biosynthesis and reduces S. aureus infection in a mouse skin-infection model. The discovery of AMY demonstrates the utility of screening complex communities against specific targets to discover small-molecule antibiotics.


September 22, 2019

Genome plasticity of agr-defective Staphylococcus aureus during clinical infection.

Therapy for bacteremia caused by Staphylococcus aureus is often ineffective, even when treatment conditions are optimal according to experimental protocols. Adapted subclones, such as those bearing mutations that attenuate agr-mediated virulence activation, are associated with persistent infection and patient mortality. To identify additional alterations in agr-defective mutants, we sequenced and assembled the complete genomes of clone pairs from colonizing and infected sites of several patients in whom S. aureus demonstrated a within-host loss of agr function. We report that events associated with agr inactivation result in agr-defective blood and nares strain pairs that are enriched in mutations compared to pairs from wild-type controls. The random distribution of mutations between colonizing and infecting strains from the same patient, and between strains from different patients, suggests that much of the genetic complexity of agr-defective strains results from prolonged infection or therapy-induced stress. However, in one of the agr-defective infecting strains, multiple genetic changes resulted in increased virulence in a murine model of bloodstream infection, bypassing the mutation of agr and raising the possibility that some changes were selected. Expression profiling correlated the elevated virulence of this agr-defective mutant to restored expression of the agr-regulated ESAT6-like type VII secretion system, a known virulence factor. Thus, additional mutations outside the agr locus can contribute to diversification and adaptation during infection by S. aureus agr mutants associated with poor patient outcomes. Copyright © 2018 Altman et al.


September 22, 2019

The opium poppy genome and morphinan production.

Morphinan-based painkillers are derived from opium poppy (Papaver somniferum L.). We report a draft of the opium poppy genome, with 2.72 gigabases assembled into 11 chromosomes with contig N50 and scaffold N50 of 1.77 and 204 megabases, respectively. Synteny analysis suggests a whole-genome duplication at ~7.8 million years ago and ancient segmental or whole-genome duplication(s) that occurred before the Papaveraceae-Ranunculaceae divergence 110 million years ago. Syntenic blocks representative of phthalideisoquinoline and morphinan components of a benzylisoquinoline alkaloid cluster of 15 genes provide insight into how this cluster evolved. Paralog analysis identified P450 and oxidoreductase genes that combined to form the STORR gene fusion essential for morphinan biosynthesis in opium poppy. Thus, gene duplication, rearrangement, and fusion events have led to evolution of specialized metabolic products in opium poppy. Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.


September 22, 2019

Comparative genomics of Czech vaccine strains of Bordetella pertussis.

Bordetella pertussis is a strictly human pathogen causing the respiratory infectious disease called whooping cough or pertussis. B. pertussis adaptation to acellular pertussis vaccine pressure has been repeatedly highlighted, but recent data indicate that adaptation of circulating strains started already in the era of the whole cell pertussis vaccine (wP) use. We sequenced the genomes of five B. pertussis wP vaccine strains isolated in the former Czechoslovakia in the pre-wP (1954-1957) and early wP (1958-1965) eras, when only limited population travel into and out of the country was possible. Four isolates exhibit a similar genome organization and form a distinct phylogenetic cluster with a geographic signature. The fifth strain is rather distinct, both in genome organization and SNP-based phylogeny. Surprisingly, despite isolation of this strain before 1966, its closest sequenced relative appears to be a recent isolate from the US. On the genome content level, the five vaccine strains contained both new and already described regions of difference. One of the new regions contains duplicated genes potentially associated with transport across the membrane. The prevalence of this region in recent isolates indicates that its spread might be associated with selective advantage leading to increased strain fitness.


September 22, 2019

Convergent evolution of complex genomic rearrangements in two fungal meiotic drive elements.

Meiotic drive is widespread in nature. The conflict it generates is expected to be an important motor for evolutionary change and innovation. In this study, we investigated the genomic consequences of two large multi-gene meiotic drive elements, Sk-2 and Sk-3, found in the filamentous ascomycete Neurospora intermedia. Using long-read sequencing, we generated the first complete and well-annotated genome assemblies of large, highly diverged, non-recombining regions associated with meiotic drive elements. Phylogenetic analysis shows that, even though Sk-2 and Sk-3 are located in the same chromosomal region, they do not form sister clades, suggesting independent origins or at least a long evolutionary separation. We conclude that they have in a convergent manner accumulated similar patterns of tandem inversions and dense repeat clusters, presumably in response to similar needs to create linkage between genes causing drive and resistance.


September 22, 2019

Genomic structural variations within five continental populations of Drosophila melanogaster.

Chromosomal structural variations (SV) including insertions, deletions, inversions, and translocations occur within the genome and can have a significant effect on organismal phenotype. Some of these effects are caused by structural variations containing genes. Large structural variations represent a significant amount of the genetic diversity within a population. We used a global sampling of Drosophila melanogaster (Ithaca, Zimbabwe, Beijing, Tasmania, and Netherlands) to represent diverse populations within the species. We used long-read sequencing and optical mapping technologies to identify SVs in these genomes. Among the five lines examined, we found an average of 2,928 structural variants within these genomes. These structural variations varied greatly in size and location, included many exonic regions, and could impact adaptation and genomic evolution. Copyright © 2018 Long et al.


September 22, 2019

Fast and inexpensive protocols for consistent extraction of high quality DNA and RNA from challenging plant and fungal samples for high-throughput SNP genotyping and sequencing applications.

Modern genotyping techniques, such as SNP analysis and genotyping by sequencing (GBS), are hampered by poor DNA quality and purity, particularly in challenging plant species, rich in secondary metabolites. We therefore investigated the utility of a pre-wash step using a buffered sorbitol solution, prior to DNA extraction using a high salt CTAB extraction protocol, in a high throughput or miniprep setting. This pre-wash appears to remove interfering metabolites, such as polyphenols and polysaccharides, from tissue macerates. We also investigated the adaptability of the sorbitol pre-wash for RNA extraction using a lithium chloride-based protocol. The method was successfully applied to a variety of tissues, including leaf, cambium and fruit of diverse plant species including annual crops, forest and fruit trees, herbarium leaf material and lyophilized fungal mycelium. We consistently obtained good yields of high purity DNA or RNA in all species tested. The protocol has been validated for thousands of DNA samples by generating high data quality in dense SNP arrays. DNA extracted from Eucalyptus spp. leaf and cambium as well as mycelium from Trichoderma spp. was readily digested with restriction enzymes and performed consistently in AFLP assays. Scaled-up DNA extractions were also suitable for long read sequencing. Successful RNA quality control and good RNA-Seq data for Eucalyptus and cashew confirms the effectiveness of the sorbitol buffer pre-wash for high quality RNA extraction.


September 22, 2019

Discovery of mcr-1-mediated colistin resistance in a highly virulent Escherichia coli lineage.

Resistance to last-line polymyxins mediated by the plasmid-borne mobile colistin resistance gene (mcr-1) represents a new threat to global human health. Here we present the complete genome sequence of an mcr-1-positive multidrug-resistant Escherichia coli strain (MS8345). We show that MS8345 belongs to serotype O2:K1:H4, has a large 241,164-bp IncHI2 plasmid that carries 15 other antibiotic resistance genes (including the extended-spectrum ß-lactamase blaCTX-M-1) and 3 putative multidrug efflux systems, and contains 14 chromosomally encoded antibiotic resistance genes. MS8345 also carries a large ColV-like virulence plasmid that has been associated with E. coli bacteremia. Whole-genome phylogeny revealed that MS8345 clusters within a discrete clade in the sequence type 95 (ST95) lineage, and MS8345 is very closely related to the highly virulent O45:K1:H4 clone associated with neonatal meningitis. Overall, the acquisition of a plasmid carrying resistance to colistin and multiple other antibiotics in this virulent E. coli lineage is concerning and might herald an era where the empirical treatment of ST95 infections becomes increasingly more difficult.IMPORTANCEEscherichia coli ST95 is a globally disseminated clone frequently associated with bloodstream infections and neonatal meningitis. However, the ST95 lineage is defined by low levels of drug resistance amongst clinical isolates, which normally provides for uncomplicated treatment options. Here, we provide the first detailed genomic analysis of an E. coli ST95 isolate that has both high virulence potential and resistance to multiple antibiotics. Using the genome, we predicted its virulence and antibiotic resistance mechanisms, which include resistance to last-line antibiotics mediated by the plasmid-borne mcr-1 gene. Finding an ST95 isolate resistant to nearly all antibiotics that also has a high virulence potential is of major clinical importance and underscores the need to monitor new and emerging trends in antibiotic resistance development in this important global lineage. Copyright © 2018 Forde et al.


September 22, 2019

Comparative analysis of blaKPC-2- and rmtB-carrying IncFII-family pKPC-LK30/pHN7A8 hybrid plasmids from Klebsiella pneumoniae CG258 strains disseminated among multiple Chinese hospitals.

We recently reported the complete sequence of a blaKPC-2- and rmtB-carrying IncFII-family plasmid p675920-1 with the pKPC-LK30/pHN7A8 hybrid structure. Comparative genomics of additional sequenced plasmids with similar hybrid structures and their prevalence in blaKPC-carrying Klebsiella pneumoniae strains from China were investigated in this follow-up study.A total of 51 blaKPC-carrying K. pneumoniae strains were isolated from 2012 to 2016 from five Chinese hospitals and genotyped by multilocus sequence typing. The blaKPC-carrying plasmids from four representative strains were sequenced and compared with p675920-1 and pCT-KPC. Plasmid transfer, carbapenemase activity determination, and bacterial antimicrobial susceptibility test were performed to characterize resistance phenotypes mediated by these plasmids. The prevalence of pCT-KPC-like plasmids in these blaKPC-carrying K. pneumoniae strains was screened by PCR.The six KPC-encoding plasmids p1068-KPC, p20049-KPC, p12139-KPC and p64917-KPC (sequenced in this study) and p675920-1 and pCT-KPC slightly differed from one another due to deletion and acquisition of various backbone and accessory regions. Two major accessory resistance regions, which included the blaKPC-2 region harboring blaKPC-2 (carbapenem resistance) and blaSHV-12 (ß-lactam resistance), and the MDR region carrying rmtB (aminoglycoside resistance), fosA3 (fosfomycin resistance), blaTEM-1B (ß-lactam resistance) and blaCTX-M-65 (ß-lactam resistance), were found in each of these six plasmids and exhibited several parallel evolution routes. The pCT-KPC-like plasmids were present in all the 51 K. pneumoniae isolates, all of which belonged to CG258.There was clonal dissemination of K. pneumoniae CG258 strains, harboring blaKPC-2- and rmtB-carrying IncFII-family pKPC-LK30/pHN7A8 hybrid plasmids, among multiple Chinese hospitals.


Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.