Identification and characterization of the WOX Gene Family revealed two WUS Clade Members associated with embryo development in Cunninghamia lanceolata.

The WUSCHEL-related homeobox (WOX) gene family is vital for plant development and stress response. In this study, we conducted a comprehensive analysis of WOX genes in Cunninghamia lanceolata (C. lanceolata) and subsequently explored the potential roles of two ClWOX genes within the WUS clade. In total, six ClWOX genes were identified through a full-length transcriptome analysis. These genes, exhibiting conserved structural and functional motifs, were assigned to the ancient clade and Modern/WUS clade, respectively, through a phylogenetic analysis. Our expression analysis indicated that these ClWOX genes were highly expressed in the middle and late developmental stages of zygotic embryos in C. lanceolata. Moreover, only ClWOX5 and ClWOX6 within the Modern/WUS clade exhibited transcriptional activity, and their expressions were also induced in response to auxin and wounding. Overexpression of ClWOX5 and ClWOX6 in Arabidopsis caused a partially sterile phenotype, resulting in a very low seed setting rate. Transcriptomic analysis revealed that expressions of many embryo-defective (EMB) genes, phytohormone-related genes, and transcription factors (TFs) were dramatically altered in ClWOX5 and ClWOX6 transgenic plants, which suggested that ClWOX5 and ClWOX6 may play specific important roles in embryo development via complex gene networks. In addition, overexpression of ClWOX5 and ClWOX6 in leaf segments promoted shoot regeneration in tobacco, indicating that ClWOX5 and ClWOX6 can promote plant regeneration and could be used to improve genetic transformation. In conclusion, these results help to elucidate the function of the WOX gene and provide a valuable basis for future studies of the developmental regulation and applications of WOX genes in C. lanceolata.

De novo assembly of Iron-Heart Cunninghamia lanceolata transcriptome and EST-SSR marker development for genetic diversity analysis.

Iron-Heart Cunninghamia lanceolata, a wild relative of Chinese fir with valuable genetic and breeding traits, has been limited in genetic studies due to a lack of genomic resources and markers. In this study, we conducted transcriptome sequencing of Iron-Heart C. lanceolata leaves using Illumina NovaSeq 6000 and performed assembly and analysis. We obtained 45,326,576 clean reads and 115,501 unigenes. Comparative analysis in five functional databases resulted in successful annotation of 26,278 unigenes, with 6,693 unigenes annotated in all databases (5.79% of the total). UniProt and Pfam databases provided annotations for 22,673 and 18,315 unigenes, respectively. Gene Ontology analysis categorized 23,962 unigenes into three categories. KEGG database alignment annotated 10,195 unigenes, classifying them into five categories: metabolism, genetic information, biological systems, cellular processes, and environmental information processing. From the unigenes, we identified 5,645 SSRs, with dinucleotides repeats being the most common (41.47%). We observed variations in repeat numbers and base compositions, with the majority of markers ranging from 12 to 29 bp in length. We randomly selected 200 primer pairs and successfully amplified 15 pairs of polymorphic SSR primers, which effectively distinguished Chinese fir plants of different origins. This study provides insights into the genetic characteristics of Iron-Heart C. lanceolata and offers a foundation for future molecular marker development, breeding programs, genetic diversity analysis, and conservation strategies.

ClNAC100 Is a NAC Transcription Factor of Chinese Fir in Response to Phosphate Starvation.

Phosphate (Pi) deficiency is one of the most limiting factors for Chinese fir growth and production. Moreover, continuous cultivation of Chinese fir for multiple generations led to the reduction of soil nutrients, which hindered the yield of Chinese fir in southern China. Although NAC (NAM, ATAF, and CUC) transcription factors (TFs) play critical roles in plant development and abiotic stress resistance, it is still unclear how they regulate the response of Chinese fir to phosphate (Pi) starvation. Based on Pi-deficient transcriptome data of Chinses fir root, we identified a NAC transcription factor with increased expression under Pi deficiency, which was obtained by PCR and named ClNAC100. RT-qPCR confirmed that the expression of ClNAC100 in the root of Chinese fir was induced by phosphate deficiency and showed a dynamic change with time. It was positively regulated by ABA and negatively regulated by JA, and ClNAC100 was highly expressed in the roots and leaves of Chinese fir. Transcriptional activation assay confirmed that ClNAC100 was a transcriptional activator. The promoter of ClNAC100 was obtained by genome walking, which was predicted to contain a large number of stress, hormone, and growth-related cis-elements. Tobacco infection was used to verify the activity of the promoter, and the core promoter was located between -1519 bp and -589 bp. We identified 18 proteins bound to the ClNAC100 promoter and 5 ClNAC100 interacting proteins by yeast one-hybrid and yeast two-hybrid, respectively. We speculated that AHL and TIFY family transcription factors, calmodulin, and E3 ubiquitin ligase in these proteins might be important phosphorus-related proteins. These results provide a basis for the further study of the regulatory mechanism and pathways of ClNAC100 under Pi starvation.

