Extreme rates and heterogeneity in insect DNA evolution

Summary DNA-DNA hybridization studies of insects, more specificallyDrosophila and cave crickets, have revealed interesting patterns of genome evolution that contrast markedly with what has been seen in other taxa, especially mammals and birds. Insect genomes are composed of sections of single-copy DNA with extreme variation in rates of evolutionary change. This variation is more extreme than between introns and exons; introns fall into the relatively conserved fraction of the genome. Attempts to calculate absolute rates of change inDrosophila DNA have all led to estimates some 5–10 times faster than those found in most vertebrates; this is true even for the more conservative part of the nuclear genome. Finally we point out that morphological similarity, chromosomal similarity, and/or ability to form interspecific hybrids is often associated with quite high levels of single-copy DNA divergence in insects as compared to mammals and birds.     

Rate accelerations in nuclear 18S rDNA of mycoheterotrophic and parasitic angiosperms

Rate variation in genes from all three genomes has been observed frequently in plant lineages with a parasitic and mycoheterotrophic mode of life. While the loss of photosynthetic ability leads to a relaxation of evolutionary constraints in genes involved in the photosynthetic apparatus, it remains to be determined how prevalent increased substitution rates are in nuclear DNA of non-photosynthetic angiosperms. In this study we infer rates of molecular evolution of 18S rDNA of all parasitic and mycoheterotorphic plant families (except Lauraceae and Polygalaceae) using relative rate tests. In several holoparasitic and mycoheterotrophic plant lineages extremely high substitution rates are observed compared to other photosynthetic angiosperms. The position and frequency of these substitutions have been identified to understand the mutation dynamics of 18S rRNA in achlorophyllous plants. Despite the presence of significantly elevated substitution rates, very few mutations occur in major functional and structural regions of the small ribosomal molecule, providing evidence that the efficiency of the translational apparatus in non-photosynthetic plants has not been affected.     

Evolution of modern birds revealed by mitogenomics: timing the radiation and origin of major orders

Mitochondrial (mt) genes and genomes are among the major sources of data for evolutionary studies in birds. This places mitogenomic studies in birds at the core of intense debates in avian evolutionary biology. Indeed, complete mt genomes are actively been used to unveil the phylogenetic relationships among major orders, whereas single genes (e.g., cytochrome c oxidase I [COX1]) are considered standard for species identification and defining species boundaries (DNA barcoding). In this investigation, we study the time of origin and evolutionary relationships among Neoaves orders using complete mt genomes. First, we were able to solve polytomies previously observed at the deep nodes of the Neoaves phylogeny by analyzing 80 mt genomes, including 17 new sequences reported in this investigation. As an example, we found evidence indicating that columbiforms and charadriforms are sister groups. Overall, our analyses indicate that by improving the taxonomic sampling, complete mt genomes can solve the evolutionary relationships among major bird groups. Second, we used our phylogenetic hypotheses to estimate the time of origin of major avian orders as a way to test if their diversification took place prior to the Cretaceous/Tertiary (K/T) boundary. Such timetrees were estimated using several molecular dating approaches and conservative calibration points. Whereas we found time estimates slightly younger than those reported by others, most of the major orders originated prior to the K/T boundary. Finally, we used our timetrees to estimate the rate of evolution of each mt gene. We found great variation on the mutation rates among mt genes and within different bird groups. COX1 was the gene with less variation among Neoaves orders and the one with the least amount of rate heterogeneity across lineages. Such findings support the choice of COX 1 among mt genes as target for developing DNA barcoding approaches in birds.     

Correlation between strand asymmetry and phylogeny in mitochondrial DNA

An evolutionary distance is introduced in order to propose an efficient and feasible procedure for phylogeny studies. Our analysis are based on the strand asymmetry property of mitochondrial DNA, but can be applied to other genomes. Comparison of our results with those reported in conventional phylogenetic trees, gives confidence about our approximation. Our findings support the hypotheses about the origin of the skew and its dependence upon evolutionary pressures, and improves previous efforts on using the strand asymmetry property of genomes for phylogeny inference. For the evolutionary distance introduced here, we observe that the more adequate technique for tree reconstructions correspond to an average link method which employs a sequential clustering algorithm.     

Clustering of Genes Coding for DNA Binding Proteins in a Regionof Atypical Evolution of the Human Genome

