Comparative population genomics in animals uncovers the determinants of genetic diversity.

Genetic diversity is the amount of variation observed between DNA sequences from distinct individuals of a given species. This pivotal concept of population genetics has implications for species health, domestication, management and conservation. Levels of genetic diversity seem to vary greatly in natural populations and species, but the determinants of this variation, and particularly the relative influences of species biology and ecology versus population history, are still largely mysterious. Here we show that the diversity of a species is predictable, and is determined in the first place by its ecological strategy. We investigated the genome-wide diversity of 76 non-model animal species by sequencing the transcriptome of two to ten individuals in each species. The distribution of genetic diversity between species revealed no detectable influence of geographic range or invasive status but was accurately predicted by key species traits related to parental investment: long-lived or low-fecundity species with brooding ability were genetically less diverse than short-lived or highly fecund ones. Our analysis demonstrates the influence of long-term life-history strategies on species response to short-term environmental perturbations, a result with immediate implications for conservation policies.

Positive selection on sperm ion channels in a brooding brittle star: consequence of life-history traits evolution.

Closely related species are key models to investigate mechanisms leading to reproductive isolation and early stages of diversification, also at the genomic level. The brittle star cryptic species complex Ophioderma longicauda encompasses the sympatric broadcast-spawning species C3 and the internal brooding species C5. Here, we used de novo transcriptome sequencing and assembly in two closely related species displaying contrasting reproductive modes to compare their genetic diversity and to investigate the role of natural selection in reproductive isolation. We reconstructed 20 146 and 22 123 genes for C3 and C5, respectively, and characterized a set of 12 229 orthologs. Genetic diversity was 1.5-2 times higher in C3 compared to C5, confirming that species with low parental investment display higher levels of genetic diversity. Forty-eight genes were the targets of positive diversifying selection during the evolution of the two species. Notably, two genes (NHE and TetraKCNG) are sperm-specific ion channels involved in sperm motility. Ancestral sequence reconstructions show that natural selection targeted the two genes in the brooding species. This may result from an adaptation to the novel environmental conditions surrounding sperm in the brooding species, either directly affecting sperm or via an increase in male/female conflict. This phenomenon could have promoted prezygotic reproductive isolation between C3 and C5. Finally, the sperm receptors to egg chemoattractants differed between C3 and C5 in the ligand-binding region. We propose that mechanisms of species-specific gamete recognition in brittle stars occur during sperm chemotaxis (sperm attraction towards the eggs), contrary to other marine invertebrates where prezygotic barriers to interspecific hybridization typically occur before sperm-egg fusion.

Population genomics of sexual and asexual lineages in fissiparous ribbon worms (Lineus, Nemertea): hybridization, polyploidy and the Meselson effect.

Comparative population genetics in asexual vs. sexual species offers the opportunity to investigate the impact of asexuality on genome evolution. Here, we analyse coding sequence polymorphism and divergence patterns in the fascinating Lineus ribbon worms, a group of marine, carnivorous nemerteans with unusual regeneration abilities, and in which asexual reproduction by fissiparity is documented. The population genomics of the fissiparous L. pseudolacteus is characterized by an extremely high level of heterozygosity and unexpectedly elevated πN /πS ratio, in apparent agreement with theoretical expectations under clonal evolution. Analysis of among-species allele sharing and read-count distribution, however, reveals that L. pseudolacteus is a triploid hybrid between Atlantic populations of L. sanguineus and L. lacteus. We model and quantify the relative impact of hybridity, polyploidy and asexuality on molecular variation patterns in L. pseudolacteus and conclude that (i) the peculiarities of L. pseudolacteus population genomics result in the first place from hybridization and (ii) the accumulation of new mutations through the Meselson effect is more than compensated by processes of heterozygosity erosion, such as gene conversion or gene copy loss. This study illustrates the complexity of the evolutionary processes associated with asexuality and identifies L. pseudolacteus as a promising model to study the first steps of polyploid genome evolution in an asexual context.

Genomic evidence for large, long-lived ancestors to placental mammals.

