The Torreya grandis genome illuminates the origin and evolution of gymnosperm-specific sciadonic acid biosynthesis.

Torreya plants produce dry fruits with assorted functions. Here, we report the 19-Gb chromosome-level genome assembly of T. grandis. The genome is shaped by ancient whole-genome duplications and recurrent LTR retrotransposon bursts. Comparative genomic analyses reveal key genes involved in reproductive organ development, cell wall biosynthesis and seed storage. Two genes encoding a C18 Δ9-elongase and a C20 Δ5-desaturase are identified to be responsible for sciadonic acid biosynthesis and both are present in diverse plant lineages except angiosperms. We demonstrate that the histidine-rich boxes of the Δ5-desaturase are crucial for its catalytic activity. Methylome analysis reveals that methylation valleys of the T. grandis seed genome harbor genes associated with important seed activities, including cell wall and lipid biosynthesis. Moreover, seed development is accompanied by DNA methylation changes that possibly fuel energy production. This study provides important genomic resources and elucidates the evolutionary mechanism of sciadonic acid biosynthesis in land plants.

Associations between picocyanobacterial ecotypes and cyanophage host genes across ocean basins and depth.

Background: Cyanophages, viruses that infect cyanobacteria, are globally abundant in the ocean's euphotic zone and are a potentially important cause of mortality for marine picocyanobacteria. Viral host genes are thought to increase viral fitness by either increasing numbers of genes for synthesizing nucleotides for virus replication, or by mitigating direct stresses imposed by the environment. The encoding of host genes in viral genomes through horizontal gene transfer is a form of evolution that links viruses, hosts, and the environment. We previously examined depth profiles of the proportion of cyanophage containing various host genes in the Eastern Tropical North Pacific Oxygen Deficient Zone (ODZ) and at the subtropical North Atlantic (BATS). However, cyanophage host genes have not been previously examined in environmental depth profiles across the oceans. Methodology: We examined geographical and depth distributions of picocyanobacterial ecotypes, cyanophage, and their viral-host genes across ocean basins including the North Atlantic, Mediterranean Sea, North Pacific, South Pacific, and Eastern Tropical North and South Pacific ODZs using phylogenetic metagenomic read placement. We determined the proportion of myo and podo-cyanophage containing a range of host genes by comparing to cyanophage single copy core gene terminase (terL). With this large dataset (22 stations), network analysis identified statistical links between 12 of the 14 cyanophage host genes examined here with their picocyanobacteria host ecotypes. Results: Picyanobacterial ecotypes, and the composition and proportion of cyanophage host genes, shifted dramatically and predictably with depth. For most of the cyanophage host genes examined here, we found that the composition of host ecotypes predicted the proportion of viral host genes harbored by the cyanophage community. Terminase is too conserved to illuminate the myo-cyanophage community structure. Cyanophage cobS was present in almost all myo-cyanophage and did not vary in proportion with depth. We used the composition of cobS phylotypes to track changes in myo-cyanophage composition. Conclusions: Picocyanobacteria ecotypes shift with changes in light, temperature, and oxygen and many common cyanophage host genes shift concomitantly. However, cyanophage phosphate transporter gene pstS appeared to instead vary with ocean basin and was most abundant in low phosphate regions. Abundances of cyanophage host genes related to nutrient acquisition may diverge from host ecotype constraints as the same host can live in varying nutrient concentrations. Myo-cyanophage community in the anoxic ODZ had reduced diversity. By comparison to the oxic ocean, we can see which cyanophage host genes are especially abundant (nirA, nirC, and purS) or not abundant (myo psbA) in ODZs, highlighting both the stability of conditions in the ODZ and the importance of nitrite as an N source to ODZ endemic LLV Prochlorococcus.

The earliest pollen-loaded insects from the Lower Permian of Russia.

Recent fossil discoveries suggest that the coevolution of insect pollinators and gymnosperms started long before the appearance of flowering plants. One of the keys to understanding the origins of pollination relationships is fossil insects with gymnosperm pollen attached to the body surface. Such fossils are exceedingly rare to find, especially from the Palaeozoic, a time when ambers with insect inclusions were absent. Here, we report compression fossils of Early Permian tillyardembiid insects (Polyneoptera) preserved with pollen on their heads, thoraces, legs and abdomens. This is the earliest finding of pollen-bearing insects, predating the previous oldest record from the Middle Jurassic by ca 120 Ma. Judging by the pollen composition, tillyardembiids visited a narrow range of host plants, including Rufloriaceae (Cordaitales). While it is impossible to say for certain whether tillyardembiids as pollen consumers contributed to pollination, a trophic specialization of this kind could be considered an evolutionary precursor of pollination mutualism.