Genome-Wide Characterization and Gene Expression Analyses of Malate Dehydrogenase (MDH) Genes in Low-Phosphorus Stress Tolerance of Chinese Fir (Cunninghamia lanceolata).

Malate dehydrogenase (MDH) genes play vital roles in developmental control and environmental stress tolerance in sessile plants by modulating the organic acid-malic acid level. However, MDH genes have not yet been characterized in gymnosperm, and their roles in nutrient deficiency are largely unexplored. In this study, 12 MDH genes were identified in Chinese fir (Cunninghamia lanceolata), namely, ClMDH-1, -2, -3, …, and -12. Chinese fir is one of the most abundant commercial timber trees in China, and low phosphorus has limited its growth and production due to the acidic soil of southern China. According to the phylogenetic analysis, MDH genes were classified into five groups, and Group 2 genes (ClMDH-7, -8, -9, and 10) were only found to be present in Chinese fir but not in Arabidopsis thaliana and Populus trichocarpa. In particular, the Group 2 MDHs also had specific functional domains-Ldh_1_N (malidase NAD-binding functional domain) and Ldh_1_C (malate enzyme C-terminal functional domain)-indicating a specific function of ClMDHs in the accumulation of malate. All ClMDH genes contained the conserved MDH gene characteristic functional domains Ldh_1_N and Ldh_1_C, and all ClMDH proteins exhibited similar structures. Twelve ClMDH genes were identified from eight chromosomes, involving fifteen ClMDH homologous gene pairs, each with a Ka/Ks ratio of <1. The analysis of cis-elements, protein interactions, and transcription factor interactions of MDHs showed that the ClMDH gene might play a role in plant growth and development, and in response to stress mechanisms. The results of transcriptome data and qRT-PCR validation based on low-phosphorus stress showed that ClMDH1, ClMDH6, ClMDH7, ClMDH2, ClMDH4, ClMDH5, ClMDH10 and ClMDH11 were upregulated under low-phosphorus stress and played a role in the response of fir to low-phosphorus stress. In conclusion, these findings lay a foundation for further improving the genetic mechanism of the ClMDH gene family in response to low-phosphorus stress, exploring the potential function of this gene, promoting the improvement of fir genetics and breeding, and improving production efficiency.

Soil nitrification and denitrification in Cunninghamia lanceolata plantations with different stand ages.

The variations in soil nitrification and denitrification processes, together with the abundances of functional microbes were investigated in Cunninghamia lanceolata plantations with different stand ages of 5, 8, 21, 27, and 40 years old. The results showed that the net nitrification rate fluctuated with increasing forest ages, with that of 8-year- and 27-year-old C. lanceolata plantations being significantly lower than other stand ages. The abundance of ammonia-oxidizing archaea (AOA) amoA in the 27-year-old plantation was significantly lower than that of the 40-year-old plantation, while there was no significant difference among the other stand ages. There was no significant difference in the abundance of AOB amoA gene, denitrifying functional genes or soil denitrification potential among different stand ages. The results of stepwise regression analysis showed that the abundance of AOA amoA gene was not significantly affected by soil physical and chemical properties. In addition, the abundance of AOB was positively associated with soil total carbon content and soil pH. The abundance of denitrifying functional genes including narG, nirK and nosZ increased with increasing soil pH. The abundance of nirK and nirS was influenced by soil total carbon. Stand age influenced soil net nitrification rate through the AOA amoA abundance. Moreover, soil denitrification potential was directly affected by stand age, or indirectly affected by stand age through soil microbial biomass carbon, soil pH and denitrifying gene abundance of narG and nirK. Compared with the denitrification process, soil nitrification and associated AOA amoA gene abundance were more sensitive to the development of C. lanceolata plantations. The rotation period sould be appropriately extended to reduce the risk of nitrogen losses resulting from soil nitrification.

