Exploring transcription factors reveals crucial members and regulatory networks involved in different abiotic stresses in Brassica napus L

Deciphering peanut complex genomes paves a way to understand its origin and domestication

Peanut (Arachis) is a key oil and protein crop worldwide with large genome. The genomes of diploid and tetraploid peanuts have been sequenced, which were compared to decipher their genome structures, evolutionary, and life secrets. Genome sequencing efforts showed that different cultivars, although Bt homeologs being more privileged in gene retention and gene expression. This subgenome bias, extended to sequence variation and point mutation, might be related to the long terminal repeat (LTR) explosions after tetraploidization, especially in At subgenomes. Except that, whole-genome sequences revealed many important genes, for example, fatty acids and triacylglycerols pathway, NBS-LRR (nucleotide-binding site-leucine-rich repeats), and seed size decision genes, were enriched after recursive polyploidization. Each ancestral polyploidy, with old ones having occurred hundreds of thousand years ago, has thousands of duplicated genes in extant genomes, contributing to genetic novelty. Notably, although full genome sequences are available, the actual At subgenome ancestor has still been elusive, highlighted with new debate about peanut origin. Although being an orphan crop lagging behind other crops in genomic resources, the genome sequencing achievement has laid a solid foundation for advancing crop enhancement and system biology research of peanut.

Comprehensive in silico analysis of the underutilized crop tef (Eragrostis tef (Zucc.) Trotter) genome reveals drought tolerance signatures

BACKGROUND

Tef (Eragrostis tef) is a C4 plant known for its tiny, nutritious, and gluten-free grains. It contains higher levels of protein, vitamins, and essential minerals like calcium (Ca), iron (Fe), copper (Cu), and zinc (Zn) than common cereals. Tef is cultivated in diverse ecological zones under diverse climatic conditions. Studies have shown that tef has great diversity in withstanding environmental challenges such as drought. Drought is a major abiotic stress severely affecting crop productivity and becoming a bottleneck to global food security. Here, we used in silico-based functional genomic analysis to identify drought-responsive genes in tef and validated their expression using quantitative RT-PCR.

RESULTS

We identified about 729 drought-responsive genes so far reported in six crop plants, including rice, wheat, maize, barley, sorghum, pearl millet, and the model plant Arabidopsis, and reported 20 genes having high-level of GO terms related to drought, and significantly enriched in several biological and molecular function categories. These genes were found to play diverse roles, including water and fluid transport, resistance to high salt, cold, and drought stress, abscisic acid (ABA) signaling, de novo DNA methylation, and transcriptional regulation in tef and other crops. Our analysis revealed substantial differences in the conserved domains of some tef genes from well-studied rice orthologs. We further analyzed the expression of sixteen tef orthologs using quantitative RT-PCR in response to PEG-induced osmotic stress.

CONCLUSIONS

The findings showed differential regulation of some drought-responsive genes in shoots, roots, or both tissues. Hence, the genes identified in this study may be promising candidates for trait improvement in crops via transgenic or gene-editing technologies.

Transcriptomic and Metabolomic Analyses Reveal the Key Genes Related to Shade Tolerance in Soybean

Soybean (Glycine max) is an important crop, rich in proteins, vegetable oils and several other phytochemicals, which is often affected by light during growth. However, the specific regulatory mechanisms of leaf development under shade conditions have yet to be understood. In this study, the transcriptome and metabolome sequencing of leaves from the shade-tolerant soybean 'Nanxiadou 25' under natural light (ND1) and 50% shade rate (SHND1) were carried out, respectively. A total of 265 differentially expressed genes (DEGs) were identified, including 144 down-regulated and 121 up-regulated genes. Meanwhile, KEGG enrichment analysis of DEGs was performed and 22 DEGs were significantly enriched in the top five pathways, including histidine metabolism, riboflavin metabolism, vitamin B6 metabolism, glycerolipid metabolism and cutin, suberine and wax biosynthesis. Among all the enrichment pathways, the most DEGs were enriched in plant hormone signaling pathways with 19 DEGs being enriched. Transcription factors were screened out and 34 differentially expressed TFs (DETFs) were identified. Weighted gene co-expression network analysis (WGCNA) was performed and identified 10 core hub genes. Combined analysis of transcriptome and metabolome screened out 36 DEGs, and 12 potential candidate genes were screened out and validated by quantitative real-time polymerase chain reaction (qRT-PCR) assay, which may be related to the mechanism of shade tolerance in soybean, such as ATP phosphoribosyl transferase (ATP-PRT2), phosphocholine phosphatase (PEPC), AUXIN-RESPONSIVE PROTEIN (IAA17), PURPLE ACID PHOSPHATASE (PAP), etc. Our results provide new knowledge for the identification and function of candidate genes regulating soybean shade tolerance and provide valuable resources for the genetic dissection of soybean shade tolerance molecular breeding.

Identification of drought-tolerant hub genes in Iranian KC-2226 genotype of Aegilops tauschii using transcriptomic analysis

