Genome-wide study of pineapple (Ananas comosus L.) bHLH transcription factors indicates that cryptochrome-interacting bHLH2 (AcCIB2) participates in flowering time regulation and abiotic stress response

Rice CRYPTOCHROME-INTERACTING BASIC HELIX-LOOP-HELIX 1-LIKE interacts with OsCRY2 and promotes flowering by upregulating Early heading date 1

Flowering time is a crucial adaptive response to seasonal variation in plants and is regulated by environmental cues such as photoperiod and temperature. In this study, we demonstrated the regulatory function of rice CRYPTOCHROME-INTERACTING BASIC HELIX-LOOP-HELIX 1-LIKE (OsCIBL1) in flowering time. Overexpression of OsCIB1L promoted flowering, whereas the oscib1l knockout mutation did not alter flowering time independent of photoperiodic conditions. Cryptochromes (CRYs) are blue light photoreceptors that enable plants to sense photoperiodic changes. OsCIBL1 interacted with OsCRY2, a member of the rice CRY family (OsCRY1a, OsCRY1b, and OsCRY2), and bound to the Early heading date 1 (Ehd1) promoter, activating the rice-specific Ehd1-Heading date 3a/RICE FLOWERING LOCUS T 1 pathway for flowering induction. Dual-luciferase reporter assays showed that the OsCIBL1-OsCRY2 complex required blue light to induce Ehd1 transcription. Natural alleles resulting from nonsynonymous single nucleotide polymorphisms in OsCIB1L and OsCRY2 may contribute to the adaptive expansion of rice cultivation areas. These results expand our understanding of the molecular mechanisms controlling rice flowering and highlight the importance of blue light-responsive genes in the geographic distribution of rice.

A natural mutation in the promoter of the aconitase gene ZjACO3 influences fruit citric acid content in jujube

Jujube (Ziziphus jujuba Mill.) is the most economically important fruit tree of the Rhamnaceae and was domesticated from wild or sour jujube (Z. jujuba Mill. var. spinosa Hu). During the process of domestication, there was a substantial reduction in the content of organic acids, particularly malate and citrate, which greatly influence the taste and nutritional value of the fruit. We previously demonstrated that ZjALMT4 is crucial for malate accumulation. However, the mechanism of citrate degradation in jujube remains poorly understood. In the present study, aconitase ZjACO3 was shown to participate in citric acid degradation in the cytoplasm through the GABA pathway. Interestingly, we discovered an E-box mutation in the ZjACO3 promoter (-484A > G; CAAGTG in sour jujube mutated to CAGGTG in cultivated jujube) that was strongly correlated with fruit citrate content; 'A' represented a high-citrate genotype and 'G' represented a low-citrate genotype. We developed and validated an ACO-based Kompetitive allele-specific PCR (KASP) marker for determining citric acid content. Yeast one-hybrid screening, transient dual-luciferase assays, and overexpression analyses showed that the transcription factor ZjbHLH113 protein directly binds to CAGGTG in the promoter of ZjACO3 in cultivated jujube plants, transcriptionally activating ZjACO3 expression, and enhancing citric acid degradation. Conversely, binding ability of the ZjbHLH113 protein to CAAGTG was weakened in sour jujube, thereby promoting citrate accumulation in the fruit. These findings will assist in elucidating the mechanism by which ZjACO3 modulates citrate accumulation in sour jujube and its cultivars.

Basic helix-loop-helix (bHLH) gene family in rye (Secale cereale L.): genome-wide identification, phylogeny, evolutionary expansion and expression analyses

BACKGROUND

Rye (Secale cereale), one of the drought and cold-tolerant crops, is an important component of the Triticae Dumortier family of Gramineae plants. Basic helix-loop-helix (bHLH), an important family of transcription factors, has played pivotal roles in regulating numerous intriguing biological processes in plant development and abiotic stress responses. However, no systemic analysis of the bHLH transcription factor family has yet been reported in rye.

RESULTS

