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

Phylogenetics of the rust fungi (Pucciniales) of South Africa, with notes on their life histories and possible origins

South Africa has an indigenous rust (Pucciniales) funga of approximately 460 species. This funga was sampled with species from as many genera as possible. The nuclear ribosomal large subunit (28S) region was amplified from samples representing 110 indigenous species, as well as the small subunit (18S) region and the cytochrome c oxidase subunit 3 (CO3) in some cases, and these were used in phylogenetic analyses. One new species is described, 12 new combinations made, six names reinstated, and two life history connections made. The life histories of this funga were summarized; it is dominated by species with contracted life histories. The majority of species are autoecious, with a small proportion being heteroecious. Of the autoecious species, many will likely be homothallic with no spermagonia. A shortened life history with homothallism allows for a single basidiospore infection to initiate a local population buildup under the prevailing unpredictable climatic conditions. Suggestions are made as to the possible origin of this funga based on the development of the modern South African flora. It is postulated that the rusts of South Africa are of relatively recent origin, consisting of three groups. Firstly, there is an African tropical element with members of the Mikronegerineae (Hemileia), the Sphaerophragmiaceae (Puccorchidium, Sphaerophragmium), and certain Uredinineae (Stomatisora). Their immediate ancestors likely occurred in the tropical forests of Africa during the Paleogene. Secondly, there is a pantropical element including the Raveneliaceae (e.g., Diorchidium, Maravalia, Ravenelia sensu lato, Uropyxis). This likely diversified during the Neogene, when the mimosoids became the dominant trees of the developing savannas. Thirdly, the Pucciniaceae invaded Africa as this continent pushed northward closing the Tethys Sea. They diversified with the development of the savannas as these become the dominant habitat in most of Africa, and are by far the largest component of the South African rust funga.

SCEP1 and SCEP2 are two new components of the synaptonemal complex central element

The synaptonemal complex (SC) is a proteinaceous structure that forms between homologous chromosomes during meiosis prophase. The SC is widely conserved across species, but its structure and roles during meiotic recombination are still debated. While the SC central region is made up of transverse filaments and central element proteins in mammals and fungi, few central element proteins have been identified in other species. Here we report the identification of two coiled-coil proteins, SCEP1 and SCEP2, that form a complex and localize at the centre of the Arabidopsis thaliana SC. In scep1 and scep2 mutants, chromosomes are aligned but not synapsed (the ZYP1 transverse filament protein is not loaded), crossovers are increased compared with the wild type, interference is lost and heterochiasmy is strongly reduced. We thus report the identification of two plant SC central elements, and homologues of these are found in all major angiosperm clades.

Functional Characterization and Cyclization Mechanism of a Diterpene Synthase Catalyzing the Skeleton Formation of Cephalotane-Type Diterpenoids

CsCTS, a new diterpene synthase from Cephalotaxus sinensis responsible for forming cephalotene, the core skeleton of cephalotane-type diterpenoids with a highly rigid 6/6/5/7 tetracyclic ring system, was functionally characterized. The stepwise cyclization mechanism is proposed mainly based on structural investigation of its derailment products, and further demonstrated through isotopic labeling experiments and density functional theory calculations. Homology modeling and molecular dynamics simulation combined with site-directed mutagenesis revealed the critical amino acid residues for the unique carbocation-driven cascade cyclization mechanism of CsCTS. Altogether, this study reports the discovery of the diterpene synthase that catalyzes the first committed step of cephalotane-type diterpenoid biosynthesis and delineates its cyclization mechanism, laying the foundation to decipher and artificially construct the complete biosynthetic pathway of this type diterpenoids.

A new fossil Acmopyle with accessory transfusion tissue and potential reproductive buds: Direct evidence for ever-wet rainforests in Eocene Patagonia

PREMISE

Acmopyle (Podocarpaceae) comprises two extant species from Oceania that are physiologically restricted to ever-wet rainforests, a confirmed fossil record based on leaf adpressions and cuticles in Australia since the Paleocene, and a few uncertain reports from New Zealand, Antarctica, and South America. We investigated fossil specimens with Acmopyle affinities from the early Eocene Laguna del Hunco site in Patagonia, Argentina.

METHODS

We studied 42 adpression leafy-shoot fossils and included them in a total evidence phylogenetic analysis.

RESULTS

Acmopyle grayae sp. nov. is based on heterophyllous leafy shoots with three distinct leaf types. Among these, bilaterally flattened leaves uniquely preserve subparallel, linear features that we interpret as accessory transfusion tissue (ATT, an extra-venous water-conducting tissue). Some apical morphologies of A. grayae shoots are compatible with the early stages of ovuliferous cone development. Our phylogenetic analysis recovers the new species in a polytomy with the two extant Acmopyle species. We report several types of insect-herbivory damage. We also transfer Acmopyle engelhardti from the middle Eocene Río Pichileufú flora to Dacrycarpus engelhardti comb. nov.

CONCLUSIONS

We confirm the biogeographically significant presence of the endangered West Pacific genus Acmopyle in Eocene Patagonia. Acmopyle is one of the most drought-intolerant genera in Podocarpaceae, possibly due to the high collapse risk of the ATT, and thus the new fossil species provides physiological evidence for the presence of an ever-wet rainforest environment at Laguna del Hunco during the Early Eocene Climatic Optimum.

An integrative view on the systematic position of the cupressophyte Cephalotaxus

We made an in-depth review of historical studies of the cupressophyte conifer genus Cephalotaxus Siebold & Zucc. with an emphasis on its systematic position. We suggest that the systematic position of the genus is better understood using an integrative approach, so the evolution of phenetic characters is discussed within the context of recent phylogenomics. We propose that the genus should be classified as a separate family Cephalotaxaceae belonging to the clade consisting of Cupressaceae, Cephalotaxaceae, and Taxaceae; the family Cephalotaxaceae is sister to the Taxaceae but not nested within the Taxaceae and is characterized by a unique set of characters including morphology, anatomy, embryology, and chemistry. The family Cephalotaxaceae shows transitional characters between the Cupressaceae and the Taxaceae; the family possesses female cones with a primary cone axis bearing 5-8 pairs of decussate bracts, which is similar to the typical female cones of the Cupressaceae, on the one hand, and may have given rise to the reduced female cone of the Taxaceae with one terminal ovule partially or completely enclosed in a fleshy aril. In parallel, the compound male cone of the Cephalotaxaceae evolved into the seemingly "simple" male cones of the Taxaceae by means of reduction, elimination, and fusion.

Diversity, Distribution, Systematics and Conservation Status of Podocarpaceae

Among conifer families, Podocarpaceae is the second largest, with amazing diversity and functional traits, and it is the dominant Southern Hemisphere conifer family. However, comprehensive studies on diversity, distribution, systematic and ecophysiological aspects of the Podocarpaceae are sparse. We aim to outline and evaluate the current and past diversity, distribution, systematics, ecophysiological adaptations, endemism, and conservation status of podocarps. We analyzed data on the diversity and distribution of living and extinct macrofossil taxa and combined it with genetic data to reconstruct an updated phylogeny and understand historical biogeography. Podocarpaceae today contains 20 genera and approximately 219 taxa (201 species, 2 subspecies, 14 varieties and 2 hybrids) placed in three clades, plus a paraphyletic group/grade of four distinct genera. Macrofossil records show the presence of more than 100 podocarp taxa globally, dominantly from the Eocene-Miocene. Australasia (New Caledonia, Tasmania, New Zealand, and Malesia) is the hotspot of living podocarps diversity. Podocarps also show remarkable adaptations from broad to scale leaves, fleshy seed cones, animal dispersal, shrubs to large trees, from lowland to alpine regions and rheophyte to a parasite (including the only parasitic gymnosperm-Parasitaxus) and a complex pattern of seed and leaf functional trait evolution.

Phylotranscriptomics and evolution of key genes for terpene biosynthesis in Pinaceae

