Genome-wide investigation of the MADS gene family and dehulling genes in tartary buckwheat (Fagopyrum tataricum)

Genome-wide identification and expression pattern analysis of MIKC-Type MADS-box genes in Chionanthus retusus, an androdioecy plant

BACKGROUND

The MADS-box gene family is widely distributed in the plant kingdom, and its members typically encoding transcription factors to regulate various aspects of plant growth and development. In particular, the MIKC-type MADS-box genes play a crucial role in the determination of floral organ development and identity recognition. As a type of androdioecy plant, Chionanthus retusus have unique gender differentiation. Manifested as male individuals with only male flowers and female individuals with only bisexual flowers. However, due to the lack of reference genome information, the characteristics of MIKC-type MADS-box genes in C. retusus and its role in gender differentiation of C. retusus remain largely unknown. Therefore, it is necessary to identify and characterize the MADS-box gene family within the genome of the C. retusus.

RESULTS

In this study, we performed a genome-wide identification and analysis of MIKC-type MADS-box genes in C. retusus (2n = 2x = 46), utilizing the latest reference genome, and studied its expression pattern in individuals of different genders. As a result, we identified a total of 61 MIKC-type MADS-box genes in C. retusus. 61 MIKC-type MADS-box genes can be divided into 12 subfamilies and distributed on 18 chromosomes. Genome collinearity analysis revealed their conservation in evolution, while gene structure, domains and motif analysis indicated their conservation in structure. Finally, based on their expression patterns in floral organs of different sexes, we have identified that CrMADS45 and CrMADS60 may potentially be involved in the gender differentiation of C. retusus.

CONCLUSIONS

Our studies have provided a general understanding of the conservation and characteristics of the MIKC-type MADS-box genes family in C. retusus. And it has been demonstrated that members of the AG subfamily, CrMADS45 and CrMADS60, may play important roles in the gender differentiation of C. retusus. This provides a reference for future breeding efforts to improve flower types in C. retusus and further investigate the role of MIKC-type MADS-box genes in gender differentiation.

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.

Genome-wide identification and characterization of TCP gene family in Dendrobium nobile and their role in perianth development

TCP is a widely distributed, essential plant transcription factor that regulates plant growth and development. An in-depth study of TCP genes in Dendrobium nobile, a crucial parent in genetic breeding and an excellent model material to explore perianth development in Dendrobium, has not been conducted. We identified 23 DnTCP genes unevenly distributed across 19 chromosomes and classified them as Class I PCF (12 members), Class II: CIN (10 members), and CYC/TB1 (1 member) based on the conserved domain and phylogenetic analysis. Most DnTCPs in the same subclade had similar gene and motif structures. Segmental duplication was the predominant duplication event for TCP genes, and no tandem duplication was observed. Seven genes in the CIN subclade had potential miR319 and -159 target sites. Cis-acting element analysis showed that most DnTCP genes contained many developmental stress-, light-, and phytohormone-responsive elements in their promoter regions. Distinct expression patterns were observed among the 23 DnTCP genes, suggesting that these genes have diverse regulatory roles at different stages of perianth development or in different organs. For instance, DnTCP4 and DnTCP18 play a role in early perianth development, and DnTCP5 and DnTCP10 are significantly expressed during late perianth development. DnTCP17, 20, 21, and 22 are the most likely to be involved in perianth and leaf development. DnTCP11 was significantly expressed in the gynandrium. Specially, MADS-specific binding sites were present in most DnTCP genes putative promoters, and two Class I DnTCPs were in the nucleus and interacted with each other or with the MADS-box. The interactions between TCP and the MADS-box have been described for the first time in orchids, which broadens our understanding of the regulatory network of TCP involved in perianth development in orchids.

Genome-Wide Analysis of the MADS-box Gene Family and Expression Analysis during Anther Development in Salvia miltiorrhiza

MADS-box genes constitute a large family of transcription factors that play important roles in plant growth and development. However, our understanding of MADS-box genes involved in anther development and male sterility in Salvia miltiorrhiza is still limited. In this study, 63 MADS-box genes were identified from the genome of the male sterility ecotype Sichuan S. miltiorrhiza (S. miltiorrhiza_SC) unevenly distributed among eight chromosomes. Phylogenetic analysis classified them into two types and 17 subfamilies. They contained 1 to 12 exons and 10 conserved motifs. Evolution analysis showed that segmental duplication was the main force for the expansion of the SmMADS gene family, and duplication gene pairs were under purifying selection. Cis-acting elements analysis demonstrated that the promoter of SmMADS genes contain numerous elements associated with plant growth and development, plant hormones, and stress response. RNA-seq showed that the expression levels of B-class and C-class SmMADS genes were highly expressed during anther development, with SmMADS11 likely playing an important role in regulating anther development and male fertility in S. miltiorrhiza_SC. Overall, this study provides a comprehensive analysis of the MADS-box gene family in S. miltiorrhiza, shedding light on their potential role in anther development and male sterility.

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.

