The MEME Suite

Comprehensive Comparative Analysis of the JAZ Gene Family in Common Wheat (Triticum aestivum) and Its D-Subgenome Donor Aegilops tauschii

JASMONATE-ZIM DOMAIN (JAZ) repressor proteins work as co-receptors in the jasmonic acid (JA) signalling pathway and are essential for plant development and environmental adaptation. Despite wheat being one of the main staple food crops, until recently, comprehensive analysis of its JAZ gene family has been limited due to the lack of complete and high-quality reference genomes. Here, using the latest reference genome, we identified 17 JAZ genes in the wheat D-genome donor Aegilops tauschii. Then, 54 TaJAZs were identified in common wheat. A systematic examination of the gene structures, conserved protein domains, and phylogenetic relationships of this gene family was performed. Five new JAZ genes were identified as being derived from tandem duplication after wheat divergence from other species. We integrated RNA-seq data and yield QTL information and found that tandemly duplicated TaJAZ genes were prone to association with spike-related traits. Moreover, 12 TaJAZ genes were located within breeding selection sweeps, including 9 tandemly duplicated ones. Haplotype variation analysis of selected JAZ genes showed significant association of TaJAZ7A and TaJAZ13A with thousand-grain weight. Our work provides a clearer picture of wheat JAZ gene evolution and puts forward the possibility of using these genes for wheat yield improvement.

Loop Catalog: a comprehensive HiChIP database of human and mouse samples

HiChIP enables cost-effective and high-resolution profiling of regulatory and structural loops. To leverage the increasing number of publicly available HiChIP datasets from diverse cell lines and primary cells, we developed the Loop Catalog (https://loopcatalog.lji.org), a web-based database featuring HiChIP loop calls for 1319 samples across 133 studies and 44 high-resolution Hi-C loop calls. We demonstrate its utility in interpreting fine-mapped GWAS variants (SNP-to-gene linking), in identifying enriched sequence motifs and motif pairs at loop anchors, and in network-level analysis of loops connecting regulatory elements (community detection). Our comprehensive catalog, spanning over 4M unique 5kb loops, along with the accompanying analysis modalities constitutes an important resource for studies in gene regulation and genome organization.

Medicago truncatula β-glucosidase 17 contributes to drought and salt tolerance through antioxidant flavonoid accumulation

Flavonoids are usually present in forms of glucosides in plants, which could be catabolized by β-glucosidase (BGLU) to form their corresponding flavonoid aglycones. In this study, we isolated three abiotic-responsive BGLU genes (MtBGLU17, MtBGLU21 and MtBGLU22) from Medicago truncatula, and found only the recombinant MtBGLU17 protein could catalyse the hydrolysis of flavonoid glycosides. The recombinant MtBGLU17 protein is active towards a variety of flavonoid glucosides, including glucosides of flavones (apigenin and luteolin), flavonols (kaempferol and quercetin), isoflavones (genistein and daidzein) and flavanone (naringenin). In particular, the recombinant MtBGLU17 protein preferentially hydrolyses flavonoid-7-O-glucosides over their corresponding 3-O-glucosides. The content of luteoin-7-O-glucoside was reduced in the MtBGLU17 overexpression plants but increased in the Tnt-1 insertional mutant lines, whereas luteoin content was increased in the MtBGLU17 overexpression plants but reduced in the Tnt-1 insertional mutant lines. Under drought and salt (NaCl) treatment, the MtBGLU17 overexpression lines showed relatively higher DPPH content, and higher CAT and SOD activity than the wild type control. These results indicated that overexpression lines of MtBGLU17 possess higher antioxidant activity and thus confer drought and salt tolerance, implying MtBGLU17 could be potentially used as a candidate gene to improve plant abiotic stress tolerance.

Epigenetic control of adaptive or homeostatic splicing during interval-training activities

Interval-training activities induce adaptive cellular changes without altering their fundamental identity, but the precise underlying molecular mechanisms are not fully understood. In this study, we demonstrate that interval-training depolarization (ITD) of pituitary cells triggers distinct adaptive or homeostatic splicing responses of alternative exons. This occurs while preserving the steady-state expression of the Prolactin and other hormone genes. The nature of these splicing responses depends on the exon's DNA methylation status, the methyl-C-binding protein MeCP2 and its associated CA-rich motif-binding hnRNP L. Interestingly, the steady expression of the Prolactin gene is also reliant on MeCP2, whose disruption leads to exacerbated multi-exon aberrant splicing and overexpression of the hormone gene transcripts upon ITD, similar to the observed hyperprolactinemia or activity-dependent aberrant splicing in Rett Syndrome. Therefore, epigenetic control is crucial for both adaptive and homeostatic splicing and particularly the steady expression of the Prolactin hormone gene during ITD. Disruption in this regulation may have significant implications for the development of progressive diseases.

EZH2 mutations in follicular lymphoma distort H3K27me3 profiles and alter transcriptional responses to PRC2 inhibition

Mutations in chromatin regulators are widespread in cancer. Among them, the histone H3 lysine 27 methyltransferase Polycomb Repressive Complex 2 (PRC2) shows distinct alterations according to tumor type. This specificity is poorly understood. Here, we model several PRC2 alterations in one isogenic system to reveal their comparative effects. Focusing then on lymphoma-associated EZH2 mutations, we show that Ezh2Y641F induces aberrant H3K27 methylation patterns even without wild-type Ezh2, which are alleviated by partial PRC2 inhibition. Remarkably, Ezh2Y641F rewires the response to PRC2 inhibition, leading to induction of antigen presentation genes. Using a unique longitudinal follicular lymphoma cohort, we further link EZH2 status to abnormal H3K27 methylation. We also uncover unexpected variability in the mutational landscape of successive biopsies, pointing to frequent co-existence of different clones and cautioning against stratifying patients based on single sampling. Our results clarify how oncogenic PRC2 mutations disrupt chromatin and transcription, and the therapeutic vulnerabilities this creates.

RAVER1 hinders lethal EMT and modulates miR/RISC activity by the control of alternative splicing

The RAVER1 protein serves as a co-factor in guiding the polypyrimidine tract-binding protein (PTBP)-dependent control of alternative splicing (AS). Whether RAVER1 solely acts in concert with PTBPs and how it affects cancer cell fate remained elusive. Here, we provide the first comprehensive investigation of RAVER1-controlled AS in cancer cell models. This reveals a pro-oncogenic role of RAVER1 in modulating tumor growth and epithelial-mesenchymal-transition (EMT). Splicing analyses and protein-association studies indicate that RAVER1 guides AS in association with other splicing regulators, including PTBPs and SRSFs. In cancer cells, one major function of RAVER1 is the stimulation of proliferation and restriction of apoptosis. This involves the modulation of AS events within the miR/RISC pathway. Disturbance of RAVER1 impairs miR/RISC activity resulting in severely deregulated gene expression, which promotes lethal TGFB-driven EMT. Among others, RAVER1-modulated splicing events affect the insertion of protein interaction modules in factors guiding miR/RISC-dependent gene silencing. Most prominently, in all three human TNRC6 proteins, RAVER1 controls AS of GW-enriched motifs, which are essential for AGO2-binding and the formation of active miR/RISC complexes. We propose, that RAVER1 is a key modulator of AS events in the miR/RISC pathway ensuring proper abundance and composition of miR/RISC effectors. This ensures balanced expression of TGFB signaling effectors and limits TGFB induced lethal EMT.

Determinant of m6A regional preference by transcriptional dynamics

N6-Methyladenosine (m6A) is the most abundant chemical modification occurring on eukaryotic mRNAs, and has been reported to be involved in almost all stages of mRNA metabolism. The distribution of m6A sites is notably asymmetric along mRNAs, with a strong preference toward the 3' terminus of the transcript. How m6A regional preference is shaped remains incompletely understood. In this study, by performing m6A-seq on chromatin-associated RNAs, we found that m6A regional preference arises during transcription. Nucleosome occupancy is remarkedly increased in the region downstream of m6A sites, suggesting an intricate interplay between m6A methylation and nucleosome-mediated transcriptional dynamics. Notably, we found a remarkable slowdown of Pol-II movement around m6A sites. In addition, inhibiting Pol-II movement increases nearby m6A methylation levels. By analyzing massively parallel assays for m6A, we found that RNA secondary structures inhibit m6A methylation. Remarkably, the m6A sites associated with Pol-II pausing tend to be embedded within RNA secondary structures. These results suggest that Pol-II pausing could affect the accessibility of m6A motifs to the methyltransferase complex and subsequent m6A methylation by mediating RNA secondary structure. Overall, our study reveals a crucial role of transcriptional dynamics in the formation of m6A regional preference.

Functional differentiation of olive PLP_deC genes: insights into metabolite biosynthesis and genetic improvement at the whole-genome level

KEY MESSAGE

This study identified 16 pyridoxal phosphate-dependent decarboxylases in olive at the whole-genome level, conducted analyses on their physicochemical properties, evolutionary relationships and characterized their activity. Group II pyridoxal phosphate-dependent decarboxylases (PLP_deC II) mediate the biosynthesis of characteristic olive metabolites, such as oleuropein and hydroxytyrosol. However, there have been no report on the functional differentiation of this gene family at the whole-genome level. This study conducted an exploration of the family members of PLP_deC II at the whole-genome level, identified 16 PLP_deC II genes, and analyzed their gene structure, physicochemical properties, cis-acting elements, phylogenetic evolution, and gene expression patterns. Prokaryotic expression and enzyme activity assays revealed that OeAAD2 and OeAAD4 could catalyze the decarboxylation reaction of tyrosine and dopa, resulting in the formation of their respective amine compounds, but it did not catalyze phenylalanine and tryptophan. Which is an important step in the synthetic pathway of hydroxytyrosol and oleuropein. This finding established the foundational data at the molecular level for studying the functional aspects of the olive PLP_deC II gene family and provided essential gene information for genetic improvement of olive.

Computational tools for plant genomics and breeding

Plant genomics and crop breeding are at the intersection of biotechnology and information technology. Driven by a combination of high-throughput sequencing, molecular biology and data science, great advances have been made in omics technologies at every step along the central dogma, especially in genome assembling, genome annotation, epigenomic profiling, and transcriptome profiling. These advances further revolutionized three directions of development. One is genetic dissection of complex traits in crops, along with genomic prediction and selection. The second is comparative genomics and evolution, which open up new opportunities to depict the evolutionary constraints of biological sequences for deleterious variant discovery. The third direction is the development of deep learning approaches for the rational design of biological sequences, especially proteins, for synthetic biology. All three directions of development serve as the foundation for a new era of crop breeding where agronomic traits are enhanced by genome design.

Human brain small extracellular vesicles contain selectively packaged, full-length mRNA

Brain cells release and take up small extracellular vesicles (sEVs) containing bioactive nucleic acids. sEV exchange is hypothesized to contribute to stereotyped spread of neuropathological changes in the diseased brain. We assess mRNA from sEVs of postmortem brain from non-diseased (ND) individuals and those with Alzheimer's disease (AD) using short- and long-read sequencing. sEV transcriptomes are distinct from those of bulk tissue, showing enrichment for genes including mRNAs encoding ribosomal proteins and transposable elements such as human-specific LINE-1 (L1Hs). AD versus ND sEVs show enrichment of inflammation-related mRNAs and depletion of synaptic signaling mRNAs. sEV mRNAs from cultured murine primary neurons, astrocytes, or microglia show similarities to human brain sEVs and reveal cell-type-specific packaging. Approximately 80% of neural sEV transcripts sequenced using long-read sequencing are full length. Motif analyses of sEV-enriched isoforms elucidate RNA-binding proteins that may be associated with sEV loading. Collectively, we show that mRNA in brain sEVs is intact, selectively packaged, and altered in disease.

Depletion of Mettl3 in cholinergic neurons causes adult-onset neuromuscular degeneration

Motor neuron (MN) demise is a hallmark of several neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). Post-transcriptional gene regulation can control RNA's fate, and defects in RNA processing are critical determinants of MN degeneration. N6-methyladenosine (m6A) is a post-transcriptional RNA modification that controls diverse aspects of RNA metabolism. To assess the m6A requirement in MNs, we depleted the m6A methyltransferase-like 3 (METTL3) in cells and mice. METTL3 depletion in embryonic stem cell-derived MNs has profound and selective effects on survival and neurite outgrowth. Mice with cholinergic neuron-specific METTL3 depletion display a progressive decline in motor behavior, accompanied by MN loss and muscle denervation, culminating in paralysis and death. Reader proteins convey m6A effects, and their silencing phenocopies METTL3 depletion. Among the m6A targets, we identified transactive response DNA-binding protein 43 (TDP-43) and discovered that its expression is under epitranscriptomic control. Thus, impaired m6A signaling disrupts MN homeostasis and triggers neurodegeneration conceivably through TDP-43 deregulation.

