Recombination drives vertebrate genome contraction

Selective and/or neutral processes may govern variation in DNA content and, ultimately, genome size. The observation in several organisms of a negative correlation between recombination rate and intron size could be compatible with a neutral model in which recombination is mutagenic for length changes. We used whole-genome data on small insertions and deletions within transposable elements from chicken and zebra finch to demonstrate clear links between recombination rate and a number of attributes of reduced DNA content. Recombination rate was negatively correlated with the length of introns, transposable elements, and intergenic spacer and with the rate of short insertions. Importantly, it was positively correlated with gene density, the rate of short deletions, the deletion bias, and the net change in sequence length. All these observations point at a pattern of more condensed genome structure in regions of high recombination. Based on the observed rates of small insertions and deletions and assuming that these rates are representative for the whole genome, we estimate that the genome of the most recent common ancestor of birds and lizards has lost nearly 20% of its DNA content up until the present. Expansion of transposable elements can counteract the effect of deletions in an equilibrium mutation model; however, since the activity of transposable elements has been low in the avian lineage, the deletion bias is likely to have had a significant effect on genome size evolution in dinosaurs and birds, contributing to the maintenance of a small genome. We also demonstrate that most of the observed correlations between recombination rate and genome contraction parameters are seen in the human genome, including for segregating indel polymorphisms. Our data are compatible with a neutral model in which recombination drives vertebrate genome size evolution and gives no direct support for a role of natural selection in this process.     

Rapid sequence turnover at an intergenic locus in Drosophila

Closely related species of Drosophila tend to have similar genome sizes. The strong imbalance in favor of small deletions relative to insertions implies that the unconstrained DNA in Drosophila is unlikely to be passively inherited from even closely related ancestors, and yet most DNA in Drosophila genomes is intergenic and potentially unconstrained. In an attempt to investigate the maintenance of this intergenic DNA, we studied the evolution of an intergenic locus on the fourth chromosome of the Drosophila melanogaster genome. This 1.2-kb locus is marked by two distinct, large insertion events: a nuclear transposition of a mitochondrial sequence and a transposition of a nonautonomous DNA transposon DNAREP1_DM. Because we could trace the evolutionary histories of these sequences, we were able to reconstruct the length evolution of this region in some detail. We sequenced this locus in all four species of the D. melanogaster species complex: D. melanogaster, D. simulans, D. sechellia, and D. mauritiana. Although this locus is similar in size in these four species, less than 10% of the sequence from the most recent common ancestor remains in D. melanogaster and all of its sister species. This region appears to have increased in size through several distinct insertions in the ancestor of the D. melanogaster species complex and has been shrinking since the split of these lineages. In addition, we found no evidence suggesting that the size of this locus has been maintained over evolutionary time; these results are consistent with the model of a dynamic equilibrium between persistent DNA loss through small deletions and more sporadic DNA gain through less frequent but longer insertions. The apparent stability of genome size in Drosophila may belie very rapid sequence turnover at intergenic loci.     

High rate of DNA loss in the Drosophila melanogaster and Drosophila virilis species groups

We recently proposed that patterns of evolution of non-LTR retrotransposable elements can be used to study patterns of spontaneous mutation. Transposition of non-LTR retrotransposable elements commonly results in creation of 5' truncated, "dead-on-arrival" copies. These inactive copies are effectively pseudogenes and, according to the neutral theory, their molecular evolution ought to reflect rates and patterns of spontaneous mutation. Maximum parsimony can be used to separate the evolution of active lineages of a non-LTR element from the fate of the "dead-on-arrival" insertions and to directly assess the relative frequencies of different types of spontaneous mutations. We applied this approach using a non-LTR element, Helena, in the Drosophila virilis group and have demonstrated a surprisingly high incidence of large deletions and the virtual absence of insertions. Based on these results, we suggested that Drosophila in general may exhibit a high rate of spontaneous large deletions and have hypothesized that such a high rate of DNA loss may help to explain the puzzling dearth of bona fide pseudogenes in Drosophila. We also speculated that variation in the rate of spontaneous deletion may contribute to the divergence of genome size in different taxa by affecting the amount of superfluous "junk" DNA such as, for example, pseudogenes or long introns. In this paper, we extend our analysis to the D. melanogaster subgroup, which last shared a common ancestor with the D. virilis group approximately 40 MYA. In a different region of the same transposable element, Helena, we demonstrate that inactive copies accumulate deletions in species of the D. melanogaster subgroup at a rate very similar to that of the D. virilis group. These results strongly suggest that the high rate of DNA loss is a general feature of Drosophila and not a peculiar property of a particular stretch of DNA in a particular species group.      

