Genome-wide investigation reveals high evolutionary rates in annual model plants

Genome-scale phylogeny and comparative genomics of the fungal order Sordariales

The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.

Detecting signals of adaptive evolution in grape plastomes with a focus on the Cretaceous-Palaeogene (K/Pg) transition

BACKGROUND AND AIMS

Although plastid genes are widely used in phylogenetic studies, signals of positive selection have been scarcely investigated in the grape family. The plastomes from 91 accessions of Vitaceae were examined to understand the extent to which positive selection is present and to identify which genes are involved. Moreover, the changes through time of genes under episodic positive selection were investigated and the hypothesis of an adaptive process following the Cretaceous-Palaeogene (K/Pg) transition about 66 million years ago was tested.

METHODS

Different codon-substitution models were used to assess pervasive and episodic positive selection events on 70 candidate plastid genes. Divergence times between lineages were estimated and stochastic character mapping analysis was used to simulate variation over time of the genes found to be under episodic positive selection.

KEY RESULTS

A total of 20 plastid genes (29 %) showed positive selection. Among them, 14 genes showed pervasive signatures of positive selection and nine genes showed episodic signatures of positive selection. In particular, four of the nine genes (psbK, rpl20, rpoB, rps11) exhibited a similar pattern showing an increase in the rate of variation close to the K/Pg transition.

CONCLUSION

Multiple analyses have shown that the grape family has experienced ancient and recent positive selection events and that the targeted genes are involved in essential functions such as photosynthesis, self-replication and metabolism. Our results are consistent with the idea that the K/Pg transition has favoured an increased rate of change in some genes. Intense environmental perturbations have influenced the rapid diversification of certain lineages, and new mutations arising on some plastid genes may have been fixed by natural selection over the course of many generations.

A phylogenetic analysis of Bromus (Poaceae: Pooideae: Bromeae) based on nuclear ribosomal and plastid data, with a focus on Bromus sect. Bromus

To investigate phylogenetic relationships among and within major lineages of Bromus, with focus on Bromus sect. Bromus, we analyzed DNA sequences from two nuclear ribosomal (ITS, ETS) and two plastid (rpl32-trnL_UAG , matK) regions. We sampled 103 ingroup accessions representing 26 taxa of B. section Bromus and 15 species of other Bromus sections. Our analyses confirm the monophyly of Bromus s.l. and identify incongruence between nuclear ribosomal and plastid data partitions for relationships within and among major Bromus lineages. Results support classification of B. pumilio and B. gracillimus within B. sect. Boissiera and B. sect. Nevskiella, respectively. These species are sister groups and are closely related to B. densus (B. sect. Mexibromus) in nrDNA trees and Bromus sect. Ceratochloa in plastid trees. Bromus sect. Bromopsis is paraphyletic. In nrDNA trees, species of Bromus sects. Bromopsis, Ceratochloa, Neobromus, and Genea plus B. rechingeri of B. sect. Bromus form a clade, in which B. tomentellus is sister to a B. sect. Genea-B. rechingeri clade. In the plastid trees, by contrast, B. sect. Bromopsis species except B. tomentosus form a clade, and B. tomentosus is sister to a clade comprising B. sect. Bromus and B. sect. Genea species. Affinities of B. gedrosianus, B. pulchellus, and B. rechingeri (members of the B. pectinatus complex), as well as B. oxyodon and B. sewerzowii, are discordant between nrDNA and plastid trees. We infer these species may have obtained their plastomes via chloroplast capture from species of B. sect. Bromus and B. sect. Genea. Within B. sect. Bromus, B. alopecuros subsp. caroli-henrici, a clade comprising B. hordeaceus and B. interruptus, and B. scoparius are successive sister groups to the rest of the section in the nrDNA phylogeny. Most relationships among the remaining species of B. sect. Bromus are unresolved in the nrDNA and plastid trees. Given these results, we infer that most B. sect. Bromus species likely diversified relatively recently. None of the subdivisional taxa proposed for Bromus sect. Bromus over the last century correspond to natural groups identified in our phylogenetic analyses except for a group including B. hordeaceus and B. interruptus.

Fitness effects of somatic mutations accumulating during vegetative growth

The unique life form of plants promotes the accumulation of somatic mutations that can be passed to offspring in the next generation, because the same meristem cells responsible for vegetative growth also generate gametes for sexual reproduction. However, little is known about the consequences of somatic mutation accumulation for offspring fitness. We evaluate the fitness effects of somatic mutations in Mimulus guttatus by comparing progeny from self-pollinations made within the same flower (autogamy) to progeny from self-pollinations made between stems on the same plant (geitonogamy). The effects of somatic mutations are evident from this comparison, as autogamy leads to homozygosity of a proportion of somatic mutations, but progeny from geitonogamy remain heterozygous for mutations unique to each stem. In two different experiments, we find consistent fitness effects of somatic mutations from individual stems. Surprisingly, several progeny groups from autogamous crosses displayed increases in fitness compared to progeny from geitonogamy crosses, likely indicating that beneficial somatic mutations occurred in some stems. These results support the hypothesis that somatic mutations accumulate during vegetative growth, but they are filtered by different forms of selection that occur throughout development, resulting in the culling of expressed deleterious mutations and the retention of beneficial mutations.