Comparison of the human and mouse genomes has revealed that significant variations in evolutionary rates exist among genomic regions and that a large part of this variation is interchromosomal. We confirm in this work, using a large collection of introns, that human chromosome 19 is the one that shows the highest divergence with respect to mouse. To search for other differences among chromosomes, we examine the distribution of gene functions in human and mouse chromosomes using the Gene Ontology definitions. We found by correspondence analysis that among the strongest clusterings of gene functions in human chromosomes is a group of genes coding for DNA binding proteins in chromosome 19. Interestingly, chromosome 19 also has a very high GC content, a feature that has been proposed to promote an opening of the chromatin, thereby facilitating binding of proteins to the DNA helix. In the mouse genome, however, a similar aggregation of genes coding for DNA binding proteins and high GC content cannot be found. This suggests that the distribution of genes coding for DNA binding proteins and the variations of the chromatin accessibility to these proteins are different in the human and mouse genomes. It is likely that the overall high synonymous and intron rates in chromosome 19 are a by-product of the high GC content of this chromosome.Department of Physiology and Molecular Biodiversity, Institut de Biologia Molecular de Barcelona, CSIC, Jordi Girona 18, 08034 Barcelona, Spain     

Coincidence,coevolution, or causation? DNA content,cellsize, and the C‐value enigma

Variation in DNA content has been largely ignored as a factor in evolution, particularly following the advent of sequence-based approaches to genomic analysis. The significant genome size diversity among organisms (more than 200000-fold among eukaryotes) bears no relationship to organismal complexity and both the origins and reasons for the clearly non-random distribution of this variation remain unclear. Several theories have been proposed to explain this 'C-value enigma' (heretofore known as the 'C-value paradox'), each of which can be described as either a mutation pressure' or 'optimal DNA' theory. Mutation pressure theories consider the large portion of non-coding DNA in eukaryotic genomes as either 'junk' or 'selfish' DNA and are important primarily in considerations of the origin of secondary DNA. Optimal DNA theories differ from mutation pressure theories by emphasizing the strong link between DNA content and cell and nuclear volumes. While mutation pressure theories generally explain this association with cell size as coincidental, the nucleoskeletal theory proposes a coevolutionary interaction between nuclear and cell volume, with DNA content adjusted adaptively following shifts in cell size. Each of these approaches to the C-value enigma is problematic for a variety of reasons and the preponderance of the available evidence instead favours the nucleotypic theory which postulates a causal link between bulk DNA amount and cell volume. Under this view, variation in DNA content is under direct selection via its impacts on cellular and organismal parameters. Until now, no satisfactory mechanism has been presented to explain this nucleotypic effect. However, recent advances in the study of cell cycle regulation suggest a possible 'gene nucleus interaction model' which may account for it. The present article provides a detailed review of the debate surrounding the C-value enigma, the various theories proposed to explain it, and the evidence in favour of a causal connection between DNA content and cell size. In addition, a new model of nucleotypic influence is developed, along with suggestions for further empirical investigation. Finally, some evolutionary implications of genome size diversity are considered, and a broadening of the traditional 'biological hierarchy' is recommended.     

Revealing hidden evolutionary capacity to cope with global change

The extent to which global change will impact the long‐term persistence of species depends on their evolutionary potential to adapt to future conditions. While the number of studies that estimate the standing levels of adaptive genetic variation in populations under predicted global change scenarios is growing all the time, few studies have considered multiple environments simultaneously and even fewer have considered evolutionary potential in multivariate context. Because conditions will not be constant, adaptation to climate change is fundamentally a multivariate process so viewing genetic variances and covariances over multivariate space will always be more informative than relying on bivariate genetic correlations between traits. A multivariate approach to understanding the evolutionary capacity to cope with global change is necessary to avoid misestimating adaptive genetic variation in the dimensions in which selection will act. We assessed the evolutionary capacity of the larval stage of the marine polychaete Galeolaria caespitosa to adapt to warmer water temperatures. Galeolaria is an important habitat‐forming species in Australia, and its earlier life‐history stages tend to be more susceptible to stress. We used a powerful quantitative genetics design that assessed the impacts of three temperatures on subsequent survival across over 30 000 embryos across 204 unique families. We found adaptive genetic variation in the two cooler temperatures in our study, but none in the warmest temperature. Based on these results, we would have concluded that this species has very little capacity to evolve to the warmest temperature. However, when we explored genetic variation in multivariate space, we found evidence that larval survival has the potential to evolve even in the warmest temperatures via correlated responses to selection across thermal environments. Future studies should take a multivariate approach to estimating evolutionary capacity to cope with global change lest they misestimate a species’ true adaptive potential.     

