Rich Genome Research Articles

Pandoravirus Celtis Illustrates the Microevolution Processes at Work in the Giant Pandoraviridae Genomes.

With genomes of up to 2.7 Mb propagated in μm-long oblong particles and initially predicted to encode more than 2000 proteins, members of the Pandoraviridae family display the most extreme features of the known viral world. The mere existence of such giant viruses raises fundamental questions about their origin and the processes governing their evolution. A previous analysis of six newly available isolates, independently confirmed by a study including three others, established that the Pandoraviridae pan-genome is open, meaning that each new strain exhibits protein-coding genes not previously identified in other family members. With an average increment of about 60 proteins, the gene repertoire shows no sign of reaching a limit and remains largely coding for proteins without recognizable homologs in other viruses or cells (ORFans). To explain these results, we proposed that most new protein-coding genes were created de novo, from pre-existing non-coding regions of the G+C rich pandoravirus genomes. The comparison of the gene content of a new isolate, pandoravirus celtis, closely related (96% identical genome) to the previously described p. quercus is now used to test this hypothesis by studying genomic changes in a microevolution range. Our results confirm that the differences between these two similar gene contents mostly consist of protein-coding genes without known homologs, with statistical signatures close to that of intergenic regions. These newborn proteins are under slight negative selection, perhaps to maintain stable folds and prevent protein aggregation pending the eventual emergence of fitness-increasing functions. Our study also unraveled several insertion events mediated by a transposase of the hAT family, 3 copies of which are found in p. celtis and are presumably active. Members of the Pandoraviridae are presently the first viruses known to encode this type of transposase.

Modeling of the GC content of the substituted bases in bacterial core genomes

BackgroundThe purpose of the present study was to examine the GC content of substituted bases (sbGC) in the core genomes of 35 bacterial species. Each species, or core genome, constituted genomes from at least 10 strains. We also wanted to explore whether sbGC for each strain was associated with the corresponding species’ core genome GC content (cgGC). We present a simple mathematical model that estimates sbGC from cgGC. The model assumes only that the estimated sbGC is a function of cgGC proportional to fixed AT→GC (α) and GC → AT (β) mutation rates. Non-linear regression was used to estimate parameters α and β from the empirical data described above.ResultsWe found that sbGC for each strain showed a non-linear association with the corresponding cgGC with a bias towards higher GC content for most core genomes (66.3% of the strains), assuming as a null-hypothesis that sbGC should be approximately equal to cgGC. The most GC rich core genomes (i.e. approximately %GC > 60), on the other hand, exhibited slightly less GC-biased sbGC than expected. The best fitted regression model indicates that GC → AT mutation rates β = (1.91 ± 0.13) p < 0.001 are approximately (1.91/0.79) = 2.42 times as high, on average, as AT→GC α = (− 0.79 ± 0.25) p < 0.001 mutation rates. Whether the observed sbGC GC-bias for all but the most GC-rich prokaryotic species is due to selection, compensating for the GC → AT mutation bias, and/or selective neutral processes is currently debated. Residual standard error was found to be σ = 0.076 indicating estimated errors of sbGC to be approximately within ±15.2% GC (95% confidence interval) for the strains of all species in the study.ConclusionNot only did our mathematical model give reasonable estimates of sbGC it also provides further support to previous observations that mutation rates in prokaryotes exhibit a universal GC → AT bias that appears to be remarkably consistent between taxa.

Virtual Genome Walking across the 32\u2009Gb Ambystoma mexicanum genome; assembling gene models and intronic sequence

