Menu

Auto Tag: Evolution

September 22, 2019

A chromosome scale assembly of the model desiccation tolerant grass Oropetium thomaeum

Oropetium thomaeum is an emerging model for desiccation tolerance and genome size evolution in grasses. A high-quality draft genome of Oropetium was recently sequenced, but the lack of a chromosome scale assembly has hindered comparative analyses and downstream functional genomics. Here, we reassembled Oropetium, and anchored the genome into ten chromosomes using Hi-C based chromatin interactions. A combination of high-resolution RNAseq data and homology-based gene prediction identified thousands of new, conserved gene models that were absent from the V1 assembly. This includes thousands of new genes with high expression across a desiccation timecourse. The sorghum and Oropetium genomes have a surprising degree of chromosome-level collinearity, and several chromosome pairs have near perfect synteny. Other chromosomes are collinear in the gene rich chromosome arms but have experienced pericentric translocations. Together, these resources will be useful for the grass comparative genomic community and further establish Oropetium as a model resurrection plant.

Read more
September 22, 2019

Periodic variation of mutation rates in bacterial genomes associated with replication timing

The causes and consequences of spatiotemporal variation in mutation rates remain to be explored in nearly all organisms. Here we examine relationships between local mutation rates and replication timing in three bacterial species whose genomes have multiple chromosomes: Vibrio fischeri, Vibrio cholerae, and Burkholderia cenocepacia Following five mutation accumulation experiments with these bacteria conducted in the near absence of natural selection, the genomes of clones from each lineage were sequenced and analyzed to identify variation in mutation rates and spectra. In lineages lacking mismatch repair, base substitution mutation rates vary in a mirrored wave-like pattern on opposing replichores of the large chromosomes of V. fischeri and V. cholerae, where concurrently replicated regions experience similar base substitution mutation rates. The base substitution mutation rates on the small chromosome are less variable in both species but occur at similar rates to those in the concurrently replicated regions of the large chromosome. Neither nucleotide composition nor frequency of nucleotide motifs differed among regions experiencing high and low base substitution rates, which along with the inferred ~800-kb wave period suggests that the source of the periodicity is not sequence specific but rather a systematic process related to the cell cycle. These results support the notion that base substitution mutation rates are likely to vary systematically across many bacterial genomes, which exposes certain genes to elevated deleterious mutational load.IMPORTANCE That mutation rates vary within bacterial genomes is well known, but the detailed study of these biases has been made possible only recently with contemporary sequencing methods. We applied these methods to understand how bacterial genomes with multiple chromosomes, like those of Vibrio and Burkholderia, might experience heterogeneous mutation rates because of their unusual replication and the greater genetic diversity found on smaller chromosomes. This study captured thousands of mutations and revealed wave-like rate variation that is synchronized with replication timing and not explained by sequence context. The scale of this rate variation over hundreds of kilobases of DNA strongly suggests that a temporally regulated cellular process may generate wave-like variation in mutation risk. These findings add to our understanding of how mutation risk is distributed across bacterial and likely also eukaryotic genomes, owing to their highly conserved replication and repair machinery. Copyright © 2018 Dillon et al.

Read more
September 22, 2019

Plasmodium vivax-like genome sequences shed new insights into Plasmodium vivax biology and evolution.

Although Plasmodium vivax is responsible for the majority of malaria infections outside Africa, little is known about its evolution and pathway to humans. Its closest genetic relative, P. vivax-like, was discovered in African great apes and is hypothesized to have given rise to P. vivax in humans. To unravel the evolutionary history and adaptation of P. vivax to different host environments, we generated using long- and short-read sequence technologies 2 new P. vivax-like reference genomes and 9 additional P. vivax-like genotypes. Analyses show that the genomes of P. vivax and P. vivax-like are highly similar and colinear within the core regions. Phylogenetic analyses clearly show that P. vivax-like parasites form a genetically distinct clade from P. vivax. Concerning the relative divergence dating, we show that the evolution of P. vivax in humans did not occur at the same time as the other agents of human malaria, thus suggesting that the transfer of Plasmodium parasites to humans happened several times independently over the history of the Homo genus. We further identify several key genes that exhibit signatures of positive selection exclusively in the human P. vivax parasites. Two of these genes have been identified to also be under positive selection in the other main human malaria agent, P. falciparum, thus suggesting their key role in the evolution of the ability of these parasites to infect humans or their anthropophilic vectors. Finally, we demonstrate that some gene families important for red blood cell (RBC) invasion (a key step of the life cycle of these parasites) have undergone lineage-specific evolution in the human parasite (e.g., reticulocyte-binding proteins [RBPs]).

