PlnTFDB: updated content and new features of the plant transcription factor database

Computational analysis and expression profiling of NAC transcription factor family involved in biotic stress response in Manihot esculenta

The Nascent polypeptide-Associated Complex (NAC) family is among the largest plant-specific TF families and plays an important role in plant growth, development, and stress responses. NAC TFs have been extensively studied in plants such as rice and Arabidopsis; however, their characterization, functions, evolution, and expression patterns in Manihot esculenta (cassava) under environmental stress remain largely unexplored. Here, we used bioinformatic analyses and biotic stress responses to investigate the physicochemical properties, chromosome location, phylogeny, gene structure, expression patterns, and cis-elements in promoter regions of the NAC TFs in cassava. We identified 119 M. esculenta NAC (MeNAC) gene families, unevenly distributed on 16 chromosomes. We investigated expression patterns of all identified MeNAC TFs under Xanthomonas axonopodis pv. manihotis (Xam) infection, strain CHN11, at different time points. Only 20 MeNAC TFs showed expression of significant bacterial resistance. Six MeNACs (MeNAC7, 26, 63, 65, 77, and 113) were selected for functional analysis. qRT-PCR assays revealed that MeNAC7, 26, 63, 65, 77, and 113 were induced in response to XamCHN11 infection and may participate in the molecular interaction of cassava and bacterial blight. Interestingly, MeNAC26, MeNAC63, MeNAC65, and MeNAC113 responded to XamCHN11 infection at 3 h post-inoculation. Furthermore, we identified 13 stress-related cis-elements in promoter regions of the MeNAC genes that are involved in diverse environmental stress responses. Phylogenetic analysis revealed that MeNAC genes with similar structures and motif distributions were grouped. This study provides valuable insights into the evolution, diversity, and characterization of MeNAC TFs. It lays the groundwork for a better understanding of their biological roles and molecular mechanisms in cassava.

Genome-wide characterization and expression profiling of FARL (FHY3/FAR1) family genes in Zea mays

A significant role of the plant is played by the transcription factor FARL, which is light signal transduction as well as plant growth and development. Despite being transposases, FARL has developed a variety of dominant biological actions in evolution and speciation. On the other hand, little is known about the Zea mays FARL protein family. This study identifies and characterizes fifteen ZmFARL genes genome-wide, and RNA sequencing data was used to profile their expression. 105 FARL proteins from five plant species were classified into five groups based on sequence alignment and phylogeny. The ZmFARL genes' exon-intron and motif distribution were conserved based on their evolutionary group. The fifteen ZmFARL genes were distributed over seven of the ten Z. mays chromosomes, although no duplication was discovered. Cis-element analysis reveals that ZmFARL genes play a variety of activities, including tissue-specific, stress- and hormone-responsive expressions. Furthermore, the results of the RNA sequencing used to profile expression showed that the genes ZmFARL2 and ZmFARL5 were much more expressed than other genes in various tissues, particularly in leaf characteristics. The identification of likely genes involved in cellular activity in Z. mays and related species will be aided by the characterization of the FARL genes.

Global transcriptome analysis identifies critical functional modules associated with multiple abiotic stress responses in microalgae Chromochloris zofingiensis

In the current study, systems biology approach was applied to get a deep insight regarding the regulatory mechanisms of Chromochloris zofingiensis under overall stress conditions. Meta-analysis was performed using p-values combination of differentially expressed genes. To identify the informative models related to stress conditions, two distinct weighted gene co-expression networks were constructed and preservation analyses were performed using medianRankand Zsummary algorithms. Moreover, functional enrichment analysis of non-preserved modules was performed to shed light on the biological performance of underlying genes in the non-preserved modules. In the next step, the gene regulatory networks between top hub genes of non-preserved modules and transcription factors were inferred using ensemble of trees algorithm. Results showed that the power of beta = 7 was the best soft-thresholding value to ensure a scale-free network, leading to the determination of 12 co-expression modules with an average size of 128 genes. Preservation analysis showed that the connectivity pattern of the six modules including the blue, black, yellow, pink, greenyellow, and turquoise changed during stress condition which defined as non-preserved modules. Examples of enriched pathways in non-preserved modules were Oxidative phosphorylation", "Vitamin B6 metabolism", and "Arachidonic acid metabolism". Constructed regulatory network between identified TFs and top hub genes of non-preserved module such as Cz06g10250, Cz03g12130 showed that some specific TFs such as C3H and SQUAMOSA promoter binding protein (SBP) specifically regulates the specific hubs. The current findings add substantially to our understanding of the stress responsive underlying mechanism of C. zofingiensis for future studies and metabolite production programs.

Integrated transcriptomic and metabolomic analysis reveals the molecular profiles of dynamic variation in Lilium brownii var. viridulum suffering from bulb rot

Lilium brownii var. viridulum, known as Longya lily, is a well-known medicinal and edible plant in China. Bulb rot is a common disease in Longya lily cultivation that severely affects the yield and quality of lilies. According field investigations, we found that different Longya lily plants in the same field had different degrees of resistance to root rot. To find the reasons leading to the difference, we performed metabolomic and transcriptomic analyses of Longya lily with different degrees of disease. The transcriptomic analyses showed that the number of differentially expressed genes increased in early and mid-stage infections (LYBH2 and LYBH3), while decreased in late-stage infection (LYBH4). A total of 2309 DEGs showed the same expression trend in diseased bulb compared healthy bulb (LYBH1). The transcription factors (TFs) analysis of DEGs showed that several common TFs, like WRKY, bHLH, AP2/ERF-ERF and MYB, were significantly activated in bulbs after decay. The metabolomic analyses showed that there were 794 differentially accumulated metabolites, and metabolites with significant changes in relative content largely were phenolic acids, followed by flavonoids and amino acids and derivatives. The combined analysis of transcriptome and metabolome indicated that phenylpropanoid biosynthesis pathway was crucial in Longya lily resistance to bulb rot. Therefore, we speculated that the different degree of resistance to bulb rot in Longya lily may be related to the transcript levels of gene and contents of metabolites in the phenylpropanoid biosynthesis pathway. Overall, these results elucidate the molecular responses of Longya lily to bulb rot and lay a theoretical foundation for breeding resistant varieties.

Identification and core gene-mining of Weighted Gene Co-expression Network Analysis-based co-expression modules related to flood resistance in quinoa seedlings

BACKGROUND

As an emerging food crop with high nutritional value, quinoa has been favored by consumers in recent years; however, flooding, as an abiotic stress, seriously affects its growth and development. Currently, reports on the molecular mechanisms related to quinoa waterlogging stress responses are lacking; accordingly, the core genes related to these processes were explored via Weighted Gene Co-expression Network Analysis (WGCNA).

RESULTS

Based on the transcriptome data, WGCNA was used to construct a co-expression network of weighted genes associated with flooding resistance-associated physiological traits and metabolites. Here, 16 closely related co-expression modules were obtained, and 10 core genes with the highest association with the target traits were mined from the two modules. Functional annotations revealed the biological processes and metabolic pathways involved in waterlogging stress, and four candidates related to flooding resistance, specifically AP2/ERF, MYB, bHLH, and WRKY-family TFs, were also identified.

CONCLUSIONS

These results provide clues to the identification of core genes for quinoa underlying quinoa waterlogging stress responses. This could ultimately provide a theoretical foundation for breeding new quinoa varieties with flooding tolerance.

Genome-wide identification and role of HSFs in antioxidant response of hot water treated zucchini fruit during cold storage

Zucchini squashes are cold-sensitive and vulnerable to chilling injury (CI) resulting from reactive oxygen species (ROS) and hot water (HW) immersing effectively reduce CI symptoms during cold storage. However, mechanism involved in reduced ROS due to HW treatment has not been characterized well. In this study, tender green zucchini fruit were treated with HW for 15 min at 45 ± 1 °C and stored for 15 d at 4 ± 1 °C and above 90 % relative humidity. Results showed substantial reduction in CI index, electrolyte leakage, malonaldehyde (MDA) contents and ROS accumulation along with increased activity of ROS-scavenging enzymes due to HW treatment. To gain insight into the molecular mechanism involved in antioxidant defense system, transcriptomic analysis revealed that heat shock factors (HSF) accumulated due to HW treatment regulated the ROS pathway during cold stress. CpHSFA4a was one of the highly expressed transcription factors (TF) due to HW treatment that regulated the transcription of ROS enzymes related genes. CpHSFA4a bind actively with heat shock element (HSE) in promoter regions of CpSOD, CpCAT, CpAPX1, CpAPX2, and CpAPX3, activated and increased the expression of these genes. In conclusion, HW treatment alleviated the CI by maintaining ROS homeostasis through CpHSFA4a mediated ROS pathway in zucchini squashes during cold storage.

In silico analysis identified bZIP transcription factors genes responsive to abiotic stress in Alfalfa (Medicago sativa L.)

BACKGROUND

Alfalfa (Medicago sativa L.) is the most cultivated forage legume around the world. Under a variety of growing conditions, forage yield in alfalfa is stymied by biotic and abiotic stresses including heat, salt, drought, and disease. Given the sessile nature of plants, they use strategies including, but not limited to, differential gene expression to respond to environmental cues. Transcription factors control the expression of genes that contribute to or enable tolerance and survival during periods of stress. Basic-leucine zipper (bZIP) transcription factors have been demonstrated to play a critical role in regulating plant growth and development as well as mediate the responses to abiotic stress in several species, including Arabidopsis thaliana, Oryza sativa, Lotus japonicus and Medicago truncatula. However, there is little information about bZIP transcription factors in cultivated alfalfa.

RESULT

