The large genome constraint hypothesis: evolution, ecology and phenotype

BACKGROUND AND AIMS: If large genomes are truly saturated with unnecessary 'junk' DNA, it would seem natural that there would be costs associated ith accumulation and replication of this excess DNA. Here we examine the available evidence to support this hypothesis, which we term the 'large genome constraint'. We examine the large genome constraint at three scales: evolution, ecology, and the plant phenotype. SCOPE: In evolution, we tested the hypothesis that plant lineages with large genomes are diversifying more slowly. We found that genera with large genomes are less likely to be highly specious -- suggesting a large genome constraint on speciation. In ecology, we found that species with large genomes are under-represented in extreme environments -- again suggesting a large genome constraint for the distribution and abundance of species. Ultimately, if these ecological and evolutionary constraints are real, the genome size effect must be expressed in the phenotype and confer selective disadvantages. Therefore, in phenotype, we review data on the physiological correlates of genome size, and present new analyses involving maximum photosynthetic rate and specific leaf area. Most notably, we found that species with large genomes have reduced maximum photosynthetic rates - again suggesting a large genome constraint on plant performance. Finally, we discuss whether these phenotypic correlations may help explain why species with large genomes are trimmed from the evolutionary tree and have restricted ecological distributions. CONCLUSION: Our review tentatively supports the large genome constraint hypothesis.     

Gene and genome duplications in vertebrates: the one-to-four (-to-eight in fish) rule and the evolution of novel gene functions

One important mechanism for functional innovation during evolution is the duplication of genes and entire genomes. Evidence is accumulating that during the evolution of vertebrates from early deuterostome ancestors entire genomes were duplicated through two rounds of duplications (the 'one-to-two-to-four' rule). The first genome duplication in chordate evolution might predate the Cambrian explosion. The second genome duplication possibly dates back to the early Devonian. Recent data suggest that later in the Devonian, the fish genome was duplicated for a third time to produce up to eight copies of the original deuterostome genome. This last duplication took place after the two major radiations of jawed vertebrate life, the ray-finned fish (Actinopterygia) and the sarcopterygian lineage, diverged. Therefore the sarcopterygian fish, which includes the coelacanth, lungfish and all land vertebrates such as amphibians, reptiles, birds and mammals, tend to have only half the number of genes compared with actinopterygian fish. Although many duplicated genes turned into pseudogenes, or even 'junk' DNA, many others evolved new functions particularly during development. The increased genetic complexity of fish might reflect their evolutionary success and diversity.     

Discovery of a large number of previously unrecognized mitochondrial pseudogenes in fish genomes

Nuclear inserted copies of mitochondrial origin (numts) vary widely among eukaryotes, with human and plant genomes harboring the largest repertoires. Numts were previously thought to be absent from fish species, but the recent release of three fish nuclear genome sequences provides the resource to obtain a more comprehensive insight into the extent of mtDNA transfer in fishes. From the sequence analyses of the genomes of Fugu rubripes, Tetraodon nigroviridis, and Danio rerio, we have identified 2, 5, and 10 recent numt integrations, respectively, which integrated into those genomes less than 0.6 million years (Myr) ago. Such results contradict the hypothesis of absence or rarity of numts in fishes, as (i) the ratio of numts to the total size of the nuclear genome in T. nigroviridis was superior to the ratio observed in several higher vertebrate species (e.g., chicken, mouse, and rat), and only surpassed by humans, and (ii) the mtDNA coverage transferred to the nuclear genome of D. rerio is exceeded only by human and mouse, within the whole range of eukaryotic genomes surveyed for numts. Additionally, 335, 336, and 471 old numts (>12.5 Myr) were detected in F. rubripes, T. nigroviridis, and D. rerio, respectively. Surprisingly, old numts are inserted preferentially into known or predicted genes, as inferred for recent numts in human. However, because in fish genomes such integrations are old, they are likely to represent evolutionary successes and they may be considered a potential important evolutionary mechanism for the enhancement of genomic coding regions.     

