Evolution and biogeography of gymnosperms

Functional and morphological evolution in gymnosperms: A portrait of implicated gene families

Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water-conducting systems, secondary metabolism, stress defense mechanisms, and small RNA-mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole-genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.

Phylogenetically distant group of terpene synthases participates in cadinene and cedrane-type sesquiterpenes accumulation in Taiwania cryptomerioides

Along with the species evolution, plants have evolved ways to produce a different collection of terpenoids to accommodate its biotic and abiotic environment, and terpene synthase (TPS) is one of the major contributors to various terpene compounds. The timber of a monotypic and relictual conifer species of Cupressace, Taiwania cryptomerioides, has excellent durability, and one of the essential factors for Taiwania to resist decay and insect pests is sesquiterpene. Compared to other conifers, Taiwania has much higher abundance of cadinene-type sesquiterpenes, and the presence of cedrene-type sesquiterpenes. To understand sesquiterpene biosynthesis in Taiwania, we functionally characterized 10 T. cryptomerioides TPSs (TcTPSs) in vivo or in planta, which could catalyze sesquiterpene formation and potentially are involved in biosynthesis of diverse sesquiterpenoids in Taiwania. The distant phylogenetic relationship and the intron loss event of TcTPSs correlate to the differentiation of chemical profile Taiwania compared to other conifers. Furthermore, we identified TcTPS3 and TcTPS12 as δ-cadinene synthase, and TcTPS6 as cedrol synthase, which demonstrates the important contributions of dynamic evolution in TPSs to the chemical diversity in plants. Combining with functional characterization and comparison of catalytic residues, we conclude at least three catalytic routes for sesquiterpene biosynthesis in this species, and the skeleton diversity has been expended in T. cryptomeriodes.

Somatic mutations substantially increase the per-generation mutation rate in the conifer Picea sitchensis

The rates and biological significance of somatic mutations have long been a subject of debate. Somatic mutations in plants are expected to accumulate with vegetative growth and time, yet rates of somatic mutations are unknown for conifers, which can reach exceptional sizes and ages. We investigated somatic mutation rates in the conifer Sitka spruce (Picea sitchensis (Bong.) Carr.) by analyzing a total of 276 Gb of nuclear DNA from the tops and bottoms of 20 old‐growth trees averaging 76 m in height. We estimate a somatic base substitution rate of 2.7 × 10⁻⁸ per base pair within a generation. To date, this is one of the highest estimated per‐generation rates of mutation among eukaryotes, indicating that somatic mutations contribute substantially to the total per‐generation mutation rate in conifers. Nevertheless, as the sampled trees are centuries old, the per‐year rate is low in comparison with nontree taxa. We argue that although somatic mutations raise genetic load in conifers, they generate important genetic variation and enable selection both among cell lineages within individual trees and among offspring.

Remarkable variation of ribosomal DNA organization and copy number in gnetophytes, a distinct lineage of gymnosperms

INTRODUCTION

Gnetophytes, comprising the genera Ephedra, Gnetum and Welwitschia, are an understudied, enigmatic lineage of gymnosperms with a controversial phylogenetic relationship to other seed plants. Here we examined the organization of ribosomal DNA (rDNA) across representative species.

METHODS

We applied high-throughput sequencing approaches to isolate and reconstruct rDNA units and to determine their intragenomic homogeneity. In addition, fluorescent in situ hybridization and Southern blot hybridization techniques were used to reveal the chromosome and genomic organization of rDNA.

KEY RESULTS

The 5S and 35S rRNA genes were separate (S-type) in Gnetum montanum, Gnetum gnemon and Welwitschia mirabilis and linked (L-type) in Ephedra altissima. There was considerable variability in 5S rDNA abundance, ranging from as few as ~4000 (W. mirabilis) to >100 000 (G. montanum) copies. A similar large variation was also observed in 5S rDNA locus numbers (two to 16 sites per diploid cell). 5S rRNA pseudogenes were interspersed between functional genes forming a single unit in E. altissima and G. montanum. Their copy number was comparable or even higher than that of functional 5S rRNA genes. In E. altissima internal transcribed spacers of 35S rDNA were long and intrinsically repetitive while in G. montanum and W. mirabilis they were short without the subrepeats.

CONCLUSIONS

Gnetophytes are distinct from other gymnosperms and angiosperms as they display surprisingly large variability in rDNA organization and rDNA copy and locus numbers between genera, with no relationship between copy numbers and genome sizes apparent. Concerted evolution of 5S rDNA units seems to have led to the amplification of 5S pseudogenes in both G. montanum and E. altissima. Evolutionary patterns of rDNA show both gymnosperm and angiosperm features underlining the diversity of the group.

Cryptic speciation in the Chinese white pine (Pinus armandii): Implications for the high species diversity of conifers in the Hengduan Mountains, a global biodiversity hotspot

Conifers are the largest and ecologically and economically most important component group of the gymnosperms. Despite their slow rate of molecular evolution, rapid and recent diversification was unexpectedly prevalent in this ancient group in the Hengduan Mountains, a world's biodiversity hotspot and gymnosperm diversity center in Southwest China. In this study, we investigated the underlying mechanisms and disentangled the interactions of geography and ecology in speciation and evolution in Pinus armandii, an important forest tree species endemic to China, by integrating analyses of population transcriptomics, population genetics and ecological niche modeling. Many lines of evidence suggest that cryptic speciation has occurred in P. armandii. During the process, geologically induced formation of Mount Gongga and other massive peaks might trigger the initial vicariance isolation of the northern and southern subdivisions, and ecologically based selection then reinforced their differentiation and local adaptation. Our ecological niche analysis and earlier reciprocal transplant experiments in P. armandii provided convincing evidences for the critical role of ecology in the process of speciation. These findings suggest that both geography and ecology contributed significantly to the abundance of very recent and rapid species divergences, which promoted the rising of the extremely high conifer diversity in the Hengduan Mountains.