Using genetics to understand the dynamics of wild primate populations

While much can be learned about primates by means of observation, the slow life history of many primates means that even decades of dedicated effort cannot illuminate long-term evolutionary processes. For example, while the size of a contemporary population can be estimated from field censuses, it is often desirable to know whether a population has been constant or changing in size over a time frame of hundreds or thousands of years. Even the nature of “a population” is open to question, and the extent to which individuals successfully disperse among defined populations is also difficult to estimate by using observational methods alone. Researchers have thus turned to genetic methods to examine the size, structure, and evolutionary histories of primate populations. Many results have been gained by study of sequence variation of maternally inherited mitochondrial DNA, but in recent years researchers have been increasingly focusing on analysis of short, highly variable microsatellite segments in the autosomal genome for a high-resolution view of evolutionary processes involving both sexes. In this review we describe some of the insights thus gained, and discuss the likely impact on this field of new technologies such as high-throughput DNA sequencing.     

Mutation and selection processes regulating short tandem repeats give rise to genetic and phenotypic diversity across species

Short tandem repeats (STRs) are units of 1–6 bp that repeat in a tandem fashion in DNA. Along with single nucleotide polymorphisms and large structural variations, they are among the major genomic variants underlying genetic, and likely phenotypic, divergence. STRs experience mutation rates that are orders of magnitude higher than other well-studied genotypic variants. Frequent copy number changes result in a wide range of alleles, and provide unique opportunities for modulating complex phenotypes through variation in repeat length. While classical studies have identified key roles of individual STR loci, the advent of improved sequencing technology, high-quality genome assemblies for diverse species, and bioinformatics methods for genome-wide STR analysis now enable more systematic study of STR variation across wide evolutionary ranges. In this review, we explore mutation and selection processes that affect STR copy number evolution, and how these processes give rise to varying STR patterns both within and across species. Finally, we review recent examples of functional and adaptive changes linked to STRs.     

Review and meta‐analysis of natural selection in mitochondrial complex I in metazoans

Variation in mitochondrial DNA is often assumed to be neutral and is used to construct the genealogical relationships among populations and species. However, if extant variation is the result of episodes of positive selection, these genealogies may be incorrect, although this information itself may provide biologically and evolutionary meaningful information. In fact, positive Darwinian selection has been detected in the mitochondrial‐encoded subunits that comprise complex I from diverse taxa with seemingly dissimilar bioenergetic life histories, but the functional implications of the selected sites are unknown. Complex I produces roughly 40% of the proton flux that is used to synthesize ATP from ADP, and a functional model based on the high‐resolution structure of complex I described a unique biomechanical apparatus for proton translocation. We reported positive selection at sites in this apparatus during the evolution of Pacific salmon, and it appeared this was also the case in published reports from other taxa, but a comparison among studies was difficult because different statistical tests were used to detect selection and oftentimes, specific sites were not reported. Here we review the literature of positive selection in mitochondrial genomes, the statistical tests used to detect selection, and the structural and functional models that are currently available to study the physiological implications of selection. We then search for signatures of positive selection among the coding mitochondrial genomes of 237 species with a common set of tests and verify that the ND5 subunit of complex I is a repeated target of positive Darwinian selection in diverse taxa. We propose a novel hypothesis to explain the results based on their bioenergetic life histories and provide a guide for laboratory and field studies to test this hypothesis.     

The evolution of long interspersed repeated DNA (L1, LINE 1) as revealed by the analysis of an ancient rodent L1 DNA family

