Diversity of 23S rRNA Genes within Individual Prokaryotic Genomes

Background

The concept of ribosomal constraints on rRNA genes is deduced primarily based on the comparison of consensus rRNA sequences between closely related species, but recent advances in whole-genome sequencing allow evaluation of this concept within organisms with multiple rRNA operons.

Methodology/Principal Findings

Using the 23S rRNA gene as an example, we analyzed the diversity among individual rRNA genes within a genome. Of 184 prokaryotic species containing multiple 23S rRNA genes, diversity was observed in 113 (61.4%) genomes (mean 0.40%, range 0.01%–4.04%). Significant (1.17%–4.04%) intragenomic variation was found in 8 species. In 5 of the 8 species, the diversity in the primary structure had only minimal effect on the secondary structure (stem versus loop transition). In the remaining 3 species, the diversity significantly altered local secondary structure, but the alteration appears minimized through complex rearrangement. Intervening sequences (IVS), ranging between 9 and 1471 nt in size, were found in 7 species. IVS in Deinococcus radiodurans and Nostoc sp. encode transposases. T. tengcongensis was the only species in which intragenomic diversity >3% was observed among 4 paralogous 23S rRNA genes.

Conclusions/Significance

These findings indicate tight ribosomal constraints on individual 23S rRNA genes within a genome. Although classification using primary 23S rRNA sequences could be erroneous, significant diversity among paralogous 23S rRNA genes was observed only once in the 184 species analyzed, indicating little overall impact on the mainstream of 23S rRNA gene-based prokaryotic taxonomy.     

Mathematical design of prokaryotic clone-based microarrays

Background  

Clone-based microarrays, on which each spot represents a random genomic fragment, are a good alternative to open reading frame-based microarrays, especially for microorganisms for which the complete genome sequence is not available. Since the generation of a genomic DNA library is a random process, it is beforehand uncertain which genes are represented. Nevertheless, the genome coverage of such an array, which depends on different variables like the insert size and the number of clones in the library, can be predicted by mathematical approaches. When applying the classical formulas that determine the probability that a certain sequence is represented in a DNA library at the nucleotide level, massive amounts of clones would be necessary to obtain a proper coverage of the genome.     

rSW-seq: Algorithm for detection of copy number alterations in deep sequencing data

Background  

Recent advances in sequencing technologies have enabled generation of large-scale genome sequencing data. These data can be used to characterize a variety of genomic features, including the DNA copy number profile of a cancer genome. A robust and reliable method for screening chromosomal alterations would allow a detailed characterization of the cancer genome with unprecedented accuracy.     

Genome size evolution at the speciation level: The cryptic species complex Brachionus plicatilis (Rotifera)

Background  

Studies on genome size variation in animals are rarely done at lower taxonomic levels, e.g., slightly above/below the species level. Yet, such variation might provide important clues on the tempo and mode of genome size evolution. In this study we used the flow-cytometry method to study the evolution of genome size in the rotifer Brachionus plicatilis, a cryptic species complex consisting of at least 14 closely related species.     

Ecological effects of cell-level processes: genome size,functional traits and regional abundance of herbaceous plant species

Background and Aims

Genome size is known to be correlated with a number of phenotypic traits associated with cell sizes and cell-division rates. Genome size was therefore used as a proxy for them in order to assess how common plant traits such as height, specific leaf area and seed size/number predict species regional abundance. In this study it is hypothesized that if there is residual correlation between genome size and abundance after these traits are partialled out, there must be additional ecological effects of cell size and/or cell-division rate.

Methods

Variation in genome size, plant traits and regional abundance were examined in 436 herbaceous species of central European flora, and relationships were sought for among these variables by correlation and path analysis.

Key Results

Species regional abundance was weakly but significantly correlated with genome size; the relationship was stronger for annuals (R² = 0·145) than for perennials (R² = 0·027). In annuals, genome size was linked to abundance via its effect on seed size, which constrains seed number and hence population growth rate. In perennials, it weakly affected (via height and specific leaf area) competitive ability. These relationships did not change qualitatively after phylogenetic correction. In both annuals and perennials there was an unresolved effect of genome size on abundance.

Conclusions

The findings indicate that additional predictors of regional abundance should be sought among variables that are linked to cell size and cell-division rate. Signals of these cell-level processes remain identifiable even at the landscape scale, and show deep differences between perennials and annuals. Plant population biology could thus possibly benefit from more systematic use of indicators of cell-level processes.     

Genome size differentiates co-occurring populations of the planktonic diatom Ditylum brightwellii (Bacillariophyta)

Background  

