The Metasequoia genome and evolutionary relationships among redwoods.

Redwood trees (Sequoioideae), including Metasequoia glyptostroboides (dawn redwood), Sequoiadendron giganteum (giant sequoia), and Sequoia sempervirens (coast redwood), are threatened and widely recognized iconic tree species. Genomic resources for redwood trees could provide clues to their evolutionary relationships. Here, we report the 8-Gb reference genome of M. glyptostroboides and a comparative analysis with two related species. More than 62% of the M. glyptostroboides genome is composed of repetitive sequences. Clade-specific bursts of long terminal repeat retrotransposons may have contributed to genomic differentiation in the three species. The chromosomal synteny between M. glyptostroboides and S. giganteum is extremely high, whereas there has been significant chromosome reorganization in S. sempervirens. Phylogenetic analysis of marker genes indicates that S. sempervirens is an autopolyploid, and more than 48% of the gene trees are incongruent with the species tree. Results of multiple analyses suggest that incomplete lineage sorting (ILS) rather than hybridization explains the inconsistent phylogeny, indicating that genetic variation among redwoods may be due to random retention of polymorphisms in ancestral populations. Functional analysis of ortholog groups indicates that gene families of ion channels, tannin biosynthesis enzymes, and transcription factors for meristem maintenance have expanded in S. giganteum and S. sempervirens, which is consistent with their extreme height. As a wetland-tolerant species, M. glyptostroboides shows a transcriptional response to flooding stress that is conserved with that of analyzed angiosperm species. Our study offers insights into redwood evolution and adaptation and provides genomic resources to aid in their conservation and management.

Assembled and annotated 26.5 Gbp coast redwood genome: a resource for estimating evolutionary adaptive potential and investigating hexaploid origin.

Sequencing, assembly, and annotation of the 26.5 Gbp hexaploid genome of coast redwood (Sequoia sempervirens) was completed leading toward discovery of genes related to climate adaptation and investigation of the origin of the hexaploid genome. Deep-coverage short-read Illumina sequencing data from haploid tissue from a single seed were combined with long-read Oxford Nanopore Technologies sequencing data from diploid needle tissue to create an initial assembly, which was then scaffolded using proximity ligation data to produce a highly contiguous final assembly, SESE 2.1, with a scaffold N50 size of 44.9 Mbp. The assembly included several scaffolds that span entire chromosome arms, confirmed by the presence of telomere and centromere sequences on the ends of the scaffolds. The structural annotation produced 118,906 genes with 113 containing introns that exceed 500 Kbp in length and one reaching 2 Mb. Nearly 19 Gbp of the genome represented repetitive content with the vast majority characterized as long terminal repeats, with a 2.9:1 ratio of Copia to Gypsy elements that may aid in gene expression control. Comparison of coast redwood to other conifers revealed species-specific expansions for a plethora of abiotic and biotic stress response genes, including those involved in fungal disease resistance, detoxification, and physical injury/structural remodeling and others supporting flavonoid biosynthesis. Analysis of multiple genes that exist in triplicate in coast redwood but only once in its diploid relative, giant sequoia, supports a previous hypothesis that the hexaploidy is the result of autopolyploidy rather than any hybridizations with separate but closely related conifer species.

Genome-wide association identifies candidate genes for drought tolerance in coast redwood and giant sequoia.

Drought is a major limitation for survival and growth in plants. With more frequent and severe drought episodes occurring due to climate change, it is imperative to understand the genomic and physiological basis of drought tolerance to be able to predict how species will respond in the future. In this study, univariate and multitrait multivariate genome-wide association study methods were used to identify candidate genes in two iconic and ecosystem-dominating species of the western USA, coast redwood and giant sequoia, using 10 drought-related physiological and anatomical traits and genome-wide sequence-capture single nucleotide polymorphisms. Population-level phenotypic variation was found in carbon isotope discrimination, osmotic pressure at full turgor, xylem hydraulic diameter, and total area of transporting fibers in both species. Our study identified new 78 new marker × trait associations in coast redwood and six in giant sequoia, with genes involved in a range of metabolic, stress, and signaling pathways, among other functions. This study contributes to a better understanding of the genomic basis of drought tolerance in long-generation conifers and helps guide current and future conservation efforts in the species.

Selective Sweeps and Polygenic Adaptation Drive Local Adaptation along Moisture and Temperature Gradients in Natural Populations of Coast Redwood and Giant Sequoia.