Pinaceae is the largest family of conifers, dominating forest ecosystems and serving as the backbone of northern, temperate and mountain forests. The terpenoid metabolism of conifers is responsive to pests, diseases, and environmental stress. Determining the phylogeny and evolution of terpene synthase genes in Pinaceae may shed light on early adaptive evolution. We used different inference methods and datasets to reconstruct the Pinaceae phylogeny based on our assembled transcriptomes. We identified the final species tree of Pinaceae by comparing and summarizing different phylogenetic trees. The genes encoding terpene synthase (TPS) and cytochrome P450 proteins in Pinaceae showed a trend of expansion compared with those in Cycas. Gene family analysis revealed that the number of TPS genes decreased while the number of P450 genes increased in loblolly pine. Expression profiles showed that TPSs and P450s were mainly expressed in leaf buds and needles, which may be the result of long-term evolution to protect these two vulnerable tissues. Our research provides insights into the phylogeny and evolution of terpene synthase genes in Pinaceae and offers some useful references for the investigation of terpenoids in conifers.

Five Species of Taxus Karyotype Based on Oligo-FISH for 5S rDNA and (AG3T3)3

As a relict plant, Taxus is used in a variety of medicinal ingredients, for instance to treat a variety of cancers. Taxus plants are difficult to distinguish from one another due to their similar morphology; indeed, some species of Taxus cytogenetic data still are unclear. Oligo-FISH can rapidly and efficiently provide insight into the genetic composition and karyotype. This is important for understanding the organization and evolution of chromosomes in Taxus species. We analysed five Taxus species using two oligonucleotide probes. (AG₃T₃)₃ signals were distributed at the chromosome ends and the centromere of five species of Taxus. The 5S rDNA signal was displayed on two chromosomes of five species of Taxus. In addition to Taxus wallichiana var. mairei, 5S rDNA signals were found proximal in the remaining four species, which signals a difference in its location. The karyotype formula of Taxus wallichiana was 2n = 2x = 24m, its karyotype asymmetry index was 55.56%, and its arm ratio was 3.0087. Taxus × media's karyotype formula was 2n = 2x = 24m, its karyotype asymmetry index was 55.09%, and its arm ratio was 3.4198. The karyotype formula of Taxus yunnanensis was 2n = 2x = 24m, its karyotype asymmetry index was 55.56%, and its arm ratio was 2.6402. The karyotype formula of Taxus cuspidate was 2n = 2x = 24m, its karyotype asymmetry index was 54.67%, its arm ratio was 3.0135, and two chromosomes exhibited the 5S rDNA signal. The karyotype formula of T. wallichiana var. mairei was 2n= 2x = 22m + 2sm, its karyotype asymmetry index was 54.33%, and its arm ratio was 2.8716. Our results provide the karyotype analysis and physical genetic map of five species of Taxus, which contributes to providing molecular cytogenetics data for Taxus.

Comprehensive collection of genes and comparative analysis of full-length transcriptome sequences from Japanese larch (Larix kaempferi) and Kuril larch (Larix gmelinii var. japonica)

BACKGROUND

Japanese larch (Larix kaempferi) is an economically important deciduous conifer species that grows in cool-temperate forests and is endemic to Japan. Kuril larch (L. gmelinii var. japonica) is a variety of Dahurian larch that is naturally distributed in the Kuril Islands and Sakhalin. The hybrid larch (L. gmelinii var. japonica × L. kaempferi) exhibits heterosis, which manifests as rapid juvenile growth and high resistance to vole grazing. Since these superior characteristics have been valued by forestry managers, the hybrid larch is one of the most important plantation species in Hokkaido. To accelerate molecular breeding in these species, we collected and compared full-length cDNA isoforms (Iso-Seq) and RNA-Seq short-read, and merged them to construct candidate gene as reference for both Larix species. To validate the results, candidate protein-coding genes (ORFs) related to some flowering signal-related genes ​were screened from the reference sequences, and the phylogenetic relationship with closely related species was elucidated.

RESULTS

Using the isoform sequencing of PacBio RS ll and the de novo assembly of RNA-Seq short-read sequences, we identified 50,690 and 38,684 ORFs in Japanese larch and Kuril larch, respectively. BUSCO completeness values were 90.5% and 92.1% in the Japanese and Kuril larches, respectively. After comparing the collected ORFs from the two larch species, a total of 19,813 clusters, comprising 22,571 Japanese larch ORFs and 22,667 Kuril larch ORFs, were contained in the intersection of the Venn diagram. In addition, we screened several ORFs related to flowering signals (SUPPRESSER OF OVEREXPRESSION OF CO1: SOC1, LEAFY: LFY, FLOWERING Locus T: FT, CONSTANCE: CO) from both reference sequences, and very similar found in other species.

CONCLUSIONS

The collected ORFs will be useful as reference sequences for molecular breeding of Japanese and Kuril larches, and also for clarifying the evolution of the conifer genome and investigating functional genomics.

A Comprehensive Evolutionary Study of Chloroplast RNA Editing in Gymnosperms: A Novel Type of G-to-A RNA Editing Is Common in Gymnosperms

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

Phylogeny and evolution of Cupressaceae: Updates on intergeneric relationships and new insights on ancient intergeneric hybridization

After the merger of the former Taxodiaceae and Cupressaceae s.s., currently the conifer family Cupressaceae (sensu lato) comprises seven subfamilies and 32 genera, most of which are important components of temperate and mountainous forests. With the exception of a recently published genus-level phylogeny of gymnosperms inferred from sequence analysis of 790 orthologs, previous phylogenetic studies of Cupressaceae were based mainly on morphological characters or a few molecular markers, and did not completely resolve the intergeneric relationships. In this study, we reconstructed a robust and well-resolved phylogeny of Cupressaceae represented by all 32 genera, using 1944 genes (Orthogroups) generated from transcriptome sequencing. Reticulate evolution analyses detected a possible ancient hybridization that occurred between ancestors of two subclades of Cupressoideae, including Microbiota-Platycladus-Tetraclinis (MPT) and Juniperus-Cupressus-Hesperocyparis-Callitropsis-Xanthocyparis (JCHCX), although both concatenation and coalescent trees are highly supported. Moreover, divergence time estimation and ancestral area reconstruction indicate that Cupressaceae very likely originated in Asia in the Triassic, and geographic isolation caused by continental separation drove the vicariant evolution of the two subfamilies Cupressoideae and Callitroideae in the northern and southern hemispheres, respectively. Evolutionary analyses of some morphological characters suggest that helically arranged linear-acicular leaves and imbricate bract-scale complexes represent ancestral states, and the shift from linear-acicular leaves to scale-like leaves was associated with the shift from helical to decussate arrangement. Our study sheds new light on phylogeny and evolutionary history of Cupressaceae, and strongly suggests that both dichotomous phylogenetic and reticulate evolution analyses be conducted in phylogenomic studies.

Recent advances on phylogenomics of gymnosperms and a new classification

Living gymnosperms comprise four major groups: cycads, Ginkgo, conifers, and gnetophytes. Relationships among/within these lineages have not been fully resolved. Next generation sequencing has made available a large number of sequences, including both plastomes and single-copy nuclear genes, for reconstruction of solid phylogenetic trees. Recent advances in gymnosperm phylogenomic studies have updated our knowledge of gymnosperm systematics. Here, we review major advances of gymnosperm phylogeny over the past 10 years and propose an updated classification of extant gymnosperms. This new classification includes three classes (Cycadopsida, Ginkgoopsida, and Pinopsida), five subclasses (Cycadidae, Ginkgoidae, Cupressidae, Pinidae, and Gnetidae), eight orders (Cycadales, Ginkgoales, Araucariales, Cupressales, Pinales, Ephedrales, Gnetales, and Welwitschiales), 13 families, and 86 genera. We also described six new tribes including Acmopyleae Y. Yang, Austrocedreae Y. Yang, Chamaecyparideae Y. Yang, Microcachrydeae Y. Yang, Papuacedreae Y. Yang, and Prumnopityeae Y. Yang, and made 27 new combinations in the genus Sabina.

Structural organization of the spongy mesophyll

