Inbreeding in Chinese fir: Insights into the adaptive growth traits of selfed progeny from mRNA, miRNA, and copy number variation.

PREMISE: The impact of inbreeding on biological processes is well documented in individuals with severe inbreeding depression. However, the biological processes influencing the adaptive growth of normal selfed individuals are unknown. Here, we aimed to investigate how inbreeding affects gene expression for adaptive growth of normal selfed seedlings from a self-fertilizing parent in Chinese fir (Cunninghamia lanceolata). METHODS: Using RNA-seq data from above- and underground tissues of abnormal and normal selfed seedlings, we analyzed GO biological processes network. We also sequenced small RNAs in the aboveground tissues and measured the copy number variations (CNV) of the hub genes. RESULTS: Phenotypic fitness analysis revealed that the normal seedlings were better adapted than their abnormal counterparts. Upregulated differentially expressed genes (DEGs) were associated with development processes, and downregulated DEGs were mainly enriched in fundamental metabolism and stress response. Results of mRNA-miRNA parallel sequencing revealed that upregulated target genes were predominantly associated with development, highlighting their crucial role in phosphorylation in signal transduction networks. We also discovered a moderate correlation (0.1328 < R2 < 0.6257) between CNV and gene expression levels for three hub genes (TMKL1, GT2, and RHY1A). CONCLUSIONS: We uncovered the key biological processes underpinning the growth of normal selfed seedlings and established the relationship between CNV and the expression levels of hub genes in selfed seedlings. Understanding the candidate genes involved in the growth of selfed seedlings will help us comprehend the genetic mechanisms behind inbreeding depression in the evolutionary biology of plants.

Characterization of hemicellulose in Cunninghamia lanceolata stem during xylogenesis.

In this study, hemicellulose was isolated from the apical, middle and basal segments of C. lanceolata stem to investigate the dynamic change of its structure during xylogenesis. Results showed that the C. lanceolata hemicellulose is mainly consisted of O-acetylgalactoglucomannan (GGM) which backbone is alternately linked by ß-d-mannopyranosyl (Manp) and ß-d-glucopyranosyl (Glcp) via (1 â†’ 4)-glycosidic bond, while the side chains are α-d-galactopyranosyl (Galp) and acetyl. In addition, 4-O-methylglucuronoarabinoxylan (GAX) is another dominant structure of C. lanceolata hemicellulose which contains a linear backbone of (1 â†’ 4)-ß-d-xylopyranosyl (Xylp) and side chains of 4-O-Me-α-d-glucuronic acid (MeGlcpA) and α-L-arabinofuranose (Araf). The thickness of the cell wall, the ratio of GGM/GAX and the molecular weight of hemicellulose were increased as the extension of growth time. The degree of glycosyl substitutions of xylan and mannan was decreased from 10.34 % (apical) to 8.38 % (basal) and from 15.63 % (apical) to 10.49 % (basal), respectively. However, the total degree of acetylation was enhanced from 0.28 (apical) to 0.37 (basal). Transcriptome analysis showed that genes (CSLA9, IRX9H1, IRX10L, IRX15L, GMGT1, TBL19, TBL25, GUX2, GUX3, GXM1, F8H1 and F8H2) related to hemicellulose biosynthesis are mainly expressed in mature part. This study is of great significance for genetic breeding and high-value utilization of C. lanceolata.

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.

Genome-Wide Association Study With Growth-Related Traits and Secondary Metabolite Contents in Red- and White-Heart Chinese Fir.

Chinese fir [Cunninghamia lanceolata (Lamb.) Hook] is an important evergreen coniferous tree species that is widely distributed in many southern provinces of China and has important economic value. The Chinese fir accounts for 1/4 and 1/3 of the total artificial forest area and stock volume, respectively. Red-heart Chinese fir is popular in the market because of its high density and red heartwood. The long-growth cycle hindered the breeding process of Chinese fir, while molecular marker-assisted breeding could accelerate it. However, Chinese fir, a perennial conifer species, has a large genome, which has not yet been published. In this study, the growth-related traits and secondary metabolite contents of red- and white-heart Chinese fir were measured and found to be different between them. There are extremely significant differences among growth-related traits (p < 0.001), but secondary metabolite contents have different correlations due to differences in chemical structure. Moreover, genotype effect analysis of the substantially correlated single nucleotide polymorphisms (SNPs) revealed that most of the loci related to each growth-related traits were different from each other, indicating a type specificity of the genes regulated different growth-related traits. Furthermore, among the loci related to secondary metabolite contents, nine loci associated with multiple metabolite phenotypes such as Marker21022_4, Marker21022_172, Marker24559_31, Marker27425_37, Marker20748_85, Marker18841_115, Marker18841_198, Marker65846_146, and Marker21486_163, suggesting the presence of pleiotropic genes. This study identified the potential SNP markers associated with secondary metabolites in Chinese fir, thus setting the basis for molecular marker-assisted selection.

