GRAS transcription factors regulate cell division planes in moss overriding the default rule.

Plant cells are surrounded by a cell wall and do not migrate, which makes the regulation of cell division orientation crucial for development. Regulatory mechanisms controlling cell division orientation may have contributed to the evolution of body organization in land plants. The GRAS family of transcription factors was transferred horizontally from soil bacteria to an algal common ancestor of land plants. SHORTROOT (SHR) and SCARECROW (SCR) genes in this family regulate formative periclinal cell divisions in the roots of flowering plants, but their roles in nonflowering plants and their evolution have not been studied in relation to body organization. Here, we show that SHR cell autonomously inhibits formative periclinal cell divisions indispensable for leaf vein formation in the moss Physcomitrium patens, and SHR expression is positively and negatively regulated by SCR and the GRAS member LATERAL SUPPRESSOR, respectively. While precursor cells of a leaf vein lacking SHR usually follow the geometry rule of dividing along the division plane with the minimum surface area, SHR overrides this rule and forces cells to divide nonpericlinally. Together, these results imply that these bacterially derived GRAS transcription factors were involved in the establishment of the genetic regulatory networks modulating cell division orientation in the common ancestor of land plants and were later adapted to function in flowering plant and moss lineages for their specific body organizations.

Specific ABA-independent tomato transcriptome reprogramming under abiotic stress combination.

Crops often have to face several abiotic stresses simultaneously, and under these conditions, the plant's response significantly differs from that observed under a single stress. However, up to the present, most of the molecular markers identified for increasing plant stress tolerance have been characterized under single abiotic stresses, which explains the unexpected results found when plants are tested under real field conditions. One important regulator of the plant's responses to abiotic stresses is abscisic acid (ABA). The ABA signaling system engages many stress-responsive genes, but many others do not respond to ABA treatments. Thus, the ABA-independent pathway, which is still largely unknown, involves multiple signaling pathways and important molecular components necessary for the plant's adaptation to climate change. In the present study, ABA-deficient tomato mutants (flacca, flc) were subjected to salinity, heat, or their combination. An in-depth RNA-seq analysis revealed that the combination of salinity and heat led to a strong reprogramming of the tomato transcriptome. Thus, of the 685 genes that were specifically regulated under this combination in our flc mutants, 463 genes were regulated by ABA-independent systems. Among these genes, we identified six transcription factors (TFs) that were significantly regulated, belonging to the R2R3-MYB family. A protein-protein interaction network showed that the TFs SlMYB50 and SlMYB86 were directly involved in the upregulation of the flavonol biosynthetic pathway-related genes. One of the most novel findings of the study is the identification of the involvement of some important ABA-independent TFs in the specific plant response to abiotic stress combination. Considering that ABA levels dramatically change in response to environmental factors, the study of ABA-independent genes that are specifically regulated under stress combination may provide a remarkable tool for increasing plant resilience to climate change.

Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation.

BACKGROUND: Regulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging. RESULTS: Here we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes. We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters. CONCLUSION: Based on these associations of (linearly) distal genes' transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

Transcriptional factors targeting in cancer stem cells for tumor modulation.

Cancer Stem Cells (CSCs) are now considered the primary "seeds" for the onset, development, metastasis, and recurrence of tumors. Despite therapeutic breakthroughs, cancer remains the leading cause of death worldwide. This is because the tumor microenvironment contains a key population of cells known as CSCs, which promote tumor aggression. CSCs are self-renewing cells that aid tumor recurrence by promoting tumor growth and persisting in patients after many traditional cancer treatments. According to reports, numerous transcription factors (TF) play a key role in maintaining CSC pluripotency and its self-renewal property. The understanding of the functions, structures, and interactional dynamics of these transcription factors with DNA has modified the hypothesis, paving the way for novel transcription factor-targeted therapies. These TFs, which are crucial and are required by cancer cells, play a vital function in the etiology of human cancer. Such CSC TFs will help with gene expression profiling, which provides crucial data for predicting the prognosis of patients. To overcome anti-cancer medication resistance and completely eradicate cancer, a potent therapy combining TFs-based CSC targets with traditional chemotherapy may be developed. In order to develop therapies that could eliminate CSCs, we here concentrated on the effect of TFs and other components of signalling pathways on cancer stemness.