Fire as a key driver of Earth's biodiversity

Many terrestrial ecosystems are fire prone, such that their composition and structure are largely due to their fire regime. Regions subject to regular fire have exceptionally high levels of species richness and endemism, and fire has been proposed as a major driver of their diversity, within the context of climate, resource availability and environmental heterogeneity. However, current fire-management practices rarely take into account the ecological and evolutionary roles of fire in maintaining biodiversity. Here, we focus on the mechanisms that enable fire to act as a major ecological and evolutionary force that promotes and maintains biodiversity over numerous spatiotemporal scales. From an ecological perspective, the vegetation, topography and local weather conditions during a fire generate a landscape with spatial and temporal variation in fire-related patches (pyrodiversity), and these produce the biotic and environmental heterogeneity that drives biodiversity across local and regional scales. There have been few empirical tests of the proposition that 'pyrodiversity begets biodiversity' but we show that biodiversity should peak at moderately high levels of pyrodiversity. Overall species richness is greatest immediately after fire and declines monotonically over time, with postfire successional pathways dictated by animal habitat preferences and varying lifespans among resident plants. Theory and data support the 'intermediate disturbance hypothesis' when mean patch species diversity is correlated with mean fire intervals. Postfire persistence, recruitment and immigration allow species with different life histories to coexist. From an evolutionary perspective, fire drives population turnover and diversification by promoting a wide range of adaptive responses to particular fire regimes. Among 39 comparisons, the number of species in 26 fire-prone lineages is much higher than that in their non-fire-prone sister lineages. Fire and its byproducts may have direct mutagenic effects, producing novel genotypes that can lead to trait innovation and even speciation. A paradigm shift aimed at restoring biodiversity-maintaining fire regimes across broad landscapes is required among the fire research and management communities. This will require ecologists and other professionals to spread the burgeoning fire-science knowledge beyond scientific publications to the broader public, politicians and media.

Codon Usage Differences among Genes Expressed in Different Tissues of Drosophila melanogaster

Codon usage patterns are affected by both mutational biases and translational selection. The frequency at which each codon is used in the genome is directly linked to the cellular concentrations of their corresponding tRNAs. Transfer RNA abundances—as well as the abundances of other potentially relevant factors, such as RNA-binding proteins—may vary across different tissues, making it possible that genes expressed in different tissues are subject to different translational selection regimes, and thus differ in their patterns of codon usage. These differences, however, are poorly understood, having been studied only in Arabidopsis, rice and human, with controversial results in human. Drosophila melanogaster is a suitable model organism to study tissue-specific codon adaptation given its large effective population size. Here, we compare 2,046 genes, each expressed specifically in one tissue of D. melanogaster. We show that genes expressed in different tissues exhibit significant differences in their patterns of codon usage, and that these differences are only partially due to differences in GC content, expression levels, or protein lengths. Remarkably, these differences are stronger when analyses are restricted to highly expressed genes. Our results strongly suggest that genes expressed in different tissues are subject to different regimes of translational selection.

Rapid cell division of Staphylococcus aureus during colonization of the human nose

BACKGROUND

Staphylococcus aureus is an important opportunistic pathogen and a commensal bacterium, thriving in the nasal cavities of 20% of the human population. Little is known about the dynamics of asymptomatic colonization and the occasional transition to infectious disease.

RESULTS

In this study, we inferred that S. aureus cells replicate every one to three hours on average while colonizing the human nose, based on two independent lines of genomic evidence. First, we collected nasal swab samples from human subjects, extracted and sequenced metagenomic DNA, and analyzed the distribution of sequencing coverage along the staphylococcal chromosome. Calibration of this data by comparison to a laboratory culture enabled measuring S. aureus cell division rates in nasal samples. Second, we applied mutation accumulation experiments paired with genome sequencing to measure spontaneous mutation rates at a genome scale. Relating these mutation rates to annual evolutionary rates confirmed that nasal S. aureus continuously pass several thousand cell divisions per year when averaged over large, globally distributed populations and over many years, corresponding to generation times of less than two hours.

CONCLUSIONS

The cell division rates we determined were higher than the fastest documented rates during fulminant disease progression (in a mouse model of systemic infection) and much higher than those previously measured in expectorated sputum from cystic fibrosis patients. This paper supplies absolute in-vivo generation times for an important bacterial commensal, indicating that colonization of the human upper respiratory tract is characterized by a highly dynamic equilibrium between bacterial growth and removal.

Different divergence events for three pairs of PEBPs in Gossypium as implied by evolutionary analysis

INTRODUCTION

The phosphatidylethanolamine-binding protein (PEBP) gene family plays a crucial role in seed germination, reproductive transformation, and other important developmental processes in plants, but its distribution in Gossypium genomes or species, evolutionary properties, and the fates of multiple duplicated genes remain unclear.

OBJECTIVES

The primary objectives of this study were to elucidate the distribution and characteristics of PEBP genes in Gossypium, as well as the evolutionary pattern of duplication and deletion, and functional differentiation of PEBPs in plants.

METHODS