Many plant leaves have two layers of photosynthetic tissue: the palisade and spongy mesophyll. Whereas palisade mesophyll consists of tightly packed columnar cells, the structure of spongy mesophyll is not well characterized and often treated as a random assemblage of irregularly shaped cells. Using micro-computed tomography imaging, topological analysis, and a comparative physiological framework, we examined the structure of the spongy mesophyll in 40 species from 30 genera with laminar leaves and reticulate venation. A spectrum of spongy mesophyll diversity encompassed two dominant phenotypes: first, an ordered, honeycomblike tissue structure that emerged from the spatial coordination of multilobed cells, conforming to the physical principles of Euler's law; and second, a less-ordered, isotropic network of cells. Phenotypic variation was associated with transitions in cell size, cell packing density, mesophyll surface-area-to-volume ratio, vein density, and maximum photosynthetic rate. These results show that simple principles may govern the organization and scaling of the spongy mesophyll in many plants and demonstrate the presence of structural patterns associated with leaf function. This improved understanding of mesophyll anatomy provides new opportunities for spatially explicit analyses of leaf development, physiology, and biomechanics.

Host Shift Speciation of the Ectomycorrhizal Genus Suillus (Suillineae, Boletales) and Biogeographic Comparison With Its Host Pinaceae

Suillus is a genus of ectomycorrhizal fungi associated almost exclusively with Pinaceae. Lack of sample collections in East Asia and unresolved basal phylogenetic relationships of the genus are the major obstacles for better understanding the Suillus evolution. A resolved phylogeny of Suillus representing global diversity was achieved by sequencing multiple nuclear ribosomal and protein coding genes and extensive samples collected in East Asia. Fungal fossils are extremely rare, and the Eocene ectomycorrhizal symbiosis (ECM) fossil of Pinus root has been widely used for calibration. This study explored an alternative calibration scenario of the ECM fossil for controversy. Ancestral host associations of Suillus were estimated by maximum likelihood and Bayesian Markov chain Monte Carlo (MCMC) analyses, inferred from current host information from root tips and field observation. Host shift speciation explains the diversification of Suillus major clades. The three basal subgenera of Suillus were inferred to be associated with Larix, and diverged in early Eocene or Upper Cretaceous. In the early Oligocene or Paleocene, subgenus Suillus diverged and switched host to Pinus subgenus Strobus, and then switched to subgenus Pinus four times. Suillus subgenus Douglasii switched host from Larix to Pseudotsuga in Oligocene or Eocene. Increased species diversity occurred in subgenus Suillus after it switched host to Pinus but no associated speciation rate shifts were detected. Ancestral biogeographic distributions of Suillus and Pinaceae were estimated under the Dispersal Extinction Cladogenesis (DEC) model. Ancestral distribution patterns of Suillus and Pinaceae are related but generally discordant. Dispersals between Eurasia and North America explain the prevalence of disjunct Suillus taxa.

Tree functional traits, forest biomass, and tree species diversity interact with site properties to drive forest soil carbon

Forests constitute important ecosystems in the global carbon cycle. However, how trees and environmental conditions interact to determine the amount of organic carbon stored in forest soils is a hotly debated subject. In particular, how tree species influence soil organic carbon (SOC) remains unclear. Based on a global compilation of data, we show that functional traits of trees and forest standing biomass explain half of the local variability in forest SOC. The effects of functional traits on SOC depended on the climatic and soil conditions with the strongest effect observed under boreal climate and on acidic, poor, coarse-textured soils. Mixing tree species in forests also favours the storage of SOC, provided that a biomass over-yielding occurs in mixed forests. We propose that the forest carbon sink can be optimised by (i) increasing standing biomass, (ii) increasing forest species richness, and (iii) choosing forest composition based on tree functional traits according to the local conditions.

Forests constitute important ecosystems in the global carbon cycle. This study investigates how tree species influence soil organic carbon using a global dataset, showing the importance of tree functional traits and forest standing biomass to optimise forest carbon sink.

Species Diversification of the Coniferous Pathogenic Fungal Genus Coniferiporia (Hymenochaetales, Basidiomycota) in Association with Its Biogeography and Host Plants

Coniferiporia, belonging to Hymenochaetaceae and now segregated from Phellinidium, is a wood-inhabiting fungal genus with three species, each having a specific geographic distribution and a strong host specificity as a forest pathogen of coniferous trees. In this study, the species diversity of Coniferiporia is further clarified with the aid of a wider sampling and multilocus-based phylogenetic analysis, which reveals a new species Coniferiporia uzbekistanensis. The molecular clock and ancestral geographic origin analyses indicate that the ancestor of Coniferiporia emerged in one of the Pinaceae and Cupressaceae, then jumped to the other plant family originated in eastern Eurasia 17.01 million years ago (Mya; 95% highest posterior density: 9.46 to 25.86 Mya), and later extended its distribution to western North America, Central Asia, and eastern Europe. Coniferiporia sulphurascens speciated on Pinaceae in eastern Eurasia 8.78 Mya (9.46 to 25.86 Mya) and then extended its distribution to western North America and eastern Europe. Coniferiporia qilianensis and C. uzbekistanensis speciated on Juniperus przewalskii in eastern Eurasia 3.67 Mya (0.36 to 8.02 Mya) and on Juniperus polycarpos in Central Asia 4.35 Mya (0.94 to 8.37 Mya), respectively. The speciation event of Coniferiporia weirii occurred 4.45 Mya (0.77 to 9.33 Mya) right after the emergence of its host, the endemic Cupressaceae species Thuja plicata, and soon after, this fungus evolved to also inhabit another endemic Cupressaceae species Calocedrus decurrens. In summary, this study for the first time unambiguously clarified and timed the adaptive evolutionary event of Coniferiporia in association with its biogeography and host plants.

Flavonoids and Phenolic Compounds From the Parasitic Gymnosperm Parasitaxus usta Endemic to New Caledonia

Parasitaxus usta (Podocarpaceae) is the only parasitic gymnosperm and endemic to New Caledonia. In this survey, 11 flavonoids and 6 phenolic compounds were isolated from the aerial parts. As for flavonoids, six flavones, apigenin 7- O-glucoside (1), luteolin (2), luteolin 7- O-glucoside (3), chrysoeriol (4), chrysoeriol 7- O-glucoside (5) and tricetin 3′- O-glucoside (6), one C-glycosylflavone, orientin (7), one flavonol, quercetin 3- O-glucoside (8), one anthocyanin, cyanidin 3- O-glucoside (9), and two biflavones, isoginkgetin (10) and agathisflavone (11) were identified by UV, liquid chromatograph–mass spectra (LC-MS), acid hydrolysis, NMR and/or HPLC comparisons with authentic samples. On the other hand, six phenolic compounds were identified as 5- O- E- p-coumaroyl quinic acid (12), 5- O- Z- p-coumaroyl quinic acid (13), 5- O- E- p-coumaroyl quinic acid methyl ester (14), 5- O- Z- p-coumaroyl quinic acid methyl ester (15), E-caffeic acid methyl ester 3- O-β-glucopyranoside (16), and Z-caffeic acid methyl ester 3- O-β-glucopyranoside (17) by UV, LC-MS and NMR. Chemical components of P usta were reported in this survey for the first time. Their chemical characters were chemotaxonomically compared with those of other Podocarpaceae species.

Comparative Analyses of 3,654 Plastid Genomes Unravel Insights Into Evolutionary Dynamics and Phylogenetic Discordance of Green Plants