In this study, 220 bHLH genes in S. cereale (ScbHLHs) were identified and named based on the chromosomal location. The evolutionary relationships, classifications, gene structures, motif compositions, chromosome localization, and gene replication events in these ScbHLH genes are systematically analyzed. These 220 ScbHLH members are divided into 21 subfamilies and one unclassified gene. Throughout evolution, the subfamilies 5, 9, and 18 may have experienced stronger expansion. The segmental duplications may have contributed significantly to the expansion of the bHLH family. To systematically analyze the evolutionary relationships of the bHLH family in different plants, we constructed six comparative genomic maps of homologous genes between rye and different representative monocotyledonous and dicotyledonous plants. Finally, the gene expression response characteristics of 22 ScbHLH genes in various biological processes and stress responses were analyzed. Some candidate genes, such as ScbHLH11, ScbHLH48, and ScbHLH172, related to tissue developments and environmental stresses were screened.

CONCLUSIONS

The results indicate that these ScbHLH genes exhibit characteristic expression in different tissues, grain development stages, and stress treatments. These findings provided a basis for a comprehensive understanding of the bHLH family in rye.

Comparative transcriptome analysis reveals candidate genes for cold stress response and early flowering in pineapple

Pineapple originates from tropical regions in South America and is therefore significantly impacted by cold stress. Periodic cold events in the equatorial regions where pineapple is grown may induce early flowering, also known as precocious flowering, resulting in monetary losses due to small fruit size and the need to make multiple passes for harvesting a single field. Currently, pineapple is one of the most important tropical fruits in the world in terms of consumption, and production losses caused by weather can have major impacts on worldwide exportation potential and economics. To further our understanding of and identify mechanisms for low-temperature tolerance in pineapple, and to identify the relationship between low-temperature stress and flowering time, we report here a transcriptomic analysis of two pineapple genotypes in response to low-temperature stress. Using meristem tissue collected from precocious flowering-susceptible MD2 and precocious flowering-tolerant Dole-17, we performed pairwise comparisons and weighted gene co-expression network analysis (WGCNA) to identify cold stress, genotype, and floral organ development-specific modules. Dole-17 had a greater increase in expression of genes that confer cold tolerance. The results suggested that low temperature stress in Dole-17 plants induces transcriptional changes to adapt and maintain homeostasis. Comparative transcriptomic analysis revealed differences in cuticular wax biosynthesis, carbohydrate accumulation, and vernalization-related gene expression between genotypes. Cold stress induced changes in ethylene and abscisic acid-mediated pathways differentially between genotypes, suggesting that MD2 may be more susceptible to hormone-mediated early flowering. The differentially expressed genes and module hub genes identified in this study are potential candidates for engineering cold tolerance in pineapple to develop new varieties capable of maintaining normal reproduction cycles under cold stress. In addition, a total of 461 core genes involved in the development of reproductive tissues in pineapple were also identified in this study. This research provides an important genomic resource for understanding molecular networks underlying cold stress response and how cold stress affects flowering time in pineapple.

Genome-Wide Characterization of bHLH Family Genes and Expression Analysis in Response to Osmotic Stress in Betula platyphylla

The bHLH family, as a superfamily of transcription factors (TFs), has special functional characteristics in plants and plays a crucial role in a plant's growth and development and helping the plant cope with various stresses. In this study, 128 bHLH family genes were screened in the birch (B. platyphylla) genome using conservative domain scan and blast analysis. These genes are clustered into 21 subfamilies based on the phylogenetic tree construction and are unevenly distributed among the 14 birch chromosomes. In all, 22 segmental duplication pairs with 27 BpbHLH genes were identified. The duplications were distributed on eight chromosomes. Analysis of gene structures and protein motifs revealed intra-group conservation of BpbHLHs. Of the BpbHLH family genes, 16 contain only one intron each. The BPChr14G06667 gene contains the most introns, that is, 19. The cis-elements, which respond to plant hormones, light, defense, and stress, were found on the promoter of BHLH family genes. As per RNA-seq data analysis, under PEG osmotic stress, most BpbHLH genes were differentially expressed, and eight were highly differentially expressed. The qRT-PCR analysis results further indicated that BPChr06G09475 was the gene with the highest expression level in leaves, roots, and stems, and that the expression of these eight genes was higher in leaves than in roots and stems and upregulated in all three tissues under osmotic stress compared to the controls. The above analysis suggests that the BpbHLH family genes have a certain biological effect under drought stress that provides a basis for molecular breeding for stress resistance in birch.

VabHLH137 promotes proanthocyanidin and anthocyanin biosynthesis and enhances resistance to Colletotrichum gloeosporioides in grapevine