Using the PEBP protein sequences in Arabidopsis thaliana as queries, blast alignment was carried out for the identification of PEBP genes in four sequenced cotton species. Using the primers designed according to the PEBP genome sequences, PEBP genes were cloned from 15 representative genomes of Gossypium genus, and the gene structure, CDS sequence, protein sequence and properties were predicted and phylogenetic analysis was performed. Taking PEBP proteins of grape as reference, grouping of orthologous gene, analysis of phylogeny and divergence of PEBPs in nine species were conducted to reconstruct the evolutionary pattern of PEBP genes in plants.

RESULTS

We identified and cloned 160 PEBPs from 15 cotton species, and the phylogenetic analysis showed that the genes could be classified into the following three subfamilies: MFT-like, FT-like and TFL1-like. There were eight single orthologous group (OG) members in each diploid and 16 double OG members in each tetraploid. An analysis of the expression and selective pressure indicated that expression divergence and strong purification selection within the same OG presented in the PEBP gene family.

CONCLUSION

An evolutionary pattern of duplication and deletion of the PEBP family in the evolutionary history of Gossypium was suggested, and three pairs of genes resulted from different whole-genome duplication events.

Mutation rate analysis via parent-progeny sequencing of the perennial peach. I. A low rate in woody perennials and a higher mutagenicity in hybrids

Mutation rates vary between species, between strains within species and between regions within a genome. What are the determinants of these forms of variation? Here, via parent–offspring sequencing of the peach we ask whether (i) woody perennials tend to have lower per unit time mutation rates compared to annuals, and (ii) hybrid strains have high mutation rates. Between a leaf from a low heterozygosity individual, derived from an intraspecific cross, to a leaf of its selfed progeny, the mutation rate is 7.77 × 10⁻⁹ point mutations per bp per generation, similar to Arabidopsis thaliana (7.0–7.4 × 10⁻⁹ point mutations per bp per generation). This suggests a low per unit time mutation rate as the generation time is much longer in peach. This is supported by our estimate of 9.48 × 10⁻⁹ point mutations per bp per generation from a 200-year-old low heterozygosity peach to its progeny. From a more highly heterozygous individual derived from an interspecific cross to its selfed progeny, the mutation rate is 1.38 × 10⁻⁸ mutations per site per generation, consistent with raised rates in hybrids. Our data thus suggest that (i) peach has an approximately order of magnitude lower mutation rate per unit time than Arabidopsis, consistent with reports of low evolutionary rates in woody perennials, and (ii) hybridization may, indeed, be associated with increased mutation rates as considered over a century ago.

Mutation rate analysis via parent-progeny sequencing of the perennial peach. II. No evidence for recombination-associated mutation

Mutation rates and recombination rates vary between species and between regions within a genome. What are the determinants of these forms of variation? Prior evidence has suggested that the recombination might be mutagenic with an excess of new mutations in the vicinity of recombination break points. As it is conjectured that domesticated taxa have higher recombination rates than wild ones, we expect domesticated taxa to have raised mutation rates. Here, we use parent–offspring sequencing in domesticated and wild peach to ask (i) whether recombination is mutagenic, and (ii) whether domesticated peach has a higher recombination rate than wild peach. We find no evidence that domesticated peach has an increased recombination rate, nor an increased mutation rate near recombination events. If recombination is mutagenic in this taxa, the effect is too weak to be detected by our analysis. While an absence of recombination-associated mutation might explain an absence of a recombination–heterozygozity correlation in peach, we caution against such an interpretation.

Multispeed genome diploidization and diversification after an ancient allopolyploidization

Hybridization and genome doubling (allopolyploidy) have led to evolutionary novelties as well as to the origin of new clades and species. Despite the importance of allopolyploidization, the dynamics of postpolyploid diploidization (PPD) at the genome level has been only sparsely studied. The Microlepidieae (MICR) is a crucifer tribe of 17 genera and c. 56 species endemic to Australia and New Zealand. Our phylogenetic and cytogenomic analyses revealed that MICR originated via an intertribal hybridization between ancestors of Crucihimalayeae (n = 8; maternal genome) and Smelowskieae (n = 7; paternal genome), both native to the Northern Hemisphere. The reconstructed ancestral allopolyploid genome (n = 15) originated probably in northeastern Asia or western North America during the Late Miocene (c. 10.6-7 million years ago) and reached the Australian mainland via long-distance dispersal. In Australia, the allotetraploid genome diverged into at least three main subclades exhibiting different levels of PPD and diversity: 1.25-fold descending dysploidy (DD) of n = 15 → n = 12 (autopolyploidy → 24) in perennial Arabidella (3 species), 1.5-fold DD of n = 15 → n = 10 in the perennial Pachycladon (11 spp.) and 2.1-3.75-fold DD of n = 15 → n = 7-4 in the largely annual crown-group genera (42 spp. in 15 genera). These results are among the first to demonstrate multispeed genome evolution in taxa descending from a common allopolyploid ancestor. It is suggested that clade-specific PPD can operate at different rates and efficacies and can be tentatively linked to life histories and the extent of taxonomic diversity.

Comprehensive analysis of correlations among codon usage bias, gene expression, and substitution rate in Arachis duranensis and Arachis ipaënsis orthologs