The Tartary buckwheat bHLH gene ALCATRAZ enables Arabidopsis thaliana silique dehiscence

Tartary buckwheat is popular because of its rich nutrients. However, the difficulty in shelling restricts food production. The gene ALCATRAZ (AtALC) plays a key role in silique dehiscence in Arabidopsis thaliana. In this study, an atalc mutant was obtained by CRISPR/Cas9, and a FtALC gene homologous to AtALC was complemented into the atalc mutant to verify its function. Phenotypic observations showed that three atalc mutant lines did not dehiscence, while ComFtALC lines recovered the dehiscence phenotype. The contents of lignin, cellulose, hemicellulose, and pectin in the siliques of all the atalc mutant lines were significantly higher than those in the wild-type and ComFtALC lines. Moreover, FtALC was found to regulate the expression of cell wall pathway genes. Finally, the interaction of FtALC with FtSHP and FtIND was verified by yeast two-hybrid, bimolecular fluorescent complimentary (BIFC) and firefly luciferase completion imaging assays (LCIs). Our findings enrich the silique regulatory network and lay the foundation for the cultivation of easily shelled tartary buckwheat varieties.

Genome-Wide Identification of the MADS-Box Gene Family during Male and Female Flower Development in Chayote (Sechium edule)

The MADS-box gene plays an important role in plant growth and development. As an important vegetable of Cucurbitaceae, chayote has great edible and medicinal value. So far, there is little molecular research on chayote, and there are no reports on the MADS-box transcription factor of chayote. In this study, the MADS-box gene family of chayote was analyzed for the first time, and a total of 70 MADS-box genes were identified, including 14 type I and 56 type II MICK MADS genes. They were randomly distributed on 13 chromosomes except for chromosome 11. The light response element, hormone response element and abiotic stress response element were found in the promoter region of 70 MADS genes, indicating that the MADS gene can regulate the growth and development of chayote, resist abiotic stress, and participate in hormone response; GO and KEGG enrichment analysis also found that SeMADS genes were mainly enriched in biological regulation and signal regulation, which further proved the important role of MADS-box gene in plant growth and development. The results of collinearity showed that segmental duplication was the main driving force of MADS gene expansion in chayote. RNA-seq showed that the expression levels of SeMADS06, SeMADS13, SeMADS26, SeMADS28, SeMADS36 and SeMADS37 gradually increased with the growth of chayote, indicating that these genes may be related to the development of root tubers of 'Tuershao'. The gene expression patterns showed that 12 SeMADS genes were specifically expressed in the male flower in 'Tuershao' and chayote. In addition, SeMADS03 and SeMADS52 may be involved in regulating the maturation of male flowers of 'Tuershao' and chayote. SeMADS21 may be the crucial gene in the development stage of the female flower of 'Tuershao'. This study laid a theoretical foundation for the further study of the function of the MADS gene in chayote in the future.

QTL mapping and candidate gene analysis for yield and grain weight/size in Tartary buckwheat

BACKGROUND

Grain weight/size influences not only grain yield (GY) but also nutritional and appearance quality and consumer preference in Tartary buckwheat. The identification of quantitative trait loci (QTLs)/genes for grain weight/size is an important objective of Tartary buckwheat genetic research and breeding programs.

RESULTS

Herein, we mapped the QTLs for GY, 1000-grain weight (TGW), grain length (GL), grain width (GW) and grain length-width ratio (L/W) in four environments using 221 recombinant inbred lines (XJ-RILs) derived from a cross of 'Xiaomiqiao × Jinqiaomai 2'. In total, 32 QTLs, including 7 for GY, 5 for TGW, 6 for GL, 11 for GW and 3 for L/W, were detected and distributed in 24 genomic regions. Two QTL clusters, qClu-1-3 and qClu-1-5, located on chromosome Ft1, were revealed to harbour 7 stable major QTLs for GY (qGY1.2), TGW (qTGW1.2), GL (qGL1.1 and qGL1.4), GW (qGW1.7 and qGW1.10) and L/W (qL/W1.2) repeatedly detected in three and above environments. A total of 59 homologues of 27 known plant grain weight/size genes were found within the physical intervals of qClu-1-3 and qClu-1-5. Six homologues, FtBRI1, FtAGB1, FtTGW6, FtMADS1, FtMKK4 and FtANT, were identified with both non-synonymous SNP/InDel variations and significantly differential expression levels between the two parents, which may play important roles in Tatary buckwheat grain weight/size control and were chosen as core candidate genes for further investigation.