Dissecting the genomic basis of local adaptation is a major goal in evolutionary biology and conservation science. Rapid changes in the climate pose significant challenges to the survival of natural populations, and the genomic basis of long-generation plant species is still poorly understood. Here, we investigated genome-wide climate adaptation in giant sequoia and coast redwood, two iconic and ecologically important tree species. We used a combination of univariate and multivariate genotype-environment association methods and a selective sweep analysis using non-overlapping sliding windows. We identified genomic regions of potential adaptive importance, showing strong associations to moisture variables and mean annual temperature. Our results found a complex architecture of climate adaptation in the species, with genomic regions showing signatures of selective sweeps, polygenic adaptation, or a combination of both, suggesting recent or ongoing climate adaptation along moisture and temperature gradients in giant sequoia and coast redwood. The results of this study provide a first step toward identifying genomic regions of adaptive significance in the species and will provide information to guide management and conservation strategies that seek to maximize adaptive potential in the face of climate change.

Genetic structure of coast redwood (Sequoia sempervirens [D. Don] Endl.) populations in and outside of the natural distribution range based on nuclear and chloroplast microsatellite markers.

Coast redwood (Sequoia sempervirens) naturally growing in southern Oregon and northern California is one of the few conifer tree species that are polyploid. Despite its unique ecological and economic importance, its population genetic structure is still insufficiently studied. To obtain additional data on its population genetic structure we genotyped 317 samples collected from populations in California (data set C) and 144 trees growing in a provenance trial in France (data set F) using 12 nuclear (five random nuclear genomic nSSRs and seven expressed sequence tag EST-SSRs) and six chloroplast (cpSSRs) microsatellite or simple sequence repeat (SSR) markers, respectively. These data sets were also used as reference to infer the origin of 147 coast redwood trees growing in Germany (data set G). Coast redwood was introduced to Europe, including Germany as an ornamental species, decades ago. Due to its fast growth and high timber quality, it could be considered as a potential commercial timber species, especially in perspective to climate warming that makes more regions in Germany suitable for its growing. The well performing trees in colder Germany could be potential frost resistant genotypes, but their genetic properties and origin are mostly unknown. Within the natural range in southern Oregon and northern California, only two relatively weak clusters were identified, one northern and one southern, separated by the San Francisco Bay. High genetic diversity, but low differentiation was found based on the 12 nuclear SSR markers for all three data sets F, C and G. We found that investigated 147 German trees represented only 37 different genotypes. They showed genetic diversity at the level less than diversity observed within the natural range in the northern or southern cluster, but more similar to the diversity observed in the southern cluster. It was difficult to assign German trees to the original single native populations using the six cpSSR markers, but rather to either the northern or southern cluster. The high number of haplotypes found in the data sets based on six cpSSR markers and low genetic differentiation based on 12 nuclear SSRs found in this study helps us study and better understand population genetic structure of this complex polyploid tree and supports the selection of potential genotypes for German forestry.

De novo transcriptome assembly of cold stressed clones of the hexaploid Sequoia sempervirens (D. Don) Endl.

Coast redwood is a very important endemic conifer timber species in Southern Oregon and Northern California in the USA. Due to its good wood properties and fast growth rate it can be considered as a prospective timber species also in other countries with similar or changing toward similar climatic conditions due to global climate warming, such as Germany. In general, it is frost sensitive and suffers from freezing temperatures. To study genetic mechanisms of frost resistance in this species and to select the most frost tolerant trees we tested 17 clones in climate control chamber experiments and generated two de novo assemblies of the coast redwood transcriptome from a pooled RNA sample using Trinity and CLC Genomic Workbench software, respectively. The hexaploid nature of the coast redwood genome makes it very challenging to successfully assemble and annotate the coast redwood transcriptome. The de novo transcriptome assembly generated by Trinity and CLC considering only reads with a minimum length of 180 bp and contigs no less than 200 bp long resulted in 634,772 and 788,464 unigenes (unique contigs), respectively.

Genotype-by-environment interaction in coast redwood outside natural distribution - search for environmental cues.

BACKGROUND: Effective matching of genotypes and environments is required for the species to reach optimal productivity and act effectively for carbon sequestration. A common garden experiment across five different environments was undertaken to assess genotype x environment interaction (GxE) of coast redwood in order to understand the performance of genotypes across environments. RESULTS: The quantitative genetic analysis discovered no GxE between investigated environments for diameter at breast height (DBH). However, no genetic component was detected at one environment possibly due to stressful conditions. The implementation of universal response function allowed for the identification of important environmental factors affecting species productivity. Additionally, this approach enabled us to predict the performance of species across the New Zealand environmental conditions. CONCLUSIONS: In combination with quantitative genetic analysis which identified genetically superior material, the URF model can directly identify the optimal geographical regions to maximize productivity. However, the finding of ideally uncorrelated climatic variables for species with narrow ecological amplitude is rather challenging, which complicates construction of informative URF model. This, along with a small number of tested environments, tended to overfit a prediction model which resulted in extreme predictions in untested environments.

Whole genome duplication in coast redwood (Sequoia sempervirens) and its implications for explaining the rarity of polyploidy in conifers.