Aegilops tauschii, as a donor of D genome to the bread wheat with a valuable source of resistance to different biotic and abiotic stresses, is used to improve the quality of wheat cultivars. Every genotype has a specific genetic content, the investigation of which can lead to the identification of useful genes such as stress tolerance genes, including drought. Therefore, 23 genotypes of Ae. tauschii were selected to evaluate their morphological and physiological traits under greenhouse conditions. Among them, a superior tolerant genotype (KC-2226) was chosen for transcriptomic analysis. Our result showed that 5007 and 3489 genes were deferentially up- and downregulated, respectively. Upregulated genes were involved in photosynthesis, glycolysis/gluconeogenesis, and amino acid biosynthesis whereas downregulated genes were often engaged in DNA synthesis, replication, repair and topological changes. The result of protein-protein interaction network analysis showed that AT1G76550 (1.46), AT1G20950 (1.42), IAR4 (1.19), and PYD2 (1.16) among upregulated genes and THY-1 (44), PCNA1 (41) and TOPII (22) among down-regulated genes had the highest interactions with other genes. In conclusion, Ae. tauschii employs elevated transcription of specific genes involved in photosynthesis, glycolysis and gluconeogenesis and amino acid biosynthesis pathways rather than genes active in DNA synthesis and repair to provide the energy needed for the plant to survive under stress conditions.

Conservation and Divergence of the Trihelix Genes in Brassica and Expression Profiles of BnaTH Genes in Brassica napus under Abiotic Stresses

Trihelix (TH) proteins are a family of plant-specific transcription factors that play a role in light response and are extensively involved in plant growth and development, as well as in various stress responses. However, the function of TH genes in Brassica napus (B. napus) remains unclear, as does the evolution and differentiation pattern of TH genes in Brassica plants. Here, we identified a total of 455 TH genes in seven species, including six Brassica species and Arabidopsis, which were grouped into five clades, GT-1, GT-2, GTγ, SH4, and SIP1, each with 69, 142, 44, 55, and 145 members, respectively. The types and distributions of motifs of the TH proteins and the structures of the TH genes are conserved in the same subgroup, and some variations in certain amino acid residues occur in B. napus when inheriting motifs from Brassica rapa (B. rapa) and Brassica oleracea (B. oleracea). Collinearity analysis revealed that the massive expansion of TH genes in tetraploid species was attributed to the hetero-tetraploidization of diploid ancestors and gene duplication events within the tetraploid species. Comparative analysis of the membership numbers of five subgroups in different species revealed that the GT-2 and SIP1 genes underwent significant expansion during evolution, possibly to support the better adaptation of plants to their environments. The differential expression of the BnaTH genes under five stresses indicates that the BnaTH genes are involved in plant responses to stresses such as drought, cold, and heat. The presence of different stress-responsive cis-elements in the upstream promoter region of the genes indicated that BnaTH genes have the potential to cope with variable environments. Meanwhile, qRT-PCR analyses also confirmed that five TH genes respond to different abiotic stresses. Our results provide information and candidates for further studies on the role of TH genes in stress resistance of B. napus.

Identification and Characterization of the HD-Zip Gene Family and Dimerization Analysis of HB7 and HB12 in Brassica napus L

Homeodomain-leucine zipper (HD-Zip) genes encode plant-specific transcription factors, which play important roles in plant growth, development, and response to environmental stress. These genes have not been fully studied in allopolyploid Brassica napus, an important kind of oil crop. In this study, 165 HD-Zip genes were identified in B. napus and classified into four subfamilies. If proteins belong to the same subfamily, they exhibit similarities in gene structure, motifs, and domain distribution patterns. BnHD-Zip genes were unevenly distributed in the An and Cn subgenomes. Whole genome triplication (WGT) events may be major mechanisms accounting for this gene family expansion. Orthologous gene analysis showed that the process of this gene family expansion was accompanied by domain loss. We further found three genes homologous to HB7 and three genes homologous to HB12, all induced by PEG, ABA, and NaCl treatment. HB7 could not form homodimers but could form heterodimers with HB12 based on yeast two-hybrid assays. The results of this study provide valuable information for further exploration of the HD-Zip gene family in B. napus.

Exploring drought-responsive crucial genes in Sorghum

Drought severely affects global food production. Sorghum is a typical drought-resistant model crop. Based on RNA-seq data for Sorghum with multiple time points and the gray correlation coefficient, this paper firstly selects candidate genes via mean variance test and constructs weighted gene differential co-expression networks (WGDCNs); then, based on guilt-by-rewiring principle, the WGDCNs and the hidden Markov random field model, drought-responsive crucial genes are identified for five developmental stages respectively. Enrichment and sequence alignment analysis reveal that the screened genes may play critical functional roles in drought responsiveness. A multilayer differential co-expression network for the screened genes reveals that Sorghum is very sensitive to pre-flowering drought. Furthermore, a crucial gene regulatory module is established, which regulates drought responsiveness via plant hormone signal transduction, MAPK cascades, and transcriptional regulations. The proposed method can well excavate crucial genes through RNA-seq data, which have implications in breeding of new varieties with improved drought tolerance.

•

We design a method that unites gene rewiring network and Markov random field model