CONCLUSIONS

Two stable major QTL clusters related to grain weight/size and six potential key candidate genes were identified by homology comparison, SNP/InDel variations and qRT‒qPCR analysis between the two parents. Our research provides valuable information for improving grain weight/size and yield in Tartary buckwheat breeding.

Dynamic transcriptome analysis suggests the key genes regulating seed development and filling in Tartary buckwheat (Fagopyrum tataricum Garetn.)

Tartary buckwheat is highly attractive for the richness of nutrients and quality, yet post-embryonic seed abortion greatly halts the yield. Seed development is crucial for determining grain yield, whereas the molecular basis and regulatory network of Tartary buckwheat seed development and filling is not well understood at present. Here, we assessed the transcriptional dynamics of filling stage Tartary buckwheat seeds at three developmental stages by RNA sequencing. Among the 4249 differentially expressed genes (DEGs), genes related to seed development were identified. Specifically, 88 phytohormone biosynthesis signaling genes, 309 TFs, and 16 expansin genes participating in cell enlargement, 37 structural genes involved in starch biosynthesis represented significant variation and were candidate key seed development genes. Cis-element enrichment analysis indicated that the promoters of differentially expressed expansin genes and starch biosynthesis genes are rich of hormone-responsive (ABA-, AUX-, ET-, and JA-), and seed growth-related (MYB, MYC and WRKY) binding sites. The expansin DEGs showed strong correlations with DEGs in phytohormone pathways and transcription factors (TFs). In total, phytohormone ABA, AUX, ET, BR and CTK, and related TFs could substantially regulate seed development in Tartary buckwheat through targeting downstream expansin genes and structural starch biosynthetic genes. This transcriptome data could provide a theoretical basis for improving yield of Tartary buckwheat.

Genome-Wide Identification of NAP1 and Function Analysis in Moso Bamboo (Phyllostachys edulis)

The nucleosome assembly protein 1 (NAP1) family is the main histone chaperone of histone H2A–H2B. To explore the function of NAP1 family genes in moso bamboo (Phyllostachys edulis), characterized by extremely rapid growth and a long flowering cycle, we originally conducted a genome-wide analysis of the PheNAP1 gene. The phylogenetic relationship, gene expression pattern, DNA methylation, and histone modification were analyzed. Eventually, 12 PheNAP1 genes were recognized from the Phyllostachys edulis genome, divided into two sorts: the NRP subfamily (four members) and the NAP subfamily (eight members). Highly conserved motifs exist in each subfamily, which are distinct between subfamilies. PheNAP1 was distributed homogeneously on 10 out of 24 chromosomes, and gene duplication contributed significantly to the enhancement of the PheNAP1 gene in the genome. Cis-acting element analysis showed that PheNAP1 family genes are involved in light, hormone, and abiotic stress responses and may play an important role in the rapid growth and flowering. PheNAP1 exhibited the highest expression level in fast-growing shoots, indicating it is closely associated with the rapid growth of moso bamboo. Besides, PheNAP1 can rescue the early-flowering phenotype of nrp1-1 nrp2-2, and it affected the expression of genes related to the flowering pathway, like BSU1, suggesting the vital role that PheNAP1 may take in the flowering process of moso bamboo. In addition, histone modification results showed that PheNAP1 could bind to phosphorylation-, acetylation-, and methylation-modified histones to further regulate gene expression. A sketch appears: that PheNAP1 can accompany histones to regulate fast-growth- and flowering-related genes in moso bamboo. The consequences of this study enrich the understanding of the epigenetic regulation mechanism of bamboo plants and lays a foundation for further studies on the role of the NAP1 gene in Phyllostachys edulis and the function of chromatin regulation in forest growth and development.

Metal Tolerance Protein Encoding Gene Family in Fagopyrum tartaricum: Genome-Wide Identification, Characterization and Expression under Multiple Metal Stresses

Metal tolerance proteins (MTP) as divalent cation transporters are essential for plant metal tolerance and homeostasis. However, the characterization and the definitive phylogeny of the MTP gene family in Fagopyrum tartaricum, and their roles in response to metal stress are still unknown. In the present study, MTP genes in Fagopyrum tartaricum were identified, and their phylogenetic relationships, structural characteristics, physicochemical parameters, as well as expression profiles under five metal stresses including Fe, Mn, Cu, Zn, and Cd were also investigated. Phylogenetic relationship analysis showed that 12 Fagopyrum tartaricum MTP genes were classified into three major clusters and seven groups. All FtMTPs had typical structural features of the MTP gene family and were predicted to be located in the cell vacuole. The upstream region of FtMTPs contained abundant cis-acting elements, implying their functions in development progress and stress response. Tissue-specific expression analysis results indicated the regulation of FtMTPs in the growth and development of Fagopyrum tataricum. Besides, the expression of most FtMTP genes could be induced by multiple metals and showed different expression patterns under at least two metal stresses. These findings provide useful information for the research of the metal tolerance mechanism and genetic improvement of Fagopyrum tataricum.

