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.)

Evolutionary analysis of MADS-box genes in buckwheat species and functional study of FdMADS28 in flavonoid metabolism

The MADS-box gene family is a transcription factor family that is widely expressed in plants. It controls secondary metabolic processes in plants and encourages the development of tissues like roots and flowers. However, the phylogenetic analysis and evolutionary model of MADS-box genes in Fagopyrum species has not been reported yet. This study identified the MADS-box genes of three buckwheat species at the whole genome level, and conducted systematic evolution and physicochemical analysis. The results showed that these genes can be divided into four subfamilies, with fragment duplication being the main way for the gene family expansion. During the domestication process from golden buckwheat to tartary buckwheat and the common buckwheat, the Ka/Ks ratio indicated that most members of the family experienced strong purification selection pressure, and with individual gene pairs experiencing positive selection. In addition, we combined the expression profile data of the MADS genes, mGWAS data, and WGCNA data to mine genes FdMADS28/48/50 that may be related to flavonoid metabolism. The results also showed that overexpression of FdMADS28 could increase rutin content by decreasing Kaempferol pathway content in hairy roots, and increase the resistance and growth of hairy roots to PEG and NaCl. This study systematically analyzed the evolutionary relationship of MADS-box genes in the buckwheat species, and elaborated on the expression patterns of MADS genes in different tissues under biotic and abiotic stresses, laying an important theoretical foundation for further elucidating their role in flavonoid metabolism.

The Characterization of G-Quadruplexes in Tobacco Genome and Their Function under Abiotic Stress

Tobacco is an ideal model plant in scientific research. G-quadruplex is a guanine-rich DNA structure, which regulates transcription and translation. In this study, the prevalence and potential function of G-quadruplexes in tobacco were systematically analyzed. In tobacco genomes, there were 2,924,271,002 G-quadruplexes in the nuclear genome, 430,597 in the mitochondrial genome, and 155,943 in the chloroplast genome. The density of the G-quadruplex in the organelle genome was higher than that in the nuclear genome. G-quadruplexes were abundant in the transcription regulatory region of the genome, and a difference in G-quadruplex density in two DNA strands was also observed. The promoter of 60.4% genes contained at least one G-quadruplex. Compared with up-regulated differentially expressed genes (DEGs), the G-quadruplex density in down-regulated DEGs was generally higher under drought stress and salt stress. The G-quadruplex formed by simple sequence repeat (SSR) and its flanking sequence in the promoter region of the NtBBX (Nitab4.5_0002943g0010) gene might enhance the drought tolerance of tobacco. This study lays a solid foundation for further research on G-quadruplex function in tobacco and other plants.

Genome-wide identification and expression analysis of NF-Y gene family in tobacco (Nicotiana tabacum L.)

Nuclear factor Y (NF-Y) gene family is an important transcription factor composed of three subfamilies of NF-YA, NF-YB and NF-YC, which is involved in plant growth, development and stress response. In this study, 63 tobacco NF-Y genes (NtNF-Ys) were identified in Nicotiana tabacum L., including 17 NtNF-YAs, 30 NtNF-YBs and 16 NtNF-YCs. Phylogenetic analysis revealed ten pairs of orthologues from tomato and tobacco and 25 pairs of paralogues from tobacco. The gene structure of NtNF-YAs exhibited similarities, whereas the gene structure of NtNF-YBs and NtNF-YCs displayed significant differences. The NtNF-Ys of the same subfamily exhibited a consistent distribution of motifs and protein 3D structure. The protein interaction network revealed that NtNF-YC12 and NtNF-YC5 exhibited the highest connectivity. Many cis-acting elements related to light, stress and hormone response were found in the promoter of NtNF-Ys. Transcriptome analysis showed that more than half of the NtNF-Y genes were expressed in all tissues, and NtNF-YB9/B14/B15/B16/B17/B29 were specifically expressed in roots. A total of 15, 12, 5, and 6 NtNF-Y genes were found to respond to cold, drought, salt, and alkali stresses, respectively. The results of this study will lay a foundation for further study of NF-Y genes in tobacco and other Solanaceae plants.

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.

Genome-wide characterization and expression analysis of MADS-box transcription factor gene family in Perilla frutescens