The plastid organelle is essential for many vital cellular processes and the growth and development of plants. The availability of a large number of complete plastid genomes could be effectively utilized to understand the evolution of the plastid genomes and phylogenetic relationships among plants. We comprehensively analyzed the plastid genomes of Viridiplantae comprising 3,654 taxa from 298 families and 111 orders and compared the genomic organizations in their plastid genomic DNA among major clades, which include gene gain/loss, gene copy number, GC content, and gene blocks. We discovered that some important genes that exhibit similar functions likely formed gene blocks, such as the psb family presumably showing co-occurrence and forming gene blocks in Viridiplantae. The inverted repeats (IRs) in plastid genomes have doubled in size across land plants, and their GC content is substantially higher than non-IR genes. By employing three different data sets [all nucleotide positions (nt123), only the first and second codon positions (nt12), and amino acids (AA)], our phylogenomic analyses revealed Chlorokybales + Mesostigmatales as the earliest-branching lineage of streptophytes. Hornworts, mosses, and liverworts forming a monophylum were identified as the sister lineage of tracheophytes. Based on nt12 and AA data sets, monocots, Chloranthales and magnoliids are successive sister lineages to the eudicots + Ceratophyllales clade. The comprehensive taxon sampling and analysis of different data sets from plastid genomes recovered well-supported relationships of green plants, thereby contributing to resolving some long-standing uncertainties in the plant phylogeny.

Gymnosperm Resprouting—A Review

Gymnosperms are generally regarded as poor resprouters, especially when compared to angiosperms and particularly following major disturbance. However, is it this clear-cut? This review investigates two main aspects of gymnosperm resprouting: (i) various papers have provided exceptions to the above generalization—how frequent are these exceptions and are there any taxonomic trends?; and (ii) assuming gymnosperms are poor resprouters are there any anatomical or physiological reasons why this is the case? Five of six non-coniferous gymnosperm genera and 24 of 80 conifer genera had at least one species with a well-developed resprouting capability. This was a wider range than would be expected from the usual observation ‘gymnosperms are poor resprouters’. All conifer families had at least three resprouting genera, except the monospecific Sciadopityaceae. Apart from the aboveground stem, buds were also recorded arising from more specialised structures (e.g., lignotubers, tubers, burls and underground stems). In some larger genera it appeared that only a relatively small proportion of species were resprouters and often only when young. The poor resprouting performance of mature plants may stem from a high proportion of apparently ‘blank’ leaf axils. Axillary meristems have been recorded in a wide range of conifer species, but they often did not form an apical dome, leaf primordia or vascular connections. Buds or meristems that did form often abscised at an early stage. While this review has confirmed that conifers do not resprout to the same degree as angiosperms, it was found that a wide diversity of gymnosperm genera can recover vegetatively after substantial disturbance. Further structural studies are needed, especially of: (i) apparently blank leaf axils and the initial development of axillary meristems; (ii) specialised regeneration structures; and (iii) why high variability can occur in the resprouting capacity within species of a single genus and within genera of the same family.

New insights into the phylogeny and evolution of Podocarpaceae inferred from transcriptomic data

Phylogenies of an increasing number of taxa have been resolved with the development of phylogenomics. However, the intergeneric relationships of Podocarpaceae, the second largest family of conifers comprising 19 genera and approximately 187 species mainly distributed in the Southern Hemisphere, have not been well disentangled in previous studies, even when genome-scale data sets were used. Here we used 993 nuclear orthologous groups (OGs) and 54 chloroplast OGs (genes), which were generated from 47 transcriptomes of Podocarpaceae and its sister group Araucariaceae, to reconstruct the phylogeny of Podocarpaceae. Our study completely resolved the intergeneric relationships of Podocarpaceae represented by all extant genera and revealed that topological conflicts among phylogenetic trees could be attributed to synonymous substitutions. Moreover, we found that two morphological traits, fleshy seed cones and flattened leaves, might be important for Podocarpaceae to adapt to angiosperm-dominated forests and thus could have promoted its species diversification. In addition, our results indicate that Podocarpaceae originated in Gondwana in the late Triassic and both vicariance and dispersal have contributed to its current biogeographic patterns. Our study provides the first robust transcriptome-based phylogeny of Podocarpaceae, an evolutionary framework important for future studies of this family.

The dynamic history of gymnosperm plastomes: Insights from structural characterization, comparative analysis, phylogenomics, and time divergence

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution.

Evolution and structural variations in chloroplast tRNAs in gymnosperms

Background

Chloroplast transfer RNAs (tRNAs) can participate in various vital processes. Gymnosperms have important ecological and economic value, and they are the dominant species in forest ecosystems in the Northern Hemisphere. However, the evolution and structural changes in chloroplast tRNAs in gymnosperms remain largely unclear.

Results

In this study, we determined the nucleotide evolution, phylogenetic relationships, and structural variations in 1779 chloroplast tRNAs in gymnosperms. The numbers and types of tRNA genes present in the chloroplast genomes of different gymnosperms did not differ greatly, where the average number of tRNAs was 33 and the frequencies of occurrence for various types of tRNAs were generally consistent. Nearly half of the anticodons were absent. Molecular sequence variation analysis identified the conserved secondary structures of tRNAs. About a quarter of the tRNA genes were found to contain precoded 3′ CCA tails. A few tRNAs have undergone novel structural changes that are closely related to their minimum free energy, and these structural changes affect the stability of the tRNAs. Phylogenetic analysis showed that tRNAs have evolved from multiple common ancestors. The transition rate was higher than the transversion rate in gymnosperm chloroplast tRNAs. More loss events than duplication events have occurred in gymnosperm chloroplast tRNAs during their evolutionary process.

Conclusions

These findings provide novel insights into the molecular evolution and biological characteristics of chloroplast tRNAs in gymnosperms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08058-3.

Comparative Repeat Profiling of Two Closely Related Conifers (Larix decidua and Larix kaempferi) Reveals High Genome Similarity With Only Few Fast-Evolving Satellite DNAs

In eukaryotic genomes, cycles of repeat expansion and removal lead to large-scale genomic changes and propel organisms forward in evolution. However, in conifers, active repeat removal is thought to be limited, leading to expansions of their genomes, mostly exceeding 10 giga base pairs. As a result, conifer genomes are largely littered with fragmented and decayed repeats. Here, we aim to investigate how the repeat landscapes of two related conifers have diverged, given the conifers' accumulative genome evolution mode. For this, we applied low-coverage sequencing and read clustering to the genomes of European and Japanese larch, Larix decidua (Lamb.) Carrière and Larix kaempferi (Mill.), that arose from a common ancestor, but are now geographically isolated. We found that both Larix species harbored largely similar repeat landscapes, especially regarding the transposable element content. To pin down possible genomic changes, we focused on the repeat class with the fastest sequence turnover: satellite DNAs (satDNAs). Using comparative bioinformatics, Southern, and fluorescent in situ hybridization, we reveal the satDNAs' organizational patterns, their abundances, and chromosomal locations. Four out of the five identified satDNAs are widespread in the Larix genus, with two even present in the more distantly related Pseudotsuga and Abies genera. Unexpectedly, the EulaSat3 family was restricted to L. decidua and absent from L. kaempferi, indicating its evolutionarily young age. Taken together, our results exemplify how the accumulative genome evolution of conifers may limit the overall divergence of repeats after speciation, producing only few repeat-induced genomic novelties.

Molecular evolution and SSRs analysis based on the chloroplast genome of Callitropsis funebris

Chloroplast genome sequences have been used to understand evolutionary events and to infer efficiently phylogenetic relationships. Callitropsis funebris (Cupressaceae) is an endemic species in China. Its phylogenetic position is controversial due to morphological characters similar to those of Cupressus, Callitropsis, and Chamaecyparis. This study used next-generation sequencing technology to sequence the complete chloroplast genome of Ca. funebris and then constructed the phylogenetic relationship between Ca. funebris and its related species based on a variety of data sets and methods. Simple sequence repeats (SSRs) and adaptive evolution analysis were also conducted. Our results showed that the monophyletic branch consisting of Ca. funebris and Cupressus tonkinensis is a sister to Cupressus, while Callitropsis is not monophyletic; Ca. nootkatensis and Ca. vietnamensis are nested in turn at the base of the monophyletic group Hesperocyparis. The statistical results of SSRs supported the closest relationship between Ca. funebris and Cupressus. By performing adaptive evolution analysis under the phylogenetic background of Cupressales, the Branch model detected three genes and the Site model detected 10 genes under positive selection; and the Branch-Site model uncovered that rpoA has experienced positive selection in the Ca. funebries branch. Molecular analysis from the chloroplast genome highly supported that Ca. funebris is at the base of Cupressus. Of note, SSR features were found to be able to shed some light on phylogenetic relationships. In short, this chloroplast genomic study has provided new insights into the phylogeny of Ca. funebris and revealed multiple chloroplast genes possibly undergoing adaptive evolution.