Anaerostipes hadrus, a butyrate-producing bacterium capable of metabolizing 5-fluorouracil

Anaerostipes hadrus (A. hadrus) is a dominant species in the human gut microbiota and considered a beneficial bacterium for producing probiotic butyrate. However, recent studies have suggested that A. hadrus may negatively affect the host through synthesizing fatty acid and metabolizing the anticancer drug 5-fluorouracil, indicating that the impact of A. hadrus is complex and unclear. Therefore, comprehensive genomic studies on A. hadrus need to be performed. We integrated 527 high-quality public A. hadrus genomes and five distinct metagenomic cohorts. We analyzed these data using the approaches of comparative genomics, metagenomics, and protein structure prediction. We also performed validations with culture-based in vitro assays. We constructed the first large-scale pan-genome of A. hadrus (n = 527) and identified 5-fluorouracil metabolism genes as ubiquitous in A. hadrus genomes as butyrate-producing genes. Metagenomic analysis revealed the wide and stable distribution of A. hadrus in healthy individuals, patients with inflammatory bowel disease, and patients with colorectal cancer, with healthy individuals carrying more A. hadrus. The predicted high-quality protein structure indicated that A. hadrus might metabolize 5-fluorouracil by producing bacterial dihydropyrimidine dehydrogenase (encoded by the preTA operon). Through in vitro assays, we validated the short-chain fatty acid production and 5-fluorouracil metabolism abilities of A. hadrus. We observed for the first time that A. hadrus can convert 5-fluorouracil to α-fluoro-β-ureidopropionic acid, which may result from the combined action of the preTA operon and adjacent hydA (encoding bacterial dihydropyrimidinase). Our results offer novel understandings of A. hadrus, exceptionally functional features, and potential applications.

IMPORTANCE

This work provides new insights into the evolutionary relationships, functional characteristics, prevalence, and potential applications of Anaerostipes hadrus.

Deciphering the differential expression patterns of yield-related negative regulators in hexaploid wheat cultivars and hybrids at different growth stages

BACKGROUND

Hexaploid bread wheat underwent a series of polyploidization events through interspecific hybridizations that conferred adaptive plasticity and resulted in duplication and neofunctionalization of major agronomic genes. The genetic architecture of polyploid wheat not only confers adaptive plasticity but also offers huge genetic diversity. However, the contribution of different gene copies (homeologs) encoded from different subgenomes (A, B, D) at different growth stages remained unexplored.

METHODS

In this study, hybrid of elite cultivars of wheat were developed via reciprocal crosses (cytoplasm swapping) and phenotypically evaluated. We assessed differential expression profiles of yield-related negative regulators in these cultivars and their F1 hybrids and identified various cis-regulatory signatures by employing bioinformatics tools. Furthermore, the preferential expression patterns of the syntenic triads encoded from A, B, and D subgenomes were assessed to decipher their functional redundancy at six different growth stages.

RESULTS

Hybrid progenies showed better heterosis such as up to 17% increase in the average number of grains and up to 50% increase in average thousand grains weight as compared to mid-parents. Based on the expression profiling, our results indicated significant dynamic transcriptional expression patterns, portraying the different homeolog-dominance at the same stage in the different cultivars and their hybrids. Albeit belonging to same syntenic triads, a dynamic trend was observed in the regulatory signatures of these genes that might be influencing their expression profiles.

CONCLUSION

These findings can substantially contribute and provide insights for the selective introduction of better cultivars into traditional and hybrid breeding programs which can be harnessed for the improvement of future wheat.

Transcriptome-informed identification and characterization of Planococcus citri cis- and trans-isoprenyl diphosphate synthase genes

Insect physiology and reproduction depend on several terpenoid compounds, whose biosynthesis is mainly unknown. One enigmatic group of insect monoterpenoids are mealybug sex pheromones, presumably resulting from the irregular coupling activity of unidentified isoprenyl diphosphate synthases (IDSs). Here, we performed a comprehensive search for IDS coding sequences of the pest mealybug Planococcus citri. We queried the available genomic and newly generated short- and long-read P. citri transcriptomic data and identified 18 putative IDS genes, whose phylogenetic analysis indicates several gene family expansion events. In vitro testing confirmed regular short-chain coupling activity with five gene products. With the candidate with highest IDS activity, we also detected low amounts of irregular coupling products, and determined amino acid residues important for chain-length preference and irregular coupling activity. This work therefore provides an important foundation for deciphering terpenoid biosynthesis in mealybugs, including the sex pheromone biosynthesis in P. citri.