Proanthocyanidins (PAs) and anthocyanins are involved in the response of plants to various environmental stresses. However, the mechanism behind defense-induced PA biosynthetic regulation is still not completely elucidated, also in grapevine. This study performed a transcriptome sequencing analysis of grape berries infected with Colletotrichum gloeosporioides to highlight the induction of the VabHLH137 factor from the basic helix-loop-helix (bHLH) XII subfamily by the fungus, which appeared to be significantly co-expressed with PA-related genes. The functional analysis of VabHLH137 overexpression and knockdown in transgenic grape calli showed that it positively regulated PA and anthocyanin biosynthesis. Moreover, VabHLH137 overexpression in the grape calli significantly increased resistance to C. gloeosporioides. A yeast one-hybrid and electrophoretic mobility shift assay revealed that VabHLH137 directly bound to the VaLAR2 promoter, enhancing its activity and interacting with VaMYBPAR, a transcriptional activator of PA biosynthesis. Furthermore, transient experiments showed that although the VabHLH137 + VaMYBPAR complex activated VaLAR2 expression, it failed to further enhance VaLAR2 expression compared to VaMYBPAR alone. The findings indicated that VabHLH137 enhanced PA biosynthesis by activating of VaLAR2 expression, providing new insight into the transcriptional regulation of defense-induced PA biosynthesis in grapevine.

Genome-wide identification and expression analysis of the bHLH gene family in passion fruit (Passiflora edulis) and its response to abiotic stress

Passion fruit is a tropical fruit crop with significant agricultural, economic and ornamental values. The growth and development of passion fruit are greatly affected by climatic conditions. In plants, the basic helix-loop-helix (bHLH) gene family plays essential roles in the floral organ and fruit development, as well as stress response. However, the characteristics and functions of the bHLH genes of passion fruit remain unclear. Here, 138 passion fruit bHLH members were identified and classified into 20 subfamilies. The structural analysis illustrated that PebHLH proteins of the specific subfamily are relatively conserved. Collinearity analysis indicated that the expansion of the PebHLH gene family mainly took place by segmental duplication, and the structural diversity of duplicated genes might contribute to their functional diversity. PebHLHs, which potentially regulate different floral organ and fruit development, were further screened out, and many of these genes were differentially expressed under various stress treatments. The co-presence of different cis-regulatory elements involved in developmental regulation, hormone and stress responses in the promoter regions of PebHLHs might be closely related to their diverse regulatory roles. Overall, this study will be helpful for further functional investigation of PebHLHs and provides clues for improvement of the passion fruit breeding.

Basic Helix-Loop-Helix Transcription Factors: Regulators for Plant Growth Development and Abiotic Stress Responses

Plant basic helix-loop-helix (bHLH) transcription factors are involved in many physiological processes, and they play important roles in the abiotic stress responses. The literature related to genome sequences has increased, with genome-wide studies on the bHLH transcription factors in plants. Researchers have detailed the functionally characterized bHLH transcription factors from different aspects in the model plant Arabidopsis thaliana, such as iron homeostasis and abiotic stresses; however, other important economic crops, such as rice, have not been summarized and highlighted. The bHLH members in the same subfamily have similar functions; therefore, unraveling their regulatory mechanisms will help us to identify and understand the roles of some of the unknown bHLH transcription factors in the same subfamily. In this review, we summarize the available knowledge on functionally characterized bHLH transcription factors according to four categories: plant growth and development; metabolism synthesis; plant signaling, and abiotic stress responses. We also highlight the roles of the bHLH transcription factors in some economic crops, especially in rice, and discuss future research directions for possible genetic applications in crop breeding.

Identification and expression analysis of pineapple sugar transporters reveal their role in the development and environmental response