Identification and core gene-mining of Weighted Gene Co-expression Network Analysis-based co-expression modules related to flood resistance in quinoa seedlings.

BACKGROUND: As an emerging food crop with high nutritional value, quinoa has been favored by consumers in recent years; however, flooding, as an abiotic stress, seriously affects its growth and development. Currently, reports on the molecular mechanisms related to quinoa waterlogging stress responses are lacking; accordingly, the core genes related to these processes were explored via Weighted Gene Co-expression Network Analysis (WGCNA). RESULTS: Based on the transcriptome data, WGCNA was used to construct a co-expression network of weighted genes associated with flooding resistance-associated physiological traits and metabolites. Here, 16 closely related co-expression modules were obtained, and 10 core genes with the highest association with the target traits were mined from the two modules. Functional annotations revealed the biological processes and metabolic pathways involved in waterlogging stress, and four candidates related to flooding resistance, specifically AP2/ERF, MYB, bHLH, and WRKY-family TFs, were also identified. CONCLUSIONS: These results provide clues to the identification of core genes for quinoa underlying quinoa waterlogging stress responses. This could ultimately provide a theoretical foundation for breeding new quinoa varieties with flooding tolerance.

A consensus genome of sika deer (Cervus nippon) and transcriptome analysis provided novel insights on the regulation mechanism of transcript factor in antler development.

BACKGROUND: Sika deer (Cervus nippon) holds significance among cervids, with three genomes recently published. However, these genomes still contain hundreds of gaps and display significant discrepancies in continuity and accuracy. This poses challenges to functional genomics research and the selection of an appropriate reference genome. Thus, obtaining a high-quality reference genome is imperative to delve into functional genomics effectively. FINDINGS: Here we report a high-quality consensus genome of male sika deer. All 34 chromosomes are assembled into single-contig pseudomolecules without any gaps, which is the most complete assembly. The genome size is 2.7G with 23,284 protein-coding genes. Comparative genomics analysis found that the genomes of sika deer and red deer are highly conserved, an approximately 2.4G collinear regions with up to 99% sequence similarity. Meanwhile, we observed the fusion of red deer's Chr23 and Chr4 during evolution, forming sika deer's Chr1. Additionally, we identified 607 transcription factors (TFs) that are involved in the regulation of antler development, including RUNX2, SOX6, SOX8, SOX9, PAX8, SIX2, SIX4, SIX6, SPI1, NFAC1, KLHL8, ZN710, JDP2, and TWST2, based on this consensus reference genome. CONCLUSIONS: Our results indicated that we acquired a high-quality consensus reference genome. That provided valuable resources for understanding functional genomics. In addition, discovered the genetic basis of sika-red hybrid fertility and identified 607 significant TFs that impact antler development.

GATA family transcription factors in alga Chlamydomonas reinhardtii.

GATA family transcription factors (GATA-TFs) are metalloproteins that regulate many metabolic pathways. These conserved proteins recognize the consensus sequence (A/T)GATA(A/G) in the promoter regions of many genes and regulate their transcription in response to environmental signals. Currently, the study of GATA-TFs is of increasing interest. GATA genes and their proteins are most actively studied in vascular plants and fungi. Based on the results of numerous studies, it has been shown that GATA factors regulate the metabolic pathways of nitrogen and carbon, and also play a major role in the processes induced by light and circadian rhythms. In algae, GATA-TFs remain poorly studied, and information about them is scattered. In this work, all known data on GATA-TFs in the unicellular green alga Chlamydomonas reinhardtii has been collected and systematized. The genome of this alga contains 12 GATA coding genes. Using the phylogenetic analysis, we identified three classes of GATA factors in C. reinhardtii according to the structure of the zinc finger domain and showed their difference from the classification of GATA factors developed on vascular plants.

Integrated transcriptomic and metabolomic analyses reveals anthocyanin biosynthesis in leaf coloration of quinoa (Chenopodium quinoa Willd.).