Evolutionary analysis of KED-rich proteins in plants.

During the course of evolution, organisms have developed genetic mechanisms in response to various environmental stresses including wounding from mechanical damage or herbivory-caused injury. A previous study of wounding response in the plant tobacco identified a unique wound-induced gene, aptly named KED due to its coding for a protein that has an unusually high content of amino acids lysine (K), glutamic acid (E) and aspartic acid (D). However, by far little is known about this intriguing gene. In this study, we investigated the evolutionary aspects of the KED-rich coding genes. We found that a consistent pattern of wound-induced KED gene expression is maintained across representative species of angiosperm and gymnosperm. KED genes can be identified in species from all groups of land plants (Embryophyta). All the KED proteins from vascular plants (Tracheophyta) including angiosperm, gymnosperm, fern and lycophyte share a conserved 19-amino acid domain near the C-terminus, whereas bryophytes (moss, liverwort and hornwort) possess KED-rich, multi-direct-repeat sequences that are distinct from the vascular plant KEDs. We detected KED-rich sequences in Charophyta species but not in Chlorophyta wherever genome sequences are available. Our studies suggest diverse and complex evolution pathways for land plant KED genes. Vascular plant KEDs exhibit high evolutionary conservation, implicating their shared function in response to wounding stress. The extraordinary enrichment of amino acids K, E and D in these groups of distinct and widely distributed proteins may reflect the structural and functional requirement for these three residues during some 600 million years of land plant evolution.

Complete mitochondrial genome of Thuja sutchuenensis and its implications on evolutionary analysis of complex mitogenome architecture in Cupressaceae.

BACKGROUND: The complex physical structure and abundant repeat sequences make it difficult to assemble the mitogenomes of seed plants, especially gymnosperms. Only approximately 33 mitogenomes of gymnosperms have been reported. However, as the most widely distributed and the second largest family among gymnosperms, Cupressaceae has only six assembled mitogenomes, including five draft mitogenomes and one complete mitogenome, which has greatly hindered the understanding of mitogenome evolution within this large family, even gymnosperms. RESULTS: In this study, we assembled and validated the complete mitogenome of Thuja sutchuenensis, with a size of 2.4 Mb. Multiple sequence units constituted its complex structure, which can be reduced to three linear contigs and one small circular contig. The analysis of repeat sequences indicated that the numbers of simple sequence repeats increased during the evolutionary history of gymnosperms, and the mitogenome of Thuja sutchuenensis harboured abundant extra-long repeats (more than 5 kb). Additionally, the longest repeat sequence identified in these seven gymnosperms also came from the mitogenome of Thuja sutchuenensis, with a length of up to 47 kb. The analysis of colinear blocks and gene clusters both revealed that the orders of mitochondrial genes within gymnosperms was not conserved. The comparative analysis showed that only four tRNAs were shared by seven gymnosperms, namely, trnD-GUC, trnE-UUC, trnI-CAU and trnY-GUA. Furthermore, four genes have undergone potential positive selection in most gymnosperm species, namely, atp8, ccmB, mttB and sdh4. CONCLUSION: We successfully assembled the second complete mitogenome within Cupressaceae and verified that it consisted of multiple sequence units. Our study also indicated that abundant long repeats may contribute to the generation of the complex conformation of the mitogenome of Thuja sutchuenensis. The investigation of Thuja sutchuenensis's mitogenome in our study provides new insight into further understanding the complex mitogenome architecture within gymnosperms.

Soil organic matter molecular composition with long-term detrital alterations is controlled by site-specific forest properties.