In plants, sugars are required for several essential functions, including growth, storage, signaling, defense and reproduction. Sugar transporters carry out the controlled movement of sugars from source (leaves) to sink (fruits and roots) tissues and determine the overall development of the plant. Various types of sugar transporter families have been described in plants, including sucrose transporters (SUC/SUT), monosaccharide transporter (MST) and SWEET (from "Sugar Will Eventually be Exported Transporters"). However, the information about pineapple sugar transporters is minimal. This study systematically identified and classified 45 MST and 4 SUC/SUT genes in the pineapple genome. We found that the expression patterns of sugar transporter genes have a spatiotemporal expression in reproductive and vegetative tissues indicating their pivotal role in reproductive growth and development. Besides, different families of sugar transporters have a diel expression pattern in photosynthetic and non-photosynthetic tissues displaying circadian rhythm associated participation of sugar transporters in the CAM pathway. Moreover, regulation of the stress-related sugar transporters during cold stress indicates their contribution to cold tolerance in pineapple. Heterologous expression (yeast complementation assays) of sugar transporters in a mutant yeast strain suggested that SUT1/2 have the ability to transport sucrose, and STP13, STP26, pGlcT-L2 and TMT4 are able to transport glucose, whereas SWEET11/13 transport both sucrose and fructose. The information provided here would help researchers further explore the underlying molecular mechanism involved in the sugar metabolism of pineapple.

Genome-wide identification and expression analysis of the bHLH gene family in cauliflower (Brassica oleracea L.)

Basic helix-loop-helix (bHLH) transcription factors (TFs) are one of the largest TF families in plant species, and they play important roles in plant growth, development and stress responses. The present study systematically identified members of the cauliflower (Brassica oleracea L.) bHLH gene family based on genomic data. Analysis of bHLH family gene numbers, evolution, collinearity, gene structures and motifs indicated that cauliflower contained 256 bHLH family genes distributed on 10 chromosomes. Most of these genes have been localized in the nucleus, and they were divided into 18 subgroups which have been relatively conserved during evolution. Promoter analysis showed that most cis-acting elements were related to MeJA and ABA. Expression analysis suggested that 14 bHLH genes may be involved in the transformation of cauliflower curd from white to purple. An expression analysis of these 14 genes in FQ136 material was performed using qRT-PCR, and 9 bHLH genes (BobHLH1, 14, 58, 61, 63, 84, 231, 239 and 243) showed significantly increased or decreased expression in cauliflower from white to purple, which suggests that these 9 genes play important roles in the accumulation of anthocyanins in cauliflower. The coexpression network of these 9 genes and anthocyanin synthesis-related key genes was analyzed using weighted gene coexpression network analysis (WGCNA). In conclusion, our observations suggested that the bHLH gene family plays an important role in the accumulation of anthocyanins in cauliflower and provide an important theoretical basis for further research on the functions of the bHLH gene family and the molecular mechanism of cauliflower coloration.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01238-9.

Genome-wide characterization and expression analysis of bHLH gene family in physic nut (Jatropha curcas L.)

The basic helix loop helix (bHLH) transcription factor perform essential roles in plant development and abiotic stress. Here, a total of 122 bHLH family members were identified from the physic nut (Jatropha curcas L.) genomic database. Chromosomal localization results showed that 120 members were located on 11 chromosomes. The phylogenetic tree manifested that the JcbHLHs could be grouped into 28 subfamilies. Syntenic analysis showed that there were 10 bHLH collinear genes among the physic nut, Arabidopsis thaliana and Oryza sativa. These genes, except JcbHLH84, were highly expressed in various tissues of the physic nut, implying a key role in plant development. Gene expression profiles showed that ten genes (especially JcbHLH33, JcbHLH45 and JcbHLH55) correspond to both salinity and drought stresses; while eight genes only respond to salinity and another eight genes only respond to drought stress. Moreover, the protein interaction network revealed that the JcbHLHs are involved in growth, development and stress signal transduction pathways. These discoveries will help to excavate several key genes may involve in salt or drought stresses and seed development, elucidate the complex transcriptional regulation mechanism of JcbHLH genes and provide the theoretical basis for stress response and genetic improvement of physic nut.

AcCIPK5, a pineapple CBL-interacting protein kinase, confers salt, osmotic and cold stress tolerance in transgenic Arabidopsis