Genome size evolution of the extant lycophytes and ferns

Ferns and lycophytes have remarkably large genomes. However, little is known about how their genome size evolved in fern lineages. To explore the origins and evolution of chromosome numbers and genome size in ferns, we used flow cytometry to measure the genomes of 240 species (255 samples) of extant ferns and lycophytes comprising 27 families and 72 genera, of which 228 species (242 samples) represent new reports. We analyzed correlations among genome size, spore size, chromosomal features, phylogeny, and habitat type preference within a phylogenetic framework. We also applied ANOVA and multinomial logistic regression analysis to preference of habitat type and genome size. Using the phylogeny, we conducted ancestral character reconstruction for habitat types and tested whether genome size changes simultaneously with shifts in habitat preference. We found that 2C values had weak phylogenetic signal, whereas the base number of chromosomes (x) had a strong phylogenetic signal. Furthermore, our analyses revealed a positive correlation between genome size and chromosome traits, indicating that the base number of chromosomes (x), chromosome size, and polyploidization may be primary contributors to genome expansion in ferns and lycophytes. Genome sizes in different habitat types varied significantly and were significantly correlated with habitat types; specifically, multinomial logistic regression indicated that species with larger 2C values were more likely to be epiphytes. Terrestrial habitat is inferred to be ancestral for both extant ferns and lycophytes, whereas transitions to other habitat types occurred as the major clades emerged. Shifts in habitat types appear be followed by periods of genomic stability. Based on these results, we inferred that habitat type changes and multiple whole-genome duplications have contributed to the formation of large genomes of ferns and their allies during their evolutionary history.     

Palaeogenomics of pterosaurs and the evolution of small genome size in flying vertebrates

The two living groups of flying vertebrates, birds and bats, both have constricted genome sizes compared with their close relatives. But nothing is known about the genomic characteristics of pterosaurs, which took to the air over 70 Myr before birds and were the first group of vertebrates to evolve powered flight. Here, we estimate genome size for four species of pterosaurs and seven species of basal archosauromorphs using a Bayesian comparative approach. Our results suggest that small genomes commonly associated with flight in bats and birds also evolved in pterosaurs, and that the rate of genome-size evolution is proportional to genome size within amniotes, with the fastest rates occurring in lineages with the largest genomes. We examine the role that drift may have played in the evolution of genome size within tetrapods by testing for correlated evolution between genome size and body size, but find no support for this hypothesis. By contrast, we find evidence suggesting that a combination of adaptation and phylogenetic inertia best explains the correlated evolution of flight and genome-size contraction. These results suggest that small genome/cell size evolved prior to or concurrently with flight in pterosaurs. We predict that, similar to the pattern seen in theropod dinosaurs, genome-size contraction preceded flight in pterosaurs and bats.     

Recombination Rate and AT-content Show Opposite Correlations in Mammalian and Other Animal Genomes

Genomes are often characterized by their genome size, recombination rate and nucleotide content. Negative correlations of recombination rate and AT-content have been reported for selected species. Here we investigate 26 completely sequenced animal genomes with respect to their genomic characteristics. After correcting for phylogenetic correlations we find significant correlations between these characters, showing that these do not evolve independently. The genomes of mammals differ when compared to other animal genomes. Correlations between AT-content and recombination rate are negative, but only for mammals. All other animals show a positive correlation for these characters. These results show that studies on the relationship between genome characteristics using single species might lead to the right conclusion, but cannot be applied as a general rule to all taxonomic groups. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

鱼类特异的基因组复制

辐鳍鱼类是脊椎动物中种类最多、分布最广的类群,其基因组大小不等。过去的观点认为,在脊椎动物进化历程中曾发生了两次基因组复制。近期的系统基因组学研究资料进一步提出,在大约350百万年,辐鳍鱼还发生了第三次基因组复制,即鱼类特异的基因组复制(fish-specificgenomeduplication,FSGD),且发生的时间正处在“物种极度丰富”的硬骨鱼谱系(真骨总目)和“物种贫乏”的谱系(辐鳍鱼纲基部的类群)出现分歧的时间点,表明FSGD与硬骨鱼物种和生物多样性的增加有关。进一步开展鱼类比较基因组学和功能基因组学研究将进一步验证FSGD这一假说。     

