Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes

Transcriptional Reactivation of Lignin Biosynthesis for the Heterologous Production of Etoposide Aglycone in Nicotiana benthamiana

Nicotiana benthamiana is a valuable plant chassis for heterologous production of medicinal plant natural products. This host is well suited for the processing of organelle-localized plant enzymes, and the conservation of the primary metabolism across the plant kingdom often provides required plant-specific precursor molecules that feed a given pathway. Despite this commonality in metabolism, limited precursor supply and/or competing host pathways can interfere with yields of heterologous products. Here, we use transient transcriptional reprogramming of endogenous N. benthamiana metabolism to drastically improve flux through the etoposide pathway derived from the medicinal plant Podophyllum spp. Specifically, coexpression of a single lignin-associated transcription factor, MYB85, with pathway genes results in unprecedented levels of heterologous product accumulation in N. benthamiana leaves: 1 mg/g dry weight (DW) of the etoposide aglycone, 35 mg/g DW (-)-deoxypodophyllotoxin, and 3.5 mg/g DW (-)-epipodophyllotoxin─up to two orders of magnitude above previously reported biosynthetic yields for the etoposide aglycone and eight times higher than what is observed for (-)-deoxypodophyllotoxin in the native medicinal plant. Unexpectedly, transient activation of lignin metabolism by transcription factor overexpression also reduces the production of undesired side products that likely result from competing N. benthamiana metabolism. Our work demonstrates that synthetic activation of lignin biosynthesis in leaf tissue is an effective strategy for optimizing the production of medicinal compounds derived from phenylpropanoid precursors in the plant chassis N. benthamiana. Furthermore, our results highlight the engineering value of MYB85, an early switch in lignin biosynthesis, for on-demand modulation of monolignol flux and support the role of MYB46 as a master regulator of lignin polymer deposition.

Trihelix Transcriptional Factor GhGT26 of Cotton Enhances Salinity Tolerance in Arabidopsis

Cotton (Gossypium hirsutum L.), the most important textile crop worldwide, often encounters abiotic stress during its growing season and its productivity is significantly limited by adverse factors. Trihelix transcription factors (also known as GT factors) are important proteins involved in the morphological development and responses to abiotic stress in plants. However, their functions and molecular mechanisms in the cotton toward abiotic stress response remain unclear. In this study, a member (GhGT26) of the cotton Trihelix family was functionally characterized in the model plant Arabidopsis. This protein containing a SANT domain belongs to the GT-1 subgroup of trihelix proteins. GhGT26 was widely expressed in tissues (with the highest level in flower) and responded to high salt and ABA treatments at the transcriptional level. Using the Arabidopsis protoplast assay system, we found that the GhGT26 protein was located in the cell nuclei. The EMSA assay revealed that the GhGT26 protein could bind to the Site1-type GT cis elements (GT-3a) and MYB elements MRE3 and MRE4. The overexpression of GhGT26 improved plant tolerance to salt stress in transgenic Arabidopsis plants. Although ABA inhibits root elongation, the statistical analysis revealed that the root lengths of GhGT26-overexpressing Arabidopsis were the same as the wild plants after ABA treatment. Our results demonstrate that GhGT26 positively regulates salt stress via ABA-independent pathways. This evidence suggests that the GhGT26 may participate in the regulation of stress tolerance in cotton.

Genome-Wide Identification of ERF Transcription Factor Family and Functional Analysis of the Drought Stress-Responsive Genes in Melilotus albus

As an important forage legume with high values in feed and medicine, Melilotus albus has been widely cultivated. The AP2/ERF transcription factor has been shown to play an important regulatory role in plant drought resistance, but it has not been reported in the legume forage crop M. albus. To digger the genes of M. albus in response to drought stress, we identified and analyzed the ERF gene family of M. albus at the genome-wide level. A total of 100 MaERF genes containing a single AP2 domain sequence were identified in this study, named MaERF001 to MaERF100, and bioinformatics analysis was performed. Collinearity analysis indicated that segmental duplication may play a key role in the expansion of the M. albus ERF gene family. Cis-acting element predictions suggest that MaERF genes are involved in various hormonal responses and abiotic stresses. The expression patterns indicated that MaERFs responded to drought stress to varying degrees. Furthermore, four up-regulated ERFs (MaERF008, MaERF037, MaERF054 and MaERF058) under drought stress were overexpressed in yeast and indicated their biological functions to confer the tolerance to drought. This work will advance the understanding of the molecular mechanisms underlying the drought response in M. albus. Further study of the promising potential candidate genes identified in this study will provide a valuable resource as the next step in functional genomics studies and improve the possibility of improving drought tolerance in M. albus by transgenic approaches.

Genome-wide identification and expression of TIFY family in cassava (Manihot esculenta Crantz)

Plant-specific TIFY [TIF(F/Y)XG] proteins serve important roles in the regulation of plant stress responses. This family encodes four subfamilies of proteins, JAZ (JASMONATE ZIM-domain), PPD (PEAPOD), ZML (Zinc-finger Inflorescence-like), and TIFY. In this work, a total of 16 JAZ, 3 PPD, 7 ZML, and 2 TIFY genes were found in cassava (Manihot esculenta Crantz) at the genome-wide level. The phylogenetics, exon-intron structure, motif organization, and conserved domains of these genes were analyzed to characterize the members of the JAZ, PPD, and ZML subfamilies. Chromosome location and synteny analyses revealed that 26 JAZ, PPD, and ZML genes were irregularly distributed across 14 of the 18 chromosomes, and 18 gene pairs were implicated in large-scale interchromosomal segmental duplication events. In addition, JAZ, PPD, and ZML gene synteny comparisons between cassava and three other plant species (Arabidopsis, Populus trichocarpa, and rice) uncovered vital information about their likely evolution. The prediction of protein interaction network and cis-acting elements reveal the function of JAZ, PPD, and ZML genes. Subsequently, expression patterns of JAZ, PPD, and ZML genes were validated by qRT-PCR as being expressed in response to osmotic, salt, and cadmium stress. Moreover, almost all JAZ subfamily genes were responsive to jasmonic acid (JA) treatment. In particular, MeJAZ1, MeJAZ13, and MeJAZ14, were highly up-regulated by three treatments, and these genes may deserve further study. This comprehensive study lays the groundwork for future research into TIFY family genes in cassava and may be valuable for genetic improvement of cassava and other related species.

Adaptive radiation in Orinus, an endemic alpine grass of the Qinghai-Tibet Plateau, based on comparative transcriptomic analysis

The species of Orinus (Poaceae) are important alpine plants with a variety of phenotypic traits and potential usages in molecular breeding toward drought-tolerant forage crops. However, the genetic basis of evolutionary adaption and diversification in the genus is still unclear. In the present study, we obtained transcriptomes for the two most divergent species, O. thoroldii and O. kokonoricus, using the Illumina platform and de novo assembly. In total, we generated 23,029 and 24,086 unigenes with N50 values of 1188 and 1203 for O. thoroldii and O. kokonoricus respectively, and identified 19,005 pairs of putative orthologs between the two species of Orinus. For these orthologs, estimations of non-synonymous/synonymous substitution rate ratios indicated that 568 pairs may be under strongly positive selection (Ka/Ks > 1), and Gene Ontogeny (GO) enrichment analysis revealed that significantly enriched pathways were in DNA repair and resistance to abiotic stress. Meanwhile, the divergence times of species between O. thoroldii and O. kokonoricus occurred 3.2 million years ago (Mya), and the recent evolutionary branch is an allotetraploid species, Cleistogenes songorica. We also detected a Ks peak of ∼0.60 for Orinus. Additionally, we identified 188 pairs of differentially expressed genes (DEGs) between the two species of Orinus, which were significantly enrich in stress resistance and lateral root development. Thus, we considered that the species diversification and evolutionary adaption of this genus was initiated by environmental selection, followed by phenotypic differentiation, finally leading to niche separation in the Qinghai-Tibet Plateau.

Genome-Wide Identification and Salt Stress Response Analysis of the bZIP Transcription Factor Family in Sugar Beet

As one of the largest transcription factor families in plants, bZIP transcription factors play important regulatory roles in different biological processes, especially in the process of stress response. Salt stress inhibits the growth and yield of sugar beet. However, bZIP-related studies in sugar beet (Beta vulgaris L.) have not been reported. This study aimed to identify the bZIP transcription factors in sugar beet and analyze their biological functions and response patterns to salt stress. Using bioinformatics, 48 BvbZIP genes were identified in the genome of sugar beet, encoding 77 proteins with large structural differences. Collinearity analysis showed that three pairs of BvbZIP genes were fragment replication genes. The BvbZIP genes were grouped according to the phylogenetic tree topology and conserved structures, and the results are consistent with those reported in Arabidopsis. Under salt stress, the expression levels of most BvbZIP genes were decreased, and only eight genes were up-regulated. GO analysis showed that the BvbZIP genes were mainly negatively regulated in stress response. Protein interaction prediction showed that the BvbZIP genes were mainly involved in light signaling and ABA signal transduction, and also played a certain role in stress responses. In this study, the structures and biological functions of the BvbZIP genes were analyzed to provide foundational data for further mechanistic studies and for facilitating the efforts toward the molecular breeding of stress-resilient sugar beet.

Characterization of B-box family genes and their expression profiles under abiotic stresses in the Melilotus albus