Plastome phylogenomics of Cephalotaxus (Cephalotaxaceae) and allied genera

BACKGROUND AND AIMS

Cephalotaxus is a paleo-endemic genus in East Asia that consists of about 7-9 conifer species. Despite its great economic and ecological importance, the relationships between Cephalotaxus and related genera, as well as the interspecific relationships within Cephalotaxus, have long been controversial, resulting in contrasting taxonomic proposals in delimitation of Cephalotaxaceae and Taxaceae. Based on plastome data, this study aims to reconstruct a robust phylogeny to infer the systematic placement and the evolutionary history of Cephalotaxus.

METHODS

A total of 11 plastomes, representing all species currently recognized in Cephalotaxus and two Torreya species, were sequenced and assembled. Combining these with previously published plastomes, we reconstructed a phylogeny of Cephalotaxaceae and Taxaceae with nearly full taxonomic sampling. Under a phylogenetic framework and molecular dating, the diversification history of Cephalotaxus and allied genera was explored.

KEY RESULTS

Phylogenetic analyses of 81 plastid protein-coding genes recovered robust relationships between Cephalotaxus and related genera, as well as providing a well-supported resolution of interspecific relationships within Cephalotaxus, Taxus, Torreya and Amentotaxus. Divergence time estimation indicated that most extant species of these genera are relatively young, although fossil and other molecular evidence consistently show that these genera are ancient plant lineages.

CONCLUSIONS

Our results justify the taxonomic proposal that recognizes Cephalotaxaceae as a monotypic family, and contribute to a clear-cut delineation between Cephalotaxaceae and Taxaceae. Given that extant species of Cephalotaxus are derived from recent divergence events associated with the establishment of monsoonal climates in East Asia and Pleistocene climatic fluctuations, they are not evolutionary relics.

Sareomycetes: more diverse than meets the eye

Since its resurrection, the resinicolous discomycete genus Sarea has been accepted as containing two species, one with black apothecia and pycnidia, and one with orange. We investigate this hypothesis using three ribosomal (nuITS, nuLSU, mtSSU) regions from and morphological examination of 70 specimens collected primarily in Europe and North America. The results of our analyses support separation of the traditional Sarea difformis s.lat. and Sarea resinae s.lat. into two distinct genera, Sarea and Zythia. Sarea as circumscribed is shown to conservatively comprise three phylospecies, with one corresponding to Sarea difformis s.str. and two, morphologically indistinguishable, corresponding to the newly combined Sarea coeloplata. Zythia is provisionally maintained as monotypic, containing only a genetically and morphologically variable Z. resinae. The new genus Atrozythia is erected for the new species A. klamathica. Arthrographis lignicola is placed in this genus on molecular grounds, expanding the concept of Sareomycetes by inclusion of a previously unknown type of asexual morph. Dating analyses using additional marker regions indicate the emergence of the Sareomycetes was roughly concurrent with the diversification of the genus Pinus, suggesting that this group of fungi emerged to exploit the newly-available resinous ecological niche supplied by Pinus or another, extinct group of conifers. Our phylogeographic studies also permitted us to study the introductions of these fungi to areas where they are not native, including Antarctica, Cape Verde, and New Zealand and are consistent with historical hypotheses of introduction.

The RGF/GLV/CLEL Family of Short Peptides Evolved Through Lineage-Specific Losses and Diversification and Yet Conserves Its Signaling Role Between Vascular Plants and Bryophytes

Short secreted plant peptides act as key signaling molecules and control a plethora of developmental and physiological processes. The ROOT GROWTH FACTOR (RGF)/GOLVEN (GLV)/CLE-Like (CLEL) family of peptides was discovered to be involved in root development in Arabidopsis thaliana. In contrast to active research efforts, which have been revealing receptors and downstream signaling components, little attention has been paid to evolutionary processes that shaped the RGF signaling system as we know it in angiosperms today. As a first step toward understanding how RGF signaling emerged and evolved, this study aimed to elucidate the phylogenetic distribution and functional conservation of RGF-like sequences. Using publicly available, genome and transcriptome data, RGF-like sequences were searched in 27 liverworts, 22 mosses, 8 hornworts, 23 lycophytes, 23 ferns, 38 gymnosperms, and 8 angiosperms. This led to the identification of more than four hundreds of RGF-like sequences in all major extant land plant lineages except for hornworts. Sequence comparisons within and between taxonomic groups identified lineage-specific characters. Notably, one of the two major RGF subgroups, represented by A. thaliana RGF6/GLV1/CLEL6, was found only in vascular plants. This subgroup, therefore, likely emerged in a common ancestor of vascular plants after its divergence from bryophytes. In bryophytes, our results infer independent losses of RGF-like sequences in mosses and hornworts. On the other hand, a single, highly similar RGF-like sequence is conserved in liverworts, including Marchantia polymorpha, a genetically tractable model species. When constitutively expressed, the M. polymorpha RGF-like sequence (MpRGF) affected plant development and growth both in A. thaliana and M. polymorpha. This suggests that MpRGF can exert known RGF-like effects and that MpRGF is under transcriptional control so that its potent activities are precisely controlled. These data suggest that RGFs are conserved as signaling molecules in both vascular plants and bryophytes and that lineage-specific diversification has increased sequence variations of RGFs. All together, our findings form a basis for further studies into RGF peptides and their receptors, which will contribute to our understandings of how peptide signaling pathways evolve.

Fruiting body form, not nutritional mode, is the major driver of diversification in mushroom-forming fungi

With ∼36,000 described species, Agaricomycetes are among the most successful groups of Fungi. Agaricomycetes display great diversity in fruiting body forms and nutritional modes. Most have pileate-stipitate fruiting bodies (with a cap and stalk), but the group also contains crust-like resupinate fungi, polypores, coral fungi, and gasteroid forms (e.g., puffballs and stinkhorns). Some Agaricomycetes enter into ectomycorrhizal symbioses with plants, while others are decayers (saprotrophs) or pathogens. We constructed a megaphylogeny of 8,400 species and used it to test the following five hypotheses regarding the evolution of morphological and ecological traits in Agaricomycetes and their impact on diversification: 1) resupinate forms are plesiomorphic, 2) pileate-stipitate forms promote diversification, 3) the evolution of gasteroid forms is irreversible, 4) the ectomycorrhizal (ECM) symbiosis promotes diversification, and 5) the evolution of ECM symbiosis is irreversible. The ancestor of Agaricomycetes was a saprotroph with a resupinate fruiting body. There have been 462 transitions in the examined morphologies, including 123 origins of gasteroid forms. Reversals of gasteroid forms are highly unlikely but cannot be rejected. Pileate-stipitate forms are correlated with elevated diversification rates, suggesting that this morphological trait is a key to the success of Agaricomycetes. ECM symbioses have evolved 36 times in Agaricomycetes, with several transformations to parasitism. Across the entire 8,400-species phylogeny, diversification rates of ectomycorrhizal lineages are no greater than those of saprotrophic lineages. However, some ECM lineages have elevated diversification rates compared to their non-ECM sister clades, suggesting that the evolution of symbioses may act as a key innovation at local phylogenetic scales.

The flattened and needlelike leaves of the pine family (Pinaceae) share a conserved genetic network for adaxial-abaxial polarity but have diverged for photosynthetic adaptation

Background

Leaves have highly diverse morphologies. However, with an evolutionary history of approximately 200 million years, leaves of the pine family are relatively monotonous and often collectively called “needles”, although they vary in length, width and cross-section shapes. It would be of great interest to determine whether Pinaceae leaves share similar morpho-physiological features and even consistent developmental and adaptive mechanisms.

Results