Protozoan genomes: gene identification and annotation

The draft sequence of several complete protozoan genomes is now available and genome projects are ongoing for a number of other species. Different strategies are being implemented to identify and annotate protein coding and RNA genes in these genomes, as well as study their genomic architecture. Since the genomes vary greatly in size, GC-content, nucleotide composition, and degree of repetitiveness, genome structure is often a factor in choosing the methodology utilised for annotation. In addition, the approach taken is dictated, to a greater or lesser extent, by the particular reasons for carrying out genome-wide analyses and the level of funding available for projects. Nevertheless, these projects have provided a plethora of material that will aid in understanding the biology and evolution of these parasites, as well as identifying new targets that can be used to design urgently required drug treatments for the diseases they cause.     

Reconstruction of ancient genome and gene order from complete microbial genome sequences

Microbial genome sequences provide us with the fossil records for inferring their origination and evolution. Assuming that current microbial genomes are the evolutionary results of ancient genomes or fragments and the neighboring genes in ancient genomes are more likely neighbors in current genomes, in this paper we proposed a paleontological algorithm and assembled the orthologous gene groups from 66 complete and current microbial genome sequences into a pseudo-ancient genome, which consists of continuous fragments of various sizes. We performed bootstrap resampling and correlation analyses and the results showed that the assembled ancient genome and fragments are statistically significant and the genes of the same fragment are inherently related and likely derived from common ancestors. This method provides a new computational tool for studying microbial genome structure and evolution.     

The repair of double-strand breaks in plants: mechanisms and consequences for genome evolution

The efficient repair of double-strand breaks (DSBs) in genomic DNA is important for the survival of all organisms. In recent years, basic mechanisms of DSB repair in somatic plant cells have been elucidated. DSBs are mainly repaired by non-homologous end-joining (NHEJ). The repair can be associated with deletions, but also insertions due to copying genomic sequences from elsewhere into the break. Species-specific differences of NHEJ have been reported and an inverse correlation of deletion size to genome size has been postulated, indicating that NHEJ might contribute significantly to evolution of genome size. DSB repair by homologous recombination (HR) might also influence genome organization. Whereas homology present in an allelic or an ectopic position is hardly used for repair, the use of homologous sequences in close proximity to the break is frequent. A 'single-strand annealing' mechanism that leads to sequence deletions between direct repeats is particularly efficient. This might explain the accumulation of single long terminal repeats of retroelements in cereal genomes. The conservative 'synthesis-dependent strand annealing' mechanism, resulting in conversions without crossovers is also prominent and seems to be significant for the evolution of tandemly arranged gene families such as resistance genes. Induction of DSBs could be used as a means for the controlled manipulation of plant genomes in an analogous way for the use of marker gene excision and site-specific integration.     

Genome size reduction in the chicken has involved massive loss of ancestral protein-coding genes

Both mean genomes size and the variance in genome size among species are smaller on average in birds (class Aves) than in the other tetrapod classes. In order to test whether loss of protein-coding genes has contributed to genome size reduction in birds, we compared the chicken genome and five mammalian genomes. Numbers of members (paralogs) were significantly lower in the chicken gene families than in the corresponding mammalian families. Phylogenetic analyses of chicken, mammal, and fish paralogs supported the hypothesis that chicken-specific loss of paralogs occurred much more frequently than mammal-specific gene duplications. Moreover, the phylogenetic analyses supported the hypothesis that a substantial majority of the paralogs lost in chicken originated from duplications prior to the most recent common ancestor of tetrapods and bony fishes. In addition to loss of paralogs, numerous gene families present in the mammalian genomes were missing in the chicken genome; over 1,000 of these families were found in bony fishes, implying presence of the family in the tetrapod ancestor. In the set of families with more members on average in the mammals than in the chicken, immune system function was associated with a greater degree of gene family size reduction in the chicken, consistent with other evidence that immune system gene families have become particularly compact in birds.     

