Genome-Wide Identification of B-Box Gene Family and Candidate Light-Related Member Analysis of Tung Tree (Vernicia fordii).

Light is one of the most important environmental factors for plant growth. In the production process of tung oil tree cultivation, due to the inappropriate growth of shading conditions, the lower branches are often dry and dead, which seriously affects the yield of tung oil trees. However, little is known about the key factors of light-induced tree photomorphogenesis. In this study, a total of 22 VfBBX family members were identified to provide a reference for candidate genes in tung tree seedlings. All members of the VfBBX family have different numbers of highly conserved B-box domains or CCT domains. Phylogenetic evolution clustered the VfBBX genes into four categories, and the highest density of members was on chromosome 6. Interspecific collinearity analysis suggested that there were six pairs of duplicate genes in VfBBX members, but the expression levels of all family members in different growth and development stages of the tung tree were significantly divergent. After different degrees of shading treatment and physiological data determination of tung tree seedlings, the differential expression level and chlorophyll synthesis genes correlation analysis revealed that VfBBX9 was a typical candidate nuclear localization transcription factor that was significantly differentially expressed in light response. This study systematically identified the VfBBX gene family and provided a reference for studying its molecular function, enhanced the theoretical basis for tung tree breeding, and identified excellent varieties.

Tung Tree (Vernicia fordii) Genome Provides A Resource for Understanding Genome Evolution and Improved Oil Production.

Tung tree (Vernicia fordii) is an economically important woody oil plant that produces tung oil rich in eleostearic acid. Here, we report a high-quality chromosome-scale genome sequence of tung tree. The genome sequence was assembled by combining Illumina short reads, Pacific Biosciences single-molecule real-time long reads, and Hi-C sequencing data. The size of tung tree genome is 1.12 Gb, with 28,422 predicted genes and over 73% repeat sequences. The V. fordii underwent an ancient genome triplication event shared by core eudicots but no further whole-genome duplication in the subsequent ca. 34.55 million years of evolutionary history of the tung tree lineage. Insertion time analysis revealed that repeat-driven genome expansion might have arisen as a result of long-standing long terminal repeat retrotransposon bursts and lack of efficient DNA deletion mechanisms. The genome harbors 88 resistance genes encoding nucleotide-binding sites; 17 of these genes may be involved in early-infection stage of Fusarium wilt resistance. Further, 651 oil-related genes were identified, 88 of which are predicted to be directly involved in tung oil biosynthesis. Relatively few phosphoenolpyruvate carboxykinase genes, and synergistic effects between transcription factors and oil biosynthesis-related genes might contribute to the high oil content of tung seed. The tung tree genome constitutes a valuable resource for understanding genome evolution, as well as for molecular breeding and genetic improvements for oil production.

Molecular mechanism of the extended oil accumulation phase contributing to the high seed oil content for the genotype of tung tree (Vernicia fordii).

BACKGROUND: Oil from seeds of the tung tree (Vernicia fordii) has unique drying properties that are industrially important. We found that the extended oil accumulation period was related to the high seed oil content at maturity among tung tree population. In order to understand the molecular mechanism underlying the high oil content in tung tree seed, Tree H and L were adopted for the further investigation, with seed oil content of about 70 and 45%, respectively. We compared the transcriptomic changes of seed at various times during oil accumulation between the two trees. RESULTS: Transcriptomes analysis revealed that many genes involved in glycolysis, fatty acid synthesis, and tri-acyl glyceride assembly still kept high expression in the late period of seed oil accumulation for Tree H only. Many genes in fatty acid degradation pathway were largely up regulated in the late period of seed oil accumulation for Tree L only. Four transcription factors related to fatty acid biosynthesis had different expression pattern in the seed oil accumulation period for the two trees. WRI1 was down regulated and kept the low expression in the late period of seed oil accumulation for the two trees. PII, LEC1 and LEC1-LIKE extended the high expression in the late period of seed oil accumulation in Tree H only. CONCLUSIONS: The continued accumulation of oil in the late period of seed oil accumulation for Tree H was associated with relatively high expression of the relevant genes in glycolysis, fatty acid synthesis and tri-acyl glyceride assembly. PII, LEC1, and LEC1-LIKE rather than WRI1 should play an important role in the oil continual accumulation in the late period of seed oil accumulation in Tree H. This study provides novel insight into the variation in seed oil content and informs plant breeding strategies to maximize oil yield.

Proteomic Analysis of Tung Tree (Vernicia fordii) Oilseeds during the Developmental Stages.