BACKGROUND: Quinoa leaves demonstrate a diverse array of colors, offering a potential enhancement to landscape aesthetics and the development of leisure-oriented sightseeing agriculture in semi-arid regions. This study utilized integrated transcriptomic and metabolomic analyses to investigate the mechanisms underlying anthocyanin synthesis in both emerald green and pink quinoa leaves. RESULTS: Integrated transcriptomic and metabolomic analyses indicated that both flavonoid biosynthesis pathway (ko00941) and anthocyanin biosynthesis pathway (ko00942) were significantly associated with anthocyanin biosynthesis. Differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were analyzed between the two germplasms during different developmental periods. Ten DEGs were verified using qRT-PCR, and the results were consistent with those of the transcriptomic sequencing. The elevated expression of phenylalanine ammonia-lyase (PAL), chalcone synthase (CHS), 4-coumarate CoA ligase (4CL) and Hydroxycinnamoyltransferase (HCT), as well as the reduced expression of flavanone 3-hydroxylase (F3H) and Flavonol synthase (FLS), likely cause pink leaf formation. In addition, bHLH14, WRKY46, and TGA indirectly affected the activities of CHS and 4CL, collectively regulating the levels of cyanidin 3-O-(3'', 6''-O-dimalonyl) glucoside and naringenin. The diminished expression of PAL, 4CL, and HCT decreased the formation of cyanidin-3-O-(6"-O-malonyl-2"-O-glucuronyl) glucoside, leading to the emergence of emerald green leaves. Moreover, the lowered expression of TGA and WRKY46 indirectly regulated 4CL activity, serving as another important factor in maintaining the emerald green hue in leaves N1, N2, and N3. CONCLUSION: These findings establish a foundation for elucidating the molecular regulatory mechanisms governing anthocyanin biosynthesis in quinoa leaves, and also provide some theoretical basis for the development of leisure and sightseeing agriculture.

Epithelial-to-mesenchymal transition and NF-kB pathways are promoted by a mutant form of DDB2, unable to bind PCNA, in UV-damaged human cells.

BACKGROUND: DNA-Damaged Binding protein 2 (DDB2) is a protein involved in the early step of Nucleotide Excision Repair. Recently, it has been reported that DDB2 is involved in epithelial-to-mesenchymal transition (EMT), key process in tumour invasiveness and metastasis formation. However, its role is not completely known. METHODS: Boyden chamber and cell adhesion assays, and ICELLigence analysis were performed to detect HEK293 adhesion and invasion. Western blotting and gelatine zymography techniques were employed to assess the EMT protein levels and MMP enzymatic activity. Immunofluorescence analysis and pull-down assays facilitated the detection of NF-kB sub-cellular localization and interaction. RESULTS: We have previously demonstrated that the loss of DDB2-PCNA binding favours genome instability, and increases cell proliferation and motility. Here, we have investigated the phenotypic and molecular EMT-like changes after UV DNA damage, in HEK293 clones stably expressing DDB2Wt protein or a mutant form unable to interact with PCNA (DDB2PCNA-), as well as in HeLa cells transiently expressing the same DDB2 constructs. Cells expressing DDB2PCNA- showed morphological modifications along with a reduced expression of E-cadherin, an increased activity of MMP-9 and an improved ability to migrate, in concomitance with a significant upregulation of EMT-associated Transcription Factors (TFs), whose expression has been reported to favour tumour invasion. We observed a higher expression of c-Myc oncogene, NF-kB, both regulating cell proliferation and metastatic process, as well as ZEB1, a TF significantly associated with tumorigenic potential and cell migratory ability. Interestingly, a novel interaction of DDB2 with NF-kB was detected and found to be increased in cells expressing the DDB2PCNA-, suggesting a direct modulation of NF-kB by DDB2. CONCLUSION: These results highlight the role of DDB2-PCNA interaction in counteracting EMT since DDB2PCNA- protein induces in HEK293 transformed cells a gain of function contributing to the acquisition of a more aggressive phenotype.

Identification of MYB gene family in medicinal tea tree Melaleuca alternifolia (Maiden and Betche) cheel and analysis of members regulating terpene biosynthesis.