•

Drought-responsive genes for five developmental stages of Sorghum are explored

•

A multilayer network reveals that Sorghum is very sensitive to pre-flowering drought

•

A drought-responsive crucial gene regulatory module is established for Sorghum

Genome-Wide Identification and Expression Analysis of the Zinc Finger Protein Gene Subfamilies under Drought Stress in Triticum aestivum

The zinc finger protein (ZFP) family is one of plants’ most diverse family of transcription factors. These proteins with finger-like structural domains have been shown to play a critical role in plant responses to abiotic stresses such as drought. This study aimed to systematically characterize Triticum aestivum ZFPs (TaZFPs) and understand their roles under drought stress. A total of 9 TaC2H2, 38 TaC3HC4, 79 TaCCCH, and 143 TaPHD were identified, which were divided into 4, 7, 12, and 14 distinct subgroups based on their phylogenetic relationships, respectively. Segmental duplication dominated the evolution of four subfamilies and made important contributions to the large-scale amplification of gene families. Syntenic relationships, gene duplications, and Ka/Ks result consistently indicate a potential strong purifying selection on TaZFPs. Additionally, TaZFPs have various abiotic stress-associated cis-acting regulatory elements and have tissue-specific expression patterns showing different responses to drought and heat stress. Therefore, these genes may play multiple functions in plant growth and stress resistance responses. This is the first comprehensive genome-wide analysis of ZFP gene families in T. aestivum to elucidate the basis of their function and resistance mechanisms, providing a reference for precise manipulation of genetic engineering for drought resistance in T. aestivum.

Gene Differential Co-Expression Networks Based on RNA-Seq: Construction and Its Applications

Gene co-expression network (GCN) becomes an increasingly important tool in omics data analysis. A great challenge for GCN construction is that the sample size is far lower than the number of genes. Traditional methods rely on considerable samples. Moreover, association signals are likely weak, nonlinear and stochastic, which are difficult to be identified among thousands of candidates. In this paper, the gray correlation coefficient (GCC) is introduced, and a novel method to construct gene differential co-expression networks (GDCNs) is proposed. Based on the GDCNs, three measures are proposed to explore informative genes. The proposed method can make full use of the information provided by a handful of samples and overcome the shortages of GCNs, which can evaluate the changes of co-expression relationships that are possibly triggered by treatments. Based on RNA-seq data of Brassica napus, GDCNs under multiple experimental conditions are constructed and investigated. It is found that the GCC-based method is very robust to data processing. The GDCNs facilitate the inference of gene functions and the identification of informative genes that are responsible for stress responsiveness. The GDCN-based approaches integrate the 'guilt by association' and the 'guilt by rewiring' rules, which provide alternative tools for omics data analysis.

Integration of the Metabolomic and Transcriptome Analysis Reveals the Remarkable Compounds of G. bicolor Young and Mature Leaves under Different Iron Nutrient Conditions

Gynura bicolor (Roxb. ex Willd.) DC. (G. bicolor) is a functional vegetable rich in iron (Fe) and widely grown in Asia (e.g., Japan and China). Because most Fe in the soil exists in the form of insoluble oxides or hydroxides, it is difficult for plants to obtain Fe from the soil. A comparative metabolomic and transcriptome study was carried out to investigate the effect of Fe deficiency on metabolite synthesis and gene expression in young and mature leaves of G. bicolor. Fe deficiency caused chlorosis and decreased the chlorophyll content in young leaves. The metabolomic results for young leaves showed that l-glutamate and 4-hydroxybutanoic acid lactone significantly increased and decreased, respectively. The transcriptome results showed that the expression levels of genes involved in ferric reduction oxidase 7 and 14-kDa proline-rich protein DC2.15-like were significantly upregulated and downregulated, respectively. However, Fe deficiency had little effect on mature leaves.

Genetic and Physiological Responses to Heat Stress in Brassica napus

Given the current rise in global temperatures, heat stress has become a major abiotic challenge affecting the growth and development of various crops and reducing their productivity. Brassica napus, the second largest source of vegetable oil worldwide, experiences a drastic reduction in seed yield and quality in response to heat. This review outlines the latest research that explores the genetic and physiological impact of heat stress on different developmental stages of B. napus with a special attention to the reproductive stages of floral progression, organogenesis, and post flowering. Several studies have shown that extreme temperature fluctuations during these crucial periods have detrimental effects on the plant and often leading to impaired growth and reduced seed production. The underlying mechanisms of heat stress adaptations and associated key regulatory genes are discussed. Furthermore, an overview and the implications of the polyploidy nature of B. napus and the regulatory role of alternative splicing in forming a priming-induced heat-stress memory are presented. New insights into the dynamics of epigenetic modifications during heat stress are discussed. Interestingly, while such studies are scarce in B. napus, opposite trends in expression of key genetic and epigenetic components have been identified in different species and in cultivars within the same species under various abiotic stresses, suggesting a complex role of these genes and their regulation in heat stress tolerance mechanisms. Additionally, omics-based studies are discussed with emphasis on the transcriptome, proteome and metabolome of B. napus, to gain a systems level understanding of how heat stress alters its yield and quality traits. The combination of omics approaches has revealed crucial interactions and regulatory networks taking part in the complex machinery of heat stress tolerance. We identify key knowledge gaps regarding the impact of heat stress on B. napus during its yield determining reproductive stages, where in-depth analysis of this subject is still needed. A deeper knowledge of heat stress response components and mechanisms in tissue specific models would serve as a stepping-stone to gaining insights into the regulation of thermotolerance that takes place in this important crop species and support future breeding of heat tolerant crops.

