Co-expression analysis aids in the identification of genes in the cuticular wax pathway in maize

A genome-wide association study identifies genes associated with cuticular wax metabolism in maize

The plant cuticle is essential in plant defense against biotic and abiotic stresses. To systematically elucidate the genetic architecture of maize (Zea mays L.) cuticular wax metabolism, 2 cuticular wax-related traits, the chlorophyll extraction rate (CER) and water loss rate (WLR) of 389 maize inbred lines, were investigated and a genome-wide association study (GWAS) was performed using 1.25 million single nucleotide polymorphisms (SNPs). In total, 57 nonredundant quantitative trait loci (QTL) explaining 5.57% to 15.07% of the phenotypic variation for each QTL were identified. These QTLs contained 183 genes, among which 21 strong candidates were identified based on functional annotations and previous publications. Remarkably, 3 candidate genes that express differentially during cuticle development encode β-ketoacyl-CoA synthase (KCS). While ZmKCS19 was known to be involved in cuticle wax metabolism, ZmKCS12 and ZmKCS3 functions were not reported. The association between ZmKCS12 and WLR was confirmed by resequencing 106 inbred lines, and the variation of WLR was significant between different haplotypes of ZmKCS12. In this study, the loss-of-function mutant of ZmKCS12 exhibited wrinkled leaf morphology, altered wax crystal morphology, and decreased C32 wax monomer levels, causing an increased WLR and sensitivity to drought. These results confirm that ZmKCS12 plays a vital role in maize C32 wax monomer synthesis and is critical for drought tolerance. In sum, through GWAS of 2 cuticular wax-associated traits, this study reveals comprehensively the genetic architecture in maize cuticular wax metabolism and provides a valuable reference for the genetic improvement of stress tolerance in maize.

Trade-offs between the accumulation of cuticular wax and jasmonic acid-mediated herbivory resistance in maize

Plants have evolved complex physical and chemical defense systems that allow them to withstand herbivory infestation. Composed of a complex mixture of very-long-chain fatty acids (VLCFAs) and their derivatives, cuticular wax constitutes the first physical line of defense against herbivores. Here, we report the function of Glossy 8 (ZmGL8), which encodes a 3-ketoacyl reductase belonging to the fatty acid elongase complex, in orchestrating wax production and jasmonic acid (JA)-mediated defenses against herbivores in maize (Zea mays). The mutation of GL8 enhanced chemical defenses by activating the JA-dependent pathway. We observed a trade-off between wax accumulation and JA levels across maize glossy mutants and 24 globally collected maize inbred lines. In addition, we demonstrated that mutants defective in cuticular wax biosynthesis in Arabidopsis thaliana and maize exhibit enhanced chemical defenses. Comprehensive transcriptomic and lipidomic analyses indicated that the gl8 mutant confers chemical resistance to herbivores by remodeling VLCFA-related lipid metabolism and subsequent JA biosynthesis and signaling. These results suggest that VLCFA-related lipid metabolism has a critical role in regulating the trade-offs between cuticular wax and JA-mediated chemical defenses.

MdBT2 regulates nitrogen-mediated cuticular wax biosynthesis via a MdMYB106-MdCER2L1 signalling pathway in apple

Cuticular waxes play important roles in plant development and the interaction between plants and their environment. Researches on wax biosynthetic pathways have been reported in several plant species. Also, wax formation is closely related to environmental condition. However, the regulatory mechanism between wax and environmental factors, especially essential mineral elements, is less studied. Here we found that nitrogen (N) played a negative role in the regulation of wax synthesis in apple. We therefore analysed wax content, composition and crystals in BTB-TAZ domain protein 2 (MdBT2) overexpressing and antisense transgenic apple seedlings and found that MdBT2 could downregulate wax biosynthesis. Furthermore, R2R3-MYB transcription factor 16-like protein (MdMYB106) interacted with MdBT2, and MdBT2 mediated its ubiquitination and degradation through the 26S proteasome pathway. Finally, HXXXD-type acyl-transferase ECERIFERUM 2-like1 (MdCER2L1) was confirmed as a downstream target gene of MdMYB106. Our findings reveal an N-mediated apple wax biosynthesis pathway and lay a foundation for further study of the environmental factors associated with wax regulatory networks in apple.

