Genome-wide identification and characterization of MADS-box family genes related to organ development and stress resistance in Brassica rapa

The Multi-Pistil Phenomenon in Higher Plants

Correct floral morphology determines the accuracy of fruit formation, which is crucial for reproductive success in higher plants. Despite this, an abnormal, multi-pistil phenotype has been observed in the flowers of many plants. In this review, we gather information on the multi-pistil phenomenon in various species and highlight potential causes, as well as possible consequences, of the trait. Our assessment of the reported multi-pistil phenotype in rice (Oryza sativa L.), wheat (Triticum aestivum L.), tomato (Solanum lycopersicum L.), Medicago, sweet cherry (Prunus avium L.), rye (Secale cereale L.), and rapeseed (Brassica napus L. and B. campestris L.) leads us to conclude that hybridization and mutation are the main factors that give rise to this phenotype. We also delve into the inheritance patterns of the multi-pistil phenotype and factors that influence this trait, such as nuclear-cytoplasmic interactions, temperature conditions, and shading. Finally, we discuss the effects of multi-pistil flowers on the yield of these plants. This analysis increases our understanding of floral development and lays the foundation for the potential utilization of the multi-pistil trait to increase seed production in crops.

Genome-wide characterization and expression profiling of MADS-box family genes during organ development and drought stress in Camelina sativa L

Camelina sativa stands out among oilseed crops due to its remarkable resilience to challenging environmental conditions and its significant potential for biodiesel production. The MADS-box transcription factors play a pivotal role in numerous biological processes within plants, encompassing growth, development, and responses to environmental stressors. In this research, by employing the BLAST, we have successfully identified 325 MADS-box genes within Camelina sativa genome. These genes were systematically categorized into two principal groups: type I (comprising Mα, Mβ, and Mγ) and type II (including MI-KCC and MIKC*) predicated upon their phylogenetic relationships, structural protein motifs, and exon-intron configurations. Our findings reveal that type II MADS-box genes have, in general, experienced a more profound expansion relative to type I genes. Specifically, the TM3 subgroup within type II MADS-box genes exhibited the highest degree of gene expansion, comprising 21 TM3 genes. The amino acid sequences encoded by these genes exhibited a length variation ranging from 150 to 820 aa. The predicted molecular weights (MW) of the CsMADS-box proteins displayed a range from 17.01 to 94.06 kDa, while the isoelectric points (pI) were observed to span from 4.13 to 10.09. Evolutionary analysis predicated on the Ka/Ks ratios indicates that the evolutionary pathway of MADS-box genes in Camelina sativa has been predominantly driven by the mechanisms of purifying selection. Moreover, an investigation of cis-acting elements has elucidated the participation of MADS-box genes in the adaptive responses to abiotic stressors. The expression profiles of six Type I and three MIKCC genes across diverse organs and under varying drought treatment conditions demonstrated that these genes are expressed in both reproductive and vegetative structures, displaying uniform expression patterns throughout several developmental phases of flowering. The expression levels of CsMADS035, CsMADS115, CsMADS131, and CsMADS181 were notably modified in reaction to drought stress conditions. The detailed annotation and comprehensive transcriptome profiling provided in this research yield essential insights into the functional roles that MADS-box genes perform in stress resistance, as well as their contributions to growth and developmental processes. This acquired understanding establishes a foundational framework for the functional characterization and potential genetic engineering initiatives pertaining to Camelina sativa, thereby augmenting the prospective application of these candidate genes.

Omics approaches to understand the MADS-box gene family in common bean (Phaseolus vulgaris L.) against drought stress

MADS-box genes are known to play important roles in diverse aspects of growth/devolopment and stress response in several plant species. However, no study has yet examined about MADS-box genes in P. vulgaris. In this study, a total of 79 PvMADS genes were identified and classified as type I and type II according to the phylogenetic analysis. While both type I and type II PvMADS classes were found to contain the MADS domain, the K domain was found to be present only in type II PvMADS proteins, in agreement with the literature. All chromosomes of the common bean were discovered to contain PvMADS genes and 17 paralogous gene pairs were identified. Only two of them were tandemly duplicated gene pairs (PvMADS-19/PvMADS-23 and PvMADS-20/PvMADS-24), and the remaining 15 paralogous gene pairs were segmentally duplicated genes. These duplications were found to play an important role in the expansion of type II PvMADS genes. Moreover, the RNAseq and RT-qPCR analyses showed the importance of PvMADS genes in response to drought stress in P. vulgaris.

Characterization and validation of TaAGL66, a gene related to fertility conversion of wheat in the presence of Aegilops kotschyi cytoplasm

MAIN CONCLUSION

TaAGL66, a MADS-box transcription factor highly expressed in fertile anthers of KTM3315A, regulates anther and/or pollen development, as well as male fertility in wheat with Aegilops kotschyi cytoplasm. Male sterility, as a string of sophisticated biological processes in higher plants, is commonly regulated by transcription factors (TFs). Among them, MADS-box TFs are mainly participated in the processes of floral organ formation and pollen development, which are tightly related to male sterility, but they have been little studied in the reproductive development in wheat. In our study, TaAGL66, a gene that was specifically expressed in spikes and highly expressed in fertile anthers, was identified by RNA sequencing and the expression profiles data of these genes, and qRT-PCR analyses, which was localized to the nucleus. Silencing of TaAGL66 under fertility condition in KTM3315A, a thermo-sensitive male sterile line with Ae. kotschyi cytoplasm, displayed severe fertility reduction, abnormal anther dehiscence, defective pollen development, decreased viability, and low seed-setting. It can be concluded that TaAGL66 plays an important role in wheat pollen development in the presence of Ae. kotschyi cytoplasm, providing new insights into the utilization of male sterility.

Identification, Evolutionary Dynamics, and Gene Expression Patterns of the ACP Gene Family in Responding to Salt Stress in Brassica Genus

Acyl carrier proteins (ACPs) have been reported to play a crucial role in responding to biotic and abiotic stresses, regulating growth and development. However, the biological function of the ACP gene family in the Brassica genus has been limited until now. In this study, we conducted a comprehensive analysis and identified a total of 120 ACP genes across six species in the Brassica genus. Among these, there were 27, 26, and 30 ACP genes in the allotetraploid B. napus, B. juncea, and B. carinata, respectively, and 14, 13, and 10 ACP genes in the diploid B. rapa, B. oleracea, and B. nigra, respectively. These ACP genes were further classified into six subclades, each containing conserved motifs and domains. Interestingly, the majority of ACP genes exhibited high conservation among the six species, suggesting that the genome evolution and polyploidization processes had relatively minor effects on the ACP gene family. The duplication modes of the six Brassica species were diverse, and the expansion of most ACPs in Brassica occurred primarily through dispersed duplication (DSD) events. Furthermore, most of the ACP genes were under purifying selection during the process of evolution. Subcellular localization experiments demonstrated that ACP genes in Brassica species are localized in chloroplasts and mitochondria. Cis-acting element analysis revealed that most of the ACP genes were associated with various abiotic stresses. Additionally, RNA-seq data revealed differential expression levels of BnaACP genes across various tissues in B. napus, with particularly high expression in seeds and buds. qRT-PCR analysis further indicated that BnaACP genes play a significant role in salt stress tolerance. These findings provide a comprehensive understanding of ACP genes in Brassica plants and will facilitate further functional analysis of these genes.

Characterization of the MIKCC-type MADS-box gene family in blueberry and its possible mechanism for regulating flowering in response to the chilling requirement

MAIN CONCLUSION

The expression peak of VcAP1.4, VcAP1.6, VcAP3.1, VcAP3.2, VcAG3, VcFLC2, and VcSVP9 coincided with the endo-dormancy release of flower buds. Additionally, GA4+7 not only increased the expression of these genes but also promoted flower bud endo-dormancy release. The MIKCC-type MADS-box gene family is involved in the regulation of flower development. A total of 109 members of the MIKCC-type MADS-box gene family were identified in blueberry. According to the phylogenetic tree, these 109 MIKCC-type MADS-box proteins were divided into 13 subfamilies, which were distributed across 40 Scaffolds. The results of the conserved motif analysis showed that among 20 motifs, motifs 1, 3, and 9 formed the MADS-box structural domain, while motifs 2, 4, and 6 formed the K-box structural domain. The presence of 66 pairs of fragment duplication events in blueberry suggested that gene duplication events contributed to gene expansion and functional differentiation. Additionally, the presence of cis-acting elements revealed that VcFLC2, VcAG3, and VcSVP9 might have significant roles in the endo-dormancy release of flower buds. Meanwhile, under chilling conditions, VcAP3.1 and VcAG7 might facilitate flower bud dormancy release. VcSEP11 might promote flowering following the release of endo-dormancy, while the elevated expression of VcAP1.7 (DAM) could impede the endo-dormancy release of flower buds. The effect of gibberellin (GA4+7) treatment on the expression pattern of MIKCC-type MADS-box genes revealed that VcAP1.4, VcAP1.6, VcAP3.1, VcAG3, and VcFLC2 might promote flower bud endo-dormancy release, while VcAP3.2, VcSEP11, and VcSVP9 might inhibit its endo-dormancy release. These results indicated that VcAP1.4, VcAP1.6, VcAP1.7 (DAM), VcAP3.1, VcAG3, VcAG7, VcFLC2, and VcSVP9 could be selected as key regulatory promoting genes for controlling the endo-dormancy of blueberry flower buds.