In silico comparative analysis reveals a mosaic conservation of genes within a novel colinear region in wheat chromosome 1AS and rice chromosome 5S

Comparative RFLP mapping has revealed extensive conservation of marker order in different grass genomes. However, microcolinearity studies at the sequence level have shown rapid genome evolution and many exceptions to colinearity. Most of these studies have focused on a limited size of genomic fragment and the extent of microcolinearity over large distances or across entire genomes remains poorly characterized in grasses. Here, we have investigated the microcolinearity between the rice genome and a total of 1,500 kb from physical BAC contigs on wheat chromosome 1AS. Using ESTs mapped in wheat chromosome bins as an additional source of physical data, we have identified 27 conserved orthologous sequences between wheat chromosome 1AS and a region of 1,210 kb located on rice chromosome 5S. Our results extend the orthology described earlier between wheat chromosome group 1S and rice chromosome 5S. Microcolinearity was found to be frequently disrupted by rearrangements which must have occurred after the divergence of wheat and rice. At the Lr10 orthologous loci, microrearrangements were due to the insertion of mobile elements, but also originated from gene movement, amplification, deletion and inversion. These mechanisms of genome evolution are at the origin of the mosaic conservation observed between the orthologous regions. Finally, in silico mapping of wheat genes identified an intragenomic colinearity between fragments from rice chromosome 1L and 5S, suggesting an ancestral segmental duplication in rice.Electronic Supplementary Material Supplementary material is available in the online version of this article at     

LTR retrotransposons and flowering plant genome size: emergence of the increase/decrease model

Long Terminal Repeat (LTR) retrotransposons are ubiquitous components of plant genomes. Because of their copy-and-paste mode of transposition, these elements tend to increase their copy number while they are active. In addition, it is now well established that the differences in genome size observed in the plant kingdom are accompanied by variations in LTR retrotransposon content, suggesting that LTR retrotransposons might be important players in the evolution of plant genome size, along with polyploidy. The recent availability of large genomic sequences for many crop species has made it possible to examine in detail how LTR retrotransposons actually drive genomic changes in plants. In the present paper, we provide a review of the recent publications that have contributed to the knowledge of plant LTR retrotransposons, as structural components of the genomes, as well as from an evolutionary genomic perspective. These studies have shown that plant genomes undergo genome size increases through bursts of retrotransposition, while there is a counteracting process that tends to eliminate the transposed copies from the genomes. This process involves recombination mechanisms that occur either between the LTRs of the elements, leading to the formation of solo-LTRs, or between direct repeats anywhere in the sequence of the element, leading to internal deletions. All these studies have led to the emergence of a new model for plant genome evolution that takes into account both genome size increases (through retrotransposition) and decreases (through solo-LTR and deletion formation). In the conclusion, we discuss this new model and present the future prospects in the study of plant genome evolution in relation to the activity of transposable elements.     

Genome size is negatively correlated with effective population size in ray-finned fish

A recent theory suggesting that genome size and complexity can increase as a passive consequence of small effective population size has generated much controversy. In this article, we demonstrate that freshwater fish species, which have smaller effective population sizes than marine fish species, have larger genomes. We show that genome size is negatively correlated with genetic variability, independent of phylogeny, body size and generation time. Genome duplication is also observed predominantly in freshwater fish. These results suggest that the raw materials of complexity originate under conditions of reduced selection efficiency.     

Nucleomorph and plastid genome sequences of the chlorarachniophyte Lotharella oceanica: convergent reductive evolution and frequent recombination in nucleomorph-bearing algae

Background