Intraspecific variations in leaf functional traits of Cunninghamia lanceolata provenances.

BACKGROUND: Studies on intra-specific variability in leaf functional traits is important to evaluate adaptation of the species to predicted climate change, and to develop long-term conservation strategy. The main objectives were to investigate the relationship between the functional traits leaves and C, N, P stoichiometry of Chinese fir from different geographical provenances and their relationship with the main environmental factors of provenance. RESULTS: In this study, we measured 12 leaf functional traits on 36-year-old Cunninghamia lanceolata trees from 13 provenances. Analysis of variance (ANOVA) was performed to examine the variability. Redundancy analysis (RA) was computed to examine the relationship between geo-climatic factors of provenance origin and leaf functional traits while Pearson's correlation coefficient was computed to assess inter-trait correlations. The results showed statistically significant differences (P < 0.01) in intraspecific leaf traits among provenances, except leaf P content. The relationships among leaf traits are consistent with the general trend observed in the leaf economic spectrum. Mean annual temperature appeared to be a key factor that influences intraspecific leaf traits variability compared to mean annual precipitation. CONCLUSION: These results provide useful insights about adaptation of leaf trait of Chinese fir in a changing climatic condition. Thus, our findings shed light on the importance of interspecific trait variability in Chinese fir and the potential effect of climate change.

The plant peptide hormone phytosulfokine promotes somatic embryogenesis by maintaining redox homeostasis in Cunninghamia lanceolata.

Somatic embryogenesis (SE) is widely used for studying the mechanisms of embryo development. However, little is known about the underlying mechanisms, especially in woody plants. Previous studies have established an SE system for Chinese fir (Cunninghamia lanceolata), but this system is genotype-dependent, which limits its application in practice. Here, we found that phytosulfokine (PSK), a plant peptide hormone, can not only increase SE efficiency, but also establish SE in recalcitrant genotypes of C. lanceolata. Proembryogenic mass (PEM) browning and determination of hydrogen peroxide (H2 O2 ) content by 2',7'-dichlorofluorescein staining indicated that a reactive oxygen species (ROS) burst occurred rapidly after PEMs were transferred to SE induction medium. Transcriptome analysis and quantitative reverse transcriptase-PCR validation showed that PSK treatment helped to maintain ROS homeostasis by decreasing the activity of peroxidases in early SE induction. This PSK-regulated redox microenvironment might be helpful to induce expression of SE-related genes like WOX2 in early SE induction. Further analyses suggested that PSK promotes SE induction in C. lanceolata partially through decreasing H2 O2 levels, which is necessary but not sufficient for SE induction in recalcitrant genotypes of C. lanceolata. Furthermore, heterologous overexpression of ClPSK in Arabidopsis led to enhanced SE induction and resistance to H2 O2 stress. Taken together, our study reveals a biological function for the plant peptide hormone PSK, extends our knowledge about SE in woody trees, and provides a valuable tool for establishing an efficient and genotype-independent SE system in C. lanceolata and other coniferous trees.

A real-time PCR for detection of pathogens of anthracnose on Chinese fir using TaqMan probe targeting ApMat gene.

BACKGROUND: Anthracnose is one of the most widespread and destructive diseases on Chinese fir. Colletotrichum cangyuanense, Colletotrichum fructicola, Colletotrichum gloeosporioides, and Colletotrichum siamense are the causal agents of anthracnose on Chinese fir. A rapid and accurate diagnosis of different pathogens is critical for the disease management. RESULTS: Phylogenetic tree and sequence similarity analysis showed that the single-locus ApMat provides superior phylogenetic information and is an appropriate target to distinguish C. cangyuanense, C. fructicola, C. gloeosporioides, and C. siamense. The real-time PCR assays with the primer sets of MQ-F/R, 1#C-F/R, YK-F/R, and WZ-F/R, and corresponding TaqMan probes of MQ-P, 1#C-P, YK-P, and WZ-P were specific for C. cangyuanense, C. fructicola, C. gloeosporioides, and C. siamense, respectively. The sensitivity tests showed that the lowest amount of gDNA that the PCRs can detect was 1 ng of genomic DNA. The validity of the assays was confirmed by directly detecting the pathogens from both the fungal culture and infected Chinese fir. CONCLUSION: These results demonstrated the potential of the TaqMan real-time PCR targeting the ApMat gene to provide rapid, specific, and reliable molecular detection of C. fructicola, C. gloeosporioides, C. siamense, and C. cangyuanense, respectively. The data also provided a reference solution for the identification of species within Colletotrichum gloeosporioides species complex (CGSC), which share similar morphological characteristics. © 2022 Society of Chemical Industry.