BACKGROUND: The tea tree (Melaleuca alternifolia) is renowned for its production of tea tree oil, an essential oil primarily composed of terpenes extracted from its shoot. MYB transcription factors, which are one of the largest TF families, play a crucial role in regulating primary and secondary metabolite synthesis. However, knowledge of the MYB gene family in M. alternifolia is limited. METHODS AND RESULTS: Here, we conducted a comprehensive genome-wide analysis of MYB genes in M. alternifolia, referred to as MaMYBs, including phylogenetic relationships, structures, promoter regions, and GO annotations. Our findings classified 219 MaMYBs into four subfamilies: one 5R-MYB, four 3R-MYBs, sixty-one MYB-related, and the remaining 153 are all 2R-MYBs. Seven genes (MYB189, MYB146, MYB44, MYB29, MYB175, MYB162, and MYB160) were linked to terpenoid synthesis based on GO annotation. Phylogenetic analysis with Arabidopsis homologous MYB genes suggested that MYB193 and MYB163 may also be involved in terpenoid synthesis. Additionally, through correlation analysis of gene expression and metabolite content, we identified 42 MYB genes associated with metabolite content. CONCLUSION: The results provide valuable insights into the importance of MYB transcription factors in essential oil production in M. alternifolia. These findings lay the groundwork for a better understanding of the MYB regulatory network and the development of novel strategies to enhance essential oil synthesis in M. alternifolia.

A novel Escherichia coli cell-based bioreporter for quantification of salicylic acid in cosmetics.

Transcription factor-based bioreporters have been extensively studied for monitoring and detecting environmental toxicants. In Escherichia coli, the multiple antibiotic resistance regulator (MarR) induces transcription upon binding to salicylic acid (SA). We generated SA-specific E. coli cell-based bioreporters utilizing the operator region of the mar operon and MarR as components of the reporter and sensing domains, respectively. Although bioreporters based on endogenous MarR and wild-type E. coli cells responded to SA, their sensitivity and selectivity were insufficient for practical sample monitoring. To improve these parameters, we genetically engineered host strains for optimal MarR expression, which enhanced the sensitivity of the biosensor to micromolar quantities of SA with increased selectivity. Under the optimized experimental conditions, the biosensor could quantify SA in environmental samples. For validation, the SA concentration in artificially contaminated SA-containing cosmetic samples was determined using the developed biosensor. Reliability assessment by comparing the concentrations determined using LC-MS/MS revealed > 90% accuracy of the bioreporters. Although bioreporters are not considered standard tools for environmental monitoring, bacterial cell-based bioreporters may serve as alternative tools owing to their affordability and simplicity. The SA biosensor developed in this study can potentially be a valuable tool for monitoring SA in environmental systems. KEY POINTS: • SA-responsive bioreporter is generated by employing mar operon system in E. coli • SA specificity and selectivity were enhanced by genetic/biochemical engineering • The novel bioreporter would be valuable for SA monitoring in environmental systems.

Comparative transcriptomic analysis of the nodulation-competent zone and inference of transcription regulatory network in silicon applied Glycine max [L.]-Merr. Roots.

KEY MESSAGE: The study unveils Si's regulatory influence by regulating DEGs, TFs, and TRs. Further bHLH subfamily and auxin transporter pathway elucidates the mechanisms enhancing root development and nodulation. Soybean is a globally important crop serving as a primary source of vegetable protein for millions of individuals. The roots of these plants harbour essential nitrogen fixing structures called nodules. This study investigates the multifaceted impact of silicon (Si) application on soybean, with a focus on root development, and nodulation employing comprehensive transcriptomic analyses and gene regulatory network. RNA sequence analysis was utilised to examine the change in gene expression and identify the noteworthy differentially expressed genes (DEGs) linked to the enhancement of soybean root nodulation and root development. A set of 316 genes involved in diverse biological and molecular pathways are identified, with emphasis on transcription factors (TFs) and transcriptional regulators (TRs). The study uncovers TF and TR genes, categorized into 68 distinct families, highlighting the intricate regulatory landscape influenced by Si in soybeans. Upregulated most important bHLH subfamily and the involvement of the auxin transporter pathway underscore the molecular mechanisms contributing to enhanced root development and nodulation. The study bridges insights from other research, reinforcing Si's impact on stress-response pathways and phenylpropanoid biosynthesis crucial for nodulation. The study reveals significant alterations in gene expression patterns associated with cellular component functions, root development, and nodulation in response to Si.

Identifying Core Genes Related to Low-Temperature Stress Resistance in Quinoa Seedlings Based on WGCNA.