The tung tree (Vernicia fordii), a non-model woody plant belonging to the Euphorbiaceae family, is a promising economic plant due to the high content of a novel high-value oil in its seeds. Many metabolic pathways are active during seed development. Oil (triacylglycerols (TAGs)) accumulates in oil bodies distributed in the endosperm cells of tung tree seeds. The relationship between oil bodies and oil content during tung tree seed development was analyzed using ultrastructural observations, which confirmed that oil accumulation was correlated with the volumes and numbers of oil bodies in the endosperm cells during three different developmental stages. For a deeper understanding of seed development, we carried out proteomic analyses. At least 144 proteins were differentially expressed during three different developmental stages. A total of 76 proteins were successfully identified using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry/mass spectrometry (MALDI-TOF/MS/MS). These proteins were grouped into 11 classes according to their functions. The major groups of differentially expressed proteins were associated with energy metabolism (25%), fatty acid metabolism (15.79%) and defense (14.47%). These results strongly suggested that a very high percentage of gene expression in seed development is dedicated to the synthesis and accumulation of TAGs.

Molecular cloning and expression analysis of a novel BCCP subunit gene from Aleurites moluccana.

Aleurites moluccana L. is grown as a roadside tree in southern China and the oil content of its seed is higher than other oil plants, such as Jatropha curcas and Camellia oleifera. A. moluccana is considered a promising energy plant because its seed oil could be used to produce biodiesel and bio-jet fuel. In addition, the bark, leaves, and kernels of A. moluccana have various medical and commercial uses. Here, a novel gene coding the biotin carboxyl carrier protein subunit (BCCP) was cloned from A. moluccana L. using the homology cloning method combined with rapid amplification of cDNA end (RACE) technology. The isolated full-length cDNA sequence (designated AM-accB) was 1188 bp, containing a 795-bp open reading frame coding for 265 amino acids. The deduced amino acid sequence of AM-accB contained a biotinylated domain located between amino acids 190 and 263. A. moluccana BCCP shows high identity at the amino acid level to its homologues in other higher plants, such as Vernicia fordii, J. curcas, and Ricinus communis (86, 77, and 70%, respectively), which all contain conserved domains for ACCase activity. The expression of the AM-accB gene during the middle stage of development and maturation in A. moluccana seeds was higher than that in early and later stages. The expression pattern of the AM-accB gene is very similar to that of the oil accumulation rate.

Fatty acid profile and unigene-derived simple sequence repeat markers in tung tree (Vernicia fordii).

Tung tree (Vernicia fordii) provides the sole source of tung oil widely used in industry. Lack of fatty acid composition and molecular markers hinders biochemical, genetic and breeding research. The objectives of this study were to determine fatty acid profiles and develop unigene-derived simple sequence repeat (SSR) markers in tung tree. Fatty acid profiles of 41 accessions showed that the ratio of α-eleostearic acid was increasing continuously with a parallel trend to the amount of tung oil accumulation while the ratios of other fatty acids were decreasing in different stages of the seeds and that α-eleostearic acid (18∶3) consisted of 77% of the total fatty acids in tung oil. Transcriptome sequencing identified 81,805 unigenes from tung cDNA library constructed using seed mRNA and discovered 6,366 SSRs in 5,404 unigenes. The di- and tri-nucleotide microsatellites accounted for 92% of the SSRs with AG/CT and AAG/CTT being the most abundant SSR motifs. Fifteen polymorphic genic-SSR markers were developed from 98 unigene loci tested in 41 cultivated tung accessions by agarose gel and capillary electrophoresis. Genbank database search identified 10 of them putatively coding for functional proteins. Quantitative PCR demonstrated that all 15 polymorphic SSR-associated unigenes were expressed in tung seeds and some of them were highly correlated with oil composition in the seeds. Dendrogram revealed that most of the 41 accessions were clustered according to the geographic region. These new polymorphic genic-SSR markers will facilitate future studies on genetic diversity, molecular fingerprinting, comparative genomics and genetic mapping in tung tree. The lipid profiles in the seeds of 41 tung accessions will be valuable for biochemical and breeding studies.

Developmental regulation of diacylglycerol acyltransferase family gene expression in tung tree tissues.