Transcriptome Profiling of Stem-Differentiating Xylem in Response to Abiotic Stresses Based on Hybrid Sequencing in Cunninghamia lanceolata.

Cunninghamia lanceolata (C. lanceolata) belongs to Gymnospermae, which are fast-growing and have desirable wood properties. However, C. lanceolata's stress resistance is little understood. To unravel the physiological and molecular regulation mechanisms under environmental stresses in the typical gymnosperm species of C. lanceolata, three-year-old plants were exposed to simulated drought stress (polyethylene glycol 8000), salicylic acid, and cold treatment at 4 °C for 8 h, 32 h, and 56 h, respectively. Regarding the physiological traits, we observed a decreased protein content and increased peroxidase upon salicylic acid and polyethylene glycol treatment. Superoxide dismutase activity either decreased or increased at first and then returned to normal under the stresses. Regarding the molecular regulation, we used both nanopore direct RNA sequencing and short-read sequencing to reveal a total of 5646 differentially expressed genes in response to different stresses, of which most had functions in lignin catabolism, pectin catabolism, and xylan metabolism, indicating that the development of stem-differentiating xylem was affected upon stress treatment. Finally, we identified a total of 51 AP2/ERF, 29 NAC, and 37 WRKY transcript factors in C. lanceolata. The expression of most of the NAC TFs increased under cold stress, and the expression of most of the WRKY TFs increased under cold and SA stress. These results revealed the transcriptomics responses in C. lanceolata to short-term stresses under this study's experimental conditions and provide preliminary clues about stem-differentiating xylem changes associated with different stresses.

Inbreeding in Chinese Fir: Insight into the Rare Self-Fertilizing Event from a Genetic View.

Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.) is a fast-growing conifer with great forestation value and prefers outcrossing with high inbreeding depression effect. Previously, we captured a special Chinese fir parent clone named as 'cx569' that lacks early inbreeding depression. In view of the fact that very little has been published about the rare self-fertilizing event in Chinese fir from a genetic view, herein, we conduct an SSR-based study on the variation of open- and self-pollinated offspring of this parent to gain a view of the rare self-fertilizing event. The results indicated that genetic diversity of self-pollinated offspring was significantly reduced by half (Ho: 0.302, vs. 0.595, p = 0.001; He: 0.274 vs. 0.512, p = 0.002) when compared to an open-pollinated set. Self-pollinated offspring also had significantly positive FIS values (FIS = 0.057, p = 0.034) with a much higher proportion of common allele (20.59% vs. 0), reflecting their heterozygote deficiency. Clustering analysis further indicated a separation of the self- and opened- pollinated groups, implying a natural preference of outcrossing for cx569. However, the cx569 still had 6% acceptance for selfing. When accepted 100% for its own pollen, the cx569 led to a genetically unique selfing group. Additionally, this selfing group seemed to be consistently homozygous at seven particular loci. These findings gave us more genetic clues to gain insight into the rare self-fertilizing event in conifer (Chinese fir).

Further Mining and Characterization of miRNA Resource in Chinese Fir (Cunninghamia lanceolata).