Bacterium-enabled transient gene activation by artificial transcription factors for resolving gene regulation in maize

Understanding gene regulatory networks is essential to elucidate developmental processes and environmental responses. Here, we studied regulation of a maize (Zea mays) transcription factor gene using designer transcription activator-like effectors (dTALes), which are synthetic Type III TALes of the bacterial genus Xanthomonas and serve as inducers of disease susceptibility gene transcription in host cells. The maize pathogen Xanthomonas vasicola pv. vasculorum was used to introduce 2 independent dTALes into maize cells to induced expression of the gene glossy3 (gl3), which encodes a MYB transcription factor involved in biosynthesis of cuticular wax. RNA-seq analysis of leaf samples identified, in addition to gl3, 146 genes altered in expression by the 2 dTALes. Nine of the 10 genes known to be involved in cuticular wax biosynthesis were upregulated by at least 1 of the 2 dTALes. A gene previously unknown to be associated with gl3, Zm00001d017418, which encodes aldehyde dehydrogenase, was also expressed in a dTALe-dependent manner. A chemically induced mutant and a CRISPR-Cas9 mutant of Zm00001d017418 both exhibited glossy leaf phenotypes, indicating that Zm00001d017418 is involved in biosynthesis of cuticular waxes. Bacterial protein delivery of dTALes proved to be a straightforward and practical approach for the analysis and discovery of pathway-specific genes in maize.

Weighted gene co-expression network analysis of nitrogen (N)-responsive genes and the putative role of G-quadruplexes in N use efficiency (NUE) in rice

Rice is an important target to improve crop nitrogen (N) use efficiency (NUE), and the identification and shortlisting of the candidate genes are still in progress. We analyzed data from 16 published N-responsive transcriptomes/microarrays to identify, eight datasets that contained the maximum number of 3020 common genes, referred to as N-responsive genes. These include different classes of transcription factors, transporters, miRNA targets, kinases and events of post-translational modifications. A Weighted gene co-expression network analysis (WGCNA) with all the 3020 N-responsive genes revealed 15 co-expression modules and their annotated biological roles. Protein-protein interaction network analysis of the main module revealed the hub genes and their functional annotation revealed their involvement in the ubiquitin process. Further, the occurrences of G-quadruplex sequences were examined, which are known to play important roles in epigenetic regulation but are hitherto unknown in N-response/NUE. Out of the 3020 N-responsive genes studied, 2298 contained G-quadruplex sequences. We compared these N-responsive genes containing G-quadruplex sequences with the 3601 genes we previously identified as NUE-related (for being both N-responsive and yield-associated). This analysis revealed 389 (17%) NUE-related genes containing G-quadruplex sequences. These genes may be involved in the epigenetic regulation of NUE, while the rest of the 83% (1811) genes may regulate NUE through genetic mechanisms and/or other epigenetic means besides G-quadruplexes. A few potentially important genes/processes identified as associated with NUE were experimentally validated in a pair of rice genotypes contrasting for NUE. The results from the WGCNA and G4 sequence analysis of N-responsive genes helped identify and shortlist six genes as candidates to improve NUE. Further, the hitherto unavailable segregation of genetic and epigenetic gene targets could aid in informed interventions through genetic and epigenetic means of crop improvement.

Identification of Key Modules and Candidate Genes for Powdery Mildew Resistance of Wheat-Agropyron cristatum Translocation Line WAT-2020-17-6 by WGCNA

As one of the serious diseases of wheat, powdery mildew (Blumeria graminis f. sp. tritici) is a long-term threat to wheat production. Therefore, it is of great significance to explore new powdery mildew-resistant genes for breeding. The wild relative species of wheat provide gene resources for resistance to powdery mildew breeding. Agropyron cristatum (2n = 4x = 28, genomes PPPP) is an important wild relative of wheat, carrying excellent genes for high yield, disease resistance, and stress resistance, which can be used for wheat improvement. To understand the molecular mechanism of powdery mildew resistance in the wheat-A. cristatum translocation line WAT2020-17-6, transcriptome sequencing was performed, and the resistance genes were analyzed by weighted gene co-expression network analysis (WGCNA). In the results, 42,845 differentially expressed genes were identified and divided into 18 modules, of which six modules were highly correlated with powdery mildew resistance. Gene ontology (GO) enrichment analysis showed that the six interested modules related to powdery mildew resistance were significantly enriched in N-methyltransferase activity, autophagy, mRNA splicing via spliceosome, chloroplast envelope, and AMP binding. The candidate hub genes of the interested modules were further identified, and their regulatory relationships were analyzed based on co-expression data. The temporal expression pattern of the 12 hub genes was verified within 96 h after powdery mildew inoculation by RT-PCR assay. In this study, we preliminarily explained the resistance mechanism of the wheat-A. cristatum translocation lines and obtained the hub candidate genes, which laid a foundation in the exploration of resistance genes in A. cristatum for powdery mildew-resistant breeding in wheat.