Diacylglycerol acyltransferases (DGAT) catalyze the final and rate-limiting step of triacylglycerol (TAG) biosynthesis in eukaryotic organisms. DGAT genes have been identified in numerous organisms. Multiple isoforms of DGAT are present in eukaryotes. We previously cloned DGAT1 and DGAT2 genes of tung tree (Vernicia fordii), whose novel seed TAGs are useful in a wide range of industrial applications. The objective of this study was to understand the developmental regulation of DGAT family gene expression in tung tree. To this end, we first cloned a tung tree gene encoding DGAT3, a putatively soluble form of DGAT that possesses 11 completely conserved amino acid residues shared among 27 DGAT3s from 19 plant species. Unlike DGAT1 and DGAT2 subfamilies, DGAT3 is absent from animals. We then used TaqMan and SYBR Green quantitative real-time PCR, along with northern and western blotting, to study the expression patterns of the three DGAT genes in tung tree tissues. Expression results demonstrate that 1) all three isoforms of DGAT genes are expressed in developing seeds, leaves and flowers; 2) DGAT2 is the major DGAT mRNA in tung seeds, whose expression profile is well-coordinated with the oil profile in developing tung seeds; and 3) DGAT3 is the major form of DGAT mRNA in tung leaves, flowers and immature seeds prior to active tung oil biosynthesis. These results suggest that DGAT2 is probably the major TAG biosynthetic isoform in tung seeds and that DGAT3 gene likely plays a significant role in TAG metabolism in other tissues. Therefore, DGAT2 should be a primary target for tung oil engineering in transgenic organisms.

Hydrophobic-domain-dependent protein-protein interactions mediate the localization of GPAT enzymes to ER subdomains.

The endoplasmic reticulum (ER) is a dynamic organelle that consists of numerous regions or 'subdomains' that have discrete morphological features and functional properties. Although it is generally accepted that these subdomains differ in their protein and perhaps lipid compositions, a clear understanding of how they are assembled and maintained has not been well established. We previously demonstrated that two diacylglycerol acyltransferase enzymes (DGAT1 and DGAT2) from tung tree (Vernicia fordii) were located in different subdomains of ER, but the mechanisms responsible for protein targeting to these subdomains were not elucidated. Here we extend these studies by describing two glycerol-3-phosphate acyltransferase-like (GPAT) enzymes from tung tree, GPAT8 and GPAT9, that both colocalize with DGAT2 in the same ER subdomains. Measurement of protein-protein interactions using the split-ubiquitin assay revealed that GPAT8 interacts with itself, GPAT9 and DGAT2, but not with DGAT1. Furthermore, mutational analysis of GPAT8 revealed that the protein's first predicted hydrophobic region, which contains an amphipathic helix-like motif, is required for interaction with DGAT2 and for DGAT2-dependent colocalization in ER subdomains. Taken together, these results suggest that the regulation and organization of ER subdomains is mediated at least in part by higher-ordered, hydrophobic-domain-dependent homo- and hetero-oligomeric protein-protein interactions.

Novel targeting signals mediate the sorting of different isoforms of the tail-anchored membrane protein cytochrome b5 to either endoplasmic reticulum or mitochondria.

Tail-anchored membrane proteins are a class of proteins that are targeted posttranslationally to various organelles and integrated by a single segment of hydrophobic amino acids located near the C terminus. Although the localization of tail-anchored proteins in specific subcellular compartments in plant cells is essential for their biological function, the molecular targeting signals responsible for sorting these proteins are not well defined. Here, we describe the biogenesis of four closely related tung (Aleurites fordii) cytochrome b5 isoforms (Cb5-A, -B, -C, and -D), which are small tail-anchored proteins that play an essential role in many cellular processes, including lipid biosynthesis. Using a combination of in vivo and in vitro assays, we show that Cb5-A, -B, and -C are targeted exclusively to the endoplasmic reticulum (ER), whereas Cb5-D is targeted specifically to mitochondrial outer membranes. Comprehensive mutational analyses of ER and mitochondrial Cb5s revealed that their C termini, including transmembrane domains (TMD) and tail regions, contained several unique physicochemical and sequence-specific characteristics that defined organelle-specific targeting motifs. Mitochondrial targeting of Cb5 was mediated by a combination of hydrophilic amino acids along one face of the TMD, an enrichment of branched beta-carbon-containing residues in the medial portion of the TMD, and a dibasic -R-R/K/H-x motif in the C-terminal tail. By contrast, ER targeting of Cb5 depended primarily upon the overall length and hydrophobicity of the TMD, although an -R/H-x-Y/F- motif in the tail was also a targeting determinant. Collectively, the results presented provide significant insight into the early biogenetic events required for entry of tail-anchored proteins into either the ER or mitochondrial targeting pathways.