In this study, we aimed to expand the current miRNA data bank of Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.) regarding its potential value for further genetic and genomic use in this species. High-throughput small RNA sequencing successfully captured 140 miRNAs from a Chinese fir selfing family harboring vigor and depressed progeny. Strikingly, 75.7% (n = 106) of these miRNAs have not been documented previously, and most (n = 105) of them belong to the novel set with 6858 putative target genes. The new datasets were then integrated with the previous information to gain insight into miRNA genetic architecture in Chinese fir. Collectively, a relatively high proportion (62%, n = 110) of novel miRNAs were found. Furthermore, we identified one MIR536 family that has not been previously documented in this species and four overlapped miRNA families (MIR159, MIR164, MIR171_1, and MIR396) from new datasets. Regarding the stability, we calculated the secondary structure free energy and found a relatively low R2 value (R2 < 0.22) between low minimal folding free energy (MFE) of pre-miRNAs and MFE of its corresponding mature miRNAs in most datasets. When in view of the conservation aspect, the phylogenetic trees showed that MIR536 and MIR159 sequences were highly conserved in gymnosperms.

Evaluation of reference genes and characterization of the MYBs in xylem radial change of Chinese fir stem.

The radial change (RC) of tree stem is the process of heartwood formation involved in complex molecular mechanism. Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.), an evergreen species, is an important fast-growing timber tree in southern China. In this study, the top four stable genes (IDH, UBC2, RCA and H2B) were selected in RC tissues of 15 years old Chinese fir stem (RC15) and the genes (H2B, 18S, TIP41 and GAPDH) were selected in RC tissues of 30 years old Chinese fir stem (RC30). The stability of the reference genes is higher in RC30 than in RC15. Sixty-one MYB transcripts were obtained on the PacBio Sequel platform from woody tissues of one 30 years old Chinese fir stem. Based on the number of MYB DNA-binding domain and phylogenetic relationships, the ClMYB transcripts contained 21 transcripts of MYB-related proteins (1R-MYB), 39 transcripts of R2R3-MYB proteins (2R-MYB), one transcript of R1R2R3-MYB protein (3R-MYB) belonged to 18 function-annotated clades and two function-unknown clades. In RC woody tissues of 30 years old Chinese fir stem, ClMYB22 was the transcript with the greatest fold change detected by both RNA-seq and qRT-PCR. Reference genes selected in this study will be helpful for further verification of transcript abundance patterns during the heartwood formation of Chinese fir.

Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir).

BACKGROUND: R2R3-MYB is a class of transcription factor crucial in regulating secondary cell wall development during wood formation. The regulation of wood formation in gymnosperm has been understudied due to its large genome size. Using Single-Molecule Real-Time sequencing, we obtained full-length transcriptomic libraries from the developmental stem of Cunninghamia lanceolata, a perennial conifer known as Chinese fir. The R2R3-MYB of C. lanceolata (hereafter named as ClMYB) associated with secondary wall development were identified based on phylogenetic analysis, expression studies and functional study on transgenic line. RESULTS: The evolutionary relationship of 52 ClMYBs with those from Arabidopsis thaliana, Eucalyptus grandis, Populus trichocarpa, Oryza sativa, two gymnosperm species, Pinus taeda, and Picea glauca were established by neighbour-joining phylogenetic analysis. A large number of ClMYBs resided in the woody-expanded subgroups that predominated with the members from woody dicots. In contrast, the woody-preferential subgroup strictly carrying the members of woody dicots contained only one candidate. The results suggest that the woody-expanded subgroup emerges before the gymnosperm/angiosperm split, while most of the woody-preferential subgroups are likely lineage-specific to woody dicots. Nine candidates shared the same subgroups with the A. thaliana orthologs, with known function in regulating secondary wall development. Gene expression analysis inferred that ClMYB1/2/3/4/5/26/27/49/51 might participate in secondary wall development, among which ClMYB1/2/5/26/27/49 were significantly upregulated in the highly lignified compression wood region, reinforcing their regulatory role associated with secondary wall development. ClMYB1 was experimentally proven a transcriptional activator that localised in the nucleus. The overexpression of ClMYB1 in Nicotiana benthamiana resulted in an increased lignin deposition in the stems. The members of subgroup S4, ClMYB3/4/5 shared the ERF-associated amphiphilic repression motif with AtMYB4, which is known to repress the metabolism of phenylpropanoid derived compounds. They also carried a core motif specific to gymnosperm lineage, suggesting divergence of the regulatory process compared to the angiosperms. CONCLUSIONS: This work will enrich the collection of full-length gymnosperm-specific R2R3-MYBs related to stem development and contribute to understanding their evolutionary relationship with angiosperm species.

Transcriptome analysis provides insights into the root response of Chinese fir to phosphorus deficiency.

