Diversification of R2R3-MYB Transcription Factors in the Tomato Family Solanaceae

In silico analysis of R2R3-MYB transcription factors in the basal eudicot model, Aquilegia coerulea

R2R3-MYBs are an important group of transcription factors that regulate crucial developmental processes across the plant kingdom; yet no comprehensive analysis of the R2R3-MYBs in the early-diverging eudicot clade of Ranunculaceae has been conducted so far. In the present study, Aquilegia coerulea is chosen to understand the extent of conservation and divergence of R2R3-MYBs as a representative of the family by analysing the genomic distribution, organization, gene structure, physiochemical properties, protein architecture, evolution and possible mode of expansion. Genome-wide analysis showed the presence of 82 putative homologues classified into 21 subgroups, based on phylogenetic analysis of full-length protein sequences. The domain has remained largely conserved across all homologues with few differences from the characterized Arabidopsis thaliana R2R3-MYBs. The topology of the phylogenetic tree remains the same when full-length protein sequences are used, indicating that the evolution of R2R3-MYBs is driven by the domain region only. This is supported by the presence of similar structures of exon-intron and conserved motifs within the same subgroup. Furthermore, comparisons of the AqcoeR2R3-MYB members with monocots and core-eudicots revealed the evolutionary expansion of a few functional clades, such as A. thaliana R2R3-MYB subgroup 6 (SG6), the upstream regulatory factors of floral pigment biosynthesis and floral color. The reconstructed evolutionary history of SG6-like genes across angiosperms highlights the occurrence of independent duplication events in the genus Aquilegia. AqcoeR2R3-MYB genes are present in all seven chromosomes of A. coerulea, most of which result from local and segmental duplications. Selection analysis of these duplicated gene pairs indicates purifying selection except one, and the physiochemical analyses of R2R3-MYBs reveal differences among the MYBs signifying their functional diversification. This study paves the way for further investigation of paralogous copies and their probable role in the evolution of different floral traits in A. coerulea. It lays the foundation for functional genomic studies of R2R3-MYBs in the basal eudicots and facilitates comparative studies among angiosperms. The work also provides a framework for deciphering novel genetic regulatory pathways that govern the diversity of floral morphology.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04119-y.

Genetic diversity analysis and marker-trait associations in Amaranthus species

Amaranth (Amaranthus spp.) is a highly nutritious, underutilized vegetable and pseudo-cereal crop. It possesses diverse abiotic stress tolerance traits, is genetically diverse and highly phenotypically plastic, making it an ideal crop to thrive in a rapidly changing climate. Despite considerable genetic diversity there is a lack of detailed characterization of germplasm or population structures. The present study utilized the DArTSeq platform to determine the genetic relationships and population structure between 188 amaranth accessions from 18 agronomically important vegetable, grain, and weedy species. A total of 74, 303 SNP alleles were generated of which 63, 821 were physically mapped to the genome of the grain species A. hypochondriacus. Population structure was inferred in two steps. First, all 188 amaranth accessions comprised of 18 species and second, only 120 A. tricolor accessions. After SNP filtering, a total of 8,688 SNPs were generated on 181 amaranth accessions of 16 species and 9,789 SNPs generated on 118 A. tricolor accessions. Both SNP datasets produced three major sub-populations (K = 3) and generate consistent taxonomic classification of the amaranth sub-genera (Amaranthus Amaranthus, Amaranthus Acnida and Amaranthus albersia), although the accessions were poorly demarcated by geographical origin and morphological traits. A. tricolor accessions were well discriminated from other amaranth species. A genome-wide association study (GWAS) of 10 qualitative traits revealed an association between specific phenotypes and genetic variants within the genome and identified 22 marker trait associations (MTAs) and 100 MTAs (P≤0.01, P≤0.001) on 16 amaranth species and 118 A.tricolor datasets, respectively. The release of SNP markers from this panel has produced invaluable preliminary genetic information for phenotyping and cultivar improvement in amaranth species.

Transcription factors evolve faster than their structural gene targets in the flavonoid pigment pathway

Dissecting the relationship between gene function and substitution rates is key to understanding genome-wide patterns of molecular evolution. Biochemical pathways provide powerful systems for investigating this relationship because the functional role of each gene is often well characterized. Here, we investigate the evolution of the flavonoid pigment pathway in the colorful Petunieae clade of the tomato family (Solanaceae). This pathway is broadly conserved in plants, both in terms of its structural elements and its MYB, basic helix–loop–helix, and WD40 transcriptional regulators, and its function has been extensively studied, particularly in model species of petunia. We built a phylotranscriptomic data set for 69 species of Petunieae to infer patterns of molecular evolution across pathway genes and across lineages. We found that transcription factors exhibit faster rates of molecular evolution (dN/dS) than their targets, with the highly specialized MYB genes evolving fastest. Using the largest comparative data set to date, we recovered little support for the hypothesis that upstream enzymes evolve slower than those occupying more downstream positions, although expression levels do predict molecular evolutionary rates. Although shifts in floral pigmentation were only weakly related to changes affecting coding regions, we found a strong relationship with the presence/absence patterns of MYB transcripts. Intensely pigmented species express all three main MYB anthocyanin activators in petals, whereas pale or white species express few or none. Our findings reinforce the notion that pathway regulators have a dynamic history, involving higher rates of molecular evolution than structural components, along with frequent changes in expression during color transitions.