It is widely assumed that our mammalian ancestors, which lived in the Cretaceous era, were tiny animals that survived massive asteroid impacts in shelters and evolved into modern forms after dinosaurs went extinct, 65 Ma. The small size of most Mesozoic mammalian fossils essentially supports this view. Paleontology, however, is not conclusive regarding the ancestry of extant mammals, because Cretaceous and Paleocene fossils are not easily linked to modern lineages. Here, we use full-genome data to estimate the longevity and body mass of early placental mammals. Analyzing 36 fully sequenced mammalian genomes, we reconstruct two aspects of the ancestral genome dynamics, namely GC-content evolution and nonsynonymous over synonymous rate ratio. Linking these molecular evolutionary processes to life-history traits in modern species, we estimate that early placental mammals had a life span above 25 years and a body mass above 1 kg. This is similar to current primates, cetartiodactyls, or carnivores, but markedly different from mice or shrews, challenging the dominant view about mammalian origin and evolution. Our results imply that long-lived mammals existed in the Cretaceous era and were the most successful in evolution, opening new perspectives about the conditions for survival to the Cretaceous-Tertiary crisis.

Mitochondrial DNA as a tool for reconstructing past life-history traits in mammals.

Reconstructing the ancestral characteristics of species is a major goal in evolutionary and comparative biology. Unfortunately, fossils are not always available and sufficiently informative, and phylogenetic methods based on models of character evolution can be unsatisfactory. Genomic data offer a new opportunity to estimate ancestral character states, through (i) the correlation between DNA evolutionary processes and species life-history traits and (ii) available reliable methods for ancestral sequence inference. Here, we assess the relevance of mitochondrial DNA--the most popular molecular marker in animals--as a predictor of ancestral life-history traits in mammals, using the order of Cetartiodactyla as a benchmark. Using the complete set of 13 mitochondrial protein-coding genes, we show that the lineage-specific nonsynonymous over synonymous substitution rate ratio (dN/dS) is closely correlated with the species body mass, longevity and age of sexual maturity in Cetartiodactyla and can be used as a marker of ancestral traits provided that the noise introduced by short branches is appropriately dealt with. Based on ancestral dN/dS estimates, we predict that the first cetartiodactyls were relatively small animals (around 20 kg). This finding is in accordance with Cope's rule and the fossil record but could not be recovered via continuous character evolution methods.

Population genomics of eusocial insects: the costs of a vertebrate-like effective population size.

The evolution of reproductive division of labour and social life in social insects has lead to the emergence of several life-history traits and adaptations typical of larger organisms: social insect colonies can reach masses of several kilograms, they start reproducing only when they are several years old, and can live for decades. These features and the monopolization of reproduction by only one or few individuals in a colony should affect molecular evolution by reducing the effective population size. We tested this prediction by analysing genome-wide patterns of coding sequence polymorphism and divergence in eusocial vs. noneusocial insects based on newly generated RNA-seq data. We report very low amounts of genetic polymorphism and an elevated ratio of nonsynonymous to synonymous changes ­ a marker of the effective population size ­ in four distinct species of eusocial insects, which were more similar to vertebrates than to solitary insects regarding molecular evolutionary processes. Moreover, the ratio of nonsynonymous to synonymous substitutions was positively correlated with the level of social complexity across ant species. These results are fully consistent with the hypothesis of a reduced effective population size and an increased genetic load in eusocial insects, indicating that the evolution of social life has important consequences at both the genomic and population levels.

The determinants of the molecular substitution process in turtles.

Neutral rates of molecular evolution vary across species, and this variation has been shown to be related to biological traits. One of the first patterns to be observed in vertebrates has been an inverse relationship between body mass (BM) and substitution rates. The effects of three major life-history traits (LHT) that covary with BM - metabolic rate, generation time and longevity (LON) - have been invoked to explain this relationship. However, most of the theoretical and empirical evidence supporting this relationship comes from endothermic vertebrates, that is, mammals and birds, in which the environmental conditions, especially temperature, do not have a direct impact on cellular and molecular biology. We analysed the variations in mitochondrial and nuclear rates of synonymous substitution across 224 turtle species and examined their correlation with two LHT (LON and BM) and two environmental variables [latitude (LAT) and habitat]. Our analyses indicate that in turtles, neutral rates of molecular evolution are hardly correlated with LON or BM. Rather, both the mitochondrial and nuclear substitution rates are significantly correlated with LAT - faster evolution in the tropics - and especially so for aquatic species. These results question the generality of the relationships reported in mammals and birds and suggest that environmental factors might be the strongest determinants of the mutation rate in ectotherms.