Polyploidy is common and an important evolutionary factor in most land plant lineages, but it is rare in gymnosperms. Coast redwood (Sequoia sempervirens) is one of just two polyploid conifer species and the only hexaploid. Evidence from fossil guard cell size suggests that polyploidy in Sequoia dates to the Eocene. Numerous hypotheses about the mechanism of polyploidy and parental genome donors have been proposed, based primarily on morphological and cytological data, but it remains unclear how Sequoia became polyploid and why this lineage overcame an apparent gymnosperm barrier to whole-genome duplication (WGD). We sequenced transcriptomes and used phylogenetic inference, Bayesian concordance analysis and paralog age distributions to resolve relationships among gene copies in hexaploid coast redwood and close relatives. Our data show that hexaploidy in coast redwood is best explained by autopolyploidy or, if there was allopolyploidy, it happened within the Californian redwood clade. We found that duplicate genes have more similar sequences than expected, given the age of the inferred polyploidization. Conflict between molecular and fossil estimates of WGD can be explained if diploidization occurred very slowly following polyploidization. We extrapolate from this to suggest that the rarity of polyploidy in gymnosperms may be due to slow diploidization in this clade.

Contrasting drought-response strategies in California redwoods.

We compared the physiology and growth of seedlings originating from different Sequoia sempervirens (D. Don.) Endl. (coast redwood) and Sequoiadendron giganteum (Lindl.) Buchh. (giant sequoia) populations subjected to progressive drought followed by a recovery period in a controlled greenhouse experiment. Our objective was to examine how multiple plant traits interact to influence the response of seedlings of each species and seed population to a single drought and recovery cycle. We measured soil and plant water status, leaf gas exchange, stem embolism and growth of control (well-watered) and drought-stressed (water withheld) seedlings from each population at the beginning, middle and end of a 6-week drought period and again 2 weeks after re-watering. The drought had a significant effect on many aspects of seedling performance, but water-stressed seedlings regained most physiological functioning by the end of the recovery period. Sequoiadendron seedlings exhibited a greater degree of isohydry (water status regulation), lower levels of stem embolism, higher biomass allocation to roots and lower sensitivity of growth to drought compared with Sequoia. Only minor intra-specific differences were observed among populations. Our results show that seedlings of the two redwood species exhibit contrasting drought-response strategies that align with the environmental conditions these trees experience in their native habitats, and demonstrate trade-offs and coordination among traits affecting plant water use, carbon gain and growth under drought.

Subpicosecond surface dynamics in genomic DNA from in vitro-grown plant species: a SERS assessment.

In this work the surface-enhanced Raman total half band widths of seven genomic DNAs from leaves of chrysanthemum (Dendranthema grandiflora Ramat.), common sundew (Drosera rotundifolia L.), edelweiss (Leontopodium alpinum Cass), Epilobium hirsutum L., Hypericum richeri ssp. transsilvanicum (Celak) Ciocârlan, rose (Rosa x hybrida L.) and redwood (Sequoia sempervirens D. Don. Endl.) have been measured. We have shown that surface-enhanced Raman spectroscopy (SERS) can be used to study the fast subpicosecond dynamics of DNA in the proximity of a metallic surface. The dependencies of the total half band widths and the global relaxation times, on the DNA molecular subgroup structure and on the type of genomic DNA, are reported. In our study, the full widths at half-maximum (FWHMs) for the SERS bands of genomic DNAs from different leaf tissues are typically in the wavenumber range from 15 to 55 cm(-1). Besides, it can be observed that molecular relaxation processes studied in this work have a global relaxation time smaller than 0.71 ps and larger than 0.19 ps. A comparison between different ranges of FT-Raman and SERS band parameters, respectively, corresponding to DNA extracted from leaf tissues is given. It is shown that the interaction between DNA and a metallic surface has the potential to lead to a shortening of the global relaxation times, as compared with molecular dynamics in solution. We have found that the surface dynamics of molecular subgroups in plant DNA is, in some cases, about two times faster than the solution dynamics of nucleic acids. This can be rationalized in a qualitative manner by invoking the complex landscape of the interaction energy between the molecule and the silver surface.

An EST dataset for Metasequoia glyptostroboides buds: the first EST resource for molecular genomics studies in Metasequoia.