Identification of WRKY transcription factors responding to abiotic stresses in Brassica napus L

A total of 278 BnWRKYs were identified and analyzed. Ectopic expression of BnWRKY149 and BnWRKY217 suggests that they function in the ABA signaling pathway. WRKY transcription factors play an important role in plant development, however, their function in Brassica napus L. abiotic stress response is still unclear. In this study, a total of 278 BnWRKY transcription factors were identified from the B. napus genome data, and they were subsequently distributed in three main groups. The protein motifs and classification of BnWRKY transcription factors were analyzed, and the locations of their corresponding encoding genes were mapped on the chromosomes of B. napus. Transcriptome analysis of rapeseed seedlings exposed to drought, salt, heat, cold and abscisic acid treatment revealed that 99 BnWRKYs responded to at least one of these stresses. The expression profiles of 12 BnWRKYs were examined with qPCR and the result coincided with RNA-seq analysis. Two genes of interest, BnWRKY149 and BnWRKY217 (homologs of AtWRKY40), were overexpressed in Arabidopsis, and the corresponding proteins were located to the nucleus. Transgene plants of BnWRKY149 and BnWRKY217 were less sensitive to ABA than Arabidopsis Col-0 plants, suggesting they might play important roles in the responses of rapeseed to abiotic stress.

Statistical Identification of Important Nodes in Biological Systems

Biological systems can be modeled and described by biological networks. Biological networks are typical complex networks with widely real-world applications. Many problems arising in biological systems can be boiled down to the identification of important nodes. For example, biomedical researchers frequently need to identify important genes that potentially leaded to disease phenotypes in animal and explore crucial genes that were responsible for stress responsiveness in plants. To facilitate the identification of important nodes in biological systems, one needs to know network structures or behavioral data of nodes (such as gene expression data). If network topology was known, various centrality measures can be developed to solve the problem; while if only behavioral data of nodes were given, some sophisticated statistical methods can be employed. This paper reviewed some of the recent works on statistical identification of important nodes in biological systems from three aspects, that is, 1) in general complex networks based on complex networks theory and epidemic dynamic models; 2) in biological networks based on network motifs; and 3) in plants based on RNA-seq data. The identification of important nodes in a complex system can be seen as a mapping from the system to the ranking score vector of nodes, such mapping is not necessarily with explicit form. The three aspects reflected three typical approaches on ranking nodes in biological systems and can be integrated into one general framework. This paper also proposed some challenges and future works on the related topics. The associated investigations have potential real-world applications in the control of biological systems, network medicine and new variety cultivation of crops.

Omics: The way forward to enhance abiotic stress tolerance in Brassica napus L

Plant abiotic stresses negative affects growth and development, causing a massive reduction in global agricultural production. Rapeseed (Brassica napus L.) is a major oilseed crop because of its economic value and oilseed production. However, its productivity has been reduced by many environmental adversities. Therefore, it is a prime need to grow rapeseed cultivars, which can withstand numerous abiotic stresses. To understand the various molecular and cellular mechanisms underlying the abiotic stress tolerance and improvement in rapeseed, omics approaches have been extensively employed in recent years. This review summarized the recent advancement in genomics, transcriptomics, proteomics, metabolomics, and their imploration in abiotic stress regulation in rapeseed. Some persisting bottlenecks have been highlighted, demanding proper attention to fully explore the omics tools. Further, the potential prospects of the CRISPR/Cas9 system for genome editing to assist molecular breeding in developing abiotic stress-tolerant rapeseed genotypes have also been explained. In short, the combination of integrated omics, genome editing, and speed breeding can alter rapeseed production worldwide.

Statistical Identification of Important Nodes in Biological Systems

Biological systems can be modeled and described by biological networks. Biological networks are typical complex networks with widely real-world applications. Many problems arising in biological systems can be boiled down to the identification of important nodes. For example, biomedical researchers frequently need to identify important genes that potentially leaded to disease phenotypes in animal and explore crucial genes that were responsible for stress responsiveness in plants. To facilitate the identification of important nodes in biological systems, one needs to know network structures or behavioral data of nodes (such as gene expression data). If network topology was known, various centrality measures can be developed to solve the problem; while if only behavioral data of nodes were given, some sophisticated statistical methods can be employed. This paper reviewed some of the recent works on statistical identification of important nodes in biological systems from three aspects, that is, 1) in general complex networks based on complex networks theory and epidemic dynamic models; 2) in biological networks based on network motifs; and 3) in plants based on RNA-seq data. The identification of important nodes in a complex system can be seen as a mapping from the system to the ranking score vector of nodes, such mapping is not necessarily with explicit form. The three aspects reflected three typical approaches on ranking nodes in biological systems and can be integrated into one general framework. This paper also proposed some challenges and future works on the related topics. The associated investigations have potential real-world applications in the control of biological systems, network medicine and new variety cultivation of crops.