Quinoa is a nutritious crop that is tolerant to extreme environmental conditions; however, low-temperature stress can affect quinoa growth, development, and quality. Considering the lack of molecular research on quinoa seedlings under low-temperature stress, we utilized a Weighted Gene Co-Expression Network Analysis to construct weighted gene co-expression networks associated with physiological indices and metabolites related to low-temperature stress resistance based on transcriptomic data. We screened 11 co-expression modules closely related to low-temperature stress resistance and selected 12 core genes from the two modules that showed the highest associations with the target traits. Following the functional annotation of these genes to determine the key biological processes and metabolic pathways involved in low-temperature stress, we identified four important transcription factors involved in resistance to low-temperature stress: gene-LOC110731664, gene-LOC110736639, gene-LOC110684437, and gene-LOC110720903. These results provide insights into the molecular genetic mechanism of quinoa under low-temperature stress and can be used to breed lines with tolerance to low-temperature stress.

The Transcription Factor MiMYB8 Suppresses Peel Coloration in Postharvest 'Guifei' Mango in Response to High Concentration of Exogenous Ethylene by Negatively Modulating MiPAL1.

Anthocyanin accumulation is regulated by specific genes during fruit ripening. Currently, peel coloration of mango fruit in response to exogenous ethylene and the underlying molecular mechanism remain largely unknown. The role of MiMYB8 on suppressing peel coloration in postharvest 'Guifei' mango was investigated by physiology detection, RNA-seq, qRT-PCR, bioinformatics analysis, yeast one-hybrid, dual-luciferase reporter assay, and transient overexpression. Results showed that compared with the control, low concentration of exogenous ethylene (ETH, 500 mg·L-1) significantly promoted peel coloration of mango fruit (cv. Guifei). However, a higher concentration of ETH (1000 mg·L-1) suppressed color transformation, which is associated with higher chlorophyll content, lower a* value, anthocyanin content, and phenylalanine ammonia-lyase (PAL) activity of mango fruit. M. indica myeloblastosis8 MiMYB8 and MiPAL1 were differentially expressed during storage. MiMYB8 was highly similar to those found in other plant species related to anthocyanin biosynthesis and was located in the nucleus. MiMYB8 suppressed the transcription of MiPAL1 by binding directly to its promoter. Transient overexpression of MiMYB8 in tobacco leaves and mango fruit inhibited anthocyanin accumulation by decreasing PAL activity and down-regulating the gene expression. Our observations suggest that MiMYB8 may act as repressor of anthocyanin synthesis by negatively modulating the MiPAL gene during ripening of mango fruit, which provides us with a theoretical basis for the scientific use of exogenous ethylene in practice.

Genome-wide identification and expression pattern analysis of R2R3-MYB transcription factor gene family involved in puerarin biosynthesis and response to hormone in Pueraria lobata var. thomsonii.

BACKGROUND: R2R3-MYB transcription factors regulate secondary metabolism, stress responses and development in various plants. Puerarin is a bioactive ingredient and most abundant secondary metabolite isolated from Pueraria lobata. The biosynthesis of puerarin proceeds via the phenylpropanoid pathway and isoflavonoids pathway, in which 9 key enzymes are involved. The expression of these structural genes is under control of specific PtR2R3-MYB genes in different plant tissues. However, how PtR2R3-MYB genes regulates structural genes in puerarin biosynthesis remains elusive. This study mined the PtR2R3-MYB genes involved in puerarin biosynthesis and response to hormone in Pueraria lobata var. thomsonii. RESULTS: A total of 209 PtR2R3-MYB proteins were identified, in which classified into 34 subgroups based on the phylogenetic topology and the classification of the R2R3-MYB superfamily in Arabidopsis thaliana. Furtherly physical and chemical characteristics, gene structure, and conserved motif analysis were also used to further analyze PtR2R3-MYBs. Combining puerarin content and RNA-seq data, speculated on the regulated puerarin biosynthesis of PtR2R3-MYB genes and structural genes, thus 21 PtR2R3-MYB genes and 25 structural genes were selected for validation gene expression and further explore its response to MeJA and GSH treatment by using qRT-PCR analysis technique. Correlation analysis and cis-acting element analysis revealed that 6 PtR2R3-MYB genes (PtMYB039, PtMYB057, PtMYB080, PtMYB109, PtMYB115 and PtMYB138) and 7 structural genes (PtHID2, PtHID9, PtIFS3, PtUGT069, PtUGT188, PtUGT286 and PtUGT297) were directly or indirectly regulation of puerarin biosynthesis in ZG11. It is worth noting that after MeJA and GSH treatment for 12-24 h, the expression changes of most candidate genes were consistent with the correlation of puerarin biosynthesis, which also shows that MeJA and GSH have the potential to mediate puerarin biosynthesis by regulating gene expression in ZG11. CONCLUSIONS: Overall, this study provides a comprehensive understanding of the PtR2R3-MYB and will paves the way to reveal the transcriptional regulation of puerarin biosynthesis and response to phytohormone of PtR2R3-MYB genes in Pueraria lobata var. thomsonii.