Mitochondrial DNA as a marker of molecular diversity: a reappraisal.

Over the last three decades, mitochondrial DNA has been the most popular marker of molecular diversity, for a combination of technical ease-of-use considerations, and supposed biological and evolutionary properties of clonality, near-neutrality and clock-like nature of its substitution rate. Reviewing recent literature on the subject, we argue that mitochondrial DNA is not always clonal, far from neutrally evolving and certainly not clock-like, questioning its relevance as a witness of recent species and population history. We critically evaluate the usage of mitochondrial DNA for species delineation and identification. Finally, we note the great potential of accumulating mtDNA data for evolutionary and functional analysis of the mitochondrial genome.

A nonhyperthermophilic common ancestor to extant life forms.

The G+C nucleotide content of ribosomal RNA (rRNA) sequences is strongly correlated with the optimal growth temperature of prokaryotes. This property allows inference of the environmental temperature of the common ancestor to all life forms from knowledge of the G+C content of its rRNA sequences. A model of sequence evolution, assuming varying G+C content among lineages and unequal substitution rates among sites, was devised to estimate ancestral base compositions. This method was applied to rRNA sequences of various species representing the major lineages of life. The inferred G+C content of the common ancestor to extant life forms appears incompatible with survival at high temperature. This finding challenges a widely accepted hypothesis about the origin of life.

Electrostatic analysis of effects of ions on the inhibition of corn root plasma membrane Mg2+-ATPase by the bivalent orthovanadate.

The microelectrophoretic mobility of corn root plasma membranes and the inhibition of the Mg+2-ATPase by vanadate were investigated under different ionic conditions. The Mg2+-ATPase was uncompetitively inhibited and a 10-fold variation of the apparent inhibition constant was observed, depending on the addition of K+ and Mg2+. The determination of the zeta potential indicated that a 5-fold decrease of the apparent inhibition constant was due to aspecific electrostatic interactions of the vanadate anion and the negative charge of the membrane. The screening and masking effects of 6 mM free Mg2+ totally abolished electrostatic interactions and allowed the direct determination of the intrinsic vanadate inhibition constant (KIi). On the other hand, a specific, non-electrostatic, effect of K+ caused a 2-fold decrease of the inhibition constant in addition to the electrostatic effect. Finally, the electrostatic analysis indicates that the Mg2+-ATPase is inhibited by the monomeric bivalent anion HVO4(2-).

Isochore evolution in mammals: a human-like ancestral structure.

Codon usage in mammals is mainly determined by the spatial arrangement of genomic G + C-content, i.e., the isochore structure. Ancestral G + C-content at third codon positions of 27 nuclear protein-coding genes of eutherian mammals was estimated by maximum-likelihood analysis on the basis of a nonhomogeneous DNA substitution model, accounting for variable base compositions among present-day sequences. Data consistently supported a human-like ancestral pattern, i.e., highly variable G + C-content among genes. The mouse genomic structure-more narrow G + C-content distribution-would be a derived state. The circumstances of isochore evolution are discussed with respect to this result. A possible relationship between G + C-content homogenization in murid genomes and high mutation rate is proposed, consistent with the negative selection hypothesis for isochore maintenance in mammals.

Detecting bottlenecks and selective sweeps from DNA sequence polymorphism.