MADS-box transcription factors are widely involved in the regulation of plant growth, developmental processes, and response to abiotic stresses. Perilla frutescens, a versatile plant, is not only used for food and medicine but also serves as an economical oil crop. However, the MADS-box transcription factor family in P. frutescens is still largely unexplored. In this study, a total of 93 PfMADS genes were identified in P. frutescens genome. These genes, including 37 Type I and 56 Type II members, were randomly distributed across 20 chromosomes and 2 scaffold regions. Type II PfMADS proteins were found to contain a greater number of motifs, indicating more complex structures and diverse functions. Expression analysis revealed that most PfMADS genes (more than 76 members) exhibited widely expression model in almost all tissues. The further analysis indicated that there was strong correlation between some MIKCC-type PfMADS genes and key genes involved in lipid synthesis and flavonoid metabolism, which implied that these PfMADS genes might play important regulatory role in the above two pathways. It was further verified that PfMADS47 can effectively mediate the regulation of lipid synthesis in Chlamydomonas reinhardtii transformants. Using cis-acting element analysis and qRT-PCR technology, the potential functions of six MIKCC-type PfMADS genes in response to abiotic stresses, especially cold and drought, were studied. Altogether, this study is the first genome-wide analysis of PfMADS. This result further supports functional and evolutionary studies of PfMADS gene family and serves as a benchmark for related P. frutescens breeding studies.

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 Identification and Expression of the Paulownia fortunei MADS-Box Gene Family in Response to Phytoplasma Infection

Paulownia witches’ broom (PaWB), caused by phytoplasmas, is the most devastating infectious disease of Paulownia. Although a few MADS-box transcription factors have been reported to be involved in the formation of PaWB, there has been little investigation into all of the MADS-box gene family in Paulownia. The objective of this study is to identify the MADS-box gene family in Paulownia fortunei on a genome-wide scale and explore their response to PaWB infection. Bioinformatics software were used for identification, characterization, subcellular localization, phylogenetic analysis, the prediction of conserved motifs, gene structures, cis-elements, and protein-protein interaction network construction. The tissue expression profiling of PfMADS-box genes was analyzed by quantitative real-time polymerase chain reaction (qRT-PCR). Transcriptome data and the protein interaction network prediction were combined to screen the genes associated with PaWB formation. We identified 89 MADS-box genes in the P. fortunei genome and categorized them into 14 subfamilies. The comprehensive analysis showed that segment duplication events had significant effects on the evolution of the PfMADS-box gene family; the motif distribution of proteins in the same subfamily are similar; development-related, phytohormone-responsive, and stress-related cis-elements were enriched in the promoter regions. The tissue expression pattern of PfMADS-box genes suggested that they underwent subfunctional differentiation. Three genes, PfMADS3, PfMADS57, and PfMADS87, might be related to the occurrence of PaWB. These results will provide a valuable resource to explore the potential functions of PfMADS-box genes and lay a solid foundation for understanding the roles of PfMADS-box genes in paulownia–phytoplasma interactions.

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.

MIKC-type MADS-box transcription factor gene family in peanut: Genome-wide characterization and expression analysis under abiotic stress

Peanut (Arachis hypogaea) is one of the most important economic crops around the world, especially since it provides vegetable oil and high-quality protein for humans. Proteins encoded by MADS-box transcription factors are widely involved in regulating plant growth and development as well as responses to abiotic stresses. However, the MIKC-type MADS-box TFs in peanut remains currently unclear. Hence, in this study, 166 MIKC-type MADS-box genes were identified in both cultivated and wild-type peanut genomes, which were divided into 12 subfamilies. We found a variety of development-, hormone-, and stress-related cis-acting elements in the promoter region of peanut MIKC-type MADS-box genes. The chromosomal distribution of peanut MADS-box genes was not random, and gene duplication contributed to the expansion of the MADS-box gene family. The interaction network of the peanut AhMADS proteins was established. Expression pattern analysis showed that AhMADS genes were specifically expressed in tissues and under abiotic stresses. It was further confirmed via the qRT-PCR technique that five selected AhMADS genes could be induced by abiotic and hormone treatments and presented different expressive profiles under various stresses. Taken together, these findings provide valuable information for the exploration of candidate genes in molecular breeding and further study of AhMADS gene functions.