B-box (BBX) proteins are one of the zinc-finger transcription factor that plays a critical role in plant development, growth, and multiple stress responses. Although BBX genes have been reported in many model organisms, no comprehensive study has yet been conducted on the BBX genes in Melilotus albus, and the biological functions of this family remain unknown. In this study, a total of 20 BBX (MaBBX) genes were identified in M. albus and were phylogenetically divided into five clades. BBX members within the same clade showed similar conserved domain, suggesting similarity of potential biological function. Analysis of MaBBX conserved motifs showed that every subfamily contained two common motifs. Distribution mapping shows that BBX proteins are nonrandomly localized in eight chromosomes. The synteny showed that most homologous gene pairs of the MaBBX gene family were amplified by segmental replication, which meant segmental replication was the main way for the MaBBX gene family to evolve. Additionally, the cis-element analysis predicted light-responsive, various hormone and stress-related elements in the promoter regions of MaBBXs. Furthermore, the expression levels of all 20 MaBBX genes were detected by qRT-PCR under salt, cold, and dark stresses in M. albus. Moreover, it was observed that 16 genes had higher expression levels after 3 h of salt treatment, 10 genes were significantly upregulated after 3 h of cold treatment, and all genes were up regulated after 3 h of dark treatment, and then appeared to decline. In addition, it was also noticed that MaBBX13 may be an important candidate for improving tolerance to abiotic stress. The prediction of protein tertiary structure showed that the tertiary structures of members of the same subfamily of MaBBX proteins were highly similar. The hypothesis exhibited that most of the MaBBX proteins were predicted to be localized to the nucleus and cytoplasm and was validated by transient expression assays of MaBBX15 in tobacco leaf epidermal cells. This study provides useful information for further investigating and researching the regulatory mechanisms of BBX family genes in response to abiotic stresses in M. albus.

RNA-seq analysis reveals key genes associated with seed germination of Fritillaria taipaiensis P.Y.Li by cold stratification

Seed dormancy is an adaptive strategy for environmental evolution. However, the molecular mechanism of the breaking of seed dormancy at cold temperatures is still unclear, and the genetic regulation of germination initiated by exposure to cold temperature requires further investigation. In the initial phase of the current study, the seed coat characteristics and embryo development of Fritillaria taipaiensis P.Y.Li at different temperatures (0°C, 4°C, 10°C & 25°C) was recorded. The results obtained demonstrated that embryo elongation and the dormancy-breaking was most significantly affected at 4°C. Subsequently, transcriptome analyses of seeds in different states of dormancy, at two stratification temperatures (4°C and 25°C) was performed, combined with weighted gene coexpression network analysis (WGCNA) and metabolomics, to explore the transcriptional regulation of seed germination in F. taipaiensis at the two selected stratification temperatures. The results showed that stratification at the colder temperature (4°C) induced an up-regulation of gene expression involved in gibberellic acid (GA) and auxin biosynthesis and the down-regulation of genes related to the abscisic acid (ABA) biosynthetic pathway. Thereby promoting embryo development and the stimulation of seed germination. Collectively, these data constitute a significant advance in our understanding of the role of cold temperatures in the regulation of seed germination in F. taipaiensis and also provide valuable transcriptomic data for seed dormancy for other non-model plant species.

Molecular basis of nitrogen starvation-induced leaf senescence

Nitrogen (N), a macronutrient, is often a limiting factor in plant growth, development, and productivity. To adapt to N-deficient environments, plants have developed elaborate N starvation responses. Under N-deficient conditions, older leaves exhibit yellowing, owing to the degradation of proteins and chlorophyll pigments in chloroplasts and subsequent N remobilization from older leaves to younger leaves and developing organs to sustain plant growth and productivity. In recent years, numerous studies have been conducted on N starvation-induced leaf senescence as one of the representative plant responses to N deficiency, revealing that leaf senescence induced by N deficiency is highly complex and intricately regulated at different levels, including transcriptional, post-transcriptional, post-translational and metabolic levels, by multiple genes and proteins. This review summarizes the current knowledge of the molecular mechanisms associated with N starvation-induced leaf senescence.

Ectopic Expression of Arabidopsis thaliana zDof1.3 in Tomato (Solanum lycopersicum L.) Is Associated with Improved Greenhouse Productivity and Enhanced Carbon and Nitrogen Use

A large collection of transgenic tomato lines, each ectopically expressing a different Arabidopsis thaliana transcription factor, was screened for variants with alterations in leaf starch. Such lines may be affected in carbon partitioning, and in allocation to the sinks. We focused on ‘L4080’, which harbored an A. thaliana zDof (DNA-binding one zinc finger) isoform 1.3 (AtzDof1.3) gene, and which had a 2−4-fold higher starch-to-sucrose ratio in source leaves over the diel (p < 0.05). Our aim was to determine whether there were associated effects on productivity. L4080 plants were altered in nitrogen (N) and carbon (C) metabolism. The N-to-C ratio was higher in six-week-old L4080, and when treated with 1/10 N, L4080 growth was less inhibited compared to the wild-type and this was accompanied by faster root elongation (p < 0.05). The six-week-old L4080 acquired 42% more dry matter at 720 ppm CO2, compared to ambient CO2 (p < 0.05), while the wild-type (WT) remained unchanged. GC-MS-TOF data showed that L4080 source leaves were enriched in amino acids compared to the WT, and at 49 DPA, fruit had 25% greater mass, higher sucrose, and increased yield (25%; p < 0.05) compared to the WT. An Affymetrix cDNA array analysis suggested that only 0.39% of the 9000 cDNAs were altered by 1.5-fold (p < 0.01) in L4080 source leaves. 14C-labeling of fruit disks identified potential differences in 14-DPA fruit metabolism suggesting that post-transcriptional regulation was important. We conclude that AtzDof1.3 and the germplasm derived therefrom, should be investigated for their ‘climate-change adaptive’ potential.

Chasing Consistency: An Update of the TCP Gene Family of Malus × Domestica

The 52 members of the Teosinte-Branched 1/Cycloidea/Proliferating Cell Factors (TCP) Transcription Factor gene family in Malus × domestica (M. × domestica) were identified in 2014 on the first genome assembly, which was released in 2010. In 2017, a higher quality genome assembly for apple was released and is now considered to be the reference genome. Moreover, as in several other species, the identified TCP genes were named based on the relative position of the genes on the chromosomes. The present work consists of an update of the TCP gene family based on the latest genome assembly of M. × domestica. Compared to the previous classification, the number of TCP genes decreased from 52 to 40 as a result of the addition of three sequences and the deduction of 15. An analysis of the intragenic identity led to the identification of 15 pairs of orthologs, shedding light on the forces that shaped the evolution of this gene family. Furthermore, a revised nomenclature system is proposed that is based both on the intragenic identity and the homology with Arabidopsis thaliana (A. thaliana) TCPs in an effort to set a common standard for the TCP classification that will facilitate any future interspecific analysis.

A multifaceted comparison between the fruit-abscission and fruit-retention cultivars in ornamental crabapple

The ornamental crabapple is a multipurpose landscaping tree that bears brilliant fruit throughout the winter. However, whether or not its fruit persists after maturation is specifically correlated to cultivar characteristics. In this work, we screened two different types that display fruit-retention ("Donald Wyman," "Red Jewel," and "Sugar Tyme") and fruit-abscission ("Radiant" and "Flame") in Northern China across the whole winter using multi-year successional records. Fruit-abscission was determined predominantly by the abscission zone established at the base of the pedicel, regardless of fruit size and pedicel length, according to the results of the comparative research. The primary physiological rationale was the accumulation of hydrolases activity (pectinesterase, cellulase, polygalacturonase, and β-glucosidase). Comparative transcriptomics further identified a number of upregulated DEGs involved in the synthesis pathways of canonical phytohormones, such as ethylene, jasmonic acid, abscisic acid, and cytokinin, as well as 12 transcription factors linked in downstream signaling in fruit-abscission cultivars. Finally, a model incorporating multi-layered modulation was proposed for the fruit abscission of ornamental crabapple. This study will serve as the foundation for the development of fruit-viewing crabapples that have an extended ornamental lifetime.

Genome-wide characterization and sequence polymorphism analyses of cysteine-rich poly comb-like protein in Glycine max

Cysteine-rich poly comb-like protein (CPP) is a member of cysteine-rich transcription factors that regulates plant growth and development. In the present work, we characterized twelve CPP transcription factors encoding genes in soybean (Glycine max). Phylogenetic analyses classified CPP genes into six clades. Sequence logos analyses between G. max and G. soja amino acid residues exhibited high conservation. The presence of growth and stress-related cis-acting elements in the upstream regions of GmCPPs highlight their role in plant development and tolerance against abiotic stress. Ka/Ks levels showed that GmCPPs experienced limited selection pressure with limited functional divergence arising from segmental or whole genome duplication events. By using the PAN-genome of soybean, a single nucleotide polymorphism was identified in GmCPP-6. To perform high throughput genotyping, a kompetitive allele-specific PCR (KASP) marker was developed. Association analyses indicated that GmCPP-6-T allele of GmCPP-6 (in exon region) was associated with higher thousand seed weight under both water regimes (well-water and water-limited). Taken together, these results provide vital information to further decipher the biological functions of CPP genes in soybean molecular breeding.

Integrated transcriptomic and metabolomic analyses revealed the regulatory mechanism of sulfur application in grain yield and protein content in wheat (Triticum aestivum L.)