Physiological and molecular characteristics associated with the anti-senescence in Camellia oleifera Abel

This study analyzed physiological and molecular characteristics associated with the resistance to aging or anti-senescence in Camellia oleifera Abel. Trees over 100 years old (ancient trees) were compared with those about 30 years old (mature trees). Total chlorophylls, chlorophyll a/b ratio, and hydrogen peroxide concentrations in ancient tree leaves were significantly higher than in their counterparts. Significantly higher activities of superoxide dismutase, peroxidase, and catalase were detected in ancient tree leaves. Nine Chl a/b-binding protein genes, 15 antioxidant enzyme genes, 21 hormone-related genes, and 301 stress-related genes were upregulated, and 42 protein-degradation genes were downregulated in ancient tree leaves. By increasing chlorophyll content and antioxidant enzyme activities and regulating the ageing-related genes expression, ancient C. oleifera leaves maintained remarkable vitality. Although further research is needed, our study may shed some light on how ancient C. oleifera trees can resist ageing and sustain their healthy growth.

Genome-Wide Characterization and Phylogenetic and Stress Response Expression Analysis of the MADS-Box Gene Family in Litchi (Litchi chinensis Sonn.)

The MADS-box protein is an important transcription factor in plants and plays an important role in regulating the plant abiotic stress response. In this study, a total of 94 MADS-box genes were predicted in the litchi genome, and these genes were widely distributed on all the chromosomes. The LcMADS-box gene family was divided into six subgroups (Mα, Mβ, Mγ, Mδ, MIKC, and UN) based on their phylogenetical relationships with Arabidopsis, and the closely linked subgroups exhibited more similarity in terms of motif distribution and intron/exon numbers. Transcriptome analysis indicated that LcMADS-box gene expression varied in different tissues, which can be divided into universal expression and specific expression. Furthermore, we further validated that LcMADS-box genes can exhibit different responses to various stresses using quantitative real-time PCR (qRT-PCR). Moreover, physicochemical properties, subcellular localization, collinearity, and cis-acting elements were also analyzed. The findings of this study provide valuable insights into the MADS-box gene family in litchi, specifically in relation to stress response. The identification of hormone-related and stress-responsive cis-acting elements in the MADS-box gene promoters suggests their involvement in stress signaling pathways. This study contributes to the understanding of stress tolerance mechanisms in litchi and highlights potential regulatory mechanisms underlying stress responses.

Genome-Wide Identification and Expression Analysis of the MADS Gene Family in Tulips (Tulipa gesneriana)

To investigate the cold response mechanism and low temperature regulation of flowering in tulips, this study identified 32 MADS-box transcription factor family members in tulips based on full-length transcriptome sequencing, named TgMADS1-TgMADS32. Phylogenetic analysis revealed that these genes can be divided into two classes: type I and type II. Structural analysis showed that TgMADS genes from different subfamilies have a similar distribution of conserved motifs. Quantitative real-time PCR results demonstrated that some TgMADS genes (e.g., TgMADS3, TgMADS15, TgMADS16, and TgMADS19) were significantly upregulated in buds and stems under cold conditions, implying their potential involvement in the cold response of tulips. In summary, this study systematically identified MADS family members in tulips and elucidated their evolutionary relationships, gene structures, and cold-responsive expression patterns, laying the foundation for further elucidating the roles of these transcription factors in flowering and the cold adaptability of tulips.

Genome-wide analysis of the MADS-box gene family in Lonicera japonica and a proposed floral organ identity model

BACKGROUND

Lonicera japonica Thunb. is widely used in traditional Chinese medicine. Medicinal L. japonica mainly consists of dried flower buds and partially opened flowers, thus flowers are an important quality indicator. MADS-box genes encode transcription factors that regulate flower development. However, little is known about these genes in L. japonica.

RESULTS

In this study, 48 MADS-box genes were identified in L. japonica, including 20 Type-I genes (8 Mα, 2 Mβ, and 10 Mγ) and 28 Type-II genes (26 MIKCc and 2 MIKC*). The Type-I and Type-II genes differed significantly in gene structure, conserved domains, protein structure, chromosomal distribution, phylogenesis, and expression pattern. Type-I genes had a simpler gene structure, lacked the K domain, had low protein structure conservation, were tandemly distributed on the chromosomes, had more frequent lineage-specific duplications, and were expressed at low levels. In contrast, Type-II genes had a more complex gene structure; contained conserved M, I, K, and C domains; had highly conserved protein structure; and were expressed at high levels throughout the flowering period. Eleven floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in L. japonica. By integrating expression pattern and protein interaction data for these genes, we developed a possible model for floral organ identity determination.

CONCLUSION

This study genome-widely identified and characterized the MADS-box gene family in L. japonica. Eleven floral homeotic MADS-box genes were identified and a possible model for floral organ identity determination was also developed. This study contributes to our understanding of the MADS-box gene family and its possible involvement in floral organ development in L. japonica.

Expression Characterization of ABCDE Class MADS-Box Genes in Brassica rapa with Different Pistil Types

MADS-box is a vital transcription factor family that functions in plant growth and development. Apart from APETALA2, all genes in the ABCDE model that explain the molecular mechanism of floral organ development belong to the MADS-box family. Carpel and ovule numbers in plants are essential agronomic traits that determine seed yield, and multilocular siliques have great potential for the development of high-yield varieties of Brassica. In this study, ABCDE genes in the MADS-box family from Brassica rapa were identified and characterized. Their tissue-specific expression patterns in floral organs and their differential expression in different pistil types of B. rapa were revealed by qRT-PCR. A total of 26 ABCDE genes were found to belong to the MADS-box family. Our proposed ABCDE model of B. rapa is consistent with that of Arabidopsis thaliana, indicating that ABCDE genes are functionally conserved. These results of qRT-PCR showed that the expression levels of class C and D genes were significantly different between the wild-type (wt) and tetracarpel (tetrac) mutant of B. rapa. Interestingly, the expression of the homologs of class E genes was imbalanced. Therefore, it is speculated that class C, D, and E genes are involved in developing the carpel and ovule of B. rapa. Our findings reveal the potential for the selection of candidate genes to improve yield traits in Brassica crops.

Molecular evolution, diversification, and expression assessment of MADS gene family in Setaria italica, Setaria viridis, and Panicum virgatum

KEY MESSAGE

This paper sheds light on the evolution and expression patterns of MADS genes in Setaria and Panicum virgatum. SiMADS51 and SiMADS64 maybe involved in the ABA-dependent pathway of drought response. The MADS gene family is a key regulatory factor family that controls growth, reproduction, and response to abiotic stress in plants. However, the molecular evolution of this family is rarely reported. Here, a total of 265 MADS genes were identified in Setaria italica (foxtail millet), Setaria viridis (green millet), and Panicum virgatum (switchgrass) and analyzed by bioinformatics, including physicochemical characteristics, subcellular localization, chromosomal position and duplicate, motif distribution, genetic structure, genetic evolvement, and expression patterns. Phylogenetic analysis was used to categorize these genes into M and MIKC types. The distribution of motifs and gene structure were similar for the corresponding types. According to a collinearity study, the MADS genes have been mostly conserved during evolution. The principal cause of their expansion is segmental duplication. However, the MADS gene family tends to shrink in foxtail millet, green millet, and switchgrass. The MADS genes were subjected to purifying selection, but several positive selection sites were also identified in three species. And most of the promoters of MADS genes contain cis-elements related to stress and hormonal response. RNA-seq and quantitative Real-time PCR (qRT-PCR) analysis also were examined. SiMADS genes expression levels are considerably changed in reaction to various treatments, following qRT-PCR analysis. This sheds fresh light on the evolution and expansion of the MADS family in foxtail millet, green millet, and switchgrass, and lays the foundation for further research on their functions.

Genome-wide association study reveals genomic loci influencing agronomic traits in Ethiopian sorghum (Sorghum bicolor (L.) Moench) landraces

Uncovering the genetic basis of agronomic traits in sorghum landraces that have adapted to various agro-climatic conditions would contribute to sorghum improvement efforts around the world. To identify quantitative trait nucleotides (QTNs) associated with nine agronomic traits in a panel of 304 sorghum accessions collected from diverse environments across Ethiopia (considered to be the center of origin and diversity), multi-locus genome-wide association studies (ML-GWAS) were performed using 79,754 high quality single nucleotide polymorphism (SNP) markers. Association analyses using six ML-GWAS models identified a set of 338 significantly (LOD ≥ 3)-associated QTNs for nine agronomic traits of sorghum accessions evaluated in two environments (E1 and E2) and their combined dataset (Em). Of these, 121 reliable QTNs, including 13 for flowering time (DF), 13 for plant height (PH), 9 for tiller number (TN), 15 for panicle weight (PWT), 30 for grain yield per panicle (GYP), 12 for structural panicle mass (SPM), 13 for hundred seed weight (HSW), 6 for grain number per panicle (GNP), and 10 for panicle exertion (PE) were consistently detected by at least three ML-GWAS methods and/or in two different environments. Notably, Ethylene responsive transcription factor gene AP2/ERF, known for regulation of plant growth, and the sorghum Terminal flower1/TF1 gene, which functions in the control of floral architecture, were identified as strong candidate genes associated with PH and HSW, respectively. This study provides an entry point for further validation studies to elucidate complex mechanisms controlling important agronomic traits in sorghum.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01381-5.