Read more
September 22, 2019

Identification of the DNA methyltransferases establishing the methylome of the cyanobacterium Synechocystis sp. PCC 6803.

DNA methylation in bacteria is important for defense against foreign DNA, but is also involved in DNA repair, replication, chromosome partitioning, and regulatory processes. Thus, characterization of the underlying DNA methyltransferases in genetically tractable bacteria is of paramount importance. Here, we characterized the methylome and orphan methyltransferases in the model cyanobacterium Synechocystis sp. PCC 6803. Single molecule real-time (SMRT) sequencing revealed four DNA methylation recognition sequences in addition to the previously known motif m5CGATCG, which is recognized by M.Ssp6803I. For three of the new recognition sequences, we identified the responsible methyltransferases. M.Ssp6803II, encoded by the sll0729 gene, modifies GGm4CC, M.Ssp6803III, encoded by slr1803, represents the cyanobacterial dam-like methyltransferase modifying Gm6ATC, and M.Ssp6803V, encoded by slr6095 on plasmid pSYSX, transfers methyl groups to the bipartite motif GGm6AN7TTGG/CCAm6AN7TCC. The remaining methylation recognition sequence GAm6AGGC is probably recognized by methyltransferase M.Ssp6803IV encoded by slr6050. M.Ssp6803III and M.Ssp6803IV were essential for the viability of Synechocystis, while the strains lacking M.Ssp6803I and M.Ssp6803V showed growth similar to the wild type. In contrast, growth was strongly diminished of the ?sll0729 mutant lacking M.Ssp6803II. These data provide the basis for systematic studies on the molecular mechanisms impacted by these methyltransferases.

Read more
September 22, 2019

Sequencing of pT5282-CTXM, p13190-KPC and p30860-NR, and comparative genomics analysis of IncX8 plasmids.

This study proposes a replicon-based scheme for typing IncX plasmids into nine separately clustering subgroups, including IncX1a, IncX1ß and IncX2-8. The complete nucleotide sequences of three IncX8 plasmids, namely pT5282-CTXM and p30860-NR from Enterobacter cloacae and p13190-KPC from Klebsiella pneumoniae, were determined and were compared with two other previously sequenced IncX8 plasmids (pCAV1043-58 and pCAV1741-16). These five plasmids possessed conserved IncX8 backbones with limited genetic variation with respect to gene content and organisation, and each of them carried one or three accessory modules that harboured resistance markers and metabolic gene clusters as well as transposons, insertion sequence (IS)-based transposition units and miniature inverted repeat transposable elements (MITEs), indicating that the relatively small IncX8 backbones were able to integrate various foreign genetic contents. The resistance genes blaCTX-M-3 and blaTEM-1 (ß-lactam resistance), blaKPC-2 (carbapenem resistance) and ?blaTEM-1, and tet(A) (tetracycline resistance) and mph(E) (macrolide resistance) were found in pT5282-CTXM, p13190-KPC and pCAV1741-16, respectively, whilst p30860-NR and pCAV1043-58 carried no resistance genes. The data presented here provide an insight into the diversification and evolution history of IncX8 plasmids. Copyright © 2018 Elsevier B.V. and International Society of Chemotherapy. All rights reserved.

Read more
September 22, 2019

Tracing genomic divergence of Vibrio bacteria in the Harveyi clade.

The mechanism of bacterial speciation remains a topic of tremendous interest. To understand the ecological and evolutionary mechanisms of speciation in Vibrio bacteria, we analyzed the genomic dissimilarities between three closely related species in the so-called Harveyi clade of the genus Vibrio, V. campbellii, V. jasicida, and V. hyugaensis The analysis focused on strains isolated from diverse geographic locations over a long period of time. The results of phylogenetic analyses and calculations of average nucleotide identity (ANI) supported the classification of V. jasicida and V. hyugaensis into two species. These analyses also identified two well-supported clades in V. campbellii; however, strains from both clades were classified as members of the same species. Comparative analyses of the complete genome sequences of representative strains from the three species identified higher syntenic coverage between genomes of V. jasicida and V. hyugaensis than that between the genomes from the two V. campbellii clades. The results from comparative analyses of gene content between bacteria from the three species did not support the hypothesis that gene gain and/or loss contributed to their speciation. We also did not find support for the hypothesis that ecological diversification toward associations with marine animals contributed to the speciation of V. jasicida and V. hyugaensis Overall, based on the results obtained in this study, we propose that speciation in Harveyi clade species is a result of stochastic diversification of local populations, which was influenced by multiple evolutionary processes, followed by extinction events.IMPORTANCE To investigate the mechanisms underlying speciation in the genus Vibrio, we provided a well-assembled reference of genomes and performed systematic genomic comparisons among three evolutionarily closely related species. We resolved taxonomic ambiguities and identified genomic features separating the three species. Based on the study results, we propose a hypothesis explaining how species in the Harveyi clade of Vibrio bacteria diversified. Copyright © 2018 American Society for Microbiology.

