Integrated small RNA, transcriptome and physiological approaches provide insight into Taxodium hybrid 'Zhongshanshan' roots in acclimation to prolonged flooding.

Although Taxodium hybrid 'Zhongshanshan' 406 (Taxodium mucronatum Tenore × Taxodium distichum; Taxodium 406) is an extremely flooding-tolerant woody plant, the physiological and molecular mechanisms underlying acclimation of its roots to long-term flooding remain largely unknown. Thus, we exposed saplings of Taxodium 406 to either non-flooding (control) or flooding for 2 months. Flooding resulted in reduced root biomass, which is in line with lower concentrations of citrate, α-ketoglutaric acid, fumaric acid, malic acid and adenosine triphosphate (ATP) in Taxodium 406 roots. Flooding led to elevated activities of pyruvate decarboxylase, alcohol dehydrogenase and lactate dehydrogenase, which is consistent with higher lactate concentration in the roots of Taxodium 406. Flooding brought about stimulated activities of superoxide dismutase and catalase and elevated reduced glutathione (GSH) concentration and GSH/oxidized glutathione, which is in agreement with reduced concentrations of O2- and H2O2 in Taxodium 406 roots. The levels of starch, soluble protein, indole-3-acetic acid, gibberellin A4 and jasmonate were decreased, whereas the concentrations of glucose, total non-structural carbohydrates, most amino acids and 1-aminocyclopropane-1-carboxylate (ACC) were improved in the roots of flooding-treated Taxodium 406. Underlying these changes in growth and physiological characteristics, 12,420 mRNAs and 42 miRNAs were significantly differentially expressed, and 886 miRNA-mRNA pairs were identified in the roots of flooding-exposed Taxodium 406. For instance, 1-aminocyclopropane-1-carboxylate synthase 8 (ACS8) was a target of Th-miR162-3p and 1-aminocyclopropane-1-carboxylate oxidase 4 (ACO4) was a target of Th-miR166i, and the downregulation of Th-miR162-3p and Th-miR166i results in the upregulation of ACS8 and ACO4, probably bringing about higher ACC content in flooding-treated roots. Overall, these results indicate that differentially expressed mRNA and miRNAs are involved in regulating tricarboxylic acid cycle, ATP production, fermentation, and metabolism of carbohydrates, amino acids and phytohormones, as well as reactive oxygen species detoxification of Taxodium 406 roots. These processes play pivotal roles in acclimation to flooding stress. These results will improve our understanding of the molecular and physiological bases underlying woody plant flooding acclimation and provide valuable insights into breeding-flooding tolerant trees.

Molecular Mechanism of Different Rooting Capacity between Two Clones of Taxodium hybrid 'Zhongshanshan'.

The conifer Taxodium hybrid 'Zhongshanshan' (T. hybrid 'Zhongshanshan') is characterized by rapid growth, strong stress resistance, and high ornamental value and has significant potential for use in afforestation, landscaping, and wood production. The main method of propagating T. hybrid 'Zhongshanshan' is tender branch cutting, but the cutting rooting abilities of different T. hybrid 'Zhongshanshan' clones differ significantly. To explore the causes of rooting ability differences at a molecular level, we analyzed the transcriptome data of cutting base and root tissues of T. hybrid 'Zhongshanshan 149' with a rooting rate of less than 5% and T. hybrid 'Zhongshanshan 118' with rooting rate greater than 60%, at the developmental time points in this study. The results indicated that differentially expressed genes between the two clones were mainly associated with copper ion binding, peroxidase, and oxidoreductase activity, response to oxidative stress, phenylpropanoid and flavonoid biosynthesis, and plant hormone signal transduction, among others. The expression pattern of ThAP2 was different throughout the development of the adventitive roots of the two clone cuttings. Therefore, this gene was selected for further study. It was shown that ThAP2 was a nuclear-localized transcription factor and demonstrated a positive feedback effect on rooting in transgenic Nicotiana benthamiana cuttings. Thus, the results of this study explain the molecular mechanism of cutting rooting and provide candidate gene resources for developing genetic breeding strategies for optimizing superior clones of T. hybrid 'Zhongshanshan'.