Plant-specific calcineurin B-like proteins (CBLs) and their interacting kinases, CBL-interacting protein kinases (CIPKs) module, are essential for dealing with various biotic and abiotic stress. The kinases (CIPKs) of this module have been well studied in several plants; however, the information about pineapple CIPKs remains limited. To understand how CIPKs function against environmental cues in pineapple, the CIPK5 gene of pineapple was cloned and characterized. The phylogenetic analyses revealed that AcCIPK5 is homologous to the CIPK12 of Arabidopsis and other plant species. Quantitative real-time PCR (qRT-PCR) analysis revealed that AcCIPK5 responds to multiple stresses, including osmotic, salt stress, heat and cold. Under optimal conditions, AcCIPK5 gets localized to the cytoplasm and cell membrane. The ectopic expression of AcCIPK5 in Arabidopsis improved the germination under osmotic and salt stress. Furthermore, AcCIPK5 positively regulated osmotic, drought, salt and cold tolerance and negatively regulated heat and fungal stress in Arabidopsis. Besides, the expression of AcCIPK impacted ABA-related genes and ROS homeostasis. Overall, the present study demonstrates that AcCIPK5 contributes to multiple stress tolerance and has the potential to be utilized in the development of stress-tolerant crops.

Integrative Analyses of Transcriptomes and Metabolomes Reveal Associated Genes and Metabolites with Flowering Regulation in Common Vetch (Vicia sativa L.)

As an important source of protein for livestock and human consumption, Vicia sativa is cultivated worldwide, but its seed production is hampered at high altitudes because of the short frost-free period. Flowering represents the transition from a vegetative to a reproductive period, and early flowering benefits plant seed production at high altitudes. However, the molecular mechanisms of flowering regulation in V. sativa remain elusive. In the present study, two V. sativa accessions with different flowering characteristics were used: Lan3 (early-flowering) was cultivated by our laboratory, and 503 (late-flowering) was selected from 222 V. sativa accessions after three years of field experiments. The shoot samples (shoot tip length = 10 cm) of these two accessions were collected 63, 70, and 77 days after sowing, and the molecular regulatory mechanism of the flowering process was identified by integrative analyses of the transcriptomes and metabolomes. Kyoto Encyclopedia of Genes and Genomes enrichment showed that the synthesis and signal transduction of plant hormone pathways were the most enriched pathways in 4274 differentially expressed genes (DEGs) and in 259 differential metabolites between Lan3 and 503. Moreover, the contents of three metabolites related to salicylic acid biosynthesis and the transcription levels of two DEGs related to salicylic acid signal transduction in Lan3 were higher than those in 503. Further verification in various accessions indicated that salicylic acid metabolism may be involved in the flowering regulation process of V. sativa. These findings provide valuable information for understanding the flowering mechanism and for promoting breeding research in V. sativa.

Integrated Analysis of Basic Helix Loop Helix Transcription Factor Family and Targeted Terpenoids Reveals Candidate AarbHLH Genes Involved in Terpenoid Biosynthesis in Artemisia argyi

Artemisia argyi is a valuable traditional medicinal plant in Asia. The essential oil from its leaves is rich in terpenoids and has been used to enhance health and well-being. In China, the market scale of industries related to A. argyi has attained tens of billions of Chinese Yuan. The basic helix-loop-helix (bHLH) family is one of the largest transcription factors families in plants that plays crucial roles in diverse biological processes and is an essential regulatory component of terpenoid biosynthesis. However, the bHLH TFs and their regulatory roles in A. argyi remain unknown. Here, 53 AarbHLH genes were identified from the transcriptome of A. argyi and were classified into 15 subfamilies based on the classification of bHLH proteins in Arabidopsis thaliana. The MEME analysis showed that the conserved motif 1 and motif 2 constituted the most conserved bHLH domain and distributed in most AarbHLH proteins. Additionally, integrated analysis of the expression profiles of AarbHLH genes and the contents of targeted terpenoids in different tissues group and JA-treated group were performed. Eleven up-regulated AarbHLHs and one down-regulated AarbHLH were screened as candidate genes that may participate in the regulation of terpenoid biosynthesis (TPS-AarbHLHs). Correlation analysis between gene expression and terpenoid contents indicated that the gene expression of these 12 TPS-AarbHLHs was significantly correlated with the content changes of 1,8-cineole or β-caryophyllene. Protein-protein interaction networks further illustrated that these TPS-AarbHLHs might be involved in terpenoid biosynthesis in A. argyi. This finding provides a basis to further investigate the regulation mechanism of AarbHLH genes in terpenoid biosynthesis, and will be helpful to improve the quality of A. argyi.

AtCRY2 Negatively Regulates the Functions of AtANN2 and AtANN3 in Drought Tolerance by Affecting Their Subcellular Localization and Transmembrane Ca2+ Flow