Read more
September 22, 2019

Analysis of the draft genome of the red seaweed Gracilariopsis chorda provides insights into genome size evolution in Rhodophyta.

Red algae (Rhodophyta) underwent two phases of large-scale genome reduction during their early evolution. The red seaweeds did not attain genome sizes or gene inventories typical of other multicellular eukaryotes. We generated a high-quality 92.1 Mb draft genome assembly from the red seaweed Gracilariopsis chorda, including methylation and small (s)RNA data. We analyzed these and other Archaeplastida genomes to address three questions: 1) What is the role of repeats and transposable elements (TEs) in explaining Rhodophyta genome size variation, 2) what is the history of genome duplication and gene family expansion/reduction in these taxa, and 3) is there evidence for TE suppression in red algae? We find that the number of predicted genes in red algae is relatively small (4,803-13,125 genes), particularly when compared with land plants, with no evidence of polyploidization. Genome size variation is primarily explained by TE expansion with the red seaweeds having the largest genomes. Long terminal repeat elements and DNA repeats are the major contributors to genome size growth. About 8.3% of the G. chorda genome undergoes cytosine methylation among gene bodies, promoters, and TEs, and 71.5% of TEs contain methylated-DNA with 57% of these regions associated with sRNAs. These latter results suggest a role for TE-associated sRNAs in RNA-dependent DNA methylation to facilitate silencing. We postulate that the evolution of genome size in red algae is the result of the combined action of TE spread and the concomitant emergence of its epigenetic suppression, together with other important factors such as changes in population size.

Read more
September 22, 2019

Genomes of ubiquitous marine and hypersaline Hydrogenovibrio, Thiomicrorhabdus and Thiomicrospira spp. encode a diversity of mechanisms to sustain chemolithoautotrophy in heterogeneous environments.

Chemolithoautotrophic bacteria from the genera Hydrogenovibrio, Thiomicrorhabdus and Thiomicrospira are common, sometimes dominant, isolates from sulfidic habitats including hydrothermal vents, soda and salt lakes and marine sediments. Their genome sequences confirm their membership in a deeply branching clade of the Gammaproteobacteria. Several adaptations to heterogeneous habitats are apparent. Their genomes include large numbers of genes for sensing and responding to their environment (EAL- and GGDEF-domain proteins and methyl-accepting chemotaxis proteins) despite their small sizes (2.1-3.1 Mbp). An array of sulfur-oxidizing complexes are encoded, likely to facilitate these organisms’ use of multiple forms of reduced sulfur as electron donors. Hydrogenase genes are present in some taxa, including group 1d and 2b hydrogenases in Hydrogenovibrio marinus and H. thermophilus MA2-6, acquired via horizontal gene transfer. In addition to high-affinity cbb3 cytochrome c oxidase, some also encode cytochrome bd-type quinol oxidase or ba3 -type cytochrome c oxidase, which could facilitate growth under different oxygen tensions, or maintain redox balance. Carboxysome operons are present in most, with genes downstream encoding transporters from four evolutionarily distinct families, which may act with the carboxysomes to form CO2 concentrating mechanisms. These adaptations to habitat variability likely contribute to the cosmopolitan distribution of these organisms.© 2018 Society for Applied Microbiology and John Wiley & Sons Ltd.

Read more
September 22, 2019

Emergence of an XDR and carbapenemase-producing hypervirulent Klebsiella pneumoniae strain in Taiwan.