Genome-wide identification and stress response analysis of BcaCPK gene family in amphidiploid Brassica carinata

BACKGROUND

Calcium-dependent protein kinases (CPKs) are crucial for recognizing and transmitting Ca2+ signals in plant cells, playing a vital role in growth, development, and stress response. This study aimed to identify and detect the potential roles of the CPK gene family in the amphidiploid Brassica carinata (BBCC, 2n = 34) using bioinformatics methods.

RESULTS

Based on the published genomic information of B. carinata, a total of 123 CPK genes were identified, comprising 70 CPK genes on the B subgenome and 53 on the C subgenome. To further investigate the homologous evolutionary relationship between B. carinata and other plants, the phylogenetic tree was constructed using CPKs in B. carinata and Arabidopsis thaliana. The phylogenetic analysis classified 123 family members into four subfamilies, where gene members within the same subfamily exhibited similar conserved motifs. Each BcaCPK member possesses a core protein kinase domain and four EF-hand domains. Most of the BcaCPK genes contain 5 to 8 introns, and these 123 BcaCPK genes are unevenly distributed across 17 chromosomes. Among these BcaCPK genes, 120 replicated gene pairs were found, whereas only 8 genes were tandem duplication, suggesting that dispersed duplication mainly drove the family amplification. The results of the Ka/Ks analysis indicated that the CPK gene family of B. carinata was primarily underwent purification selection in evolutionary selection. The promoter region of most BcaCPK genes contained various stress-related cis-acting elements. qRT-PCR analysis of 12 selected CPK genes conducted under cadmium and salt stress at various points revealed distinct expression patterns among different family members in response to different stresses. Specifically, the expression levels of BcaCPK2.B01a, BcaCPK16.B02b, and BcaCPK26.B02 were down-regulated under both stresses, whereas the expression levels of other members were significantly up-regulated under at least one stress.

CONCLUSION

This study systematically identified the BcaCPK gene family in B. carinata, which contributes to a better understanding the CPK genes in this species. The findings also serve as a reference for analyzing stress responses, particularly in relation to cadmium and salt stress in B. carinata.

DNA-PK and ATM drive phosphorylation signatures that antagonistically regulate cytokine responses to herpesvirus infection or DNA damage

The DNA-dependent protein kinase, DNA-PK, is an essential regulator of DNA damage repair. DNA-PK-driven phosphorylation events and the activated DNA damage response (DDR) pathways are also components of antiviral intrinsic and innate immune responses. Yet, it is not clear whether and how the DNA-PK response differs between these two forms of nucleic acid stress-DNA damage and DNA virus infection. Here, we define DNA-PK substrates and the signature cellular phosphoproteome response to DNA damage or infection with the nuclear-replicating DNA herpesvirus, HSV-1. We establish that DNA-PK negatively regulates the ataxia-telangiectasia-mutated (ATM) DDR kinase during viral infection. In turn, ATM blocks the binding of DNA-PK and the nuclear DNA sensor IFI16 to viral DNA, thereby inhibiting cytokine responses. However, following DNA damage, DNA-PK enhances ATM activity, which is required for IFN-β expression. These findings demonstrate that the DDR autoregulates cytokine expression through the opposing modulation of DDR kinases.

Evolutionary Conservation in Protein-Protein Interactions and Structures of the Elongator Sub-Complex ELP456 from Arabidopsis and Yeast

The Elongator complex plays a pivotal role in the wobble uridine modification of the tRNA anticodon. Comprising two sets of six distinct subunits, namely, Elongator proteins (ELP1-ELP6) and associated proteins, the holo-Elongator complex demonstrates remarkable functional and structural conservation across eukaryotes. However, the precise details of the evolutionary conservation of the holo-Elongator complex and its individual sub-complexes (i.e., ELP123; ELP456) in plants remain limited. In this study, we conducted an in vivo analysis of protein-protein interactions among Arabidopsis ELP4, ELP5, and ELP6 proteins. Additionally, we predicted their structural configurations and performed a comparative analysis with the structure of the yeast Elp456 sub-complex. Protein-protein interaction analysis revealed that AtELP4 interacts with AtELP6 but not directly with AtELP5. Furthermore, we found that the Arabidopsis Elongator-associated protein, Deformed Roots and Leaves 1 (DRL1), did not directly bind to AtELP proteins. The structural comparison of the ELP456 sub-complex between Arabidopsis and yeast demonstrated high similarity, encompassing the RecA-ATPase fold and the positions of hydrogen bonds, despite their relatively low sequence homology. Our findings suggest that Arabidopsis ELP4, ELP5, and ELP6 proteins form a heterotrimer, with ELP6 serving as a bridge, indicating high structural conservation between the ELP456 sub-complexes from Arabidopsis and yeast.

Genome mining of Mycobacterium tuberculosis: targeting SufD as a novel drug candidate through in silico characterization and inhibitor screening

Tuberculosis (TB) stands as the second most fatal infectious disease globally, causing 1.3 million deaths in 2022. The resurgence of TB and the alarming rise of antibiotic resistance demand urgent call to develop novel antituberculosis drugs. Despite concerted efforts to control TB, the disease persists and spreads rapidly on a global scale. Targeting stress response pathways in Mycobacterium tuberculosis (Mtb) has become imperative to achieve complete eradication. This study employs subtractive genomics to identify and prioritize potential drug targets among the hypothetical proteins of Mtb, focusing on indispensable pathways. Amongst 177 essential hypothetical proteins, 152 were nonhomologous to human. These proteins participated in 34 pathways, and a 20-fold enrichment of SUF pathway genes led to its selection as a target pathway. Fe-S clusters are fundamental, widely distributed protein cofactors involved in vital cellular processes. The survival of Mtb in a hypoxic environment relies on the iron-sulfur (Fe-S) cluster biogenesis pathway for the repair of damaged Fe-S clusters. It also protects pathogen against drugs, ensuring controlled iron utilization and contributing to drug resistance. In Mtb, six proteins of Fe-S cluster assembly pathway are encoded by the suf operon. The present study was focused on SufD because of its role in iron acquisition and prevention of Fenton reaction. The research further delves into the in silico characterization of SufD, utilizing bioinformatics tools for sequence and structure based analysis. The protein's structural features, including the identification of conserved regions, motifs, and 3D structure prediction enhanced functional annotation. Target based virtual screening of compounds from the ChEMBL database resulted in 12 inhibitors with best binding affinities. Drug likeness and ADMET profiling of potential inhibitors identified promising compounds with favorable drug-like properties. The study also involved cloning in SUMO-pRSF-Duet1 expression vector, overexpression, and purification of recombinant SufD from E. coli BL21 (DE3) cells. Optimization of expression conditions resulted in soluble production, and subsequent purification highlighting the efficacy of the SUMO fusion system for challenging Mtb proteins in E. coli. These findings provide valuable insights into pharmacological targets for future experimental studies, holding promise for the development of targeted therapy against Mtb.

C-terminally encoded peptide-like genes are associated with the development of primary root at qRL16.1 in soybean

Root architecture traits are belowground traits that harness moisture and nutrients from the soil and are equally important to above-ground traits in crop improvement. In soybean, the root length locus qRL16.1 was previously mapped on chromosome 16. The qRL16.1 has been characterized by transcriptome analysis of roots in near-isogenic lines (NILs), gene expression analysis in a pair of lines contrasting with alleles of qRL16.1, and differential gene expression analysis in germplasm accessions contrasting with root length. Two candidate genes, Glyma.16g108500 and Glyma.16g108700, have shown relatively higher expression in longer root accessions than in shorter rooting accessions. The C-terminal domain of Glyma.16g108500 and Glyma.16g108700 is similar to the conserved domain of C-terminally encoded peptides (CEPs) that regulate root length and nutrient response in Arabidopsis. Two polymorphisms upstream of Glyma.16g108500 showed a significant association with primary root length and total root length traits in a germplasm set. Synthetic peptide assay with predicted CEP variants of Glyma.16g108500 and Glyma.16g108700 demonstrated their positive effect on primary root length. The two genes are root-specific in the early stage of soybean growth and showed differential expression only in the primary root. These genes will be useful for improving soybean to develop a deep and robust root system to withstand low moisture and nutrient regimes.

RNA-Binding Protein Motifs Predict microRNA Secretion and Cellular Retention in Hypothalamic and Other Cell Types