Phylogenomics resolves the deep phylogeny of seed plants and indicates partial convergent or homoplastic evolution between Gnetales and angiosperms

After decades of molecular phylogenetic studies, the deep phylogeny of gymnosperms has not been resolved, and the phylogenetic placement of Gnetales remains one of the most controversial issues in seed plant evolution. To resolve the deep phylogeny of seed plants and to address the sources of phylogenetic conflict, we conducted a phylotranscriptomic study with a sampling of all 13 families of gymnosperms and main lineages of angiosperms. Multiple datasets containing up to 1 296 042 sites across 1308 loci were analysed, using concatenation and coalescence approaches. Our study generated a consistent and well-resolved phylogeny of seed plants, which places Gnetales as sister to Pinaceae and thus supports the Gnepine hypothesis. Cycads plus Ginkgo is sister to the remaining gymnosperms. We also found that Gnetales and angiosperms have similar molecular evolutionary rates, which are much higher than those of other gymnosperms. This implies that Gnetales and angiosperms might have experienced similar selective pressures in evolutionary histories. Convergent molecular evolution or homoplasy is partially responsible for the phylogenetic conflicts in seed plants. Our study provides a robustly reconstructed backbone phylogeny that is important for future molecular and morphological studies of seed plants, in particular gymnosperms, in the light of evolution.

A transcriptome-based resolution for a key taxonomic controversy in Cupressaceae

BACKGROUND AND AIMS

Rapid evolutionary divergence and reticulate evolution may result in phylogenetic relationships that are difficult to resolve using small nucleotide sequence data sets. Next-generation sequencing methods can generate larger data sets that are better suited to solving these puzzles. One major and long-standing controversy in conifers concerns generic relationships within the subfamily Cupressoideae (105 species, approx. 1/6 of all conifers) of Cupressaceae, in particular the relationship between Juniperus, Cupressus and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) clade. Here we attempt to resolve this question using transcriptome-derived data.

METHODS

Transcriptome sequences of 20 species from Cupressoideae were collected. Using MarkerMiner, single-copy nuclear (SCN) genes were extracted. These were applied to estimate phylogenies based on concatenated data, species trees and a phylogenetic network. We further examined the effect of alternative backbone topologies on downstream analyses, including biogeographic inference and dating analysis.

RESULTS

Based on the 73 SCN genes (>200 000 bp total alignment length) we considered, all tree-building methods lent strong support for the relationship (HCX, (Juniperus, Cupressus)); however, strongly supported conflicts among individual gene trees were also detected. Molecular dating suggests that these three lineages shared a most recent common ancestor approx. 60 million years ago (Mya), and that Juniperus and Cupressus diverged about 56 Mya. Ancestral area reconstructions (AARs) suggest an Asian origin for the entire clade, with subsequent dispersal to North America, Europe and Africa.

CONCLUSIONS

Our analysis of SCN genes resolves a controversial phylogenetic relationship in the Cupressoideae, a major clade of conifers, and suggests that rapid evolutionary divergence and incomplete lineage sorting probably acted together as the source for conflicting phylogenetic inferences between gene trees and between our robust results and recently published studies. Our updated backbone topology has not substantially altered molecular dating estimates relative to previous studies; however, application of the latest AAR approaches has yielded a clearer picture of the biogeographic history of Cupressoideae.

Pronounced genetic differentiation in Fokienia hodginsii revealed by simple sequence repeat markers

Fokienia hodginsii is a Tertiary relict conifer of the monotypic genus Fokienia (Cupressaceae s.l.). Currently, the species is distributed in southern China, northern Vietnam, and northern Laos and listed as a "near threatened" species by the IUCN. In this study, a total of 427 individuals of F. hodginsii were sampled from China and Vietnam to characterize its genetic diversity and population differentiation. Based on the profiles of 12 simple sequence repeat (SSR) markers, we observed a high level of genetic diversity in F. hodginsii at the species level (H e =0.635), albeit slightly lower than that of its sister species Chamaecyparis obtusa. Signals of bottleneck events were detected in the populations GXDMS, GXHJ, V-PXB, and V-HB, probably due to Pleistocene glaciations or overexploitation in recent years. Pronounced genetic differentiation (F st   = 0.157) was found in this species. The inbreeding index (F is  = 0.176 ± 0.024) indicated that F. hodginsii has a mixed mating system. Significant correlation was found between the pairwise genetic differentiation and geographic distance (r = 0.882, p = 0.01), suggesting that genetic differentiation among the populations follows the model of isolation by distance (IBD). STRUCTURE analysis and principal coordinate analysis revealed that these populations were divided into four groups: the western China group located mainly in the Yunnan-Guizhou Plateau, the central China group located mostly in the Luoxiao Mountains and Nanling Mountains, the eastern China group located in the Wuyi Mountains and the Vietnam group containing two populations in Vietnam. The different terrains and elevations of populations may be the most likely factors leading to the differentiation between the western China group and the central China group, while the geographic isolation caused by the lack of appropriate habitats may greatly contribute to the differentiation between the central China group and the eastern China group. Based on the results, some conservation suggestions for this species are provided, such as establishing seed orchards and multiple nature reserves.

Biogeography of Coptis Salisb. (Ranunculales, Ranunculaceae, Coptidoideae), an Eastern Asian and North American genus

BACKGROUND