Characterization of Cuticular Wax in Tea Plant and Its Modification in Response to Low Temperature

Cuticular wax ubiquitously covers the outer layer of plants and protects them against various abiotic and biotic stresses. Nevertheless, the characteristics of cuticular wax and its role in cold resistance in tea plants remain unclear. In our study, cuticular wax from different tissues, cultivars, and leaves during different spatio-temporal growth stages were characterized and compared in tea plants. The composition, distribution pattern, and structural profile of cuticular wax showed considerable tissue specificity, particularly in petals and seeds. During the spatial development of tea leaves, total wax content increased from the first to fifth leaf in June, while a decreasing pattern was observed in September. Additionally, the total wax content and number of wax compounds were enhanced, and the wax composition significantly varied with leaf growth from June to September. Ten cultivars showed considerable differences in total wax content and composition, such as the predominance of saturated fatty acids and primary alcohols in SYH and HJY cultivars, respectively. Correlation analysis suggested that n-hexadecanoic acid is positively related to cold resistance in tea plants. Further transcriptome analysis from cold-sensitive AJBC, cold-tolerant CYQ, and EC 12 cultivars indicated that the inducible expression of wax-related genes was associated with the cold tolerance of different cultivars in response to cold stress. Our results revealed the characterization of cuticular wax in tea plants and provided new insights into its modification in cold tolerance.

Stocky1, a Novel Gene Involved in Maize Seedling Development and Cuticle Integrity

The cuticle is the plant’s outermost layer that covers the surfaces of aerial parts. This structure is composed of a variety of aliphatic molecules and is well-known for its protective role against biotic and abiotic stresses in plants. Mutants with a permeable cuticle show developmental defects such as organ fusions and altered seed germination and viability. In this study, we identified a novel maize mutant, stocky1, with unique features: lethal at the seedling stage, and showing a severely dwarfed phenotype, due to a defective cuticle. For the first time, the mutant was tentatively mapped to chromosome 5, bin 5.04. The mutant phenotype investigated in this work has the potential to contribute to the elucidation of the role of the cuticle during plant development. The possibility of controlling this trait is of relevance in the context of climate change, as it may contribute to tolerance to abiotic stresses.

Transcriptome and Resequencing Analyses Provide Insight into Differences in Organic Acid Accumulation in Two Pear Varieties

Fruit acidity is one of the main determinants of fruit flavor and a target trait in fruit breeding. However, the genomic mechanisms governing acidity variation among different pear varieties remain poorly understood. In this study, two pear varieties with contrasting organic acid levels, ‘Dangshansuli’ (low-acidity) and ‘Amute’ (high-acidity), were selected, and a combination of transcriptome and population genomics analyses were applied to characterize their patterns of gene expression and genetic variation. Based on RNA-seq data analysis, differentially expressed genes (DEGs) involved in organic acid metabolism and accumulation were identified. Weighted correlation network analysis (WGCNA) revealed that nine candidate TCA (tricarboxylic acid)-related DEGs and three acid transporter-related DEGs were located in three key modules. The regulatory networks of the above candidate genes were also predicted. By integrating pear resequencing data, two domestication-related genes were found to be upregulated in ‘Amute’, and this trend was further validated for other pear varieties with high levels of organic acid, suggesting distinct selective sweeps during pear dissemination and domestication. Collectively, this study provides insight into organic acid differences related to expression divergence and domestication in two pear varieties, pinpointing several candidate genes for the genetic manipulation of acidity in pears.

Gene-Metabolite Network Analysis Revealed Tissue-Specific Accumulation of Therapeutic Metabolites in Mallotus japonicus