The Role of Chloroplast Gene Expression in Plant Responses to Environmental Stress

Chloroplasts are plant organelles that carry out photosynthesis, produce various metabolites, and sense changes in the external environment. Given their endosymbiotic origin, chloroplasts have retained independent genomes and gene-expression machinery. Most genes from the prokaryotic ancestors of chloroplasts were transferred into the nucleus over the course of evolution. However, the importance of chloroplast gene expression in environmental stress responses have recently become more apparent. Here, we discuss the emerging roles of the distinct chloroplast gene expression processes in plant responses to environmental stresses. For example, the transcription and translation of psbA play an important role in high-light stress responses. A better understanding of the connection between chloroplast gene expression and environmental stress responses is crucial for breeding stress-tolerant crops better able to cope with the rapidly changing environment.

Transcriptomic analysis reveals the changes of energy production and AsA-GSH cycle in oat embryos during seed ageing

Deterioration during seed storage generally causes seed vigour declining. However, the mechanism of deterioration occurred still not clear. Seeds and embryos of oat (Avena sativa L.) were selected to analyze the relation of physiological and metabolic reactions with DEGs by using RNA-seq. Oat seed vigour declined during seeds aged 0 day (CK), 16 days (CD16) and 32 days (CD32). The changes of MDA and H₂O₂ contents, antioxidant enzymes activities of APX, DHAR, MDHAR and GR related with AsA-GSH cycle in embryos illustrated that seed vigour declined to the minimum at CD32. Transcriptomic analysis showed a total of 11335 and 8274 DEGs were identified at CD16 and CD32 compared with CK respectively, of which 4070 were overlapped. When seed vigour declined to the moderate level (CD16), the accumulation of H₂O₂ caused by the inhibition of complex I in ETC could be alleviated with AsA-GSH cycle. RNA-seq and qRT-PCR results both showed alternative oxidase in alternate respiratory pathway was upregulated which would maintain seed respiration. However, as seed vigour was at the lowest level (CD32), blocked ETC caused by down-regulation of complex III, including Ubiquinol-cytochrome C reductase complex 14kD subunit and Ubiquinol-cytochrome C reductase, UQCRX/QCR9 like, were more seriously and H₂O₂ scavenging was limited by the inactive AsA-GSH cycle. It could be suggested that the function of AsA-GSH would play a key role for regulating the physiological responses of ETC in embryos during seed ageing. These results would provide an insight into embryo for the transcriptomic information during oat seed ageing.

Screening of abiotic stress-responsive cotton genes using a cotton full-length cDNA overexpressing Arabidopsis library

Cotton (Gossypium hirsutum L.) is a major crop and the main source of natural fiber worldwide. Because various abiotic and biotic stresses strongly influence cotton fiber yield and quality, improved stress resistance of this crop plant is urgently needed. In this study, we used Gateway technology to construct a normalized full-length cDNA overexpressing (FOX) library from upland cotton cultivar ZM12 under various stress conditions. The library was transformed into Arabidopsis to produce a cotton-FOX-Arabidopsis library. Screening of this library yielded 6,830 transgenic Arabidopsis lines, of which 757 were selected for sequencing to ultimately obtain 659 cotton ESTs. GO and KEGG analyses mapped most of the cotton ESTs to plant biological process, cellular component, and molecular function categories. Next, 156 potential stress-responsive cotton genes were identified from the cotton-FOX-Arabidopsis library under drought, salt, ABA, and other stress conditions. Four stress-related genes identified from the library, designated as GhCAS, GhAPX, GhSDH, and GhPOD, were cloned from cotton complementary DNA, and their expression patterns under stress were analyzed. Phenotypic experiments indicated that overexpression of these cotton genes in Arabidopsis affected the response to abiotic stress. The method developed in this study lays a foundation for high-throughput cloning and rapid identification of cotton functional genes.

Integrated transcriptomics and metabolomics analysis of catechins, caffeine and theanine biosynthesis in tea plant (Camellia sinensis) over the course of seasons

BACKGROUND

Catechins, caffeine, and theanine as three important metabolites in the tea leaves play essential roles in the formation of specific taste and shows potential health benefits to humans. However, the knowledge on the dynamic changes of these metabolites content over seasons, as well as the candidate regulatory factors, remains largely undetermined.

RESULTS

An integrated transcriptomic and metabolomic approach was used to analyze the dynamic changes of three mainly metabolites including catechins, caffeine, and theanine, and to explore the potential influencing factors associated with these dynamic changes over the course of seasons. We found that the catechins abundance was higher in Summer than that in Spring and Autumn, and the theanine abundance was significantly higher in Spring than that in Summer and Autumn, whereas caffeine exhibited no significant changes over three seasons. Transcriptomics analysis suggested that genes in photosynthesis pathway were significantly down-regulated which might in linkage to the formation of different phenotypes and metabolites content in the tea leaves of varied seasons. Fifty-six copies of nine genes in catechins biosynthesis, 30 copies of 10 genes in caffeine biosynthesis, and 12 copies of six genes in theanine biosynthesis were detected. The correlative analysis further presented that eight genes can be regulated by transcription factors, and highly correlated with the changes of metabolites abundance in tea-leaves.