Numerous studies have favored dispersal (colonization) over vicariance (past fragmentation) events to explain eastern Asian-North American distribution patterns. In plants, however the disjunction between eastern Asia and western North America has been rarely examined using the integration of phylogenetic, molecular dating, and biogeographical methods. Meanwhile, the biogeographic patterns within eastern Asia remain poorly understood. The goldthread genus Coptis Salisb. includes 15 species disjunctly distributed in North America, Japan, mainland China, and Taiwan. We present a dated phylogeny for Coptis under the optimal clock model and infer its historical biogeography by comparing different biogeographic models.

RESULTS

The split of Coptis and Xanthorhiza Marshall occurred in the middle Miocene (ca. 15.47 Ma). Coptis started their diversification in the early late Miocene (ca. 9.55 Ma). A late Miocene vicariance event resulted in the eastern Asian and western North American disjunction in the genus. Within eastern Asia, dispersals from mainland Asia to Japan and from Japan to Taiwan occurred at ca. 4.85 Ma and at ca. 1.34 Ma, respectively.

CONCLUSIONS

Our analyses provide evidence that both vicariance and dispersal events have played important roles in shaping the current distribution and endemism of Coptis, likely resulting from eustatic sea-level changes, mountain formation processes and an increasing drier and cooler climate from the middle Miocene onwards.

A genome for gnetophytes and early evolution of seed plants

Gnetophytes are an enigmatic gymnosperm lineage comprising three genera, Gnetum, Welwitschia and Ephedra, which are morphologically distinct from all other seed plants. Their distinctiveness has triggered much debate as to their origin, evolution and phylogenetic placement among seed plants. To increase our understanding of the evolution of gnetophytes, and their relation to other seed plants, we report here a high-quality draft genome sequence for Gnetum montanum, the first for any gnetophyte. By using a novel genome assembly strategy to deal with high levels of heterozygosity, we assembled >4 Gb of sequence encoding 27,491 protein-coding genes. Comparative analysis of the G. montanum genome with other gymnosperm genomes unveiled some remarkable and distinctive genomic features, such as a diverse assemblage of retrotransposons with evidence for elevated frequencies of elimination rather than accumulation, considerable differences in intron architecture, including both length distribution and proportions of (retro) transposon elements, and distinctive patterns of proliferation of functional protein domains. Furthermore, a few gene families showed Gnetum-specific copy number expansions (for example, cellulose synthase) or contractions (for example, Late Embryogenesis Abundant protein), which could be connected with Gnetum's distinctive morphological innovations associated with their adaptation to warm, mesic environments. Overall, the G. montanum genome enables a better resolution of ancestral genomic features within seed plants, and the identification of genomic characters that distinguish Gnetum from other gymnosperms.

Phylogenetic Classification of Seed Plants of Taiwan

BACKGROUND