The "living fossil" Metasequoia glyptostroboides Hu et Cheng, commonly known as dawn redwood or Chinese redwood, is the only living species in the genus and is valued for its essential oil and crude extracts that have great potential for anti-fungal activity. Despite its paleontological significance and economical value as a rare relict species, genomic resources of Metasequoia are very limited. In order to gain insight into the molecular mechanisms behind the formation of reproductive buds and the transition from vegetative phase to reproductive phase in Metasequoia, we performed sequencing of expressed sequence tags from Metasequoia vegetative buds and female buds. By using the 454 pyrosequencing technology, a total of 1,571,764 high-quality reads were generated, among which 733,128 were from vegetative buds and 775,636 were from female buds. These EST reads were clustered and assembled into 114,124 putative unique transcripts (PUTs) with an average length of 536 bp. The 97,565 PUTs that were at least 100 bp in length were functionally annotated by a similarity search against public databases and assigned with Gene Ontology (GO) terms. A total of 59 known floral gene families and 190 isotigs involved in hormone regulation were captured in the dataset. Furthermore, a set of PUTs differentially expressed in vegetative and reproductive buds, as well as SSR motifs and high confidence SNPs, were identified. This is the first large-scale expressed sequence tags ever generated in Metasequoia and the first evidence for floral genes in this critically endangered deciduous conifer species.

Small RNAs of Sequoia sempervirens during rejuvenation and phase change.

In this work, the population of small RNAs (sRNAs) was studied in the gymnosperm Sequoia sempervirens during phase changes, specifically in the juvenile, adult and rejuvenated plants obtained in vitro. The potential target genes of Sequoia sRNAs were predicted through bioinformatics. Rejuvenation is a pivotal process in woody plants that enables them to regain their growth potential, which results in the recovery of physiologic and molecular characteristics that were lost when the juveniles mature into adult plants. The results from the five repeated graftings of juvenile, adult and rejuvenated plants in vitro showed that sRNAs could be classified into structural RNAs (Group I), small interfering RNAs (Group II), annotated microRNAs (Group III, and unannotated sRNAs (Group IV). The results indicate that only 573 among 15,485,415 sRNAs (Groups III and IV) had significantly different expression patterns associated with rejuvenation and phase change. A total of 215 sRNAs exhibited up-regulated expression patterns in adult shoots, and 358 sRNAs were down-regulated. Expression profiling and prediction of possible target genes of these unique small RNAs indicate possible functions in the control of photosynthetic efficiency and rooting competence abundance during plant rejuvenation. Moreover, the increase in SsmiR156 and decrease in SsmiR172 during plant rejuvenation suggested that these two microRNAs extensively affect phase transition.

DNA methylation and genome rearrangement characteristics of phase change in cultured shoots of Sequoia sempervirens.

Epigenetic machinery regulates the expression of individual genes and plays a crucial role in globally shaping and maintaining developmental patterning. We studied the extent of DNA methylation in the nucleus, mitochondrion and chloroplast in cultured Sequoia sempervirens (coast redwood) adult, juvenile and rejuvenated shoots by measuring the ratio of methylcytosine to total cytosine using high-performance liquid chromatography (HPLC). We also analyzed nuclear DNA (nuDNA) polymorphisms of different shoot types by methylation-sensitive amplified fragment length polymorphism (MSAP) and Southern blot analysis. The extent of nuDNA methylation was greater in the adult vegetative than juvenile and rejuvenated shoots (8% vs 6.5-7.5%). In contrast, the proportion of methylcytosine was higher in mitochondrial DNA (mDNA) of juvenile and rejuvenated shoots than adult shoots (6.6% vs 7.8-8.2%). MSAP and Southern blot analyses identified three MSAP fragments which could be applied as phase-specific molecular markers. We also found nuclear genome and mtDNA rearrangement may be as important as DNA methylation status during the phase change. Our findings strongly suggest that DNA methylation and genome rearrangement may affect the dynamic tissue- and cell type-specific changes that determine the developmental phase of S. sempervirens shoots.

Clones identification of Sequoia sempervirens (D. Don) Endl. in Chile by using PCR-RAPDs technique.

A protocol of polymerase chain reaction-random amplified polymorphic DNAs (PCR-RAPDs) was established to analyse the gene diversity and genotype identification for clones of Sequoia sempervirens (D. Don) Endl. in Chile. Ten (out of 34) clones from introduction trial located in Voipir-Villarrica, Chile, were studied. The PCR-RAPDs technique and a modified hexadecyltrimethylammonium bromide (CTAB) protocol were used for genomic DNA extraction. The PCR tests were carried out employing 10-mer random primers. The amplification products were detected by electrophoresis in agarose gels. Forty nine polymorphic bands were obtained with the selected primers (BG04, BF07, BF12, BF13, and BF14) and were ordered according to their molecular size. The genetic similarity between samples was calculated by the Jaccard index and a dendrogram was constructed using a cluster analysis of unweighted pair group method using arithmetic averages (UPGMA). Of the primers tested, 5 (out of 60) RAPD primers were selected for their reproducibility and high polymorphism. A total of 49 polymorphic RAPD bands were detected out of 252 bands. The genetic similarity analysis demonstrates an extensive genetic variability between the tested clones and the dendrogram depicts the genetic relationships among the clones, suggesting a geographic relationship. The results indicate that the RAPD markers permitted the identification of the assayed clones, although they are derived from the same geographic origin.