Large repeat rich genomes present challenges for assembly using short read technologies. The 32 Gb axolotl genome is estimated to contain ~19 Gb of repetitive DNA making an assembly from short reads alone effectively impossible. Indeed, this model species has been sequenced to 20× coverage but the reads could not be conventionally assembled. Using an alternative strategy, we have assembled subsets of these reads into scaffolds describing over 19,000 gene models. We call this method Virtual Genome Walking as it locally assembles whole genome reads based on a reference transcriptome, identifying exons and iteratively extending them into surrounding genomic sequence. These assemblies are then linked and refined to generate gene models including upstream and downstream genomic, and intronic, sequence. Our assemblies are validated by comparison with previously published axolotl bacterial artificial chromosome (BAC) sequences. Our analyses of axolotl intron length, intron-exon structure, repeat content and synteny provide novel insights into the genic structure of this model species. This resource will enable new experimental approaches in axolotl, such as ChIP-Seq and CRISPR and aid in future whole genome sequencing efforts. The assembled sequences and annotations presented here are freely available for download from https://tinyurl.com/y8gydc6n. The software pipeline is available from https://github.com/LooseLab/iterassemble.

The complete mitochondrial genome sequences of two Isospora species (Eimeriidae, Eucoccidiorida, Coccidiasina, Apicomplexa) causing coccidiosis in superb glossy starlings, Lamprotornis superbus (Aves: Sturnidae)

Complete mitochondrial genomes are reported for two Isospora species causing systemic coccidiosis in Superb Glossy Starlings (Aves: Sturnidae). The A/T rich (34.7% G/C) genomes were 6223 bp in length for Isospora greineri and 6217 bp for Isospora superbusi. Each encoded 3 protein-coding genes, (COI, COIII and CytB) plus 18 LSU and 14 SSU rDNA fragments. Arrangement of protein- and rRNA-coding regions was identical to known Eimeria sp. mt genomes; start codon usage was conventional. The mitochondrial genome structures of Isospora and Eimeria species are conserved and reflect the close phylogenetic association between these eimeriid genera of apicomplexan parasites.

Elongation factor Tu is a multifunctional and processed moonlighting protein

Many bacterial moonlighting proteins were originally described in medically, agriculturally, and commercially important members of the low G + C Firmicutes. We show Elongation factor Tu (Ef-Tu) moonlights on the surface of the human pathogens Staphylococcus aureus (SaEf-Tu) and Mycoplasma pneumoniae (MpnEf-Tu), and the porcine pathogen Mycoplasma hyopneumoniae (MhpEf-Tu). Ef-Tu is also a target of multiple processing events on the cell surface and these were characterised using an N-terminomics pipeline. Recombinant MpnEf-Tu bound strongly to a diverse range of host molecules, and when bound to plasminogen, was able to convert plasminogen to plasmin in the presence of plasminogen activators. Fragments of Ef-Tu retain binding capabilities to host proteins. Bioinformatics and structural modelling studies indicate that the accumulation of positively charged amino acids in short linear motifs (SLiMs), and protein processing promote multifunctional behaviour. Codon bias engendered by an A + T rich genome may influence how positively-charged residues accumulate in SLiMs.

Comparative Genomics of Bacillus amyloliquefaciens Strains Reveals a Core Genome with Traits for Habitat Adaptation and a Secondary Metabolites Rich Accessory Genome.

The Gram positive, non-pathogenic endospore-forming soil inhabiting prokaryote Bacillus amyloliquefaciens is a plant growth-promoting rhizobacterium. Bacillus amyloliquefaciens processes wide biocontrol abilities and numerous strains have been reported to suppress diverse bacterial, fungal and fungal-like pathogens. Knowledge about strain level biocontrol abilities is warranted to translate this knowledge into developing more efficient biocontrol agents and bio-fertilizers. Ever-expanding genome studies of B. amyloliquefaciens are showing tremendous increase in strain-specific new secondary metabolite clusters which play key roles in the suppression of pathogens and plant growth promotion. In this report, we have used genome mining of all sequenced B. amyloliquefaciens genomes to highlight species boundaries, the diverse strategies used by different strains to promote plant growth and the diversity of their secondary metabolites. Genome composition of the targeted strains suggest regions of genomic plasticity that shape the structure and function of these genomes and govern strain adaptation to different niches. Our results indicated that B. amyloliquefaciens: (i) suffer taxonomic imprecision that blurs the debate over inter-strain genome diversity and dynamics, (ii) have diverse strategies to promote plant growth and development, (iii) have an unlocked, yet to be delimited impressive arsenal of secondary metabolites and products, (iv) have large number of so-called orphan gene clusters, i.e., biosynthetic clusters for which the corresponding metabolites are yet unknown, and (v) have a dynamic pan genome with a secondary metabolite rich accessory genome.