The relationship between evolutionary rates and gene expression in model plant orthologs is well documented. However, little is known about the relationships between gene expression and evolutionary trends in Arachis orthologs. We identified 7,435 one-to-one orthologs, including 925 single-copy and 6,510 multiple-copy sequences in Arachis duranensis and Arachis ipaënsis. Codon usage was stronger for shorter polypeptides, which were encoded by codons with higher GC contents. Highly expressed coding sequences had higher codon usage bias, GC content, and expression breadth. Additionally, expression breadth was positively correlated with polypeptide length, but there was no correlation between gene expression and polypeptide length. Inferred selective pressure was also negatively correlated with both gene expression and expression breadth in all one-to-one orthologs, while positively but non-significantly correlated with gene expression in sequences with signatures of positive selection. Gene expression levels and expression breadth were significantly higher for single-copy genes than for multiple-copy genes. Similarly, the gene expression and expression breadth in sequences with signatures of purifying selection were higher than those of sequences with positive selective signatures. These results indicated that gene expression differed between single-copy and multiple-copy genes as well as sequences with signatures of positive and purifying selection.

Evolutionary Analyses of GRAS Transcription Factors in Angiosperms

GRAS transcription factors (TFs) play critical roles in plant growth and development such as gibberellin and mycorrhizal signaling. Proteins belonging to this gene family contain a typical GRAS domain in the C-terminal sequence, whereas the N-terminal region is highly variable. Although, GRAS genes have been characterized in a number of plant species, their classification is still not completely resolved. Based on a panel of eight representative species of angiosperms, we identified 29 orthologous groups or orthogroups (OGs) for the GRAS gene family, suggesting that at least 29 ancestor genes were present in the angiosperm lineage before the "Amborella" evolutionary split. Interestingly, some taxonomic groups were missing members of one or more OGs. The gene number expansion usually observed in transcription factors was not observed in GRAS while the genome triplication ancestral to the eudicots (γ hexaploidization event) was detectable in a limited number of GRAS orthogroups. We also found conserved OG-specific motifs in the variable N-terminal region. Finally, we could regroup OGs in 17 subfamilies for which names were homogenized based on a literature review and described 5 new subfamilies (DLT, RAD1, RAM1, SCLA, and SCLB). This study establishes a consistent framework for the classification of GRAS members in angiosperm species, and thereby a tool to correctly establish the orthologous relationships of GRAS genes in most of the food crops in order to facilitate any subsequent functional analyses in the GRAS gene family. The multi-fasta file containing all the sequences used in our study could be used as database to perform diagnostic BLASTp to classify GRAS genes from other non-model species.

Patterns of Substitution Rate Variation at Many Nuclear Loci in Two Species Trios in the Brassicaceae Partitioned with ANOVA

There are marked variations among loci and among lineages in rates of nucleotide substitution. The generation time hypothesis (GTH) is a neutral explanation for substitution rate heterogeneity that has genomewide application, predicting that species with shorter generation times accumulate DNA sequence substitutions faster than species with longer generation times do since faster genome replication provides more opportunities for mutations to occur and reach fixation by genetic drift. Relatively few studies have rigorously evaluated the GTH in plants, and there are numerous alternative hypotheses for plant substitution rate variation. One major challenge has been finding pairs of closely related plant species with contrasting generation times and appropriate outgroup taxa that all also have DNA sequence data for numerous loci. To test for causes of rate variation, we obtained sequence data for 256 genes for Arabidopsis thaliana, normally reproducing every year, and the biennial Arabidopsis lyrata with three closely related outgroup taxa (Brassica rapa, Capsella grandiflora, and Neslia paniculata) as well as the biennial Brassica oleracea and the annual B. rapa lineage with the outgroup N. paniculata. A sign test indicated that more loci than expected by chance have faster rates of substitution on the branch leading to the annual than to the perennial for one three-species trio but not another. Tajima's 1D and 2D tests, and a likelihood ratio test that incorporated saturation correction, rejected rate homogeneity for up to 26 genes (up to 14 genes when correcting for multiple tests), consistently showing faster rates for the annual lineage in the Arabidopsis species trio. ANOVA showed significant rate heterogeneity between the Arabidopsis and Brassica species trios (about 6 % of rate variation) and among loci (about 26-32 % of rate variation). The lineage-by-locus interaction which would be caused by locus- and lineage-specific natural selection explained about 13 % of substitution rate variation in one ANOVA model using substitution rates from genes partitioned into odd and even codons but was not a significant effect without partitioned genes. Annual/perennial lineage and species trio by annual/perennial lineage each explained about 1 % of substitution rate variation.

(Why) Does Evolution Favour Embryogenesis?

Complex multicellular organisms typically possess life cycles in which zygotes (formed by gamete fusion) and meiosis occur. Canonical animal embryogenesis describes development from zygote to birth. It involves polarisation of the egg/zygote, asymmetric cell divisions, establishment of axes, symmetry breaking, formation of organs, and parental nutrition (at least in early stages). Similar developmental patterns have independently evolved in other eukaryotic lineages, including land plants and brown algae. The question arises whether embryo-like structures and associated developmental processes recurrently emerge because they are local optima of the evolutionary landscape. To understand which evolutionary principles govern complex multicellularity, we need to analyse why and how similar processes evolve convergently - von Baer's and Haeckel's phylotypic stage revisited in other phyla.