B-type response regulator hst1 controls salinity tolerance in rice by regulating transcription factors and antioxidant mechanisms

Salinity is a serious environmental problem that limits plant yield in almost half of the agricultural fields. The hitomebore salt tolerant 1(hst1) is a mutant B-type response regulator gene that was reported to improve salinity tolerance in the 'YNU31-2-4' (YNU) genotype. The sister line (SL) is salt-sensitive, and the nearest genomic relative of the YNU plant has the OsRR22 gene, which is the non-mutant form of the hst1 gene. Biochemical and comprehensive transcriptome analysis of SL and YNU plants was performed to clarify the salinity tolerance mechanism(s) mediated by the hst1 gene. The hst1 gene reduced Na⁺ ions, lipid peroxidation, and H₂O₂ content, and improve proline and antioxidant enzymes activities under salt stress. Various transporter and gene-specific transcriptional regulator genes up-regulated in presence of the hst1 gene under saline conditions, identifying that post-stress transcription factors (OsbHLH056, OsH43, OsGRAS29, and OsMADS1) contributed to improved salinity tolerance in YNU plants. Specifically, OsSalT, miR156, and OsLPT1.16 genes were up-regulated, while upstream (OsHKs and OsHPs) and downstream regulators of the OsRR22 gene were down-regulated in YNU plants under saline conditions. Notably, the transcription factors reprogramming, upstream and downstream genes, indicate that these pathways are transcriptionally regulated by the hst1 gene. The findings of the regulatory role of the hst1 gene on plant transcriptome provide a greater understanding of hst1-mediated salt tolerance in rice plants. This knowledge will contribute to understanding the salinity tolerance mechanisms in rice and the evolution of salt-tolerant crops with the ability to withstand higher salinity to ensure food security during climate change.

Genome-Wide Analysis of the Mads-Box Transcription Factor Family in Solanum melongena

The MADS-box transcription factors are known to be involved in several aspects of plant growth and development, especially in floral organ specification. However, little is known in eggplant. Here, 120 eggplant MADS-box genes were identified and categorized into type II (MIKCC and MIKC*) and type I (Mα, Mβ, and Mγ) subfamilies based on phylogenetic relationships. The exon number in type II SmMADS-box genes was greater than that in type I SmMADS-box genes, and the K-box domain was unique to type II MADS-box TFs. Gene duplication analysis revealed that segmental duplications were the sole contributor to the expansion of type II genes. Cis-elements of MYB binding sites related to flavonoid biosynthesis were identified in three SmMADS-box promoters. Flower tissue-specific expression profiles showed that 46, 44, 38, and 40 MADS-box genes were expressed in the stamens, stigmas, petals, and pedicels, respectively. In the flowers of SmMYB113-overexpression transgenic plants, the expression levels of 3 SmMADS-box genes were co-regulated in different tissues with the same pattern. Correlation and protein interaction predictive analysis revealed six SmMADS-box genes that might be involved in the SmMYB113-regulated anthocyanin biosynthesis pathway. This study will aid future studies aimed at functionally characterizing important members of the MADS-box gene family.

Genome-wide identification of MADS-box gene family in orchardgrass and the positive role of DgMADS114 and DgMADS115 under different abiotic stress

Abiotic stress, a major factor limit growth and productivity of major crops. Orchardgrass is one of the most important cool-season forage grasses in the world, and it is highly tolerant to abiotic stress. The MADS-box transcription factor family is one of the largest families in plants, and it plays vital roles in multiple biological processes. However, MADS-box transcription factors in orchardgrass, especially those involved in abiotic stress, have not yet been elucidated. Here, 123 DgMADS-box members were identified in orchardgrass and a detailed overview has been presented. Syntenic analysis indicated that the expansion of the DgMADS-box genes in orchardgrass is mainly dependent on tandem duplication events. Some DgMADS-box genes were induced by multiple abiotic stresses, indicating that these genes may play critical regulatory roles in orchardgrass response to various abiotic stresses. Heterologous expression showed that DgMADS114 and DgMADS115 could enhance stress tolerance of transgenic Arabidopsis, as revealed by longer root length or higher survival rates under PEG, NaCl, ABA, and heat stress. The results of this study provide a scientific basis for clarifying the functional characterization of MADS-box genes in orchardgrass in response to environmental stress can be further used to improve forages and crops via breeding programs.

Functional Characterization of the Cystine-Rich-Receptor-like Kinases (CRKs) and Their Expression Response to Sclerotinia sclerotiorum and Abiotic Stresses in Brassica napus

Cysteine-rich receptor-like kinases (CRKs) are transmembrane proteins that bind to the calcium ion to regulate stress-signaling and plant development-related pathways, as indicated by several pieces of evidence. However, the CRK gene family hasn’t been inadequately examined in Brassica napus. In our study, 27 members of the CRK gene family were identified in Brassica napus, which are categorized into three phylogenetic groups and display synteny relationship to the Arabidopsis thaliana orthologs. All the CRK genes contain highly conserved N-terminal PKINASE domain; however, the distribution of motifs and gene structure were variable conserved. The functional divergence analysis between BnaCRK groups indicates a shift in evolutionary rate after duplication events, demonstrating that BnaCRKs might direct a specific function. RNA-Seq datasets and quantitative real-time PCR (qRT-PCR) exhibit the complex expression profile of the BnaCRKs in plant tissues under multiple stresses. Nevertheless, BnaA06CRK6-1 and BnaA08CRK8 from group B were perceived to play a predominant role in the Brassica napus stress signaling pathway in response to drought, salinity, and Sclerotinia sclerotiorum infection. Insights gained from this study improve our knowledge about the Brassica napus CRK gene family and provide a basis for enhancing the quality of rapeseed.

Genome-wide analysis of MADS-box gene family in kiwifruit (Actinidia chinensis var. chinensis) and their potential role in floral sex differentiation

Kiwifruit (Actinidia chinensis Planch.) is a functionally dioecious plant, which displays diverse morphology in male and female flowers. MADS-box is an ancient and huge gene family that plays a key role in plant floral organ differentiation. In this study, we have identified 89 MADS-box genes from A. chinensis Red 5 genome. These genes are distributed on 26 chromosomes and are classified into type I (21 genes) and type II (68 genes). Overall, type II AcMADS-box genes have more complex structures than type I with more exons, protein domains, and motifs, indicating that type II genes may have more diverse functions. Gene duplication analysis showed that most collinearity occurred in type II AcMADS-box genes, which was consistent with a large number of type II genes. Analysis of cis-acting elements in promoters showed that AcMADS-box genes are mainly associated with light and phytohormone responsiveness. The expression profile of AcMADS-box genes in different tissues showed that most genes were highly expressed in flowers. Further, the qRT-PCR analysis of the floral organ ABCDE model-related genes in male and female flowers revealed that AcMADS4, AcMADS56, and AcMADS70 were significantly expressed in female flowers. It indicated that those genes may play an important role in the sex differentiation of kiwifruit. This work provided a comprehensive analysis of the AcMADS-box genes and may help facilitate our understanding of the sex differentiation regulatory mechanism in kiwifruit.

Genome-wide identification and expression analysis of MIKCC genes in rose provide insight into their effects on flower development

The MIKC^C-type gene family plays important roles in plant growth, development, and tolerance of biotic and abiotic stress, especially during floral organ differentiation. However, there have been no studies of MIKC^C-type genes in rose, and functional differentiation of family members has not been explored. In this study, we identified 42 MIKC^C-type genes in rose, classified the genes into 12 subfamilies, and constructed a phylogenetic tree. We performed expression analysis of these genes, and found that expression patterns correlated with the predicted subfamily, indicating that the features of MIKC^C-type genes were broadly conserved during evolution. Collinear analysis of MIKC^C genes among Rosaceae species confirmed the occurrence of whole genome duplications (WGD) and revealed some species-specific MIKC^C genes. Transcriptome analysis showed that the expression of some MIKC^C-type genes responded to low temperatures (4°C, 24 h) during flower organ differentiation. These conserved, duplicated, and novel expression patterns of MIKC^C-type genes may have facilitated the adaptation of rose to various internal and external environmental changes. The results of this study provide a theoretical basis for future functional analysis of the MIKC^C genes in rose and investigation of the evolutionary pattern of the MIKC^C gene family in the Rosaceae genome.

Transcription-associated metabolomic profiling reveals the critical role of frost tolerance in wheat

BACKGROUND

Low temperature is a crucial stress factor of wheat (Triticum aestivum L.) and adversely impacts on plant growth and grain yield. Multi-million tons of grain production are lost annually because crops lack the resistance to survive in winter. Particularlly, winter wheat yields was severely damaged under extreme cold conditions. However, studies about the transcriptional and metabolic mechanisms underlying cold stresses in wheat are limited so far.

RESULTS