A Model for Genome Size Evolution

We present a model for genome size evolution that takes into account both local mutations such as small insertions and small deletions, and large chromosomal rearrangements such as duplications and large deletions. We introduce the possibility of undergoing several mutations within one generation. The model, albeit minimalist, reveals a non-trivial spontaneous dynamics of genome size: in the absence of selection, an arbitrary large part of genomes remains beneath a finite size, even for a duplication rate 2.6-fold higher than the rate of large deletions, and even if there is also a systematic bias toward small insertions compared to small deletions. Specifically, we show that the condition of existence of an asymptotic stationary distribution for genome size non-trivially depends on the rates and mean sizes of the different mutation types. We also give upper bounds for the median and other quantiles of the genome size distribution, and argue that these bounds cannot be overcome by selection. Taken together, our results show that the spontaneous dynamics of genome size naturally prevents it from growing infinitely, even in cases where intuition would suggest an infinite growth. Using quantitative numerical examples, we show that, in practice, a shrinkage bias appears very quickly in genomes undergoing mutation accumulation, even though DNA gains and losses appear to be perfectly symmetrical at first sight. We discuss this spontaneous dynamics in the light of the other evolutionary forces proposed in the literature and argue that it provides them a stability-related size limit below which they can act.     

DNA indels in coding regions reveal selective constraints on protein evolution in the human lineage

Background  

Insertions and deletions of DNA segments (indels) are together with substitutions the major mutational processes that generate genetic variation. Here we focus on recent DNA insertions and deletions in protein coding regions of the human genome to investigate selective constraints on indels in protein evolution.     

A phylogenetic analysis of indel dynamics in the cotton genus

Genome size evolution is a dynamic process involving counterbalancing mechanisms whose actions vary across lineages and over time. Whereas the primary mechanism of expansion, transposable element (TE) amplification, has been widely documented, the evolutionary dynamics of genome contraction have been less thoroughly explored. To evaluate the relative impact and evolutionary stability of the mechanisms that affect genome size, we conducted a phylogenetic analysis of indel rates for 2 genomic regions in 4 Gossypium genomes: the 2 coresident genomes (A(T) and D(T)) of tetraploid cotton and its model diploid progenitors, Gossypium arboreum (A) and Gossypium raimondii (D). We determined the rates of sequence gain or loss along each branch, partitioned by mechanism, and how these changed during species divergence. In general, there has been a propensity toward growth of the diploid genomes and contraction in the polyploid. Most of the size difference between the diploid species occurred prior to polyploid divergence and was largely attributable to TE amplification in the A/A(T) genome. After separating from the true parents of the polyploid genomes, both diploid genomes experienced slower sequence gain than in the ancestor, due to fewer TE insertions in the A genome and a combination of increased deletions and decreased TE insertions in the D genome. Both genomes of the polyploid displayed increased rates of deletion and decreased rates of insertion, leading to a rate of near stasis in D(T) and overall contraction in A(T) resulting in polyploid genome contraction. As expected, TE insertions contributed significantly to the genome size differences; however, intrastrand homologous recombination, although rare, had the most significant impact on the rate of deletion. Small indel data for the diploids suggest the possibility of a bias as the smaller genomes add less or delete more sequence through small indels than do the larger genomes, whereas data for the polyploid suggest increased sequence turnover in general (both as small deletions and small insertions). Illegitimate recombination, although not demonstrated to be a dominant mechanism of genome size change, was biased in the polyploid toward deletions, which may provide a partial explanation of polyploid genomic downsizing.     

A molecular clock based on the expansion of gene families.

There is evidence to suggest that eukaryotic genomes are subject to frequent insertions and deletions of non-coding DNA. This may lead to a gradual increase or decrease in genome size, or to a dynamic equilibrium in which the overall size remains constant. We argue, however, that there is a bias favouring an accumulation of non-coding DNA in the proximity of genes. Such bias causes a progressive change in genome structure regardless of whether the overall genome size increases, decreases or remains constant. We show that this change may serve as a 'molecular clock', supplementing that provided by nucleotide substitution rates.     

Evolution of genome size: multilevel selection, mutation bias or dynamical chaos?