Biological classification, the hierarchical arrangement of scientific names of organisms, constitutes the core infrastructure of biological databases. For an efficient management of biological databases, adopting a stable and universal biological classification system is crucial. Currently in Taiwan Biodiversity Information Facility (TaiBIF; http://taibif.tw/ ), the national portal website that integrates Taiwan's biodiversity information databases, angiosperms are arranged according to Cronquist's System of Classification, which is not compatible with current trend of the Angiosperm Phylogeny Group (APG) classification. To consolidate the function and management of the database, TaiBIF is moving to adopt the APG IV classification and Christenhusz et al. (Phytotaxa 19:55-70, 2011)'s classification of gymnosperms, which we summarize as the Phylogenetic Classification of Seed Plants of Taiwan.

RESULTS

The Phylogenetic Classification of Seed Plants of Taiwan places gymnosperms in five families [vs. eight families in the Flora of Taiwan (FOT)] and angiosperms in 210 families (vs. 193 families in FOT). Three FOT gymnosperm families are synonymized in current treatment. Of the 210 APG IV families, familial circumscriptions of 114 families are identical with FOT and 50 families are recircumscription of FOT, with 46 families newly added. Of the 29 FOT families not included in current classification, two families are excluded and 27 families are synonymized.

CONCLUSIONS

The adoption of the Phylogenetic Classification of Seed Plants of Taiwan in TaiBIF will provide better service and efficient management of the nation's biodiversity information databases.

Single-Copy Genes as Molecular Markers for Phylogenomic Studies in Seed Plants

Phylogenetic relationships among seed plant taxa, especially within the gymnosperms, remain contested. In contrast to angiosperms, for which several genomic, transcriptomic and phylogenetic resources are available, there are few, if any, molecular markers that allow broad comparisons among gymnosperm species. With few gymnosperm genomes available, recently obtained transcriptomes in gymnosperms are a great addition to identifying single-copy gene families as molecular markers for phylogenomic analysis in seed plants. Taking advantage of an increasing number of available genomes and transcriptomes, we identified single-copy genes in a broad collection of seed plants and used these to infer phylogenetic relationships between major seed plant taxa. This study aims at extending the current phylogenetic toolkit for seed plants, assessing its ability for resolving seed plant phylogeny, and discussing potential factors affecting phylogenetic reconstruction. In total, we identified 3,072 single-copy genes in 31 gymnosperms and 2,156 single-copy genes in 34 angiosperms. All studied seed plants shared 1,469 single-copy genes, which are generally involved in functions like DNA metabolism, cell cycle, and photosynthesis. A selected set of 106 single-copy genes provided good resolution for the seed plant phylogeny except for gnetophytes. Although some of our analyses support a sister relationship between gnetophytes and other gymnosperms, phylogenetic trees from concatenated alignments without 3rd codon positions and amino acid alignments under the CAT + GTR model, support gnetophytes as a sister group to Pinaceae. Our phylogenomic analyses demonstrate that, in general, single-copy genes can uncover both recent and deep divergences of seed plant phylogeny.

Contrasting Rates of Molecular Evolution and Patterns of Selection among Gymnosperms and Flowering Plants

The majority of variation in rates of molecular evolution among seed plants remains both unexplored and unexplained. Although some attention has been given to flowering plants, reports of molecular evolutionary rates for their sister plant clade (gymnosperms) are scarce, and to our knowledge differences in molecular evolution among seed plant clades have never been tested in a phylogenetic framework. Angiosperms and gymnosperms differ in a number of features, of which contrasting reproductive biology, life spans, and population sizes are the most prominent. The highly conserved morphology of gymnosperms evidenced by similarity of extant species to fossil records and the high levels of macrosynteny at the genomic level have led scientists to believe that gymnosperms are slow-evolving plants, although some studies have offered contradictory results. Here, we used 31,968 nucleotide sites obtained from orthologous genes across a wide taxonomic sampling that includes representatives of most conifers, cycads, ginkgo, and many angiosperms with a sequenced genome. Our results suggest that angiosperms and gymnosperms differ considerably in their rates of molecular evolution per unit time, with gymnosperm rates being, on average, seven times lower than angiosperm species. Longer generation times and larger genome sizes are some of the factors explaining the slow rates of molecular evolution found in gymnosperms. In contrast to their slow rates of molecular evolution, gymnosperms possess higher substitution rate ratios than angiosperm taxa. Finally, our study suggests stronger and more efficient purifying and diversifying selection in gymnosperm than in angiosperm species, probably in relation to larger effective population sizes.

Global Analysis of Small RNA Dynamics during Seed Development of Picea glauca and Arabidopsis thaliana Populations Reveals Insights on their Evolutionary Trajectories

While DNA methylation carries genetic signals and is instrumental in the evolution of organismal complexity, small RNAs (sRNAs), ~18-24 ribonucleotide (nt) sequences, are crucial mediators of methylation as well as gene silencing. However, scant study deals with sRNA evolution via featuring their expression dynamics coupled with species of different evolutionary time. Here we report an atlas of sRNAs and microRNAs (miRNAs, single-stranded sRNAs) produced over time at seed-set of two major spermatophytes represented by populations of Picea glauca and Arabidopsis thaliana with different seed-set duration. We applied diverse profiling methods to examine sRNA and miRNA features, including size distribution, sequence conservation and reproduction-specific regulation, as well as to predict their putative targets. The top 27 most abundant miRNAs were highly overlapped between the two species (e.g., miR166,-319 and-396), but in P. glauca, they were less abundant and significantly less correlated with seed-set phases. The most abundant sRNAs in libraries were deeply conserved miRNAs in the plant kingdom for Arabidopsis but long sRNAs (24-nt) for P. glauca. We also found significant difference in normalized expression between populations for population-specific sRNAs but not for lineage-specific ones. Moreover, lineage-specific sRNAs were enriched in the 21-nt size class. This pattern is consistent in both species and alludes to a specific type of sRNAs (e.g., miRNA, tasiRNA) being selected for. In addition, we deemed 24 and 9 sRNAs in P. glauca and Arabidopsis, respectively, as sRNA candidates targeting known adaptive genes. Temperature had significant influence on selected gene and miRNA expression at seed development in both species. This study increases our integrated understanding of sRNA evolution and its potential link to genomic architecture (e.g., sRNA derivation from genome and sRNA-mediated genomic events) and organismal complexity (e.g., association between different sRNA expression and their functionality).

NorWood: a gene expression resource for evo-devo studies of conifer wood development

The secondary xylem of conifers is composed mainly of tracheids that differ anatomically and chemically from angiosperm xylem cells. There is currently no high-spatial-resolution data available profiling gene expression during wood formation for any coniferous species, which limits insight into tracheid development. RNA-sequencing data from replicated, high-spatial-resolution section series throughout the cambial and woody tissues of Picea abies were used to generate the NorWood.conGenIE.org web resource, which facilitates exploration of the associated gene expression profiles and co-expression networks. Integration within PlantGenIE.org enabled a comparative regulomics analysis, revealing divergent co-expression networks between P. abies and the two angiosperm species Arabidopsis thaliana and Populus tremula for the secondary cell wall (SCW) master regulator NAC Class IIB transcription factors. The SCW cellulose synthase genes (CesAs) were located in the neighbourhoods of the NAC factors in A. thaliana and P. tremula, but not in P. abies. The NorWood co-expression network enabled identification of potential SCW CesA regulators in P. abies. The NorWood web resource represents a powerful community tool for generating evo-devo insights into the divergence of wood formation between angiosperms and gymnosperms and for advancing understanding of the regulation of wood development in P. abies.

β-Ginkgotides: Hyperdisulfide-constrained peptides from Ginkgo biloba

Hyperdisulfide-constrained peptides are distinguished by their high stability and diverse functions. Thus far, these peptides have been reported from animals only but their occurrence in plants are rare. Here, we report the discovery, synthesis and characterization of a hyperdisulfide-constrained peptides family of approximately 2 kDa, β-ginkgotides (β-gB1 and β-gB2) from Ginkgo biloba. Proteomic analysis showed β-ginkgotides contain 18‒20 amino acids, of which 16 residues form a conserved six-cysteine core with a highly clustered cysteine spacing of C‒CC‒C‒CC, an arrangement that has not been reported in cysteine-rich peptides. Disulfide mapping revealed a novel disulfide connectivity of CysI‒IV, CysII‒VI and CysIII‒V. Oxidative folding of synthetic β-gB1 to the native form was obtained in 70% yield. The synthetic β-gB1 displays a compact structure with no regular secondary structural elements, as determined by NMR spectroscopy. Transcriptomic analysis showed precursor βgb1 has a four-domain architecture and revealed an additional 76 β-ginkgotide-like peptides in 59 different gymnosperms, but none in angiosperms. Phylogenetic clustering analysis demonstrated β-ginkgotides belong to a new cysteine-rich peptide family. β-Ginkgotide is resistant to thermal, chemical and proteolytic degradation. Together, β-ginkgotides represent the first-in-class hyperdisulfide-constrained peptide family from plants with a novel scaffold that could be useful for engineering metabolically stable peptidyl therapeutics.

Baptism by fire: the pivotal role of ancient conflagrations in evolution of the Earth's flora

Fire became a defining feature of the Earth's processes as soon as land plants evolved 420 million years ago and has played a major role in shaping the composition and physiognomy of many ecosystems ever since. However, there remains a general lack of appreciation of the place of fire in the origin, evolution, ecology and conservation of the Earth's biodiversity. We review the literature on the presence of fire throughout the Earth's history following the evolution of land plants and examine the evidence for the origin and evolution of adaptive functional traits, biomes and major plant groups in relation to fire. We show that: (1) fire activities have fluctuated throughout geological time due to variations in climate, and more importantly in atmospheric oxygen, as these affected fuel levels and flammability; (2) fire promoted the early evolution and spread of major terrestrial plant groups; (3) fire has shaped the floristics, structure and function of major global biomes; and (4) fire has initiated and maintained the evolution of a wide array of fire-adapted functional traits since the evolution of land plants. We conclude that fire has been a fundamental agent of natural selection on terrestrial plants throughout the history of life on the Earth's land surface. We suggest that a paradigm shift is required to reassess ecological and evolutionary theories that exclude a role for fire, and also there is a need to review fire-suppression policies on ecosystem management and biodiversity conservation in global fire-prone regions.

Insights into the Existence of Isomeric Plastomes in Cupressoideae (Cupressaceae)

The cypress family (Cupressaceae) possesses highly rearranged plastomes that lack a pair of large inverted repeats typically found in land plants. A few cypress species have been reported to contain isomeric plastomes, but whether the existence of isomeric plastomes is ubiquitous in the family remains to be investigated with a broader taxon sampling. In this study, we sequenced the complete plastomes of ten species in Cupressoideae, the largest cypress subfamily. Cupressoideae showed relatively accelerated rates of substitutions at both nonsynonymous and synonymous sites as compared with other subfamilies of Cupressaceae. Our PCR and read mapping analyses together suggested the existence of isomeric plastomes in eight of the ten sequenced Cupressoideae species. The isomeric plastomes were also detected in 176 individuals from nine wild populations of four Cupressoideae species. Within Calocedrus macrolepis, we discovered a new type of isomeric plastomes that was likely derived from homologous recombination mediated by an 11-bp repeat. We conclude that isomeric plastomes are commonly present in Cupressoideae, thereby contributing to increased plastomic complexity.

Multiple measures could alleviate long-branch attraction in phylogenomic reconstruction of Cupressoideae (Cupressaceae)

Long-branch attraction (LBA) is a major obstacle in phylogenetic reconstruction. The phylogenetic relationships among Juniperus (J), Cupressus (C) and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) subclades of Cupressoideae are controversial. Our initial analyses of plastid protein-coding gene matrix revealed both J and C with much longer stem branches than those of HCX, so their sister relationships may be attributed to LBA. We used multiple measures including data filtering and modifying, evolutionary model selection and coalescent phylogenetic reconstruction to alleviate the LBA artifact. Data filtering by strictly removing unreliable aligned regions and removing substitution saturation genes and rapidly evolving sites could significantly reduce branch lengths of subclades J and C and recovered a relationship of J (C, HCX). In addition, using coalescent phylogenetic reconstruction could elucidate the LBA artifact and recovered J (C, HCX). However, some valid methods for other taxa were inefficient in alleviating the LBA artifact in J-C-HCX. Different strategies should be carefully considered and justified to reduce LBA in phylogenetic reconstruction of different groups. Three subclades of J-C-HCX were estimated to have experienced ancient rapid divergence within a short period, which could be another major obstacle in resolving relationships. Furthermore, our plastid phylogenomic analyses fully resolved the intergeneric relationships of Cupressoideae.