In this study, 14,466 differentially expressed genes (DEGs) were obtained between wild-type and cold-sensitive mutants, of which 5278 DEGs were acquired after cold treatment. 88 differential accumulated metabolites (DAMs) were detected, including P-coumaroyl putrescine of alkaloids, D-proline betaine of mino acids and derivativ, Chlorogenic acid of the Phenolic acids. The comprehensive analysis of metabolomics and transcriptome showed that the cold resistance of wheat was closely related to 13 metabolites and 14 key enzymes in the flavonol biosynthesis pathway. The 7 enhanced energy metabolites and 8 up-regulation key enzymes were also compactly involved in the sucrose and amino acid biosynthesis pathway. Moreover, quantitative real-time PCR (qRT-PCR) revealed that twelve key genes were differentially expressed under cold, indicating that candidate genes POD, Tacr7, UGTs, and GSTU6 which were related to cold resistance of wheat.

CONCLUSIONS

In this study, we obtained the differentially expressed genes and differential accumulated metabolites in wheat under cold stress. Using the DEGs and DAMs, we plotted regulatory pathway maps of the flavonol biosynthesis pathway, sucrose and amino acid biosynthesis pathway related to cold resistance of wheat. It was found that candidate genes POD, Tacr7, UGTs and GSTU6 are related to cold resistance of wheat. This study provided valuable molecular information and new genetic engineering clues for the further study on plant resistance to cold stress.

Genome-wide identification, phylogenetic and expression pattern analysis of MADS-box family genes in foxtail millet (Setaria italica)

Foxtail millet (Setaria italica) is rich in nutrients and extremely beneficial to human health. We identified and comprehensively analyzed 89 MADS-box genes in the foxtail millet genome. According to the classification of MADS-box genes in Arabidopsis thaliana and rice, the SiMADS-box genes were divided into M-type (37) and MIKC-type (52). During evolution, the differentiation of MIKC-type MADS-box genes occurred before that of monocotyledons and dicotyledons. The SiMADS-box gene structure has undergone much differentiation, and the number of introns in the MIKC-type subfamily is much greater than that in the M-type subfamily. Analysis of gene duplication events revealed that MIKC-type MADS-box gene segmental duplication accounted for the vast majority of gene duplication events, and MIKC-type MADS-box genes played a major role in the amplification of SiMADS-box genes. Collinearity analysis showed highest collinearity between foxtail millet and maize MADS-box genes. Analysis of tissue-specific expression showed that SiMADS-box genes are highly expressed throughout the grain-filling process. Expression analysis of SiMADS-box genes under eight different abiotic stresses revealed many stress-tolerant genes, with induced expression of SiMADS33 and SiMADS78 under various stresses warranting further attention. Further, some SiMADS-box proteins may interact under external stress. This study provides insights for MADS-box gene mining and molecular breeding of foxtail millet in the future.

The MADS transcription factor GhFYF is involved in abiotic stress responses in upland cotton (Gossypium hirsutum L.)

Cotton is an important textile industry raw material crops, which plays a critical role in the development of society. MADS transcription factors (TFs) play a key role about the flowering time, flower development, and abiotic stress responses in plants, but little is known about their functions on abiotic stress in cotton. In this study, a MIKC^C subfamily gene from cotton, GhFYF (FOREVER YOUNG FLOWER), was isolated and characterized. Our data showed that GhFYF localized to the nucleus. A β-glucuronidase (GUS) activity assay revealed that the promoter of GhFYF was mainly expressed in the flower and seed of ProGhFYF::GUS transgenic A. thaliana plants. The GUS staining of flowers and seeds was deepened after drought, salt treatment, and the expression level of the GUS gene and corresponding stress genes AtERD10, AtAnnexin1 are up-regulated in the inflorescence. Overexpression GhFYF in A. thaliana could promote the seed germination and growth under different salt concentrations, and determin the proline content. Yeast two-hybrid (Y2H) assays showed that GhFYF interacted with the HAD-like protein GhGPP2, which has responds to abiotic stress. Our findings indicate that GhFYF is involved in abiotic stress responses, especially for salt stress. This work establishes a solid foundation for further functional analysis of the GhFYF gene in cotton.

A rice transcription factor, OsMADS57, positively regulates high salinity tolerance in transgenic Arabidopsis thaliana and Oryza sativa plants

MADS-box transcription factors (TFs) play indispensable roles in various aspects of plant growth, development as well as in response to environmental stresses. Several MADS-box genes have been reported to be involved in the salt tolerance in different plant species. However, the role of the transcription factor OsMADS57 under salinity stress is still unknown. Here, the results of this study showed that OsMADS57 was mainly expressed in roots and leaves of rice plants (Oryza sativa). Gene expression pattern analysis revealed that OsMADS57 was induced by NaCl. Overexpression of OsMADS57 in both Arabidopsis thaliana (A. thaliana) and rice could improve their salt tolerance, which was demonstrated by higher germination rates, longer root length and better growth status of overexpression plants than wild type (WT) under salinity conditions. In contrast, RNA interference (RNAi) lines of rice showed more sensitivity towards salinity. Moreover, less reactive oxygen species (ROS) accumulated in OsMADS57 overexpressing lines when exposed to salt stress, as measured by 3, 3'-diaminobenzidine (DAB) or nitroblue tetrazolium (NBT) staining. Further experiments exhibited that overexpression of OsMADS57 in rice significantly increased the tolerance ability of plants to oxidative damage under salt stress, mainly by increasing the activities of antioxidative enzymes such as superoxide dismutase (SOD) and peroxidase (POD), reducing malonaldehyde (MDA) content and improving the expression of stress-related genes. Taken together, these results demonstrated that OsMADS57 plays a positive role in enhancing salt tolerance by activating the antioxidant system.

Genome-wide identification of the MADS-box transcription factor family in autotetraploid cultivated alfalfa (Medicago sativa L.) and expression analysis under abiotic stress

Background

Alfalfa, the “queen of forage”, is the most extensively cultivated forage legume in the world. The development and yield of alfalfa are seriously limited by abiotic stress. MADS-box transcription factors are one of the largest gene families and play a pivotal role in plant development and abiotic stress. However, little is known regarding the MADS-box transcription factors in autotetraploid cultivated alfalfa.

Results

In the present study, we identified 120 MsMADS-box genes in the alfalfa genome. Phylogenetic analysis indicated that 75 type-I MsMADS-box genes were classified into the Mα, Mβ, and Mγ subgroups, and 45 type-II MsMADS-box genes were classified into 11 subgroups. The promoter region of MsMADS-box genes containing several hormone and stress related elements. Chromosomal location analysis revealed that 117 MsMADS-box genes were unevenly distributed on 32 chromosomes, and the remaining three genes were located on unmapped scaffolds. A total of nine pairs of segmental duplications and four groups of tandem duplications were found. Expression analysis showed that MsMADS-box genes were differentially expressed in various tissues and under abiotic stresses. qRT-PCR analysis revealed that the expression profiles of eight selected MsMADS-box genes were distinct under various stresses.

Conclusions

In this study, MsMADS-box genes were identified in the cultivated alfalfa genome based on autotetraploid level, and further confirmed by Gene Ontology (GO) analysis, phylogenetic analysis, sequence features and expression analysis. Taken together, these findings will provide clues for further study of MsMADS-box functions and alfalfa molecular breeding.

Our study is the first to systematically identify and characterize the MADS-box transcription factors in autotetraploid cultivated alfalfa (Medicago sativa L.), and eight MsMADS-box genes were significantly involved in response to various stresses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07911-9.

Genome-Wide Analysis of MADS-Box Genes in Foxtail Millet (Setaria italica L.) and Functional Assessment of the Role of SiMADS51 in the Drought Stress Response

MADS-box transcription factors play vital roles in multiple biological processes in plants. At present, a comprehensive investigation into the genome-wide identification and classification of MADS-box genes in foxtail millet (Setaria italica L.) has not been reported. In this study, we identified 72 MADS-box genes in the foxtail millet genome and give an overview of the phylogeny, chromosomal location, gene structures, and potential functions of the proteins encoded by these genes. We also found that the expression of 10 MIKC-type MADS-box genes was induced by abiotic stresses (PEG-6000 and NaCl) and exogenous hormones (ABA and GA), which suggests that these genes may play important regulatory roles in response to different stresses. Further studies showed that transgenic Arabidopsis and rice (Oryza sativa L.) plants overexpressing SiMADS51 had reduced drought stress tolerance as revealed by lower survival rates and poorer growth performance under drought stress conditions, which demonstrated that SiMADS51 is a negative regulator of drought stress tolerance in plants. Moreover, expression of some stress-related genes were down-regulated in the SiMADS51-overexpressing plants. The results of our study provide an overall picture of the MADS-box gene family in foxtail millet and establish a foundation for further research on the mechanisms of action of MADS-box proteins with respect to abiotic stresses.

A MADS-Box Gene CiMADS43 Is Involved in Citrus Flowering and Leaf Development through Interaction with CiAGL9