BACKGROUND: Phosphorus is one of the essential elements for plant growth and development, but available phosphorus (Pi) content in many soil types is low. As a fast-growing tree species for timber production, Chinese fir is in great demand of Pi, and the lack of Pi in soil restricts the increase of productivity of Chinese fir plantation. Root morphology and the synthesis and secretion of organic acids play an important role in the uptake of phosphorus, but the molecular mechanisms of Chinese fir root responses to Pi deficiency are largely unexplored. In this study, seedlings of Yang 061 clone were grown under three Pi supply levels (0, 5 and 10 mg·L-1 P) and morphological attributes, organic acid content and enzyme activity were measured. The transcriptome data of Chinese fir root system were obtained and the expression levels of phosphorus responsive genes and organic acid synthesis related genes on citric acid and glyoxylate cycle pathway were determined. RESULTS: We annotated 50,808 Unigenes from the transcriptome of Chinese fir roots. Among differentially expressed genes, seven genes of phosphate transporter family and 17 genes of purple acid phosphatase family were up-regulated by Pi deficiency, two proteins of SPX domain were up-regulated and one was down-regulated. The metabolic pathways of the citric acid and glyoxylate cycle pathway were mapped, and the expression characteristics of the related Unigenes under different phosphorus treatments were analyzed. The genes involved in malic acid and citric acid synthesis were up-regulated, and the activities of the related enzymes were significantly enhanced under long-term Pi stress. The contents of citric acid and malic acid in the roots of Chinese fir increased after 30 days of Pi deficiency. CONCLUSION: Chinese fir roots showed increased expression of genes related with phosphorus starvation, citrate and malate synthesis genes, increased content of organic acids, and enhanced activities of related enzymes under Pi deficiency. The results provide a new insight for revealing the molecular mechanism of adaption to Pi deficiency and the pathway of organic acid synthesis in Chinese fir roots.

Drought stress introduces growth, physiological traits and ecological stoichiometry changes in two contrasting Cunninghamia lanceolata cultivars planted in continuous-plantation soils.

BACKGROUND: The decrease in Cunninghamia lanceolata (Lamb.) production on continuously planted soil is an essential problem. In this study, two-year-old seedlings of two cultivars (a normal cultivar, NC, and a super cultivar, SC) were grown in two types of soil (not planted (NP) soil; continuously planted (CP) soil) with three watering regimes, and the interactive effects on plant growth and physiological traits were investigated in a greenhouse experiment. The water contents of the soil in the control (CK) (normal water content), medium water content (MWC) and low water content (LWC) treatments reached 75-80 %, 45-50 % and 20-25 % of the field water capacity, respectively. RESULTS: The results indicated that the CP soil had a negative effect on growth and physiological traits and that the LWC treatment caused even more severe and comprehensive negative effects. In both cultivars, the CP soil significantly decreased the height increment (HI), basal diameter increment (DI), dry matter accumulation (DMA), net photosynthetic rate (Pn), total chlorophyll content (TChl), carotenoid content (Caro) and photosynthetic nitrogen use efficiency (PNUE). Compared to the NP soil, the CP soil also decreased the proline and soluble protein contents, nitrogen use efficiency (NUE) and phosphorus use efficiency (PUE) and increased the nitrogen:phosphorus ratio in roots, stems and leaves. The LWC treatment decreased growth and photosynthesis, changed ecological stoichiometry, induced oxidative stress, promoted water use efficiency and damaged chloroplast ultrastructure. Significant increases in ascorbate peroxidase (APX), peroxidase (POD), soluble protein and proline contents were found in the LWC treatment. Compared with the NC, the SC was more tolerant to the CP soil and water stress, as indicated by the higher levels of DMA, Pn, and WUE. After exposure to the CP soil and watering regimes, the decreases in biomass accumulation and gas exchange were more pronounced. CONCLUSIONS: The combination of drought and CP soil may have detrimental effects on C. lanceolata growth, and low water content enhances the impacts of CP soil stress on C. lanceolata seedlings. The superiority of the SC over the NC is significant in Chinese fir plantation soil. Therefore, continuously planted soil can be utilized to cultivate improved varieties of C. lanceolata and maintain water capacity. This can improve their growth and physiological performance to a certain extent.

Half-leaf width symmetric distribution reveals buffering strategy of Cunninghamia lanceolata.