Genome-Wide Analysis of the MADS-Box Gene Family in Maize: Gene Structure, Evolution, and Relationships

The MADS-box gene family is one of the largest families in plants and plays an important roles in floral development. The MADS-box family includes the SRF-like domain and K-box domain. It is considered that the MADS-box gene family encodes a DNA-binding domain that is generally related to transcription factors, and plays important roles in regulating floral development. Our study identified 211 MADS-box protein sequences in the Zea mays proteome and renamed all the genes based on the gene annotations. All the 211 MADS-box protein sequences were coded by 98 expressed genes. Phylogenetic analysis of the MADS-box genes showed that all the family members were categorized into five subfamilies: MIKC-type, Mα, Mβ, Mγ, and Mδ. Gene duplications are regarded as products of several types of errors during the period of DNA replication and reconstruction; in our study all the 98 MADS-box genes contained 22 pairs of segmentally duplicated events which were distributed on 10 chromosomes. We compared expression data in different tissues from the female spikelet, silk, pericarp aleurone, ear primordium, leaf zone, vegetative meristem, internode, endosperm crown, mature pollen, embryo, root cortex, secondary root, germination kernels, primary root, root elongation zone, and root meristem. According to analysis of gene ontology pathways, we found a total of 41 pathways in which MADS-box genes in maize are involved. All the studies we conducted provided an overview of MADS-box gene family members in maize and showed multiple functions as transcription factors. The related research of MADS-box domains has provided the theoretical basis of MADS-box domains for agricultural applications.

Identification of bHLH genes through genome-wide association study and antisense expression of ZjbHLH076/ZjICE1 influence tolerance to low temperature and salinity in Zoysia japonica

Abiotic stress greatly affects plant growth and developmental processes, resulting in poor productivity. A variety of basic helix-loop-helix (bHLH) transcription factors (TFs) that play important roles in plant abiotic stress response pathways have been identified. However, bHLH proteins of Zoysia japonica, one of the warm-season turfgrasses, have not been widely studied. In this study, 141 bHLH genes (ZjbHLHs) were identified and classified into 22 subfamilies. The ZjbHLHs were mapped on 19 chromosomes except for Chr17 and one pair of the tandemly arrayed genes was identified on Chr06. Also, the co-linearity of ZjbHLHs was found to have been driven mostly by segmental duplication events. The subfamily IIIb genes of our present interest, possessed various stress responsive cis-elements in their promoters. ZjbHLH076/ZjICE1, a MYC-type bHLH TF in subfamily IIIb was analyzed by overexpression and its loss-of-function via overexpressing a short ZjbHLH076/ZjICE1 fragment in the antisense direction. The overexpression of ZjbHLH076/ZjICE1 enhanced the tolerance to cold and salinity stress in the transgenic Z. japonica plants. However, the anti-sense expression of ZjbHLH076/ZjICE1 showed sensitive to these abiotic stresses. These results suggest that ZjbHLH076/ZjICE1 would be a promising candidate for the molecular breeding program to improve the abiotic stress tolerance of Z. japonica.

Tartary buckwheat database (TBD): an integrative platform for gene analysis of and biological information on Tartary buckwheat

Rice, wheat, corn, and potatoes are four crops that provide a daily source of nutrition for humans, but there are many problems that have been found with these crops. First, they lack amino acids and minerals which are necessary for balanced nutrition, and they also are grown very widely and as monocultures, which increases the risk of the human food system being destroyed by climate change. Thus, by introducing coarse cereals with good characteristics, we can enrich human food resources, realize agricultural diversification, improve dietary structure, and mitigate risks. Tartary buckwheat (Fagopyrum tataricum) is a widely cultivated edible and medicinal crop with unique nutritional and excellent economic value. It contains flavonoids, such as rutin and quercetin, which are not found in cereal crops. Rutin is a major flavonoid that can enhance blood flow and aid in the use of vitamin C and the production of collagen. In addition, such antioxidants have been shown to effectively reduce cholesterol levels, blood clots, and hypertension, particularly for the prevention of inflammatory liver injury (Middleton et al., 2000; Lee et al., 2013; Suzuki et al., 2014; Huang et al., 2016; Nishimura et al., 2016). Meanwhile, Tartary buckwheat can tolerate poor climate and acidic soils containing high amounts of aluminum, which is toxic to other crops (Wang et al., 2015). The self-pollination of Tartary buckwheat has resulted in a decrease in genomic heterozygosity, which is valuable for breeding and a stable production trait (Wang and Campbell, 2007). Therefore, Tartary buckwheat is an important minor crop, which is expected to become the target of many breeding efforts in the future.