Forest ecosystems are important global soil carbon (C) reservoirs, but their capacity to sequester C is susceptible to climate change factors that alter the quantity and quality of C inputs. To better understand forest soil C responses to altered C inputs, we integrated three molecular composition published data sets of soil organic matter (SOM) and soil microbial communities for mineral soils after 20 years of detrital input and removal treatments in two deciduous forests: Bousson Forest (BF), Harvard Forest (HF), and a coniferous forest: H.J. Andrews Forest (HJA). Soil C turnover times were estimated from radiocarbon measurements and compared with the molecular-level data (based on nuclear magnetic resonance and specific analysis of plant- and microbial-derived compounds) to better understand how ecosystem properties control soil C biogeochemistry and dynamics. Doubled aboveground litter additions did not increase soil C for any of the forests studied likely due to long-term soil priming. The degree of SOM decomposition was higher for bacteria-dominated sites with higher nitrogen (N) availability while lower for the N-poor coniferous forest. Litter exclusions significantly decreased soil C, increased SOM decomposition state, and led to the adaptation of the microbial communities to changes in available substrates. Finally, although aboveground litter determined soil C dynamics and its molecular composition in the coniferous forest (HJA), belowground litter appeared to be more influential in broadleaf deciduous forests (BH and HF). This synthesis demonstrates that inherent ecosystem properties regulate how soil C dynamics change with litter manipulations at the molecular-level. Across the forests studied, 20 years of litter additions did not enhance soil C content, whereas litter reductions negatively impacted soil C concentrations. These results indicate that soil C biogeochemistry at these temperate forests is highly sensitive to changes in litter deposition, which are a product of environmental change drivers.

Transcriptome analysis provides insights into light condition effect on paclitaxel biosynthesis in yew saplings.

BACKGROUND: Taxus is a rare gymnosperm plant that is the sole producer of the anticancer drug paclitaxel. The growth and development of Taxus is affected by environmental factors such as light. However, little is known about how light conditions affect growth and metabolic processes, especially paclitaxel biosynthesis. RESULTS: In this study, we applied three different light conditions to Taxus chinensis young saplings and investigated the physiological response and gene expression. Our observations showed that exposure to high light led to oxidative stress, caused photoinhibition, and damaged the photosynthetic systems in T. chinensis. The paclitaxel content in T. chinensis leaves was significantly decreased after the light intensity increased. Transcriptomic analysis revealed that numerous genes involved in paclitaxel biosynthesis and phenylpropanoid metabolic pathways were downregulated under high light. We also analyzed the expression of JA signaling genes, bHLH, MYB, AP2/ERF transcription factors, and the CYP450 families that are potentially related to paclitaxel biosynthesis. We found that several CYP450s, MYB and AP2/ERF genes were induced by high light. These genes may play an important role in tolerance to excessive light or heat stress in T. chinensis. CONCLUSIONS: Our study elucidates the molecular mechanism of the effects of light conditions on the growth and development of T. chinensis and paclitaxel biosynthesis, thus facilitating the artificial regeneration of Taxus and enhancing paclitaxel production.

Study of the Pollen Grain Metabolome under Deposition of Nitrogen and Phosphorus in Taxus baccata L. and Juniperus communis L.

Nitrogen plays an important role in both quantitative and qualitative aspects of plant reproduction, including pollen grain compounds and seed production. Recent studies have pointed out that pollen grains produced by male plants of T. baccata and J. communis subjected to a long period of fertilizer supplementation have lower in vitro germination ability and higher nitrogen content. To gain molecular insights into these observations, we conducted GC-MS analysis of both species to characterize the metabolomes of dry, mature pollen grains, which allowed for the identification and quantification of more than 200 metabolites. The results demonstrated that fertilizer supplementation impacts the relative content of 14 metabolites in J. communis (9 downregulated and 5 upregulated) and 21 in T. baccata (6 downregulated and 15 upregulated). Although plants showed little similarity in patterns, in metabolite profiles, both up and down fold-changes were observed. This is the first report on the gymnosperm pollen grain metabolomic profile and changes induced by long-term nitrogen and phosphorus supplementation. Pollen grains produced by fertilizer-supplemented male individuals had significantly lower relative content of linolenic acid, 5,6-dihydrouracil, maltotriose, galactonic acid, D-xylulose, and glycerol-α-phosphate but higher content of sorbitol, glucosamine, and 1,5-anhydro-D-glucitol as well as n-acetyl-d-hexosamine, dimethyl phthalate, glycine, galactose-6-phosphate, D-fructose-6-phosphate, pyroglutamic acid, and 3-(3-hydroxyphenyl)-3-hydroxypropionic acid. Thus, in pollen grain samples earlier shown to have different germination abilities, the presence of different metabolites indicates a significant environmental impact on the quality of gymnosperm pollen grains.