Mallotus japonicus is a valuable traditional medicinal plant in East Asia for applications as a gastrointestinal drug. However, the molecular components involved in the biosynthesis of bioactive metabolites have not yet been explored, primarily due to a lack of omics resources. In this study, we established metabolome and transcriptome resources for M. japonicus to capture the diverse metabolite constituents and active transcripts involved in its biosynthesis and regulation. A combination of untargeted metabolite profiling with data-dependent metabolite fragmentation and metabolite annotation through manual curation and feature-based molecular networking established an overall metabospace of M. japonicus represented by 2129 metabolite features. M. japonicus de novo transcriptome assembly showed 96.9% transcriptome completeness, representing 226,250 active transcripts across seven tissues. We identified specialized metabolites biosynthesis in a tissue-specific manner, with a strong correlation between transcripts expression and metabolite accumulations in M. japonicus. The correlation- and network-based integration of metabolome and transcriptome datasets identified candidate genes involved in the biosynthesis of key specialized metabolites of M. japonicus. We further used phylogenetic analysis to identify 13 C-glycosyltransferases and 11 methyltransferases coding candidate genes involved in the biosynthesis of medicinally important bergenin. This study provides comprehensive, high-quality multi-omics resources to further investigate biological properties of specialized metabolites biosynthesis in M. japonicus.

Identification of biological pathway and process regulators using sparse partial least squares and triple-gene mutual interaction

Identification of biological process- and pathway-specific regulators is essential for advancing our understanding of regulation and formation of various phenotypic and complex traits. In this study, we applied two methods, triple-gene mutual interaction (TGMI) and Sparse Partial Least Squares (SPLS), to identify the regulators of multiple metabolic pathways in Arabidopsis thaliana and Populus trichocarpa using high-throughput gene expression data. We analyzed four pathways: (1) lignin biosynthesis pathway in A. thaliana and P. trichocarpa; (2) flavanones, flavonol and anthocyannin biosynthesis in A. thaliana; (3) light reaction pathway and Calvin cycle in A. thaliana. (4) light reaction pathway alone in A. thaliana. The efficiencies of two methods were evaluated by examining the positive known regulators captured, the receiver operating characteristic (ROC) curves and the area under ROC curves (AUROC). Our results showed that TGMI is in general more efficient than SPLS in identifying true pathway regulators and ranks them to the top of candidate regulatory gene lists, but the two methods are to some degree complementary because they could identify some different pathway regulators. This study identified many regulators that potentially regulate the above pathways in plants and are valuable for genetic engineering of these pathways.

Conserved co-functional network between maize and Arabidopsis aid in the identification of seed defective genes in maize

BACKGROUND

The biological pathways related to Arabidopsis seed development have been well studied and functional genes involved in it have been discovered. However, functional studies about maize seed development were more limited compared to Arabidopsis.

OBJECTIVE

Therefore, transferring knowledge from Arabidopsis into maize would facilitate functional studies about maize seed development.

METHOD

In this study, public transcriptome data of the two species related to seed development were obtained. Co-expression network in each species was compared by integrating orthology information.

RESULTS

This conserved co-functional network contained 4510 maize and 4808 Arabidopsis genes, respectively. Most of these genes were expressed in throughout embryo, early or later endosperm/seed. These conserved co-functional genes were significantly enriched for members of PPR protein family, which was consistent with that PPR proteins play an important role in maize seed development. Spatial-temporally co-functional genes were discovered in the seed coat and embryo. Furthermore, 66 well-studied genes involved in Arabidopsis seed development were co-functional with 319 maize genes and one maize gene (GRMZM2G036050) was further confirmed using an EMS-induced seed defective mutant by bulked segregating RNA sequencing (BSR) analysis.

CONCLUSIONS

Altogether, these results showed the potential of this approach to support functional studies in maize seed development by transferring knowledge from Arabidopsis.

The FUSED LEAVES1-ADHERENT1 regulatory module is required for maize cuticle development and organ separation