In the present study, 237 bZIP genes were identified in alfalfa from publicly available sequencing data. Multiple sequence alignments showed the presence of intact bZIP motifs in the identified sequences. Based on previous phylogenetic analyses in A. thaliana, alfalfa bZIPs were similarly divided and fell into 10 groups. The physico-chemical properties, motif analysis and phylogenetic study of the alfalfa bZIPs revealed high specificity within groups. The differential expression of alfalfa bZIPs in a suite of tissues indicates that bZIP genes are specifically expressed at different developmental stages in alfalfa. Similarly, expression analysis in response to ABA, cold, drought and salt stresses, indicates that a subset of bZIP genes are also differentially expressed and likely play a role in abiotic stress signaling and/or tolerance. RT-qPCR analysis on selected genes further verified these differential expression patterns.

CONCLUSIONS

Taken together, this work provides a framework for the future study of bZIPs in alfalfa and presents candidate bZIPs involved in stress-response signaling.

In-Depth Characterization of bZIP Genes in the Context of Endoplasmic Reticulum (ER) Stress in Brassica campestris ssp. chinensis

Numerous studies have been conducted to investigate the genomic characterization of bZIP genes and their involvement in the cellular response to endoplasmic reticulum (ER) stress. These studies have provided valuable insights into the coordinated cellular response to ER stress, which is mediated by bZIP transcription factors (TFs). However, a comprehensive and systematic investigations regarding the role of bZIP genes and their involvement in ER stress response in pak choi is currently lacking in the existing literature. To address this knowledge gap, the current study was initiated to elucidate the genomic characteristics of bZIP genes, gain insight into their expression patterns during ER stress in pak choi, and investigate the protein-to-protein interaction of bZIP genes with the ER chaperone BiP. In total, 112 members of the BcbZIP genes were identified through a comprehensive genome-wide analysis. Based on an analysis of sequence similarity, gene structure, conserved domains, and responsive motifs, the identified BcbZIP genes were categorized into 10 distinct subfamilies through phylogenetic analysis. Chromosomal location and duplication events provided insight into their genomic context and evolutionary history. Divergence analysis estimated their evolutionary history with a predicted divergence time ranging from 0.73 to 80.71 million years ago (MYA). Promoter regions of the BcbZIP genes were discovered to exhibit a wide variety of cis-elements, including light, hormone, and stress-responsive elements. GO enrichment analysis further confirmed their roles in the ER unfolded protein response (UPR), while co-expression network analysis showed a strong relationship of BcbZIP genes with ER-stress-responsive genes. Moreover, gene expression profiles and protein-protein interaction with ER chaperone BiP further confirmed their roles and capacity to respond to ER stress in pak choi.

Evaluation of aphid resistance on different rose cultivars and transcriptome analysis in response to aphid infestation

BACKGROUND

The rose is one of the most important ornamental flowers in the world for its aesthetic beauty but can be attacked by many pests such as aphids. Aphid infestation causes tremendous damage on plant tissues leading to harmed petals and leaves. Rose cultivars express different levels of resistance to aphid infestation yet the information remains unclear. Not only that, studies about the transcriptional analysis on defending mechanisms against aphids in rose are limited so far.

RESULTS

In this study, the aphid resistance of 20 rose cultivars was evaluated, and they could be sorted into six levels based on the number ratio of aphids. And then, a transcriptome analysis was conducted after aphid infestation in one high resistance (R, Harmonie) and one highly susceptibility (S, Carefree Wonder) rose cultivar. In open environment the majority of rose cultivars had the highest aphid number at May 6th or May 15th in 2020 and the resistance to infestation could be classified into six levels. Differential expression analysis revealed that there were 1,626 upregulated and 767 downregulated genes in the R cultivar and 481 upregulated and 63 downregulated genes in the S cultivar after aphid infestation. Pathway enrichment analysis of the differentially expressed genes revealed that upregulated genes in R and S cultivars were both enriched in defense response, biosynthesis of secondary metabolites (phenylpropanoid, alkaloid, and flavonoid), carbohydrate metabolism (galactose, starch, and sucrose metabolism) and lipid processing (alpha-linolenic acid and linolenic acid metabolism) pathways. In the jasmonic acid metabolic pathway, linoleate 13S-lipoxygenase was specifically upregulated in the R cultivar, while genes encoding other crucial enzymes, allene oxide synthase, allene oxide cyclase, and 12-oxophytodienoate reductase were upregulated in both cultivars. Transcription factor analysis and transcription factor binding search showed that WRKY transcription factors play a pivotal role during aphid infestation in the R cultivar.

CONCLUSIONS

Our study indicated the potential roles of jasmonic acid metabolism and WRKY transcription factors during aphid resistance in rose, providing clues for future research.

Revealing the specific regulations of nitric oxide on the postharvest ripening and senescence of bitter melon fruit

Bitter melon fruit is susceptible to yellowing, softening, and rotting under room-temperature storage conditions, resulting in reduced commercial value. Nitric oxide (NO) is an important signaling molecule and plays a crucial role in regulating the fruit postharvest quality. In this study, we investigated the effects of NO treatment on changes in sensory and firmness of bitter melon fruit during postharvest storage. Moreover, transcriptomic, metabolomic, and proteomic analyses were performed to elucidate the regulatory mechanisms through which NO treatment delays the ripening and senescence of bitter melon fruit. Our results show that differentially expressed genes (DEGs) were involved in fruit texture (CSLE, β-Gal, and PME), plant hormone signal transduction (ACS, JAR4, and AUX28), and fruit flavor and aroma (SUS2, LOX, and GDH2). In addition, proteins differentially abundant were associated with fruit texture (PLY, PME, and PGA) and plant hormone signal transduction (PBL15, JAR1, and PYL9). Moreover, NO significantly increased the abundance of key enzymes involved in the phenylpropanoid biosynthetic pathway, thus enhancing the disease resistance and alleviating softening of bitter melon fruit. Finally, differential metabolites mainly included phenolic acids, terpenoids, and flavonoids. These results provide a theoretical basis for further studies on the physiological changes associated with postharvest ripening and senescence of bitter melon fruit.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-023-00110-y.

Trancriptome data mining in combination with co-expression network analysis identifies the functional modules and critical regulators in Hordeum vulgare L. in response to cold stress

Cold stress, as an abiotic stress, is one of the most limiting factors which pose a great threat to the plant's productivity. To understand the transcriptional regulation and connectivity pattern of genes involved in barley cold stress responses, co-expression network analysis was performed based on the global transcriptome profiling. The microarray datasets related to cold stress treatments were retrieved from the Gene Expression Omnibus (GEO) and Array express databases. Four microarray datasets related to cold stress-responsive transcriptome in barley were included in our study. Gene co-expression analysis was constructed using WGCNA method. Module-Trait Relationships (MTR) analysis and hub genes determination and validation were carried out. Finally, transcription factor and kinase regulatory networks were Inferred using machine learning algorithm. The co-expression modules were determined using beta index = 10. In total 13 co-expressed modules were identified with an average size of 153 genes. Functional enrichment based on gene ontology (GO) showed that each of the stress related significant modules were enriched in different biological processes. Annotation of significant modules identifies some TFs and Kinases such as ethylene-responsive transcription factor 1-like, transcription factor PCL1-like, transcription factor MYC2, WRKY, serine/threonine-protein kinase PBL7, and receptor-like protein kinase At2g42960 were contributed in barley cold stress response. Our analysis highlighted the functional importance of ABA signaling pathway, ROS signaling, defensive and protective proteins, degrading protein, Ca2+ related signaling, ribosome-mediated translation and etc. in responding of barley to cold stress condition. The current findings add substantially to our understanding of the cold responsive underlying mechanism of barley which can serve in future studies and breeding programs.

Genome-Wide Analysis of the Liriodendron chinense Hsf Gene Family under Abiotic Stress and Characterization of the LcHsfA2a Gene

Heat shock factors (Hsfs) play a crucial role in plant defense processes. However, the distribution and functional characteristics of Hsf genes in the relict plant Liriodendron chinense are still unclear. In this study, a total of 19 LcHsfs were identified and divided into three separate subgroups, comprising 10 LcHsfA, 7 LcHsfB, and 2 LcHsfC genes, respectively, based on their phylogenetic tree and the presence/absence of conserved protein domains. Whole-genome duplication and segmental duplication led to an expansion of the LhHsf gene family. The promoters of LcHsf genes are enriched for different types of cis-acting elements, including hormone responsive and abiotic-stress-responsive elements. The expression of LcHsfA3, LcHsfA4b, LcHsfA5, LcHsfB1b, and LcHsfB2b increased significantly as a result of both cold and drought treatments. LcHsfA2a, LcHsfA2b, and LcHsfA7 act as important genes whose expression levels correlate strongly with the expression of the LcHsp70, LcHsp110, and LcAPX genes under heat stress. In addition, we found that transiently transformed 35S:LcHsfA2a seedlings showed significantly lower levels of hydrogen peroxide (H2O2) after heat stress and showed a stronger thermotolerance. This study sheds light on the possible functions of LcHsf genes under abiotic stress and identifies potentially useful genes to target for molecular breeding, in order to develop more stress-resistant varieties.

Transcriptomic and physiological analysis of the response of Spirodela polyrrhiza to sodium nitroprusside

BACKGROUND

Spirodela polyrrhiza is a simple floating aquatic plant with great potential in synthetic biology. Sodium nitroprusside (SNP) stimulates plant development and increases the biomass and flavonoid content in some plants. However, the molecular mechanism of SNP action is still unclear.

RESULTS

To determine the effect of SNP on growth and metabolic flux in S. polyrrhiza, the plants were treated with different concentrations of SNP. Our results showed an inhibition of growth, an increase in starch, soluble protein, and flavonoid contents, and enhanced antioxidant enzyme activity in plants after 0.025 mM SNP treatment. Differentially expressed transcripts were analysed in S. polyrrhiza after 0.025 mM SNP treatment. A total of 2776 differentially expressed genes (1425 upregulated and 1351 downregulated) were identified. The expression of some genes related to flavonoid biosynthesis and NO biosynthesis was upregulated, while the expression of some photosynthesis-related genes was downregulated. Moreover, SNP stress also significantly influenced the expression of transcription factors (TFs), such as ERF, BHLH, NAC, and WRKY TFs.