Red Algal Phylogenomics Provides a Robust Framework for Inferring Evolution of Key Metabolic Pathways.

Red algae comprise an anciently diverged, species-rich phylum with morphologies that span unicells to large seaweeds. Here, leveraging a rich red algal genome and transcriptome dataset, we used 298 single-copy orthologous nuclear genes from 15 red algal species to erect a robust multi-gene phylogeny of Rhodophyta. This tree places red seaweeds (Bangiophyceae and Florideophyceae) at the base of the mesophilic red algae with the remaining non-seaweed mesophilic lineages forming a well-supported sister group. The early divergence of seaweeds contrasts with the evolution of multicellular land plants and brown algae that are nested among multiple, unicellular or filamentous sister lineages. Using this novel perspective on red algal evolution, we studied the evolution of the pathways for isoprenoid biosynthesis. This analysis revealed losses of the mevalonate pathway on at least three separate occasions in lineages that contain Cyanidioschyzon, Porphyridium, and Chondrus. Our results establish a framework for in-depth studies of the origin and evolution of genes and metabolic pathways in Rhodophyta.

TRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes.

High-throughput genome sequencing continues to grow the need for rapid, accurate genome annotation and tRNA genes constitute the largest family of essential, ever-present non-coding RNA genes. Newly developed tRNAscan-SE 2.0 has advanced the state-of-the-art methodology in tRNA gene detection and functional prediction, captured by rich new content of the companion Genomic tRNA Database. Previously, web-server tRNA detection was isolated from knowledge of existing tRNAs and their annotation. In this update of the tRNAscan-SE On-line resource, we tie together improvements in tRNA classification with greatly enhanced biological context via dynamically generated links between web server search results, the most relevant genes in the GtRNAdb and interactive, rich genome context provided by UCSC genome browsers. The tRNAscan-SE On-line web server can be accessed at http://trna.ucsc.edu/tRNAscan-SE/.

Staphylococcus aureus Transcriptome Architecture: From Laboratory to Infection-Mimicking Conditions.

Staphylococcus aureus is a major pathogen that colonizes about 20% of the human population. Intriguingly, this Gram-positive bacterium can survive and thrive under a wide range of different conditions, both inside and outside the human body. Here, we investigated the transcriptional adaptation of S. aureus HG001, a derivative of strain NCTC 8325, across experimental conditions ranging from optimal growth in vitro to intracellular growth in host cells. These data establish an extensive repertoire of transcription units and non-coding RNAs, a classification of 1412 promoters according to their dependence on the RNA polymerase sigma factors SigA or SigB, and allow identification of new potential targets for several known transcription factors. In particular, this study revealed a relatively low abundance of antisense RNAs in S. aureus, where they overlap only 6% of the coding genes, and only 19 antisense RNAs not co-transcribed with other genes were found. Promoter analysis and comparison with Bacillus subtilis links the small number of antisense RNAs to a less profound impact of alternative sigma factors in S. aureus. Furthermore, we revealed that Rho-dependent transcription termination suppresses pervasive antisense transcription, presumably originating from abundant spurious transcription initiation in this A+T-rich genome, which would otherwise affect expression of the overlapped genes. In summary, our study provides genome-wide information on transcriptional regulation and non-coding RNAs in S. aureus as well as new insights into the biological function of Rho and the implications of spurious transcription in bacteria.

Upstream AUGs and upstream ORFs can regulate the downstream ORF in Plasmodium falciparum.