Nucleomorphs are residual nuclei derived from eukaryotic endosymbionts in chlorarachniophyte and cryptophyte algae. The endosymbionts that gave rise to nucleomorphs and plastids in these two algal groups were green and red algae, respectively. Despite their independent origin, the chlorarachniophyte and cryptophyte nucleomorph genomes share similar genomic features such as extreme size reduction and a three-chromosome architecture. This suggests that similar reductive evolutionary forces have acted to shape the nucleomorph genomes in the two groups. Thus far, however, only a single chlorarachniophyte nucleomorph and plastid genome has been sequenced, making broad evolutionary inferences within the chlorarachniophytes and between chlorarachniophytes and cryptophytes difficult. We have sequenced the nucleomorph and plastid genomes of the chlorarachniophyte Lotharella oceanica in order to gain insight into nucleomorph and plastid genome diversity and evolution.

Results

The L. oceanica nucleomorph genome was found to consist of three linear chromosomes totaling ~610 kilobase pairs (kbp), much larger than the 373 kbp nucleomorph genome of the model chlorarachniophyte Bigelowiella natans. The L. oceanica plastid genome is 71 kbp in size, similar to that of B. natans. Unexpectedly long (~35 kbp) sub-telomeric repeat regions were identified in the L. oceanica nucleomorph genome; internal multi-copy regions were also detected. Gene content analyses revealed that nucleomorph house-keeping genes and spliceosomal intron positions are well conserved between the L. oceanica and B. natans nucleomorph genomes. More broadly, gene retention patterns were found to be similar between nucleomorph genomes in chlorarachniophytes and cryptophytes. Chlorarachniophyte plastid genomes showed near identical protein coding gene complements as well as a high level of synteny.

Conclusions

We have provided insight into the process of nucleomorph genome evolution by elucidating the fine-scale dynamics of sub-telomeric repeat regions. Homologous recombination at the chromosome ends appears to be frequent, serving to expand and contract nucleomorph genome size. The main factor influencing nucleomorph genome size variation between different chlorarachniophyte species appears to be expansion-contraction of these telomere-associated repeats rather than changes in the number of unique protein coding genes. The dynamic nature of chlorarachniophyte nucleomorph genomes lies in stark contrast to their plastid genomes, which appear to be highly stable in terms of gene content and synteny.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-374) contains supplementary material, which is available to authorized users.     

A new look at genome size,evolutionary duration and genetic variation in salamanders

Salamanders have some of the largest genomes among vertebrates, and also some of the lowest reported levels of genetic diversity. Paedomorphs, in particular, have the largest genomes on average among urodela, and display exceptionally low levels of nucleotide and protein variation. Here, we address the question of genetic variation in relation to genome size in eight different salamander families. Using the rag1 gene as a probe for evolutionary rates, we found that rates of substitution are exceptionally low in obligate paedomorphs (neotenes) and other salamander species. Substitution rates in some cases are as low as those reported for cartilaginous fish, which have the slowest mutation rates recorded so far in vertebrates. Confirming and extending an earlier study, we also found that genome size is correlated with phylogenetic age in Plethodontidae, indicating a more general trend in genome size evolution in urodela. The Plethodontidae, furthermore, display much higher levels of genetic variance than the obligate neotene families, consistent with greater habitat heterogeneity in terrestrial salamanders. Finally, we present the first direct evidence of a gene, rag1, whose substitution rate is negatively associated with genome size. Based on these and other observations, we propose a hypothesis according to which mutation rates in nuclear genes tend to increase as genome size decreases during the course of vertebrate evolution.     

Modelling the evolution of genomes with integrated external and internal functions