Cellular microRNAs (miRNAs) can be selectively secreted or retained, adding another layer to their critical role in regulating human health and disease. To date, select RNA-binding proteins (RBPs) have been proposed to be a mechanism underlying miRNA localization, but the overall relevance of RBPs in systematic miRNA sorting remains unclear. This study profiles intracellular and small extracellular vesicles' (sEVs) miRNAs in NPY-expressing hypothalamic neurons. These findings were corroborated by the publicly available sEV and intracellular miRNA profiles of white and brown adipocytes, endothelium, liver, and muscle from various databases. Using experimentally determined binding motifs of 93 RBPs, our enrichment analysis revealed that sEV-originating miRNAs contained significantly different RBP motifs than those of intracellularly retained miRNAs. Multiple RBP motifs were shared across cell types; for instance, RBM4 and SAMD4 are significantly enriched in neurons, hepatocytes, skeletal muscle, and endothelial cells. Homologs of both proteins physically interact with Argonaute1/2 proteins, suggesting that they play a role in miRNA sorting. Machine learning modelling also demonstrates that significantly enriched RBP motifs could predict cell-specific preferential miRNA sorting. Non-optimized machine learning modeling of the motifs using Random Forest and Naive Bayes in all cell types except WAT achieved an area under the receiver operating characteristic (ROC) curve of 0.77-0.84, indicating a high predictive accuracy. Given that the RBP motifs have a significant predictive power, these results underscore the critical role that RBPs play in miRNA sorting within mammalian cells and reinforce the importance of miRNA sequencing in preferential localization. For the future development of small RNA therapeutics, considering these RBP-RNA interactions could be crucial to maximize delivery effectiveness and minimize off-target effects.

A bird's-eye view: exploration of the flavin-containing monooxygenase superfamily in common wheat

The Flavin Monooxygenase (FMO) gene superfamily in plants is involved in various processes most widely documented for its involvement in auxin biosynthesis, specialized metabolite biosynthesis, and plant microbial defense signaling. The roles of FMOs in defense signaling and disease resistance have recently come into focus as they may present opportunities to increase immune responses in plants including leading to systemic acquired resistance, but are not well characterized. We present a comprehensive catalogue of FMOs found in genomes across vascular plants and explore, in depth, 170 wheat TaFMO genes for sequence architecture, cis-acting regulatory elements, and changes due to Transposable Element insertions. A molecular phylogeny separates TaFMOs into three clades (A, B, and C) for which we further report gene duplication patterns, and differential rates of homoeologue expansion and retention among TaFMO subclades. We discuss Clade B TaFMOs where gene expansion is similarly seen in other cereal genomes. Transcriptome data from various studies point towards involvement of subclade B2 TaFMOs in disease responses against both biotrophic and necrotrophic pathogens, substantiated by promoter element analysis. We hypothesize that certain TaFMOs are responsive to both abiotic and biotic stresses, providing potential targets for enhancing disease resistance, plant yield and other important agronomic traits. Altogether, FMOs in wheat and other crop plants present an untapped resource to be exploited for improving the quality of crops.

Genome-Wide Identification and Analysis of the EIN3/EIL Transcription Factor Gene Family in Doubled Haploid (DH) Poplar

Ethylene (ET) is an important phytohormone that regulates plant growth, development and stress responses. The ethylene-insensitive3/ethylene-insensitive3-like (EIN3/EIL) transcription factor family, as a key regulator of the ET signal transduction pathway, plays an important role in regulating the expression of ET-responsive genes. Although studies of EIN3/EIL family members have been completed in many species, their role in doubled haploid (DH) poplar derived from another culture of diploid Populus simonii × P. nigra (donor tree, DT) remains ambiguous. In this study, a total of seven EIN3/EIL gene family members in the DH poplar genome were identified. Basic physical and chemical property analyses of these genes were performed, and these proteins were predicted to be localized to the nucleus. According to the phylogenetic relationship, EIN3/EIL genes were divided into two groups, and the genes in the same group had a similar gene structure and conserved motifs. The expression patterns of EIN3/EIL genes in the apical buds of different DH poplar plants were analyzed based on transcriptome data. At the same time, the expression patterns of PsnEIL1, PsnEIN3, PsnEIL4 and PsnEIL5 genes in different tissues of different DH plants were detected via RT-qPCR, including the apical buds, young leaves, functional leaves, xylem, cambium and roots. The findings presented above indicate notable variations in the expression levels of PsnEIL genes across various tissues of distinct DH plants. Finally, the PsnEIL1 gene was overexpressed in DT, and the transgenic plants showed a dwarf phenotype, indicating that the PsnEIL1 gene was involved in regulating the growth and development of poplar. In this study, the EIN3/EIL gene family of DH poplar was analyzed and functionally characterized, which provides a theoretical basis for the future exploration of the EIN3/EIL gene function.

The BAS chromatin remodeler determines brassinosteroid-induced transcriptional activation and plant growth in Arabidopsis

Brassinosteroid (BR) signaling leads to the nuclear accumulation of the BRASSINAZOLE-RESISTANT 1 (BZR1) transcription factor, which plays dual roles in activating or repressing the expression of thousands of genes. BZR1 represses gene expression by recruiting histone deacetylases, but how it activates transcription of BR-induced genes remains unclear. Here, we show that BR reshapes the genome-wide chromatin accessibility landscape, increasing the accessibility of BR-induced genes and reducing the accessibility of BR-repressed genes in Arabidopsis. BZR1 physically interacts with the BRAHMA-associated SWI/SNF (BAS)-chromatin-remodeling complex on the genome and selectively recruits the BAS complex to BR-activated genes. Depletion of BAS abrogates the capacities of BZR1 to increase chromatin accessibility, activate gene expression, and promote cell elongation without affecting BZR1's ability to reduce chromatin accessibility and expression of BR-repressed genes. Together, these data identify that BZR1 recruits the BAS complex to open chromatin and to mediate BR-induced transcriptional activation of growth-promoting genes.

Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

AMPpred-MFA: An Interpretable Antimicrobial Peptide Predictor with a Stacking Architecture, Multiple Features, and Multihead Attention

Antimicrobial peptides (AMPs) are small molecular polypeptides that can be widely used in the prevention and treatment of microbial infections. Although many computational models have been proposed to help identify AMPs, a high-performance and interpretable model is still lacking. In this study, new benchmark data sets are collected and processed, and a stacking deep architecture named AMPpred-MFA is carefully designed to discover and identify AMPs. Multiple features and a multihead attention mechanism are utilized on the basis of a bidirectional long short-term memory (LSTM) network and a convolutional neural network (CNN). The effectiveness of AMPpred-MFA is verified through five independent tests conducted in batches. Experimental results show that AMPpred-MFA achieves a state-of-the-art performance. The visualization interpretability analyses and ablation experiments offer a further understanding of the model behavior and performance, validating the importance of our feature representation and stacking architecture, especially the multihead attention mechanism. Therefore, AMPpred-MFA can be considered a reliable and efficient approach to understanding and predicting AMPs.

Rapid profiling of transcription factor-cofactor interaction networks reveals principles of epigenetic regulation

Transcription factor (TF)-cofactor (COF) interactions define dynamic, cell-specific networks that govern gene expression; however, these networks are understudied due to a lack of methods for high-throughput profiling of DNA-bound TF-COF complexes. Here we describe the Cofactor Recruitment (CoRec) method for rapid profiling of cell-specific TF-COF complexes. We define a lysine acetyltransferase (KAT)-TF network in resting and stimulated T cells. We find promiscuous recruitment of KATs for many TFs and that 35% of KAT-TF interactions are condition specific. KAT-TF interactions identify NF-κB as a primary regulator of acutely induced H3K27ac. Finally, we find that heterotypic clustering of CBP/P300-recruiting TFs is a strong predictor of total promoter H3K27ac. Our data supports clustering of TF sites that broadly recruit KATs as a mechanism for widespread co-occurring histone acetylation marks. CoRec can be readily applied to different cell systems and provides a powerful approach to define TF-COF networks impacting chromatin state and gene regulation.

LogoMotif: A Comprehensive Database of Transcription Factor Binding Site Profiles in Actinobacteria