ITS and trnH-psbA as Efficient DNA Barcodes to Identify Threatened Commercial Woody Angiosperms from Southern Brazilian Atlantic Rainforests

The Araucaria Forests in southern Brazil are part of the Atlantic Rainforest, a key hotspot for global biodiversity. This habitat has experienced extensive losses of vegetation cover due to commercial logging and the intense use of wood resources for construction and furniture manufacturing. The absence of precise taxonomic tools for identifying Araucaria Forest tree species motivated us to test the ability of DNA barcoding to distinguish species exploited for wood resources and its suitability for use as an alternative testing technique for the inspection of illegal timber shipments. We tested three cpDNA regions (matK, trnH-psbA, and rbcL) and nrITS according to criteria determined by The Consortium for the Barcode of Life (CBOL). The efficiency of each marker and selected marker combinations were evaluated for 30 commercially valuable woody species in multiple populations, with a special focus on Lauraceae species. Inter- and intraspecific distances, species discrimination rates, and ability to recover species-specific clusters were evaluated. Among the regions and different combinations, ITS was the most efficient for identifying species based on the 'best close match' test; similarly, the trnH-psbA + ITS combination also demonstrated satisfactory results. When combining trnH-psbA + ITS, Maximum Likelihood analysis demonstrated a more resolved topology for internal branches, with 91% of species-specific clusters. DNA barcoding was found to be a practical and rapid method for identifying major threatened woody angiosperms from Araucaria Forests such as Lauraceae species, presenting a high confidence for recognizing members of Ocotea. These molecular tools can assist in screening those botanical families that are most targeted by the timber industry in southern Brazil and detecting certain species protected by Brazilian legislation and could be a useful tool for monitoring wood exploitation.

Comparative Genomics of NAC Transcriptional Factors in Angiosperms: Implications for the Adaptation and Diversification of Flowering Plants

NAC proteins constitute one of the largest groups of plant-specific transcription factors and are known to play essential roles in various developmental processes. They are also important in plant responses to stresses such as drought, soil salinity, cold, and heat, which adversely affect growth. The current knowledge regarding the distribution of NAC proteins in plant lineages comes from relatively small samplings from the available data. In the present study, we broadened the number of plant species containing the NAC family origin and evolution to shed new light on the evolutionary history of this family in angiosperms. A comparative genome analysis was performed on 24 land plant species, and NAC ortholog groups were identified by means of bidirectional BLAST hits. Large NAC gene families are found in those species that have experienced more whole-genome duplication events, pointing to an expansion of the NAC family with divergent functions in flowering plants. A total of 3,187 NAC transcription factors that clustered into six major groups were used in the phylogenetic analysis. Many orthologous groups were found in the monocot and eudicot lineages, but only five orthologous groups were found between P. patens and each representative taxa of flowering plants. These groups were called basal orthologous groups and likely expanded into more recent taxa to cope with their environmental needs. This analysis on the angiosperm NAC family represents an effort to grasp the evolutionary and functional diversity within this gene family while providing a basis for further functional research on vascular plant gene families.

Multiple glacial refugia lead to genetic structuring and the potential for reproductive isolation in a herbaceous plant

PREMISE OF THE STUDY

Glacial cycles have influenced the genetic structure of many species. In addition to facilitating genetic divergence, isolation in multiple glacial refugia may have contributed to the development of genetic incompatibility and reproductive isolation. We examined the phylogeography of Campanulastrum americanum, a monocarpic herbaceous plant that exhibits strong intraspecific reproductive isolation, to determine whether the current genetic structure reflects a history of multiple glacial refugia.

METHODS

Chloroplast loci and nuclear RAD sequencing were used to characterize the range-wide phylogeography of C. americanum, in order to determine locations of potential glacial refugia and recolonization routes. Potential locations of refugia during the Last Glacial Maximum were also identified using ecological niche modeling.

KEY RESULTS

Together, the chloroplast and nuclear phylogenies found support for three geographically structured, genetically divergent lineages, among which gene flow appears to be restricted. The distribution of these lineages indicates that C. americanum survived the Last Glacial Maximum in at least three refugia located in the Appalachians and on the Atlantic and Gulf coasts. The ecological niche model also supported the existence of multiple refugia.

CONCLUSIONS

The isolation of populations of C. americanum in multiple refugia has led to a degree of phylogeographic structure greater than that found in most previously studied plants in eastern North America, which may be attributable to its short generation time. Reproductively isolated populations of C. americanum belong to divergent lineages, which suggests that survival in multiple glacial refugia contributed to the development of reproductive isolation in this species.

Plastid phylogenomics of the cool-season grass subfamily: clarification of relationships among early-diverging tribes