Identification and Functional Analysis of ThADH1 and ThADH4 Genes Involved in Tolerance to Waterlogging Stress in Taxodium hybrid 'Zhongshanshan 406'.

The Taxodium hybrid 'Zhongshanshan 406' (T. hybrid 'Zhongshanshan 406') [Taxodium mucronatum Tenore × Taxodium distichum (L.). Rich] has an outstanding advantage in flooding tolerance and thus has been widely used in wetland afforestation in China. Alcohol dehydrogenase genes (ADHs) played key roles in ethanol metabolism to maintain energy supply for plants in low-oxygen conditions. Two ADH genes were isolated and characterized-ThADH1 and ThADH4 (GenBank ID: AWL83216 and AWL83217-basing on the transcriptome data of T. hybrid 'Zhongshanshan 406' grown under waterlogging stress. Then the functions of these two genes were investigated through transient expression and overexpression. The results showed that the ThADH1 and ThADH4 proteins both fall under ADH III subfamily. ThADH1 was localized in the cytoplasm and nucleus, whereas ThADH4 was only localized in the cytoplasm. The expression of the two genes was stimulated by waterlogging and the expression level in roots was significantly higher than those in stems and leaves. The respective overexpression of ThADH1 and ThADH4 in Populus caused the opposite phenotype, while waterlogging tolerance of the two transgenic Populus significantly improved. Collectively, these results indicated that genes ThADH1 and ThADH4 were involved in the tolerance and adaptation to anaerobic conditions in T. hybrid 'Zhongshanshan 406'.

Comparative chloroplast genomics of the genus Taxodium.

BACKGROUND: Chloroplast (cp) genome information would facilitate the development and utilization of Taxodium resources. However, cp genome characteristics of Taxodium were poorly understood. RESULTS: We determined the complete cp genome sequences of T. distichum, T. mucronatum, and T. ascendens. The cp genomes are 131,947 bp to 132,613 bp in length, encode 120 genes with the same order, and lack typical inverted repeat (IR) regions. The longest small IR, a 282 bp trnQ-containing IR, were involved in the formation of isomers. Comparative analysis of the 3 cp genomes showed that 91.57% of the indels resulted in the periodic variation of tandem repeat (TR) motifs and 72.46% single nucleotide polymorphisms (SNPs) located closely to TRs, suggesting a relationship between TRs and mutational dynamics. Eleven hypervariable regions were identified as candidates for DNA barcode development. Hypothetical cp open reading frame 1(Ycf1) was the only one gene that has an indel in coding DNA sequence, and the indel is composed of a long TR. When extended to cupressophytes, ycf1 genes have undergone a universal insertion of TRs accompanied by extreme length expansion. Meanwhile, ycf1 also located in rearrangement endpoints of cupressophyte cp genomes. All these characteristics highlight the important role of repeats in the evolution of cp genomes. CONCLUSIONS: This study added new evidence for the role of repeats in the dynamics mechanism of cp genome mutation and rearrangement. Moreover, the information of TRs and hypervariable regions would provide reliable molecular resources for future research focusing on the infrageneric taxa identification, phylogenetic resolution, population structure and biodiversity for the genus Taxodium and Cupressophytes.

An Integrated Transcriptome and Proteome Analysis Reveals Putative Regulators of Adventitious Root Formation in Taxodium 'Zhongshanshan'.

Adventitious root (AR) formation from cuttings is the primary manner for the commercial vegetative propagation of trees. Cuttings is also the main method for the vegetative reproduction of Taxodium 'Zhongshanshan', while knowledge of the molecular mechanisms regulating the processes is limited. Here, we used mRNA sequencing and an isobaric tag for relative and absolute quantitation-based quantitative proteomic (iTRAQ) analysis to measure changes in gene and protein expression levels during AR formation in Taxodium 'Zhongshanshan'. Three comparison groups were established to represent the three developmental stages in the AR formation process. At the transcript level, 4743 genes showed an expression difference in the comparison groups as detected by RNA sequencing. At the protein level, 4005 proteins differed in their relative abundance levels, as indicated by the quantitative proteomic analysis. A comparison of the transcriptome and proteome data revealed regulatory aspects of metabolism during AR formation and development. In summary, hormonal signal transduction is different at different developmental stages during AR formation. Other factors related to carbohydrate and energy metabolism and protein degradation and some transcription factor activity levels, were also correlated with AR formation. Studying the identified genes and proteins will provide further insights into the molecular mechanisms controlling AR formation.

High-density genetic map construction and quantitative trait loci identification for growth traits in (Taxodium distichum var. distichum × T. mucronatum) × T. mucronatum.

BACKGROUND: 'Zhongshanshan' is the general designation for the superior interspecific hybrid clones of Taxodium species, which is widely grown for economic and ecological purposes in southern China. Growth is the priority objective in 'Zhongshanshan' tree improvement. A high-density linkage map is vital to efficiently identify key quantitative trait loci (QTLs) that affect growth. RESULTS: In total, 403.16 Gb of data, containing 2016,336 paired-end reads, was obtained after preprocessing. The average sequencing depth was 28.49 in T. distichum var. distichum, 25.18 in T. mucronatum, and 11.12 in each progeny. In total, 524,662 high-quality SLAFs were detected, of which 249,619 were polymorphic, and 6166 of the polymorphic markers met the requirements for use in constructing a genetic map. The final map harbored 6156 SLAF markers on 11 linkage groups, and was 1137.86 cM in length, with an average distance of 0.18 cM between adjacent markers. Separate QTL analyses of traits in different years by CIM detected 7 QTLs. While combining multiple-year data, 13 QTLs were detected by ICIM. 5 QTLs were repeatedly detected by the two methods, and among them, 3 significant QTLs (q6-2, q4-2 and q2-1) were detected in at least two traits. Bioinformatic analysis discoveried a gene annotated as a leucine-rich repeat receptor-like kinase gene within q4-2. CONCLUSIONS: This map is the most saturated one constructed in a Taxodiaceae species to date, and would provide useful information for future comparative mapping, genome assembly, and marker-assisted selection.

Inferences of population structure and demographic history for Taxodium distichum, a coniferous tree in North America, based on amplicon sequencing analysis.

PREMISE OF THE STUDY: Studies of natural genetic variation can elucidate the genetic basis of phenotypic variation and the past population structure of species. Our study species, Taxodium distichum, is a unique conifer that inhabits the flood plains and swamps of North America. Morphological and ecological differences in two varieties, T. distichum var. distichum (bald cypress) and T. distichum var. imbricarium (pond cypress), are well known, but little is known about the level of genetic differentiation between the varieties and the demographic history of local populations. METHODS: We analyzed nucleotide polymorphisms at 47 nuclear loci from 96 individuals collected from the Mississippi River Alluvial Valley (MRAV), and Gulf Coastal populations in Texas, Louisiana, and Florida using high-throughput DNA sequencing. Standard population genetic statistics were calculated, and demographic parameters were estimated using a composite-likelihood approach. KEY RESULTS: Taxodium distichum in North America can be divided into at least three genetic groups, bald cypress in the MRAV and Texas, bald cypress in Florida, and pond cypress in Florida. The levels of genetic differentiation among the groups were low but significant. Several loci showed the signatures of positive selection, which might be responsible for local adaptation or varietal differentiation. CONCLUSIONS: Bald cypress was genetically differentiated into two geographical groups, and the boundary was located between the MRAV and Florida. This differentiation could be explained by population expansion from east to west. Despite the overlap of the two varieties' ranges, they were genetically differentiated in Florida. The estimated demographic parameters suggested that pond cypress split from bald cypress during the late Miocene.

Transcriptome analysis of the Taxodium 'Zhongshanshan 405' roots in response to salinity stress.

Taxodium 'Zhongshanshan' is an interspecies hybrid of Taxodium distichum and Taxodium mucronatum, and has been widely planted in southeastern China. It has great ecological and economic potential. However, the scant genomic resources in genus Taxodium have greatly hindered further exploration of its underlying salinity-tolerance mechanism. To understand the genetic basis of its salt tolerance, high-throughput sequencing of mRNA (RNA-Seq) was used to analyze transcriptome changes of 'Zhongshanshan 405' clone roots treated with NaCl stress. After de novo assembly, 70,312 unigenes were achieved, and 41,059 of them were annotated. 9038 differentially expressed genes (DEGs) were identified among the treatments, and 7959 DEGs were found between salt-stressed roots and control, with 489 up-regulated and 570 down-regulated shared by all of the treatments. Genes related to transport, signal transductions as well as undescribed transcripts were among those DEGs in response to salt stress. Gene ontology classification analysis revealed that salt stress-related categories including 'oxidoreductase activity', 'metal ion binding', and 'membrane' were highly enriched among these DEGs. Moreover, the gene expression pattern of 12 unigenes revealed by quantitative real-time polymerase chain reaction (qRT-PCR) confirmed the RNA-Seq data. Our study not only provided the large-scale assessment of transcriptome resources of Taxodium but also guidelines for probing the molecular mechanism underlying 'Zhongshanshan' salt tolerance.

[Response in cuttings of Taxodium hybrid 'Zhongshanshan' and their parents to drought and re-hydration.] / ä¸­å±±æ‰åŠå ¶çˆ¶æ¯æœ¬å¹¼è‹—å¯¹å¹²æ—±èƒè¿«å’Œå¤æ°´çš„å“åº”.

In this study, three different strains of Taxodium hybrid 'Zhongshanshan' varieties [T. hybrid 302 (T. distichum♀×T. mucronatum♂), T. hybrid 407 (T. mucronatum♀×T. distichum♂), T. hybrid 118 (T. hybrid 302 ♀×T. macronatum ♂)] and their parents, T. distichum and T. mucronatum, were applied to investigate the response of photosynthetic characteristics, antioxidant enzyme systems and morphological characteristics to drought stress and recovery. The results indicated that as drought days were prolonged, all plants' net photosynthetic rate (Pn) decreased, while proline accumulated. Meanwhile, the antioxidases functioned to eliminate malonaldehyde toxicity. On the 8th day, the decrease of Pn of T. distichum was the biggest, T. hybrid 118 plants showed the highest water use efficiency and the smallest MDA content, while T. macronatum plants increased the activity of superoxide dismutase and content of proline. After rewatering for 2 days, all these parameters showed signs of recovery, and the T. hybrid 118 plants showed the fastest recovery rate since their Pn and proline content had recovered for 74.4% and 60.2%, respectively. Then after recovered for 9 days, all tested parameters had almost restored to equivalent levels of CK plants. The total biomass of T. hybrid 118 plants was not affected, while the ratio of root to shoot was significantly (P<0.05) increased. The drought-resistance capacity ranged as T. macronatum>T. hybrid 118>T. hybrid 407>T. hybrid 302>T. distichum plants. In conclusion, the backcross generation T. hybrid 118 plants largely inherited the drought resistance of T. mucronatum, and the result would be instrumental in breeding and popularization of drought-resistant hybrid varieties.

Evolutionary rate variation in two conifer species, Taxodium distichum (L.) Rich. var. distichum (baldcypress) and Cryptomeria japonica (Thunb. ex L.f.) D. Don (Sugi, Japanese cedar).

With the advance of sequencing technologies, large-scale data of expressed sequence tags and full-length cDNA sequences have been reported for several conifer species. Comparative analyses of evolutionary rates among diverse taxa provide insights into taxon-specific molecular evolutionary features and into the origin of variation in evolutionary rates within genomes and between species. Here, we estimated evolutionary rates in two conifer species, Taxodium distichum and Cryptomeria japonica, to illuminate the molecular evolutionary features of these species, using hundreds of genes and employing Chamaecyparis obtusa as an outgroup. Our results show that the mutation rates based on synonymous substitution rates (dS) of T. distichum and C. japonica are approximately 0.67 × 10(-9) and 0.59 × 10(-9)/site/year, respectively, which are 15-25 times lower than those of annual angiosperms. We found a significant positive correlation between dS and GC3. This implies that a local mutation bias, such as context dependency of the mutation bias, exists within the genomes of T. distichum and C. japonica, and/or that selection acts on synonymous sites in these species. In addition, the means of the ratios of synonymous to nonsynonymous substitution rate in the two species are almost the same, suggesting that the average intensity of functional constraint is constant between the lineages. Finally, we tested the possibility of positive selection based on the site model, and detected one candidate gene for positive selection.

Characterization of the genome of bald cypress.

BACKGROUND: Bald cypress (Taxodium distichum var. distichum) is a coniferous tree of tremendous ecological and economic importance. It is a member of the family Cupressaceae which also includes cypresses, redwoods, sequoias, thujas, and junipers. While the bald cypress genome is more than three times the size of the human genome, its 1C DNA content is amongst the smallest of any conifer. To learn more about the genome of bald cypress and gain insight into the evolution of Cupressaceae genomes, we performed a Cot analysis and used Cot filtration to study Taxodium DNA. Additionally, we constructed a 6.7 genome-equivalent BAC library that we screened with known Taxodium genes and select repeats. RESULTS: The bald cypress genome is composed of 90% repetitive DNA with most sequences being found in low to mid copy numbers. The most abundant repeats are found in fewer than 25,000 copies per genome. Approximately 7.4% of the genome is single/low-copy DNA (i.e., sequences found in 1 to 5 copies). Sequencing of highly repetitive Cot clones indicates that most Taxodium repeats are highly diverged from previously characterized plant repeat sequences. The bald cypress BAC library consists of 606,336 clones (average insert size of 113 kb) and collectively provides 6.7-fold genome equivalent coverage of the bald cypress genome. Macroarray screening with known genes produced, on average, about 1.5 positive clones per probe per genome-equivalent. Library screening with Cot-1 DNA revealed that approximately 83% of BAC clones contain repetitive sequences iterated 103 to 104 times per genome. CONCLUSIONS: The BAC library for bald cypress is the first to be generated for a conifer species outside of the family Pinaceae. The Taxodium BAC library was shown to be useful in gene isolation and genome characterization and should be an important tool in gymnosperm comparative genomics, physical mapping, genome sequencing, and gene/polymorphism discovery. The single/low-copy (SL) component of bald cypress is 4.6 times the size of the Arabidopsis genome. As suggested for other gymnosperms, the large amount of SL DNA in Taxodium is likely the result of divergence among ancient repeat copies and gene/pseudogene duplication.

Contrasting patterns of DNA variation in natural populations of two related conifers, Cryptomeria japonica and Taxodium distichum (Cupressaceae sensu lato).

We investigated DNA variation within and between two closely related conifers, Cryptomeria japonica and Taxodium distichum, at nuclear loci encoding ferredoxin, glutamyl-tRNA reductase, lycopene beta cyclase, and phosphoribosylanthranilate transferase. Average nucleotide diversity at silent sites was estimated to be 0.0035 (SE 0.0012) in C. japonica and 0.0058 (SE 0.0006) in T. distichum. One population in C. japonica was differentiated from the others but generally there was not much differentiation among populations or varieties within the two species. However, the two species seemed to differ in frequency spectra of DNA polymorphisms. Excesses of intermediate-frequency variants were found in C. japonica, whereas excesses of both rare and high-frequency variants were found in T. distichum, which suggested different histories of population structures in the two species. Deviations from the standard neutral expectations in DNA polymorphisms were found by applications of neutrality tests. The results show that actions of selection to respective loci seem to differ between the two species, indicating differences of interaction among evolutionary factors.

DNA variation in a conifer, Cryptomeria japonica (Cupressaceae sensu lato).

We investigated the nucleotide variation of a conifer, Cryptomeria japonica, and the divergence between this species and its closest relative, Taxodium distichum, at seven nuclear loci (Acl5, Chi1, Ferr, GapC, HemA, Lcyb, and Pat). Samples of C. japonica were collected from three areas, Kantou-Toukai, Hokuriku, and Iwate. No apparent geographic differentiation was found among these samples. However, the frequency spectrum of the nucleotide polymorphism revealed excesses of intermediate-frequency variants, which suggests that the population was not panmictic and a constant size in the past. The average nucleotide diversity, pi, for silent sites was 0.00383. However, values of pi for silent sites vary among loci. Comparisons of polymorphism to divergence among loci (the HKA test) showed that the polymorphism at the Acl5 locus was significantly lower. We also observed a nearly significant excess of replacement polymorphisms at the Lcyb locus. These results suggested possibilities of natural selection acting at some of the loci. Intragenic recombination was detected only once at the Chi1 locus and was not detected at the other loci. The low level of population recombination rate, 4Nr, seemed to be due to both low level of recombination, r, and small population size, N.