In the past two years, new data on conceptual aspects of the evolution of eukaryotic genome size have appeared, including the adaptivity of genome enlargement, the mechanisms of genome size change and the relation of genome size to organismal complexity. New data on the hypotheses of "selfish DNA" and "mutational equilibrium" have been recently obtained. A relationship is emerging between the intragenomic distribution of noncoding DNA and differential gene expression, which suggests that noncoding DNA is involved in epigenetic organization of the genome and organismal complexity. The standpoint of dynamical chaos, which integrates multilevel selection and mutation biases, may provide a framework for studying the evolution of genome size.     

Deletional bias and the evolution of bacterial genomes

Although bacteria increase their DNA content through horizontal transfer and gene duplication, their genomes remain small and, in particular, lack nonfunctional sequences. This pattern is most readily explained by a pervasive bias towards higher numbers of deletions than insertions. When selection is not strong enough to maintain them, genes are lost in large deletions or inactivated and subsequently eroded. Gene inactivation and loss are particularly apparent in obligate parasites and symbionts, in which dramatic reductions in genome size can result not from selection to lose DNA, but from decreased selection to maintain gene functionality. Here we discuss the evidence showing that deletional bias is a major force that shapes bacterial genomes.     

Breakage-fusion-bridge Cycles and Large Insertions Contribute to the Rapid Evolution of Accessory Chromosomes in a Fungal Pathogen

Chromosomal rearrangements are a major driver of eukaryotic genome evolution, affecting speciation, pathogenicity and cancer progression. Changes in chromosome structure are often initiated by mis-repair of double-strand breaks in the DNA. Mis-repair is particularly likely when telomeres are lost or when dispersed repeats misalign during crossing-over. Fungi carry highly polymorphic chromosomal complements showing substantial variation in chromosome length and number. The mechanisms driving chromosome polymorphism in fungi are poorly understood. We aimed to identify mechanisms of chromosomal rearrangements in the fungal wheat pathogen Zymoseptoria tritici. We combined population genomic resequencing and chromosomal segment PCR assays with electrophoretic karyotyping and resequencing of parents and offspring from experimental crosses to show that this pathogen harbors a highly diverse complement of accessory chromosomes that exhibits strong global geographic differentiation in numbers and lengths of chromosomes. Homologous chromosomes carried highly differentiated gene contents due to numerous insertions and deletions. The largest accessory chromosome recently doubled in length through insertions totaling 380 kb. Based on comparative genomics, we identified the precise breakpoint locations of these insertions. Nondisjunction during meiosis led to chromosome losses in progeny of three different crosses. We showed that a new accessory chromosome emerged in two viable offspring through a fusion between sister chromatids. Such chromosome fusion is likely to initiate a breakage-fusion-bridge (BFB) cycle that can rapidly degenerate chromosomal structure. We suggest that the accessory chromosomes of Z. tritici originated mainly from ancient core chromosomes through a degeneration process that included BFB cycles, nondisjunction and mutational decay of duplicated sequences. The rapidly evolving accessory chromosome complement may serve as a cradle for adaptive evolution in this and other fungal pathogens.     

Species-specific double-strand break repair and genome evolution in plants

Even closely related eukaryotic species may differ drastically in genome size. While insertion of retroelements represents a major source of genome enlargement, the mechanism mediating species- specific deletions is fairly obscure. We analyzed the formation of deletions during double-strand break (DSB) repair in Arabidopsis thaliana and tobacco, two dicotyledonous plant species differing >20-fold in genome size. DSBs were induced by the rare cutting restriction endonuclease I-SCE:I and deletions were identified by loss of function of a negative selectable marker gene containing an I-SCE:I site. Whereas the partial use of micro-homologies in junction formation was similar in both species, in tobacco 40% of the deletions were accompanied by insertions. No insertions could be detected in Arabidopsis , where larger deletions were more frequent, indicating a putative inverse correlation between genome size and the average length of deletions. Such a correlation has been postulated before by a theoretical study on the evolution of related insect genomes and our study now identifies a possible molecular cause for the phenomenon, indicating that species-specific differences in DSB repair might indeed influence genome evolution.     

Plant genome evolution: lessons from comparative genomics at the DNA level

Angiosperm genomes show tremendous variability in genome size and chromosome number. Nevertheless, comparative genetic mapping has revealed genome collinearity of closely related species. Sequence-based comparisons were used to assess the conservation of gene arrangements. Numerous small rearrangements, insertions/deletions, duplications, inversions and translocations have been detected. Importantly, comparative sequence analyses have unambiguously shown micro-collinearity of distantly related plant species. Duplications and subsequent gene loss have been identified as a particular important factor in the evolution of plant genomes.     