CONCLUSIONS

Taken together, these findings provide novel insights into the mechanisms of underlying the SNP stress response in S. polyrrhiza and show that the metabolic flux of fixed CO2 is redirected into the starch synthesis and flavonoid biosynthesis pathways after SNP treatment.

Functional and Transcriptome Analysis Reveal Specific Roles of Dimocarpus longan DlRan3A and DlRan3B in Root Hair Development, Reproductive Growth, and Stress Tolerance

Ran GTPases play essential roles in plant growth and development. Our previous studies revealed the nuclear localization of DlRan3A and DlRan3B proteins and proposed their functional redundancy and distinction in Dimocarpus longan somatic embryogenesis, hormone, and abiotic stress responses. To further explore the possible roles of DlRan3A and DlRan3B, gene expression analysis by qPCR showed that their transcripts were both more abundant in the early embryo and pulp in longan. Heterologous expression of DlRan3A driven by its own previously cloned promoter led to stunted growth, increased root hair density, abnormal fruits, bigger seeds, and enhanced abiotic stress tolerance. Conversely, constitutive promoter CaMV 35S (35S)-driven expression of DlRan3A, 35S, or DlRan3B promoter-controlled expression of DlRan3B did not induce the alterations in growth phenotype, while they rendered different hypersensitivities to abiotic stresses. Based on the transcriptome profiling of longan Ran overexpression in tobacco plants, we propose new mechanisms of the Ran-mediated regulation of genes associated with cell wall biosynthesis and expansion. Also, the transgenic plants expressing DlRan3A or DlRan3B genes controlled by 35S or by their own promoter all exhibited altered mRNA levels of stress-related and transcription factor genes. Moreover, DlRan3A overexpressors were more tolerant to salinity, osmotic, and heat stresses, accompanied by upregulation of oxidation-related genes, possibly involving the Ran-RBOH-CIPK network. Analysis of a subset of selected genes from the Ran transcriptome identified possible cold stress-related roles of brassinosteroid (BR)-responsive genes. The marked presence of genes related to cell wall biosynthesis and expansion, hormone, and defense responses highlighted their close regulatory association with Ran.

Protein interaction mapping reveals widespread targeting of development-related host transcription factors by phytoplasma effectors

Phytoplasmas are pathogenic bacteria that reprogram plant host development for their own benefit. Previous studies have characterized a few different phytoplasma effector proteins that destabilize specific plant transcription factors. However, these are only a small fraction of the potential effectors used by phytoplasmas; therefore, the molecular mechanisms through which phytoplasmas modulate their hosts require further investigation. To obtain further insights into the phytoplasma infection mechanisms, we generated a protein-protein interaction network between a broad set of phytoplasma effectors and a large, unbiased collection of Arabidopsis thaliana transcription factors and transcriptional regulators. We found widespread, but specific, interactions between phytoplasma effectors and host transcription factors, especially those related to host developmental processes. In particular, many unrelated effectors target specific sets of TCP transcription factors, which regulate plant development and immunity. Comparison with other host-pathogen protein interaction networks shows that phytoplasma effectors have unusual targets, indicating that phytoplasmas have evolved a unique and unusual infection strategy. This study contributes a rich and solid data source that guides further investigations of the functions of individual effectors, as demonstrated for some herein. Moreover, the dataset provides insights into the underlying molecular mechanisms of phytoplasma infection.

Transcriptome Expression Profiling Reveals the Molecular Response to Salt Stress in Gossypium anomalum Seedlings

Some wild cotton species are remarkably tolerant to salt stress, and hence represent valuable resources for improving salt tolerance of the domesticated allotetraploid species Gossypium hirsutum L. Here, we first detected salt-induced stress changes in physiological and biochemical indexes of G. anomalum, a wild African diploid cotton species. Under 350 mmol/L NaCl treatment, the photosynthetic parameters declined significantly, whereas hydrogen peroxide (H2O2) and malondialdehyde (MDA) contents increased. Catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD) activity and proline (PRO) content also significantly increased, reaching peak values at different stages of salt stress. We used RNA-Seq to characterize 15,476 differentially expressed genes in G. anomalum roots after 6, 12, 24, 72, and 144 h of salt stress. Gene Ontology enrichment analysis revealed these genes to be related to sequence-specific DNA and iron ion binding and oxidoreductase, peroxidase, antioxidant, and transferase activity; meanwhile, the top enriched pathways from the Kyoto Encyclopedia of Genes and Genomes database were plant hormone signal transduction, phenylpropanoid biosynthesis, fatty acid degradation, carotenoid biosynthesis, zeatin biosynthesis, starch and sucrose metabolism, and MAPK signaling. A total of 1231 transcription factors were found to be expressed in response to salt stress, representing ERF, MYB, WRKY, NAC, C2H2, bZIP, and HD-ZIP families. Nine candidate genes were validated by quantitative real-time PCR and their expression patterns were found to be consistent with the RNA-Seq data. These data promise to significantly advance our understanding of the molecular response to salt stress in Gossypium spp., with potential value for breeding applications.

GhCNGC13 and 32 Act as Critical Links between Growth and Immunity in Cotton

Cyclic nucleotide-gated ion channels (CNGCs) remain poorly studied in crop plants, most of which are polyploid. In allotetraploid Upland cotton (Gossypium hirsutum), silencing GhCNGC13 and 32 impaired plant growth and shoot apical meristem (SAM) development, while triggering plant autoimmunity. Both growth hormones (indole-3-acetic acid and gibberellin) and stress hormones (abscisic acid, salicylic acid, and jasmonate) increased, while leaf photosynthesis decreased. The silenced plants exhibited an enhanced resistance to Botrytis cinerea; however, Verticillium wilt resistance was weakened, which was associated with LIPOXYGENASE2 (LOX2) downregulation. Transcriptomic analysis of silenced plants revealed 4835 differentially expressed genes (DEGs) with functional enrichment in immunity and photosynthesis. These DEGs included a set of transcription factors with significant over-representation in the HSF, NAC, and WRKY families. Moreover, numerous members of the GhCNGC family were identified among the DEGs, which may indicate a coordinated action. Collectively, our results suggested that GhCNGC13 and 32 functionally link to photosynthesis, plant growth, and plant immunity. We proposed that GhCNGC13 and 32 play a critical role in the "growth-defense tradeoff" widely observed in crops.

Role of bZIP Transcription Factors in Response to NaCl Stress in Tamarix ramosissima under Exogenous Potassium (K+)

Salt stress is a significant environmental factor affecting plant growth and development, with NaCl stress being one of the most common types of salt stress. The halophyte, Tamarix ramosissima Ledeb (T. ramosissima), is frequently utilized for the afforestation of saline-alkali soils. Indeed, there has been limited research and reports by experts and scholars on the regulatory mechanisms of basic leucine zipper (bZIP) genes in T. ramosissima when treated with exogenous potassium (K+) to alleviate the effects of NaCl stress. This study focused on the bZIP genes in T. ramosissima roots under NaCl stress with additional KCl applied. We identified key candidate genes and metabolic pathways related to bZIP and validated them through quantitative real-time PCR (qRT-PCR). The results revealed that under NaCl stress with additional KCl applied treatments at 0 h, 48 h, and 168 h, based on Pfam protein domain prediction and physicochemical property analysis, we identified 20 related bZIP genes. Notably, four bZIP genes (bZIP_2, bZIP_6, bZIP_16, and bZIP_18) were labeled with the plant hormone signal transduction pathway, showing a predominant up-regulation in expression levels. The results suggest that these genes may mediate multiple physiological pathways under NaCl stress with additional KCl applied at 48 h and 168 h, enhancing signal transduction, reducing the accumulation of ROS, and decreasing oxidative damage, thereby enhancing the tolerance of T. ramosissima to NaCl stress. This study provides gene resources and a theoretical basis for further breeding of salt-tolerant Tamarix species and the involvement of bZIP transcription factors in mitigating NaCl toxicity.

Identification of candidate genes associated with less-photosensitive anthocyanin phenotype using an EMS mutant (pind) in eggplant (Solanum melongena L.)

Eggplant (Solanum melongena L.) is a highly nutritious and economically important vegetable crop. However, the fruit peel of eggplant often shows poor coloration owing to low-light intensity during cultivation, especially in the winter. The less-photosensitive varieties produce anthocyanin in low light or even dark conditions, making them valuable breeding materials. Nevertheless, genes responsible for anthocyanin biosynthesis in less-photosensitive eggplant varieties are not characterized. In this study, an EMS mutant, named purple in the dark (pind), was used to identify the key genes responsible for less-photosensitive coloration. Under natural conditions, the peel color and anthocyanin content in pind fruits were similar to that of wildtype '14-345'. The bagged pind fruits were light purple, whereas those of '14-345' were white; and the anthocyanin content in the pind fruit peel was significantly higher than that in '14-345'. Genetic analysis revealed that the less-photosensitive trait was controlled by a single dominant gene. The candidate gene was mapped on chromosome 10 in the region 7.72 Mb to 11.71 Mb. Thirty-five differentially expressed genes, including 12 structural genes, such as CHS, CHI, F3H, DFR, ANS, and UFGT, and three transcription factors MYB113, GL3, and TTG2, were identified in pind using RNA-seq. Four candidate genes EGP21875 (myb domain protein 113), EGP21950 (unknown protein), EGP21953 (CAAX amino-terminal protease family protein), and EGP21961 (CAAX amino-terminal protease family protein) were identified as putative genes associated with less-photosensitive anthocyanin biosynthesis in pind. These findings may clarify the molecular mechanisms underlying less-photosensitive anthocyanin biosynthesis in eggplant.