MADS-box genes are involved in various developmental processes including vegetative development, flower architecture, flowering, pollen formation, seed and fruit development. However, the function of most MADS-box genes and their regulation mechanism are still unclear in woody plants compared with model plants. In this study, a MADS-box gene (CiMADS43) was identified in citrus. Phylogenetic and sequence analysis showed that CiMADS43 is a GOA-like Bsister MADS-box gene. It was localized in the nucleus and as a transcriptional activator. Overexpression of CiMADS43 promoted early flowering and leaves curling in transgenic Arabidopsis. Besides, overexpression or knockout of CiMADS43 also showed leaf curl phenotype in citrus similar to that of CiMADS43 overexpressed in Arabidopsis. Protein–protein interaction found that a SEPALLATA (SEP)-like protein (CiAGL9) interacted with CiMADS43 protein. Interestingly, CiAGL9 also can bind to the CiMADS43 promoter and promote its transcription. Expression analysis also showed that these two genes were closely related to seasonal flowering and the development of the leaf in citrus. Our findings revealed the multifunctional roles of CiMADS43 in the vegetative and reproductive development of citrus. These results will facilitate our understanding of the evolution and molecular mechanisms of MADS-box genes in citrus.

Identification of transcription factors interacting with a 1274 bp promoter of MaPIP1;1 which confers high-level gene expression and drought stress Inducibility in transgenic Arabidopsis thaliana

BACKGROUND

Drought stress can severely affect plant growth and crop yield. The cloning and identification of drought-inducible promoters would be of value for genetically-based strategies to improve resistance of crops to drought.

RESULTS

Previous studies showed that the MaPIP1;1 gene encoding an aquaporin is involved in the plant drought stress response. In this study, the promoter pMaPIP1;1, which lies 1362 bp upstream of the MaPIP1;1 transcriptional initiation site, was isolated from the banana genome..And the transcription start site(A) is 47 bp before the ATG. To functionally validate the promoter, various lengths of pMaPIP1;1 were deleted and fused to GUS to generate pMaPIP1;1::GUS fusion constructs that were then transformed into Arabidopsis to generate four transformants termed M-P1, M-P2, M-P3 and M-P4.Mannitol treatment was used to simulate drought conditions. All four transformants reacted well to mannitol treatment. M-P2 (- 1274 bp to - 1) showed the highest transcriptional activity among all transgenic Arabidopsis tissues, indicating that M-P2 was the core region of pMaPIP1;1. This region of the promoter also confers high levels of gene expression in response to mannitol treatment. Using M-P2 as a yeast one-hybrid bait, 23 different transcription factors or genes that interacted with MaPIP1;1 were screened. In an dual luciferase assay for complementarity verification, the transcription factor MADS3 positively regulated MaPIP1;1 transcription when combined with the banana promoter. qRT-PCR showed that MADS3 expression was similar in banana leaves and roots under drought stress. In banana plants grown in 45% soil moisture to mimic drought stress, MaPIP1;1 expression was maximized, which further demonstrated that the MADS3 transcription factor can synergize with MaPIP1;1.

CONCLUSIONS

Together our results revealed that MaPIP1;1 mediates molecular mechanisms associated with drought responses in banana, and will expand our understanding of how AQP gene expression is regulated. The findings lay a foundation for genetic improvement of banana drought resistance.

GsMAS1 Encoding a MADS-box Transcription Factor Enhances the Tolerance to Aluminum Stress in Arabidopsis thaliana

The MADS-box transcription factors (TFs) are essential in regulating plant growth and development, and conferring abiotic and metal stress resistance. This study aims to investigate GsMAS1 function in conferring tolerance to aluminum stress in Arabidopsis. The GsMAS1 from the wild soybean BW69 line encodes a MADS-box transcription factor in Glycine soja by bioinformatics analysis. The putative GsMAS1 protein was localized in the nucleus. The GsMAS1 gene was rich in soybean roots presenting a constitutive expression pattern and induced by aluminum stress with a concentration-time specific pattern. The analysis of phenotypic observation demonstrated that overexpression of GsMAS1 enhanced the tolerance of Arabidopsis plants to aluminum (Al) stress with larger values of relative root length and higher proline accumulation compared to those of wild type at the AlCl3 treatments. The genes and/or pathways regulated by GsMAS1 were further investigated under Al stress by qRT-PCR. The results indicated that six genes resistant to Al stress were upregulated, whereas AtALMT1 and STOP2 were significantly activated by Al stress and GsMAS1 overexpression. After treatment of 50 μM AlCl3, the RNA abundance of AtALMT1 and STOP2 went up to 17-fold and 37-fold than those in wild type, respectively. Whereas the RNA transcripts of AtALMT1 and STOP2 were much higher than those in wild type with over 82% and 67% of relative expression in GsMAS1 transgenic plants, respectively. In short, the results suggest that GsMAS1 may increase resistance to Al toxicity through certain pathways related to Al stress in Arabidopsis.

Genome-Wide Identification and Analysis of the Cytochrome B5 Protein Family in Chinese Cabbage (Brassica rapa L. ssp. Pekinensis)

Cytochrome B5 (CB5) family proteins play an important role in various oxidation/reduction reactions in cells as the electron donor and are involved in a variety of biotic and abiotic stress processes. However, the function of the CB5s in Brassica rapa is still unclear. In this study, we carried out genome-wide identification, characterization, and expression analysis of BrCB5s in different tissues under adversities and stresses. It was identified that fifteen BrCB5s were distributed on different chromosomes, which were classified into seven groups (A-G) according to its phylogenetic relationship. Phylogenetic analysis of the CB5 protein sequences from six species showed that the BrCB5s conduct a close evolutionary process with the CB5s of Arabidopsis thaliana and far from those of Oryza sativa. Protein interaction analysis showed that 40 interaction patterns were predicted including two Sucrose Transporter 4 subfamily proteins (SUT 4) and Fatty Acid Hydroxylase 2 protein (FAH 2) can interact with most members of BrCB5s. The expression profile analysis indicated that BrCB5s were differentially expressed in different tissues, and the transcript abundances were significantly different under various abiotic stresses and plant hormone treatments. Our study provides a basis for a better understanding of the characteristics and biological functions of the CB5 family genes in Chinese cabbage during plant development, especially in plant responses to multiple stresses.

Genome-wide analysis of the H3K27me3 epigenome and transcriptome in Brassica rapa

BACKGROUND

Genome-wide maps of histone modifications have been obtained for several plant species. However, most studies focus on model systems and do not enforce FAIR data management principles. Here we study the H3K27me3 epigenome and associated transcriptome of Brassica rapa, an important vegetable cultivated worldwide.

FINDINGS

We performed H3K27me3 chromatin immunoprecipitation followed by high-throughput sequencing and transcriptomic analysis by 3'-end RNA sequencing from B. rapa leaves and inflorescences. To analyze these data we developed a Reproducible Epigenomic Analysis pipeline using Galaxy and Jupyter, packaged into Docker images to facilitate transparency and reuse. We found that H3K27me3 covers roughly one-third of all B. rapa protein-coding genes and its presence correlates with low transcript levels. The comparative analysis between leaves and inflorescences suggested that the expression of various floral regulatory genes during development depends on H3K27me3. To demonstrate the importance of H3K27me3 for B. rapa development, we characterized a mutant line deficient in the H3K27 methyltransferase activity. We found that braA.clf mutant plants presented pleiotropic alterations, e.g., curly leaves due to increased expression and reduced H3K27me3 levels at AGAMOUS-like loci.

CONCLUSIONS

We characterized the epigenetic mark H3K27me3 at genome-wide levels and provide genetic evidence for its relevance in B. rapa development. Our work reveals the epigenomic landscape of H3K27me3 in B. rapa and provides novel genomics datasets and bioinformatics analytical resources. We anticipate that this work will lead the way to further epigenomic studies in the complex genome of Brassica crops.

Genome wide survey, evolution and expression analysis of PHD finger genes reveal their diverse roles during the development and abiotic stress responses in Brassica rapa L

BACKGROUND

Plant homeodomain (PHD) finger proteins are widely present in all eukaryotes and play important roles in chromatin remodeling and transcriptional regulation. The PHD finger can specifically bind a number of histone modifications as an "epigenome reader", and mediate the activation or repression of underlying genes. Many PHD finger genes have been characterized in animals, but only few studies were conducted on plant PHD finger genes to this day. Brassica rapa (AA, 2n = 20) is an economically important vegetal, oilseed and fodder crop, and also a good model crop for functional and evolutionary studies of important gene families among Brassica species due to its close relationship to Arabidopsis thaliana.

RESULTS

We identified a total of 145 putative PHD finger proteins containing 233 PHD domains from the current version of B. rapa genome database. Gene ontology analysis showed that 67.7% of them were predicted to be located in nucleus, and 91.3% were predicted to be involved in protein binding activity. Phylogenetic, gene structure, and additional domain analyses clustered them into different groups and subgroups, reflecting their diverse functional roles during plant growth and development. Chromosomal location analysis showed that they were unevenly distributed on the 10 B. rapa chromosomes. Expression analysis from RNA-Seq data showed that 55.7% of them were constitutively expressed in all the tested tissues or organs with relatively higher expression levels reflecting their important housekeeping roles in plant growth and development, while several other members were identified as preferentially expressed in specific tissues or organs. Expression analysis of a subset of 18 B. rapa PHD finger genes under drought and salt stresses showed that all these tested members were responsive to the two abiotic stress treatments.

CONCLUSIONS