Genome-wide analysis of the NF-Y gene family and their roles in relation to fruit development in Tartary buckwheat (Fagopyrum tataricum)

Nuclear factor Y (NF-Y) is a heterotrimeric transcription factor playing crucial roles in various biological process in plant. However, thorough research on NF-Y gene family of Tartary buckwheat (Fagopyrum tataricum) is little. In this study, 38 FtNF-Y genes (12 FtNF-YAs, 17 FtNF-YBs, and 9 FtNF-YCs) were identified and renamed on the basis of their subfamily and chromosomal location. Their gene structure, genomic mapping, motif composition, conserved domain, phylogenetic relationships, cis-acting elements and gene expression were investigated. Illustration of gene structures and conserved domains of FtNF-Ys revealed their functional conservation and specificity. Construction of phylogenetic trees of NF-Ys in Tartary buckwheat, Arabidopsis, tomato, rice and banana, allowed us to predict functional similarities among NF-Ys from different species. Gene expression analysis displayed that twenty-four FtNF-Ys were expressed in all the tissues and the transcript levels of them were different, suggesting their function varieties. Moreover, expression profiles of twenty FtNF-Ys along five different fruit development stages acquired by real-time quantitative PCR (RT-qPCR) demonstrated distinct abundance diversity at different stages, providing some clues of potential fruit development regulators. Our study could provide helpful reference information for further function characterization of FtNF-Ys and for the fruit quality enhancement of Tartary buckwheat.

Genome-wide analysis and characterization of GRAS family in switchgrass

Panicum virgatum, a model plant of cellulosic ethanol conversion, not only has high large biomass and strong adaptability to soil, but also grows well in marginal soil and has the advantage of improving saline-alkali soil. GRAS transcription factor gene family play important roles in individual environment adaption, and these vital functions has been proved in several plants, however, the research of GRAS in the development of switchgrass (Panicum virgatum) were limited. A comprehensive study was investigated to explore the relationship between GRAS gene family and resistance. According to the phylogenetic analysis, a total of 144 GRAS genes were identified and renamed which were classified into eight subfamilies. Chromosome distribution, tandem and segmental repeats analysis indicated that gene duplication events contributed a lot to the expansion of GRAS genes in the switchgrass genome. Sixty-six GRAS genes in switchgrass were identified as having orthologous genes with rice through gene duplication analysis. Most of these GRAS genes contained zero or one intron, and closely related genes in evolution shared similar motif composition. Interaction networks were analyzed including DELLA and ten interaction proteins that were primarily involved in gibberellin acid mediated signaling. Notably, online analysis indicated that the promoter regions of the identified PvGRAS genes contained many cis-elements including light responsive elements, suggesting that PvGRAS might involve in light signal cross-talking. This work provides key insights into resistance and bioavailability in switchgrass and would be helpful to further study the function of GRAS and GRAS-mediated signal transduction pathways.

Genome-wide identification and expression analysis of the trehalose-6-phosphate synthase (TPS) gene family in cucumber (Cucumis sativus L.)

Trehalose-6-phosphate synthase (TPS) is significant in the growth, development and stress resistance of plants. We identified the cucumber TPS family and its physicochemical properties, domains, gene structures, evolutionary relationships, gene locations, cis-acting elements, conserved motifs, and expression patterns using bioinformatics. Our results uncovered seven CsTPS genes in the cucumber genome and named CsTPS1-CsTPS7 according to their locations in the chromosomes. Seven CsTPS genes were randomly distributed in six cucumber chromosomes. Domain analysis showed that the TPS and TPP domains exist in all CsTPSs, and an additional hydrolase-3 domain exist in CsTPS3, CsTPS5 and CsTPS6. Phylogenetic analysis showed that TPS proteins from Arabidopsis, rice, soybean, and cucumber were divided into two subfamilies (Class I and Class II) and they were further divided into seven subgroups. TPS proteins from Arabidopsis and cucumber were grouped together, suggesting a close evolutionary relationship. Gene structure analysis indicated that most Class I genes contained 16-17 introns, while Class II genes (except CsTPS7) had two introns. Motif analysis showed that Class II genes had 10 complete conserved motifs, while Class I genes lacked motif 8 and motif 9. Furthermore, CsTPS genes possessed numerous cis-acting elements related to stress, hormone, and light response in the promoter regions. GO analysis indicated multiple functions for the CsTPS proteins. Expression analysis of CsTPS genes in different tissues found that they were expressed in roots, stems and leaves, with the highest expression levels in roots. The expression analysis of CsTPSs under different treatments showed that CsTPS genes may participate in the response to abiotic stress, plant hormones and sugar treatments.