BACKGROUND: Leaf length and width could be a functioning relationship naturally as plant designs. Single-vein leaves have the simplest symmetrical distribution and structural design, which means that fast-growing single-vein species could interpret the scheme more efficiently. The distribution of leaf length and width can be modulated for better adaptation, providing an informative perspective on the various operational strategies in an emergency, while this mechanism is less clear. Here we selected six age groups of Cunninghamia lanceolata pure forests, including saplings, juveniles, mature, and old-growth trees. We pioneered a tapering model to describe half-leaf symmetric distribution with mathematical approximation based on every measured leaf along developmental sequence, and evaluated the ratio of leaf basal part length to total length (called tipping leaf length ratio). RESULTS: The tipping leaf length ratio varied among different tree ages. That means the changes of tipping leaf length ratio and leaf shape are a significant but less-noticed reflection of trees tradeoff strategies at different growth stages. For instance, there exhibited relatively low ratio during sapling and juvenile, then increased with increasing age, showing the highest value in their maturity, and finally decreased on mature to old-growth transition. The tipping leaf length ratio serves as a cost-benefit ratio, thus the subtle changes in the leaf symmetrical distribution within individuals reveal buffering strategy, indicating the selection for efficient design of growth and hydraulic in their developmental sequences. CONCLUSIONS: Our model provides a physical explanation of varied signatures for tree operations in hydraulic buffering through growth stages, and the buffering strategy revealed from leaf distribution morphologically provides evidence on the regulation mechanism of leaf biomechanics, hydraulics and physiologies. Our insight contributes greatly to plant trait modeling, policy and management, and will be of interest to some scientists and policy makers who are involved in climate change, ecology and environment protection, as well as forest ecology and management.

The Transcriptome of Cunninghamia lanceolata male/female cone reveal the association between MIKC MADS-box genes and reproductive organs development.

BACKGROUND: Cunninghamia lanceolata (Chinese fir), a member of the conifer family Cupressaceae, is one of the most popular cultivated trees for wood production in China. Continuous research is being performed to improve C. lanceolata breeding values. Given the high rate of seed abortion (one of the reasons being the failure of ovule and pollen development) in C. lanceolata, the proper formation of female/male cones could theoretically increase the number of offspring in future generations. MIKC MADS-box genes are well-known for their roles in the flower/cone development and comprise the typical/atypical floral development model for both angiosperms and gymnosperms. RESULTS: We performed a transcriptomic analysis to find genes differentially expressed between female and male cones at a single, carefully determined developmental stage, focusing on the MIKC MADS-box genes. We finally obtained 47 unique MIKC MADS-box genes from C. lanceolata and divided these genes into separate branches. 27 out of the 47 MIKC MADS-box genes showed differential expression between female and male cones, and most of them were not expressed in leaves. Out of these 27 genes, most B-class genes (AP3/PI) were up-regulated in the male cone, while TM8 genes were up-regulated in the female cone. Then, with no obvious overall preference for AG (class C + D) genes in female/male cones, it seems likely that these genes are involved in the development of both cones. Finally, a small number of genes such as GGM7, SVP, AGL15, that were specifically expressed in female/male cones, making them candidate genes for sex-specific cone development. CONCLUSIONS: Our study identified a number of MIKC MADS-box genes showing differential expression between female and male cones in C. lanceolata, illustrating a potential link of these genes with C. lanceolata cone development. On the basis of this, we postulated a possible cone development model for C. lanceolata. The gene expression library showing differential expression between female and male cones shown here, can be used to discover unknown regulatory networks related to sex-specific cone development in the future.

Genome survey of Chinese fir (Cunninghamia lanceolata): Identification of genomic SSRs and demonstration of their utility in genetic diversity analysis.