Whole plastid genomes are being sequenced rapidly from across the green plant tree of life, and phylogenetic analyses of these are increasing resolution and support for relationships that have varied among or been unresolved in earlier single- and multi-gene studies. Pooideae, the cool-season grass lineage, is the largest of the 12 grass subfamilies and includes important temperate cereals, turf grasses and forage species. Although numerous studies of the phylogeny of the subfamily have been undertaken, relationships among some 'early-diverging' tribes conflict among studies, and some relationships among subtribes of Poeae have not yet been resolved. To address these issues, we newly sequenced 25 whole plastomes, which showed rearrangements typical of Poaceae. These plastomes represent 9 tribes and 11 subtribes of Pooideae, and were analysed with 20 existing plastomes for the subfamily. Maximum likelihood (ML), maximum parsimony (MP) and Bayesian inference (BI) robustly resolve most deep relationships in the subfamily. Complete plastome data provide increased nodal support compared with protein-coding data alone at nodes that are not maximally supported. Following the divergence of Brachyelytrum, Phaenospermateae, Brylkinieae-Meliceae and Ampelodesmeae-Stipeae are the successive sister groups of the rest of the subfamily. Ampelodesmeae are nested within Stipeae in the plastome trees, consistent with its hybrid origin between a phaenospermatoid and a stipoid grass (the maternal parent). The core Pooideae are strongly supported and include Brachypodieae, a Bromeae-Triticeae clade and Poeae. Within Poeae, a novel sister group relationship between Phalaridinae and Torreyochloinae is found, and the relative branching order of this clade and Aveninae, with respect to an Agrostidinae-Brizinae clade, are discordant between MP and ML/BI trees. Maximum likelihood and Bayesian analyses strongly support Airinae and Holcinae as the successive sister groups of a Dactylidinae-Loliinae clade.

Dissecting Molecular Evolution in the Highly Diverse Plant Clade Caryophyllales Using Transcriptome Sequencing

Many phylogenomic studies based on transcriptomes have been limited to “single-copy” genes due to methodological challenges in homology and orthology inferences. Only a relatively small number of studies have explored analyses beyond reconstructing species relationships. We sampled 69 transcriptomes in the hyperdiverse plant clade Caryophyllales and 27 outgroups from annotated genomes across eudicots. Using a combined similarity- and phylogenetic tree-based approach, we recovered 10,960 homolog groups, where each was represented by at least eight ingroup taxa. By decomposing these homolog trees, and taking gene duplications into account, we obtained 17,273 ortholog groups, where each was represented by at least ten ingroup taxa. We reconstructed the species phylogeny using a 1,122-gene data set with a gene occupancy of 92.1%. From the homolog trees, we found that both synonymous and nonsynonymous substitution rates in herbaceous lineages are up to three times as fast as in their woody relatives. This is the first time such a pattern has been shown across thousands of nuclear genes with dense taxon sampling. We also pinpointed regions of the Caryophyllales tree that were characterized by relatively high frequencies of gene duplication, including three previously unrecognized whole-genome duplications. By further combining information from homolog tree topology and synonymous distance between paralog pairs, phylogenetic locations for 13 putative genome duplication events were identified. Genes that experienced the greatest gene family expansion were concentrated among those involved in signal transduction and oxidoreduction, including a cytochrome P450 gene that encodes a key enzyme in the betalain synthesis pathway. Our approach demonstrates a new approach for functional phylogenomic analysis in nonmodel species that is based on homolog groups in addition to inferred ortholog groups.

A generation-time effect on the rate of molecular evolution in bacteria

Molecular evolutionary rate varies significantly among species and a strict global molecular clock has been rejected across the tree of life. Generation time is one primary life-history trait that influences the molecular evolutionary rate. Theory predicts that organisms with shorter generation times evolve faster because of the accumulation of more DNA replication errors per unit time. Although the generation-time effect has been demonstrated consistently in plants and animals, the evidence of its existence in bacteria is lacking. The bacterial phylum Firmicutes offers an excellent system for testing generation-time effect because some of its members can enter a dormant, nonreproductive endospore state in response to harsh environmental conditions. It follows that spore-forming bacteria would--with their longer generation times--evolve more slowly than their nonspore-forming relatives. It is therefore surprising that a previous study found no generation-time effect in Firmicutes. Using a phylogenetic comparative approach and leveraging on a large number of Firmicutes genomes, we found sporulation significantly reduces the genome-wide spontaneous DNA mutation rate and protein evolutionary rate. Contrary to the previous study, our results provide strong evidence that the evolutionary rates of bacteria, like those of plants and animals, are influenced by generation time.

Genomic analysis of NAC transcription factors in banana (Musa acuminata) and definition of NAC orthologous groups for monocots and dicots

Identifying the molecular mechanisms underlying tolerance to abiotic stresses is important in crop breeding. A comprehensive understanding of the gene families associated with drought tolerance is therefore highly relevant. NAC transcription factors form a large plant-specific gene family involved in the regulation of tissue development and responses to biotic and abiotic stresses. The main goal of this study was to set up a framework of orthologous groups determined by an expert sequence comparison of NAC genes from both monocots and dicots. In order to clarify the orthologous relationships among NAC genes of different species, we performed an in-depth comparative study of four divergent taxa, in dicots and monocots, whose genomes have already been completely sequenced: Arabidopsis thaliana, Vitis vinifera, Musa acuminata and Oryza sativa. Due to independent evolution, NAC copy number is highly variable in these plant genomes. Based on an expert NAC sequence comparison, we propose forty orthologous groups of NAC sequences that were probably derived from an ancestor gene present in the most recent common ancestor of dicots and monocots. These orthologous groups provide a curated resource for large-scale protein sequence annotation of NAC transcription factors. The established orthology relationships also provide a useful reference for NAC function studies in newly sequenced genomes such as M. acuminata and other plant species.