Genomic gigantism: DNA loss is slow in mountain grasshoppers

Several studies have shown DNA loss to be inversely correlated with genome size in animals. These studies include a comparison between Drosophila and the cricket, Laupala, but there has been no assessment of DNA loss in insects with very large genomes. Podisma pedestris, the brown mountain grasshopper, has a genome over 100 times as large as that of Drosophila and 10 times as large as that of Laupala. We used 58 paralogous nuclear pseudogenes of mitochondrial origin to study the characteristics of insertion, deletion, and point substitution in P. pedestris and Italopodisma. In animals, these pseudogenes are "dead on arrival"; they are abundant in many different eukaryotes, and their mitochondrial origin simplifies the identification of point substitutions accumulated in nuclear pseudogene lineages. There appears to be a mononucleotide repeat within the 643-bp pseudogene sequence studied that acts as a strong hot spot for insertions or deletions (indels). Because the data for other insect species did not contain such an unusual region, hot spots were excluded from species comparisons. The rate of DNA loss relative to point substitution appears to be considerably and significantly lower in the grasshoppers studied than in Drosophila or Laupala. This suggests that the inverse correlation between genome size and the rate of DNA loss can be extended to comparisons between insects with large or gigantic genomes (i.e., Laupala and Podisma). The low rate of DNA loss implies that in grasshoppers, the accumulation of point mutations is a more potent force for obscuring ancient pseudogenes than their loss through indel accumulation, whereas the reverse is true for Drosophila. The main factor contributing to the difference in the rates of DNA loss estimated for grasshoppers, crickets, and Drosophila appears to be deletion size. Large deletions are relatively rare in Podisma and Italopodisma.     

DEVELOPMENTAL CORRELATES OF GENOME SIZE IN PLETHODONTID SALAMANDERS AND THEIR IMPLICATIONS FOR GENOME EVOLUTION

We present an analysis of the evolutionary relationship between genome size (C-value, mass of DNA per haploid nucleus) and developmental rate using observations of limb regeneration in salamanders of the family Plethodontidae. Rates of growth and differentiation of regenerating limbs are reported for 27 plethodontid species whose C-values range from 14 to 76 picograms. A phylogenetic analysis employing Felsenstein's method of independent contrasts indicates that rate of differentiation is inversely proportional to genome size, although we have not identified any statistically significant association between genome size and the growth rate of regenerating tissue. Our results are consistent with an interpretation that genome size may place a limit on the maximum rate of regeneration attainable in plethodontid salamanders. The implications of our findings for the “junk DNA,” “nucleotypic DNA,” “selfish DNA,” and “skeletal DNA” hypotheses of genome evolution are discussed.     

Short-range template switching in great ape genomes explored using pair hidden Markov models

Many complex genomic rearrangements arise through template switch errors, which occur in DNA replication when there is a transient polymerase switch to an alternate template nearby in three-dimensional space. While typically investigated at kilobase-to-megabase scales, the genomic and evolutionary consequences of this mutational process are not well characterised at smaller scales, where they are often interpreted as clusters of independent substitutions, insertions and deletions. Here we present an improved statistical approach using pair hidden Markov models, and use it to detect and describe short-range template switches underlying clusters of mutations in the multi-way alignment of hominid genomes. Using robust statistics derived from evolutionary genomic simulations, we show that template switch events have been widespread in the evolution of the great apes’ genomes and provide a parsimonious explanation for the presence of many complex mutation clusters in their phylogenetic context. Larger-scale mechanisms of genome rearrangement are typically associated with structural features around breakpoints, and accordingly we show that atypical patterns of secondary structure formation and DNA bending are present at the initial template switch loci. Our methods improve on previous non-probabilistic approaches for computational detection of template switch mutations, allowing the statistical significance of events to be assessed. By specifying realistic evolutionary parameters based on the genomes and taxa involved, our methods can be readily adapted to other intra- or inter-species comparisons.     

Reptiles: a group of transition in the evolution of genome size and of the nucleotypic effect