Independent Evolution of the MYB Family in Brown Algae

Myeloblastosis (MYB) proteins represent one of the largest families of eukaryotic transcription factors and regulate important processes in growth and development. Studies on MYBs have mainly focused on animals and plants; however, comprehensive analysis across other supergroups such as SAR (stramenopiles, alveolates, and rhizarians) is lacking. This study characterized the structure, evolution, and expression of MYBs in four brown algae, which comprise the biggest multicellular lineage of SAR. Subfamily 1R-MYB comprised heterogeneous proteins, with fewer conserved motifs found outside the MYB domain. Unlike the SHAQKY subgroup of plant 1R-MYB, THAQKY comprised the largest subgroup of brown algal 1R-MYBs. Unlike the expansion of 2R-MYBs in plants, brown algae harbored more 3R-MYBs than 2R-MYBs. At least ten 2R-MYBs, fifteen 3R-MYBs, and one 6R-MYB orthologs existed in the common ancestor of brown algae. Phylogenetic analysis showed that brown algal MYBs had ancient origins and a diverged evolution. They showed strong affinity with stramenopile species, while not with red algae, green algae, or animals, suggesting that brown algal MYBs did not come from the secondary endosymbiosis of red and green plastids. Sequence comparison among all repeats of the three types of MYB subfamilies revealed that the repeat of 1R-MYBs showed higher sequence identity with the R3 of 2R-MYBs and 3R-MYBs, which supports the idea that 1R-MYB was derived from loss of the first and second repeats of the ancestor MYB. Compared with other species of SAR, brown algal MYB proteins exhibited a higher proportion of intrinsic disordered regions, which might contribute to multicellular evolution. Expression analysis showed that many MYB genes are responsive to different stress conditions and developmental stages. The evolution and expression analyses provided a comprehensive analysis of the phylogeny and functions of MYBs in brown algae.

Genome-wide identification and genetic characterization of the CaMYB family and its response to five types of heavy metal stress in hot pepper (Capsicum annuum cv. CM334)

MYB proteins play a crucial role in plant growth and development and stress responses. In this study, 160 members of the MYB gene family from the pepper genome database were used to analyze gene structures, chromosome localization, collinearity, genetic affinity and expression in response to heavy metals. The results identified R2R3-MYB members and further phylogenetically classified them into 35 subgroups based on highly conserved gene structures and motifs. Collinearity analysis showed that segmental duplication events played a crucial role in the functional expansion of the CaMYB gene family by intraspecific collinearity, and at least 12 pairs of CaMYB genes existed between species prior to the differentiation between monocots and dicots. Moreover, the upstream CaMYB genes were mainly localized to the phytohormone elements ABRE and transcription factor elements MYB and MYC. Further analysis revealed that MYB transcription factors were closely associated with a variety of abiotic stress-related proteins (e.g., MAC-complex and SKIP). Under the stress of five metal ions, Cd²⁺, Cu²⁺, Pb²⁺, Zn²⁺, and Fe³⁺, the expression levels of some CaMYB family genes were upregulated. Of these genes, pairing homologous 1 (PH-1), PH-13, and PH-15 in the roots of Capsicum annuum were upregulated to the greatest extent, indicating that these three MYB family members are particularly sensitive to these five metals. This study provides a theoretical reference for the analysis of the molecular regulatory mechanism of MYB family genes in mediating the response to heavy metals in plants. This study reveals the mode of interaction between MYB and a variety of abiotic stress proteins and clarifies the biological functions of CaMYB family members in the regulation of heavy metal stress.

Characterization of the Liriodendron Chinense MYB Gene Family and Its Role in Abiotic Stress Response