DNA methylome of the 20-gigabase Norway spruce genome

DNA methylation plays important roles in many biological processes, such as silencing of transposable elements, imprinting, and regulating gene expression. Many studies of DNA methylation have shown its essential roles in angiosperms (flowering plants). However, few studies have examined the roles and patterns of DNA methylation in gymnosperms. Here, we present genome-wide high coverage single-base resolution methylation maps of Norway spruce (Picea abies) from both needles and somatic embryogenesis culture cells via whole genome bisulfite sequencing. On average, DNA methylation levels of CG and CHG of Norway spruce were higher than most other plants studied. CHH methylation was found at a relatively low level; however, at least one copy of most of the RNA-directed DNA methylation pathway genes was found in Norway spruce, and CHH methylation was correlated with levels of siRNAs. In comparison with needles, somatic embryogenesis culture cells that are used for clonally propagating spruce trees showed lower levels of CG and CHG methylation but higher level of CHH methylation, suggesting that like in other species, these culture cells show abnormal methylation patterns.

Vesicles Are Persistent Features of Different Plastids

Peripheral vesicles in plastids have been observed repeatedly, primarily in proplastids and developing chloroplasts, in which they are suggested to function in thylakoid biogenesis. Previous observations of vesicles in mature chloroplasts have mainly concerned low temperature pretreated plants occasionally treated with inhibitors blocking vesicle fusion. Here, we show that such vesicle-like structures occur not only in chloroplasts and proplastids, but also in etioplasts, etio-chloroplasts, leucoplasts, chromoplasts and even transforming desiccoplasts without any specific pretreatment. Observations are made both in C3 and C4 species, in different cell types (meristematic, epidermis, mesophyll, bundle sheath and secretory cells) and different organs (roots, stems, leaves, floral parts and fruits). Until recently not much focus has been given to the idea that vesicle transport in chloroplasts could be mediated by proteins, but recent data suggest that the vesicle system of chloroplasts has similarities with the cytosolic coat protein complex II system. All current data taken together support the idea of an ongoing, active and protein-mediated vesicle transport not only in chloroplasts but also in other plastids, obviously occurring regardless of chemical modifications, temperature and plastid developmental stage.