Based on a detailed morpho-anatomical study of leaves from all 11 Pinaceae genera, we particularly investigated the expression patterns of adaxial-abaxial polarity genes in two types of leaves (needlelike and flattened) and compared their photosynthetic capacities. We found that the two types of leaves share conserved spatial patterning of vasculatures and genetic networks for adaxial-abaxial polarity, although they display different anatomical structures in the mesophyll tissue differentiation and distribution direction. In addition, the species with needlelike leaves exhibited better photosynthetic capacity than the species with flattened leaves.

Conclusions

Our study provides the first evidence for the existence of a conserved genetic module controlling adaxial-abaxial polarity in the development of different Pinaceae leaves.

Beyond Wallace: a new lineage of Chrysorthenches (Lepidoptera: Yponomeutoidea: Glyphipterigidae) reveals a journey tracking its host-plants, Podocarpus (Pinopsida: Podocarpaceae)

A northward trans-Wallacean radiation is demonstrated for Chrysorthenches, a member of the Orthenches group. Here we review Chrysorthenches and allied genera resulting in a generic transfer of Diathryptica callibrya to Chrysorthenches and two new congeners: C. muraseaeSohn & Kobayashisp. nov. from Japan and C. smaragdinaSohnsp. nov. from Thailand. We review morphological characters of Chrysorthenches and allied genera, and find polyphyly of Diathryptica and the association of the Orthenches-group with Glyphipterigidae. These findings were supported in a maximum likelihood phylogeny of DNA barcodes from ten yponomeutoids. We analysed 30 morphological characters for 12 species of Chrysorthenches, plus one outgroup, via a cladistic approach. The resulting cladogram redefined two pre-existing Chrysorthenches species-groups and identified one novel lineage: the C. callibrya species-group. We review the host associations between Chrysorthenches and Podocarpaceae, based on mapping the working phylogenies. Our review suggests that ancestral Chrysorthenches colonized Podocarpus and later shifted to other podocarp genera. Biogeographical patterns of Chrysorthenches show that they evolved long after the Podocarpaceae radiation. Disjunctive trans-Wallacean distribution of the C. callibrya species-group is possibly related to the tracking of their host-plants and the complicated geological history of the island-arc system connecting Australia and East Asia.

The Mitogenome of Norway Spruce and a Reappraisal of Mitochondrial Recombination in Plants

Plant mitogenomes can be difficult to assemble because they are structurally dynamic and prone to intergenomic DNA transfers, leading to the unusual situation where an organelle genome is far outnumbered by its nuclear counterparts. As a result, comparative mitogenome studies are in their infancy and some key aspects of genome evolution are still known mainly from pregenomic, qualitative methods. To help address these limitations, we combined machine learning and in silico enrichment of mitochondrial-like long reads to assemble the bacterial-sized mitogenome of Norway spruce (Pinaceae: Picea abies). We conducted comparative analyses of repeat abundance, intergenomic transfers, substitution and rearrangement rates, and estimated repeat-by-repeat homologous recombination rates. Prompted by our discovery of highly recombinogenic small repeats in P. abies, we assessed the genomic support for the prevailing hypothesis that intramolecular recombination is predominantly driven by repeat length, with larger repeats facilitating DNA exchange more readily. Overall, we found mixed support for this view: Recombination dynamics were heterogeneous across vascular plants and highly active small repeats (ca. 200 bp) were present in about one-third of studied mitogenomes. As in previous studies, we did not observe any robust relationships among commonly studied genome attributes, but we identify variation in recombination rates as a underinvestigated source of plant mitogenome diversity.

The Origin and Evolution of Plastid Genome Downsizing in Southern Hemispheric Cypresses (Cupressaceae)

Plastome downsizing is rare in photosynthetic seed plants. However, a large-scale study of five cupressophyte families (conifers II) indicated that the plastomes of some Cupressaceous genera are notably reduced and compact. Here, we enriched taxon sampling in Cupressaceae by adding plastomes of ten previously unreported genera to determine the origin, evolution, and consequences of plastome reduction in this family. We discovered that plastome downsizing is specific to Callitroideae (a Southern Hemispheric subfamily). Their plastomes are the smallest, encode the fewest plastid genes, and contain the fewest GC-end codons among Cupressaceae. We show that repeated tRNA losses and shrinkage of intergenic spacers together contributed to the plastome downsizing in Callitroideae. Moreover, our absolute nucleotide substitution rate analyses suggest relaxed functional constraints in translation-related plastid genes (clpP, infA, rpl, and rps), but not in photosynthesis- or transcription-related ones, of Callitris (the most diverse genus in Callitroideae). We hypothesize that the small and low-GC plastomes of Callitroideae emerged ca. 112-75 million years ago as an adaptation to increased competition with angiosperms on the Gondwana supercontinent. Our findings highlight Callitroideae as another case of plastome downsizing in photosynthetic seed plant lineages.

Chemical Analysis of Pollen by FT-Raman and FTIR Spectroscopies

Pollen studies are important for the assessment of present and past environment, including biodiversity, sexual reproduction of plants and plant-pollinator interactions, monitoring of aeroallergens, and impact of climate and pollution on wild communities and cultivated crops. Although information on chemical composition of pollen is of importance in all of those research areas, pollen chemistry has been rarely measured due to complex and time-consuming analyses. Vibrational spectroscopies, coupled with multivariate data analysis, have shown great potential for rapid chemical characterization, identification and classification of pollen. This study, comprising 219 species from all principal taxa of seed plants, has demonstrated that high-quality Raman spectra of pollen can be obtained by Fourier transform (FT) Raman spectroscopy. In combination with Fourier transform infrared spectroscopy (FTIR), FT-Raman spectroscopy is obtaining comprehensive information on pollen chemistry. Presence of all the main biochemical constituents of pollen, such as proteins, lipids, carbohydrates, carotenoids and sporopollenins, have been identified and detected in the spectra, and the study shows approaches to measure relative and absolute content of these constituents. The results show that FT-Raman spectroscopy has clear advantage over standard dispersive Raman measurements, in particular for measurement of pollen samples with high pigment content. FT-Raman spectra are strongly biased toward chemical composition of pollen wall constituents, namely sporopollenins and pigments. This makes Raman spectra complementary to FTIR spectra, which over-represent chemical constituents of the grain interior, such as lipids and carbohydrates. The results show a large variability in pollen chemistry for families, genera and even congeneric species, revealing wide range of reproductive strategies, from storage of nutrients to variation in carotenoids and phenylpropanoids. The information on pollen's chemical patterns for major plant taxa should be of outstanding value for various studies in plant biology and ecology, including aerobiology, palaeoecology, forensics, community ecology, plant-pollinator interactions, and climate effects on plants.

A new macrofossil ephedroid plant with unusual bract morphology from the Lower Cretaceous Jiufotang Formation of northeastern China

BACKGROUND

The evolution of the Jehol Biota of western Liaoning in China includes three phases, initiation in the Dabeigou phase, radiation in the Yixian phase, and decline in the Jiufotang phase. Numerous ephedroid macrofossils were reported from the Lower Cretaceous Yixian Formation. However, so far none has been found in the younger Jiufotang Formation (ca. 120.3 Ma) of western Liaoning.

RESULTS

Here we report a new species Jianchangia verticillata gen. et sp. nov. with unusual morphology from the Lower Cretaceous of the Jiufotang Formation, Lamadong Village, Jianchang County, Liaoning. This species is the first record of gnetophytes from the Jiufotang Formation. It is similar to other ephedroid species from the Yixian Formation in possessing linear leaves with parallel veins, jointed shoots with swollen nodes and longitudinally furrowed internodes, and ovulate cones possessing two whorls of bracts enclosing two chlamydosperms, but differs from all known species by the ovulate cone having multiple fine linear verticillate bracts.

CONCLUSIONS

This study expands our knowledge about the diversity of early gnetophytes in the Lower Cretaceous, and demonstrates the lineage continuity of gnetophytes from the Yixian Formation to the younger Jiufotang Formation.