Using multi-locus allelic sequence data to estimate genetic divergence among four Lilium (Liliaceae) cultivars

Next Generation Sequencing (NGS) may enable estimating relationships among genotypes using allelic variation of multiple nuclear genes simultaneously. We explored the potential and caveats of this strategy in four genetically distant Lilium cultivars to estimate their genetic divergence from transcriptome sequences using three approaches: POFAD (Phylogeny of Organisms from Allelic Data, uses allelic information of sequence data), RAxML (Randomized Accelerated Maximum Likelihood, tree building based on concatenated consensus sequences) and Consensus Network (constructing a network summarizing among gene tree conflicts). Twenty six gene contigs were chosen based on the presence of orthologous sequences in all cultivars, seven of which also had an orthologous sequence in Tulipa, used as out-group. The three approaches generated the same topology. Although the resolution offered by these approaches is high, in this case there was no extra benefit in using allelic information. We conclude that these 26 genes can be widely applied to construct a species tree for the genus Lilium.

The timetable for allopolyploidy in flowering plants

BACKGROUND

Our understanding of the processes and dynamics of allopolyploid speciation, the long-term consequences of ploidal change, and the genetic and chromosomal changes in new emerged allopolyploids has substantially increased during the past few decades. Yet we remain uncertain about the time since lineage divergence when two taxa are capable of spawning such entities. Indeed, the matter has seemed intractable. Knowledge of the window of opportunity for allopolyploid production is very important because it provides temporal insight into a key evolutionary process, and a temporal reference against which other modes of speciation may be measured.

SCOPE

This Viewpoint paper reviews and integrates published information on the crossability of herbaceous species and the fertility of their hybrids in relation to species' divergence times. Despite limitations in methodology and sampling, the estimated times to hybrid sterility are somewhat congruent across disparate lineages. Whereas the waiting time for hybrid sterility is roughly 4-5 million years, the waiting time for cross-incompatibility is roughly 8-10 million years, sometimes considerably more. Strict allopolyploids may be formed in the intervening time window. The progenitors of several allopolyploids diverged between 4 and 6 million years before allopolyploid synthesis, as expected. This is the first study to propose a general temporal framework for strict allopolyploidy. This Viewpoint paper hopefully will stimulate interest in studying the tempo of speciation and the tempo of reproductive isolation in general.

Differences in evolution rates among eudicotyledon species observed by analysis of protein divergence

Genome evolution rates can vary considerably among plants. In particular, a correlation has often been reported between the evolution rate and annual/perennial habit, possibly associated with differences in generation time. For example, among the rosid species whose genome is fully sequenced, Vitis vinifera, a perennial species, was shown to have the genome that evolved the slowest. In order to extend knowledge of evolution rates to the asterid clade, one of the two major clades of core eudicotyledonous, the protein evolution rates in three asterid species, one perennial (Coffea canephora) and two annual species (Solanum lycopersicum and Mimulus guttatus), were investigated and compared with V. vinifera. Significant differences were observed among these species, and the proteins that evolved the most slowly were those of V. vinifera. Among the species belonging to the asterid clade, C. canephora appears to have evolved more slowly than the others. These findings are consistent with a correlation between perennial habit and slow evolution rates. The C. canephora genome seems to be an appropriate model for paleogenomic studies of asterids.

A universal trend among proteomes indicates an oily last common ancestor