Diatoms are one of the most species-rich groups of eukaryotic microbes known. Diatoms are also the only group of eukaryotic micro-algae with a diplontic life history, suggesting that the ancestral diatom switched to a life history dominated by a duplicated genome. A key mechanism of speciation among diatoms could be a propensity for additional stable genome duplications. Across eukaryotic taxa, genome size is directly correlated to cell size and inversely correlated to physiological rates. Differences in relative genome size, cell size, and acclimated growth rates were analyzed in isolates of the diatom Ditylum brightwellii. Ditylum brightwellii consists of two main populations with identical 18s rDNA sequences; one population is distributed globally at temperate latitudes and the second appears to be localized to the Pacific Northwest coast of the USA. These two populations co-occur within the Puget Sound estuary of WA, USA, although their peak abundances differ depending on local conditions.     

False positive PCR detection of <Emphasis Type="Italic">Tropheryma whipplei</Emphasis> in the saliva of healthy people

Background  

Tropheryma whipplei, the agent of Whipple's disease (WD), has been recently isolated and the genomes of two isolates have been fully sequenced. Previous diagnosis tools for the diagnosis of the disease used sequence analysis of the 16S rRNA gene. Using this target gene, the high percentage of detection of the bacterium in saliva of healthy people was in contrast to the negative results obtained with specific target genes. The aim of our study was to compare previously published primers targeting the 16S rRNA gene to real-time PCR with Taqman* probes targeting specific repeat genes only found in the genome of T. whipplei in a series of 57 saliva from healthy people.     

Genome increase as a clock for the origin and evolution of life

Background  

The size of non-redundant functional genome can be an indicator of biological complexity of living organisms. Several positive feedback mechanisms including gene cooperation and duplication with subsequent specialization may result in the exponential growth of biological complexity in macro-evolution.     

Bayesian DNA copy number analysis

Background  

Some diseases, like tumors, can be related to chromosomal aberrations, leading to changes of DNA copy number. The copy number of an aberrant genome can be represented as a piecewise constant function, since it can exhibit regions of deletions or gains. Instead, in a healthy cell the copy number is two because we inherit one copy of each chromosome from each our parents.     

Evolution of ribosomal DNA-derived satellite repeat in tomato genome

Background  

Tandemly repeated DNA, also called as satellite DNA, is a common feature of eukaryotic genomes. Satellite repeats can expand and contract dramatically, which may cause genome size variation among genetically-related species. However, the origin and expansion mechanism are not clear yet and needed to be elucidated.     

Large-scale polymorphism of heterochromatic repeats in the DNA of <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Background  

The composition of the individual eukaryote's genome and its variation within a species remain poorly defined. Even for a sequenced genome such as that of the model plant Arabidopsis thaliana accession Col-0, the large arrays of heterochromatic repeats are incompletely sequenced, with gaps of uncertain size persisting in them.     

Comparative assessment of performance and genome dependence among phylogenetic profiling methods

Background  

The rapidly increasing speed with which genome sequence data can be generated will be accompanied by an exponential increase in the number of sequenced eukaryotes. With the increasing number of sequenced eukaryotic genomes comes a need for bioinformatic techniques to aid in functional annotation. Ideally, genome context based techniques such as proximity, fusion, and phylogenetic profiling, which have been so successful in prokaryotes, could be utilized in eukaryotes. Here we explore the application of phylogenetic profiling, a method that exploits the evolutionary co-occurrence of genes in the assignment of functional linkages, to eukaryotic genomes.     

SeeGH – A software tool for visualization of whole genome array comparative genomic hybridization data

Background  

Array comparative genomic hybridization (CGH) is a technique which detects copy number differences in DNA segments. Complete sequencing of the human genome and the development of an array representing a tiling set of tens of thousands of DNA segments spanning the entire human genome has made high resolution copy number analysis throughout the genome possible. Since array CGH provides signal ratio for each DNA segment, visualization would require the reassembly of individual data points into chromosome profiles.     

Genome size reduction can trigger rapid phenotypic evolution in invasive plants

Background and Aims

The study of rapid evolution in invasive species has highlighted the fundamental role played by founder events, emergence of genetic novelties through recombination and rapid response to new selective pressures. However, whether rapid adaptation of introduced species can be driven by punctual changes in genome organization has received little attention. In plants, variation in genome size, i.e. variation in the amount of DNA per monoploid set of chromosomes through loss or gain of repeated DNA sequences, is known to influence a number of physiological, phenological and life-history features. The present study investigated whether change in genome size has contributed to the evolution of greater potential of vegetative growth in invasive populations of an introduced grass.

Methods

The study was based on the recent demonstration that invasive genotypes of reed canarygrass (Phalaris arundinacea) occurring in North America have emerged from recombination between introduced European strains. The genome sizes of more than 200 invasive and native genotypes were measured and their genome size was related to their phenotypic traits measured in a common glasshouse environment. Population genetics data were used to infer phylogeographical relationships between study populations, and the evolutionary history of genome size within the study species was inferred.