BackgroundUpstream open reading frames (uORFs) and upstream AUGs (uAUGs) can regulate the translation of downstream ORFs. The AT rich genome of Plasmodium falciparum, due to the higher AT content of start and stop codons, has the potential to give rise to a large number of uORFs and uAUGs that may affect expression of their flanking ORFs.MethodsA bioinformatics approach was used to detect uATGs associated with different genes in the parasite. To study the effect of some of these uAUGs on the expression of the downstream ORF, promoters and 5′ leaders containing uAUGs and uORFs were cloned upstream of a luciferase reporter gene. Luciferase assays were carried out in transient transfection experiments to assess the effects of uAUGs and mutations on reporter expression.ResultsThe average number of uATGs and uORFs seen in P. falciparum coding sequences (CDS) is expectedly high compared to other less biased genomes. Certain genes, including the var gene family contain the maximum number of uATGs and uORFs in the parasite. They possess ~5 times more uORFs and ~4.5 times more uAUGs within 100 bases upstream of the start codons than other CDS of the parasite. A 60 bp upstream region containing three ORFs and five ATGs from var gene PF3D7_0400100 and a gene of unknown function (PF3D7_0517100) when cloned upstream of the luciferase start codon, driven by the hsp86 promoter, resulted in loss of luciferase activity. This was restored when all the ATGs present in the −60 bp were mutated to TTGs. Point mutations in the ATGs showed that even one AUG was sufficient to repress the luciferase gene.ConclusionsOverall, this work indicates that the P. falciparum genome has a large number of uATGs and uORFs that can repress the expression of flanking ORFs. The role of AUGs in translation initiation suggests that this repression is mediated by preventing the translation initiation complex from reaching the main AUG of the downstream ORF. How the P. falciparum ribosome is able to bypass these uAUGs and uORFs for highly expressed genes remains a question for future research.Electronic supplementary materialThe online version of this article (doi:10.1186/s12936-015-1040-5) contains supplementary material, which is available to authorized users.

Effects of Genic Base Composition on Growth Rate in G+C-rich Genomes.

The source and significance of the wide variation in the genomic base composition of bacteria have been a matter of continued debate. Although the variation was originally attributed to a strictly neutral process, i.e., species-specific differences in mutational patterns, recent genomic comparisons have shown that bacteria with G+C-rich genomes experience a mutational bias toward A+T. This difference between the mutational input to a genome and its overall base composition suggests the action of natural selection. Here, we examine if selection acts on G+C contents in Caulobacter crescentus and Pseudomonas aeruginosa, which both have very G+C-rich genomes, by testing whether the expression of gene variants that differ only in their base compositions at synonymous sites affects cellular growth rates. In C. crescentus, expression of the more A+T-rich gene variants decelerated growth, indicating that selection on genic base composition is, in part, responsible for the high G+C content of this genome. In contrast, no comparable effect was observed in P. aeruginosa, which has similarly high genome G+C contents. Selection for increased genic G+C-contents in C. crescentus acts independently of the species-specific codon usage pattern and represents an additional selective force operating in bacterial genomes.

Highlight--Stacking the Deck: Heritage and Environment Shape Nucleotide Composition in Prokaryotes

The basic text of all life is written in a simple alphabet. Adenine, thymine, guanine, and cytosine form the units that allow DNA to write the books of Earth’s elaborately complex lifeforms.Curiously though, despite doing a seemingly equivalent job at conveying information, the prevalence of certain base pairs (G-C or A-T) in the genome is mysteriously variable, across both phyla and environment. Currently, it’s unclear why this should be. That makes nucleotide content a fascinating area ripe study, says Ruth Hershberg, an assistant professor at the Technion-Israel Institute of Technology in Haifa. It’s a highly variable, easy-to-measure trait. One that, she says, “affects every nucleotide in the gene. We really don’t know why it’s so variable, or what determines it, so that makes it interesting.” It might also be important too, she says. Mutations tend to move in a G-C ! A-T direction, meaning G-C rich genomes could potentially have higher mutation rates. Nucleotide content might also affect codon usage and regulatory genes, she says. Hershberg and colleagues, including Erin Reichenberger, Gail Rosen, and Uri Hershberg (her brother) have recently begun to unpick why nucleotide content varies so widely, publishing their findings online April 9 in Genome Biology Evolution (Reichenberger et al. 2015).

Targeting and function of proteins mediating translation initiation in organelles of Plasmodium falciparum.