Sulfur fertilizers play an important role in increasing the yield and improving the dough quality of bread wheat, but their regulatory mechanism remains unclear. In this study, 0 kg·ha^-1 (S0) and 60 kg·ha^-1 (S60) of sulfur were applied on the anthesis date; subsequently, immature wheat grains at 8, 13, and 18 days post-anthesis (DPA) were subjected to integrated transcriptomic and metabolomic analyses to investigate the changes in the gene/metabolite activity in a typical strong-gluten wheat, Gaoyou2018 (GY2018). Our data show that the S60 treatment could significantly increase the grain yield and grain protein content by 13.2 and 3.6%, respectively. The transcriptomic analysis revealed that 10,694 differentially expressed genes (DEGs) were induced by S60 from 8 to 18 DPA when compared with their corresponding no-sulfur controls, and most DEGs were mainly involved in lipid metabolism and amino acid metabolism pathways. Ninety-seven MYB transcription factors (TFs) were identified as responsive to the S60 treatment; of these, 66 showed significantly differential expression at 13 DPA, and MYB118 might participate in the process of sulfur metabolism by regulating glucosinolate synthesis. In total, 542 significantly enriched differentially expressed (DE) metabolites (DEMs) were identified following the S60 treatment, which mainly included secondary metabolites, carbohydrates, and amino acids. Several metabolites (e.g., glutathione, sucrose, GDP-alpha-D-glucose, and amino acids) exhibited altered abundances following the S60 treatment. The combination of transcriptomic and metabolomic analyses highlighted the important role of amino acid metabolism (especially cysteine, methionine, and glutathione metabolism) and starch and sucrose metabolism pathways after S60 application. Our results provide valuable information enhancing our understanding of the molecular mechanism of the response to sulfur and provide useful clues for grain protein quality formation and yield improvement in bread wheat.

The membrane associated NAC transcription factors ANAC060 and ANAC040 are functionally redundant in the inhibition of seed dormancy in Arabidopsis thaliana

The NAC family of transcription factors is involved in plant development and various biotic and abiotic stresses. The Arabidopsis thaliana ANAC genes ANAC060, ANAC040, and ANAC089 are highly homologous based on protein and nucleotide sequence similarity. These three genes are predicted to be membrane bound transcription factors (MTFs) containing a conserved NAC domain, but divergent C-terminal regions. The anac060 mutant shows increased dormancy when compared with the wild type. Mutations in ANAC040 lead to higher seed germination under salt stress, and a premature stop codon in ANAC089 Cvi allele results in seeds exhibiting insensitivity to high concentrations of fructose. Thus, these three homologous MTFs confer distinct functions, although all related to germination. To investigate whether the differences in function are caused by a differential spatial or temporal regulation, or by differences in the coding sequence (CDS), we performed swapping experiments in which the promoter and CDS of the three MTFs were exchanged. Seed dormancy and salt and fructose sensitivity analyses of transgenic swapping lines in mutant backgrounds showed that there is functional redundancy between ANAC060 and ANAC040, but not between ANAC060 and ANAC089.

Genome-Wide Analysis of the ERF Family and Identification of Potential Genes Involved in Fruit Ripening in Octoploid Strawberry

Ethylene response factors (ERFs) belonging to the APETALA2/ERF superfamily acted at the end of the ethylene signaling pathway, and they were found to play important roles in plant growth and development. However, the information of ERF genes in strawberry and their involvement in fruit ripening have been limited. Here, a total of 235 ERF members were identified from 426 AP2/ERF genes at octoploid strawberry genome level and classified into six subgroups according to their sequence characteristics and phylogenetic relationship. Conserved motif and gene structure analysis supported the evolutionary conservation of FaERFs. Syntenic analysis showed that four types of duplication events occurred during the expansion of FaERF gene family. Of these, WGD/segmental duplication played a major role. Transcriptomic data of FaERF genes during fruit ripening and in response to abscisic acid screened one activator (FaERF316) and one repressor (FaERF118) that were involved in fruit ripening. Transcriptional regulation analysis showed some transcription factors related to ripening such as ABI4, TCP15, and GLK1 could bind to FaERF316 or FaERF118 promoters, while protein-protein interaction analysis displayed some proteins associated with plant growth and development could interact with FaERF118 or FaERF316. These results suggested that FaERF118 and FaERF316 were potential genes to regulate strawberry ripening. In summary, the present study provides the comprehensive and systematic information on FaERF family evolution and gains insights into FaERF's potential regulatory mechanism in strawberry ripening.

Genome-Wide Analysis of the SPL Gene Family and Expression Analysis during Flowering Induction in Prunus × yedoensis 'Somei-yoshino'

SQUAMOSA Promoter-Binding Protein-Like (SPL) genes encode plant-specific transcription factors which bind to the SQUAMOSA promoter of the MADS-box genes to regulate its expression. It plays important regulatory roles in floral induction and development, fertility, light signals and hormonal transduction, and stress response in plants. In this study, 32 PySPL genes with complete SBP (squamosa promoter binding protein) conserved domain were identified from the genome of Prunus × yedoensis ‘Somei-yoshino’ and analyzed by bioinformatics. 32 PySPLs were distributed on 13 chromosomes, encoding 32 PySPL proteins with different physical and chemical properties. The phylogenetic tree constructed with Arabidopsis thaliana and Oryza sativa can be divided into 10 subtribes, indicating PySPLs of different clusters have different biological functions. The conserved motif prediction showed that the number and distribution of motifs on each PySPL is varied. The gene structure analysis revealed that PySPLs harbored exons ranging from 2 to 10. The predictive analysis of acting elements showed that the promoter of PySPLs contain a large number of light-responsive elements, as well as response elements related to hormone response, growth and development and stress response. The analysis of the PySPLs expressions in flower induction and flower organs based on qRT-PCR showed that PySPL06/22 may be the key genes of flower development, PySPL01/06 and PySPL22 may play a role in the development of sepal and pistil, respectively. The results provide a foundation for the study of SPL transcription factors of Prunus × yedoensis ‘Somei-yoshino’ and provide more reference information of the function of SPL gene in flowering.

PIF4 enhances DNA binding of CDF2 to co-regulate target gene expression and promote Arabidopsis hypocotyl cell elongation

How specificity is conferred within gene regulatory networks is an important problem in biology. The basic helix-loop-helix PHYTOCHROME-INTERACTING FACTORs (PIFs) and single zinc-finger CYCLING DOF FACTORs (CDFs) mediate growth responses of Arabidopsis to light and temperature. We show that these two classes of transcription factor (TF) act cooperatively. CDF2 and PIF4 are temporally and spatially co-expressed, they interact to form a protein complex and act in the same genetic pathway to promote hypocotyl cell elongation. Furthermore, PIF4 substantially strengthens genome-wide occupancy of CDF2 at a subset of its target genes. One of these, YUCCA8, encodes an auxin biosynthesis enzyme whose transcription is increased by PIF4 and CDF2 to contribute to hypocotyl elongation. The binding sites of PIF4 and CDF2 in YUCCA8 are closely spaced, and in vitro PIF4 enhances binding of CDF2. We propose that this occurs by direct protein interaction and because PIF4 binding alters DNA conformation. Thus, we define mechanisms by which PIF and CDF TFs cooperate to achieve regulatory specificity and promote cell elongation in response to light.

Genome-Wide Identification and Expression Analysis of WRKY Transcription Factors in Akebiatrifoliata: A Bioinformatics Study

WRKY transcription factors have been found in most plants and play an important role in regulating organ growth and disease response. Outlining the profile of WRKY genes is a very useful project for studying morphogenesis and resistance formation. In the present study, a total of 63 WRKY genes consisting of 13 class I, 41 class II, and 9 class III genes were identified from the newly published A. trifoliata genome, of which 62 were physically distributed on all 16 chromosomes. Structurally, two AkWRKY genes (AkWRKY6 and AkWRKY52) contained four domains, and AkWRKY17 lacked the typical heptapeptide structure. Evolutionarily, 42, 16, and 5 AkWRKY genes experienced whole genome duplication (WGD) or fragmentation, dispersed duplication, and tandem duplication, respectively; 28 Ka/Ks values of 30 pairs of homologous genes were far lower than 1, while those of orthologous gene pairs between AkWRKY41 and AkWRKY52 reached up to 2.07. Transcriptome analysis showed that many of the genes were generally expressed at a low level in 12 fruit samples consisting of three tissues, including rind, flesh, and seeds, at four developmental stages, and interaction analysis between AkWRKY and AkNBS genes containing W-boxes suggested that AkWRKY24 could play a role in plant disease resistance by positively regulating AkNBS18. In summary, the WRKY gene family of A. trifoliata was systemically characterized for the first time, and the data and information obtained regarding AkWRKY could be very useful in further theoretically elucidating the molecular mechanisms of plant development and response to pathogens and practically improving favorable traits such as disease resistance.

Genome-Wide Identification and Characterization of G2-Like Transcription Factor Genes in Moso Bamboo (Phyllostachys edulis)