A coalescence-based maximum-likelihood method is presented that aims to (i) detect diversity-reducing events in the recent history of a population and (ii) distinguish between demographic (e.g., bottlenecks) and selective causes (selective sweep) of a recent reduction of genetic variability. The former goal is achieved by taking account of the distortion in the shape of gene genealogies generated by diversity-reducing events: gene trees tend to be more star-like than under the standard coalescent. The latter issue is addressed by comparing patterns between loci: demographic events apply to the whole genome whereas selective events affect distinct regions of the genome to a varying extent. The maximum-likelihood approach allows one to estimate the time and strength of diversity-reducing events and to choose among competing hypotheses. An application to sequence data from an African population of Drosophila melanogaster shows that the bottleneck hypothesis is unlikely and that one or several selective sweeps probably occurred in the recent history of this population.

Effects of Elevated Sucrose-Phosphate Synthase Activity on Photosynthesis, Assimilate Partitioning, and Growth in Tomato (Lycopersicon esculentum var UC82B).

The expression of a sucrose-phosphate synthase (SPS) gene from maize (Zea mays, a monocotyledon) in tomato (Lycopersicon esculentum, a dicotyledon) resulted in marked increases in extractable SPS activity in the light and the dark. Diurnal modulation of the native tomato SPS activity was found. However, when the maize enzyme was present the tomato leaf cells were unable to regulate its activation state. No detrimental effects were observed and total dry matter production was unchanged. However, carbon allocation within the plants was modified such that in shoots it increased, whereas in roots it decreased. There was, therefore, a change in the shoot:root dry weight ratio favoring the shoot. This was positively correlated with increased SPS activity in leaves. SPS was a major determinant of the amount of starch in leaves as well as sucrose. There was a strong positive correlation between the ratio of sucrose to starch and SPS activity in leaves. Therefore, SPS activity is a major determinant of the partitioning of photosynthetically fixed carbon in the leaf and in the whole plant. The photosynthetic rate in air was not significantly increased as a result of elevated leaf SPS activity. However, the light- and CO2-saturated rate of photosynthesis was increased by about 20% in leaves expressing high SPS. In addition, the temporary enhancement of the photosynthetic rate following brief exposures to low light was increased in the high SPS plants relative to controls. We conclude that the level of SPS in the leaves plays a pivotal role in carbon partitioning. Furthermore, high SPS levels have the potential to boost photosynthetic rates under favorable conditions.

Markov-modulated Markov chains and the covarion process of molecular evolution.

The covarion (or site specific rate variation, SSRV) process of biological sequence evolution is a process by which the evolutionary rate of a nucleotide/amino acid/codon position can change in time. In this paper, we introduce time-continuous, space-discrete, Markov-modulated Markov chains as a model for representing SSRV processes, generalizing existing theory to any model of rate change. We propose a fast algorithm for diagonalizing the generator matrix of relevant Markov-modulated Markov processes. This algorithm makes phylogeny likelihood calculation tractable even for a large number of rate classes and a large number of states, so that SSRV models become applicable to amino acid or codon sequence datasets. Using this algorithm, we investigate the accuracy of the discrete approximation to the Gamma distribution of evolutionary rates, widely used in molecular phylogeny. We show that a relatively large number of classes is required to achieve accurate approximation of the exact likelihood when the number of analyzed sequences exceeds 20, both under the SSRV and among site rate variation (ASRV) models.

Mouse biodiversity in the genomic era.

Comparative genomics has developed by comparison of distantly related genomes, for which the link between the reported evolutionary changes and species development/physiology/ecology is not obvious. It is argued that the mouse (genus Mus) is an optimal model for microevolutionary genomics in vertebrates. This is because the mouse genome sequence, physical and genetic map have been completed, because mouse genetics, morpho-anatomy, pathology, behavior and ecology are well-studied, and because the Mus genus is a diverse, well- documented taxon, allowing comparative studies at the level of individual, population, subspecies, and species. The potential of the interaction between mouse genome and mouse biodiversity is illustrated by recent studies of speciation in the house mouse Mus musculus, and studies about the evolution of isochores, the peculiar pattern of GC-content variation across mammalian genomes.

Molecular phylogeny of Eubacteria: a new multiple tree analysis method applied to 15 sequence data sets questions the monophyly of gram-positive bacteria.