All aerial epidermal cells in land plants are covered by the cuticle, an extracellular hydrophobic layer that provides protection against abiotic and biotic stresses and prevents organ fusion during development. Genetic and morphological analysis of the classic maize adherent1 (ad1) mutant was combined with genome-wide binding analysis of the maize MYB transcription factor FUSED LEAVES1 (FDL1), coupled with transcriptional profiling of fdl1 mutants. We show that AD1 encodes an epidermally-expressed 3-KETOACYL-CoA SYNTHASE (KCS) belonging to a functionally uncharacterized clade of KCS enzymes involved in cuticular wax biosynthesis. Wax analysis in ad1 mutants indicates that AD1 functions in the formation of very-long-chain wax components. We demonstrate that FDL1 directly binds to CCAACC core motifs present in AD1 regulatory regions to activate its expression. Over 2000 additional target genes of FDL1, including many involved in cuticle formation, drought response and cell wall organization, were also identified. Our results identify a regulatory module of cuticle biosynthesis in maize that is conserved across monocots and eudicots, and highlight previously undescribed factors in lipid metabolism, transport and signaling that coordinate organ development and cuticle formation.

ZmSRL5 is involved in drought tolerance by maintaining cuticular wax structure in maize

Cuticular wax is a natural barrier on terrestrial plant organs, which protects plants from damages caused by a variety of stresses. Here, we report the identification and functional characterization of a cuticular-wax-related gene, Zea mays L. SEMI-ROLLED LEAF 5 (ZmSRL5). The loss-of-function mutant srl5, which was created by a 3,745 bp insertion in the first intron that led to the premature transcript, exhibited abnormal wax crystal morphology and distribution, which, in turn, caused the pleiotropic phenotypes including increased chlorophyll leaching and water loss rate, decreased leaf temperature, sensitivity to drought, as well as semi-rolled mature leaves. However, total wax amounts showed no significant difference between wild type and semi-rolled leaf5 (srl5) mutant. The phenotype of srl5 was confirmed through the generation of two allelic mutants using CRISPR/Cas9. ZmSRL5 encodes a CASPARIAN-STRIP-MEMBRANE-DOMAIN-LIKE (CASPL) protein located in plasma membrane, and highly expressed in developing leaves. Further analysis showed that the expressions of the most wax related genes were not affected or slightly altered in srl5. This study, thus, primarily uncovers that ZmSRL5 is required for the structure formation of the cuticular wax and could increase the drought tolerance by maintaining the proper cuticular wax structure in maize.

Synthetic Biology of Cannabinoids and Cannabinoid Glucosides in Nicotiana benthamiana and Saccharomyces cerevisiae

Phytocannabinoids are a group of plant-derived metabolites that display a wide range of psychoactive as well as health-promoting effects. The production of pharmaceutically relevant cannabinoids relies on extraction and purification from cannabis (Cannabis sativa) plants yielding the major constituents, Δ9-tetrahydrocannabinol and cannabidiol. Heterologous biosynthesis of cannabinoids in Nicotiana benthamiana or Saccharomyces cerevisiae may provide cost-efficient and rapid future production platforms to acquire pure and high quantities of both the major and the rare cannabinoids as well as novel derivatives. Here, we used a meta-transcriptomic analysis of cannabis to identify genes for aromatic prenyltransferases of the UbiA superfamily and chalcone isomerase-like (CHIL) proteins. Among the aromatic prenyltransferases, CsaPT4 showed CBGAS activity in both N. benthamiana and S. cerevisiae. Coexpression of selected CsaPT pairs and of CHIL proteins encoding genes with CsaPT4 did not affect CBGAS catalytic efficiency. In a screen of different plant UDP-glycosyltransferases, Stevia rebaudiana SrUGT71E1 and Oryza sativa OsUGT5 were found to glucosylate olivetolic acid, cannabigerolic acid, and Δ9-tetrahydrocannabinolic acid. Metabolic engineering of N. benthamiana for production of cannabinoids revealed intrinsic glucosylation of olivetolic acid and cannabigerolic acid. S. cerevisiae was engineered to produce olivetolic acid glucoside and cannabigerolic acid glucoside.

Gene co-expression network analysis to identify critical modules and candidate genes of drought-resistance in wheat

AIM

To establish a gene co-expression network for identifying principal modules and hub genes that are associated with drought resistance mechanisms, analyzing their mechanisms, and exploring candidate genes.

METHODS AND FINDINGS