A comparison between genome size and some phenotypic parameters, such as developmental length and metabolic rate, showed in reptiles a nucleotypic correlation similar to the one observed in birds and mammals. Indeed, like homeotherms, reptiles exhibit a highly significant, inverse correlation of genome size with metabolic rate but unlike amphibians, no relationship with developmental length. Several lines of evidence suggest that these nucleotypic correlations are influenced by body temperature, which also affects the guanine + cytosine nuclear percentage, and that they play an important role in the adaptation of these amniotes. However, the reptilian suborders exhibit differences in the quantitative and compositional characters of the genome that do not completely correspond to differences in the phenotypic parameters commonly involved in the nucleotypic effect. Thus, additional factors could have influenced genome size in this class. These data could be explained with the model of Hartl and Petrov, who observed an inverse correlation between genome size, non-coding portion of the genome and rate of DNA loss and hypothesized a strong role for different spectra of spontaneous insertions and deletions (indels) in the variations of genome size. It is thus reasonable to surmise that variations in the reptilian genome were initially influenced by different indels spectra typical of the diverse lineages, possibly related to different chromosome compartmentalizations. The consequent size increases or decreases would have influenced various morphological and functional cell parameters, and through these some phenotypic characteristics of the whole organism, especially the metabolic rate, very important for environmental adaptation and thus subject to natural selection. Through this "nucleotypic" bond, natural selection would also have controlled genome size variations.     

An evolutionary correlate of genome size change in plethodontid salamanders

Variation in the amount of nuclear DNA, the C-value, does not correlate with differences in morphological complexity. There are two classes of explanations for this observation, which is known as the ''C-value paradox''. The quantity of DNA may serve a ''nucleotypic'' function that is positively selected. Alternatively, large genomes may consist of junk DNA, which increases until it negatively affects fitness. Attempts to resolve the C-value paradox focus on the link between genome size and fitness. This link is usually sought in life history traits, particularly developmental rates. I examined the relationship among two life history traits, egg size and embryonic developmental time and genome size, in 15 species of plethodontid salamanders. Surprisingly, there is no correlation between egg size and developmental time, a relationship included in models of life history evolution. However, genome size is positively correlated with embryonic developmental time, a result that is robust with respect to many sources of variation in the data. Without information on the targets of natural selection it is not possible with these data to distinguish between nucleotypic and junk DNA explanations for the C-value paradox.     

Phylogenetic inference and comparative evolution of a complex microsatellite and its flanking regions in carnivores

We sequenced locus Mel 08, with complex short repetitive motifs, in 24 carnivore species belonging to five different families in order to explore mutational changes in the region in the context of locus and species evolution. This non-coding locus includes up to four different parts or repetitive motifs showing size variability. The variability consists of repeat additions and deletions; substitutions, insertions and/or deletions creating interruptions in the repeat; and substitutions, insertions and deletions in the flanking regions. The locus has different repeat expansions in different carnivore subfamilies. We hypothesize that the complexity of this locus is due to a high mutation rate at an ancestral DNA sequence and, thus, prompts the emergence of repeats at mutational hotspots. High levels of homoplasy were evident, with nine electromorphs representing 28 haplotypes never shared across species. The variability in flanking regions was informative for phylogenetic inference and their evolutionary content. Tree topologies were congruent with relevant hypotheses on current conflicts in carnivore phylogenies, such as: (i) the monophyly of Lutrinae, (ii) the paraphyly of Mustelinae, (iii) the basal position of the Eurasian badger, Meles meles , in the Mustelidae, (iv) the classification of skunks as a separate family, Mephitidae, and (v) the placement of the red panda, Ailurus fulgens , as a monotypic family, Ailuridae, at a basal position in the Musteloidea.     

A Model of Genetic Search for Beneficial Mutations: Estimating the Constructive Capacities of Mutagenesis

We attempted to answer the following question: What evolutionary conditions are required to generate novel genetic modules? Our broad formulation of the problem allows us to simultaneously consider such issues as the relationship between the stage of "genetic search" and the rate of adaptive evolution; the theoretical limits to the generative capacities of spontaneous mutagenesis; and the correlation between genome organization and evolvability. We show that adaptive evolution is feasible only when the mutation rate is fine-tuned to a specific range of values and the structures of the genome and genes are optimized in a certain way. Our quantitative analysis has demonstrated that the rate of evolution of novelty depends on several parameters, such as genome size, the length of a module, the size of the adjacent nonfunctional DNA spacers, and the mutation rate at various genomic scales. We evaluated the efficiency of some mechanisms that increase evolvability: bias in the spectrum of mutation rates towards small mutations, and the availability and size of nonfunctional DNA spacers. We show that the probability of successful duplication and insertion of a copy of a functional module increases by several orders of magnitude depending on the length of the spacers flanking the module. We infer that the adaptive evolution of multicellular organisms has become feasible because of the abundance of nonfunctional DNA spacers, particularly introns, in the genome. We also discuss possible reasons underlying evolutionary retention of the mechanisms that increase evolvability.