Genome-Wide Investigation of Major Enzyme-Encoding Genes in the Flavonoid Metabolic Pathway in Tartary Buckwheat (Fagopyrum tataricum)

Key enzymes play a vital role in plant growth and development. However, the evolutionary relationships between genes encoding key enzymes in the metabolic pathway of Tartary buckwheat flavonoids are poorly understood. Based on the published Tartary buckwheat genome sequence and related Tartary buckwheat transcriptome data, 48 key enzyme-encoding genes involved in flavonoid metabolism were screened from the Tartary buckwheat genome in this study; the chromosome localization, gene structure and promoter elements of these enzyme-encoding gene were also investigated. Gene structure analysis revealed relatively conserved 5' exon sequences among the 48 genes, indicating that the structural diversity of key enzyme-encoding genes is low in Tartary buckwheat. Through promoter analysis, these key enzyme-encoding genes were found to contain a large number of light-response elements and hormone-response elements. In addition, some genes could bind MYB transcription factors, participating in the regulation of flavonoid biosynthesis. The transcription level of the 48 key enzyme-encoding gene varied greatly among tissues. In this study, we identified 48 key enzyme-encoding genes involved in flavonoid metabolic pathways, and elucidated the structure, evolution and tissue-specific expression patterns of these genes. These results lay a foundation for further understanding the functional characteristics and evolutionary relationships of key enzyme-encoding genes involved in the flavonoid metabolic pathway in Tartary buckwheat.

Mapping QTLs for 1000-grain weight and genes controlling hull type using SNP marker in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum), an important pseudocereal crop, has high economic value due to its nutritional and medicinal properties. However, dehulling of Tartary buckwheat is difficult owing to its thick and tough hull, which has greatly limited the development of the Tartary buckwheat processing industry. The construction of high-resolution genetic maps serves as a basis for identifying quantitative trait loci (QTLs) and qualitative trait genes for agronomic traits. In this study, a recombinant inbred lines (XJ-RILs) population derived from a cross between the easily dehulled Rice-Tartary type and Tartary buckwheat type was genotyped using restriction site-associated DNA (RAD) sequencing to construct a high-density SNP genetic map. Furthermore, QTLs for 1000-grain weight (TGW) and genes controlling hull type were mapped in multiple environments.

RESULTS

In total, 4151 bin markers comprising 122,185 SNPs were used to construct the genetic linkage map. The map consisted of 8 linkage groups and covered 1444.15 cM, with an average distance of 0.35 cM between adjacent bin markers. Nine QTLs for TGW were detected and distributed on four loci on chromosome 1 and 4. A major locus detected in all three trials was mapped in 38.2-39.8 cM region on chromosome 1, with an LOD score of 18.1-37.0, and explained for 23.6-47.5% of the phenotypic variation. The genes controlling hull type were mapped to chromosome 1 between marker Block330 and Block331, which was closely followed by the major locus for TGW. The expression levels of the seven candidate genes controlling hull type present in the region between Block330 and Block336 was low during grain development, and no significant difference was observed between the parental lines. Six non-synonymous coding SNPs were found between the two parents in the region.

CONCLUSIONS

We constructed a high-density SNP genetic map for the first time in Tartary buckwheat. The mapped major loci controlling TGW and hull type will be valuable for gene cloning and revealing the mechanism underlying grain development and easy dehulling, and marker-assisted selection in Tartary buckwheat.

Genome-wide identification of genes involved in heterotrimeric G-protein signaling in Tartary buckwheat (Fagopyrum tataricum) and their potential roles in regulating fruit development

Tartary buckwheat (Fagopyrum tataricum Gaertn.) is an economical crop with excellent edible, nutritional, and medicinal values. However, the production of Tartary buckwheat is very low and it is urgent to breed high-yield varieties for satisfying the increasing market demand. Heterotrimeric G-protein signaling involves in the regulation of agronomical traits and fruit development in several plant species. In this study, fifteen genes involved in G-protein signaling were characterized in Tartary buckwheat and their potential roles in fruit development were revealed by expression analysis. The exon-intron organization and conserved motif of Tartary buckwheat G-protein signaling genes were similar to those in other dicot plants. All these genes were ubiquitously and differently expressed in five tissues. The expression patterns of Tartary buckwheat G-protein signaling genes in fruit suggested they may play important roles in the fruit at early development stage, which was supported by meta-analysis of G-protein signaling genes' expression in the fruits from different species. Furthermore, we found the expression of G-protein signaling genes in fruit showed high correlation with 178 transcription factors, which indicated a transcriptional regulatory loop moderating G-protein signaling genes' expression during fruit development. This paper provides new insights into the physiological functions of G-protein signaling in fruit.