The malaria parasite Plasmodium falciparum has two translationally active organelles - the apicoplast and mitochondrion, which import nuclear-encoded translation factors to mediate protein synthesis. Initiation of translation is a complex step wherein initiation factors (IFs) act in a regulated manner to form an initiation complex. We identified putative organellar IFs and investigated the targeting, structure and function of IF1, IF2 and IF3 homologues encoded by the parasite nuclear genome. A single PfIF1 is targeted to the apicoplast. Apart from its critical ribosomal interactions, PfIF1 also exhibited nucleic-acid binding and melting activities and mediated transcription anti-termination. This suggests a prominent ancillary function for PfIF1 in destabilisation of DNA and RNA hairpin loops encountered during transcription and translation of the A+T rich apicoplast genome. Of the three putative IF2 homologues, only one (PfIF2a) was an organellar protein with mitochondrial localisation. We additionally identified an IF3 (PfIF3a) that localised exclusively to the mitochondrion and another protein, PfIF3b, that was apicoplast targeted. PfIF3a exhibited ribosome anti-association activity, and monosome splitting by PfIF3a was enhanced by ribosome recycling factor (PfRRF2) and PfEF-G(Mit). These results fill a gap in our understanding of organellar translation in Plasmodium, which is the site of action of several anti-malarial compounds.

The complete mitochondrial genome sequence of an Isospora sp. (Eimeriidae, Eucoccidiorida, Coccidiasina, Apicomplexa) causing systemic coccidiosis in domestic Canaries (Serinus canaria Linn.)

We report a complete mitochondrial genome sequence for an Isospora sp. causing systemic coccidiosis in canaries, Serinus canaria. The A + T rich (65.2%) genome was 6216 bp in length and possessed 3 protein-coding genes, (COI; COIII and CytB), 19 LSU and 14 SSU rDNA fragments, including 1 newly identified putative LSU fragment. Arrangement of coding regions was identical to that of available Eimeria sp. mt genomes and start codon usage for protein-coding genes was conventional. The similar mitochondrial genome sequences and structures of Isospora and Eimeria species confirm the close relationship between these eimeriid genera of apicomplexan parasites.

Homopolymer tract organization in the human malarial parasite Plasmodium falciparum and related Apicomplexan parasites.

BackgroundHomopolymeric tracts, particularly poly dA.dT, are enriched within the intergenic sequences of eukaryotic genomes where they appear to act as intrinsic regulators of nucleosome positioning. A previous study of the incomplete genome of the human malarial parasite Plasmodium falciparum reports a higher than expected enrichment of poly dA.dT tracts, far above that anticipated even in this highly AT rich genome. Here we report an analysis of the relative frequency, length and spatial arrangement of homopolymer tracts for the complete P. falciparum genome, extending this analysis to twelve additional genomes of Apicomplexan parasites important to human and animal health. In addition, using nucleosome-positioning data available for P. falciparum, we explore the correlation of poly dA.dT tracts with nucleosome-positioning data over key expression landmarks within intergenic regions.ResultsWe describe three apparent lineage-specific patterns of homopolymeric tract organization within the intergenic regions of these Apicomplexan parasites. Moreover, a striking pattern of enrichment of overly long poly dA.dT tracts in the intergenic regions of Plasmodium spp. uniquely extends into protein coding sequences. There is a conserved spatial arrangement of poly dA.dT immediately flanking open reading frames and over predicted core promoter sites. These key landmarks are all relatively depleted in nucleosomes in P. falciparum, as would be expected for poly dA.dT acting as nucleosome exclusion sequences.ConclusionsPrevious comparative studies of homopolymer tract organization emphasize evolutionary diversity; this is the first report of such an analysis within a single phylum. Our data provide insights into the evolution of homopolymeric tracts and the selective pressures at play in their maintenance and expansion.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-848) contains supplementary material, which is available to authorized users.

Reduced representation approaches to interrogate genome diversity in large repetitive plant genomes