Three new species of Calepitrimerus (Acari: Eriophyidae) from China.

The genus Calepitrimerus (Acari: Eriophyoidea) comprise more than 90 extant species. Most species infest angiosperms, while only five were reported from gymnosperms. Herein, we describe and illustrate three new species of the genus Calepitrimerus from China: two from angiosperms, C. huayanlingense sp. nov. on Ficus erecta (Moraceae), C. cryptocaryaer sp. nov. on Cryptocarya concinna (Lauraceae), and one from a gymnosperm, C. liuanense sp. nov. on Cryptomeria japonica (Cupressaceae). All new species were found vagrant on lower surfaces of their hosts leaves, albeit inducing no apparent symptom. A list of Calepitrimerus species in China is tabulated.

Two new species of Pharaxonothinae beetles (Coleoptera: Erotylidae) inhabiting cones of the cycad Ceratozamia santillanii Prez-Farr. amp; Vovides (Cycadales: Zamiaceae) in Mexico.

Two beetles, Ceratophila (Vovidesa) chipjonesi, new species, and Pharaxonotha perezi, new species, (Coleoptera: Erotylidae: Pharaxonothinae) are described from the male cones of the New World cycad Ceratozamia santillanii Prez-Farr. Vovides (Cycadales: Zamiaceae). Morphological analysis and keys are provided for distinguishing these from related beetles inhabiting Ceratozamia in Mexico.

Phylotranscriptomics reveal the spatio-temporal distribution and morphological evolution of Macrozamia, an Australian endemic genus of Cycadales.

BACKGROUND AND AIMS: Cycads are regarded as an ancient lineage of living seed plants, and hold important clues to understand the early evolutionary trends of seed plants. The molecular phylogeny and spatio-temporal diversification of one of the species-rich genera of cycads, Macrozamia, have not been well reconstructed. METHODS: We analysed a transcriptome dataset of 4740 single-copy nuclear genes (SCGs) of 39 Macrozamia species and two outgroup taxa. Based on concatenated (maximum parsimony, maximum likelihood) and multispecies coalescent analyses, we first establish a well-resolved phylogenetic tree of Macrozamia. To identify cyto-nuclear incongruence, the plastid protein coding genes (PCGs) from transcriptome data are extracted using the software HybPiper. Furthermore, we explore the biogeographical history of the genus and shed light on the pattern of floristic exchange between three distinct areas of Australia. Six key diagnostic characters are traced on the phylogenetic framework using two comparative methods, and infra-generic classification is investigated. KEY RESULTS: The tree topologies of concatenated and multi-species coalescent analyses of SCGs are mostly congruent with a few conflicting nodes, while those from plastid PCGs show poorly supported relationships. The genus contains three major clades that correspond to their distinct distributional areas in Australia. The crown group of Macrozamia is estimated to around 11.80 Ma, with a major expansion in the last 5-6 Myr. Six morphological characters show homoplasy, and the traditional phenetic sectional division of the genus is inconsistent with this current phylogeny. CONCLUSIONS: This first detailed phylogenetic investigation of Macrozamia demonstrates promising prospects of SCGs in resolving phylogenetic relationships within cycads. Our study suggests that Macrozamia, once widely distributed in Australia, underwent major extinctions because of fluctuating climatic conditions such as cooling and mesic biome disappearance in the past. The current close placement of morphologically distinct species in the phylogenetic tree may be related to neotenic events that occurred in the genus.

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.

Mega-sized pericentromeric blocks of simple telomeric repeats and their variants reveal patterns of chromosome evolution in ancient Cycadales genomes.