Phylogenetic relationships between the major eubacterial phyla were studied using the sequences of 15 homologous bacterial genes. Neither the classical concatenation strategy nor a new multiple tree analysis method involving statistical tests of the inferred phylogenetic relationships provided any (solid) conclusions about eubacterial phylogeny; no pairs of eubacterial phyla proved to be closer to each other in the 15 reconstructed trees than would be expected for trees with random topologies. The phylogeny of Eubacteria therefore appears to be tightly bush-like. Moreover, results from both concatenation and multiple tree analysis raise doubts concerning the monophyly of the so-called Gram-positive bacteria phylum, since the monophyly hypothesis is no more strongly supported by data than its alternatives. It is noteworthy that the structural bases for the Gram-positive phenotype are not incompatible with the hypothesis of independent emergence of this character at two different times.

A test of the hypothesis of an autopolyploid vs. allopolyploid origin for a tetraploid lineage: application to the genus barbus (Cyprinidae)

A new method is described for determination of the origin of polyploid lineages. It tests the hypothesis that a tetraploid lineage originated via autopolyploidization vs. allopolyploidization. The method is based on the hypothesis that, in the case of autopolyploidy, any genetic marker in the first tetraploid ancestor is represented by two copies (one for each homoeologous chromosome of the haploid complement), whereas in allopolyploidy some markers absent from one of the hybridizing species will display one copy at most. The model requires knowledge of the phylogeny (topology and branch lengths) of a sample of species descending from the same tetraploidization event, together with the number of homoeologous copies present in each species for a set of neutral markers. The likelihood of a given proportion of the markers being present in both homoeologous chromosome pairs of the ancestral tetraploid is expressed as a function of the deletion rate of a marker. In the case of an autopolyploid origin, this proportion equals one. A likelihood-ratio test was carried out to test this hypothesis. The method was used to examine five microsatellite loci in eight species of Barbus (sensu lato). Assuming the validity of the hypotheses on phylogenetic relationships and evolutionary rates, the test rejects the possibility that European tetraploid barbs originated through autopolyploidy. This is the first test that can reject autopolyploidy, and it would appear particularly useful for phylogenetic studies in taxa where hybridization is known and where, consequently, undetected reticulate evolution may impair phylogenetic reconstruction.

Reference-free transcriptome assembly in non-model animals from next-generation sequencing data.

Next-generation sequencing (NGS) technologies offer the opportunity for population genomic study of non-model organisms sampled in the wild. The transcriptome is a convenient and popular target for such purposes. However, designing genetic markers from NGS transcriptome data requires assembling gene-coding sequences out of short reads. This is a complex task owing to gene duplications, genetic polymorphism, alternative splicing and transcription noise. Typical assembling programmes return thousands of predicted contigs, whose connection to the species true gene content is unclear, and from which SNP definition is uneasy. Here, the transcriptomes of five diverse non-model animal species (hare, turtle, ant, oyster and tunicate) were assembled from newly generated 454 and Illumina sequence reads. In two species for which a reference genome is available, a new procedure was introduced to annotate each predicted contig as either a full-length cDNA, fragment, chimera, allele, paralogue, genomic sequence or other, based on the number of, and overlap between, blast hits to the appropriate reference. Analyses showed that (i) the highest quality assemblies are obtained when 454 and Illumina data are combined, (ii) typical de novo assemblies include a majority of irrelevant cDNA predictions and (iii) assemblies can be appropriately cleaned by filtering contigs based on length and coverage. We conclude that robust, reference-free assembly of thousands of genes from transcriptomic NGS data is possible, opening promising perspectives for transcriptome-based population genomics in animals. A Galaxy pipeline implementing our best-performing assembling strategy is provided.

Coevolution within and between genes.

Interacting biological systems do not evolve independently, as exemplified many times at the cellular, organismal and ecosystem levels. Biological molecules interact tightly, and should therefore coevolve as well. Here we review the literature about molecular coevolution, between residues within RNAs or proteins, and between proteins. A panel of methodological and bioinformatic approaches have been developed to address this issue, yielding contrasting results: a strong coevolutionary signal is detected in RNA stems, whereas proteins show only moderate, uneasy to interpret departure from the independence hypothesis. The reasons for this discrepancy are discussed.