CONCLUSION

Sunshine intensity as a key factor can affect photosynthesis of tea plants, further affect the expression of major Transcription factors (TFs) and structural genes in, and finally resulted in the various amounts of catechins, caffeine and theaine in tea-leaves over three seasons. These findings provide new insights into abundance and influencing factors of metabolites of tea in different seasons, and further our understanding in the formation of flavor, nutrition and medicinal function.

MYB43 in Oilseed Rape (Brassica napus) Positively Regulates Vascular Lignification, Plant Morphology and Yield Potential but Negatively Affects Resistance to Sclerotinia sclerotiorum

Arabidopsis thaliana MYB43 (AtMYB43) is suggested to be involved in cell wall lignification. PtrMYB152, the Populus orthologue of AtMYB43, is a transcriptional activator of lignin biosynthesis and vessel wall deposition. In this research, MYB43 genes from Brassica napus (rapeseed) and its parental species B. rapa and B. oleracea were molecularly characterized, which were dominantly expressed in stem and other vascular organs and showed responsiveness to Sclerotinia sclerotiorum infection. The BnMYB43 family was silenced by RNAi, and the transgenic rapeseed lines showed retardation in growth and development with smaller organs, reduced lodging resistance, fewer silique number and lower yield potential. The thickness of the xylem layer decreased by 28%; the numbers of sclerenchymatous cells, vessels, interfascicular fibers, sieve tubes and pith cells in the whole cross section of the stem decreased by 28%, 59%, 48%, 34% and 21% in these lines, respectively. The contents of cellulose and lignin decreased by 17.49% and 16.21% respectively, while the pectin content increased by 71.92% in stems of RNAi lines. When inoculated with S. sclerotiorum, the lesion length was drastically decreased by 52.10% in the stems of transgenic plants compared with WT, implying great increase in disease resistance. Correspondingly, changes in the gene expression patterns of lignin biosynthesis, cellulose biosynthesis, pectin biosynthesis, cell cycle, SA- and JA-signals, and defensive pathways were in accordance with above phenotypic modifications. These results show that BnMYB43, being a growth-defense trade-off participant, positively regulates vascular lignification, plant morphology and yield potential, but negatively affects resistance to S. sclerotiorum. Moreover, this lignification activator influences cell biogenesis of both lignified and non-lignified tissues of the whole vascular organ.

Comparative Transcriptome Analysis of Two Aegilops tauschii with Contrasting Drought Tolerance by RNA-Seq

As the diploid progenitor of common wheat, Aegilops tauschii is considered to be a valuable resistance source to various biotic and abiotic stresses. However, little has been reported concerning the molecular mechanism of drought tolerance in Ae. tauschii. In this work, the drought tolerance of 155 Ae. tauschii accessions was firstly screened on the basis of their coleoptile lengths under simulated drought stress. Subsequently, two accessions (XJ002 and XJ098) with contrasting coleoptile lengths were selected and intensively analyzed on rate of water loss (RWL) as well as physiological characters, confirming the difference in drought tolerance at the seedling stage. Further, RNA-seq was utilized for global transcriptome profiling of the two accessions seedling leaves under drought stress conditions. A total of 6969 differentially expressed genes (DEGs) associated with drought tolerance were identified, and their functional annotations demonstrated that the stress response was mediated by pathways involving alpha-linolenic acid metabolism, starch and sucrose metabolism, peroxisome, mitogen-activated protein kinase (MAPK) signaling, carbon fixation in photosynthetic organisms, and glycerophospholipid metabolism. In addition, DEGs with obvious differences between the two accessions were intensively analyzed, indicating that the expression level of DEGs was basically in alignment with the physiological changes of Ae. tauschii under drought stress. The results not only shed fundamental light on the regulatory process of drought tolerance in Ae. tauschii, but also provide a new gene resource for improving the drought tolerance of common wheat.

Arabidopsis GDSL1 overexpression enhances rapeseed Sclerotinia sclerotiorum resistance and the functional identification of its homolog in Brassica napus