The genomes that organisms transmit between generations contain information about different kinds of functions. The genome with the "best" mix and number of genes for these functions is the one that natural selection favours. Here I introduce a new way to model simple organisms with genes for external and internal functions, and use it to study the evolution of genome size. The external functions are exemplified by resource use and the internal functions by mutation control (repair). It is shown that even with a suitable proportion of genes for mutation control, the genomes in the organisms do not forever incorporate genes that increase resource use. Instead they evolve towards an optimal genome of limited size. The optimal proportion of genes for mutation control is shown to have an upper limit given by the ease with which transmission accuracy is improved by adding extra genes for this purpose to the genome. The model illustrates how natural selection on genomes integrates systems for the transmission of genetic information with systems relating to the external adaptation of the organism. It also opens up for other, more detailed theoretical investigations of genome functions.     

Mutation rate and genome reduction in endosymbiotic and free-living bacteria

Genome reduction has been considered the hallmark of endosymbiotic bacteria, such as endocellular mutualists or obligatory pathogens until it was found exactly the same in several free-living bacteria. In endosymbiotic bacteria genome reduction is mainly attributed to degenerative processes due to small population size. These cannot affect the free-living bacteria with reduced genomes because they are known to have very large population sizes. It has been proposed that selection for simplification drove genome reduction in these free-living bacteria. For at least one of them (Prochlorococcus), genome reduction is associated with accelerated evolution and we suggest an alternative hypothesis based on increase in mutation rate as the primary cause of genome reduction in free-living bacteria.     

Structural organization of the human gastrointestinal glutathione peroxidase (GPX2) promoter and 3'-nontranscribed region: transcriptional response to exogenous redox agents

We describe a new family of repetitive elements, named Mimo, from the mosquito Culex pipiens. Structural characteristics of these elements fit well with those of miniature inverted-repeat transposable elements (MITEs), which are ubiquitous and highly abundant in plant genomes. The occurrence of Mimo in C. pipiens provides new evidence that MITEs are not restricted to plant genomes, but may be widespread in arthropods as well. The copy number of Mimo elements in C. pipiens (1000 copies in a 540 Mb genome) supports the hypothesis that there is a positive correlation between genome size and the magnitude of MITE proliferation. In contrast to most MITE families described so far, members of the Mimo family share a high sequence conservation, which may reflect a recent amplification history in this species. In addition, we found that Mimo elements are a frequent nest for other MITE-like elements, suggesting that multiple and successive MITE transposition events have occurred very recently in the C. pipiens genome. Despite evidence for recent mobility of these MITEs, no element has been found to encode a protein; therefore, we do not know how they have transposed and have spread in the genome. However, some sequence similarities in terminal inverted-repeats suggest a possible filiation of some of these mosquito MITEs with pogo-like DNA transposons.     

Sampling,sequencing and the SAD

The Species Abundance Distribution (SAD) is a common metric for characterizing macroscopic ecological communities. Recently, this metric has been applied to analysis of microbial communities as well. However, as compared to macroscopic communities, sampling of microscopic communities is different. In particular, most microbial communities are studied using sequencing techniques. These techniques have known biases that result in certain taxa being detected more often than others, even if the taxa are present in the sample at equivalent abundances. There are, for example, amplification biases that result in some sequences being amplified more than others. Likewise, differences in genome size across organisms can result in different numbers of reads from different taxa, again resulting in biased detection. A number of bioinformatics methods have been devised to account for biases in sequencing data, allowing for more accurate estimates of relative taxon abundances. However, because the sampling process itself is affected by biased detection, and because sampling (and under-sampling in particular) can influence the shape of the SAD, it is possible that, even when corrected for through re-scaling, detection biases can affect SAD predictions from sequencing data. To test this hypothesis, we construct a simulation model of the sampling process, focusing on biased detection in shotgun sequencing that arises from genome size differences across microbial taxa. Interestingly, we find that, although genome size itself does not impact SAD predictions, predictions can vary depending on the range of genome sizes that are represented in a community, as well as how genome size is distributed (i.e., whether the majority of species have small versus large genomes). Our results suggest that care should be taken when comparing SADs across environments, particularly when those environments might have taxa with different genome size distributions. Furthermore, our results indicate that relatively deep sequencing might be required to avoid drawing spurious inferences about ecological differences across microbial communities.