G2-like (GLK) transcription factors contribute significantly and extensively in regulating chloroplast growth and development in plants. This study investigated the genome-wide identification, phylogenetic relationships, conserved motifs, promoter cis-elements, MCScanX, divergence times, and expression profile analysis of PeGLK genes in moso bamboo (Phyllostachys edulis). Overall, 78 putative PeGLKs (PeGLK1–PeGLK78) were identified and divided into 13 distinct subfamilies. Each subfamily contains members displaying similar gene structure and motif composition. By synteny analysis, 42 orthologous pairs and highly conserved microsynteny between regions of GLK genes across moso bamboo and maize were found. Furthermore, an analysis of the divergence times indicated that PeGLK genes had a duplication event around 15 million years ago (MYA) and a divergence happened around 38 MYA between PeGLK and ZmGLK. Tissue-specific expression analysis showed that PeGLK genes presented distinct expression profiles in various tissues, and many members were highly expressed in leaves. Additionally, several PeGLKs were significantly up-regulated under cold stress, osmotic stress, and MeJA and GA treatment, implying that they have a likelihood of affecting abiotic stress and phytohormone responses in plants. The results of this study provide a comprehensive understanding of the moso bamboo GLK gene family, as well as elucidating the potential functional characterization of PeGLK genes.

Comparative transcriptome and metabolome analysis reveal key regulatory defense networks and genes involved in enhanced salt tolerance of Actinidia (kiwifruit)

The Actinidia (kiwifruit) is an emerging fruit plant that is severely affected by salt stress in northern China. Plants have evolved several signaling network mechanisms to cope with the detrimental effects of salt stress. To date, no reported work is available on metabolic and molecular mechanisms involved in kiwifruit salt tolerance. Therefore, the present study aims to decipher intricate adaptive responses of two contrasting salt tolerance kiwifruit species Actinidia valvata [ZMH (an important genotype), hereafter referred to as R] and Actinidia deliciosa ['Hayward' (an important green-fleshed cultivar), hereafter referred to as H] under 0.4% (w/w) salt stress for time courses of 0, 12, 24, and 72 hours (hereafter refered to as h) by combined transcriptome and metabolome analysis. Data revealed that kiwifruit displayed specific enrichment of differentially expressed genes (DEGs) under salt stress. Interestingly, roots of R plants showed a differential expression pattern for up-regulated genes. The KEGG pathway analysis revealed the enrichment of DEGs related to plant hormone signal transduction, glycine metabolism, serine and threonine metabolism, glutathione metabolism, and pyruvate metabolism in the roots of R under salt stress. The WGCNA resulted in the identification of five candidate genes related to glycine betaine (GB), pyruvate, total soluble sugars (TSS), and glutathione biosynthesis in kiwifruit. An integrated study of transcriptome and metabolome identified several genes encoding metabolites involved in pyruvate metabolism. Furthermore, several genes encoding transcription factors were mainly induced in R under salt stress. Functional validation results for overexpression of a candidate gene betaine aldehyde dehydrogenase (AvBADH, R_transcript_80484) from R showed significantly improved salt tolerance in Arabidopsis thaliana (hereafter referred to as At) and Actinidia chinensis ['Hongyang' (an important red-fleshed cultivar), hereafter referred to as Ac] transgenic plants than in WT plants. All in all, salt stress tolerance in kiwifruit roots is an intricate regulatory mechanism that consists of several genes encoding specific metabolites.

Plant-DTI: Extending the landscape of TF protein and DNA interaction in plants by a machine learning-based approach

As a sessile organism, plants hold elaborate transcriptional regulatory systems that allow them to adapt to variable surrounding environments. Current understanding of plant regulatory mechanisms is greatly constrained by limited knowledge of transcription factor (TF)-DNA interactions. To mitigate this problem, a Plant-DTI predictor (Plant DBD-TFBS Interaction) was developed here as the first machine-learning model that covered the largest experimental datasets of 30 plant TF families, including 7 plant-specific DNA binding domain (DBD) types, and their transcription factor binding sites (TFBSs). Plant-DTI introduced a novel TFBS feature construction, called TFBS base-preference, which enhanced the specificity of TFBS to DBD types. The proposed model showed better predictive performance with the TFBS base-preference than the simple binary representation. Plant-DTI was validated with 22 independent ChIP-seq datasets. It accurately predicted the measured DBD-TFBS pairs along with their TFBS motifs, and effectively predicted interactions of other TFs containing similar DBD types. Comparing to the existing state-of-art methods, Plant-DTI prediction showed a figure of merit in sensitivity and specificity with respect to the position weight matrix (PWM) and TSPTFBS methods. Finally, the proposed Plant-DTI model helped to fill the knowledge gap in the regulatory mechanisms of the cassava sucrose synthase 1 gene (MeSUS1). Plant-DTI predicted MeERF72 as a regulator of MeSUS1 in consistence with the yeast one-hybrid (Y1H) experiment. Taken together, Plant-DTI would help facilitate the prediction of TF-TFBS and TF-target gene (TG) interactions, thereby accelerating the study of transcriptional regulatory systems in plant species.

Genome-wide identification, structural analysis and expression profiles of short internodes related sequence gene family in quinoa

Based on the whole genome data information of Chenopodium quinoa Willd, the CqSRS gene family members were systematically identified and analyzed by bioinformatics methods, and the responses of CqSRS genes to NaCl (100 mmol/L), salicylic acid (200 umol/L) and low temperature (4°C) were detected by qRT-PCR. The results showed that a total of 10 SHI related sequence genes were identified in quinoa, and they were distributed on 9 chromosomes, and there were four pairs of duplicated genes. The number of amino acids encoded ranged from 143 aa to 370 aa, and the isoelectric point ranged from 4.81 to 8.90. The secondary structure was mainly composed of random coil (Cc). Most of the SRS gene encoding proteins were located in the cytoplasm (5 CqSRS). Phylogenetic analysis showed that the CqSRS genes were divided into three groups, and the gene structure showed that the number of exons of CqSRS was between two-five. Promoter analysis revealed that there are a total of 44 elements related to plant hormone response elements, light response elements, stress response elements and tissue-specific expression in the upstream regin of the gene. Protein interaction showed that all 10 CqSRS proteins appeared in the known protein interaction network diagram in Arabidopsis. Expression profile analysis showed that CqSRS genes had different expression patterns, and some genes had tissue-specific expression. qRT-PCR showed that all SRS family genes responded to ABA、NaCl、drought and low-temperature treatments, but the expression levels of different CqSRS genes were significantly different under various stresses. This study lays a foundation for further analyzed the function of CqSRS genes.

CisCross: A gene list enrichment analysis to predict upstream regulators in Arabidopsis thaliana