Annexins are a multifunctional class of calcium-binding proteins in plants, and their physiological functions and regulation in response to drought stress remain to be elucidated. Here, we found that AtANN2 and AtANN3 conferred to drought tolerance under short-day and long-day conditions, respectively. Under their functional photoperiod, AtANN2 and AtANN3 gene expression was enhanced in the mannitol-treated roots, and their encoded proteins were rapidly targeted to the plasma membrane, and mediated significant Ca²⁺ flows across the plasma membrane. Cryptochromes as photoreceptors can not only sense the photoperiod and regulate ion channels on the plasma membrane to influence ion flow but also induce downstream physiological responses. AtCRY2 repressed the functions of AtANN2 and AtANN3 by affecting their plasma membrane localization and inhibited AtANN2- and AtANN3-dependent transmembrane Ca²⁺ flow in response to drought stress. Taken together, these results uncover a mechanism linking Annexins-AtCRY2 to transmembrane Ca²⁺ flow and resulting in enhanced drought tolerance in Arabidopsis.

Genome-wide identification and expression analysis of the bHLH transcription factor family and its response to abiotic stress in foxtail millet (Setaria italica L.)

Background

Members of the basic helix-loop-helix (bHLH) transcription factor family perform indispensable functions in various biological processes, such as plant growth, seed maturation, and abiotic stress responses. However, the bHLH family in foxtail millet (Setaria italica), an important food and feed crop, has not been thoroughly studied.

Results

In this study, 187 bHLH genes of foxtail millet (SibHLHs) were identified and renamed according to the chromosomal distribution of the SibHLH genes. Based on the number of conserved domains and gene structure, the SibHLH genes were divided into 21 subfamilies and two orphan genes via phylogenetic tree analysis. According to the phylogenetic tree, the subfamilies 15 and 18 may have experienced stronger expansion in the process of evolution. Then, the motif compositions, gene structures, chromosomal spread, and gene duplication events were discussed in detail. A total of sixteen tandem repeat events and thirty-eight pairs of segment duplications were identified in bHLH family of foxtail millet. To further investigate the evolutionary relationship in the SibHLH family, we constructed the comparative syntenic maps of foxtail millet associated with representative monocotyledons and dicotyledons species. Finally, the gene expression response characteristics of 15 typical SibHLH genes in different tissues and fruit development stages, and eight different abiotic stresses were analysed. The results showed that there were significant differences in the transcription levels of some SibHLH members in different tissues and fruit development stages, and different abiotic stresses, implying that SibHLH members might have different physiological functions.

Conclusions

In this study, we identified 187 SibHLH genes in foxtail millet and further analysed the evolution and expression patterns of the encoded proteins. The findings provide a comprehensive understanding of the bHLH family in foxtail millet, which will inform further studies on the functional characteristics of SibHLH genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08095-y.

Genome-wide identification and expression analysis of the bHLH transcription factor family and its response to abiotic stress in sorghum [Sorghum bicolor (L.) Moench]

Background

Basic helix-loop-helix (bHLH) is a superfamily of transcription factors that is widely found in plants and animals, and is the second largest transcription factor family in eukaryotes after MYB. They have been shown to be important regulatory components in tissue development and many different biological processes. However, no systemic analysis of the bHLH transcription factor family has yet been reported in Sorghum bicolor.

Results

We conducted the first genome-wide analysis of the bHLH transcription factor family of Sorghum bicolor and identified 174 SbbHLH genes. Phylogenetic analysis of SbbHLH proteins and 158 Arabidopsis thaliana bHLH proteins was performed to determine their homology. In addition, conserved motifs, gene structure, chromosomal spread, and gene duplication of SbbHLH genes were studied in depth. To further infer the phylogenetic mechanisms in the SbbHLH family, we constructed six comparative syntenic maps of S. bicolor associated with six representative species. Finally, we analyzed the gene-expression response and tissue-development characteristics of 12 typical SbbHLH genes in plants subjected to six different abiotic stresses. Gene expression during flower and fruit development was also examined.

Conclusions

This study is of great significance for functional identification and confirmation of the S. bicolor bHLH superfamily and for our understanding of the bHLH superfamily in higher plants.

Supplementary Information

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