Technology and software improvements in the last decade now provide methodologies to access the genome sequence of not only a single accession, but also multiple accessions of plant species. This provides a means to interrogate species diversity at the genome level. Ample diversity among accessions in a collection of species can be found, including single-nucleotide polymorphisms, insertions and deletions, copy number variation and presence/absence variation. For species with small, non-repetitive rich genomes, re-sequencing of query accessions is robust, highly informative, and economically feasible. However, for species with moderate to large sized repetitive-rich genomes, technical and economic barriers prevent en masse genome re-sequencing of accessions. Multiple approaches to access a focused subset of loci in species with larger genomes have been developed, including reduced representation sequencing, exome capture and transcriptome sequencing. Collectively, these approaches have enabled interrogation of diversity on a genome scale for large plant genomes, including crop species important to worldwide food security.

The role of large pedigrees in an era of high-throughput sequencing.

Rare variation is the current frontier in human genetics. The large pedigree design is practical, efficient, and well-suited for investigating rare variation. In large pedigrees, specific rare variants that co-segregate with a trait will occur in sufficient numbers so that effects can be measured, and evidence for association can be evaluated, by making use of methods that fully use the pedigree information. Evidence from linkage analysis can focus investigation, both reducing the multiple testing burden and expanding the variants that can be evaluated and followed up, as recent studies have shown. The large pedigree design requires only a small fraction of the sample size needed to identify rare variants of interest in population-based designs, and many highly suitable, well-understood, and available statistical and computational tools already exist. Samples consisting of large pedigrees with existing rich phenotype and genome scan data should be prime candidates for high-throughput sequencing in the search of the determinants of complex traits.

Insertion site preference of Mu, Tn5, and Tn7 transposons

BackgroundTransposons, segments of DNA that can mobilize to other locations in a genome, are often used for insertion mutagenesis or to generate priming sites for sequencing of large DNA molecules. For both of these uses, a transposon with minimal insertion bias is desired to allow complete coverage with minimal oversampling.FindingsThree transposons, Mu, Tn5, and Tn7, were used to generate insertions in the same set of fosmids containing Candida glabrata genomic DNA. Tn7 demonstrates markedly less insertion bias than either Mu or Tn5, with both Mu and Tn5 biased toward sequences containing guanosine (G) and cytidine (C). This preference of Mu and Tn5 yields less uniform spacing of insertions than for Tn7, in the adenosine (A) and thymidine (T) rich genome of C. glabrata (39% GC).ConclusionsIn light of its more uniform distribution of insertions, Tn7 should be considered for applications in which insertion bias is deleterious.

The SAR11 Group of Alpha-Proteobacteria Is Not Related to the Origin of Mitochondria

Although free living, members of the successful SAR11 group of marine alpha-proteobacteria contain a very small and A+T rich genome, two features that are typical of mitochondria and related obligate intracellular parasites such as the Rickettsiales. Previous phylogenetic analyses have suggested that Candidatus Pelagibacter ubique, the first cultured member of this group, is related to the Rickettsiales+mitochondria clade whereas others disagree with this conclusion. In order to determine the evolutionary position of the SAR11 group and its relationship to the origin of mitochondria, we have performed phylogenetic analyses on the concatenation of 24 proteins from 5 mitochondria and 71 proteobacteria. Our results support that SAR11 group is not the sistergroup of the Rickettsiales+mitochondria clade and confirm that the position of this group in the alpha-proteobacterial tree is strongly affected by tree reconstruction artefacts due to compositional bias. As a consequence, genome reduction and bias toward a high A+T content may have evolved independently in the SAR11 species, which points to a different direction in the quest for the closest relatives to mitochondria and Rickettsiales. In addition, our analyses raise doubts about the monophyly of the newly proposed Pelagibacteraceae family.

Plasmodium falciparum has rare correlation properties

A plot of the correlation function of a given DNA sequence has certain characteristic features common to almost all organisms. One common feature is that the correlation values at distances that are multiples of three is higher than correlation values at other distances. Because of this such a correlation plot can be divided into two or three curves with different scalings. P. falciparum has a rare correlation property which is probably unique. I have analyzed genomes of many bacteria, fungi and protozoa and found that P. falciparum is the only organism whose DNA sequence correlation plot can be divided into four curves with different scalings. This property is neither shared by other species of the Plasmodium genus nor by other AT rich genomes. This could be a hint that the DNA sequence of P. falciparum has undergone certain rare mutational events.