Our results reveal that the PHD finger genes play diverse roles in plant growth and development, and can serve as a source of candidate genes for genetic engineering and improvement of Brassica crops against abiotic stresses. This study provides valuable information and lays the foundation for further functional determination of PHD finger genes across the Brassica species.

Candidate genes expression profiling during wilting in chickpea caused by Fusarium oxysporum f. sp. ciceris race 5

Chickpea production may be seriously threatened by Fusarium wilt, a disease caused by the soil-borne fungus Fusarium oxysporum f. sp. ciceris. F. oxysporum race 5 is the most important race in the Mediterranean basin. Recently, the region responsible for resistance race 5 has been delimited within a region on chromosome 2 that spans 820 kb. To gain a better understanding of this genomic region, we used a transcriptomic approach based on quantitative real-time PCR to analyze the expression profiles of 22 selected candidate genes. We used a pair of near-isogenic lines (NILs) differing in their sensitivity to Fusarium race 5 (resistant vs susceptible) to monitor the transcriptional changes over a time-course experiment (24, 48, and 72 hours post inoculation, hpi). Qualitative differences occurred during the timing of regulation. A cluster of 12 genes were induced by the resistant NIL at 24 hpi, whereas a second cluster contained 9 genes induced by the susceptible NIL at 48 hpi. Their possible functions in the molecular defence of chickpea is discussed. Our study provides new insight into the molecular defence against Fusarium race 5 and demonstrates that development of NILs is a rich resource to facilitate the detection of candidate genes. The new genes regulated here may be useful against other Fusarium races.

Genome-Wide Identification, Gene Structure and Expression Analysis of the MADS-Box Gene Family Indicate Their Function in the Development of Tobacco (Nicotiana tabacum L.)

MADS-box genes play a pivotal role in various processes, including floral and seed development, controlling flowering time, regulation of fruits ripening, and respond to abiotic and biotic stressors in planta. Tobacco (Nicotiana tabacum) has been widely used as a model plant for analyzing the gene function, however, there has been less information on the regulation of flowering, and the associated genes. In the present study, a total of 168 NtMADS-box genes were identified from tobacco, and their phylogenetic relationship, chromosome locations, and gene structures were further analyzed. NtMADS-box genes can be clustered into four sub-families of Mα, Mγ, MIKC*, and MIKCC. A total of 111 NtMADS-box genes were distributed on 20 chromosomes, and 57 NtMADS-box genes were located on the unanchored scaffolds due to the complex and incomplete assembly of the tobacco genome. Expression profiles of NtMADS-box genes by microarray from 23 different tissues indicated that members in different NtMADS-box gene subfamilies might play specific roles in the growth and flower development, and the transcript levels of 24 NtMADS-box genes were confirmed by quantitative real-time PCR. Importantly, overexpressed NtSOC1/NtMADS133 could promote early flowering and dwarfism in transgenic tobacco plants. Therefore, our findings provide insights on the characterization of NtMADS-box genes to further study their functions in plant development.

Identification and evolutionary characterization of salt-responsive transcription factors in the succulent halophyte Suaeda fruticosa

Transcription factors are key regulatory elements that affect gene expression in response to specific signals, including environmental stresses such as salinity. Halophytes are specialized plants that have the ability to complete their life cycle in saline environments. In this study we have identified and characterized the evolutionary relationships of putative transcription factors (TF) in an obligate succulent halophyte, Suaeda fruticosa, that are involved in conferring salt tolerance. Using RNA-seq data we have analyzed the expression patterns of certain TF families, predicted protein-protein interactions, and analyzed evolutionary trajectories to elucidate their possible roles in salt tolerance. We have detected the top differentially expressed (DE) transcription factor families (MYB, CAMTA, MADS-box and bZIP) that show the most pronounced response to salinity. The majority of DE genes in the four aforementioned TF families cluster together on TF phylogenetic trees, which suggests common evolutionary origins and trajectories. This research represents the first comprehensive TF study of a leaf succulent halophyte including their evolutionary relationships with TFs in other halophyte and salt-senstive plants. These findings provide a foundation for understanding the function of salt-responsive transcription factors in salt tolerance and associated gene regulation in plants.

Na+ compartmentation strategy of Chinese cabbage in response to salt stress

Na⁺/H⁺ antiporter (NHX), responsible for counter-transport of Na⁺ and H⁺ across membranes (Na⁺ compartmentalization), plays a central role in plant salt-tolerance. In order to explore the Na⁺ compartmentalization modes and salt tolerance strategy in Chinese cabbage (Brassica rapa L. ssp. pekinensis), the seedlings of a salt-susceptible cabbage cultivar (Kuaicai 38) and a salt-tolerant cabbage cultivar (Qingmaye) were exposed to 100-400 mM NaCl for 30 days. Both of these cultivars showed a gradual decrease in fresh weight and water content and an increase in root-shoot ratio with the increasing NaCl-treatment concentration. The distribution of Na⁺ in these two cultivars was similar, with the green leaves showing the highest Na⁺ content, followed by inflated midribs, stems, and roots. The Na⁺ concentration in the apoplast was higher than that in the protoplast of the leaves. The expression levels of BrNHX1-1 and BrNHX1-2 in the leaves of Qingmaye were the highest among all BrNHX members, and increased after salt treatment. However, only BrNHX1-1 was expressed in Kuaicai 38. These results indicate that Na⁺ compartmentation into vacuoles is the major salt-adaptation strategy in Chinese cabbage. Coordinated overexpression of BrNHX1-1 and BrNHX1-2 may confer greater salt-tolerance for Chinese cabbage.

Genome-wide analysis of spatiotemporal gene expression patterns during floral organ development in Brassica rapa

Flowering is a key agronomic trait that directly influences crop yield and quality and serves as a model system for elucidating the molecular basis that controls successful reproduction, adaptation, and diversification of flowering plants. Adequate knowledge of continuous series of expression data from the floral transition to maturation is lacking in Brassica rapa. To unravel the genome expression associated with the development of early small floral buds (< 2 mm; FB2), early large floral buds (2-4 mm; FB4), stamens (STs) and carpels (CPs), transcriptome profiling was carried out with a Br300K oligo microarray. The results showed that at least 6848 known nonredundant genes (30% of the genes of the Br300K) were differentially expressed during the floral transition from vegetative tissues to maturation. Functional annotation of the differentially expressed genes (DEGs) (fold change ≥ 5) by comparison with a close relative, Arabidopsis thaliana, revealed 6552 unigenes (4579 upregulated; 1973 downregulated), including 131 Brassica-specific and 116 functionally known floral Arabidopsis homologs. Additionally, 1723, 236 and 232 DEGs were preferentially expressed in the tissues of STs, FB2, and CPs. These DEGs also included 43 transcription factors, mainly AP2/ERF-ERF, NAC, MADS-MIKC, C2H2, bHLH, and WRKY members. The differential gene expression during flower development induced dramatic changes in activities related to metabolic processes (23.7%), cellular (22.7%) processes, responses to the stimuli (7.5%) and reproduction (1%). A relatively large number of DEGs were observed in STs and were overrepresented by photosynthesis-related activities. Subsequent analysis via semiquantitative RT-PCR, histological analysis performed with in situ hybridization of BrLTP1 and transgenic reporter lines (BrLTP promoter::GUS) of B. rapa ssp. pekinensis supported the spatiotemporal expression patterns. Together, these results suggest that a temporally and spatially regulated process of the selective expression of distinct fractions of the same genome leads to the development of floral organs. Interestingly, most of the differentially expressed floral transcripts were located on chromosomes 3 and 9. This study generated a genome expression atlas of the early floral transition to maturation that represented the flowering regulatory elements of Brassica rapa.

Genome-Wide Analysis of the miRNA-mRNAs Network Involved in Cold Tolerance in Populus simonii × P. nigra

Background: Cold tolerance is important for plants’ geographical distribution and survival in extreme seasonal variations of climate. However, Populus simonii × P. nigra shows wide adaptability and strong cold resistance. Transcriptional and post-transcriptional regulation of cold-responsive genes is crucial for cold tolerance in plants. To understand the roles of regulatory RNAs under cold induction in Populus simonii × P. nigra, we constructed cDNA and small RNA libraries from leaf buds treated or not with −4 °C for 8 h for analysis. Results: Through high-throughput sequencing and differential expression analysis, 61 miRNAs and 1229 DEGs were identified under cold induction condition in Populus simonii × P. nigra. The result showed that miR167a, miR1450, miR319a, miR395b, miR393a-5p, miR408-5p, and miR168a-5p were downregulated, whereas transcription level of miR172a increased under the cold treatment. Thirty-one phased-siRNA were also obtained (reads ≥ 4) and some of them proceeded from TAS3 loci. Analysis of the differentially expressed genes (DEGs) showed that transcription factor genes such as Cluster-15451.2 (putative MYB), Cluster-16493.29872 (putative bZIP), Cluster-16493.29175 (putative SBP), and Cluster-1378.1 (putative ARF) were differentially expressed in cold treated and untreated plantlets of Populus simonii × P. nigra. Integrated analysis of miRNAs and transcriptome showed miR319, miR159, miR167, miR395, miR390, and miR172 and their target genes, including MYB, SBP, bZIP, ARF, LHW, and ATL, were predicted to be involved in ARF pathway, SPL pathway, DnaJ related photosystem II, and LRR receptor kinase, and many of them are known to resist chilling injury. Particularly, a sophisticated regulatory model including miRNAs, phasiRNAs, and targets of them was set up. Conclusions: Integrated analysis of miRNAs and transcriptome uncovered the complicated regulation of the tolerance of cold in Populus simonii × P. nigra. MiRNAs, phasiRNAs, and gene-encoded transcription factors were characterized at a whole genome level and their expression patterns were proved to be complementary. This work lays a foundation for further research of the pathway of sRNAs and regulatory factors involved in cold tolerance.