Evolutionary loss of light-harvesting proteins Lhcb6 and Lhcb3 in major land plant groups--break-up of current dogma

Photosynthesis in plants and algae relies on the coordinated function of photosystems (PS) I and II. Their efficiency is augmented by finely-tuned light-harvesting proteins (Lhcs) connected to them. The most recent Lhcs (in evolutionary terms), Lhcb6 and Lhcb3, evolved during the transition of plants from water to land and have so far been considered to be an essential characteristic of land plants. We used single particle electron microscopy and sequence analysis to study architecture and composition of PSII supercomplex from Norway spruce and related species. We have found that there are major land plant families that lack functional lhcb6 and lhcb3 genes, which notably changes the organization of PSII supercomplexes. The Lhcb6 and Lhcb3 proteins have been lost in the gymnosperm genera Picea and Pinus (family Pinaceae) and Gnetum (Gnetales). We also revealed that the absence of these proteins in Norway spruce modifies the PSII supercomplex in such a way that it resembles its counterpart in the alga Chlamydomonas reinhardtii, an evolutionarily older organism. Our results break a deep-rooted concept of Lhcb6 and Lhcb3 proteins being the essential characteristic of land plants, and beg the question of what the evolutionary benefit of their loss could be.

A high frequency of allopolyploid speciation in the gymnospermous genus Ephedra and its possible association with some biological and ecological features

The origin and evolution of polyploids have been studied extensively in angiosperms and ferns but very rarely in gymnosperms. With the exception of three species of conifers, all natural polyploid species of gymnosperms belong to Ephedra, in which more than half of the species show polyploid cytotypes. Here, we investigated the origin and evolution of polyploids of Ephedra distributed in the Qinghai–Tibetan Plateau (QTP) and neighbouring areas. Flow cytometry (FCM) was used to measure the ploidy levels of the sampled species that are represented by multiple individuals from different populations, and then, two single‐copy nuclear genes (LFY and DDB2) and two chloroplast DNA fragments were used to unravel the possible origins and maternal donors of the polyploids. The results indicate that the studied polyploid species are allopolyploids, and suggest that allotetraploidy is a dominant mode of speciation in Ephedra. The high percentage of polyploids in the genus could be related to some of its biological attributes such as vegetative propagation, a relatively high rate of unreduced gamete formation, and a small genome size relative to most other gymnosperms. Significant ecological divergences between allotetraploids and their putative progenitors were detected by PCAs and anova and Tukey's tests, with the exception of E. saxatilis. The overlap of geographical distributions and ecological niches of some diploid species could have provided opportunities for interspecific hybridization and allopolyploid speciation.

Does the legacy of historical biogeography shape current invasiveness in pines?

Why are some introduced species more successful at establishing and spreading than others? Until now, characteristics of extant species have been intensively investigated to answer this question. We propose to gain new insights on species invasiveness by exploring the long-term biogeographic and evolutionary history of lineages. We exemplify our approach using one of the best-studied invasive plant genera, Pinus. We notably estimated the historical biogeography of pines and the rates of trait evolution in pines. These estimates were analysed with regard to species invasiveness status. The results revealed that currently invasive species belong to lineages that were particularly successful at colonizing new regions in the past. We also showed that highly mobile lineages had faster rates of niche evolution, but that these rates are poor proxies for species adaptive potential in invaded regions (estimated by niche shift among native and invaded regions). In summary, working at the interface of ecology, historical biogeography and evolutionary history offers stimulating perspectives to improve our understanding of the drivers of invasion success.

Sequenced genomes and rapidly emerging technologies pave the way for conifer evolutionary developmental biology

Conifers, Ginkgo, cycads and gnetophytes comprise the four groups of extant gymnosperms holding a unique position of sharing common ancestry with the angiosperms. Comparative studies of gymnosperms and angiosperms are the key to a better understanding of ancient seed plant morphologies, how they have shifted over evolution to shape modern day species, and how the genes governing these morphologies have evolved. However, conifers and other gymnosperms have been notoriously difficult to study due to their long generation times, inaccessibility to genetic experimentation and unavailable genome sequences. Now, with three draft genomes from spruces and pines, rapid advances in next generation sequencing methods for genome wide expression analyses, and enhanced methods for genetic transformation, we are much better equipped to address a number of key evolutionary questions relating to seed plant evolution. In this mini-review we highlight recent progress in conifer developmental biology relevant to evo-devo questions. We discuss how genome sequence data and novel techniques might allow us to explore genetic variation and naturally occurring conifer mutants, approaches to reduce long generation times to allow for genetic studies in conifers, and other potential upcoming research avenues utilizing current and emergent techniques. Results from developmental studies of conifers and other gymnosperms in comparison to those in angiosperms will provide information to trace core molecular developmental control tool kits of ancestral seed plants, but foremost they will greatly improve our understanding of the biology of conifers and other gymnosperms in their own right.

Evidence of intense chromosomal shuffling during conifer evolution