Sclerotinia stem rot (SSR) caused by Sclerotinia sclerotiorum is a devastating disease of rapeseed (Brassica napus L.). To date, the genetic mechanisms of rapeseed' interactions with S. sclerotiorum are not fully understood, and molecular-based breeding is still the most effective control strategy for this disease. Here, Arabidopsis thaliana GDSL1 was characterized as an extracellular GDSL lipase gene functioning in Sclerotinia resistance. Loss of AtGDSL1 function resulted in enhanced susceptibility to S. sclerotiorum. Conversely, overexpression of AtGDSL1 in B. napus enhanced resistance, which was associated with increased reactive oxygen species (ROS) and salicylic acid (SA) levels, and reduced jasmonic acid levels. In addition, AtGDSL1 can cause an increase in lipid precursor phosphatidic acid levels, which may lead to the activation of downstream ROS/SA defence-related pathways. However, the rapeseed BnGDSL1 with highest sequence similarity to AtGDSL1 had no effect on SSR resistance. A candidate gene association study revealed that only one AtGDSL1 homolog from rapeseed, BnaC07g35650D (BnGLIP1), significantly contributed to resistance traits in a natural B. napus population, and the resistance function was also confirmed by a transient expression assay in tobacco leaves. Moreover, genomic analyses revealed that BnGLIP1 locus was embedded in a selected region associated with SSR resistance during the breeding process, and its elite allele type belonged to a minor allele in the population. Thus, BnGLIP1 is the functional equivalent of AtGDSL1 and has a broad application in rapeseed S. sclerotiorum-resistance breeding.

Transcriptional analyses reveal the molecular mechanism governing shade tolerance in the invasive plant Solidago canadensis

Solidago canadensis is an invasive plant that is capable of adapting to variable light conditions. To elucidate the shade tolerance mechanism in S. canadensis at the molecular level, transcriptome analyses were performed for leaves growing under natural light and three shade level conditions. Many differentially expressed genes (DEGs) were found in the comparative analysis, including those involved in photosynthesis, antioxidant, and secondary metabolism of phenol- and flavonoid-related pathways. Most genes encoding proteins involved in photosynthesis, such as photosystem I reaction center subunit (Psa), photosystem II core complex protein (Psb), and light-harvesting chlorophyll protein (Lhca and Lhcb), and reactive oxygen species (ROS) scavenging-related enzymes, such as superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), were upregulated with the shade levels. Furthermore, most of the DEGs related to secondary metabolite synthesis were also upregulated in the shade conditions. Our study indicates that S. canadensis can respond to shade stress by modulating the expression of several photosynthesis-related, free radical scavenging-related, and secondary metabolism-related genes; thus, this species has the ability to adapt to different light conditions.

Transcriptome analysis reveals key genes involved in the regulation of nicotine biosynthesis at early time points after topping in tobacco (Nicotiana tabacum L.)

BACKGROUND

Nicotiana tabacum is an important economic crop. Topping, a common agricultural practice employed with flue-cured tobacco, is designed to increase leaf nicotine contents by increasing nicotine biosynthesis in roots. Many genes are found to be differentially expressed in response to topping, particularly genes involved in nicotine biosynthesis, but comprehensive analyses of early transcriptional responses induced by topping are not yet available. To develop a detailed understanding of the mechanisms regulating nicotine biosynthesis after topping, we have sequenced the transcriptomes of Nicotiana tabacum roots at seven time points following topping.

RESULTS

Differential expression analysis revealed that 4830 genes responded to topping across all time points. Amongst these, nine gene families involved in nicotine biosynthesis and two gene families involved in nicotine transport showed significant changes during the immediate 24 h period following topping. No obvious preference to the parental species was detected in the differentially expressed genes (DEGs). Significant changes in transcript levels of nine genes involved in nicotine biosynthesis and phytohormone signal transduction were validated by qRT-PCR assays. 549 genes encoding transcription factors (TFs), found to exhibit significant changes in gene expression after topping, formed 15 clusters based on similarities of their transcript level time-course profiles. 336 DEGs involved in phytohormone signal transduction, including genes functionally related to the phytohormones jasmonic acid, abscisic acid, auxin, ethylene, and gibberellin, were identified at the earliest time point after topping.

CONCLUSIONS

Our research provides the first detailed analysis of the early transcriptional responses to topping in N. tabacum, and identifies excellent candidates for further detailed studies concerning the regulation of nicotine biosynthesis in tobacco roots.

Comprehensive analysis of bZIP transcription factors uncovers their roles during dimorphic floret differentiation and stress response in Cleistogenes songorica

BACKGROUND

Transcription factors act as important regulators of transcription networks. Basic leucine zipper (bZIP) transcription factors have been shown to be involved in multiple biological processes in plants. However, no information is available for the bZIP family in Cleistogenes songorica, which is an important xerophytic and allotetraploid grass in desert grasslands.

RESULTS

In this study, 86 CsbZIPs were identified in the allotetraploid C. songorica genome. For location analysis, CsbZIPs were distributed evenly across two subgenomes of C. songorica. Phylogenetic tree analysis among three species indicated that CsbZIPs were evolutionarily more closely related to OsbZIPs than AtbZIPs. Syntenic and phylogenetic analyses confirmed that the CsbZIPs were mainly expanded by whole-genome duplication events. Furthermore, it was determined that rice and C. songorica might have undergone purified selection during their long evolutionary history by calculating the Ks values and Ka/Ks ratios of orthologous gene pairs. By analysing the expression patterns of CsbZIPs in different tissues and under abiotic stresses, 21 CsbZIP genes were differentially expressed between chasmogamous (CH) and cleistogamous (CL) flowers, including two FLOWERING LOCUS D (FD) genes. In shoots and roots, 79.1 and 87.2% of the CsbZIP genes, respectively, displayed transcript changes under at least one stress treatment, such as heat, cold, drought and salt. Strikingly, 17 common CsbZIP genes showed differential expression under stress response and during CL flowering. Co-expression network, GO annotation and real-time quantitative reverse transcription PCR (qRT-PCR) analyses revealed a close relationship between CL flowering-associated genes and abiotic stress-related genes.