Genome-wide identification and classification of MIKC-type MADS-box genes in Streptophyte lineages and expression analyses to reveal their role in seed germination of orchid

BACKGROUND

MADS-box genes play crucial roles in plant floral organ formation and plant reproductive development. However, there is still no information on genome-wide identification and classification of MADS-box genes in some representative plant species. A comprehensive investigation of MIKC-type genes in the orchid Dendrobium officinale is still lacking.

RESULTS

Here we conducted a genome-wide analysis of MADS-box proteins from 29 species. In total, 1689 MADS-box proteins were identified. Two types of MADS-box genes, termed type I and II, were found in land plants, but not in liverwort. The SQUA, DEF/GLO, AG and SEP subfamilies existed in all the tested flowering plants, while SQUA was absent in the gymnosperm Ginkgo biloba, and no genes of the four subfamilies were found in a charophyte, liverwort, mosses, or lycophyte. This strongly corroborates the notion that clades of floral organ identity genes led to the evolution of flower development in flowering plants. Nine subfamilies of MIKC^C genes were present in two orchids, D. officinale and Phalaenopsis equestris, while the TM8, FLC, AGL15 and AGL12 subfamilies may be lost. In addition, the four clades of floral organ identity genes in both orchids displayed a conservative and divergent expression pattern. Only three MIKC-type genes were induced by cold stress in D. officinale while 15 MIKC-type genes showed different levels of expression during seed germination.

CONCLUSIONS

MIKC-type genes were identified from streptophyte lineages, revealing new insights into their evolution and development relationships. Our results show a novel role of MIKC-type genes in seed germination and provide a useful clue for future research on seed germination in orchids.

Genome-wide characterization of ALDH Superfamily in Brassica rapa and enhancement of stress tolerance in heterologous hosts by BrALDH7B2 expression

Aldehyde dehydrogenase (ALDH) carries out oxidation of toxic aldehydes using NAD⁺/NADP⁺ as cofactors. In the present study, we performed a genome-wide identification and expression analysis of genes in the ALDH gene family in Brassica rapa. A total of 23 ALDH genes in the superfamily have been identified according to the classification of ALDH Gene Nomenclature Committee (AGNC). They were distributed unevenly across all 10 chromosomes. All the 23 Brassica rapa ALDH (BrALDH) genes exhibited varied expression patterns during treatments with abiotic stress inducers and hormonal treatments. The relative expression profiles of ALDH genes in B. rapa showed that they are predominantly expressed in leaves and stem suggesting their function in the vegetative tissues. BrALDH7B2 showed a strong response to abiotic stress and hormonal treatments as compared to other ALDH genes; therefore, it was overexpressed in heterologous hosts, E. coli and yeast to study its possible function under abiotic stress conditions. Over-expression of BrALDH7B2 in heterologous systems, E. coli and yeast cells conferred significant tolerance to abiotic stress treatments. Results from this work demonstrate that BrALDH genes are a promising and untapped genetic resource for crop improvement and could be deployed further in the development of drought and salinity tolerance in B. rapa and other economically important crops.

Molecular markers based on sequence variation in BoFLC1.C9 for characterizing early- and late-flowering cabbage genotypes

BACKGROUND

Cabbage (Brassica oleracea var. capitata) is popular worldwide for consumption as a leafy vegetable. Premature flowering is triggered by low temperature, and deteriorates quality of cabbage as vegetable. In general, growers prefer late-flowering varieties to assure good quality compact head. Here, we report BoFLC1.C9 as a gene with clear sequence variation between cabbage lines with different flowering times, and proposed as molecular marker to characterize early- and late-flowering cabbage lines.

RESULTS

We identified sequence variation of 67 bp insertions in intron 2, which were contributed in flowering time variation between two inbred lines through rapid down-regulation of the BoFLC1.C9 gene in early-flowering line compared to late-flowering one upon vernalization. One set of primer 'F7R7' proposed as marker, of which was explained with 83 and 80% of flowering time variation in 141 F₂ individuals and 20 commercial lines, respectively.

CONCLUSIONS

This F7R7 marker could be used as genetic tools to characterize flowering time variation and to select as well to develop early- and late-flowering cabbage cultivars.

Asymmetric birth and death of type I and type II MADS-box gene subfamilies in the rubber tree facilitating laticifer development

The rubber tree (Hevea brasiliensis Muell. Arg.) is a rubber producing crop and contains specialized laticifers. MADS-box genes are a family of transcription factor genes that regulate plant development, especially floral organ and gametophyte development. 97 MADS-box genes were identified in the rubber tree through transcriptomes and genome mining. 93.8% of the genes were mapped onto the genome scaffolds in correspondence to the coverage (93.8%) of current version of sequenced genome. Phylogenetic analysis indicates that type II MADS-box genes have been more actively duplicated than their orthologous genes in Arabidopsis and rice, so that most (70, 72.2%) of the MADS-box genes in the rubber tree belong to type II subfamily. This is a high percentage compared to those in Arabidopsis (43.7%) and rice (56.8%). Moreover, 69 out of 70 type II genes in the rubber tree are transcribed, and they are mostly predominantly expressed in flowers, but some genes are predominantly expressed in laticifers, suggesting their roles in both flower and laticifer development. The number of type I genes in the rubber tree is only 27 (27.8%), a much smaller number compared to their orthologous genes in Arabidopsis (56.3%) and rice (43.2%). At the same time, most of the type I genes (55.6%, 15) in the rubber tree are silent and are probably pseudogenes. The high birth rate and low death rate of type II genes and low birth rate and high death rate of type I genes may corresponds to special developmental requirements in the rubber tree, e.g. the development of laticifer system for biosynthesis of cis-polyisoprene, the rubber. Moreover, atypical MIKC* factors (e.g. HbMADS1 in S-clade, and HbMADS20 in P-clade) are identified. These genes are diverged to typical MIKC* genes in sequences and facilitate functions required in laticifer development and rubber biosynthesis, which is not necessary in Arabidopsis and rice.

Genome wide analysis of MADS-box gene family in Brassica oleracea reveals conservation and variation in flower development

BACKGROUND

MADS-box genes play important roles in vegetative growth and reproductive development and are essential for the correct development of plants (particularly inflorescences, flowers, and fruits). However, this gene family has not been identified nor their functions analyzed in Brassica oleracea.

RESULTS

In this study, we performed a whole-genome survey of the complete set of MADS-box genes in B. oleracea. In total, 91 MADS-box transcription factors (TFs) were identified and categorized as type I (Mα, Mβ, Mγ) and type II (MIKCC, MIKC*) groups according to the phylogeny and gene structure analysis. Among these genes, 59 were randomly distributed on 9 chromosomes, while the other 23 were assigned to 19 scaffolds and 9 genes from NCBI had no location information. Both RNA-sequencing and quantitative real-time-PCR analysis suggested that MIKC genes had more active and complex expression patterns than M type genes and most type II genes showed high flowering-related expression profiles. Additional quantitative real-time-PCR analysis of pedicel and four flower whorls revealed that the structure of the B.oleracea MIKC genes was conserved, but their homologues showed variable expression patterns compared to those in Arabidopsis thaliana.

CONCLUSION

This paper gives a detailed overview of the BolMADS genes and their expression patterns. The results obtained in this study provide useful information for understanding the molecular regulation of flower development and further functional characterization of MADS-box genes in B. oleracea.

CaMADS, a MADS-box transcription factor from pepper, plays an important role in the response to cold, salt, and osmotic stress

MADS-box family transcription factors play essential roles in the growth and development of plants, and some MADS-box genes have also been found to participate in the stress-response. At present, little information regarding stress-related MADS-box genes is available in pepper. We isolated a MADS-box transcription factor gene from Capsicum annuum, which we named CaMADS. CaMADS expression is induced by low and high temperature, salt, and osmotic stress, and also by abscisic acid (ABA), salicylic acid (SA), methyl-jasmonic acid (MeJA), and CaCl_2. To understand the function of CaMADS in the abiotic stress response, we generated pepper plants in which CaMADS expression was down-regulated using VIGS (Virus-induced Gene Silencing), and also transgenic Arabidopsis plants overexpressing CaMADS. We found that CaMADS-down-regulated seedlings were more seriously injured than WT after cold, NaCl, and mannitol treatment, and showed increased electrolyte leakage, malondialdehyde (MDA) levels, and lower chlorophyll content. CaMADS-overexpressing Arabidopsis plants were more tolerant to these stresses than WT, and showed significantly high survival rates and lower H₂O₂ and super oxide radical contents after cold treatment. CaMADS-overexpressing plants had higher germination rates and percentages of green cotyledons following NaCl and mannitol treatment. Root lengths and fresh weight in CaMADS-overexpressing plants were also significantly longer and higher, respectively, than in WT plants. Taken together, our results suggest that CaMADS functions as a positive stress-responsive transcription factor in the cold, salt, and osmotic stress signaling pathways.