The Expansion and Diversification of Pentatricopeptide Repeat RNA-Editing Factors in Plants

The RNA-binding pentatricopeptide repeat (PPR) family comprises hundreds to thousands of genes in most plants, but only a few dozen in algae, indicating massive gene expansions during land plant evolution. The nature and timing of these expansions has not been well defined due to the sparse sequence data available from early-diverging land plant lineages. In this study, we exploit the comprehensive OneKP datasets of over 1000 transcriptomes from diverse plants and algae toward establishing a clear picture of the evolution of this massive gene family, focusing on the proteins typically associated with RNA editing, which show the most spectacular variation in numbers and domain composition across the plant kingdom. We characterize over 2 250 000 PPR motifs in over 400 000 proteins. In lycophytes, polypod ferns, and hornworts, nearly 10% of expressed protein-coding genes encode putative PPR editing factors, whereas they are absent from algae and complex-thalloid liverworts. We show that rather than a single expansion, most land plant lineages with high numbers of editing factors have continued to generate novel sequence diversity. We identify sequence variations that imply functional differences between PPR proteins in seed plants versus non-seed plants and variations we propose to be linked to seed-plant-specific editing co-factors. Finally, using the sequence variations across the datasets, we develop a structural model of the catalytic DYW domain associated with C-to-U editing and identify a clade of unique DYW variants that are strong candidates as U-to-C RNA-editing factors, given their phylogenetic distribution and sequence characteristics.

Target Nuclear and Off-Target Plastid Hybrid Enrichment Data Inform a Range of Evolutionary Depths in the Orchid Genus Epidendrum

Universal angiosperm enrichment probe sets designed to enrich hundreds of putatively orthologous nuclear single-copy loci are increasingly being applied to infer phylogenetic relationships of different lineages of angiosperms at a range of evolutionary depths. Studies applying such probe sets have focused on testing the universality and performance of the target nuclear loci, but they have not taken advantage of off-target data from other genome compartments generated alongside the nuclear loci. Here we do so to infer phylogenetic relationships in the orchid genus Epidendrum and closely related genera of subtribe Laeliinae. Our aims are to: 1) test the technical viability of applying the plant anchored hybrid enrichment (AHE) method (Angiosperm v.1 probe kit) to our focal group, 2) mine plastid protein coding genes from off-target reads; and 3) evaluate the performance of the target nuclear and off-target plastid loci in resolving and supporting phylogenetic relationships along a range of taxonomical depths. Phylogenetic relationships were inferred from the nuclear data set through coalescent summary and site-based methods, whereas plastid loci were analyzed in a concatenated partitioned matrix under maximum likelihood. The usefulness of target and flanking non-target nuclear regions and plastid loci was assessed through the estimation of their phylogenetic informativeness. Our study successfully applied the plant AHE probe kit to Epidendrum, supporting the universality of this kit in angiosperms. Moreover, it demonstrated the feasibility of mining plastome loci from off-target reads generated with the Angiosperm v.1 probe kit to obtain additional, uniparentally inherited sequence data at no extra sequencing cost. Our analyses detected some strongly supported incongruences between nuclear and plastid data sets at shallow divergences, an indication of potential lineage sorting, hybridization, or introgression events in the group. Lastly, we found that the per site phylogenetic informativeness of the ycf1 plastid gene surpasses that of all other plastid genes and several nuclear loci, making it an excellent candidate for assessing phylogenetic relationships at medium to low taxonomic levels in orchids.

The complete mitochondrial genome of Taxus cuspidata (Taxaceae): eight protein-coding genes have transferred to the nuclear genome

BACKGROUND

Gymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms.

RESULTS

Here, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species.

CONCLUSIONS

We have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

Comparative Genomic Analysis Reveals the Mechanism Driving the Diversification of Plastomic Structure in Taxaceae Species

Inverted repeat (IR) regions in the plastomes from land plants induce homologous recombination, generating isomeric plastomes. While the plastomes of Taxaceae species often lose one of the IR regions, considerable isomeric plastomes were created in Taxaceae species with a hitherto unclarified mechanism. To investigate the detailed mechanism underpinning the IR-independent genesis of plastomic diversity, we sequenced four Taxaceae plastomes, including Taxus cuspidata Siebold & Zuccarini, Taxus fauna Nan Li & R. R. Mill, and two individuals of Taxus wallichiana Zuccarini. Then we compared these structures with those of previously reported Taxaceae plastomes. Our analysis identified four distinct plastome forms that originated from the rearrangements of two IR-flanking inverted fragments. The presence of isomeric plastomes was then verified in T. cuspidata individuals. Both rearrangement analyses and phylogenetic results indicated that Taxaceae were separated into two clades, one including Taxus and Pseudotaxus and another formed by Amentotaxus and Torreya. Our reconstructed scenario suggests that the minimum number of inversion events required for the transformation of the plastome of Cephalotaxus oliveri Masters into the diversified Taxaceae plastomes ranged from three to six. To sum up, our study reveals a distinct pattern and the mechanism driving the structural diversification of Taxaceae plastomes, which will advance our understanding of the maintenance of plastomic diversity and complexity in conifers.

The complete chloroplast genome of Microcycas calocoma (Miq.) A. DC. (Zamiaceae, Cycadales) and evolution in Cycadales

Cycadales is an extant group of seed plants occurring in subtropical and tropical regions comprising putatively three families and 10 genera. At least one complete plastid genome sequence has been reported for all of the 10 genera except Microcycas, making it an ideal plant group to conduct comprehensive plastome comparisons at the genus level. This article reports for the first time the plastid genome of Microcycas calocoma. The plastid genome has a length of 165,688 bp with 134 annotated genes including 86 protein-coding genes, 47 non-coding RNA genes (39 tRNA and eight rRNA) and one pseudogene. Using global sequence variation analysis, the results showed that all cycad genomes share highly similar genomic profiles indicating significant slow evolution and little variation. However, identity matrices coinciding with the inverted repeat regions showed fewer similarities indicating that higher polymorphic events occur at those sites. Conserved non-coding regions also appear to be more divergent whereas variations in the exons were less discernible indicating that the latter comprises more conserved sequences. Phylogenetic analysis using 81 concatenated protein-coding genes of chloroplast (cp) genomes, obtained using maximum likelihood and Bayesian inference with high support values (>70% ML and = 1.0 BPP), confirms that Microcycas is closest to Zamia and forms a monophyletic clade with Ceratozamia and Stangeria. While Stangeria joined the Neotropical cycads Ceratozamia, Zamia and Microcyas, Bowenia grouped with the Southern Hemisphere cycads Encephalartos, Lepidozamia and Macrozamia. All Cycas species formed a distinct clade separated from the other genera. Dioon, on the other hand, was outlying from the rest of Zamiaceae encompassing two major clades—the Southern Hemisphere cycads and the Neotropical cycads. Analysis of the whole cp genomes in phylogeny also supports that the previously recognized family—Stangeriaceae—which contained Bowenia and Stangeria, is not monophyletic. Thus, the cp genome topology obtained in our study is congruent with other molecular phylogenies recognizing only a two-family classification (Cycadaceae and Zamiaceae) within extant Cycadales.

Evolutionary history of a relict conifer, Pseudotaxus chienii (Taxaceae), in south-east China during the late Neogene: old lineage, young populations

BACKGROUND AND AIMS

Many monotypic gymnosperm lineages in south-east China paradoxically remain in relict status despite long evolutionary histories and ample opportunities for allopatric speciation, but this paradox has received little attention and has yet to be resolved. Here, we address this issue by investigating the evolutionary history of a relict conifer, Pseudotaxus chienii (Taxaceae).

METHODS

DNA sequences from two chloroplast regions and 14 nuclear loci were obtained for 134 samples. The demographic history was inferred and the contribution of isolation by environment (IBE) in patterning genetic divergence was compared with that of isolation by distance (IBD).

KEY RESULTS