Summary All modern mammals contain a distinctive, highly repeated (⩾50,000 members) family of long interspersed repeated DNA called the L1 (LINE 1) family. While the modern L1 families were derived from a common ancestor that predated the mammalian radiation ∼80 million years ago, most of the members of these families were generated within the last 5 million years. However, recently we demonstrated that modern murine (Old World rats and mice) genomes share an older long interspersed repeated DNA family that we called Lx. Here we report our analysis of the DNA sequence of Lx family members and the relationship of this family to the modern L1 families in mouse and rat. The extent of DNA sequence divergence between Lx members indicates that the Lx amplification occurred about 12 million years ago, around the time of the murine radiation. Parsimony analysis revealed that Lx elements were ancestral to both the modern rat and mouse L1 families. However, we found that few if any of the evolutionary intermediates between the Lx and the modern L1 families were extensively amplified. Because the modern L1 families have evolved under selective pressure, the evolutionary intermediates must have been capable of replication. Therefore, replicationcompetent L1 elements can reside in genomes without undergoing extensive amplification. We discuss the bearing of our findings on the evolution of L1 DNA elements and the mammalian genome.     

A bird's-eye view of the C-value enigma: genome size, cell size, and metabolic rate in the class aves

For half a century, variation in genome size (C-value) has been an unresolved puzzle in evolutionary biology. While the initial "C-value paradox" was solved with the discovery of noncoding DNA, a much more complex "C-value enigma" remains. The present study focuses on one aspect of this puzzle, namely the small genome sizes of birds. Significant negative correlations are reported between resting metabolic rate and both C-value and erythrocyte size. Cell size is positively correlated with both nucleus size and C-value in birds, as in other vertebrates. These findings shed light on the constraints acting on genome size in birds and illustrate the importance of interactions among various levels of the biological hierarchy, ranging from the subchromosomal to the ecological. Following from a discussion of the mechanistic bases of the correlations reported and the processes by which birds achieved and/or maintain small genomes, a pluralistic approach to the C-value enigma is recommended.     

A meta‐analysis of the strength and nature of cytoplasmic genetic effects

Genetic variation in cytoplasmic genomes (i.e. the mitochondrial genome in animals, and the combined mitochondrial and chloroplast genomes in plants) was traditionally assumed to accumulate under a neutral equilibrium model. This view has, however, come under increasing challenge from studies that have experimentally linked cytoplasmic genetic effects to the expression of life history phenotypes. Such results suggest that genetic variance located within the cytoplasm might be of evolutionary importance and potentially involved in shaping population evolutionary trajectories. As a step towards assessing this assertion, here we conduct a formal meta‐analytic review to quantitatively assess the extent to which cytoplasmic genetic effects contribute to phenotypic expression across animal and plant kingdoms. We report that cytoplasmic effect sizes are generally moderate in size and associated with variation across a range of factors. Specifically, cytoplasmic effects on morphological traits are generally larger than those on life history or metabolic traits. Cytoplasmic effect sizes estimated at the between‐species scale (via interspecies mix‐and‐matching of cytoplasmic and nuclear genomes) are larger than those at the within‐species scale. Furthermore, cytoplasmic effects tied to epistatic interactions with the nuclear genome tend to be stronger than additive cytoplasmic effects, at least when restricting the data set to gonochorous animal species. Our results thus confirm that cytoplasmic genetic variation is commonly tied to phenotypic expression across plants and animals, implicate the cytoplasmic–nuclear interaction as a key unit on which natural selection acts and generally suggest that the genetic variation that lies within the cytoplasm is likely to be entwined in adaptive evolutionary processes.     

Similarity-Based Analysis of Allele Frequency Distribution among Multiple Populations Identifies Adaptive Genomic Structural Variants

Structural variants have a considerable impact on human genomic diversity. However, their evolutionary history remains mostly unexplored. Here, we developed a new method to identify potentially adaptive structural variants based on a similarity-based analysis that incorporates genotype frequency data from 26 populations simultaneously. Using this method, we analyzed 57,629 structural variants and identified 576 structural variants that show unusual population differentiation. Of these putatively adaptive structural variants, we further showed that 24 variants are multiallelic and overlap with coding sequences, and 20 variants are significantly associated with GWAS traits. Closer inspection of the haplotypic variation associated with these putatively adaptive and functional structural variants reveals deviations from neutral expectations due to: 1) population differentiation of rapidly evolving multiallelic variants, 2) incomplete sweeps, and 3) recent population-specific negative selection. Overall, our study provides new methodological insights, documents hundreds of putatively adaptive variants, and introduces evolutionary models that may better explain the complex evolution of structural variants.     

Multiple cellular mechanisms prevent chromosomal rearrangements involving repetitive DNA