Genome-Wide Identification and Expression Analysis of the MADS-Box Gene Family in Sweet Potato [Ipomoea batatas (L.) Lam]

MADS-box gene, one of the largest transcription factor families in plants, is a class of transcription factors widely present in eukaryotes. It plays an important role in plant growth and development and participates in the growth and development of flowers and fruits. Sweet potato is the seventh most important food crop in the world. Its tuberous roots, stems, and leaves contain a large number of proteins, lipids, carotenoids, anthocyanins, conjugated phenolic acids, and minerals, which have high edible, forage, and medicinal value, and is also an important energy crop. At present, MADS-box genes in sweet potato have rarely been reported, and there has been no study on the genome-wide identification and classification of MADS-box genes in Ipomoea batatas. This study provided the first comprehensive analysis of sweet potato MADS-box genes. We identified 95 MADS-box genes, analyzed the structure and protein of sweet potato MADS-box genes, and categorized them based on phylogenetic analysis with Arabidopsis MADS-box proteins. Chromosomal localization indicated an unequal number of MADS-box genes in all 14 chromosomes except LG3, with more than 10 MADS-box genes located on chromosomes LG7, LG11, and LG15. The MADS domain and core motifs of the sweet potato MADS-box genes were identified by motif analysis. We identified 19 MADS-box genes with collinear relationships and analyzed duplication events. Cis-acting elements, such as light-responsive, auxin-responsive, drought-inducible, and MeJA-responsive elements, were found in the promoter region of the MADS-box genes in sweet potato, which further indicates the basis of MADS-box gene regulation in response to environmental changes and hormones. RNA-seq suggested that sweet potato MADS-box genes exhibit tissue-specific expression patterns, with 34 genes highly expressed in sweet potato flowers and fruits, and 19 genes highly expressed in the tuberous root, pencil root, or fibrous root. qRT-PCR again validated the expression levels of the 10 genes and found that IbMADS1, IbMADS18, IbMADS19, IbMADS79, and IbMADS90 were highly expressed in the tuberous root or fibrous root, and IbMADS18, IbMADS31, and IbMADS83 were highly expressed in the fruit. In this study, the molecular basis of MADS-box genes of sweet potato was analyzed from various angles. The effects of MADS-box genes on the growth and development of sweet potato were investigated, which may provide a certain theoretical basis for molecular breeding of sweet potato.

Genome-wide identification and expression analysis of NtbHLH gene family in tobacco (Nicotiana tabacum) and the role of NtbHLH86 in drought adaptation

The bHLH transcription factors play pivotal roles in plant growth and development, production of secondary metabolites and responses to various environmental stresses. Although the bHLH genes have been well studied in model plant species, a comprehensive investigation of the bHLH genes is required for tobacco with newly obtained high-quality genome. In the present study, a total of 309 NtbHLH genes were identified and can be divided into 23 subfamilies. The conserved amino acids which are essential for their function were predicted for the NtbHLH proteins. Moreover, the NtbHLH genes were conserved during evolution through analyzing the gene structures and conserved motifs. A total of 265 NtbHLH genes were localized in the 24 tobacco chromosomes while the remained 44 NtbHLH genes were mapped to the scaffolds due to the complexity of tobacco genome. Moreover, transcripts of NtbHLH genes were obviously tissue-specific expressed from the gene-chip data from 23 tobacco tissues, and expressions of 20 random selected NtbHLH genes were further confirmed by quantitative real-time PCR, indicating their potential functions in the plant growth and development. Importantly, overexpressed NtbHLH86 gene confers improve drought tolerance in tobacco indicating that it might be involved in the regulation of drought stress. Therefore, our findings here provide a valuable information on the characterization of NtbHLH genes and further investigation of their functions in tobacco.

The Euscaphis japonica genome and the evolution of malvids