Liriodendron chinense (Lchi) is a Magnoliaceae plant, which is a basic angiosperm left behind by the Pleistocene and mainly distributed in the south of the Yangtze River. Liriodendron hybrids has good wood properties and is widely used in furniture and in other fields. It is not clear if they can adapt to different environmental conditions, such as drought and high and low temperatures, and the molecular mechanisms for this adaptation are unknown. Among plant transcription factors (TFs), the MYB gene family is one of the largest and is often involved in stress or adversity response signaling, growth, and development. Therefore, studying the role of MYBTFs in regulating abiotic stress signaling, growth, and development in Lchi is helpful to promote afforestation in different environments. In our research, a genome-wide analysis of the LchiMYB gene family was performed, including the phylogenetic relationship tree, gene exon-intron structure, collinearity, and chromosomal position. According to the evolutionary tree, 190 LchiMYBs were divided into three main branches. LchiMYBs were evenly distributed across 19 chromosomes, with their collinearity, suggesting that segment duplication events may have contributed to LchiMYB gene expansion. Transcriptomes from eight tissues, 11 stages of somatic embryogenesis, and leaves after cold, heat, and drought stress were used to analyze the function of the MYB gene family. The results of tissue expression analysis showed that most LchiMYB genes regulated bark, leaf, bud, sepal, stigma, and stamen development, as well as the four important stages (ES3, ES4, ES9, and PL) of somatic embryogenesis. More than 60 LchiMYBs responded to heat, cold, and drought stress; some of which underwent gene duplication during evolution. LchiMYB3 was highly expressed under all three forms of stress, while LchiMYB121 was strongly induced by both cold and heat stress. Eight genes with different expression patterns were selected and verified by quantitative real-time PCR (qRT-PCR) experiments. The results suggested that these LchiMYBs may regulate Lchi growth development and resistance to abiotic stress. This study shows the cross-regulatory function of LchiMYBs in the growth and development, asexual reproduction, and abiotic resistance of Lchi. This information will prove pivotal to directing further studies on the biological function of Lchi MYBTFs in genetic improvement and abiotic stress response.

Genome-Wide Identification and Capsaicinoid Biosynthesis-Related Expression Analysis of the R2R3-MYB Gene Family in Capsicum annuum L

Capsaicinoids are naturally specialized metabolites in pepper and are the main reason that Capsicum fruits have a pungent smell. During the synthesis of capsaicin, MYB transcription factors play key regulatory roles. In particular, R2R3-MYB subfamily genes are the most important members of the MYB family and are critical candidate factors in capsaicinoid biosynthesis. The 108 R2R3-MYB genes in pepper were identified in this study and all are shown to have two highly conserved MYB binding domains. Phylogenetic and structural analyses clustered CaR2R3-MYB genes into seven groups. Interspecies collinearity analysis found that the R2R3-MYB family contains 16 duplicated gene pairs and the highest gene density is on chromosome 00 and 03. The expression levels of CaR2R3-MYB differentially expressed genes (DEGs) and capsaicinoid-biosynthetic genes (CBGs) in fruit development stages were obtained via RNA-seq and quantitative polymerase chain reaction (qRT-PCR). Co-expression analyses reveal that highly expressed CaR2R3-MYB genes are co-expressed with CBGs during early stages of pericarp and placenta development processes. It is speculated that six candidate CaR2R3-MYB genes are involved in regulating the synthesis of capsaicin and dihydrocapsaicin. This study is the first systematic analysis of the CaR2R3-MYB gene family and provided references for studying their molecular functions. At the same time, these results also laid the foundation for further research on the capsaicin characteristics of CaR2R3-MYB genes in pepper.

Convergent Evolution at the Pathway Level: Predictable Regulatory Changes during Flower Color Transitions

The predictability of evolution, or whether lineages repeatedly follow the same evolutionary trajectories during phenotypic convergence remains an open question of evolutionary biology. In this study, we investigate evolutionary convergence at the biochemical pathway level and test the predictability of evolution using floral anthocyanin pigmentation, a trait with a well-understood genetic and regulatory basis. We reconstructed the evolution of floral anthocyanin content across 28 species of the Andean clade Iochrominae (Solanaceae) and investigated how shifts in pigmentation are related to changes in expression of seven key anthocyanin pathway genes. We used phylogenetic multivariate analysis of gene expression to test for phenotypic and developmental convergence at a macroevolutionary scale. Our results show that the four independent losses of the ancestral pigment delphinidin involved convergent losses of expression of the three late pathway genes (F3'5'h, Dfr, and Ans). Transitions between pigment types affecting floral hue (e.g., blue to red) involve changes to the expression of branching genes F3'h and F3'5'h, while the expression levels of early steps of the pathway are strongly conserved in all species. These patterns support the idea that the macroevolution of floral pigmentation follows predictable evolutionary trajectories to reach convergent phenotype space, repeatedly involving regulatory changes. This is likely driven by constraints at the pathway level, such as pleiotropy and regulatory structure.

Parallel Seed Color Adaptation during Multiple Domestication Attempts of an Ancient New World Grain