Comparative cellular, physiological and transcriptome analyses reveal the potential easy dehulling mechanism of rice-tartary buckwheat (Fagopyrum Tararicum)

BACKGROUND

Tartary buckwheat has gained popularity in the food marketplace due to its abundant nutrients and high bioactive flavonoid content. However, its difficult dehulling process has severely restricted its food processing industry development. Rice-tartary buckwheat, a rare local variety, is very easily dehulled, but the cellular, physiological and molecular mechanisms responsible for this easy dehulling remains largely unclear.

RESULTS

In this study, we integrated analyses of the comparative cellular, physiological, transcriptome, and gene coexpression network to insight into the reason that rice-tartary buckwheat is easy to dehull. Compared to normal tartary buckwheat, rice-tartary buckwheat has significantly brittler and thinner hull, and thinner cell wall in hull sclerenchyma cells. Furthermore, the cellulose, hemicellulose, and lignin contents of rice-tartary buckwheat hull were significantly lower than those in all or part of the tested normal tartary buckwheat cultivars, respectively, and the significant difference in cellulose and hemicellulose contents between rice-tartary buckwheat and normal tartary buckwheat began at 10 days after pollination (DAP). Comparative transcriptome analysis identified a total of 9250 differentially expressed genes (DEGs) between the rice- and normal-tartary buckwheat hulls at four different development stages. Weighted gene coexpression network analysis (WGCNA) of all DEGs identified a key module associated with the formation of the hull difference between rice- and normal-tartary buckwheat. In this specific module, many secondary cell wall (SCW) biosynthesis regulatory and structural genes, which involved in cellulose and hemicellulose biosynthesis, were identified as hub genes and displayed coexpression. These identified hub genes of SCW biosynthesis were significantly lower expression in rice-tartary buckwheat hull than in normal tartary buckwheat at the early hull development stages. Among them, the expression of 17 SCW biosynthesis relative-hub genes were further verified by quantitative real-time polymerase chain reaction (qRT-PCR).

CONCLUSIONS

Our results showed that the lower expression of SCW biosynthesis regulatory and structural genes in rice-tartary buckwheat hull in the early development stages contributes to its easy dehulling by reducing the content of cell wall chemical components, which further effects the cell wall thickness of hull sclerenchyma cells, and hull thickness and mechanical strength.

Basic helix-loop-helix (bHLH) gene family in Tartary buckwheat (Fagopyrum tataricum): Genome-wide identification, phylogeny, evolutionary expansion and expression analyses

Tartary buckwheat (Fagopyrum tataricum) a kind of edible and medicinal plant, is of great nutritional value. It is difficult to remove the hull of Tartary buckwheat fruit and breeding new easy-dehulled varieties has been one of the major breeding objectives. The bHLH gene family plays a vital role in plant growth and fruit dehiscence. In order to improve Tartary buckwheat breeding, we need to study the bHLH gene family for excavating genes with potential regulation of fruit development and dehiscence. Here, 164 Fagopyrum tataricum bHLH (FtbHLH) genes were identified. Analyses of gene structure and motif composition illustrate that the members of specific FtbHLH subfamily are relatively conserved. Synteny and phylogenetic analyses of bHLH genes in Tartary buckwheat and other plants lay a foundation for further exploring the evolutionary characteristic of the FtbHLH genes (FtbHLHs). qRT-PCR experiments showed that FtbHLHs expression patterns were different in plant organs, indicating that they may perform diverse functions. In addition, some genes that potentially regulate flower and fruit development and easy dehulling were screened out. Overall, this study will be helpful for further analyzing the biological function of FtbHLHs and provides clues for improving the genetic breeding and economic value of the Tartary buckwheat.

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.

Genome-wide identification, expression analysis and functional study of the GRAS gene family in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

GRAS are plant-specific transcription factors that play important roles in plant growth and development. Although the GRAS gene family has been studied in many plants, there has been little research on the GRAS genes of Tartary buckwheat (Fagopyrum tataricum), which is an important crop rich in rutin. The recently published whole genome sequence of Tartary buckwheat allows us to study the characteristics and expression patterns of the GRAS gene family in Tartary buckwheat at the genome-wide level.

RESULTS