Metabolite Profiling and Transcriptome Analysis Unveil the Mechanisms of Red-Heart Chinese Fir [Cunninghamia lanceolata (Lamb.) Hook] Heartwood Coloration.

Red-heart Chinese fir (Cunninghamia lanceolata) has the advantages of high density and attractive color, making it popular in the market. To date, most studies about stems of woody plants have only been reported at the cytological level because of few living cells. In this study, the xylem was successfully partitioned into three effective sampling areas: sapwood, transition zone, and heartwood. Secondary metabolites, cell survival, and differentially expressed genes in the three sampling areas were, respectively, investigated. First, we identified the phenylpropanoid and flavonoid pathways closely related to color. Based on the chemical structure of secondary metabolites in pathways, two notable directions had been found. Luteolin's glycosylation products might be the key substances that regulated the color of heartwood in red-heart Chinese fir because of the 1,000-fold difference between red-heart and white-heart. We also found pinocembrin and pinobanksin in Chinese fir, which were rarely reported before. At the cytological level, we believed that the transition zone of red-heart Chinese fir was a critical region for color production because of the fewer living ray parenchyma cells. In addition, transcriptome and quantitative reverse transcription PCR (qRT-PCR) proved that genes regulating the entire phenylpropanoid pathway, upstream of the flavonoid pathway, and some glycosyltransferases were significantly upregulated in the transition zone of red-heart and then colored the heartwood by increasing metabolites. This is the first report on the color-related secondary metabolites regulated by differential genes in red-heart Chinese fir. This study will broaden our knowledge on the effects of metabolites on coloring woody plant xylems.

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.

Somatic Embryogenesis of Immature Cunninghamia lanceolata (Lamb.) Hook Zygotic Embryos.

Two efficient somatic embryogenesis systems were developed in Chinese fir, the most important conifer for industrial wood production in China. Three development stages (cleavage polyembryony, dominant embryo, and precotyledon) of immature embryos derived from 25 genotypes of open-pollinated mother trees were used as initial explants. Cleavage polyembryony-stage embryos with a 12.44% induction rate was the most embryogenic response stage. The highest frequency of embryogenic callus (13.86%) induction was obtained from DCR medium supplemented with 1.5 mg L-1 2,4-dichlorophenoxyacetic acid (2,4-D) and 0.3 mg L-1 kinetin (KN). An average of 53.33 early somatic embryos were produced from approximately 0.2 g (fresh weight) embryogenic callus after 2 weeks of incubation on medium supplemented with 50 µmol L-1 abscisic acid (ABA) and 100 g L-1 polyethylene glycol (PEG) 6000. About 53% dominant embryos have an embryogenic response after a 6-week cultivation on medium supplemented with 1.0-2.0 mg L-1 benzyladenine (BA), 0.2 mg L-1 naphthylacetic acid (NAA) or 2,4-D, and 0.004 mg L-1 thidiazuron (TDZ). After three successive transfer cultures on medium supplemented with 1.5 mg L-1 BA, 0.2 mg L-1 NAA, and 0.004 mg L-1 TDZ, 4.49-16.51% of the embryos developed into somatic embryos.

Global Reprogramming of Transcription in Chinese Fir (Cunninghamia lanceolata) during Progressive Drought Stress and after Rewatering.

Chinese fir (Cunninghamia lanceolata), an evergreen conifer, is the most commonly grown afforestation species in southeast China due to its rapid growth and good wood qualities. To gain a better understanding of the drought-signalling pathway and the molecular metabolic reactions involved in the drought response, we performed a genome-wide transcription analysis using RNA sequence data. In this study, Chinese fir plantlets were subjected to progressively prolonged drought stress, up to 15 d, followed by rewatering under controlled environmental conditions. Based on observed morphological changes, plantlets experienced mild, moderate, or severe water stress before rehydration. Transcriptome analysis of plantlets, representing control and mild, moderate, and severe drought-stress treatments, and the rewatered plantlets, identified several thousand genes whose expression was altered in response to drought stress. Many genes whose expression was tightly coupled to the levels of drought stress were identified, suggesting involvement in Chinese fir drought adaptation responses. These genes were associated with transcription factors, signal transport, stress kinases, phytohormone signalling, and defence/stress response. The present study provides the most comprehensive transcriptome resource and the first dynamic transcriptome profiles of Chinese fir under drought stress. The drought-responsive genes identified in this study could provide further information for understanding the mechanisms of drought tolerance in Chinese fir.

Ectopic expression of a poplar APETALA3-like gene in tobacco causes early flowering and fast growth.