Carbapenemase-producing Klebsiella pneumoniae causes high mortality owing to the limited therapeutic options available. Here, we investigated an emergent carbapenem-resistant K. pneumoniae strain with hypervirulence found among KPC-2-producing strains in Taiwan.KPC-producing K. pneumoniae strains were collected consecutively from clinical specimens at the Taipei Veterans General Hospital between January 2012 and December 2014. Capsular types and the presence of rmpA/rmpA2 were analysed, and PFGE and MLST performed using these strains. The strain positive for rmpA/rmpA2 was tested in an in vivo mouse lethality study to verify its virulence and subjected to WGS to delineate its genomic features.A total of 62 KPC-2-producing K. pneumoniae strains were identified; all of these belonged to ST11 and capsular genotype K47. One strain isolated from a fatal case with intra-abdominal abscess (TVGHCRE225) harboured rmpA and rmpA2 genes. This strain was resistant to tigecycline and colistin, in addition to carbapenems, and did not belong to the major cluster in PFGE. TVGHCRE225 exhibited high in vivo virulence in the mouse lethality experiment. WGS showed that TVGHCRE225 acquired a novel hybrid virulence plasmid harbouring a set of virulence genes (iroBCDN, iucABCD, rmpA and rmpA2, and iutA) compared with the classic ST11 KPC-2-producing strain.We identified an XDR ST11 KPC-2-producing K. pneumoniae strain carrying a hybrid virulent plasmid in Taiwan. Active surveillance focusing on carbapenem-resistant hypervirulent K. pneumoniae strains is necessary, as the threat to human health is imminent.

Read more
September 22, 2019

Transcriptome analysis of Neisseria gonorrhoeae during natural infection reveals differential expression of antibiotic resistance determinants between men and women.

Neisseria gonorrhoeae is a bacterial pathogen responsible for the sexually transmitted infection gonorrhea. Emergence of antimicrobial resistance (AMR) of N. gonorrhoeae worldwide has resulted in limited therapeutic choices for this infection. Men who seek treatment often have symptomatic urethritis; in contrast, gonococcal cervicitis in women is usually minimally symptomatic, but may progress to pelvic inflammatory disease. Previously, we reported the first analysis of gonococcal transcriptome expression determined in secretions from women with cervical infection. Here, we defined gonococcal global transcriptional responses in urethral specimens from men with symptomatic urethritis and compared these with transcriptional responses in specimens obtained from women with cervical infections and in vitro-grown N. gonorrhoeae isolates. This is the first comprehensive comparison of gonococcal gene expression in infected men and women. RNA sequencing analysis revealed that 9.4% of gonococcal genes showed increased expression exclusively in men and included genes involved in host immune cell interactions, while 4.3% showed increased expression exclusively in women and included phage-associated genes. Infected men and women displayed comparable antibiotic-resistant genotypes and in vitro phenotypes, but a 4-fold higher expression of the Mtr efflux pump-related genes was observed in men. These results suggest that expression of AMR genes is programed genotypically and also driven by sex-specific environments. Collectively, our results indicate that distinct N. gonorrhoeae gene expression signatures are detected during genital infection in men and women. We propose that therapeutic strategies could target sex-specific differences in expression of antibiotic resistance genes.IMPORTANCE Recent emergence of antimicrobial resistance of Neisseria gonorrhoeae worldwide has resulted in limited therapeutic choices for treatment of infections caused by this organism. We performed global transcriptomic analysis of N. gonorrhoeae in subjects with gonorrhea who attended a Nanjing, China, sexually transmitted infection (STI) clinic, where antimicrobial resistance of N. gonorrhoeae is high and increasing. We found that N. gonorrhoeae transcriptional responses to infection differed in genital specimens taken from men and women, particularly antibiotic resistance gene expression, which was increased in men. These sex-specific findings may provide a new approach to guide therapeutic interventions and preventive measures that are also sex specific while providing additional insight to address antimicrobial resistance of N. gonorrhoeae. Copyright © 2018 Nudel et al.

Read more
September 22, 2019

Chromosomally encoded mcr-5 in colistin non-susceptible Pseudomonas aeruginosa.

Whole genome sequencing (WGS) of historical Pseudomonas aeruginosa clinical isolates identified a chromosomal copy of mcr-5 within a Tn3-like transposon in P. aeruginosa MRSN 12280. The isolate was non-susceptible to colistin by broth microdilution and genome analysis revealed no mutations known to confer colistin resistance. To the best of our knowledge, this is the first report of mcr in colistin non-susceptible P. aeruginosa.

Read more
September 22, 2019