In this study, 47 GRAS genes of Tartary buckwheat were identified and divided into 10 subfamilies: LISCL, HAM, DELLA, SCR, PAT1, SCL4/7, LAS, SHR, SCL3, and DLT. FtGRAS genes were unevenly distributed on 8 chromosomes, and members of the same subfamily contained similar gene structures and motif compositions. Some FtGRAS genes may have been produced by gene duplications; tandem duplication contributed more to the expansion of the GRAS gene family in Tartary buckwheat. Real-time PCR showed that the transcription levels of FtGRAS were significantly different in different tissues and fruit development stages, implying that FtGRAS might have different functions. Furthermore, an increase in fruit weight was induced by exogenous paclobutrazol, and the transcription level of the DELLA subfamily member FtGRAS22 was significantly upregulated during the whole fruit development stage. Therefore, FtGRAS22 may be a potential target for molecular breeding or genetic editing.

CONCLUSIONS

Collectively, this systematic analysis lays a foundation for further study of the functional characteristics of GRAS genes and for the improvement of Tartary buckwheat crops.

Genome-wide identification, phylogeny, evolutionary expansion and expression analyses of bZIP transcription factor family in tartaty buckwheat

Background

In reported plants, the bZIP family is one of the largest transcription factor families. bZIP genes play roles in the light signal, seed maturation, flower development, cell elongation, seed accumulation protein, abiotic and biological stress and other biological processes. While, no detailed identification and genome-wide analysis of bZIP family genes in Fagopyum talaricum (tartary buckwheat) has previously been published. The recently reported genome sequence of tartary buckwheat provides theoretical basis for us to study and discuss the characteristics and expression of bZIP genes in tartary buckwheat based on the whole genome.

Results

In this study, 96 FtbZIP genes named from FtbZIP1 to FtbZIP96 were identified and divided into 11 subfamilies according to their genetic relationship with 70 bZIPs of A. thaliana. FtbZIP genes are not evenly distributed on the chromosomes, and we found tandem and segmental duplication events of FtbZIP genes on 8 tartary buckwheat chromosomes. According to the results of gene and motif composition, FtbZIP located in the same group contained analogous intron/exon organizations and motif composition. By qRT-PCR, we quantified the expression of FtbZIP members in stem, root, leaf, fruit, and flower and during fruit development. Exogenous ABA treatment increased the weight of tartary buckwheat fruit and changed the expressions of FtbZIP genes in group A.

Conclusions

Through our study, we identified 96 FtbZIP genes in tartary buckwheat and synthetically further analyzed the structure composition, evolution analysis and expression pattern of FtbZIP proteins. The expression pattern indicates that FtbZIP is important in the course of plant growth and development of tartary buckwheat. Through comprehensively analyzing fruit weight and FtbZIP genes expression after ABA treatment and endogenous ABA content of tartary buckwheat fruit, ABA may regulate downstream gene expression by regulating the expression of FtPinG0003523300.01 and FtPinG0003196200.01, thus indirectly affecting the fruit development of tartary buckwheat. This will help us to further study the function of FtbZIP genes in the tartary buckwheat growth and improve the fruit of tartary buckwheat.

Electronic supplementary material

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

Genome-wide analysis of the NAC transcription factor family in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

The NAC (NAM, ATAF1/2, and CUC2) transcription factor family represents a group of large plant-specific transcriptional regulators, participating in plant development and response to external stress. However, there is no comprehensive study on the NAC genes of Tartary buckwheat (Fagopyrum tataricum), a large group of extensively cultivated medicinal and edible plants. The recently published Tartary buckwheat genome permits us to explore all the FtNAC genes on a genome-wide basis.

RESULTS

In the present study, 80 NAC (FtNAC) genes of Tartary buckwheat were obtained and named uniformly according to their distribution on chromosomes. Phylogenetic analysis of NAC proteins in both Tartary buckwheat and Arabidopsis showed that the FtNAC proteins are widely distributed in 15 subgroups with one subgroup unclassified. Gene structure analysis found that multitudinous FtNAC genes contained three exons, indicating that the structural diversity in Tartary buckwheat NAC genes is relatively low. Some duplication genes of FtNAC have a conserved structure that was different from others, indicating that these genes may have a variety of functions. By observing gene expression, we found that FtNAC genes showed abundant differences in expression levels in various tissues and at different stages of fruit development.

CONCLUSIONS

In this research, 80 NAC genes were identified in Tartary buckwheat, and their phylogenetic relationships, gene structures, duplication, global expression and potential roles in Tartary buckwheat development were studied. Comprehensive analysis will be useful for a follow-up study of functional characteristics of FtNAC genes and for the development of high-quality Tartary buckwheat varieties.