Unraveling the involvement of WRKY TFs in regulating plant disease defense signaling.

MAIN CONCLUSION: This review article explores the intricate role, regulation, and signaling mechanisms of WRKY TFs in response to biotic stress, particularly emphasizing their pivotal role in the trophism of plant-pathogen interactions. Transcription factors (TFs) play a vital role in governing both plant defense and development by controlling the expression of various downstream target genes. Early studies have shown the differential expression of certain WRKY transcription factors by microbial infections. Several transcriptome-wide studies later demonstrated that diverse sets of WRKYs are significantly activated in the early stages of viral, bacterial, and fungal infections. Furthermore, functional investigations indicated that overexpression or silencing of certain WRKY genes in plants can drastically alter disease symptoms as well as pathogen multiplication rates. Hence the new aspects of pathogen-triggered WRKY TFs mediated regulation of plant defense can be explored. The already recognized roles of WRKYs include transcriptional regulation of defense-related genes, modulation of hormonal signaling, and participation in signal transduction pathways. Some WRKYs have been shown to directly bind to pathogen effectors, acting as decoys or resistance proteins. Notably, the signaling molecules like salicylic acid, jasmonic acid, and ethylene which are associated with plant defense significantly increase the expression of several WRKYs. Moreover, induction of WRKY genes or heightened WRKY activities is also observed during ISR triggered by the beneficial microbes which protect the plants from subsequent pathogen infection. To understand the contribution of WRKY TFs towards disease resistance and their exact metabolic functions in infected plants, further studies are required. This review article explores the intrinsic transcriptional regulation, signaling mechanisms, and hormonal crosstalk governed by WRKY TFs in plant disease defense response, particularly emphasizing their specific role against different biotrophic, hemibiotrophic, and necrotrophic pathogen infections.

The pivotal role of MYB transcription factors in plant disease resistance.

MAIN CONCLUSION: MYB transcription factors are essential for diverse biology processes in plants. This review has focused on the potential molecular actions of MYB transcription factors in plant immunity. Plants possess a variety of molecules to defend against disease. Transcription factors (TFs) serve as gene connections in the regulatory networks controlling plant growth and defense against various stressors. As one of the largest TF families in plants, MYB TFs coordinate molecular players that modulate plant defense resistance. However, the molecular action of MYB TFs in plant disease resistance lacks a systematic analysis and summary. Here, we describe the structure and function of the MYB family in the plant immune response. Functional characterization revealed that MYB TFs often function either as positive or negative modulators towards different biotic stressors. Moreover, the MYB TF resistance mechanisms are diverse. The potential molecular actions of MYB TFs are being analyzed to uncover functions by controlling the expression of resistance genes, lignin/flavonoids/cuticular wax biosynthesis, polysaccharide signaling, hormone defense signaling, and the hypersensitivity response. MYB TFs have a variety of regulatory modes that fulfill pivotal roles in plant immunity. MYB TFs regulate the expression of multiple defense genes and are, therefore, important for increasing plant disease resistance and promoting agricultural production.

MIMIC: an optimization method to identify cell type-specific marker panel for cell sorting.