Analysis of the Gli-D2 locus identifies a genetic target for simultaneously improving the breadmaking and health-related traits of common wheat.

Gliadins are a major component of wheat seed proteins. However, the complex homoeologous Gli-2 loci (Gli-A2, -B2 and -D2) that encode the a-gliadins in commercial wheat are still poorly understood. Here we analyzed the Gli-D2 locus of Xiaoyan 81 (Xy81), a winter wheat cultivar. A total of 421.091 kb of the Gli-D2 sequence was assembled from sequencing multiple bacterial artificial clones, and 10 a-gliadin genes were annotated. Comparative genomic analysis showed that Xy81 carried only eight of the a-gliadin genes of the D genome donor Aegilops tauschii, with two of them each experiencing a tandem duplication. A mutant line lacking Gli-D2 (DLGliD2) consistently exhibited better breadmaking quality and dough functionalities than its progenitor Xy81, but without penalties in other agronomic traits. It also had an elevated lysine content in the grains. Transcriptome analysis verified the lack of Gli-D2 a-gliadin gene expression in DLGliD2. Furthermore, the transcript and protein levels of protein disulfide isomerase were both upregulated in DLGliD2 grains. Consistent with this finding, DLGliD2 had increased disulfide content in the flour. Our work sheds light on the structure and function of Gli-D2 in commercial wheat, and suggests that the removal of Gli-D2 and the gliadins specified by it is likely to be useful for simultaneously enhancing the end-use and health-related traits of common wheat. Because gliadins and homologous proteins are widely present in grass species, the strategy and information reported here may be broadly useful for improving the quality traits of diverse cereal crops.© 2018 The Authors The Plant Journal © 2018 John Wiley & Sons Ltd.

Read more
September 22, 2019

Human copy number variants are enriched in regions of low mappability.

Copy number variants (CNVs) are known to affect a large portion of the human genome and have been implicated in many diseases. Although whole-genome sequencing (WGS) can help identify CNVs, most analytical methods suffer from limited sensitivity and specificity, especially in regions of low mappability. To address this, we use PopSV, a CNV caller that relies on multiple samples to control for technical variation. We demonstrate that our calls are stable across different types of repeat-rich regions and validate the accuracy of our predictions using orthogonal approaches. Applying PopSV to 640 human genomes, we find that low-mappability regions are approximately 5 times more likely to harbor germline CNVs, in stark contrast to the nearly uniform distribution observed for somatic CNVs in 95 cancer genomes. In addition to known enrichments in segmental duplication and near centromeres and telomeres, we also report that CNVs are enriched in specific types of satellite and in some of the most recent families of transposable elements. Finally, using this comprehensive approach, we identify 3455 regions with recurrent CNVs that were missing from existing catalogs. In particular, we identify 347 genes with a novel exonic CNV in low-mappability regions, including 29 genes previously associated with disease.

Read more
September 22, 2019

Genomic signatures of mitonuclear coevolution across populations of Tigriopus californicus.

The copepod Tigriopus californicus shows extensive population divergence and is becoming a model for understanding allopatric differentiation and the early stages of speciation. Here, we report a high-quality reference genome for one population (~190?megabases across 12 scaffolds, and ~15,500 protein-coding genes). Comparison with other arthropods reveals 2,526 genes presumed to be specific to T. californicus, with an apparent proliferation of genes involved in ion transport and receptor activity. Beyond the reference population, we report re-sequenced genomes of seven additional populations, spanning the continuum of reproductive isolation. Populations show extreme mitochondrial DNA divergence, with higher levels of amino acid differentiation than observed in other taxa. Across the nuclear genome, we find elevated protein evolutionary rates and positive selection in genes predicted to interact with mitochondrial DNA and the proteins and RNA it encodes in multiple pathways. Together, these results support the hypothesis that rapid mitochondrial evolution drives compensatory nuclear evolution within isolated populations, thereby providing a potentially important mechanism for causing intrinsic reproductive isolation.

Read more
September 22, 2019

A synthetic-diploid benchmark for accurate variant-calling evaluation.

Existing benchmark datasets for use in evaluating variant-calling accuracy are constructed from a consensus of known short-variant callers, and they are thus biased toward easy regions that are accessible by these algorithms. We derived a new benchmark dataset from the de novo PacBio assemblies of two fully homozygous human cell lines, which provides a relatively more accurate and less biased estimate of small-variant-calling error rates in a realistic context.

Read more

Subscribe for blog updates:

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.