Although recent advances have been gained on genome evolution in angiosperm lineages, virtually nothing is known about karyotype evolution in the other group of seed plants, the gymnosperms. Here, we used high-density gene-based linkage mapping to compare the karyotype structure of two families of conifers (the most abundant group of gymnosperms) separated around 290 Ma: Pinaceae and Cupressaceae. We propose for the first time a model based on the fusion of 20 ancestral chromosomal blocks that may have shaped the modern karyotpes of Pinaceae (with n = 12) and Cupressaceae (with n = 11). The considerable difference in modern genome organization between these two lineages contrasts strongly with the remarkable level of synteny already reported within the Pinaceae. It also suggests a convergent evolutionary mechanism of chromosomal block shuffling that has shaped the genomes of the spermatophytes.

Exploring the Formation of a Disjunctive Pattern between Eastern Asia and North America Based on Fossil Evidence from Thuja (Cupressaceae)

Thuja, a genus of Cupressaceae comprising five extant species, presently occurs in both East Asia (3 species) and North America (2 species) and has a long fossil record from Paleocene to Pleistocene in the Northern Hemisphere. Two distinct hypotheses have been proposed to account for the origin and present distribution of this genus. Here we recognize and describe T. sutchuenensis Franch., a new fossil Thuja from the late Pliocene sediments of Zhangcun, Shanxi, North China, based on detailed comparisons with all living species and other fossil ones, integrate the global fossil records of this genus plotted in a set of paleomaps from different time intervals, which show that Thuja probably first appeared at high latitudes of North America in or before the Paleocene. This genus reached Greenland in the Paleocene, then arrived in eastern Asia in the Miocene via the land connection between East Asia and western North America. In the late Pliocene, it migrated into the interior of China. With the Quaternary cooling and drying, Thuja gradually retreated southwards to form today's disjunctive distribution between East Asia and North America.

Species delimitation, genetic diversity and population historical dynamics of Cycas diannanensis (Cycadaceae) occurring sympatrically in the Red River region of China

Delimitating species boundaries could be of critical importance when evaluating the species' evolving process and providing guidelines for conservation genetics. Here, species delimitation was carried out on three endemic and endangered Cycas species with resembling morphology and overlapped distribution range along the Red River (Yuanjiang) in China: Cycas diananensis Z. T. Guan et G. D. Tao, Cycas parvula S. L. Yang and Cycas multiovula D. Y. Wang. A total of 137 individuals from 15 populations were genotyped by using three chloroplastic (psbA-trnH, atpI-atpH, and trnL-rps4) and two single copy nuclear (RPB1 and SmHP) DNA sequences. Basing on the carefully morphological comparison and cladistic haplotype aggregation (CHA) analysis, we propose all the populations as one species, with the rest two incorporated into C. diannanensis. Genetic diversity and structure analysis of the conflated C. diannanensis revealed this species possessed a relative lower genetic diversity than estimates of other Cycas species. The higher genetic diversity among populations and relative lower genetic diversity within populations, as well as obvious genetic differentiation among populations inferred from chloroplastic DNA (cpDNA) suggested a recent genetic loss within this protected species. Additionally, a clear genetic structure of C. diannanensis corresponding with geography was detected based on cpDNA, dividing its population ranges into "Yuanjiang-Nanhun" basin and "Ejia-Jiepai" basin groups. Demographical history analyses based on combined cpDNA and one nuclear DNA (nDNA) SmHP both showed the population size of C. diannanensis began to decrease in Quaternary glaciation with no subsequent expansion, while another nDNA RPB1 revealed a more recent sudden expansion after long-term population size contraction, suggesting its probable bottleneck events in history. Our findings offer grounded views for clarifying species boundaries of C. diannanensis when determining the conservation objectives. For operational guidelines, the downstream populations which occupy high and peculiar haplotypes should be given prior in-situ conservation. In addition, ex-situ conservation and reintroduction measures for decades of generations are supplemented for improving the population size and genetic diversity of the endemic and endangered species.

Ancient nuclear plastid DNA in the yew family (taxaceae)

Plastid-to-nucleus DNA transfer provides a rich genetic resource to the complexity of plant nuclear genome architecture. To date, the evolutionary route of nuclear plastid DNA (nupt) remain unknown in conifers. We have sequenced the complete plastomes of two yews, Amentotaxus formosana and Taxus mairei (Taxaceae of coniferales). Our comparative genomic analyses recovered an evolutionary scenario for plastomic reorganization from ancestral to extant plastomes in the three sampled Taxaceae genera, Amentotaxus, Cephalotaxus, and Taxus. Specific primers were designed to amplify nonsyntenic regions between ancestral and extant plastomes, and 12.6 kb of nupts were identified based on phylogenetic analyses. These nupts have significantly accumulated GC-to-AT mutations, reflecting a nuclear mutational environment shaped by spontaneous deamination of 5-methylcytosin. The ancestral initial codon of rps8 is retained in the T. nupts, but its corresponding extant codon is mutated and requires C-to-U RNA-editing. These findings suggest that nupts can help recover scenarios of the nucleotide mutation process. We show that the Taxaceae nupts we retrieved may have been retained because the Cretaceous and they carry information of both ancestral genomic organization and nucleotide composition, which offer clues for understanding the plastome evolution in conifers.

Molecular phylogenetics and evolutionary history of sect. Quinquefoliae (Pinus): implications for Northern Hemisphere biogeography