Having DNA-binding profiles for a sufficient number of genome-encoded transcription factors (TFs) opens up the perspectives for systematic evaluation of the upstream regulators for the gene lists. Plant Cistrome database, a large collection of TF binding profiles detected using the DAP-seq method, made it possible for Arabidopsis. Here we re-processed raw DAP-seq data with MACS2, the most popular peak caller that leads among other ones according to quality metrics. In the benchmarking study, we confirmed that the improved collection of TF binding profiles supported a more precise gene list enrichment procedure, and resulted in a more relevant ranking of potential upstream regulators. Moreover, we consistently recovered the TF binding profiles that were missing in the previous collection of DAP-seq peak sets. We developed the CisCross web service (https://plamorph.sysbio.ru/ciscross/) that gives more flexibility in the analysis of potential upstream TF regulators for Arabidopsis thaliana genes.

MYB Transcription Factors Becoming Mainstream in Plant Roots

The function of the root system is crucial for plant survival, such as anchoring plants, absorbing nutrients and water from the soil, and adapting to stress. MYB transcription factors constitute one of the largest transcription factor families in plant genomes with structural and functional diversifications. Members of this superfamily in plant development and cell differentiation, specialized metabolism, and biotic and abiotic stress processes are widely recognized, but their roles in plant roots are still not well characterized. Recent advances in functional studies remind us that MYB genes may have potentially key roles in roots. In this review, the current knowledge about the functions of MYB genes in roots was summarized, including promoting cell differentiation, regulating cell division through cell cycle, response to biotic and abiotic stresses (e.g., drought, salt stress, nutrient stress, light, gravity, and fungi), and mediate phytohormone signals. MYB genes from the same subfamily tend to regulate similar biological processes in roots in redundant but precise ways. Given their increasing known functions and wide expression profiles in roots, MYB genes are proposed as key components of the gene regulatory networks associated with distinct biological processes in roots. Further functional studies of MYB genes will provide an important basis for root regulatory mechanisms, enabling a more inclusive green revolution and sustainable agriculture to face the constant changes in climate and environmental conditions.

Comparative genomics and evolutionary analysis of plant CNGCs

Comparative genomics and computational biology offer powerful research tools for studying evolutionary mechanisms of organisms, and the identification and characterization of conserved/distant genes and gene families. The plant CNGC gene family encodes evolutionary conserved ion channel proteins involved in important signaling pathways and biological functions. The fundamental ideas and standard procedures for genome-wide identification and evolutionary analysis of plant cyclic nucleotide-gated ion channels employing various software, tools, and online servers have been discussed. In particular, this developed method focused on practical procedures involving the comparative analysis of paralogs and orthologs of CNGC genes in different plant species at different levels including phylogenetic analysis, nomenclature and classification, gene structure, molecular protein evolution, and duplication events as mechanisms of gene family expansion and synteny.

Contrasting Roles of Ethylene Response Factors in Pathogen Response and Ripening in Fleshy Fruit

Fleshy fruits are generally hard and unpalatable when unripe; however, as they mature, their quality is transformed by the complex and dynamic genetic and biochemical process of ripening, which affects all cell compartments. Ripening fruits are enriched with nutrients such as acids, sugars, vitamins, attractive volatiles and pigments and develop a pleasant taste and texture and become attractive to eat. Ripening also increases sensitivity to pathogens, and this presents a crucial problem for fruit postharvest transport and storage: how to enhance pathogen resistance while maintaining ripening quality. Fruit development and ripening involve many changes in gene expression regulated by transcription factors (TFs), some of which respond to hormones such as auxin, abscisic acid (ABA) and ethylene. Ethylene response factor (ERF) TFs regulate both fruit ripening and resistance to pathogen stresses. Different ERFs regulate fruit ripening and/or pathogen responses in both fleshy climacteric and non-climacteric fruits and function cooperatively or independently of other TFs. In this review, we summarize the current status of studies on ERFs that regulate fruit ripening and responses to infection by several fungal pathogens, including a systematic ERF transcriptome analysis of fungal grey mould infection of tomato caused by Botrytis cinerea. This deepening understanding of the function of ERFs in fruit ripening and pathogen responses may identify novel approaches for engineering transcriptional regulation to improve fruit quality and pathogen resistance.

The apple GARP family gene MdHHO3 regulates the nitrate response and leaf senescence

The regulation of plant gene expression by nitrate is a complex regulatory process. Here, we identified 90 GARP family genes in apples by genome-wide analysis. As a member of the GARP gene family, the expression of MdHHO3 (Malus domestica HYPERSENSITIVITY TO LOW PHOSPHATE-ELICITED PRIMARY ROOT SHORTENING1 HOMOLOG 3) is upregulated under N (nitrogen) supply. The results of DNA-binding site analysis and electrophoretic mobility shift assays (EMSA) showed that MdHHO3 binds to the motif-containing GAATC. Furthermore, MdHHO3 binds to its promoter sequence and inhibits its activity. In addition, the overexpression of MdHHO3 in apple calli resulted in less accumulation of nitrate in 35S:MdHHO3-GFP calli and downregulated the expression of the nitrate transport-related genes but upregulated the expression of the nitrate assimilation-related genes. Similarly, the expression of the nitrate transport-related genes was downregulated and the expression of the nitrate assimilation-related genes was upregulated in MdHHO3 overexpression Arabidopsis and tobacco plants. Interaction experiments showed that MdHHO3 could bind to the promoter MdNRT2.1 (NITRATE TRANSPORTER 2.1) and negatively regulate its expression. Moreover, the exposure of MdHHO3-overexpressing Arabidopsis and tobacco to nitrate deficiency resulted in an early senescence phenotype as compared to the WT plants. These results show that MdHHO3 can not only negatively regulate nitrate accumulation in response to nitrate but also promote early leaf senescence under nitrate deficiency. This information may be useful to further reveal the mechanism of the nitrate response and demonstrates that nitrate deficiency induces leaf senescence in apples.

Genome-Wide Identification and Expression Analysis of Dof Transcription Factors in Lotus (Nelumbo nucifera Gaertn.)

Lotus (Nelumbo nucifera Gaertn.) is a traditional Chinese aquatic flower with high ornamental and economic value, but water salinity seriously affects lotus cultivation and distribution. The Dof transcription factors (TFs) play a crucial function in the regulatory network of growth and defense in plants. However, no systematic investigations of the Dof TFs in lotus have been performed. In this study, comprehensive searches of the lotus genome yielded 29 potential NnDofs. We carried out a series of standardized analyses, which include physical properties, multiple sequence alignment, phylogenetic analysis, gene structure, motif composition, cis-acting element prediction, chromosome distribution, and synteny analysis. The results showed that segment duplication probably caused the NnDofs gene family expansion. The potential functions of NnDofs in lotus development and stress conditions are speculated by promoter analysis. Furthermore, a complete expression investigation of NnDofs utilizing an RNA-seq atlas and quantitative real-time polymerase chain reaction (qRT-PCR) was performed. The majority of the NnDofs exhibit tissue-specific expression patterns, and many genes have been identified as being extremely sensitive to salt stressors. Overall, this study is the first to report a genome-wide assessment of the Dof family in lotus, and the findings offer vital insights for prospective functional studies on lotus salinity stress.

The NAC transcription factor ANAC017 regulates aluminum tolerance by regulating the cell wall-modifying genes

Aluminum (Al) toxicity is one of the key factors limiting crop production in acid soils; however, little is known about its transcriptional regulation in plants. In this study, we characterized the role of a NAM, ATAF1/2, and cup-shaped cotyledon 2 (NAC) transcription factors (TFs), ANAC017, in the regulation of Al tolerance in Arabidopsis (Arabidopsis thaliana). ANAC017 was localized in the nucleus and exhibited constitutive expression in the root, stem, leaf, flower, and silique, although its expression and protein accumulation were repressed by Al stress. Loss of function of ANAC017 enhanced Al tolerance when compared with wild-type Col-0 and was accompanied by lower root and root cell wall Al content. Furthermore, both hemicellulose and xyloglucan content decreased in the anac017 mutants, indicating the possible interaction between ANAC017 and xyloglucan endotransglucosylase/hydrolase (XTH). Interestingly, the expression of XTH31, which is responsible for xyloglucan modification, was downregulated in the anac017 mutants regardless of Al supply, supporting the possible interaction between ANAC017 and XTH31. Yeast one-hybrid, dual-luciferase reporter assay, and chromatin immunoprecipitation-quantitative PCR analysis revealed that ANAC017 positively regulated the expression of XTH31 through directly binding to the XTH31 promoter region, and overexpression of XTH31 in the anac017 mutant background rescued its Al-tolerance phenotype. In conclusion, we identified that the tTF ANAC017 acts upstream of XTH31 to regulate Al tolerance in Arabidopsis.

Rapid sequence and functional diversification of a miRNA superfamily targeting calcium signaling components in seed plants

MicroRNA (miRNA)-directed posttranscriptional gene silencing (miR-PTGS) is an integral component of gene regulatory networks governing plant development and responses to the environment. The sequence homology between Sly-miR4376, a miRNA common to Solanaceae and reported to target autoinhibited Ca²⁺ -ATPase 10 (ACA10) messenger RNA (mRNA) in tomato, and Arabidopsis miR391 (Ath-miR391), previously annotated as a nonconserved member of the deeply conserved miR390 family, has prompted us to revisit the function of Ath-miR391, as well as its regulatory conservation. A combination of genetic, molecular, and bioinformatic analyses revealed a hidden conservation for miR-PTGS of ACA10 homologs in spermatophytes. We found that the Arabidopsis ACA10 mRNA undergoes miR391-directed cleavage in vivo. Furthermore, transgenic overexpression of miR391 recapitulated the compact inflorescence (cif) phenotypes characteristic of ACA10 loss-of-function mutants, due to miR391-directed PTGS of ACA10. Significantly, comprehensive data mining revealed robust evidence for widespread PTGS of ACA10 homologs directed by a superfamily of related miRNAs sharing a conserved sequence core. Intriguingly, the ACA-targeting miRNAs in Poaceae also direct PTGS for calmodulin-like proteins which are putative Ca²⁺ sensors. The PTGS of ACA10 homologs is therefore directed by a miRNA superfamily that is of ancient origin and has undergone rapid sequence diversification associated with functional innovation.

Intrinsically disordered protein, DNA binding with one finger transcription factor (OsDOF27) implicates thermotolerance in yeast and rice

Intrinsically disorder regions or proteins (IDRs or IDPs) constitute a large subset of the eukaryotic proteome, which challenges the protein structure-function paradigm. These IDPs lack a stable tertiary structure, yet they play a crucial role in the diverse biological process of plants. This study represents the intrinsically disordered nature of a plant-specific DNA binding with one finger transcription factor (DOF-TF). Here, we have investigated the role of OsDOF27 and characterized it as an intrinsically disordered protein. Furthermore, the molecular role of OsDOF27 in thermal stress tolerance has been elucidated. The qRT-PCR analysis revealed that OsDOF27 was significantly upregulated under different abiotic stress treatments in rice, particularly under heat stress. The stress-responsive transcript induction of OsDOF27 was further correlated with enriched abiotic stress-related cis-regulatory elements present in its promoter region. The in vivo functional analysis of the potential role of OsDOF27 in thermotolerance was further studied in yeast and in planta. Ectopic expression of OsDOF27 in yeast implicates thermotolerance response. Furthermore, the rice transgenic lines with overexpressing OsDOF27 revealed a positive role in mitigating heat stress tolerance. Collectively, our results evidently show the intrinsically disorderedness in OsDOF27 and its role in thermal stress response in rice.

Genome-Wide Identification of R2R3-MYB Transcription Factor and Expression Analysis under Abiotic Stress in Rice

The myeloblastosis (MYB) family comprises a large group of transcription factors (TFs) that has a variety of functions. Among them, the R2R3-MYB type of proteins are the largest group in plants, which are involved in controlling various biological processes such as plant growth and development, physiological metabolism, defense, and responses to abiotic and biotic stresses. In this study, bioinformatics was adopted to conduct genome-wide identification of the R2R3-MYB TFs in rice. We identified 190 MYB TFs (99 R2R3-MYBs), which are unevenly distributed on the 12 chromosomes of rice. Based on the phylogenetic clustering and protein sequence characteristics, OsMYBs were classified into five subgroups, and 59.6% of the Os2R_MYB genes contained two introns. Analysis of cis-acting elements in the 2000 bp upstream region of Os2R_MYB genes showed that all Os2R_MYB genes contained plant hormones-related or stress-responsive elements since 91.9%, 79.8%, 79.8%, and 58.6% of Os2R_MYB genes contain ABRE, TGACG, CGTCA, and MBS motifs, respectively. Protein–protein network analysis showed that the Os2R_MYBs were involved in metabolic process, biosynthetic process, and tissue development. In addition, some genes showed a tissue-specific or developmental-stage-specific expression pattern. Moreover, the transcription levels of 20 Os2R_MYB genes under polyethylene glycol (PEG) and cadmium chloride (CdCl₂) stress inducers were dissected by qRT-PCR. The results indicated genes with an altered expression upon PEG or CdCl₂ stress induction. These results potentially supply a basis for further research on the role that Os2R_MYB genes play in plant development and stress responses.

Pleiotropic ZmICE1 Is an Important Transcriptional Regulator of Maize Endosperm Starch Biosynthesis

Starch, the major component of cereal grains, affects crop yield and quality and is widely used in food and industrial applications. The biosynthesis of maize starch is a complex process involving a series of functional enzymes. However, the sophisticated regulatory mechanisms of starch biosynthetic genes have not been fully elaborated. The basic/helix-loop-helix (bHLH) transcription factors are widely distributed in eukaryotes and participate in many physiological processes. In this study, 202 bHLH encoding genes were identified in the maize genome by Blast method. ZmICE1 gene, which belongs to the ICE subfamily of the bHLH family, was obtained and expressed mainly in maize filling endosperm and co-expressed with 14 starch biosynthesis genes. Based on the comparative analyses across different plant species, we revealed that the gene structures and protein domains of the ICE subfamily were conserved between monocots and dicots, suggesting their functional conservation feature. Yeast activation and subcellular localization assays suggested that ZmICE1 had transcriptional activation activity and localized in the nucleus. Yeast one-hybrid assays confirmed that ZmICE1 could directly bind to the promoters of ZmSSIIa and ZmGBSSI. Transient gene expression analysis in maize endosperm revealed that ZmICE1 positively regulated the expression of ZmSSIIa, but inhibited the expression of ZmGBSSI. Our results indicated that ZmICE1 could function as a regulator of maize starch biosynthesis.

Homeobox 27, a Homeodomain Transcription Factor, Confers Tolerances to CMV by Associating with Cucumber Mosaic Virus 2b Protein

Transcription factors (TFs) play an important role in plant development; however, their role during viral infection largely remains unknown. The present study was designed to uncover the role transcription factors play in Cucumber mosaic virus (CMV) infection. During the screening of an Arabidopsis thaliana (Col-0) transcription factor library, using the CMV 2b protein as bait in the yeast two-hybrid system, the 2b protein interacted with Homeobox protein 27 (HB27). HB27 belongs to the zinc finger homeodomain family and is known to have a regulatory role in flower development, and responses to biotic and abiotic stress. The interaction between CMV 2b and HB27 proteins was further validated using in planta (bimolecular fluorescence complementation assay) and in vitro far-Western blotting (FWB) methods. In the bimolecular fluorescence complementation assay, these proteins reconstituted YFP fluorescence in the nucleus and the cytoplasmic region as small fluorescent dots. In FWB, positive interaction was detected using bait anti-MYC antibody on the target HB27-HA protein. During CMV infection, upregulation (~3-fold) of the HB27 transcript was observed at 14 days post-infection (dpi) in A. thaliana plants, and expression declined to the same as healthy plants at 21 dpi. To understand the role of the HB27 protein during CMV infection, virus accumulation was determined in HB27-overexpressing (HB27 OE) and knockout mutants. In HB27-overexpressing lines, infected plants developed mild symptoms, accumulating a lower virus titer at 21 dpi compared to wild-type plants. Additionally, knockout HB27 mutants had more severe symptoms and a higher viral accumulation than wild-type plants. These results indicate that HB27 plays an important role in the regulation of plant defense against plant virus infection.

The Transcription Factor MYB37 Positively Regulates Photosynthetic Inhibition and Oxidative Damage in Arabidopsis Leaves Under Salt Stress

MYB transcription factors (TFs) mediate plant responses and defenses to biotic and abiotic stresses. The effects of overexpression of MYB37, an R2R3 MYB subgroup 14 transcription factors in Arabidopsis thaliana, on chlorophyll content, chlorophyll fluorescence parameters, reactive oxygen species (ROS) metabolism, and the contents of osmotic regulatory substances were studied under 100 mM NaCl stress. Compared with the wild type (Col-0), MYB37 overexpression significantly alleviated the salt stress symptoms in A. thaliana plants. Chlorophyll a (Chl a) and chlorophyll b (Chl b) contents were significantly decreased in OE-1 and OE-2 than in Col-0. Particularly, the Chl a/b ratio was also higher in OE-1 and OE-2 than in Col-0 under NaCl stress. However, MYB37 overexpression alleviated the degradation of chlorophyll, especially Chl a. Salt stress inhibited the activities of PSII and PSI in Arabidopsis leaves, but did not affect the activity of PSII electron donor side oxygen-evolving complex (OEC). MYB37 overexpression increased photosynthesis in Arabidopsis by increasing PSII and PSI activities. MYB37 overexpression also promoted the transfer of electrons from Q _A to Q _B on the PSII receptor side of Arabidopsis under NaCl stress. Additionally, MYB37 overexpression increased Y(II) and Y(NPQ) of Arabidopsis under NaCl stress and decreased Y(NO). These results indicate that MYB37 overexpression increases PSII activity and regulates the proportion of energy dissipation in Arabidopsis leaves under NaCl stress, thus decreasing the proportion of inactivated reaction centers. Salt stress causes excess electrons and energy in the photosynthetic electron transport chain of Arabidopsis leaves, resulting in the release of reactive oxygen species (ROS), such as superoxide anion and hydrogen peroxide, leading to oxidative damage. Nevertheless, MYB37 overexpression reduced accumulation of malondialdehyde in Arabidopsis leaves under NaCl stress and alleviated the degree of membrane lipid peroxidation caused by ROS. Salt stress also enhanced the accumulation of soluble sugar (SS) and proline (Pro) in Arabidopsis leaves, thus reducing salt stress damage to plants. Salt stress also degraded soluble protein (SP). Furthermore, the accumulation of osmoregulation substances SS and Pro in OE-1 and OE-2 was not different from that in Col-0 since MYB37 overexpression in Arabidopsis OE-1, and OE-2 did not significantly affect plants under NaCl stress. However, SP content was significantly higher in OE-1 and OE-2 than in Col-0. These results indicate that MYB37 overexpression can alleviate the degradation of Arabidopsis proteins under NaCl stress, promote plant growth and improve salt tolerance.

Genome-wide identification, phylogenetic analysis, and expression profiles of trihelix transcription factor family genes in quinoa (Chenopodium quinoa Willd.) under abiotic stress conditions

Background

The trihelix family of transcription factors plays essential roles in the growth, development, and abiotic stress response of plants. Although several studies have been performed on the trihelix gene family in several dicots and monocots, this gene family is yet to be studied in Chenopodium quinoa (quinoa).

Results

In this study, 47 C. quinoa trihelix (CqTH) genes were in the quinoa genome. Phylogenetic analysis of the CqTH and trihelix genes from Arabidopsis thaliana and Beta vulgaris revealed that the genes were clustered into five subfamilies: SIP1, GTγ, GT1, GT2, and SH4. Additionally, synteny analysis revealed that the CqTH genes were located on 17 chromosomes, with the exception of chromosomes 8 and 11, and 23 pairs of segmental duplication genes were detected. Furthermore, expression patterns of 10 CqTH genes in different plant tissues and at different developmental stages under abiotic stress and phytohormone treatment were examined. Among the 10 genes, CqTH02, CqTH25, CqTH18, CqTH19, CqTH25, CqTH31, and CqTH36, were highly expressed in unripe achenes 21 d after flowering and in mature achenes compared with other plant tissues. Notably, the 10 CqTH genes were upregulated in UV-treated leaves, whereas CqTH36 was consistently upregulated in the leaves under all abiotic stress conditions.

Conclusions

The findings of this study suggest that gene duplication could be a major driver of trihelix gene evolution in quinoa. These findings could serve as a basis for future studies on the roles of CqTH transcription factors and present potential genetic markers for breeding stress-resistant and high-yielding quinoa varieties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08726-y.

Global Integrated Genomic and Transcriptomic Analyses of MYB Transcription Factor Superfamily in C3 Model Plant Oryza sativa (L.) Unravel Potential Candidates Involved in Abiotic Stress Signaling

Plant transcription factors (TFs) are significant players in transcriptional regulations, signal transduction, and constitute an integral part of signaling networks. MYB TFs are major TF superfamilies that play pivotal roles in regulation of transcriptional reprogramming, physiological processes, and abiotic stress (AbS) responses. To explore the understanding of MYB TFs, genome and transcriptome-wide identification was performed in the C3 model plant, Oryza sativa (OsMYB). This study retrieved 114 OsMYB TFs that were computationally analyzed for their expression profiling, gene organization, cis-acting elements, and physicochemical properties. Based on the microarray datasets, six OsMYB genes which were sorted out and identified by a differential expression pattern were noted in various tissues. Systematic expression profiling of OsMYB TFs showed their meta-differential expression of different AbS treatments, spatio-temporal gene expression of various tissues and their growth in the field, and gene expression profiling in responses to phytohormones. In addition, the circular ideogram of OsMYB genes in related C4 grass plants conferred the gene synteny. Protein–protein interactions of these genes revealed the molecular crosstalk of OsMYB TFs. Transcriptional analysis (qPCR) of six OsMYB players in response to drought and salinity stress suggested the involvement in individual and combined AbS responses. To decipher how these OsMYB play functional roles in AbS dynamics, further research is a prerequisite.

Genome-Wide Characterization and Analysis of the bHLH Transcription Factor Family in Suaeda aralocaspica, an Annual Halophyte With Single-Cell C4 Anatomy

Basic helix-loop-helix (bHLH) transcription factors play important roles in plant growth, development, metabolism, hormone signaling pathways, and responses to abiotic stresses. However, comprehensive genomic and functional analyses of bHLH genes have not yet been reported in desert euhalophytes. Suaeda aralocaspica, an annual C4 halophyte without Kranz anatomy, presents high photosynthetic efficiency in harsh natural habitats and is an ideal plant for identifying transcription factors involved in stress resistance. In this study, 83 bHLH genes in S. aralocaspica were identified and categorized into 21 subfamilies based on conserved motifs, gene structures, and phylogenetic analysis. Functional annotation enrichment revealed that the majority of SabHLHs were enriched in Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways involved in the response to stress conditions, as transcription factors. A number of cis-acting elements related to plant hormones and stress responses were also predicted in the promoter regions of SabHLHs, which were confirmed by expression analysis under various abiotic stress conditions (NaCl, mannitol, low temperature, ABA, GA3, MeJA, and SA); most were involved in tolerance to drought and salinity. SabHLH169 (076) protein localized in the nucleus was involved in transcriptional activity, and gene expression could be affected by different light qualities. This study is the first comprehensive analysis of the bHLH gene family in S. aralocaspica. These data will facilitate further characterization of their molecular functions in the adaptation of desert plants to abiotic stress.

UPLC-MS/MS Profile Combined With RNA-Seq Reveals the Amino Acid Metabolism in Zanthoxylum bungeanum Leaves Under Drought Stress

Zanthoxylum bungeanum leaves have a unique taste and incomparable nutritional value and hence are popular as a food item and traditional medicine in China. However, the studies on the metabolites in Z. bungeanum leaves are quite limited, especially for amino acids. Therefore, this study explored the amino acid component in Z. bungeanum leaves and also the accumulation mechanism under drought stress in two Z. bungeanum cultivars using the widely targeted metabolome combined with transcriptome analysis. A total of 56 amino acids and their derivatives were identified in Z. bungeanum leaves, including eight essential amino acids. The total amino acid content with most individual amino acids increased under progressive drought stress. More differentially accumulated amino acids (DAAs) and differentially expressed genes (DEGs) were found in FJ (Z. bungeanum cv. ‘Fengjiao’) than in HJ (Z. bungeanum cv. ‘Hanjiao’). The orthogonal projections to latent structures discriminant analysis identified nine and seven indicator DAAs in FJ and HJ leaves, respectively. The weighted gene co-expression network analysis (WGCNA) showed that the green module was significantly correlated with most indicator DAAs and revealed the important role of FBA3, DELTA-OAT, PROC, and 15 transcription factor genes in regulating the amino acid synthesis. Furthermore, the correlation analysis and redundancy analysis (RDA) identified four candidate synthesis genes (ASNS, AK, ASPS, and PK) in amino acid biosynthesis pathway. This study provided useful information for the development of Z. bungeanum leaves in food and nutrition industry and also laid the foundations for future molecular breeding.

Less is more: natural variation disrupting a miR172 gene at the di locus underlies the recessive double-flower trait in peach (P. persica L. Batsch)

BACKGROUND

With the domestication of ornamental plants, artificial selective pressure favored the propagation of mutations affecting flower shape, and double-flower varieties are now readily available for many species. In peach two distinct loci control the double-flower phenotype: the dominant Di2 locus, regulated by the deletion of the binding site for miR172 in the euAP2 PETALOSA gene Prupe.6G242400, and the recessive di locus, of which the underlying factor is still unknown.

RESULTS

Based on its genomic location a candidate gene approach was used to identify genetic variants in a diverse panel of ornamental peach accessions and uncovered three independent mutations in Prupe.2G237700, the gene encoding the transcript for microRNA miR172d: a ~5.0 Kb LTR transposable element and a ~1.2 Kb insertion both positioned upstream of the sequence encoding the pre-miR172d within the transcribed region of Prupe.2G237700, and a ~9.5 Kb deletion encompassing the whole gene sequence. qRT-PCR analysis confirmed that expression of pre-miR172d was abolished in di/di genotypes homozygous for the three variants.

CONCLUSIONS

Collectively, PETALOSA and the mutations in micro-RNA miR172d identified in this work provide a comprehensive collection of the genetic determinants at the base of the double-flower trait in the peach germplasms.

The MYB transcription factor MYB103 acts upstream of TRICHOME BIREFRINGENCE-LIKE27 in regulating aluminum sensitivity by modulating the O-acetylation level of cell wall xyloglucan in Arabidopsis thaliana

Modification of the O-acetylation level of xyloglucan (XyG) appears to affect aluminum (Al) sensitivity in Arabidopsis by modulating its binding capacity to Al. However, the transcriptional regulation of this process remains largely unknown. In our previous studies, we found that the expression of TRICHOME BIREFRINGENCE-LIKE27 (TBL27), which is responsible for the O-acetylation of XyG, was downregulated under Al stress. In the present study, we showed that the expression of an R2R3-type transcription factor-encoding gene, MYB103, was also inhibited by Al exposure and exhibited a co-expression pattern with TBL27 in roots and siliques, suggesting a potential link between MYB103 and TBL27. The loss of function of MYB103 resulted in increased Al sensitivity, as indicated by more inhibited root growth and elevated root Al content compared with the wild type. Moreover, we also detected increased Al accumulation in the root cell wall and the hemicellulose fraction, which was attributed to the changes in the O-acetylation level of XyG rather than the XyG content itself. In addition, further analysis revealed that MYB103 positively activated TBL27 expression by directly binding to the TBL27 promoter region, and TBL27 overexpression in the myb103 mutant rescued the Al-sensitive phenotype of the mutant to the wild-type level. Taken together, we conclude that MYB103 acts upstream of TBL27 to positively regulate Al resistance by modulating the O-acetylation of the cell wall XyG.

Genome-Wide Analysis of the TCP Gene Family and Their Expression Pattern Analysis in Tea Plant (Camellia sinensis)

TEOSINTE BRANCHED1/CYCLOIDEA/PCF (TCP) transcription factors TEOSINTE BRANCHED1/CYCLOIDEA/PCF have been suggested to control the cell growth and proliferation in meristems and lateral organs. A total of 37 CsTCP genes were identified and divided into two classes, class I (PCF, group 1) and class II (CIN CYC/TB1, groups 2, and 3). The residues of TEOSINTE BRANCHED1/CYCLOIDEA/PCF of Camellia sinensis (Tea plant) (CsTCP) proteins between class I and class II were definitely different in the loop, helix I, and helix II regions; however, eighteen conserved tandem was found in bHLH. There are a large number of CsTCP homologous gene pairs in three groups. Additionally, most CsTCP proteins have obvious differences in motif composition. The results illuminated that CsTCP proteins in different groups are supposed to have complementary functions, whereas those in the same class seem to display function redundancies. There is no relationship between the number of CsTCP gene members and genome size, and the CsTCP gene family has only expanded since the divergence of monocots and eudicots. WGD/segmental duplication played a vital role in the expansion of the CsTCP gene family in tea plant, and the CsTCP gene family has expanded a lot. Most CsTCP genes of group 1 are more widely and non-specifically expressed, and the CsTCP genes of group 2 are mainly expressed in buds, flowers, and leaves. Most genes of group 1 and some genes of group 2 were up-/downregulated in varying degrees under different stress, CsTCP genes of group 3 basically do not respond to stress. TCP genes involved in abiotic stress response mostly belong to PCF group. Some CsTCP genes may have the same function as the homologous genes in Arabidopsis, but there is functional differentiation.

Genome-Wide Identification and Characterization of GARP Transcription Factor Gene Family Members Reveal Their Diverse Functions in Tea Plant (Camellia sinensis)

Golden2, ARR-B, Psr1 (GARP) proteins are plant-specific transcription factors that play vital and diverse roles in plants. However, systematic research on the GARP gene family in plants, including tea plant (Camellia sinensis), is scarce. In this study, a total of 69 GARP genes were identified and characterized from the tea plant genome based on the B-motif sequence signature. The CsGARP genes were clustered into five subfamilies: PHR1/PHL1, KAN, NIGT1/HRS1/HHO, GLK and ARR-B subfamilies. The phylogenetic relationships, gene structures, chromosomal locations, conserved motifs and regulatory cis-acting elements of the CsGARP family members were comprehensively analyzed. The expansion of CsGARP genes occurred via whole-genome duplication/segmental duplication, proximal duplication, and dispersed duplication under purifying selective pressure. The expression patterns of the CsGARP genes were systematically explored from various perspectives: in different tissues during different seasons; in different leaf color stages of tea plant; under aluminum treatment and nitrogen treatment; and in response to abiotic stresses such as cold, drought and salt and to biotic stress caused by Acaphylla theae. The results demonstrate that CsGARP family genes are ubiquitously expressed and play crucial roles in the regulation of growth and development of tea plant and the responses to environmental stimuli. Collectively, these results not only provide valuable information for further functional investigations of CsGARPs in tea plant but also contribute to broadening our knowledge of the functional diversity of GARP family genes in plants.

Identification and expression analysis of GARP superfamily genes in response to nitrogen and phosphorus stress in Spirodela polyrhiza

BACKGROUND

GARP transcription factors perform critical roles in plant development and response to environmental stimulus, especially in the phosphorus (P) and nitrogen (N) sensing and uptake. Spirodela polyrhiza (giant duckweed) is widely used for phytoremediation and biomass production due to its rapid growth and efficient N and P removal capacities. However, there has not yet been a comprehensive analysis of the GRAP gene family in S. polyrhiza.

RESULTS

We conducted a comprehensive study of GRAP superfamily genes in S. polyrhiza. First, we investigated 35 SpGARP genes which have been classified into three groups based on their gene structures, conserved motifs, and phylogenetic relationship. Then, we identified the duplication events, performed the synteny analysis, and calculated the K_a/K_s ratio in these SpGARP genes. The regulatory and co-expression networks of SpGARPs were further constructed using cis-acting element analysis and weighted correlation network analysis (WGCNA). Finally, the expression pattern of SpGARP genes were analyzed using RNA-seq data and qRT-PCR, and several NIGT1 transcription factors were found to be involved in both N and P starvation responses.

CONCLUSIONS

The study provides insight into the evolution and function of GARP superfamily in S. polyrhiza, and lays the foundation for the further functional verification of SpGARP genes.

Characterization of Expression and Epigenetic Features of Core Genes in Common Wheat

The availability of multiple wheat genome sequences enables us to identify core genes and characterize their genetic and epigenetic features, thereby advancing our understanding of their biological implications within individual plant species. It is, however, largely understudied in wheat. To this end, we reanalyzed genome sequences from 16 different wheat varieties and identified 62,299 core genes. We found that core and non-core genes have different roles in subgenome differentiation. Meanwhile, according to their expression profiles, these core genes can be classified into genes related to tissue development and stress responses, including 3376 genes highly expressed in both spikelets and at high temperatures. After associating with six histone marks and open chromatin, we found that these core genes can be divided into eight sub-clusters with distinct epigenomic features. Furthermore, we found that ca. 51% of the expressed transcription factors (TFs) were marked with both H3K27me3 and H3K4me3, indicative of the bivalency feature, which can be involved in tissue development through the TF-centered regulatory network. Thus, our study provides a valuable resource for the functional characterization of core genes in stress responses and tissue development in wheat.

Genome-wide identification, phylogeny and expression analysis of the SPL gene family and its important role in salt stress in Medicago sativa L

BACKGROUND

SQUAMOSA promoter-binding protein-like (SPL) transcription factors are widely present in plants and are involved in signal transduction, the stress response and development. The SPL gene family has been characterized in several model species, such as A. thaliana and G. max. However, there is no in-depth analysis of the SPL gene family in forage, especially alfalfa (Medicago sativa L.), one of the most important forage crops worldwide.

RESULT

In total, 76 putative MsSPL genes were identified in the alfalfa genome with an uneven distribution. Based on their identity and gene structure, these MsSPLs were divided into eight phylogenetic groups. Seventy-three MsSPL gene pairs arose from segmental duplication events, and the MsSPLs on the four subgenomes of individual chromosomes displayed high collinearity with the corresponding M. truncatula genome. The prediction of the cis-elements in the promoter regions of the MsSPLs detected two copies of ABA (abscisic acid)-responsive elements (ABREs) on average, implying their potential involvement in alfalfa adaptation to adverse environments. The transcriptome sequencing of MsSPLs in roots and leaves revealed that 54 MsSPLs were expressed in both tissues. Upon salt treatment, three MsSPLs (MsSPL17, MsSPL23 and MsSPL36) were significantly regulated, and the transcription level of MsSPL36 in leaves was repressed to 46.6% of the control level.

CONCLUSION

In this study, based on sequence homology, we identified 76 SPL genes in the alfalfa. The SPLs with high identity shared similar gene structures and motifs. In total, 71.1% (54 of 76) of the MsSPLs were expressed in both roots and leaves, and the majority (74.1%) preferred underground tissues to aerial tissues. MsSPL36 in leaves was significantly repressed under salt stress. These findings provide comprehensive information regarding the SPB-box gene family for improve alfalfa tolerance to high salinity.

Genome-wide identification and expression analysis of the R2R3-MYB gene family in tobacco (Nicotiana tabacum L.)

BACKGROUND

The R2R3-MYB transcription factor is one of the largest gene families in plants and involved in the regulation of plant development, hormone signal transduction, biotic and abiotic stresses. Tobacco is one of the most important model plants. Therefore, it will be of great significance to investigate the R2R3-MYB gene family and their expression patterns under abiotic stress and senescence in tobacco.

RESULTS

A total of 174 R2R3-MYB genes were identified from tobacco (Nicotiana tabacum L.) genome and were divided into 24 subgroups based on phylogenetic analysis. Gene structure (exon/intron) and protein motifs were especially conserved among the NtR2R3-MYB genes, especially members within the same subgroup. The NtR2R3-MYB genes were distributed on 24 tobacco chromosomes. Analysis of gene duplication events obtained 3 pairs of tandem duplication genes and 62 pairs of segmental duplication genes, suggesting that segmental duplications is the major pattern for R2R3-MYB gene family expansion in tobacco. Cis-regulatory elements of the NtR2R3-MYB promoters were involved in cellular development, phytohormones, environmental stress and photoresponsive. Expression profile analysis showed that NtR2R3-MYB genes were widely expressed in different maturity tobacco leaves, and however, the expression patterns of different members appeared to be diverse. The qRT-PCR analysis of 15 NtR2R3-MYBs confirmed their differential expression under different abiotic stresses (cold, salt and drought), and notably, NtMYB46 was significantly up-regulated under three treatments.

CONCLUSIONS

In summary, a genome-wide identification, evolutionary and expression analysis of R2R3-MYB gene family in tobacco were conducted. Our results provided a solid foundation for further biological functional study of NtR2R3-MYB genes in tobacco.

Identification of TCP family in moso bamboo (Phyllostachys edulis) and salt tolerance analysis of PheTCP9 in transgenic Arabidopsis

Bioinformatic analysis of moso bamboo TEOSINTE BRANCHED 1, CYCLOIDEA, and PROLIFERATING CELL FACTORS (TCP) transcription factors reveals their conservation and variation as well as the probable biological functions in abiotic stress response. Overexpressing PheTCP9 in Arabidopsis thaliana illustrates it may exhibit a new vision in different aspects of response to salt stress. Plant specific TCPs play important roles in plant growth, development and stress response, but studies of TCP in moso bamboo are limited. Therefore, in this study, a total of 40 TCP genes (PheTCP1 ~ 40) were identified and characterized from moso bamboo genome and divided into three different subfamilies, namely, 7 in TEOSINTE BRANCHED 1 / CYCLOIDEA (TB1/CYC), 14 in CINCINNATA (CIN) and 19 in PROLIFERATING CELL FACTOR (PCF). Subsequently, we analyzed the gene structures and conserved domain of these genes and found that the members from the same subfamilies exhibited similar exon/intron distribution patterns. Selection pressure and gene duplication analysis results indicated that PheTCP genes underwent strong purification selection during evolution. There were many cis-elements related to phytohermone and stress responsive existing in the upstream promoter regions of PheTCP genes, such as ABRE, CGTCA-motif and ARE. Subcellular localization experiments showed that PheTCP9 was a nuclear localized protein. As shown by β-glucuronidase (GUS) activity, the promoter of PheTCP9 was significantly indicated by salt stress. PheTCP9 was significantly induced in the roots, stems and leaves of moso bamboo. It was also significantly induced by NaCl solution. Overexpressing PheTCP9 increased the salt tolerance of transgenic Arabidopsis. Meanwhile, H₂O₂ and malondialdehyde (MDA) contents were significantly lower in PheTCP9 over expression (OE) transgenic Arabidopsis than WT. Catalase (CAT) activity, K⁺/Na⁺ ratio as well as CAT2 expression level was also much improved in transgenic Arabidopsis than WT under salt conditions. In addition, PheTCP9 OE transgenic Arabidopsis held higher survival rates of seedlings than WT under NaCl conditions. These results showed the positive regulation functions of PheTCP9 in plants under salt conditions.

Systematic analysis and expression profiles of TCP gene family in Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn.) revealed the potential function of FtTCP15 and FtTCP18 in response to abiotic stress

Background

As transcription factors, the TCP genes are considered to be promising targets for crop enhancement for their responses to abiotic stresses. However, information on the systematic characterization and functional expression profiles under abiotic stress of TCPs in Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn.) is limited.

Results

In this study, we identified 26 FtTCPs and named them according to their position on the chromosomes. Phylogenetic tree, gene structure, duplication events, and cis-acting elements were further studied and syntenic analysis was conducted to explore the bioinformatic traits of the FtTCP gene family. Subsequently, 12 FtTCP genes were selected for expression analysis under cold, dark, heat, salt, UV, and waterlogging (WL) treatments by qRT-PCR. The spatio-temporal specificity, correlation analysis of gene expression levels and interaction network prediction revealed the potential function of FtTCP15 and FtTCP18 in response to abiotic stresses. Moreover, subcellular localization confirmed that FtTCP15 and FtTCP18 localized in the nucleus function as transcription factors.

Conclusions

In this research, 26 TCP genes were identified in Tartary buckwheat, and their structures and functions have been systematically explored. Our results reveal that the FtTCP15 and FtTCP18 have special cis-elements in response to abiotic stress and conserved nature in evolution, indicating they could be promising candidates for further functional verification under multiple abiotic stresses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08618-1.