Actinobacteria undergo a complex multicellular life cycle and produce a wide range of specialized metabolites, including the majority of the antibiotics. These biological processes are controlled by intricate regulatory pathways, and to better understand how they are controlled we need to augment our insights into the transcription factor binding sites. Here, we present LogoMotif (https://logomotif.bioinformatics.nl), an open-source database for characterized and predicted transcription factor binding sites in Actinobacteria, along with their cognate position weight matrices and hidden Markov models. Genome-wide predictions of binding site locations in Streptomyces model organisms are supplied and visualized in interactive regulatory networks. In the web interface, users can freely access, download and investigate the underlying data. With this curated collection of actinobacterial regulatory interactions, LogoMotif serves as a basis for binding site predictions, thus providing users with clues on how to elicit the expression of genes of interest and guide genome mining efforts.

Genome-wide identification of HMT gene family explores BpHMT2 enhancing selenium accumulation and tolerance in Broussonetia papyrifera

Broussonetia papyrifera, a valuable feed resource, is known for its fast growth, wide adaptability, high protein content and strong selenium enrichment capacity. Selenomethionine (SeMet), the main selenium form in selenium fortification B. papyrifera, is safe for animals and this enhances its nutritional value as a feed resource. However, the molecular mechanisms underlying SeMet synthesis remain unclear. This study identified three homocysteine S-methyltransferase genes from the B. papyrifera genome. The phylogenetic tree demonstrated that BpHMTs were divided into two classes, and BpHMT2 in the Class 2-D subfamily evolved earlier and possesses more fundamental functions. On the basis of the correlation between gene expression levels and selenium content, BpHMT2 was identified as a key candidate gene associated with selenium tolerance. Subcellular localization experiments confirmed the targeting of BpHMT2 in nucleus, cell membrane and chloroplasts. Moreover, three BpHMT2 overexpression Arabidopsis thaliana lines were confirmed to enhance plant selenium tolerance and SeMet accumulation. Overall, our finding provides insights into the molecular mechanisms of selenium metabolism in B. papyrifera, highlighting the potential role of BpHMT2 in SeMet synthesis. This research contributes to our understanding of selenium-enriched feed resources, with increased SeMet content contributing to the improved nutritional value of B. papyrifera as a feed resource.

Genome-wide identification and in-silico expression analysis of CCO gene family in sunflower (Helianthus annnus) against abiotic stress

Carotenoid cleavage oxygenases (CCOs) enzymes play an important role in plant growth and development by producing a wide array of apocarotenoids and their derivatives. These compounds are vital for colouring flowers and fruits and synthesizing plant hormones such as abscisic acid and strigolactones. Despite their importance, the gene family responsible for CCO enzymes in sunflowers has not been identified. In this study, we identify the CCO genes of the sunflower plant to fill this knowledge gap. Phylogenetic and synteny analysis indicated that the Helianthus annnus CCO (HaCCO) genes were conserved in different plant species and they could be divided into three subgroups based on their conserved domains. Analysis using MEME tool and multiple sequence alignment identified conserved motifs in the HaCCO gene sequence. Cis-regulatory elements (CREs) analysis of the HaCCO genes indicated the presence of various responsive elements related to plant hormones, development, and responses to both biotic and abiotic stresses. This implies that these genes may respond to plant hormones, developmental cues, and drought stress, offering potential applications in the development of more resistant crops. Genes belonging to the 9-cis-epoxy carotenoid dioxygenases (NCED) subgroups predominantly exhibited chloroplast localization, whereas the genes found in other groups are primarily localized in the cytoplasm. These 21 identified HaCCOs were regulated by 60 miRNAs, indicating the crucial role of microRNAs in gene regulation in sunflowers. Gene expression analysis under drought stress revealed significant up-regulation of HaNCED16 and HaNCED19, genes that are pivotal in ABA hormone biosynthesis. During organ-specific gene expression analysis, HaCCD12 and HaCCD20 genes exhibit higher activity in leaves, indicating a potential role in leaf pigmentation. This study provides a foundation for future research on the regulation and functions of the CCO gene family in sunflower and beyond. There is potential for developing molecular markers that could be employed in breeding programs to create new sunflower lines resistant to biotic and abiotic stresses.

Flowering genes identification, network analysis, and database construction for 837 plants

Flowering is one of the most important biological phenomena in the plant kingdom, which not only has important ecological significance, but also has substantial horticultural ornamental value. In this study, we undertook an exhaustive review of the advancements in our understanding of plant flowering genes. We delved into the identification and conducted comparative analyses of flowering genes across virtually all sequenced angiosperm plant genomes. Furthermore, we established an extensive angiosperm flowering atlas, encompassing a staggering 183 720 genes across eight pathways, along with 10 155 ABCDE mode genes, which play a pivotal role in plant flowering regulation. Through the examination of expression patterns, we unveiled the specificities of these flowering genes. An interaction network between flowering genes of the ABCDE model and their corresponding upstream genes offered a blueprint for comprehending their regulatory mechanisms. Moreover, we predicted the miRNA and target genes linked to the flowering processes of each species. To culminate our efforts, we have built a user-friendly web interface, named the Plant Flowering-time Gene Database (PFGD), accessible at http://pfgd.bio2db.com/. We firmly believe that this database will serve as a cornerstone in the global research community, facilitating the in-depth exploration of flowering genes in the plant kingdom. In summation, this pioneering endeavor represents the first comprehensive collection and comparative analysis of flowering genes in plants, offering valuable resources for the study of plant flowering genetics.

Characterization of the Fatty Acyl-CoA Reductase (FAR) Gene Family and Its Response to Abiotic Stress in Rice (Oryza sativa L.)

Fatty acyl-CoA reductase (FAR) is an important NADPH-dependent enzyme that can produce primary alcohol from fatty acyl-CoA or fatty acyl-carrier proteins as substrates. It plays a pivotal role in plant growth, development, and stress resistance. Herein, we performed genome-wide identification and expression analysis of FAR members in rice using bioinformatics methods. A total of eight OsFAR genes were identified, and the OsFARs were comprehensively analyzed in terms of phylogenetic relationships, duplication events, protein motifs, etc. The cis-elements of the OsFARs were predicted to respond to growth and development, light, hormones, and abiotic stresses. Gene ontology annotation analysis revealed that OsFAR proteins participate in biological processes as fatty acyl-CoA reductase during lipid metabolism. Numerous microRNA target sites were present in OsFARs mRNAs. The expression analysis showed that OsFARs were expressed at different levels during different developmental periods and in various tissues. Furthermore, the expression levels of OsFARs were altered under abiotic stresses, suggesting that FARs may be involved in abiotic stress tolerance in rice. The findings presented here serve as a solid basis for further exploring the functions of OsFARs.

TbsP and TrmB jointly regulate gapII to influence cell development phenotypes in the archaeon Haloferax volcanii

Microbial cells must continually adapt their physiology in the face of changing environmental conditions. Archaea living in extreme conditions, such as saturated salinity, represent important examples of such resilience. The model salt-loving organism Haloferax volcanii exhibits remarkable plasticity in its morphology, biofilm formation, and motility in response to variations in nutrients and cell density. However, the mechanisms regulating these lifestyle transitions remain unclear. In prior research, we showed that the transcriptional regulator, TrmB, maintains the rod shape in the related species Halobacterium salinarum by activating the expression of enzyme-coding genes in the gluconeogenesis metabolic pathway. In Hbt. salinarum, TrmB-dependent production of glucose moieties is required for cell surface glycoprotein biogenesis. Here, we use a combination of genetics and quantitative phenotyping assays to demonstrate that TrmB is essential for growth under gluconeogenic conditions in Hfx. volcanii. The ∆trmB strain rapidly accumulated suppressor mutations in a gene encoding a novel transcriptional regulator, which we name trmB suppressor, or TbsP (a.k.a. "tablespoon"). TbsP is required for adhesion to abiotic surfaces (i.e., biofilm formation) and maintains wild-type cell morphology and motility. We use functional genomics and promoter fusion assays to characterize the regulons controlled by each of TrmB and TbsP, including joint regulation of the glucose-dependent transcription of gapII, which encodes an important gluconeogenic enzyme. We conclude that TrmB and TbsP coregulate gluconeogenesis, with downstream impacts on lifestyle transitions in response to nutrients in Hfx. volcanii.

Changes in chromatin accessibility and transcriptional landscape induced by HDAC inhibitors in TP53 mutated patient-derived colon cancer organoids

Here we present the generation and characterization of patient-derived organoids (PDOs) from colorectal cancer patients. PDOs derived from two patients with TP53 mutations were tested with two different HDAC inhibitors (SAHA and NKL54). Cell death induction, transcriptome, and chromatin accessibility changes were analyzed. HDACIs promote the upregulation of low expressed genes and the downregulation of highly expressed genes. A similar differential effect is observed at the level of chromatin accessibility. Only SAHA is a potent inducer of cell death, which is characterized by the upregulation of BH3-only genes BIK and BMF. Up-regulation of BIK is associated with increased accessibility in an intronic region that has enhancer properties. SAHA, but not NKL54, also causes downregulation of BCL2L1 and decreases chromatin accessibility in three distinct regions of the BCL2L1 locus. Both inhibitors upregulate the expression of innate immunity genes and members of the MHC family. In summary, our exploratory study indicates a mechanism of action for SAHA and demonstrate the low efficacy of NKL54 as a single agent for apoptosis induction, using two PDOs. These observations need to be validated in a larger cohort of PDOs.

Genomic characterization and expression profiling of the lytic polysaccharide monooxygenases AA9 family in thermophilic fungi Thermothelomyces fergusii in response to carbon source media

Enhancing enzyme activity and stability in biomass degradation can improve substrate saccharification and, increases biorefinery efficiency. For the first time, we identified 20 lytic polysaccharide monooxygenases (LPMOs) AA9 genes in the genome of Thermothelomyces fergusii. Our results showed that TfAA9 was categorized into LPMOs1, LPMOs2, and LPMOs3 subgroups based on protein diversity. Protein- 3D structure analysis showed strong interactions between Myceliophthora thermophila AA9 proteins and 17 TfAA9 proteins. Gene ontology analysis indicated a high enrichment of cellulase activity in TfAA9 genes. KEGG pathways analysis revealed the role of TfAA9 proteins in the endohydrolysis of 1,4-beta-D-glucosidic linkages in cellulose. Numerous TfAA9s gene transcripts were up-regulated on avicel, cellobiose, and glucose, with a higher proportion on avicel. Protein concentration, endoglucanase, and cellulase activity were also boosted on avicel. However, limited fungal biomass was observed on avicel, despite the abundance of AA9 LPMOs in the T. fergusii genome. These findings expand our understanding of fungal AA9 genes and their role in lignocellulolytic degradation. The disparity between biomass and enzymatic activity suggests screening TfAA9 genes for highly active enzymes and redundant genes via heterologous expression. In short, functional characterization of these genes could contribute to improving the saccharification process of industrial raw materials.

Genome-Wide Identification of the SPP/SPPL Gene Family and BnaSPPL4 Regulating Male Fertility in Rapeseed (Brassica napus L.)

Signal peptide peptidase (SPP) and its homologs, signal peptide peptidase-like (SPPL) proteases, are members of the GxGD-type aspartyl protease family, which is widespread in plants and animals and is a class of transmembrane proteins with significant biological functions. SPP/SPPLs have been identified; however, the functions of SPP/SPPL in rapeseed (Brassica napus L.) have not been reported. In this study, 26 SPP/SPPLs were identified in rapeseed and categorized into three groups: SPP, SPPL2, and SPPL3. These members mainly contained the Peptidase_A22 and PA domains, which were distributed on 17 out of 19 chromosomes. Evolutionary analyses indicated that BnaSPP/SPPLs evolved with a large number of whole-genome duplication (WGD) events and strong purifying selection. Members are widely expressed and play a key role in the growth and development of rapeseed. The regulation of rapeseed pollen fertility by the BnaSPPL4 gene was further validated through experiments based on bioinformatics analysis, concluding that BnaSPPL4 silencing causes male sterility. Cytological observation showed that male infertility caused by loss of BnaSPPL4 gene function occurs late in the mononucleate stage due to microspore dysplasia.

Premeiotic 24-nt phasiRNAs are present in the Zea genus and unique in biogenesis mechanism and molecular function

Reproductive phasiRNAs are broadly present in angiosperms and play crucial roles in sustaining male fertility. While the premeiotic 21-nt phasiRNAs and meiotic 24-nt phasiRNA pathways have been extensively studied in maize (Zea mays) and rice (Oryza sativa), a third putative category of reproductive phasiRNAs-named premeiotic 24-nt phasiRNAs-have recently been reported in barley (Hordeum vulgare) and wheat (Triticum aestivum). To determine whether premeiotic 24-nt phasiRNAs are also present in maize and related species and begin to characterize their biogenesis and function, we performed a comparative transcriptome and degradome analysis of premeiotic and meiotic anthers from five maize inbred lines and three teosinte species/subspecies. Our data indicate that a substantial subset of the 24-nt phasiRNA loci in maize and teosinte are already highly expressed at premeiotic phase. The premeiotic 24-nt phasiRNAs are similar to meiotic 24-nt phasiRNAs in genomic origin and dependence on DCL5 for biogenesis, however, premeiotic 24-nt phasiRNAs are unique in that they are likely (i) not triggered by microRNAs, (ii) not loaded by AGO18 proteins, and (iii) not capable of mediating cis-cleavage. In addition, we also observed a group of premeiotic 24-nt phasiRNAs in rice using previously published data. Together, our results indicate that the premeiotic 24-nt phasiRNAs constitute a unique class of reproductive phasiRNAs and are present more broadly in the grass family (Poaceae) than previously known.

Genome-wide identification and expression analysis of TIFY genes under MeJA, cold and PEG-induced drought stress treatment in Dendrobium huoshanense

The TIFY family consists of plant-specific genes that regulates multiple plant functions, including developmental and defense responses. Here, we performed a comprehensive genomic analysis of TIFY genes in Dendrobium huoshanense. Our analysis encompassed their phylogenetic relationships, gene structures, chromosomal distributions, promoter regions, and patterns of collinearity. A total of 16 DhTIFY genes were identified, and classified into distinct clusters named JAZ, PPD, ZIM, and TIFY based on their phylogenetic relationship. These DhTIFYs exhibited an uneven distribution across 7 chromosomes. The expansion of the DhTIFY gene family appears to have been significantly influenced by whole-genome and segmental duplication events. The ratio of non-synonymous to synonymous substitutions (Ka/Ks) implies that the purifying selection has been predominant, maintaining a constrained functional diversification after duplication events. Gene structure analysis indicated that DhTIFYs exhibited significant structural variation, particularly in terms of gene organization and intron numbers. Moreover, numerous cis-acting elements related to hormone signaling, developmental processes, and stress responses were identified within the promoter regions. Subsequently, qRT-PCR experiments demonstrated that the expression of DhTIFYs is modulated in response to MeJA (Methyl jasmonate), cold, and drought treatment. Collectively, these results enhance our understanding of the functional dynamics of TIFY genes in D. huoshanense and may pinpoint potential candidates for detailed examination of the biological roles of TIFY genes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01442-9.

AIRE relies on Z-DNA to flag gene targets for thymic T cell tolerization

AIRE is an unconventional transcription factor that enhances the expression of thousands of genes in medullary thymic epithelial cells and promotes clonal deletion or phenotypic diversion of self-reactive T cells1-4. The biological logic of AIRE's target specificity remains largely unclear as, in contrast to many transcription factors, it does not bind to a particular DNA sequence motif. Here we implemented two orthogonal approaches to investigate AIRE's cis-regulatory mechanisms: construction of a convolutional neural network and leveraging natural genetic variation through analysis of F1 hybrid mice5. Both approaches nominated Z-DNA and NFE2-MAF as putative positive influences on AIRE's target choices. Genome-wide mapping studies revealed that Z-DNA-forming and NFE2L2-binding motifs were positively associated with the inherent ability of a gene's promoter to generate DNA double-stranded breaks, and promoters showing strong double-stranded break generation were more likely to enter a poised state with accessible chromatin and already-assembled transcriptional machinery. Consequently, AIRE preferentially targets genes with poised promoters. We propose a model in which Z-DNA anchors the AIRE-mediated transcriptional program by enhancing double-stranded break generation and promoter poising. Beyond resolving a long-standing mechanistic conundrum, these findings suggest routes for manipulating T cell tolerance.

An alfalfa MYB-like transcriptional factor MsMYBH positively regulates alfalfa seedling drought resistance and undergoes MsWAV3-mediated degradation

Drought is a major threat to alfalfa (Medicago sativa L.) production. The discovery of important alfalfa genes regulating drought response will facilitate breeding for drought-resistant alfalfa cultivars. Here, we report a genome-wide association study of drought resistance in alfalfa. We identified and functionally characterized an MYB-like transcription factor gene (MsMYBH), which increases the drought resistance in alfalfa. Compared with the wild-types, the biomass and forage quality were enhanced in MsMYBH overexpressed plants. Combined RNA-seq, proteomics and chromatin immunoprecipitation analysis showed that MsMYBH can directly bind to the promoters of MsMCP1, MsMCP2, MsPRX1A and MsCARCAB to improve their expression. The outcomes of such interactions include better water balance, high photosynthetic efficiency and scavenge excess H2O2 in response to drought. Furthermore, an E3 ubiquitin ligase (MsWAV3) was found to induce MsMYBH degradation under long-term drought, via the 26S proteasome pathway. Furthermore, variable-number tandem repeats in MsMYBH promoter were characterized among a collection of germplasms, and the variation is associated with promoter activity. Collectively, our findings shed light on the functions of MsMYBH and provide a pivotal gene that could be leveraged for breeding drought-resistant alfalfa. This discovery also offers new insights into the mechanisms of drought resistance in alfalfa.

Arabidopsis BBX14 is involved in high light acclimation and seedling development

The development of photosynthetically competent seedlings requires both light and retrograde biogenic signaling pathways. The transcription factor GLK1 functions at the interface between these pathways and receives input from the biogenic signal integrator GUN1. BBX14 was previously identified, together with GLK1, in a core module that mediates the response to high light (HL) levels and biogenic signals, which was studied by using inhibitors of chloroplast development. Our chromatin immunoprecipitation-Seq experiments revealed that BBX14 is a direct target of GLK1, and RNA-Seq analysis suggests that BBX14 may function as a regulator of the circadian clock. In addition, BBX14 plays a role in chlorophyll biosynthesis during early onset of light. Knockout of BBX14 results in a long hypocotyl phenotype dependent on a retrograde signal. Furthermore, the expression of BBX14 and BBX15 during biogenic signaling requires GUN1. Investigation of the role of BBX14 and BBX15 in GUN-type biogenic (gun) signaling showed that the overexpression of BBX14 or BBX15 caused de-repression of CA1 mRNA levels, when seedlings were grown on norflurazon. Notably, transcripts of the LHCB1.2 marker are not de-repressed. Furthermore, BBX14 is required to acclimate plants to HL stress. We propose that BBX14 is an integrator of biogenic signals and that BBX14 is a nuclear target of retrograde signals downstream of the GUN1/GLK1 module. However, we do not classify BBX14 or BBX15 overexpressors as gun mutants based on a critical evaluation of our results and those reported in the literature. Finally, we discuss a classification system necessary for the declaration of new gun mutants.

Genomic data and ecological niche modeling reveal an unusually slow rate of molecular evolution in the Cretaceous Eupteleaceae

Living fossils are evidence of long-term sustained ecological success. However, whether living fossils have little molecular changes remains poorly known, particularly in plants. Here, we have introduced a novel method that integrates phylogenomic, comparative genomic, and ecological niche modeling analyses to investigate the rate of molecular evolution of Eupteleaceae, a Cretaceous relict angiosperm family endemic to East Asia. We assembled a high-quality chromosome-level nuclear genome, and the chloroplast and mitochondrial genomes of a member of Eupteleaceae (Euptelea pleiosperma). Our results show that Eupteleaceae is most basal in Ranunculales, the earliest-diverging order in eudicots, and shares an ancient whole-genome duplication event with the other Ranunculales. We document that Eupteleaceae has the slowest rate of molecular changes in the observed angiosperms. The unusually low rate of molecular evolution of Eupteleaceae across all three independent inherited genomes and genes within each of the three genomes is in association with its conserved genome architecture, ancestral woody habit, and conserved niche requirements. Our findings reveal the evolution and adaptation of living fossil plants through large-scale environmental change and also provide new insights into early eudicot diversification.

Sorghum SbGhd7 is a major regulator of floral transition and directly represses genes crucial for flowering activation

Molecular genetic understanding of flowering time regulation is crucial for sorghum development. GRAIN NUMBER, PLANT HEIGHT AND HEADING DATE 7 (SbGhd7) is one of the six classical loci conferring photoperiod sensitivity of sorghum flowering. However, its functions remain poorly studied. The molecular functions of SbGhd7 were characterized. The gene regulatory network controlled by SbGhd7 was constructed and validated. The biological roles of SbGhd7 and its major targets were studied. SbGhd7 overexpression (OE) completely prevented sorghum flowering. Additionally, we show that SbGhd7 is a major negative regulator of flowering, binding to the promoter motif TGAATG(A/T)(A/T/C) and repressing transcription of the major florigen FLOWERING LOCUS T 10 (SbFT10) and floral activators EARLY HEADING DATE (SbEhd1), FLAVIN-BINDING, KELCH REPEAT, F-BOX1 (SbFKF1) and EARLY FLOWERING 3 (SbELF3). Reinforcing the direct effect of SbGhd7, SbEhd1 OE activated the promoters of three functional florigens (SbFT1, SbFT8 and SbFT10), dramatically accelerating flowering. Our studies demonstrate that SbGhd7 is a major repressor of sorghum flowering by directly and indirectly targeting genes for flowering activation. The mechanism appears ancient. Our study extends the current model of floral transition regulation in sorghum and provides a framework for a comprehensive understanding of sorghum photoperiod response.

Identification of a viral gene essential for the genome replication of a domesticated endogenous virus in ichneumonid parasitoid wasps

Thousands of endoparasitoid wasp species in the families Braconidae and Ichneumonidae harbor "domesticated endogenous viruses" (DEVs) in their genomes. This study focuses on ichneumonid DEVs, named ichnoviruses (IVs). Large quantities of DNA-containing IV virions are produced in ovary calyx cells during the pupal and adult stages of female wasps. Females parasitize host insects by injecting eggs and virions into the body cavity. After injection, virions rapidly infect host cells which is followed by expression of IV genes that promote the successful development of wasp offspring. IV genomes consist of two components: proviral segment loci that serve as templates for circular dsDNAs that are packaged into capsids, and genes from an ancestral virus that produce virions. In this study, we generated a chromosome-scale genome assembly for Hyposoter didymator that harbors H. didymator ichnovirus (HdIV). We identified a total of 67 HdIV loci that are amplified in calyx cells during the wasp pupal stage. We then focused on an HdIV gene, U16, which is transcribed in calyx cells during the initial stages of replication. Sequence analysis indicated that U16 contains a conserved domain in primases from select other viruses. Knockdown of U16 by RNA interference inhibited virion morphogenesis in calyx cells. Genome-wide analysis indicated U16 knockdown also inhibited amplification of HdIV loci in calyx cells. Altogether, our results identified several previously unknown HdIV loci, demonstrated that all HdIV loci are amplified in calyx cells during the pupal stage, and showed that U16 is required for amplification and virion morphogenesis.

Populus MYC2 orchestrates root transcriptional reprogramming of defence pathway to impair Laccaria bicolor ectomycorrhizal development

The jasmonic acid (JA) signalling pathway plays an important role in the establishment of the ectomycorrhizal symbiosis. The Laccaria bicolor effector MiSSP7 stabilizes JA corepressor JAZ6, thereby inhibiting the activity of Populus MYC2 transcription factors. Although the role of MYC2 in orchestrating plant defences against pathogens is well established, its exact contribution to ECM symbiosis remains unclear. This information is crucial for understanding the balance between plant immunity and symbiotic relationships. Transgenic poplars overexpressing or silencing for the two paralogues of MYC2 transcription factor (MYC2s) were produced, and their ability to establish ectomycorrhiza was assessed. Transcriptomics and DNA affinity purification sequencing were performed. MYC2s overexpression led to a decrease in fungal colonization, whereas its silencing increased it. The enrichment of terpene synthase genes in the MYC2-regulated gene set suggests a complex interplay between the host monoterpenes and fungal growth. Several root monoterpenes have been identified as inhibitors of fungal growth and ECM symbiosis. Our results highlight the significance of poplar MYC2s and terpenes in mutualistic symbiosis by controlling root fungal colonization. We identified poplar genes which direct or indirect control by MYC2 is required for ECM establishment. These findings deepen our understanding of the molecular mechanisms underlying ECM symbiosis.

Genome-wide identification and evolutionary analysis of the AP2/EREBP, COX and LTP genes in Zea mays L. under drought stress

AP2 (APETALA2)/EREBP (ethylene-responsive element-binding protein), cytochrome c oxidase (COX) and nonspecific lipid transfer proteins (LTP) play important roles in the response to drought stress. This is the first study to identify the COX gene in Zea mays L. via genome-wide analysis. The qRT‒PCR results indicated that AP2/EREBP, COX and LTP were downregulated, with fold changes of 0.84, 0.53 and 0.31, respectively, after 12 h of drought stress. Genome-wide analysis identified 78 AP2/EREBP, 6 COX and 10 LTP genes in Z. mays L. Domain analysis confirmed the presence of the AP2 domain, Cyt_c_Oxidase_Vb domain and nsLTP1 in the AP2/EREBP, COX and LTP proteins, respectively. The AP2/EREBP protein family (AP2) includes five different domain types: the AP2/ERF domain, the EREBP-like factor (EREBP), the ethylene responsive factor (ERF), the dehydration responsive element binding protein (DREB) and the SHN SHINE. Synteny analysis of the AP2/EREBP, COX and LTP genes revealed collinearity orthologous relationships in O. sativa, H. vulgare and A. thaliana. AP2/EREBP genes were found on the 10 chromosomes of Z. mays L. COX genes were found on chromosomes 1, 3, 4, 5, 7 and 8. LTP genes were found on chromosomes 1, 3, 6, 8, 9 and 10. In the present study, the Ka/Ks ratios of the AP2/EREBP paralogous pairs indicated that the AP2/EREBP genes were influenced primarily by purifying selection, which indicated that the AP2/EREBP genes received strong environmental pressure during evolution. The Ka/Ks ratios of the COX-3/COX-4 paralogous pairs indicate that the COX-3/COX-4 genes were influenced primarily by Darwinian selection (driving change). For the LTP genes, the Ka/Ks ratios of the LTP-1/LTP-10, LTP-5/LTP-3 and LTP-4/LTP-8 paralogous pairs indicate that these genes were influenced primarily by purifying selection, while the Ka/Ks ratios of the LTP-2/LTP-6 paralogous pairs indicate that these genes were influenced primarily by Darwinian selection. The duplication time of the AP2/EREBP paralogous gene pairs in Z. mays L. ranged from approximately 9.364 to 100.935 Mya. The duplication time of the COX-3/COX-4 paralogous gene pair was approximately 5.217 Mya. The duplication time of the LTP paralogous gene pairs ranged from approximately 19.064 to 96.477 Mya. The major focus of research is to identify the genes that are responsible for drought stress tolerance to improve maize for drought stress tolerance. The results of the present study will improve the understanding of the functions of the AP2/EREBP, COX and LTP genes in response to drought stress.

The Evolution, Expression Patterns, and Domestication Selection Analysis of the Annexin Gene Family in the Barley Pan-Genome

Plant annexins constitute a conserved protein family that plays crucial roles in regulating plant growth and development, as well as in responses to both biotic and abiotic stresses. In this study, a total of 144 annexin genes were identified in the barley pan-genome, comprising 12 reference genomes, including cultivated barley, landraces, and wild barley. Their chromosomal locations, physical-chemical characteristics, gene structures, conserved domains, and subcellular localizations were systematically analyzed to reveal the certain differences between wild and cultivated populations. Through a cis-acting element analysis, co-expression network, and large-scale transcriptome analysis, their involvement in growth, development, and responses to various stressors was highlighted. It is worth noting that HvMOREXann5 is only expressed in pistils and anthers, indicating its crucial role in reproductive development. Based on the resequencing data from 282 barley accessions worldwide, genetic variations in thefamily were investigated, and the results showed that 5 out of the 12 identified HvMOREXanns were affected by selection pressure. Genetic diversity and haplotype frequency showed notable reductions between wild and domesticated barley, suggesting that a genetic bottleneck occurred on the annexin family during the barley domestication process. Finally, qRT-PCR analysis confirmed the up-regulation of HvMOREXann7 under drought stress, along with significant differences between wild accessions and varieties. This study provides some insights into the genome organization and genetic characteristics of the annexin gene family in barley at the pan-genome level, which will contribute to better understanding its evolution and function in barley and other crops.

NSD2 is a requisite subunit of the AR/FOXA1 neo-enhanceosome in promoting prostate tumorigenesis

The androgen receptor (AR) is a ligand-responsive transcription factor that binds at enhancers to drive terminal differentiation of the prostatic luminal epithelia. By contrast, in tumors originating from these cells, AR chromatin occupancy is extensively reprogrammed to drive hyper-proliferative, metastatic, or therapy-resistant phenotypes, the molecular mechanisms of which remain poorly understood. Here, we show that the tumor-specific enhancer circuitry of AR is critically reliant on the activity of Nuclear Receptor Binding SET Domain Protein 2 (NSD2), a histone 3 lysine 36 di-methyltransferase. NSD2 expression is abnormally gained in prostate cancer cells and its functional inhibition impairs AR trans-activation potential through partial off-loading from over 40,000 genomic sites, which is greater than 65% of the AR tumor cistrome. The NSD2-dependent AR sites distinctly harbor a chimeric AR-half motif juxtaposed to a FOXA1 element. Similar chimeric motifs of AR are absent at the NSD2-independent AR enhancers and instead contain the canonical palindromic motifs. Meta-analyses of AR cistromes from patient tumors uncovered chimeric AR motifs to exclusively participate in tumor-specific enhancer circuitries, with a minimal role in the physiological activity of AR. Accordingly, NSD2 inactivation attenuated hallmark cancer phenotypes that were fully reinstated upon exogenous NSD2 re-expression. Inactivation of NSD2 also engendered increased dependency on its paralog NSD1, which independently maintained AR and MYC hyper-transcriptional programs in cancer cells. Concordantly, a dual NSD1/2 PROTAC degrader, called LLC0150, was preferentially cytotoxic in AR-dependent prostate cancer as well as NSD2-altered hematologic malignancies. Altogether, we identify NSD2 as a novel subunit of the AR neo-enhanceosome that wires prostate cancer gene expression programs, positioning NSD1/2 as viable paralog co-targets in advanced prostate cancer.

Omics-based construction of regulatory variants can be applied to help decipher pig liver-related traits

Genetic variants can influence complex traits by altering gene expression through changes to regulatory elements. However, the genetic variants that affect the activity of regulatory elements in pigs are largely unknown, and the extent to which these variants influence gene expression and contribute to the understanding of complex phenotypes remains unclear. Here, we annotate 90,991 high-quality regulatory elements using acetylation of histone H3 on lysine 27 (H3K27ac) ChIP-seq of 292 pig livers. Combined with genome resequencing and RNA-seq data, we identify 28,425 H3K27ac quantitative trait loci (acQTLs) and 12,250 expression quantitative trait loci (eQTLs). Through the allelic imbalance analysis, we validate two causative acQTL variants in independent datasets. We observe substantial sharing of genetic controls between gene expression and H3K27ac, particularly within promoters. We infer that 46% of H3K27ac exhibit a concomitant rather than causative relationship with gene expression. By integrating GWAS, eQTLs, acQTLs, and transcription factor binding prediction, we further demonstrate their application, through metabolites dulcitol, phosphatidylcholine (PC) (16:0/16:0) and published phenotypes, in identifying likely causal variants and genes, and discovering sub-threshold GWAS loci. We provide insight into the relationship between regulatory elements and gene expression, and the genetic foundation for dissecting the molecular mechanism of phenotypes.

Unraveling the role of Xist in X chromosome inactivation: insights from rabbit model and deletion analysis of exons and repeat A

X chromosome inactivation (XCI) is a process that equalizes the expression of X-linked genes between males and females. It relies on Xist, continuously expressed in somatic cells during XCI maintenance. However, how Xist impacts XCI maintenance and its functional motifs remain unclear. In this study, we conducted a comprehensive analysis of Xist, using rabbits as an ideal non-primate model. Homozygous knockout of exon 1, exon 6, and repeat A in female rabbits resulted in embryonic lethality. However, X∆ReAX females, with intact X chromosome expressing Xist, showed no abnormalities. Interestingly, there were no significant differences between females with homozygous knockout of exons 2-5 and wild-type rabbits, suggesting that exons 2, 3, 4, and 5 are less important for XCI. These findings provide evolutionary insights into Xist function.

The transcription of the main gene associated with Treacher-Collins syndrome (TCOF1) is regulated by G-quadruplexes and cellular nucleic acid binding protein (CNBP)

Treacle ribosome biogenesis factor 1 (TCOF1) is responsible for about 80% of mandibular dysostosis (MD) cases. We have formerly identified a correlation between TCOF1 and CNBP (CCHC-type zinc finger nucleic acid binding protein) expression in human mesenchymal cells. Given the established role of CNBP in gene regulation during rostral development, we explored the potential for CNBP to modulate TCOF1 transcription. Computational analysis for CNBP binding sites (CNBP-BSs) in the TCOF1 promoter revealed several putative binding sites, two of which (Hs791 and Hs2160) overlap with putative G-quadruplex (G4) sequences (PQSs). We validated the folding of these PQSs measuring circular dichroism and fluorescence of appropriate synthetic oligonucleotides. In vitro studies confirmed binding of purified CNBP to the target PQSs (both folded as G4 and unfolded) with Kd values in the nM range. ChIP assays conducted in HeLa cells chromatin detected the CNBP binding to TCOF1 promoter. Transient transfections of HEK293 cells revealed that Hs2160 cloned upstream SV40 promoter increased transcription of downstream firefly luciferase reporter gene. We also detected a CNBP-BS and PQS (Dr2393) in the zebrafish TCOF1 orthologue promoter (nolc1). Disrupting this G4 in zebrafish embryos by microinjecting DNA antisense oligonucleotides complementary to Dr2393 reduced the transcription of nolc1 and recapitulated the craniofacial anomalies characteristic of Treacher Collins Syndrome. Both cnbp overexpression and Morpholino-mediated knockdown in zebrafish induced nolc1 transcription. These results suggest that CNBP modulates the transcriptional expression of TCOF1 through a mechanism involving G-quadruplex folding/unfolding, and that this regulation is active in vertebrates as distantly related as bony fish and humans. These findings may have implications for understanding and treating MD.