Thousands of plants have been selected as crops; yet, only a few are fully domesticated. The lack of adaptation to agroecological environments of many crop plants with few characteristic domestication traits potentially has genetic causes. Here, we investigate the incomplete domestication of an ancient grain from the Americas, amaranth. Although three grain amaranth species have been cultivated as crop for millennia, all three lack key domestication traits. We sequenced 121 crop and wild individuals to investigate the genomic signature of repeated incomplete adaptation. Our analysis shows that grain amaranth has been domesticated three times from a single wild ancestor. One trait that has been selected during domestication in all three grain species is the seed color, which changed from dark seeds to white seeds. We were able to map the genetic control of the seed color adaptation to two genomic regions on chromosomes 3 and 9, employing three independent mapping populations. Within the locus on chromosome 9, we identify an MYB-like transcription factor gene, a known regulator for seed color variation in other plant species. We identify a soft selective sweep in this genomic region in one of the crop species but not in the other two species. The demographic analysis of wild and domesticated amaranths revealed a population bottleneck predating the domestication of grain amaranth. Our results indicate that a reduced level of ancestral genetic variation did not prevent the selection of traits with a simple genetic architecture but may have limited the adaptation of complex domestication traits.

A novel R3 MYB transcriptional repressor associated with the loss of floral pigmentation in Iochroma

Losses of floral pigmentation represent one of the most common evolutionary transitions in flower color, yet the genetic basis for these changes has been elucidated in only a handful of cases. Here we used crossing studies, bulk-segregant RNA sequencing, phylogenetic analyses and functional tests to identify the gene(s) responsible for the transition to white flowers in Iochroma loxense. Crosses between I. loxense and its blue-flowered sister species, I. cyaneum, suggested that a single locus controls the flower color difference and that the white allele causes a nearly complete loss of pigmentation. Examining sequence variation across phenotypic pools from the crosses, we found that alleles at a novel R3 MYB transcription factor were tightly associated with flower color variation. This gene, which we term MYBL1, falls into a class of MYB transcriptional repressors and, accordingly, higher expression of this gene is associated with downregulation of multiple anthocyanin pigment pathway genes. We confirmed the repressive function of MYBL1 through stable transformation of Nicotiana. The mechanism underlying the evolution of white flowers in I. loxense differs from that uncovered in previous studies, pointing to multiple mechanisms for achieving fixed transitions in flower color intensity.

Identification and Characterization of the Diverse Stress-Responsive R2R3-RMYB Transcription Factor from Hibiscus sabdariffa L

Various regulatory proteins play a fundamental role to manage the healthy plant growth under stress conditions. Differential display reverse transcriptase PCR and random amplification of cDNA ends (RACE) was used to explore the osmotic stress-responsive transcripts. We identified and characterized the salt stress-responsive R2R3 type RMYB transcription factor from Hibiscus sabdariffa which has an open reading frame of 690 bp, encoding 229 long chain amino acids. In silico analysis confirmed the conserved R2 and R3 domain as well as an NLS-1 localization site. The deduced amino acids of RMYB shared 83, 81, 80, 79, 72, 71, and 66% homology with Arabidopsis thaliana, Glycine max, Oryza sativa, Zea maize, Malus domestica, Populus tremula × Populus alba, and Medicago sativa specific MYB family, respectively. We observed the gene upregulation in stem, leaf, and root tissue in response to abiotic stress. Furthermore, RMYB gene was cloned into plant expression vector under CaMV35S promoter and transformed to Gossypium hirsutum: a local cotton cultivar. Overexpression of RMYB was observed in transgenic plants under abiotic stresses which further suggests its regulatory role in response to stressful conditions. The RMYB transcription factor-overexpressing in transgenic cotton plants may be used as potential agent for the development of stress tolerant crop cultivars.

Stress-Mediated cis-Element Transcription Factor Interactions Interconnecting Primary and Specialized Metabolism in planta

Plant specialized metabolites are being used worldwide as therapeutic agents against several diseases. Since the precursors for specialized metabolites come through primary metabolism, extensive investigations have been carried out to understand the detailed connection between primary and specialized metabolism at various levels. Stress regulates the expression of primary and specialized metabolism genes at the transcriptional level via transcription factors binding to specific cis-elements. The presence of varied cis-element signatures upstream to different stress-responsive genes and their transcription factor binding patterns provide a prospective molecular link among diverse metabolic pathways. The pattern of occurrence of these cis-elements (overrepresentation/common) decipher the mechanism of stress-responsive upregulation of downstream genes, simultaneously forming a molecular bridge between primary and specialized metabolisms. Though many studies have been conducted on the transcriptional regulation of stress-mediated primary or specialized metabolism genes, but not much data is available with regard to cis-element signatures and transcription factors that simultaneously modulate both pathway genes. Hence, our major focus would be to present a comprehensive analysis of the stress-mediated interconnection between primary and specialized metabolism genes via the interaction between different transcription factors and their corresponding cis-elements. In future, this study could be further utilized for the overexpression of the specific transcription factors that upregulate both primary and specialized metabolism, thereby simultaneously improving the yield and therapeutic content of plants.