Three genetic clusters were identified. Approximate Bayesian computation analyses showed that the three clusters diverged in the late Pliocene (~3.68 Ma) and two admixture events were detected. Asymmetric gene flow and similar population divergence times (~ 3.74 Ma) were characterized using the isolation with migration model. Neither IBD nor IBE contributed significantly to genetic divergence, and the contribution of IBE was much smaller than that of IBD.

CONCLUSIONS

These results suggest that several monotypic relict gymnosperm lineages like P. chienii in south-east China did not remain in situ and undiversified for millions of years. On the contrary, they have been evolving and the extant populations have become established more recently, having insufficient time to speciate. Our findings provide a new perspective for understanding the formation and evolution of the relict gymnosperm flora of China as well as of the Sino-Japanese Flora.

Transcriptome analysis of two radiated Cycas species and the subsequent species delimitation of the Cycas taiwaniana complex

PREMISE

Cycas is an important gymnosperm genus, and the most diverse of all cycad genera. The C. taiwaniana complex of species are morphologically similar and difficult to distinguish due to a lack of genomic resources.

METHODS

We characterized the transcriptomes of two closely related and endangered Cycas species endemic to Hainan, China: C. hainanensis and C. changjiangensis. Three single-copy nuclear genes in the C. taiwaniana complex were sequenced based on these transcriptomes, enabling us to evaluate the species boundaries using the multispecies coalescent method implemented in the Bayesian Phylogenetics and Phylogeography program.

RESULTS

We obtained 68,184 and 81,561 unigenes for C. changjiangensis and C. hainanensis, respectively. We identified six positively selected genes that are mainly involved in stimulus responses, suggesting that environmental adaptation may have played an important role in the relatively recent divergence of these species. The similar K S distribution peaks at 1.0 observed for the paralogs in the two species indicate a common whole-genome duplication event. Our species delimitation analysis indicated that the C. taiwaniana complex consists of three distinct species, which correspond to the previously reported morphological differences.

DISCUSSION

Our study provides valuable genetic resources for Cycas species and guidance for the taxonomic treatment of the C. taiwaniana complex, as well as new insights into evolution of species within Cycas.

A Target Capture-Based Method to Estimate Ploidy From Herbarium Specimens

Whole genome duplication (WGD) events are common in many plant lineages, but the ploidy status and possible occurrence of intraspecific ploidy variation are unknown for most species. Standard methods for ploidy determination are chromosome counting and flow cytometry approaches. While flow cytometry approaches typically use fresh tissue, an increasing number of studies have shown that recently dried specimens can be used to yield ploidy data. Recent studies have started to explore whether high-throughput sequencing (HTS) data can be used to assess ploidy levels by analyzing allelic frequencies from single copy nuclear genes. Here, we compare different approaches using a range of yam (Dioscorea) tissues of varying ages, drying methods and quality, including herbarium tissue. Our aims were to: (1) explore the limits of flow cytometry in estimating ploidy level from dried samples, including herbarium vouchers collected between 1831 and 2011, and (2) optimize a HTS-based method to estimate ploidy by considering allelic frequencies from nuclear genes obtained using a target-capture method. We show that, although flow cytometry can be used to estimate ploidy levels from herbarium specimens collected up to fifteen years ago, success rate is low (5.9%). We validated our HTS-based estimates of ploidy using 260 genes by benchmarking with dried samples of species of known ploidy (Dioscorea alata, D. communis, and D. sylvatica). Subsequently, we successfully applied the method to the 85 herbarium samples analyzed with flow cytometry, and successfully provided results for 91.7% of them, comprising species across the phylogenetic tree of Dioscorea. We also explored the limits of using this HTS-based approach for identifying high ploidy levels in herbarium material and the effects of heterozygosity and sequence coverage. Overall, we demonstrated that ploidy diversity within and between species may be ascertained from historical collections, allowing the determination of polyploidization events from samples collected up to two centuries ago. This approach has the potential to provide insights into the drivers and dynamics of ploidy level changes during plant evolution and crop domestication.

Battistuzzi F, Kumar S. A Machine Learning Method for Detecting Autocorrelation of Evolutionary Rates in Large Phylogenies

New species arise from pre-existing species and inherit similar genomes and environments. This predicts greater similarity of the tempo of molecular evolution between direct ancestors and descendants, resulting in autocorrelation of evolutionary rates in the tree of life. Surprisingly, molecular sequence data have not confirmed this expectation, possibly because available methods lack the power to detect autocorrelated rates. Here, we present a machine learning method, CorrTest, to detect the presence of rate autocorrelation in large phylogenies. CorrTest is computationally efficient and performs better than the available state-of-the-art method. Application of CorrTest reveals extensive rate autocorrelation in DNA and amino acid sequence evolution of mammals, birds, insects, metazoans, plants, fungi, parasitic protozoans, and prokaryotes. Therefore, rate autocorrelation is a common phenomenon throughout the tree of life. These findings suggest concordance between molecular and nonmolecular evolutionary patterns, and they will foster unbiased and precise dating of the tree of life.

Ancient duplications and grass-specific transposition influenced the evolution of LEAFY transcription factor genes

The LFY transcription factor gene family are important in the promotion of cell proliferation and floral development. Understanding their evolution offers an insight into floral development in plant evolution. Though a promiscuous transition intermediate and a gene duplication event within the LFY family had been identified previously, the early evolutionary path of this family remained elusive. Here, we reconstructed the LFY family phylogeny using maximum-likelihood and Bayesian inference methods incorporating LFY genes from all major lineages of streptophytes. The well-resolved phylogeny unveiled a high-confidence duplication event before the functional divergence of types I and II LFY genes in the ancestry of liverworts, mosses and tracheophytes, supporting sub-functionalization of an ancestral promiscuous gene. The identification of promiscuous genes in Osmunda suggested promiscuous LFY genes experienced an ancient transient duplication. Genomic synteny comparisons demonstrated a deep genomic positional conservation of LFY genes and an ancestral lineage-specific transposition activity in grasses.

Araucaria angustifolia chloroplast genome sequence and its relation to other Araucariaceae

Araucaria angustifolia is endemic to southern Brazil. Known as Brazilian pine, A. angustifolia is the only native conifer species with economic and social relevance in this country. Due to massive exploitation, it has suffered a significant population decline and currently is classified as critically endangered. This encouraged the scientific community to investigate genetic features in Brazilian pine to increase resources for management and preservation. In this work, RNA-Seq data was used to determine the complete nucleotide sequence of the A. angustifolia chloroplast genome (cpDNA). The cpDNA is 146,203 bp in length and contains 122 genes, including 80 protein-coding genes, 5 ribosomal RNA genes, and 37 tRNA genes. Coding regions comprise 45.02%, 4.96% correspond to rRNAs and tRNAs, and 50.02% of the genome encompasses non-coding regions. Genes found in the inverted repeat (IR) are present as single copy, with exception of the rrn5 and trnI-CAU loci. The typical LSC, SSC, IRa and IRb organization reported in several land-plant groups is not present in A. angustifolia cpDNA. Phylogenetic analyses using Bayesian and Maximum Likelihood methods clustered A. angustifolia in the Araucariaceae family, with A. heterophylla and A. columnaris as congeneric species. The screening of A. angustifolia cpDNA reveled 100 SSRs, 14 of them corresponding to tetrapolymer loci.

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.

Reconstructing the ancestral flower of extant angiosperms: the 'war of the whorls' is heating up

The origin of the angiosperm flower is a long-standing problem of botany and evolutionary biology. One widely accepted milestone towards solving it is the reconstruction of the ancestral flower of extant angiosperms, here termed 'AFEA'. A recent approach employing novel methods gave results that were not anticipated. Most notably the reconstructed phyllotaxis of AFEA soon was criticized and sparked a heated debate in the literature. To better explain, clarify, and perhaps cool the debate, we first summarize the results of previous attempts to reconstruct AFEA and contrast them with the more recent, controversial prediction of its structure. We then outline the major arguments made by contrasting parties in the recent debate. Finally, we discuss two key topics, the molecular mechanism of phyllotaxis and the role of gene regulatory networks during flower development and evolution, that may help to clarify the issue in the intermediate future.