Repetitive DNA is present in the eukaryotic genome in the form of segmental duplications, tandem and interspersed repeats, and satellites. Repetitive sequences can be beneficial by serving specific cellular functions (e.g. centromeric and telomeric DNA) and by providing a rapid means for adaptive evolution. However, such elements are also substrates for deleterious chromosomal rearrangements that affect fitness and promote human disease. Recent studies analyzing the role of nuclear organization in DNA repair and factors that suppress non-allelic homologous recombination (NAHR) have provided insights into how genome stability is maintained in eukaryotes. In this review, we outline the types of repetitive sequences seen in eukaryotic genomes and how recombination mechanisms are regulated at the DNA sequence, cell organization, chromatin structure, and cell cycle control levels to prevent chromosomal rearrangements involving these sequences.     

Liberating genetic variance through sex

Genetic variation in fitness is the fundamental prerequisite for adaptive evolutionary change. If there is no variation in survival and reproduction or if this variation has no genetic basis, then the composition of a population will not evolve over time. Consequently, the factors influencing genetic variation in fitness have received close attention from evolutionary biologists. One key factor is the mode of reproduction. Indeed, it has long been thought that sex enhances fitness variation and that this explains the ubiquity of sexual reproduction among eukaryotes. Nevertheless, theoretical studies have demonstrated that sex need not always increase genetic variation in fitness. In particular, if fitness interactions among beneficial alleles (epistasis) are positive, sex can reduce genetic variance in fitness. Empirical data have been sorely needed to settle the issue of whether sex does enhance fitness variation. A recent flurry of studies[1-4] has demonstrated that sex and recombination do dramatically increase genetic variation in fitness and consequently the rate of adaptive evolution. Interpreted in light of evolutionary theory, these studies rule out positive in these experiments epistasis as a major source of genetic associations. Further studies are needed, however, to tease apart other possible sources.     

Comparative analysis of codon usage patterns in chloroplast genomes of the Asteraceae family

Codon usage bias (CUB) is an important evolutionary feature in a genome and has been widely documented from prokaryotes to eukaryotes. However, the significance of CUB in the Asteraceae family has not been well understood, with no Asteraceae species having been analyzed for this characteristic. Here, we use bioinformatics approaches to comparatively analyze the general patterns and influencing factors of CUB in five Asteraceae chloroplast (cp) genomes. The results indicated that the five genomes had similar codon usage patterns, showing a strong bias towards a high representation of NNA and NNT codons. Neutrality analysis showed that these cp genomes had a narrow GC distribution and no significant correlation was observed between GC₁₂ and GC₃. Parity Rule 2 (PR2) plot analysis revealed that purines were used more frequently than pyrimidines. Effective number of codons (ENc)-plot analysis showed that most genes followed the parabolic line of trajectory, but several genes with low ENc values lying below the expected curve were also observed. Furthermore, correspondence analysis of relative synonymous codon usage (RSCU) yielded a first axis that explained only a partial amount of variation of codon usage. These findings suggested that both natural selection and mutational bias contributed to codon bias, while selection was the major force to shape the codon usage in these Asteraceae cp genomes. Our study, which is the first to investigate codon usage patterns in Asteraceae plastomes, will provide helpful information about codon distribution and variation in these species, and also shed light on the genetic and evolutionary mechanisms of codon biology within this family.     

In silico predicted robustness of viroids RNA secondary structures. I. The effect of single mutations

Viroids are plant subviral pathogens whose genomes are constituted by a single-stranded and covalently closed small RNA molecule that does not encode for any protein. Despite this genomic simplicity, they are able of inducing devastating symptoms in susceptible plants. Most of the 29 described viroid species fold into a rodlike or quasi-rodlike structure, whereas a few of them fold as branched structures. The shape of these RNA structures is perhaps one of the most characteristic properties of viroids and sometimes is considered their only phenotype. Here we use RNA thermodynamic secondary structure prediction algorithms to compare the mutational robustness of all viroid species. After characterizing the statistical properties of the distribution of mutational effects on structure stability and the wideness of neutral neighborhood for each viroid species, we show an evolutionary trend toward increased structural robustness during viroid radiation, giving support to the adaptive value of robustness. Differences in robustness among the 2 viroid families can be explained by the larger fragility of branched structures compared with the rodlike ones. We also show that genomic redundancy can contribute to the robustness of these simple RNA genomes.