A MADS-box gene, designated PtAP3, was isolated from a floral bud cDNA library derived from Populus tomentosa. Analysis by multiple alignments of both nucleotide and amino acid sequences, together with phylogenetic analysis, revealed that PtAP3 is an ortholog of Arabidopsis AP3. Analysis of RNA extracts from vegetative and reproductive tissues of P. tomentosa by RT-PCR indicated that PtAP3 is expressed in roots, stems, leaves and vegetative and floral buds. Notably, the expression of PtAP3 fluctuated during floral bud development between September and February with differences between male and female buds. In the former, a gradual down-regulation during this period, interrupted by a slight up-regulation in December, was followed by a sharper up-regulation on February. In developing female floral buds, expression was stable from September to November, sharply up-regulated in December, and then gradually down-regulated until February. The functional role of PtAP3 was investigated in transgenic tobacco plants. Of 25 transformants, nine displayed an earlier flowering phenotype compared with the wild type plants. Furthermore, transgenic tobacco had faster growth and more leaves than untransformed controls. The traits proved to be heritable between the T0 and T1 generations. Our results demonstrate a regulatory role of the PtAP3 gene during plant flowering and growth and suggest that the gene may be an interesting target for genetic modification to induce early flowering in plants.

Bidirectionalization of a methyl jasmonate-inducible plant promoter.

A unidirectional promoter can be transformed into a bidirectional module by artificial methods. Here we report the bidirectionalization of the methyl jasmonate (MeJA)-inducible PtDrl02 promoter derived from poplar [(Populus tomentosa × P. bolleana) × P. tomentosa] in planta. Construction of the bidirectional PtDrl02 promoter (designated as mPtDrl02) was rapidly achieved by introducing a minimal 35S promoter in the opposite orientation to the 5' end of PtDrl02. ß-Glucuronidase (GUS) and green fluorescent protein (GFP) reporter genes were also interchangeably linked to the 3'-and 5'-ends of mPtDrl02 to produce GFP/mPtDrl02/GUS and GUS/mPtDrl02/GFP vectors, respectively. Using the Agrobacterium-mediated transient expression approach, we demonstrated that the mPtDrl02 module was able to drive gene (GUS and GFP) expression in both orientations simultaneously. Furthermore, the cooperative and concurrent activity from both directions of the mPtDrl02 module was demonstrated following MeJA induction. To our knowledge, this is the first report of an artificial MeJA-responsive bidirectional promoter in perennial plants.

Functional identification and regulation of the PtDrl02 gene promoter from triploid white poplar.

The PtDrl02 gene belongs to the TIR-NBS gene family in triploid white poplar (Populus tomentosa x P. bolleana) x P. tomentosa. Its expression pattern displays tissue-specificity, and the transcript level can be induced by wounding, methyl jasmonate (MeJA), and salicylic acid (SA). To understand the regulatory mechanism controlling PtDrl02 gene expression, we functionally characterized the PtDrl02 promoter region. Using the beta-glucuronidase as a reporter, we found that the PtDrl02 promoter directed gene expression mainly in the aerial parts of the plants and was confined to the cortex tissues of leaf veins, petioles, stems, and stem piths, showing a typical tissue-specific expression pattern. Deletion analysis revealed two positive regulatory regions (-985 to -669 and -669 to -467) responsible for the basal activity of the PtDrl02 promoter. Impressively, the sequence from -669 to -467 was shown to contain cis-element (s) responding to wounding and MeJA, while the promoter region between -244 and 0 could individually display wounding-responsiveness, and the fragment from -467 to -244 was required for SA- and NaCl-inducible expression of the PtDrl02 promoter. Additionally, it was found that the -985 to -669 sequence was the ABA-responding promoter fragment. These results suggested that the PtDrl02 promoter was modulated by multiple cis-regulatory elements in distinct and complex patterns to regulate PtDrl02 gene expression. Our study also suggested that the PtDrl02 gene 5' untranslated region, as well as a Populus WRKY transcription factor, PtWRKY1, was involved in the regulation of PtDrl02 promoter activities.

Functional analysis of 5' untranslated region of a TIR-NBS-encoding gene from triploid white poplar.

Genome-wide analyses have identified a set of TIR-NBS-encoding genes in plants. However, the molecular mechanism underlying the expression of these genes is still unknown. In this study, we presented a TIR-NBS-encoding gene, PtDrl02, that displayed a low level of tissue-specific expression in a triploid white poplar [(Populus tomentosa x P. bolleana) x P. tomentosa], and analyzed the effects of the 5' untranslated region (UTR) on gene expression. The 5' UTR sequence repressed the reporter activity of beta-glucuronidase (GUS) gene under PtDrl02 promoter by 113.5-fold with a staining ratio of 2.97% in the transgenic tobacco plants. Quantitative RT-PCR assays revealed that the 5' UTR sequence decreased the transcript level of the GUS reporter gene by 13.3-fold, implying a regulatory role of 5' UTR in transcription and/or mRNA destabilization. The comparison of GUS activity with the transcript abundance indicated that the 5' UTR sequence decreased the translation efficiency of target gene by 88.3%. Additionally, the analysis of the transgenic P-985/UTRDelta/GUS plants showed that both the exon1 sequence and the leading intron within the 5' UTR region were responsible for the regulation of gene expression. Our results suggested a negative effect of the 5' UTR of PtDrl02 gene on gene expression.