Multi-omics data allow us to select a small set of informative markers for the discrimination of specific cell types and study of cellular heterogeneity. However, it is often challenging to choose an optimal marker panel from the high-dimensional molecular profiles for a large amount of cell types. Here, we propose a method called Mixed Integer programming Model to Identify Cell type-specific marker panel (MIMIC). MIMIC maintains the hierarchical topology among different cell types and simultaneously maximizes the specificity of a fixed number of selected markers. MIMIC was benchmarked on the mouse ENCODE RNA-seq dataset, with 29 diverse tissues, for 43 surface markers (SMs) and 1345 transcription factors (TFs). MIMIC could select biologically meaningful markers and is robust for different accuracy criteria. It shows advantages over the standard single gene-based approaches and widely used dimensional reduction methods, such as multidimensional scaling and t-SNE, both in accuracy and in biological interpretation. Furthermore, the combination of SMs and TFs achieves better specificity than SMs or TFs alone. Applying MIMIC to a large collection of 641 RNA-seq samples covering 231 cell types identifies a panel of TFs and SMs that reveal the modularity of cell type association networks. Finally, the scalability of MIMIC is demonstrated by selecting enhancer markers from mouse ENCODE data. MIMIC is freely available at https://github.com/MengZou1/MIMIC.

Integrated miRNA-mRNA analysis reveals candidate miRNA family regulating arbuscular mycorrhizal symbiosis of Poncirus trifoliata.

Over 70% land plants live in mutualistic symbiosis with arbuscular mycorrhizal (AM) fungi, and maintenance of symbiosis requires transcriptional and post-transcriptional regulation. The former has been widely studied, whereas the latter mediated by symbiotic microRNAs (miRNAs) remains obscure, especially in woody plants. Here, we performed high-throughput sequencing of the perennial woody citrus plant Poncirus trifoliata and identified 3750 differentially expressed genes (DEGs) and 42 miRNAs (DEmiRs) upon AM fungal colonization. By analyzing cis-regulatory elements in the promoters of the DEGs, we predicted 329 key AM transcription factors (TFs). A miRNA-mRNA regulatory network was then constructed by integrating these data. Several candidate miRNA families of P. trifoliata were identified whose members target known symbiotic genes, such as miR167h-AMT2;3 and miR156e-EXO70I, or key TFs, such as miR164d-NAC and miR477a-GRAS, thus are involved in AM symbiotic processes of fungal colonization, arbuscule development, nutrient exchange and phytohormone signaling. Finally, analysis of selected miRNA family revealed that a miR159b conserved in mycorrhizal plant species and a Poncirus-specific miR477a regulate AM symbiosis. The role of miR477a was likely to target GRAS family gene RAD1 in citrus plants. Our results not only revealed that miRNA-mRNA network analysis, especially miRNA-TF analysis, is effective in identifying miRNA family regulating AM symbiosis, but also shed light on miRNA-mediated post-transcriptional regulation of AM symbiosis in woody citrus plants.

MsWRKY33 increases alfalfa (Medicago sativa L.) salt stress tolerance through altering the ROS scavenger via activating MsERF5 transcription.

Alfalfa (Medicago sativa L.) is considered to be the most important forage crop on a global scale. Nevertheless, soil salinity significantly decreases productivity, seriously threatening food security worldwide. One viable strategy is to explore salt stress-responsive factors and elucidate their underlying molecular mechanism, and utilize them in further alfalfa breeding. In the present study, we designated MsWRKY33 as a representative salt stress-responsive factor preferentially expressed in alfalfa roots and leaves. Subsequently, it was demonstrated that MsWRKY33 was localized in the cell nucleus, and functioned as a transcriptional activator of the W-box element. Transgenic alfalfa overexpressing MsWRKY33 displayed enhanced salt stress tolerance and antioxidant activities with no significant difference in other agronomic traits. Transcriptome profiling of MsWRKY33 transgenic alfalfa under control and salt treatment unveiled significantly altered expression of reactive oxygen species (ROS) scavenger genes in transgenic alfalfa. Subsequent examination revealed that MsWRKY33 binded to the promoter of MsERF5, activating its expression and consequently fine-tuning the ROS-scavenging enzyme activity. Furthermore, MsWRKY33 interacted with the functional fragment of MsCaMBP25, which participates in Ca2+ signaling transduction. Collectively, this research offers new insight into the molecular mechanism of alfalfa salt stress tolerance and highlights the potential utility of MsWRKY33 in alfalfa breeding.