Key Results

Invasive genotypes had a smaller genome than European native genotypes from which they are derived. This smaller genome size had phenotypic effects that increased the species'' invasive potential, including a higher early growth rate, due to a negative relationship between genome size and rate of stem elongation. Based on inferred phylogeographical relationships of invasive and native populations, evolutionary models were consistent with a scenario of genome reduction by natural selection during the invasion process, rather than a scenario of stochastic change.

Conclusions

Punctual reduction in genome size could cause rapid changes in key phenotypic traits that enhance invasive ability. Although the generality of genome size variation leading to phenotypic evolution and the specific genomic mechanisms involved are not known, change in genome size may constitute an important but previously under-appreciated mechanism of rapid evolutionary change that may promote evolutionary novelties over short time scales.Key words: Biological invasion, evolutionary models, genome size, Phalaris arundinacea, quantile regression, relative growth rate, rapid evolution     

Draft genome sequence of the Daphnia pathogen Octosporea bayeri: insights into the gene content of a large microsporidian genome and a model for host-parasite interactions

Background  

The highly compacted 2.9-Mb genome of Encephalitozoon cuniculi placed the microsporidia in the spotlight, encoding a mere 2,000 proteins and a highly reduced suite of biochemical pathways. This extreme level of reduction is not universal across the microsporidia, with genomes known to vary up to sixfold in size, suggesting that some genomes may harbor a gene content that is not as reduced as that of Enc. cuniculi. In this study, we present an in-depth survey of the large genome of Octosporea bayeri, a pathogen of Daphnia magna, with an estimated genome size of 24 Mb, in order to shed light on the organization and content of a large microsporidian genome.     

Gene family evolution: an in-depth theoretical and simulation analysis of non-linear birth-death-innovation models

Background  

The size distribution of gene families in a broad range of genomes is well approximated by a generalized Pareto function. Evolution of ensembles of gene families can be described with Birth, Death, and Innovation Models (BDIMs). Analysis of the properties of different versions of BDIMs has the potential of revealing important features of genome evolution.     

Evolutionary constraints on yeast protein size

Background  

Despite a strong evolutionary pressure to reduce genome size, proteins vary in length over a surprisingly wide range also in very compact genomes. Here we investigated the evolutionary forces that act on protein size in the yeast Saccharomyces cerevisiae utilizing a system-wide bioinformatics approach. Data on yeast protein size was compared to global experimental data on protein expression, phenotypic pleiotropy, protein-protein interactions, protein evolutionary rate and biochemical classification.     

Accelerated evolution associated with genome reduction in a free-living prokaryote

Background

Three complete genomes of Prochlorococcus species, the smallest and most abundant photosynthetic organism in the ocean, have recently been published. Comparative genome analyses reveal that genome shrinkage has occurred within this genus, associated with a sharp reduction in G+C content. As all examples of genome reduction characterized so far have been restricted to endosymbionts or pathogens, with a host-dependent lifestyle, the observed genome reduction in Prochlorococcus is the first documented example of such a process in a free-living organism.

Results

Our results clearly indicate that genome reduction has been accompanied by an increased rate of protein evolution in P. marinus SS120 that is even more pronounced in P. marinus MED4. This acceleration has affected every functional category of protein-coding genes. In contrast, the 16S rRNA gene seems to have evolved clock-like in this genus. We observed that MED4 and SS120 have lost several DNA-repair genes, the absence of which could be related to the mutational bias and the acceleration of amino-acid substitution.

Conclusions

We have examined the evolutionary mechanisms involved in this process, which are different from those known from host-dependent organisms. Indeed, most substitutions that have occurred in Prochlorococcus have to be selectively neutral, as the large size of populations imposes low genetic drift and strong purifying selection. We assume that the major driving force behind genome reduction within the Prochlorococcus radiation has been a selective process favoring the adaptation of this organism to its environment. A scenario is proposed for genome evolution in this genus.     

QTL mapping of agronomic traits in tef [<Emphasis Type="Italic">Eragrostis tef</Emphasis>(Zucc) Trotter]

Background  

Tef [Eragrostis tef (Zucc.) Trotter] is the major cereal crop in Ethiopia. Tef is an allotetraploid with a base chromosome number of 10 (2n = 4× = 40) and a genome size of 730 Mbp. The goal of this study was to identify agronomically important quantitative trait loci (QTL) using recombinant inbred lines (RIL) derived from an inter-specific cross between E. tef and E. pilosa (30-5).     

GenomeMatcher: A graphical user interface for DNA sequence comparison

Background  

The number of available genome sequences is increasing, and easy-to-use software that enables efficient comparative analysis is needed.