Despite progresses in ancestral protein sequence reconstruction, much needs to be unraveled about the nature of the putative last common ancestral proteome that served as the prototype of all extant lifeforms. Here, we present data that indicate a steady decline (oil escape) in proteome hydrophobicity over species evolvedness (node number) evident in 272 diverse proteomes, which indicates a highly hydrophobic (oily) last common ancestor (LCA). This trend, obtained from simple considerations (free from sequence reconstruction methods), was corroborated by regression studies within homologous and orthologous protein clusters as well as phylogenetic estimates of the ancestral oil content. While indicating an inherent irreversibility in molecular evolution, oil escape also serves as a rare and universal reaction-coordinate for evolution (reinforcing Darwin's principle of Common Descent), and may prove important in matters such as (i) explaining the emergence of intrinsically disordered proteins, (ii) developing composition- and speciation-based "global" molecular clocks, and (iii) improving the statistical methods for ancestral sequence reconstruction.

Increased phenotypic plasticity to climate may have boosted the invasion success of polyploid Centaurea stoebe

Phenotypic plasticity may allow organisms to cope with altered environmental conditions as e.g. after the introduction into a new range. In particular polyploid organisms, containing more than two sets of chromosomes, may show high levels of plasticity, which could in turn increase their environmental tolerance and invasiveness. Here, we studied the role of phenotypic plasticity in the invasion of Centaurea stoebe (Asteraceae), which in the native range in Europe occurs as diploids and tetraploids, whereas in the introduced range in North America so far only tetraploids have been found. In a common garden experiment at two sites in the native range, we grew half-sibs of the three geo-cytotypes (native European diploids, European tetraploids and invasive North American tetraploids) from a representative sample of 27 populations. We measured the level and the adaptive significance of phenotypic plasticity in eco-physiological and life-history traits in response to the contrasting climatic conditions at the two study sites as well as three different soil conditions in pots, simulating the most crucial abiotic differences between the native and introduced range. European tetraploids showed increased levels of phenotypic plasticity as compared to diploids in response to the different climatic conditions in traits associated with rapid growth and fast phenological development. Moreover, we found evidence for adaptive plasticity in these traits, which suggests that increased plasticity may have contributed to the invasion success of tetraploid C. stoebe by providing an advantage under the novel climatic conditions. However, in invasive tetraploids phenotypic plasticity was similar to that of native tetraploids, indicating no evolution of increased plasticity during invasions. Our findings provide the first empirical support for increased phenotypic plasticity associated with polyploids, which may contribute to their success as invasive species in novel environments.

Slow but not low: genomic comparisons reveal slower evolutionary rate and higher dN/dS in conifers compared to angiosperms

BACKGROUND

Comparative genomics can inform us about the processes of mutation and selection across diverse taxa. Among seed plants, gymnosperms have been lacking in genomic comparisons. Recent EST and full-length cDNA collections for two conifers, Sitka spruce (Picea sitchensis) and loblolly pine (Pinus taeda), together with full genome sequences for two angiosperms, Arabidopsis thaliana and poplar (Populus trichocarpa), offer an opportunity to infer the evolutionary processes underlying thousands of orthologous protein-coding genes in gymnosperms compared with an angiosperm orthologue set.

RESULTS

Based upon pairwise comparisons of 3,723 spruce and pine orthologues, we found an average synonymous genetic distance (dS) of 0.191, and an average dN/dS ratio of 0.314. Using a fossil-established divergence time of 140 million years between spruce and pine, we extrapolated a nucleotide substitution rate of 0.68 × 10(-9) synonymous substitutions per site per year. When compared to angiosperms, this indicates a dramatically slower rate of nucleotide substitution rates in conifers: on average 15-fold. Coincidentally, we found a three-fold higher dN/dS for the spruce-pine lineage compared to the poplar-Arabidopsis lineage. This joint occurrence of a slower evolutionary rate in conifers with higher dN/dS, and possibly positive selection, showcases the uniqueness of conifer genome evolution.

CONCLUSIONS

Our results are in line with documented reduced nucleotide diversity, conservative genome evolution and low rates of diversification in conifers on the one hand and numerous examples of local adaptation in conifers on the other hand. We propose that reduced levels of nucleotide mutation in large and long-lived conifer trees, coupled with large effective population size, were the main factors leading to slow substitution rates but retention of beneficial mutations.

Genome evolution in diploid and tetraploid Coffea species as revealed by comparative analysis of orthologous genome segments

Sequence comparison of orthologous regions enables estimation of the divergence between genomes, analysis of their evolution and detection of particular features of the genomes, such as sequence rearrangements and transposable elements. Despite the economic importance of Coffea species, little genomic information is currently available. Coffea is a relatively young genus that includes more than one hundred diploid species and a single tetraploid species. Three Coffea orthologous regions of 470-900 kb were analyzed and compared: both subgenomes of allotetraploid Coffea arabica (contributed by the diploid species Coffea eugenioides and Coffea canephora) and the genome of diploid C. canephora. Sequence divergence was calculated on global alignments or on coding and non-coding sequences separately. A search for transposable elements detected 43 retrotransposons and 198 transposons in the sequences analyzed. Comparative insertion analysis made it possible to locate 165 TE insertions in the phylogenetic tree of the three genomes/subgenomes. In the tetraploid C. arabica, a homoeologous non-reciprocal transposition (HNRT) was detected and characterized: a 50 kb region of the C. eugenioides derived subgenome replaced the C. canephora derived counterpart. Comparative sequence analysis on three Coffea genomes/subgenomes revealed almost perfect gene synteny, low sequence divergence and a high number of shared transposable elements. Compared to the results of similar analysis in other genera (Aegilops/Triticum and Oryza), Coffea genomes/subgenomes appeared to be dramatically less diverged, which is consistent with the relatively recent radiation of the Coffea genus. Based on nucleotide substitution frequency, the HNRT was dated at 10,000-50,000 years BP, which is also the most recent estimation of the origin of C. arabica.

Sequencing and analysis of plastid genome in mycoheterotrophic orchid Neottia nidus-avis

Plastids are the semiautonomous organelles that possess their own genome inherited from the cyanobacterial ancestor. The primary function of plastids is photosynthesis so the structure and evolution of plastid genomes are extensively studied in photosynthetic plants. In contrast, little is known about the plastomes of nonphotosynthetic species. In higher plants, plastid genome sequences are available for only three strictly nonphotosynthetic species, the liverwort Aneura mirabilis and two flowering plants, Epifagus virginiana and Rhizanthella gardneri. We report here the complete sequence of a plastid genome of nonphotosynthetic mycoheterotrophic orchid Neottia nidus-avis, determined using 454 pyrosequencing technology. It was found to be reduced in both genome size and gene content; this reduction is however not as drastic as in the other nonphotosynthetic orchid, R. gardneri. Neottia plastome lacks all genes encoding photosynthetic proteins, RNA polymerase subunits but retains most genes of translational apparatus. Those genes that are retained have an increased rate of both synonymous and nonsynonymous substitutions but do not exhibit relaxation of purifying selection either in Neottia or in Rhizanthella.