Climatic changes and tectonic events in the Cenozoic have greatly influenced the evolution and geographic distribution of the temperate flora. Such consequences should be most evident in plant groups that are ancient, widespread, and diverse. As one of the most widespread genera of trees, Pinus provides a good model for investigating the history of species diversification and biogeographic disjunction in the Northern Hemisphere. In this study, we reconstructed the phylogeny and investigated the evolutionary and biogeographic history of sect. Quinquefoliae (Pinus), a species-rich lineage disjunctly distributed in Asia, Europe and North America, based on complete taxon sampling and by using nine DNA fragments from chloroplast (cp), mitochondrial (mt) and nuclear genomes. The monophyly of the three subsections, Krempfianae, Gerardianae, and Strobus, is well-supported by cpDNA and nuclear gene phylogenies. However, neither subsect. Gerardianae nor subsect. Strobus forms a monophyletic group in the mtDNA phylogeny, in which sect. Quinquefoliae was divided into two major clades, one consisting of the North American and northeastern Asian species as well as the European P. peuce of subsect. Strobus, and the other comprising the remaining Eurasian species belonging to three subsections. The significant topological incongruence among the gene trees, in conjunction with divergence time estimation and ancestral area reconstruction, indicates that both ancient and relatively recent introgressive hybridization events occurred in the evolution of sect. Quinquefoliae, particularly in northeastern Asia and northwestern North America. In addition, the phylogenetic analysis suggests that the species of subsect. Strobus from subtropical eastern Asia and neighboring areas may have a single origin, although species non-monophyly is very widespread in the nuclear gene trees. Moreover, our study seems to support a Tethyan origin of sect. Quinquefoliae given the distributions and phylogenetic positions of subsects. Krempfianae and Gerardianae, and also highlights the importance of active mountain buildings and climatic changes during the Late Neogene in shaping the species diversity and geographic distribution of Pinus.

Development and evolution of the female gametophyte and fertilization process in Welwitschia mirabilis (Welwitschiaceae)

PREMISE OF THE STUDY

The female gametophyte of Welwitschia has long been viewed as highly divergent from other members of the Gnetales and, indeed, all other seed plants. However, the formation of female gametes and the process of fertilization have never been observed.

METHODS

Standard histological techniques were applied to study gametophyte development and the fertilization process in Welwitschia.

KEY RESULTS

In Welwitschia, fertilization events occur when pollen tubes with binucleate sperm cells grow down through the nucellus and encounter prothallial tubes, free nuclear tubular extensions of the micropylar end of the female gametophyte that grow up through the nucellus. Entry of a binucleate sperm cell into a vacuolate prothallial tube appears to stimulate the rapid coagulation of cytoplasm around a single female nucleus, which differentiates into an egg cell. One sperm nucleus enters the female gamete, while the second sperm nucleus remains outside and ultimately degenerates. Only a single fertilization event occurs per mating pair of pollen tube and prothallial tube.

CONCLUSIONS

Welwitschia lacks the gnetalean pattern of regular double fertilization, as found in Ephedra and Gnetum, involving sperm from a single pollen tube to yield two zygotes. Moreover, an analysis of character evolution indicates that the female gametophyte of Welwitschia is highly apomorphic both among seed plants, and specifically within Gnetales, but also shares several key synapomorphies with its sister taxon Gnetum. Finally, the biological role of prothallial tubes in Welwitschia is examined from the perspectives of gamete competition and kin conflict.

Phylogeny and divergence times of gymnosperms inferred from single-copy nuclear genes

Phylogenetic reconstruction is fundamental to study evolutionary biology and historical biogeography. However, there was not a molecular phylogeny of gymnosperms represented by extensive sampling at the genus level, and most published phylogenies of this group were constructed based on cytoplasmic DNA markers and/or the multi-copy nuclear ribosomal DNA. In this study, we use LFY and NLY, two single-copy nuclear genes that originated from an ancient gene duplication in the ancestor of seed plants, to reconstruct the phylogeny and estimate divergence times of gymnosperms based on a complete sampling of extant genera. The results indicate that the combined LFY and NLY coding sequences can resolve interfamilial relationships of gymnosperms and intergeneric relationships of most families. Moreover, the addition of intron sequences can improve the resolution in Podocarpaceae but not in cycads, although divergence times of the cycad genera are similar to or longer than those of the Podocarpaceae genera. Our study strongly supports cycads as the basal-most lineage of gymnosperms rather than sister to Ginkgoaceae, and a sister relationship between Podocarpaceae and Araucariaceae and between Cephalotaxaceae-Taxaceae and Cupressaceae. In addition, intergeneric relationships of some families that were controversial, and the relationships between Taxaceae and Cephalotaxaceae and between conifers and Gnetales are discussed based on the nuclear gene evidence. The molecular dating analysis suggests that drastic extinctions occurred in the early evolution of gymnosperms, and extant coniferous genera in the Northern Hemisphere are older than those in the Southern Hemisphere on average. This study provides an evolutionary framework for future studies on gymnosperms.

Mitochondrial DNA capture and divergence in Pinus provide new insights into the evolution of the genus

The evolution of the mitochondrial (mt) genome is far from being fully understood. Systematic investigations into the modes of inheritance, rates and patterns of recombination, nucleotide substitution, and structural changes in the mt genome are still lacking in many groups of plants. In this study, we sequenced >11kbp mtDNA segments from multiple accessions of 36 pine species to characterize the evolutionary patterns of mtDNA in the genus Pinus. We found extremely low substitution rates and complex repetitive sequences scattered across different genome regions, as well as chimeric structures that were probably generated by multiple intergenomic recombinations. The mtDNA-based phylogeny of the genus differed from that based on chloroplast and nuclear DNA in the placement of several groups of species. Such discordances suggest a series of mtDNA capture events during past range shifts of the pine species and that both vertical and horizontal inheritance are implicated in the evolution of mtDNA in Pinus. MtDNA dating revealed that most extant lineages of the genus originated during Oligocene-Miocene radiation and subgenus Strobus diversified earlier than subgenus Pinus. Our findings illustrate a reticular evolutionary pathway for the mt genome through capture and recombination in the genus Pinus, and provide new insights into the evolution of the genus.