Genome-wide transcriptome analysis reveals the molecular mechanism of high temperature-induced floral abortion in Litchi chinensis

Background

Warm winter and hot spring attributed to global warming affected floral development and may induce floral abortion, resulted in poor flowering in litchi. To identify genes potentially involved in litchi floral abortion, six RNA-sequencing (RNA-Seq) libraries of the developing panicles (DPs) under low temperature (LT) conditions and the shrinking panicles (SPs) under high temperature (HT) conditions were constructed.

Results

3.07–8.97 × 10⁶ clean reads were generated. Digital expression of the DPs with that of the SPs was compared. As a result, 1320 up-regulated and 981 down-regulated differentially expressed genes (DEGs) were identified. From the enriched GO-term, 54 temperature responsive DEGs, 23 hormone homeostasis- or biosynthesis-related DEGs, 137 hormone signal transduction or responsive DEGs, and 18 flowering-related DEGs were identified. Partial Least Squares Structural Equation Modeling (PLS-SEM) analysis indicated that the effects of hormone-related DEGs on NACs, MYBs, WRKYs were stronger than that on flowering-related DEGs. Expression pattern analysis of the inflorescence in ‘Nuomici’ and ‘Huaizhi’ under LT and HT conditions showed that genes homologous to AIL6 (LcAIL6), LHY (LcLHY), MED16 (LcMED16), SKIP20 (LcSKIP20), POD20 (LcPOD20) in the two cultivars had similar expression trends.

Conclusion

This study elucidated the transcriptome in the HT-induced floral abortion and identified key genes involved in the process. NACs, MYBs, WRKYs may act as central players involved in the HT-induced floral abortion underlying hormonal control. Increased transcript in LcLHY, LcMED16, LcSKIP20, LcPOD20 and decreased transcript in LcAIL6 might be related to the inhibition of floral development. Our studies provided potential genes for the future molecular breeding of new cultivars that can reduce floral abortion under warm climates, and a novel clue to reveal the relationship of biological processes based on the RNA-seq data using PLS-SEM.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5493-8) contains supplementary material, which is available to authorized users.

Identification of Genes Involved in Low Temperature-Induced Senescence of Panicle Leaf in Litchi chinensis

Warm winters and hot springs may promote panicle leaf growing and repress floral development. To identify genes potentially involved in litchi panicle leaf senescence, eight RNA-sequencing (RNA-Seq) libraries of the senescing panicle leaves under low temperature (LT) conditions and the developing panicle leaves under high temperature (HT) conditions were constructed. For each library, 4.78–8.99 × 10⁶ clean reads were generated. Digital expression of the genes was compared between the senescing and developing panicle leaves. A total of 6477 upregulated differentially expressed genes (DEGs) (from developing leaves to senescing leaves), and 6318 downregulated DEGs were identified, 158 abscisic acid (ABA)-, 68 ethylene-, 107 indole-3-acetic acid (IAA)-, 27 gibberellic acid (GA)-, 68 cytokinin (CTK)-, 37 salicylic acid (SA)-, and 23 brassinolide (BR)-related DEGs. Confirmation of the RNA-Seq data by quantitative real-time PCR (qRT-PCR) analysis suggested that expression trends of the 10 candidate genes using qRT-PCR were similar to those revealed by RNA-Seq, and a significantly positive correlation between the obtained data from qRT-PCR and RNA-Seq were found, indicating the reliability of our RNA-Seq data. The present studies provided potential genes for the future molecular breeding of new cultivars that can induce panicle leaf senescence and reduce floral abortion under warm climates.

Identification of new regulators through transcriptome analysis that regulate anthocyanin biosynthesis in apple leaves at low temperatures

Anthocyanin pigments play many roles in plants, including providing protection against biotic and abiotic stresses. To identify new regulatory genes in apple (Malus domestica) that may be involved in regulating low temperature induced anthocyanin biosynthesis, we performed RNA-seq analysis of leaves from the ‘Gala’ apple cultivar following exposure to a low temperature (16 °C). A visible red color appeared on the upper leaves and the anthocyanin content increased significantly after the low temperature treatment. Genes from the flavonoid biosynthesis pathway were significantly enriched among the differentially expressed genes, and the expression of several transcription factors was shown by WGCNA (weighted gene co-expression network analysis) to correlate with anthocyanin accumulation, including members of the MYB, MADS, WRKY, WD40, Zinc Finger and HB-ZIP families. Three MYB transcription factors (MdMYB12, MdMYB22 and MdMYB114), which had several CBF/DREB response elements in their promoters, were significantly induced by low temperature exposure and their expression also correlated highly with anthocyanin accumulation. We hypothesize that they may act as regulators of anthocyanin biosynthesis and be regulated by CBF/DREB transcription factors in apple leaves under low temperature conditions. The analyses presented here provide insights into the molecular mechanisms underlying anthocyanin accumulation during low temperature exposure.

miR824/AGAMOUS-LIKE16 Module Integrates Recurring Environmental Heat Stress Changes to Fine-Tune Poststress Development

Plant development is continually fine-tuned based on environmental factors. How environmental perturbations are integrated into the developmental programs and how poststress adaptation is regulated remains an important topic to dissect. Vegetative to reproductive phase change is a very important developmental transition that is complexly regulated based on endogenous and exogenous cues. Proper timing of flowering is vital for reproductive success. It has been shown previously that AGAMOUS LIKE 16 (AGL16), a MADS-box transcription factor negatively regulates flowering time transition through FLOWERING LOCUS T (FT), a central downstream floral integrator. AGL16 itself is negatively regulated by the microRNA miR824. Here we present a comprehensive molecular analysis of miR824/AGL16 module changes in response to mild and recurring heat stress. We show that miR824 accumulates gradually in response to heat due to the combination of transient transcriptional induction and posttranscriptional stability. miR824 induction requires heat shock cis-elements and activity of the HSFA1 family and HSFA2 transcription factors. Parallel to miR824 induction, its target AGL16 is decreased, implying direct causality. AGL16 posttranscriptional repression during heat stress, however, is more complex, comprising of a miRNA-independent, and a miR824-dependent pathway. We also show that AGL16 expression is leaf vein-specific and overlaps with miR824 (and FT) expression. AGL16 downregulation in response to heat leads to a mild derepression of FT. Finally, we present evidence showing that heat stress regulation of miR824/AGL16 is conserved within Brassicaceae. In conclusion, due to the enhanced post-transcriptional stability of miR824, stable repression of AGL16 is achieved following heat stress. This may serve to fine-tune FT levels and alter flowering time transition. Stress-induced miR824, therefore, can act as a "posttranscriptional memory factor" to extend the acute impact of environmental fluctuations in the poststress period.

A Large-Scale Genome-Wide Association Analyses of Ethiopian Sorghum Landrace Collection Reveal Loci Associated With Important Traits

The eastern Africa region, Ethiopia and its surroundings, is considered as the center of origin and diversity for sorghum, and has contributed to global sorghum genetic improvement. The germplasm from this region harbors enormous genetic variation for various traits but little is known regarding the genetic architecture of most traits. Here, 1425 Ethiopian landrace accessions were phenotyped under field conditions for presence or absence of awns, panicle compactness and shape, panicle exsertion, pericarp color, glume cover, plant height and smut resistance under diverse environmental conditions in Ethiopia. In addition, F1 hybrids obtained from a subset of 1341 accessions crossed to an A1 cytoplasmic male sterile line, ATx623, were scored for fertility/sterility reactions. Subsequently, genotyping-by-sequencing generated a total of 879,407 SNPs from which 72,190 robust SNP markers were selected after stringent quality control (QC). Pairwise distance-based hierarchical clustering identified 11 distinct groups. Of the genotypes assigned to either one of the 11 sub-populations, 65% had high ancestry membership coefficient with the likelihood of more than 0.60 and the remaining 35% represented highly admixed accessions. A genome-wide association study (GWAS) identified loci and SNPs associated with aforementioned traits. GWAS based on compressed mixed linear model (CMLM) identified SNPs with significant association (FDR ≤ 0.05) to the different traits studied. The percentage of total phenotypic variation explained with significant SNPs across traits ranged from 2 to 43%. Candidate genes showing significant association with different traits were identified. The sorghum bHLH transcription factor, ABORTED MICROSPORES was identified as a strong candidate gene conditioning male fertility. Notably, sorghum CLAVATA1 receptor like kinase, known for regulation of plant growth, and the ETHYLENE RESPONSIVE TRANSCRIPTION FACTOR gene RAP2-7, known to suppress transition to flowering, were significantly associated with plant height. In addition, the YELLOW SEED1 like MYB transcription factor and TANNIN1 showed strong association with pericarp color validating previous observations. Overall, the genetic architecture of natural variation representing the complex Ethiopian sorghum germplasm was established. The study contributes to the characterization of genes and alleles controlling agronomic traits, and will serve as a source of markers for molecular breeding.