Chinese fir (Cunninghamia lanceolata) is an important coniferous species that accounts for 20-30% of the total commercial timber production in China. Though traditional breeding of Chinese fir has achieved remarkable success, molecular-assisted breeding has made little progress due to limited availability of genomic information. In this study, a survey of Chinese fir genome was performed using the Illumina HiSeq Xten sequencing platform. K-mer analysis indicated that Chinese fir has a large genome of approximately 11.6 Gb with 74.89% repetitive elements and is highly heterozygous. Meanwhile, its genome size was estimated to be 13.2 Gb using flow cytometry. A total of 778.02 Gb clean reads were assembled into 10,982,272 scaffolds with an N50 of 1.57 kb. In total, 362,193 SSR loci were detected with a frequency of 13.18 kb. Dinucleotide repeats were the most abundant (up to 73.6% of the total SSRs), followed by trinucleotide and tetranucleotide repeats. Forty-six polymorphic pairs were developed, and 298 alleles were successfully amplified from 199 Chinese fir clones. The average PIC value was 0.53, indicating that the identified genomic SSR (gSSR) markers have a high degree of polymorphism. In addition, these breeding resources were divided into three groups, and a limited gene flow existed among these inferred groups.

Transcriptomic Revelation of Phenolic Compounds Involved in Aluminum Toxicity Responses in Roots of Cunninghamia lanceolata (Lamb.) Hook.

: Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.) is one of the most important coniferous evergreen tree species in South China due to its desirable attributes of fast growth and production of strong and hardy wood. However, the yield of Chinese fir is often inhibited by aluminum (Al) toxicity in acidic soils of South China. Understanding the molecular mechanisms of Chinese fir root responses to Al toxicity might help to further increase its productivity. Here we used the Illumina Hiseq4000 platform to carry out transcriptome analysis of Chinese fir roots subjected to Al toxicity conditions. A total of 88.88 Gb of clean data was generated from 12 samples and assembled into 105,732 distinct unigenes. The average length and N50 length of these unigenes were 839 bp and 1411 bp, respectively. Among them, 58362 unigenes were annotated through searches of five public databases (Nr: NCBI non-redundant protein sequences, Swiss-Prot: A manually annotated and reviewed protein sequence database, GO: Gene Ontology, KOG/COG: Clusters of Orthologous Groups of proteins, and KEGG: the Kyoto Encyclopedia of Genes and Genomes database), which led to association of unigenes with 44 GO terms. Plus, 1615 transcription factors (TFs) were functionally classified. Then, differentially expressed genes (DEGs, |log2(fold change)| ≥ 1 and FDR ≤ 0.05) were identified in comparisons labelled TC1 (CK-72 h/CK-1 h) and TC2 (Al-72 h/Al-1 h). A large number of TC2 DEGs group were identified, with most being down-regulated under Al stress, while TC1 DEGs were primarily up-regulated. Combining GO, KEGG, and MapMan pathway analysis indicated that many DEGs are involved in primary metabolism, including cell wall metabolism and lipid metabolism, while other DEGs are associated with signaling pathways and secondary metabolism, including flavonoids and phenylpropanoids metabolism. Furthermore, TFs identified in TC1 and TC2 DEGs represented 21 and 40 transcription factor families, respectively. Among them, expression of bHLH, C2H2, ERF, bZIP, GRAS, and MYB TFs changed considerably under Al stress, which suggests that these TFs might play crucial roles in Chinese fir root responses to Al toxicity. These differentially expressed TFs might act in concert with flavonoid and phenylpropanoid pathway genes in fulfilling of key roles in Chinese fir roots responding to Al toxicity.

Genetic Characterization of Chinese fir from Six Provinces in Southern China and Construction of a Core Collection.

Large ex situ germplasm collections of plants generally contain significant diversity. A set of 700 well-conserved Chinese fir (Cunninghamia lanceolata (Lamb.) Hook) clones from six provinces in southern China in the ex situ gene bank of Longshan State Forest, was analyzed using 21 simple sequence repeat markers, with the aim of assessing the genetic diversity of these germplasm resources. Genetic analysis revealed extensive genetic variation among the accessions, with an average of 8.31 alleles per locus and a mean Shannon index of 1.331. Excluding loci with null alleles, we obtained a low level of genetic differentiation among provinces, consistent with the interpopulation genetic variation (1%). Three clusters were identified by STRUCTURE, which did not match the individuals' geographical provenances. Ten traits related to growth and wood properties were quantified in these individuals, and there was substantial variation in all traits across individuals, these provide a potential source of variation for genetic improvement of the Chinese fir. Screening large collections for multiple-trait selective breeding programs is laborious and expensive; a core collection of 300 accessions, representative of the germplasm, was established, based on genotypic and phenotypic data. The identified small, but diverse, collections will be useful for further genome-wide association studies.