CONCLUSIONS

BZIP TFs were comprehensively analysed and identified in allotetraploid C. songorica. Our results provide insights into the evolutionary history of the bZIP family in C. songorica and provide abiotic stress-responsive and CL-associated candidate CsbZIP genes for potential applications in the genetic improvement of plants.

Differential Expression of Circular RNAs in Polytocous and Monotocous Uterus during the Reproductive Cycle of Sheep

Simple Summary

The uterus is an important reproductive organ that provides nutrition and place for embryonic development. In this study, we identified circular RNAs by deep sequencing and analyzed their expression in the uteri of polytocous and monotocous sheep (FecB++) during follicular and luteal phases. Gene Ontology (GO) and KEGG enrichment analyses revealed that the source genes of these differential circular RNAs (circRNAs) were mainly enriched in reproductive hormone- and energy metabolism-related pathways. These results provide information on the molecular mechanisms of sheep prolificacy.

Abstract

CircRNA plays important roles in cell proliferation, differentiation, autophagy and apoptosis during development. However, there are few reports on circRNAs related to livestock reproduction. In this study, we identified circRNAs by deep sequencing and analyzed their expression in the uteri of polytocous and monotocous sheep (FecB++) during follicular and luteal phases. There were 147 and 364 circRNAs with differential expression in the follicular and luteal phases, respectively. GO and KEGG enrichment analysis was performed for the host genes of the circRNAs to predict the functions of differentially expressed circRNAs. These source genes were mainly involved in the estrogen signaling pathway, TGFβ signaling pathway, GnRH signaling pathway, oxytocin signaling pathway, pentose phosphate pathway, and starch and sucrose metabolism related to reproduction and energy metabolism. CircRNA expression patterns were validated by RT-qPCR. Our findings provide a solid foundation for the identification and characterization of key important circRNAs involved in reproduction.

Identification of candidate genes controlling oil content by combination of genome-wide association and transcriptome analysis in the oilseed crop Brassica napus

BACKGROUND

Increasing seed oil content is one of the most important targets for rapeseed (Brassica napus) breeding. However, genetic mechanisms of mature seed oil content in Brassica napus (B. napus) remain little known. To identify oil content-related genes, a genome-wide association study (GWAS) was performed using 588 accessions.

RESULTS

High-throughput genome resequencing resulted in 385,692 high-quality single nucleotide polymorphism (SNPs) with a minor allele frequency (MAF) > 0.05. We identified 17 loci that were significantly associated with seed oil content, among which 12 SNPs were distributed on the A3 (11 loci) and A1 (one loci) chromosomes, and five novel significant SNPs on the C5 (one loci) and C7 (four loci) chromosomes, respectively. Subsequently, we characterized differentially expressed genes (DEGs) between the seeds and silique pericarps on main florescences and primary branches of extremely high- and low-oil content accessions (HO and LO). A total of 64 lipid metabolism-related DEGs were identified, 14 of which are involved in triacylglycerols (TAGs) biosynthesis and assembly. Additionally, we analyzed differences in transcription levels of key genes involved in de novo fatty acid biosynthesis in the plastid, TAGs assembly and lipid droplet packaging in the endoplasmic reticulum (ER) between high- and low-oil content B. napus accessions.

CONCLUSIONS

The combination of GWAS and transcriptome analyses revealed seven candidate genes located within the confidence intervals of significant SNPs. Current findings provide valuable information for facilitating marker-based breeding for higher seed oil content in B. napus.

Genome-wide identification and evolutionary analysis of TGA transcription factors in soybean

The gain of function in genes and gene families is a continuous process and is a key factor in understanding gene and genome evolution in plants. TGACG-Binding (TGA) transcription factors (TFs) have long been known for their essential roles in plant defence in Arabidopsis, but their roles in legume symbiosis are yet to be explored. Here, we identified a total of 25 TGA (named GmTGA1-GmTGA25) genes in soybean. Through phylogenetic analysis, we discovered a clade of GmTGA proteins that appear to be legume-specific. Among them, two GmTGAs were unique by possessing the autophagy sequence in their proteins, while the third one was an orphan gene in soybean. GmTGAs were structurally different from AtTGAs, and their expression patterns also differed with the dominant expression of AtTGAs and GmTGAs in aerial and underground parts, respectively. Moreover, twenty-five GmTGAs showed a strong correlation among the gene expression in roots, nodules, and root hairs. The qRT-PCR analysis results revealed that among 15 tested GmTGAs, six were induced and four were suppressed by rhizobia inoculation, while 11 of these GmTGAs were induced by high nitrate. Our findings suggested the important roles of GmTGAs in symbiotic nodulation and in response to nitrogen availability in soybean.