Malvids is one of the largest clades of rosids, includes 58 families and exhibits remarkable morphological and ecological diversity. Here, we report a high-quality chromosome-level genome assembly for Euscaphis japonica, an early-diverging species within malvids. Genome-based phylogenetic analysis suggests that the unstable phylogenetic position of E. japonica may result from incomplete lineage sorting and hybridization event during the diversification of the ancestral population of malvids. Euscaphis japonica experienced two polyploidization events: the ancient whole genome triplication event shared with most eudicots (commonly known as the γ event) and a more recent whole genome duplication event, unique to E. japonica. By resequencing 101 samples from 11 populations, we speculate that the temperature has led to the differentiation of the evergreen and deciduous of E. japonica and the completely different population histories of these two groups. In total, 1012 candidate positively selected genes in the evergreen were detected, some of which are involved in flower and fruit development. We found that reddening and dehiscence of the E. japonica pericarp and long fruit-hanging time promoted the reproduction of E. japonica populations, and revealed the expression patterns of genes related to fruit reddening, dehiscence and abscission. The key genes involved in pentacyclic triterpene synthesis in E. japonica were identified, and different expression patterns of these genes may contribute to pentacyclic triterpene diversification. Our work sheds light on the evolution of E. japonica and malvids, particularly on the diversification of E. japonica and the genetic basis for their fruit dehiscence and abscission.

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.

Bioinformatics and Expression Analysis of IDA-Like Genes Reveal Their Potential Functions in Flower Abscission and Stress Response in Tobacco (Nicotiana tabacum L.)

The inflorescence deficient in abscission-like (IDL) genes have been shown to play critical roles in floral organ abscission, lateral root formation and various stress responses in Arabidopsis. The IDL gene family has been characterized in a number of plant species, while limited information is available about IDL genes of tobacco. In the current study, 15 NtIDL members were identified in the tobacco genome, and were classified into six groups together with IDL members from other species. Evolution analysis suggested that the NtIDL members form group VI might have originated from duplication events. Notably, NtIDL06 shared high similarities with AtIDA in the EPIP sequence, and its encoding gene was highly expressed in the abscission zone of flowers at late developmental stages, implying that NtIDL06 might regulate tobacco flower abscission. In addition, the results from cis-elements analysis of promoters and expression after stress treatments suggested that NtIDL members might be involved in various stress responses of tobacco. The results from this study provide information for further functional analysis related to flower abscission and stress responses of NtIDL genes.

Genome-wide identification and expression analysis of the MADS-box transcription factor family in Camellia sinensis

The MADS-box genes are an important class of transcription factors and play critical roles in flower development. However, the functions of these genes in the economically important drinking plant, Camellia sinensis, are still not reported. Here, an evolutionary analysis of tea MADS-box genes was performed at whole genome level. A total of 83 MADS-box genes were identified in tea, and their gene structures and expression patterns were further analyzed. The tea MADS-box genes were classified into Mα (26), Mβ (12), Mγ (9), MIKC* (7), and MIKCC (29) clade according to their phylogenetic relationship with Arabidopsis thaliana. Several cis-elements were identified in the promoter regions of the CsMADS genes that are important in regulating growth, development, light responses, and the response to several stresses. Most CsMADS genes display clear different expression patterns in different organs and different species of tea plant. The expression of CsMADS genes can be regulated by abiotic stresses and phytohormone treatment. Our results lay the foundation for future research on the function of CsMADS genes and beneficial for improving tea agricultural traits in the future.

The Phoebe genome sheds light on the evolution of magnoliids

Lauraceae includes the genus Phoebe, and the family is linked to the evolution of magnoliids. We sequenced the genome of Phoebe bournei Nanmu. The assembled genome size was 989.19 Mb, with a contig N50 value of 2.05 Mb. A total of 28,198 protein-coding genes were annotated in P. bournei. Whole-genome duplication (WGD) analysis showed that Lauraceae has experienced two WGD events; the older WGD event occurred just before the divergence of Lauraceae and Magnoliales, and the more recent WGD was shared by all lineages of Lauraceae. The phylogenetic tree showed that magnoliids form a sister clade to monocots and eudicots. We also identified 63 MADS-box genes, including AGL12-like genes that may be related to the regulation of P. bournei roots and FIN219-like genes encoding GH3 proteins, which are involved in photomorphogenesis. SAUR50-like genes involved in light signal-mediated pedicel or stem development were also identified. Four ATMYB46- and three PtrEPSP-homologous genes related to lignin biosynthesis were identified. These genes may be associated with the formation of straight trunks in P. bournei. Overall, the P. bournei reference genome provides insight into the origin, evolution, and diversification of Phoebe and other magnoliids.