Simple telomeric repeats composed of six to seven iterating nucleotide units are important sequences typically found at the ends of chromosomes. Here we analyzed their abundance and homogeneity in 42 gymnosperm (29 newly sequenced), 29 angiosperm (one newly sequenced), and eight bryophytes using bioinformatics, conventional cytogenetic and molecular biology approaches to explore their diversity across land plants. We found more than 10 000-fold variation in the amounts of telomeric repeats among the investigated taxa. Repeat abundance was positively correlated with increasing intragenomic sequence heterogeneity and occurrence at non-telomeric positions, but there was no correlation with genome size. The highest abundance/heterogeneity was found in the gymnosperm genus Cycas (Cycadaceae), in which megabase-sized blocks of telomeric repeats (i.e., billions of copies) were identified. Fluorescent in situ hybridization experiments using variant-specific probes revealed canonical Arabidopsis-type telomeric TTTAGGG repeats at chromosome ends, while pericentromeric blocks comprised at least four major telomeric variants with decreasing abundance: TTTAGGG>TTCAGGG >TTTAAGG>TTCAAGG. Such a diversity of repeats was not found in the sister cycad family Zamiaceae or in any other species analyzed. Using immunocytochemistry, we showed that the pericentromeric blocks of telomeric repeats overlapped with histone H3 serine 10 phosphorylation signals. We show that species of Cycas have amplified their telomeric repeats in centromeric and telomeric positions on telocentric chromosomes to extraordinary high levels. The ancestral chromosome number reconstruction suggests their occurrence is unlikely to be the product of ancient Robertsonian chromosome fusions. We speculate as to how the observed chromosome dynamics may be associated with the diversification of cycads.

The Dynamic Changes of Alternative Electron Flows upon Transition from Low to High Light in the Fern Cyrtomium fortune and the Gymnosperm Nageia nagi.

In photosynthetic organisms except angiosperms, an alternative electron sink that is mediated by flavodiiron proteins (FLVs) plays the major role in preventing PSI photoinhibition while cyclic electron flow (CEF) is also essential for normal growth under fluctuating light. However, the dynamic changes of FLVs and CEF has not yet been well clarified. In this study, we measured the P700 signal, chlorophyll fluorescence, and electrochromic shift spectra in the fern Cyrtomium fortune and the gymnosperm Nageia nagi. We found that both species could not build up a sufficient proton gradient (∆pH) within the first 30 s after light abruptly increased. During this period, FLVs-dependent alternative electron flow was functional to avoid PSI over-reduction. This functional time of FLVs was much longer than previously thought. By comparison, CEF was highly activated within the first 10 s after transition from low to high light, which favored energy balancing rather than the regulation of a PSI redox state. When FLVs were inactivated during steady-state photosynthesis, CEF was re-activated to favor photoprotection and to sustain photosynthesis. These results provide new insight into how FLVs and CEF interact to regulate photosynthesis in non-angiosperms.

Projections of leaf turgor loss point shifts under future climate change scenarios.

Predicting the consequences of climate change is of utmost importance to mitigate impacts on vulnerable ecosystems; plant hydraulic traits are particularly useful proxies for predicting functional disruptions potentially occurring in the near future. This study assessed the current and future regional patterns of leaf water potential at turgor loss point (Ψtlp ) by measuring and projecting the Ψtlp of 166 vascular plant species (159 angiosperms and 7 gymnosperms) across a large climatic range spanning from alpine to Mediterranean areas in NE Italy. For angiosperms, random forest models predicted a consistent shift toward more negative values in low-elevation areas, whereas for gymnosperms the pattern was more variable, particularly in the alpine sector (i.e., Alps and Prealps). Simulations were also developed to evaluate the number of threatened species under two Ψtlp plasticity scenarios (low vs. high plasticity), and it was found that in the worst-case scenario approximately 72% of the angiosperm species and 68% of gymnosperms within a location were at risk to exceed their physiological plasticity. The different responses to climate change by specific clades might produce reassembly in natural communities, undermining the resilience of natural ecosystems to climate change.

Independent recruitment of glycosyltransferase family 61 members for xylan substitutions in conifers.