42 data sets including PRJNA380841 and PRJNA369686 were used to construct the co-expression network through weighted gene co-expression network analysis (WGCNA). A total of 1,896,897,901 (284.30 Gb) clean reads and 35,021 differentially expressed genes (DEGs) were obtained from 42 samples. Functional enrichment analysis indicated that photosynthesis, DNA replication, glycolysis/gluconeogenesis, starch and sucrose metabolism, arginine and proline metabolism, and cell cycle were significantly influenced by drought stress. Furthermore, the DEGs with similar expression patterns, detected by K-means clustering, were grouped into 29 clusters. Genes involved in the modules, such as dark turquoise, yellow, and brown, were found to be appreciably linked with drought resistance. Twelve central, greatly correlated genes in stage-specific modules were subsequently confirmed and validated at the transcription levels, including TraesCS7D01G417600.1 (PP2C), TraesCS5B01G565300.1 (ERF), TraesCS4A01G068200.1 (HSP), TraesCS2D01G033200.1 (HSP90), TraesCS6B01G425300.1 (RBD), TraesCS7A01G499200.1 (P450), TraesCS4A01G118400.1 (MYB), TraesCS2B01G415500.1 (STK), TraesCS1A01G129300.1 (MYB), TraesCS2D01G326900.1 (ALDH), TraesCS3D01G227400.1 (WRKY), and TraesCS3B01G144800.1 (GT).

CONCLUSIONS

Analyzing the response of wheat to drought stress during different growth stages, we have detected three modules and 12 hub genes that are associated with drought resistance mechanisms, and five of those genes are newly identified for drought resistance. The references provided by these modules will promote the understanding of the drought-resistance mechanism. In addition, the candidate genes can be used as a basis of transgenic or molecular marker-assisted selection for improving the drought resistance and increasing the yields of wheat.

Non-homology-based prediction of gene functions in maize (Zea mays ssp. mays)

Advances in genome sequencing and annotation have eased the difficulty of identifying new gene sequences. Predicting the functions of these newly identified genes remains challenging. Genes descended from a common ancestral sequence are likely to have common functions. As a result, homology is widely used for gene function prediction. This means functional annotation errors also propagate from one species to another. Several approaches based on machine learning classification algorithms were evaluated for their ability to accurately predict gene function from non-homology gene features. Among the eight supervised classification algorithms evaluated, random-forest-based prediction consistently provided the most accurate gene function prediction. Non-homology-based functional annotation provides complementary strengths to homology-based annotation, with higher average performance in Biological Process GO terms, the domain where homology-based functional annotation performs the worst, and weaker performance in Molecular Function GO terms, the domain where the accuracy of homology-based functional annotation is highest. GO prediction models trained with homology-based annotations were able to successfully predict annotations from a manually curated "gold standard" GO annotation set. Non-homology-based functional annotation based on machine learning may ultimately prove useful both as a method to assign predicted functions to orphan genes which lack functionally characterized homologs, and to identify and correct functional annotation errors which were propagated through homology-based functional annotations.

Genome-Wide Analysis of Glycoside Hydrolase Family 1 β-glucosidase Genes in Brassica rapa and Their Potential Role in Pollen Development

Glycoside hydrolase family 1 (GH1) β-glucosidases (BGLUs) are encoded by a large number of genes, and are involved in many developmental processes and stress responses in plants. Due to their importance in plant growth and development, genome-wide analyses have been conducted in model plants (Arabidopsis and rice) and maize, but not in Brassica species, which are important vegetable crops. In this study, we systematically analyzed B. rapaBGLUs (BrBGLUs), and demonstrated the involvement of several genes in pollen development. Sixty-four BrBGLUs were identified in Brassica databases, which were anchored onto 10 chromosomes, with 10 tandem duplications. Phylogenetic analysis revealed that 64 genes were classified into 10 subgroups, and each subgroup had relatively conserved intron/exon structures. Clustering with Arabidopsis BGLUs (AtBGLUs) facilitated the identification of several important subgroups for flavonoid metabolism, the production of glucosinolates, the regulation of abscisic acid (ABA) levels, and other defense-related compounds. At least six BrBGLUs might be involved in pollen development. The expression of BrBGLU10/AtBGLU20, the analysis of co-expressed genes, and the examination of knocked down Arabidopsis plants strongly suggests that BrBGLU10/AtBGLU20 has an indispensable function in pollen development. The results that are obtained from this study may provide valuable information for the further understanding of β-glucosidase function and Brassica breeding, for nutraceuticals-rich Brassica crops.