A comparative full-length transcriptomic resource provides insight into the perennial monocarpic mass flowering

Perennial monocarpic mass flowering represents as a key developmental innovation in flowering time diversity in several biological and economical essential families, such as the woody bamboos and the shrubby Strobilanthes. However, molecular and genetic mechanisms underlying this important biodiversity remain poorly investigated. Here, we generated a full-length transcriptome resource incorporated into the BlueOmics database (http://blueomics.iflora.cn) for two Strobilanthes species, which feature contrasting flowering time behaviors. Using about 112 and 104 Gb Iso-seq reads together with ~185 and ~75 Gb strand-specific RNA seq data, we annotated 80 971 and 79 985 non-redundant full-length transcripts for the perennial polycarpic Strobilanthes tetrasperma and the perennial monocarpic Strobilanthes biocullata, respectively. In S. tetrasperma, we identified 8794 transcripts showing spatiotemporal expression in nine tissues. In leaves and shoot apical meristems at two developmental stages, 977 and 1121 transcripts were differentially accumulated in S. tetrasperma and S. biocullata, respectively. Interestingly, among the 33 transcription factors showing differential expression in S. tetrasperma but without differential expression in S. biocullata, three were involved potentially in the photoperiod and circadian-clock pathway of flowering time regulation (FAR1 RELATED SEQUENCE 12, FRS12; NUCLEAR FACTOR Y A1, NFYA1; PSEUDO-RESPONSE REGULATOR 5, PRR5), hence provides an important clue in deciphering the flowering diversity mechanisms. Our data serve as a key resource for further dissection of molecular and genetic mechanisms underpinning key biological innovations, here, the perennial monocarpic mass flowering.

Comparative Transcriptome Sequencing and Endogenous Phytohormone Content of Annual Grafted Branches of Zelkova schneideriana and Its Dwarf Variety HenTianGao

Zelkova schneideriana is a fast-growing tree species endemic to China. Recent surveys and reports have highlighted a continued decline in its natural populations; therefore, it is included in the Red List of Threatened Species by The International Union for Conservation of Nature. A new variety "HenTianGao" (H) has been developed with smaller plant height, slow growth, and lower branching points. In this study, we attempted to understand the differences in plant height of Z. schneideriana (J) and its dwarf variety H. We determined the endogenous hormone content in the annual grafted branches of both J and H. J exhibited higher gibberellic acid (GA)-19 and trans-Zeatin (tZ) levels, whereas H had higher levels of indole-3-acetic acid (IAA) catabolite 2-oxindole-3-acetic acid (OxIAA), IAA-Glu conjugate, and jasmonic acid (JA) (and its conjugate JA-Ile). The transcriptome comparison showed differential regulation of 20,944 genes enriched in growth and development, signaling, and metabolism-related pathways. The results show that the differential phytohormone level (IAA, JA, tZ, and GA) was consistent with the expression of the genes associated with their biosynthesis. The differences in relative OxIAA, IAA-Glu, GA19, trans-Zeatin, JA, and JA-Ile levels were linked to changes in respective signaling-related genes. We also observed significant differences in the expression of cell size, number, proliferation, cell wall biosynthesis, and remodeling-related genes in J and H. The differences in relative endogenous hormone levels, expression of biosynthesis, and signaling genes provide a theoretical basis for understanding the plant height differences in Z. schneideriana.

Comparative transcriptome analysis of gene responses of salt-tolerant and salt-sensitive rice cultivars to salt stress

Salt stress is one unfavorable factor of global climate change that adversely affects rice plant growth and yield. To identify novel salt-tolerant genes and new varieties of salt-tolerant rice, a better understanding of the molecular regulation mechanism of salt tolerance in rice is needed. In this study we used transcriptome analyses to examine changes in gene expression of salt-tolerant and salt-sensitive rice plants. The salt-tolerant cultivar HH11 and salt-sensitive cultivar IR29 were treated with 200 mM NaCl solution for 0 h, 6 h, 24 h and 48 h at the three leaf stage. Physiological parameters and transcriptome were measured and analyzed after each treatment. Activity of SOD and POD, as well as the MDA and protein content of the two rice cultivars generally increased with increasing time of exposure to NaCl. Meanwhile, the APX activity first increased, then decreased in both cultivars, with maximum values seen at 6 h for IR29 and at 24 h for HH11. The GR and GPX activity of HH11 were stronger than that of IR29 in response to salt stress. The H2O2 content first increased at 0-6 h, then decreased at 6-24 h, and then increased again at 24-48 h under salt stress. Compared with IR29, SOD, POD and APX activity of HH11 was more sluggish in response to salt stress, reaching the maximum at 24 h or 48 h. The MDA, H2O2 and proline content of HH11 was lower than that of IR29 under salt stress. Relative to untreated HH11 plants (0 h) and those exposed to salt for 6 h, 24 h, and 48 h (H0-H6, H0-H24 and H0-H48), 7462, 6363 and 6636, differentially expressed genes (DEGs), respectively, were identified. For IR29, the respective total DEGs were 7566, 6075 and 6136. GO and KEGG enrichment analysis showed that metabolic pathways related to antioxidative responses and osmotic balance played vital roles in salt stress tolerance. Sucrose and starch metabolism, in addition to flavonoid biosynthesis and glutathione metabolism, showed positive responses to salt stress. Expression of two SPS genes (LOC_Os01g69030 and LOC_Os08g20660) and two GST genes (LOC_Os06g12290 and LOC_Os10g38740) was up-regulated in both HH11 and IR29, whereas expression of LOC_Os09g12660, a glucose-1-phosphate adenylyltransferase gene, and two SS genes (LOC_Os04g17650 and LOC_Os04g24430) was up-regulated differential expression in HH11. The results showed that HH11 had more favorable adjustment in antioxidant and osmotic activity than IR29 upon exposure to salt stress, and highlighted candidate genes that could play roles in the function and regulation mechanism of salt tolerance in rice.

Transcriptional deciphering of the metabolic pathways associated with the bioactive ingredients of wolfberry species with different quality characteristics

BACKGROUND

Wolfberry is rich in carotenoids, flavonoids, vitamins, alkaloids, betaines and other bioactive ingredients. For over 2,000 years, wolfberry has been used in China as a medicinal and edible plant resource. Nevertheless, the content of bioactive ingredients varies by cultivars, resulting in uneven quality across wolfberry cultivars and species. To date, research has revealed little about the underlying molecular mechanism of the metabolism of flavonoids, carotenoids, and other bioactive ingredients in wolfberry.

RESULTS

In this context, the transcriptomes of the Lycium barbarum L. cultivar 'Ningqi No. 1' and Lycium chinense Miller were compared during the fruit maturity stage using the Illumina NovaSeq 6000 sequencing platform, and subsequently, the changes of the gene expression profiles in two types of wolfberries were analysed. In total, 256,228,924 clean reads were obtained, and 8817 differentially expressed genes (DEGs) were identified, then assembled by Basic Local Alignment Search Tool (BLAST) similarity searches and annotated using Gene Ontology (GO), Clusters of Orthologous Groups of proteins (KOG), and the Kyoto Encyclopedia of Genes and Genomes (KEGG). By combining these transcriptome data with data from the PubMed database, 36 DEGs related to the metabolism of bioactive ingredients and implicated in the metabolic pathway of carotenoids, flavonoids, terpenoids, alkaloids, vitamins, etc., were identified. In addition, among the 9 differentially expressed transcription factors, LbAPL, LbPHL11 and LbKAN4 have raised concerns. The protein physicochemical properties, structure prediction and phylogenetic analysis indicated that LbAPL and LbPHL11 may be good candidate genes involved in regulating the flavonoid metabolism pathway in wolfberry.

CONCLUSIONS

This study provides preliminary evidence for the differences in bioactive ingredient content at the transcription level among different wolfberry species, as well as a research and theoretical basis for the screening, cloning and functional analysis of key genes involved in the metabolism of bioactive ingredients in wolfberry.

De novo sequencing, assembly, and characterization of Asparagus racemosus transcriptome and analysis of expression profile of genes involved in the flavonoid biosynthesis pathway

Asparagus racemosus is known for its diverse content of secondary metabolites, i.e., saponins, alkaloids, and a wide range of flavonoids. Flavonoids, including phenols and polyphenols, have a significant role in plant physiology and are synthesized in several tissues. Despite the diverse role of flavonoids, genetic information is limited for flavonoid biosynthesis pathways in A. racemosus. The current study explores full-scale functional genomics information of A. racemosus by de novo transcriptome sequencing using Illumina paired-end sequencing technology to elucidate the genes involved in flavonoid biosynthesis pathways. The de novo assembly of high-quality paired-end reads resulted in ∼2.3 million high-quality reads with a pooled transcript of 45,647 comprising ∼76 Mb transcriptome with a mean length (bp) of 1,674 and N50 of 1,868bp. Furthermore, the coding sequence (CDS) prediction analysis from 45,647 pooled transcripts resulted in 45,444 CDS with a total length and mean length of 76,398,686 and 1,674, respectively. The Gene Ontology (GO) analysis resulted in a high number of CDSs assigned to 25,342 GO terms, which grouped the predicted CDS into three main domains, i.e., Biological Process (19,550), Molecular Function (19,873), and Cellular Component (14,577). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database was used to categorize 6,353 CDS into 25 distinct biological pathway categories, in which the majority of mapped CDS were shown to be related to translation (645), followed by signal transduction (532), carbohydrate metabolism (524), folding, sorting, and degradation (522). Among these, only ∼64 and 14 CDSs were found to be involved in the phenylpropanoid and flavonoid biosynthesis pathways, respectively. Quantitative Real-time PCR was used to check the expression profile of fourteen potential flavonoid biosynthesis pathway genes. The qRT-PCR analysis result matches the transcriptome sequence data validating the Illumina sequence results. Moreover, a large number of genes associated with the flavonoids biosynthesis pathway were found to be upregulated under the induction of methyl jasmonate. The present-day study on transcriptome sequence data of A. racemosus can be utilized for characterizing genes involved in flavonoid biosynthesis pathways and for functional genomics analysis in A. racemosus using the reverse genetics approach (CRISPR/Cas9 technology).

SlWRKY16 and SlWRKY31 of tomato, negative regulators of plant defense, involved in susceptibility activation following root-knot nematode Meloidogyne javanica infection

The involvement of WRKY transcription factors in plant-nematode interactions, and in particular, how these WRKYs participate in regulating the complex morphological and physiological changes occurring after nematode infection, are the topic of active research. We characterized the functional role of the unstudied tomato WRKY genes SlWRKY16 and SlWRKY31 in regulating tomato roots' response to infection by the root-knot nematode Meloidogyne javanica. Using promoter-GUS reporter gene fusions and qRT-PCR, we show that both SlWRKYs are predominantly expressed during the first half of the parasitic life stages, when feeding-site induction and construction occur. Expression of SlWRKY16 increased sharply 15 days after inoculation, whereas SlWRKY31 was already induced earlier, but reached its maximum expression at this time. Both genes were downregulated at the mature female stage. To determine biological function, we produced transgenic lines overexpressing SlWRKY16 and SlWRKY31 in tomato hairy roots. Overexpression of both genes resulted in enhanced M. javanica infection, reflected by increased galling occurrence and reproduction. Expression profiling of marker genes responsive to defense-associated phytohormones indicated reductions in salicylic acid defense-related PR-1 and jasmonic acid defense-related PI in inoculated roots overexpressing SlWRK16 and SlWRKY31, respectively. Our results suggest that SlWRKY16 and SlWRKY31 function as negative regulators of plant immunity induced upon nematode infection.

Integrated Metabolomic and Transcriptomic Analysis of the Quinoa Seedling Response to High Relative Humidity Stress

Quinoa is of great interest because it is cold- and drought-resistant; however, little research has been performed on quinoa under high relative humidity (RH) stress. In this study, quinoa seedlings of a highly HR-resistant variety ("Dianli-439") and a sensitive variety ("Dianli-969") were subjected to morphological and physiological measurements and metabolome and transcriptome analyses to investigate their response to high RH stress. In total, 1060 metabolites were detected, and lipids and flavonoids were the most abundant, with 173 and 167 metabolites, respectively. In total, 13,095 differentially expressed genes were identified, and the results showed that abscisic acid, auxin, and jasmonic-acid-related genes involved in plant hormone signaling may be involved in the response of quinoa seedlings to high RH stress. The analysis of the transcription factors revealed that the AP2/ERF family may also play an important role in the response to high RH stress. We identified the possible regulatory mechanisms of the hormone signaling pathways under high RH stress. Our findings can provide a basis for the selection and identification of highly resistant quinoa varieties and the screening of the metabolite-synthesis- and gene-regulation-related mechanisms in quinoa in response to RH stress.

Genome-Wide Identification and Expression Analysis of C3H Zinc Finger Family in Potato (Solanum tuberosum L.)

Transcription factors containing a CCCH structure (C3H) play important roles in plant growth and development, and their stress response, but research on the C3H gene family in potato has not been reported yet. In this study, we used bioinformatics to identify 50 C3H genes in potato and named them StC3H-1 to StC3H-50 according to their location on chromosomes, and we analyzed their physical and chemical properties, chromosome location, phylogenetic relationship, gene structure, collinearity relationship, and cis-regulatory element. The gene expression pattern analysis showed that many StC3H genes are involved in potato growth and development, and their response to diverse environmental stresses. Furthermore, RT-qPCR data showed that the expression of many StC3H genes was induced by high temperatures, indicating that StC3H genes may play important roles in potato response to heat stress. In addition, Some StC3H genes were predominantly expressed in the stolon and developing tubers, suggesting that these StC3H genes may be involved in the regulation of tuber development. Together, these results provide new information on StC3H genes and will be helpful for further revealing the function of StC3H genes in the heat stress response and tuber development in potato.

Metabolic characteristics of self-pollinated wheat seed under red and blue light during early development

MAIN CONCLUSION

Blue light has a greater effect on jasmonic acid and flavonoid accumulation in wheat seeds than red light; blue light reduces starch synthesis and the size of starch granules and seeds. This study sought to elucidate the effects of blue and red light on seed metabolism to provide important insights regarding the role of light quality in regulating seed growth and development. We used combined multi-omics analysis to investigate the impact of red and blue light (BL) on the induction of secondary metabolite accumulation in the hexaploid wheat Dianmai 3 after pollination. Flavonoids and alkaloids were the most differentially abundant metabolites detected under different treatments. Additionally, we used multi-omics and weighted correlation network analysis to screen multiple candidate genes associated with jasmonic acid (JA) and flavonoids. Expression regulatory networks were constructed based on RNA-sequencing data and their potential binding sites. The results revealed that BL had a greater effect on JA and flavonoid accumulation in wheat seeds than red light. Furthermore, BL reduced starch synthesis and stunted the size of starch granules and seeds. Collectively, these findings clarify the role of BL in the metabolic regulation of early seed development in wheat.

Recent Advances in Molecular Mechanism and Breeding Utilization of Brown Planthopper Resistance Genes in Rice: An Integrated Review

Over half of the world's population relies on rice as their staple food. The brown planthopper (Nilaparvata lugens Stål, BPH) is a significant insect pest that leads to global reductions in rice yields. Breeding rice varieties that are resistant to BPH has been acknowledged as the most cost-effective and efficient strategy to mitigate BPH infestation. Consequently, the exploration of BPH-resistant genes in rice and the development of resistant rice varieties have become focal points of interest and research for breeders. In this review, we summarized the latest advancements in the localization, cloning, molecular mechanisms, and breeding of BPH-resistant rice. Currently, a total of 70 BPH-resistant gene loci have been identified in rice, 64 out of 70 genes/QTLs were mapped on chromosomes 1, 2, 3, 4, 6, 8, 10, 11, and 12, respectively, with 17 of them successfully cloned. These genes primarily encode five types of proteins: lectin receptor kinase (LecRK), coiled-coil-nucleotide-binding-leucine-rich repeat (CC-NB-LRR), B3-DNA binding domain, leucine-rich repeat domain (LRD), and short consensus repeat (SCR). Through mediating plant hormone signaling, calcium ion signaling, protein kinase cascade activation of cell proliferation, transcription factors, and miRNA signaling pathways, these genes induce the deposition of callose and cell wall thickening in rice tissues, ultimately leading to the inhibition of BPH feeding and the formation of resistance mechanisms against BPH damage. Furthermore, we discussed the applications of these resistance genes in the genetic improvement and breeding of rice. Functional studies of these insect-resistant genes and the elucidation of their network mechanisms establish a strong theoretical foundation for investigating the interaction between rice and BPH. Furthermore, they provide ample genetic resources and technical support for achieving sustainable BPH control and developing innovative insect resistance strategies.

PearMODB: a multiomics database for pear (Pyrus) genomics, genetics and breeding study

Pear (Pyrus ssp.) belongs to Rosaceae and is an important fruit tree widely cultivated around the world. Currently, challenges to cope with the burgeoning sets of multiomics data are rapidly increasing. Here, we constructed the Pear Multiomics Database (PearMODB) by integrating genome, transcriptome, epigenome and population variation data, and aimed to provide a portal for accessing and analyzing pear multiomics data. A variety of online tools were built including gene search, BLAST, JBrowse, expression heatmap, synteny analysis and primer design. The information of DNA methylation sites and single-nucleotide polymorphisms can be retrieved through the custom JBrowse, providing an opportunity to explore the genetic polymorphisms linked to phenotype variation. Moreover, different gene families involving transcription factors, transcription regulators and disease resistance (nucleotide-binding site leucine-rich repeat) were identified and compiled for quick search. In particular, biosynthetic gene clusters (BGCs) were identified in pear genomes, and specialized webpages were set up to show detailed information of BGCs, laying a foundation for studying metabolic diversity among different pear varieties. Overall, PearMODB provides an important platform for pear genomics, genetics and breeding studies. Database URL http://pearomics.njau.edu.cn.

Transcriptomic profile of lettuce seedlings (Lactuca sativa) response to microalgae extracts used as biostimulant agents

To reduce the use of chemical fertilizers and maximize agricultural yields, the use of microalgae extracts as biostimulants has recently attracted significant attention due to their favourable impact on both plant growth and their ability to induce tolerance towards environmental stressors. Lettuce (Lactuca sativa) is one of the most important fresh vegetables that often requires applications of chemical fertilizers to increase quality and productivity. Accordingly, the purpose of this study was to analyse the transcriptome reprogramming of lettuce (L. sativa) seedlings in response to either Chlorella vulgaris or Scenedesmus quadricauda extracts by applying an RNAseq approach. Differential gene expression analysis revealed that the core gene set that responded to microalgal treatments in a species-independent manner includes 1330 clusters, 1184 of which were down-regulated and 146 up-regulated, clearly suggesting that the repression of gene expression is the main effect of algal treatments. The deregulation of 7197 transcripts in the C. vulgaris treated seedlings compared to control samples (LsCv vs. LsCK) and 7118 transcripts in the S. quadricauda treated seedlings compared to control samples (LsSq vs. LsCK) were counted. Although the number of deregulated genes turned out to be similar between the algal treatments, the level of deregulation was higher in LsCv versus LsCK than in LsSq versus LsCK. In addition, 2439 deregulated transcripts were observed in the C. vulgaris treated seedlings compared to S. quadricauda treated samples (LsCv vs. LsSq comparison) suggesting that a specific transcriptomic profile was induced by the single algal extracts. 'Plant hormone signal transduction' category includes a very elevated number of DEGs, many of them specifically indicating that C. vulgaris actives both genes involved in the auxin biosynthesis and transduction pathways, whereas S. quadricauda up-regulates genes implicated in the cytokinin biosynthesis pathway. Finally, algal treatments induced the deregulation of genes encoding small hormone-like molecules that are known to act alone or by interacting with major plant hormones. In conclusion, this study offers the groundwork to draw up a list of putative gene targets with the aim of lettuce genetic improvement that will allow a limited or even null use of synthetic fertilizers and pesticides in the management of this crop.

Widening the landscape of transcriptional regulation of green algal photoprotection

Availability of light and CO₂, substrates of microalgae photosynthesis, is frequently far from optimal. Microalgae activate photoprotection under strong light, to prevent oxidative damage, and the CO₂ Concentrating Mechanism (CCM) under low CO₂, to raise intracellular CO₂ levels. The two processes are interconnected; yet, the underlying transcriptional regulators remain largely unknown. Employing a large transcriptomic data compendium of Chlamydomonas reinhardtii's responses to different light and carbon supply, we reconstruct a consensus genome-scale gene regulatory network from complementary inference approaches and use it to elucidate transcriptional regulators of photoprotection. We show that the CCM regulator LCR1 also controls photoprotection, and that QER7, a Squamosa Binding Protein, suppresses photoprotection- and CCM-gene expression under the control of the blue light photoreceptor Phototropin. By demonstrating the existence of regulatory hubs that channel light- and CO₂-mediated signals into a common response, our study provides an accessible resource to dissect gene expression regulation in this microalga.

Comparative Transcriptome and Widely Targeted Metabolome Analysis Reveals the Molecular Mechanism of Powdery Mildew Resistance in Tomato

Powdery mildew is a serious problem in tomato production; therefore, the PM-resistant tomato inbred line, '63187', and the susceptible tomato variety, 'Moneymaker (MM)', were used as experimental materials for the combined analysis of transcriptome and widely targeted metabolome on tomato leaves at 0 h post inoculation (hpi), 12 hpi, and 48 hpi. The results indicated that 276 genes were expressed in all treatments, and the K-means cluster analysis showed that these genes were divided into eight classes in '63187' and ten classes in 'MM'. KEGG enrichment showed that amino acid metabolism, signal transduction, energy metabolism, and other secondary metabolites biosynthesis pathways were significantly enriched. Interestingly, the analysis of WRKY family transcription factors (TFs) showed that the expression of four TFs in '63187' increased with no obvious change in 'MM'; and the expression of one TF in 'MM' increased with no obvious change in '63187'. The combined analysis revealed that both phenylpropanoid biosynthesis and flavonoid biosynthesis pathways were enriched in '63187' and 'MM'. In '63187', six metabolites involved in this pathway were downregulated, and four genes were highly expressed, while in 'MM', three metabolites were upregulated, four metabolites were downregulated, and ten genes were highly expressed. These metabolites and genes might be candidates for PM resistance or susceptibility in subsequent studies. These results provide favorable molecular information for the study of the different resistances of tomatoes to PM, and they provide a basis for the breeding of tomato varieties resistant to PM.

Reciprocal inhibition of expression between RAV1 and BES1 modulates plant growth and development in Arabidopsis

RAV1 (Related to ABI3/VP1) is a plant-specific B3 and AP2 domain-containing transcription factor that acts as a negative regulator of growth in many plant species. The expression of RAV1 is downregulated by brassinosteroids (BRs); large-scale transcriptome analyses have shown that the expression of RAV1 was previously targeted by BRI1-EMS-SUPPRESOR1 (BES1) and BRASSINAZOLE-RESISTANT1 (BZR1), which are critical transcription factors for the BR-signaling process. Using RAV1-overexpressing transgenic plants, we showed that RAV1 overexpression reduced the BR signaling capacity, resulting in the downregulation of BR biosynthetic genes and BES1 expression. Furthermore, we demonstrated that BES1, not BZR1, is directly bound to the RAV1 promoter and repressed RAV1 expression, and vice versa; RAV1 is also bound to the BES1 promoter and repressed BES1 expression. This mutual inhibition was specific to RAV1 and BES1 because RAV1 exhibited binding activity to the BZR1 promoter but did not repress BZR1 expression. We observed that constitutively activated BR signaling phenotypes in bes1-D were attenuated by the repression of endogenous BES1 expression in transgenic bes1-D plants overexpressing RAV1. RNA-sequencing analysis of RAV1-overexpressing transgenic plants and bes1-D mutant plants revealed differentially expressed genes by RAV1 and BES1 and genes that were oppositely co-regulated by RAV1 and BES1. RAV1 and BES1 regulated different transcriptomes but co-regulated a specific set of genes responsible for the balance between growth and defense. These results suggested that the mutual inhibitory transcriptional activities of RAV1 and BES1 provide fine regulatory mechanisms for plant growth and development.

Regulation of the regulators: Transcription factors controlling biosynthesis of plant secondary metabolites during biotic stresses and their regulation by miRNAs

Biotic stresses threaten to destabilize global food security and cause major losses to crop yield worldwide. In response to pest and pathogen attacks, plants trigger many adaptive cellular, morphological, physiological, and metabolic changes. One of the crucial stress-induced adaptive responses is the synthesis and accumulation of plant secondary metabolites (PSMs). PSMs mitigate the adverse effects of stress by maintaining the normal physiological and metabolic functioning of the plants, thereby providing stress tolerance. This differential production of PSMs is tightly orchestrated by master regulatory elements, Transcription factors (TFs) express differentially or undergo transcriptional and translational modifications during stress conditions and influence the production of PSMs. Amongst others, microRNAs, a class of small, non-coding RNA molecules that regulate gene expression post-transcriptionally, also play a vital role in controlling the expression of many such TFs. The present review summarizes the role of stress-inducible TFs in synthesizing and accumulating secondary metabolites and also highlights how miRNAs fine-tune the differential expression of various stress-responsive transcription factors during biotic stress.

Overexpression of the cotton trihelix transcription factor GhGT23 in Arabidopsis mediates salt and drought stress tolerance by binding to GT and MYB promoter elements in stress-related genes

Cotton (Gossypium hirsutum L.) is the world's most economically valuable textile crop. However, cotton plants are often subjected to numerous abiotic stresses that can dramatically limit yield. Trihelix transcription factors (TTFs) play important roles in abiotic stress responses in many plant species, and efforts to better understand their roles in cotton abiotic stress responses are ongoing. In this study, a member of the cotton TTF family (GhGT23) was functionally characterized. This protein contains a SANT domain and is a member of the SIP subfamily of TTF proteins. GhGT23 was significantly (p < 0.05) and highly expressed in cotton fiber compared to relatively low expression in other tissues. A significant (p < 0.05) increase in GhGT23 expression occurred in cotton seedlings within 12 hours of drought, salt, and ABA exposure. The GhGT23 protein localized in the nucleus but exhibited no signs of transactivation activity. GhGT23 overexpression in Arabidopsis conferred enhanced drought and salt stress tolerance. The expression of stress-related genes was higher in transgenic Arabidopsis expressing GhGT23 than in wild-type plants subjected to salt stress. The results of electrophoretic mobility shift assay revealed that GhGT23 could bind to the GT cis-elements GT-1Box (Box II), GT2-Box, GT3-Box, GT-3a (Site1-type), GT-3b, and Box as well as the MYB cis-elements MBS1 and MRE4. Our results demonstrate that GhGT23 positively regulates salt and drought stress responses, possibly by enhancing the expression of stress-related genes.

Transcriptome dynamics provides insights into divergences of the photosynthesis pathway between Saccharum officinarum and Saccharum spontaneum

Saccharum spontaneum and Saccharum officinarum contributed to the genetic background of modern sugarcane cultivars. Saccharum spontaneum has shown a higher net photosynthetic rate and lower soluble sugar than S. officinarum. Here, we analyzed 198 RNA-sequencing samples to investigate the molecular mechanisms for the divergences of photosynthesis and sugar accumulation between the two Saccharum species. We constructed gene co-expression networks based on differentially expressed genes (DEGs) both for leaf developmental gradients and diurnal rhythm. Our results suggested that the divergence of sugar accumulation may be attributed to the enrichment of major carbohydrate metabolism and the oxidative pentose phosphate pathway. Compared with S. officinarum, S. spontaneum DEGs showed a high enrichment of photosynthesis and contained more complex regulation of photosynthesis-related genes. Noticeably, S. spontaneum lacked gene interactions with sulfur assimilation stimulated by photorespiration. In S. spontaneum, core genes related to clock and photorespiration displayed a sensitive regulation by the diurnal rhythm and phase-shift. Small subunit of Rubisco (RBCS) displayed higher expression in the source tissues of S. spontaneum. Additionally, it was more sensitive under a diurnal rhythm, and had more complex gene networks than that in S. officinarum. This indicates that the differential regulation of RBCS Rubisco contributed to photosynthesis capacity divergence in both Saccharum species.

Global Responses of Autopolyploid Sugarcane Badila (Saccharum officinarum L.) to Drought Stress Based on Comparative Transcriptome and Metabolome Profiling

Sugarcane (Saccharum spp. hybrid) is frequently affected by seasonal drought, which causes substantial declines in quality and yield. To understand the drought resistance mechanisms of S. officinarum, the main species of modern sugarcane, at a molecular level, we carried out a comparative analysis of transcriptome and metabolome profiling of the sugarcane variety Badila under drought stress (DS). Compared with control group (CG) plants, plants exposed to DS had 13,744 (6663 up-regulated and 7081 down-regulated) differentially expressed genes (DEGs). GO and KEGG analysis showed that the DEGs were enriched in photosynthesis-related pathways and most DEGs had down-regulated expression. Moreover, the chlorophyll content, photosynthesis (Photo), stomatal conductance (Cond), intercellular carbon dioxide concentration (Ci) and transpiration rate (Trmmol) were sharply decreased under DS. These results indicate that DS has a significant negative influence on photosynthesis in sugarcane. Metabolome analysis identified 166 (37 down-regulated and 129 up-regulated) significantly regulated metabolites (SRMs). Over 50% of SRMs were alkaloids, amino acids and their derivatives, and lipids. The five most significantly enriched KEGG pathways among SRMs were Aminoacyl-tRNA biosynthesis, 2-Oxocarboxylic acid metabolism, Biosynthesis of amino acids, Phenylalanine metabolism, and Arginine and proline metabolism (p < 0.05). Comparing CG with DS for transcriptome and metabolome profiling (T_CG/DS and M_CG/DS, respectively), we found three of the same KEGG-enriched pathways, namely Biosynthesis of amino acids, Phenylalanine metabolism and Arginine and proline metabolism. The potential importance of Phenylalanine metabolism and Arginine and proline metabolism was further analyzed for response to DS in sugarcane. Seven SRMs (five up-regulated and two down-regulated) and 60 DEGs (17 up-regulated and 43 down-regulated) were enriched in Phenylalanine metabolism under DS, of which novel.31261, Sspon.04G0008060-1A, Sspon.04G0008060-2B and Sspon.04G0008060-3C were significantly correlated with 7 SRMs. In Arginine and proline metabolism, eight SRMs (seven up-regulated and one down-regulated) and 63 DEGs (32 up-regulated and 31 down-regulated) were enriched, of which Sspon.01G0026110-1A (OAT) and Sspon.03G0002750-3D (P5CS) were strongly associated with proline (r > 0.99). These findings present the dynamic changes and possible molecular mechanisms of Phenylalanine metabolism as well as Arginine and proline metabolism under DS and provide a foundation for future research and sugarcane improvement.

Transcriptome analysis reveals the fruit color variation in Ailanthus altissima

Anthocyanins are responsible for the intensity of color in plants; however, the systematic mechanisms underlying the color differences in the fruit of Ailanthus altissima remain unknown. Therefore, this study aims to analyze the transcriptomes of the white and red fruit of A. altissima by screening and validating the key genes involved in flavonoid and anthocyanin biosynthesis. Samples of A. altissima fruit were collected 30, 45, and 60 days after flowering, and their pigment and sugar content were determined. The anthocyanin content was significantly higher in red than in white fruits. Transcriptome analysis was also performed on the fruit samples, 73,807 unigenes were assembled and annotated to seven databases. Twenty-one co-expressed modules were identified via weighted gene co-expression network analysis, of which two were associated with flavonoids and anthocyanins. Furthermore, in three growth stages, 126, 30, and 124 differentially expressed genes were screened between white and red fruit. Genes involved in flavonoid and anthocyanin metabolism were identified. AaDFR (A. altissima bifunctional dihydroflavonol 4-reductase/flavanone 4-reductase) and AaANS (A. altissima anthocyanidin synthase) were associated with flavonoid and anthocyanin metabolism. Members of the AaDFR and AaANS families were also identified, and their basic physicochemical characteristics, conserved domains, motif compositions, phylogenetics, and expression levels were analyzed. The overexpression of AaDFR and AaANS in transgenic Arabidopsis significantly increased the content of seed and foliar flavonoids and anthocyanins. The study elucidated the different mechanisms underlying fruit color development and provided insight into A. altissima plants breeding with commercially desirable properties.

Characterization and expression analysis of basic leucine zipper (bZIP) transcription factors responsive to chilling injury in peach fruit

BACKGROUND

Peach (Prunus persica L.) is prone to chilling injury as exhibited by inhibition of the ethylene production, failure in softening, and the manifestation of internal browning. The basic leucine zipper (bZIP) transcription factors play an essential role in regulatory networks that control many processes associated with physiological, abiotic and biotic stress responses in fruits. Formerly, the underlying molecular and regulatory mechanism of (bZIP) transcription factors responsive to chilling injury in peach fruit is still elusive.

METHODS AND RESULTS

In the current experiment, the solute peach 'Zhongyou Peach No. 13' was used as the test material and cold storage at low temperature (4 °C). It was found that long-term low-temperature storage induced the production of ethylene, the hardness of the pulp decreased, and the low temperature also induced ABA accumulation. The changes of ABA and ethylene in peach fruits during low-temperature storage were clarified. Since the bZIP transcription factor is involved in the regulation of downstream pathways of ABA signals, 47 peach bZIP transcription factor family genes were identified through bioinformatics analysis. Further based on RT-qPCR analysis, 18 PpbZIP genes were discovered to be expressed in refrigerated peach fruits. Among them, the expression of PpbZIP23 and PpbZIP25 was significantly reduced during the refrigeration process, the promoter analysis of these genes found that this region contains the MYC/MYB/ABRES binding element, but not the DRES/CBFS element, indicating that the expression may be regulated by the ABA-dependent cold induction pathway, thereby responding to chilling injury in peach fruit.

CONCLUSIONS

Over investigation will provide new insights for further postharvest protocols related to molecular changes during cold storage and will prove a better cope for chilling injury.

Genome-Wide Characterization and Analysis of R2R3-MYB Genes Related to Fruit Ripening and Stress Response in Banana (Musa acuminata L. AAA Group, cv. 'Cavendish')

MYB is an important type of transcription factor in eukaryotes. It is widely involved in a variety of biological processes and plays a role in plant morphogenesis, growth and development, primary and secondary metabolite synthesis, and other life processes. In this study, bioinformatics methods were used to identify the R2R3-MYB transcription factor family members in the whole Musa acuminata (DH-Pahang) genome, one of the wild ancestors of banana. A total of 280 MaMYBs were obtained, and phylogenetic analysis indicated that these MaMYBs could be classified into 33 clades with MYBs from Arabidopsis thaliana. The amino acid sequences of the R2 and R3 Myb-DNA binding in all MaMYB protein sequences were quite conserved, especially Arg-12, Arg-13, Leu-23, and Leu-79. Distribution mapping results showed that 277 MaMYBs were localized on the 11 chromosomes in the Musa acuminata genome. The MaMYBs were distributed unevenly across the 11 chromosomes. More than 40.0% of the MaMYBs were located in collinear fragments, and segmental duplications likely played a key role in the expansion of the MaMYBs. Moreover, the expression profiles of MaMYBs in different fruit development and ripening stages and under various abiotic and biotic stresses were investigated using available RNA-sequencing data to obtain fruit development, ripening-specific, and stress-responsive candidate genes. Weighted gene co-expression network analysis (WGCNA) was used to analyze transcriptome data of banana from the above 11 samples. We found MaMYBs participating in important metabolic biosynthesis pathways in banana. Collectively, our results represent a comprehensive genome-wide study of the MaMYB gene family, which should be helpful in further detailed studies on MaMYBs functions related to fruit development, postharvest ripening, and the seedling response to stress in an important banana cultivar.

Orchid NAC Transcription Factors: A Focused Analysis of CUPULIFORMIS Genes

Plant transcription factors are involved in different developmental pathways. NAC transcription factors (No Apical Meristem, Arabidopsis thaliana Activating Factor, Cup-shaped Cotyledon) act in various processes, e.g., plant organ formation, response to stress, and defense mechanisms. In Antirrhinum majus, the NAC transcription factor CUPULIFORMIS (CUP) plays a role in determining organ boundaries and lip formation, and the CUP homologs of Arabidopsis and Petunia are involved in flower organ formation. Orchidaceae is one of the most species-rich families of angiosperms, known for its extraordinary diversification of flower morphology. We conducted a transcriptome and genome-wide analysis of orchid NACs, focusing on the No Apical Meristem (NAM) subfamily and CUP genes. To check whether the CUP homologs could be involved in the perianth formation of orchids, we performed an expression analysis on the flower organs of the orchid Phalaenopsis aphrodite at different developmental stages. The expression patterns of the CUP genes of P. aphrodite suggest their possible role in flower development and symmetry establishment. In addition, as observed in other species, the orchid CUP1 and CUP2 genes seem to be regulated by the microRNA, miR164. Our results represent a preliminary study of NAC transcription factors in orchids to understand the role of these genes during orchid flower formation.

Natural variation of Alfin-like family affects seed size and drought tolerance in rice

The Alfin-like (AL) family is a group of small plant-specific transcriptional factors involved in abiotic stresses in dicotyledon. In an early study, we found an AL gene in rice that was associated with grain yield under drought stress. However, little information is known about the AL family in rice. In this study, AL genes in the rice genome were identified, and the OsAL proteins were found to locate in the nucleus and have no transcriptional self-activation activity. The expression of the OsALs was regulated by different environmental stimulations and plant hormones. Association and domestication analysis revealed that natural variation of most OsALs was significantly associated with yield traits, drought resistance and divergence in grain size in indica and japonica rice varieties. Hap1 of OsAL7.1 and Hap7 of OsAL11 were favorable haplotypes of seed weight and germination under osmotic stress. Furthermore, osal7.1 and osal11 mutants have larger seeds and are more sensitive to abscisic acid and mannitol during germination stage. Overexpressing of OsAL7.1 and OsAL11 in rice weakened the tolerance to drought in the adult stage. Thus, our work provides informative knowledge for exploring and harnessing haplotype diversity of OsALs to improve yield stability under drought stress.

Comparative transcriptomic profiling of peach and nectarine cultivars reveals cultivar-specific responses to chilled postharvest storage

INTRODUCTION

Peach (Prunus persica (L.) Batsch,) and nectarine fruits (Prunus persica (L.) Batsch, var nectarine), are characterized by a rapid deterioration at room temperature. Therefore, cold storage is widely used to delay fruit post-harvest ripening and extend fruit commercial life. Physiological disorders, collectively known as chilling injury, can develop typically after 3 weeks of low-temperature storage and affect fruit quality.

METHODS

A comparative transcriptomic analysis was performed to identify regulatory pathways that develop before chilling injury symptoms are detectable using next generation sequencing on the fruits of two contrasting cultivars, one peach (Sagittaria) and one nectarine, (Big Top), over 14 days of postharvest cold storage.

RESULTS

There was a progressive increase in the number of differentially expressed genes between time points (DEGs) in both cultivars. More (1264) time point DEGs were identified in 'Big Top' compared to 'Sagittaria' (746 DEGs). Both cultivars showed a downregulation of pathways related to photosynthesis, and an upregulation of pathways related to amino sugars, nucleotide sugar metabolism and plant hormone signal transduction with ethylene pathways being most affected. Expression patterns of ethylene related genes (including biosynthesis, signaling and ERF transcription factors) correlated with genes involved in cell wall modification, membrane composition, pathogen and stress response, which are all involved later during storage in development of chilling injury.

DISCUSSION

Overall, the results show that common pathways are activated in the fruit of 'Big Top' nectarine and 'Sagittaria' peach in response to cold storage but include also differences that are cultivar-specific responses.

Dehydration response in Stylosanthes scabra: Transcriptional, biochemical, and physiological modulations

Stylosanthes scabra, popularly known as stylo, is native to the Brazilian Caatinga semiarid region and stands out as a drought-tolerant shrub forage crop. This work provides information about the plant response during the first 48 h of water deficit, followed by a rehydration treatment. Besides root transcriptomics data, 13 physiological or biochemical parameters were scrutinized. Additionally, RNA-Seq annotated transcripts not associated with the "Viridiplantae" clade were taxonomically categorized. It was found that S. scabra quickly perceives and recovers from the oscillations of the imposed water regime. Physiologically, mechanisms that minimize evapotranspiration or protect the photosynthetic apparatus stood out. Biochemically, it was found that the root tissue invests in synthesizing compounds that can act as osmolytes (proline and sugars), emphasizing the importance of osmoregulation to water deficit acclimation. Consistently, transcriptome and qPCR analyses showed that a set of enriched biological processes with upregulated (UR) transcripts were involved in protective functions against reactive oxygen species or encoding enzymes of important metabolic pathways, which might contribute to S. scabra response to water deficit. Additionally, several UR kinases and transcription factors were identified. Finally, in an innovative approach, some naturally occurring microbial groups (such as Schizosaccharomyces, Bradyrhizobium, etc.) were identified in the S. scabra roots. This study reveals insights into the physiological, biochemical, and molecular mechanisms underlying the S. scabra response to water deficit and provides candidate genes that may be useful in developing drought-tolerant crop varieties through biotechnological applications.

Detection and validation of cis-regulatory motifs in osmotic stress-inducible synthetic gene switches via computational and experimental approaches

Synthetic cis -regulatory modules can improve our understanding of gene regulatory networks. We applied an ensemble approach for de novo cis motif discovery among the promoters of 181 drought inducible differentially expressed soybean (Glycine max L.) genes. A total of 43 cis motifs were identified in promoter regions of all gene sets using the binding site estimation suite of tools (BEST). Comparative analysis of these motifs revealed similarities with known cis -elements found in PLACE database and led to the discovery of cis -regulatory motifs that were not yet implicated in drought response. Compiled with the proposed synthetic promoter design rationale, three synthetic assemblies were constructed by concatenating multiple copies of drought-inducible cis motifs in a specific order with inter-motif spacing using random bases and placed upstream of 35s minimal core promoter. Each synthetic module substituted 35S promoter in pBI121 and pCAMBIA3301 to drive glucuronidase expression in soybean hairy roots and Arabidopsis thaliana L. Chimeric soybean seedlings and 3-week-old transgenic Arabidopsis plants were treated with simulated with different levels of osmotic stress. Histochemical staining of transgenic soybean hairy roots and Arabidopsis displayed drought-inducible GUS activity of synthetic promoters. Fluorometric assay and expression analysis revealed that SP2 is the better manual combination of cis -elements for stress-inducible expression. qRT-PCR results further demonstrated that designed synthetic promoters are not tissue-specific and thus active in different parts upon treatment with osmotic stress in Arabidopsis plants. This study provides tools for transcriptional upgradation of valuable crops against drought stress and adds to the current knowledge of synthetic biology.

Genome-Wide Identification and Salt Stress Response Analysis of the bZIP Transcription Factor Family in Sugar Beet

As one of the largest transcription factor families in plants, bZIP transcription factors play important regulatory roles in different biological processes, especially in the process of stress response. Salt stress inhibits the growth and yield of sugar beet. However, bZIP-related studies in sugar beet (Beta vulgaris L.) have not been reported. This study aimed to identify the bZIP transcription factors in sugar beet and analyze their biological functions and response patterns to salt stress. Using bioinformatics, 48 BvbZIP genes were identified in the genome of sugar beet, encoding 77 proteins with large structural differences. Collinearity analysis showed that three pairs of BvbZIP genes were fragment replication genes. The BvbZIP genes were grouped according to the phylogenetic tree topology and conserved structures, and the results are consistent with those reported in Arabidopsis. Under salt stress, the expression levels of most BvbZIP genes were decreased, and only eight genes were up-regulated. GO analysis showed that the BvbZIP genes were mainly negatively regulated in stress response. Protein interaction prediction showed that the BvbZIP genes were mainly involved in light signaling and ABA signal transduction, and also played a certain role in stress responses. In this study, the structures and biological functions of the BvbZIP genes were analyzed to provide foundational data for further mechanistic studies and for facilitating the efforts toward the molecular breeding of stress-resilient sugar beet.

Transcriptome data reveal gene clusters and key genes in pepper response to heat shock

Climate change and global warming pose a great threat to plant growth and development as well as crop productivity. To better study the genome-wide gene expression under heat, we performed a time-course (0.5 to 24 h) transcriptome analysis in the leaf and root of 40-day-old pepper plants under 40°C as well as in control plants. Clustering analysis (K-means) showed that the expression of 29,249 genes can be grouped into 12 clusters with distinct expression dynamics under stress. Gene ontology (GO) enrichment analysis and transcription factor (TF) identification were performed on the clusters with certain expression patterns. Comparative analysis between the heat-treated and control plants also identified differentially expressed genes (DEGs), which showed the largest degree of change at 24 h. Interestingly, more DEGs were identified in the root than in the leaf. Moreover, we analyzed the gene expression of 25 heat shock factor genes (HSFs) in pepper after heat stress, identified five of these HSFs that responded to heat stress, and characterized the role of these genes in heat-tolerant (17CL30) and heat-susceptible (05S180) pepper lines. The findings of this study improve our understanding of the genome-wide heat stress response in pepper.

BGFD: an integrated multi-omics database of barley gene families

BACKGROUND

A gene family comprises a group of genes with similar functional domains that play various roles in plant growth, development, and responses to environmental stimuli. Barley (Hordeum vulgare L.) is the fourth most cultivated cereal crop worldwide, and it is an important model species for genetic studies. Systematic identification and annotation of gene families are key for studies of molecular function and evolutionary history.

RESULTS

We constructed a multi-omics database containing 5593 genes of 77 gene families called the Barley Gene Family Database (BGFD: http://barleygfdb.com ). BGFD is a free, user-friendly, and web-accessible platform that provides data on barley family genes. BGFD provides intuitive visual displays to facilitate studies of the physicochemical properties, gene structure, phylogenetic relationships, and motif organization of genes. Massive multi-omics datasets have been acquired and processed to generate an atlas of expression pattern profiles and genetic variation in BGFD. The platform offers several practical toolkits to conduct searches, browse, and employ BLAST functions, and the data are downloadable.

CONCLUSIONS

BGFD will aid research on the domestication and adaptive evolution of barley; it will also facilitate the screening of candidate genes and exploration of important agronomic traits in barley.

Genome-wide identification and characterization of AP2/ERF gene superfamily during flower development in Actinidia eriantha

BACKGROUND

As one of the largest transcription factor families in plants, AP2/ERF gene superfamily plays important roles in plant growth, development, fruit ripening and biotic and abiotic stress responses. Despite the great progress has been made in kiwifruit genomic studies, little research has been conducted on the AP2/ERF genes of kiwifruit. The increasing kiwifruit genome resources allowed us to reveal the tissue expression profiles of AP2/ERF genes in kiwifruit on a genome-wide basis.

RESULTS

In present study, a total of 158 AP2/ERF genes in A. eriantha were identified. All genes can be mapped on the 29 chromosomes. Phylogenetic analysis divided them into four main subfamilies based on the complete protein sequences. Additionally, our results revealed that the same subfamilies contained similar gene structures and conserved motifs. Ka/Ks calculation indicated that AP2/ERF gene family was undergoing a strong purifying selection and the evolutionary rates were slow. RNA-seq showed that the AP2/ERF genes were expressed differently in different flower development stages and 56 genes were considered as DEGs among three contrasts. Moreover, qRT-PCR suggested partial genes showed significant expressions as well, suggesting they could be key regulators in flower development in A. eriantha. In addition, two genes (AeAP2/ERF061, AeAP2/ERF067) had abundant transcription level based on transcriptomes, implying that they may play a crucial role in plant flower development regulation and flower tissue forming.

CONCLUSIONS

We identified AP2/ERF genes and demonstrated their gene structures, conserved motifs, and phylogeny relationships of AP2/ERF genes in two related species of kiwifruit, A. eriantha and A. chinensis, and their potential roles in flower development in A. eriantha. Such information would lay the foundation for further functional identification of AP2/ERF genes involved in kiwifruit flower development.