MAIN CONCLUSION: Several pine members of the gymnosperm-specific GT61 clades were demonstrated to be arabinosyltransferases and xylosyltransferases catalyzing the transfer of 2-O-Araf, 3-O-Araf and 2-O-Xyl side chains onto xylooligomer acceptors, indicating their possible involvement in Araf and Xyl substitutions of xylan in pine. Xylan in conifer wood is substituted at O-2 with methylglucuronic acid (MeGlcA) as well as at O-3 with arabinofuranose (Araf), which differs from xylan in dicot wood that is typically decorated with MeGlcA but not Araf. Currently, glycosyltransferases responsible for conifer xylan arabinosylation have not been identified. Here, we investigated the roles of pine glycosyltransferase family 61 (GT61) members in xylan substitutions. Biochemical characterization of four pine wood-associated GT61 members showed that they exhibited three distinct glycosyltransferase activities involved in xylan substitutions. Two of them catalyzed the addition of 2-O-α-Araf or 3-O-α-Araf side chains onto xylooligomer acceptors and thus were named Pinus taeda xylan 2-O-arabinosyltransferase 1 (PtX2AT1) and 3-O-arabinosyltransferase 1 (PtX3AT1), respectively. Two other pine GT61 members were found to be xylan 2-O-xylosyltransferases (PtXYXTs) adding 2-O-ß-Xyl side chains onto xylooligomer acceptors. Furthermore, sequential reactions with PtX3AT1 and the PtGUX1 xylan glucuronyltransferase demonstrated that PtX3AT1 could efficiently arabinosylate glucuronic acid (GlcA)-substituted xylooligomers and likewise, PtGUX1 was able to add GlcA side chains onto 3-O-Araf-substituted xylooligomers. Phylogenetic analysis revealed that PtX2AT1, PtX3AT1 and PtXYXTs resided in three gymnosperm-specific GT61 clades that are separated from the grass-expanded GT61 clade harboring xylan 3-O-arabinosyltransferases and 2-O-xylosyltransferases, suggesting that they might have been recruited independently for xylan substitutions in gymnosperms. Together, our findings have established several pine GT61 members as xylan 2-O- and 3-O-arabinosyltransferases and 2-O-xylosyltransferases and they indicate that pine xylan might also be substituted with 2-O-Araf and 2-O-Xyl side chains.

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.

Evolution of complex genome architecture in gymnosperms.

Gymnosperms represent an ancient lineage that diverged from early spermatophytes during the Devonian. The long fossil records and low diversity in living species prove their complex evolutionary history, which included ancient radiations and massive extinctions. Due to their ultra-large genome size, the whole-genome assembly of gymnosperms has only generated in the past 10 years and is now being further expanded into more taxonomic representations. Here, we provide an overview of the publicly available gymnosperm genome resources and discuss their assembly quality and recent findings in large genome architectures. In particular, we describe the genomic features most related to changes affecting the whole genome. We also highlight new realizations relative to repetitive sequence dynamics, paleopolyploidy, and long introns. Based on the results of relevant genomic studies of gymnosperms, we suggest additional efforts should be made toward exploring the genomes of medium-sized (5-15 gigabases) species. Lastly, more comparative analyses among high-quality assemblies are needed to understand the genomic shifts and the early species diversification of seed plants.

Nutrient allocation strategies of four conifers from semiarid to extremely arid environments.

Although the contents of limiting elements in plants, such as nitrogen (N) and phosphorus (P), have been widely studied from subtropical to humid-temperate zones, the strategies used by coniferous species to allocation N and P in arid and semiarid forests remain unclear. In this study, samples of 545 leaves, 194 twigs, and 78 fine roots were collected from four coniferous species (Pinus tabuliformis, Picea wilsonii, Juniperus przewalskii, and Picea crassifolia) of three genera (Pinus, Picea, and Juniperus) in the northeastern Tibetan Plateau, and the contents of C, N, and P were analyzed. Two key parameters, namely the allometric exponent and coefficient of variation, were calculated to illustrate the relative investment of plants to N and P uptake and plasticity (variation of N and P), respectively. The contents of N and P and the N:P ratios were the highest in leaves, but their plasticity was the lowest. This confirmed the hypothesis that the leaves of coniferous species have a high content of limiting nutrients and homeostasis. At the regional level, the allometric exponent of N and P in leaves was 0.68, 0.74 in twigs, and 0.78 in fine roots, which is consistent with the results on a global scale. Thus, this invariant allometric relationship suggests the existence of an important mechanism that constrains the allocation of plant nutrients across broad environmental gradients. However, the allocation strategies for N and P shifted with the species, climate, and soil nutrients. Namely: their preferred nutrient uptake was P when the trees had a better nutritional status (semiarid environments, mean annual precipitations (MAP) > 300 mm), but the investment of N was strengthened when the habitat conditions become more severe (extremely arid environments, MAP <100 mm). Thus, our results can provide a novel perspective to understand the strategies of plant nutrient uptake in arid and semiarid forests.