The MEME Suite

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.

Genome-Wide Identification, Expression Analysis under Abiotic Stress and Co-Expression Analysis of MATE Gene Family in Torreya grandis

The multidrug and toxin efflux (MATE) family participates in numerous biological processes and plays important roles in abiotic stress responses. However, information about the MATE family genes in Torreya grandis remains unclear. In this study, our genome-wide investigation identified ninety MATE genes in Torreya grandis, which were divided into five evolutionary clades. TgMATE family members are located on eleven chromosomes, and a total of thirty TgMATEs exist in tandem duplication. The promoter analysis showed that most TgMATEs contain the cis-regulatory elements associated with stress and hormonal responses. In addition, we discovered that most TgMATE genes responded to abiotic stresses (aluminum, drought, high temperatures, and low temperatures). Weighted correlation network analysis showed that 147 candidate transcription factor genes regulated the expression of 14 TgMATE genes, and it was verified through a double-luciferase assay. Overall, our findings offer valuable information for the characterization of the TgMATE gene mechanism in responding to abiotic stress and exhibit promising prospects for the stress tolerance breeding of Torreya grandis.

Identification of WRKY gene family in Dioscorea opposita Thunb. reveals that DoWRKY71 enhanced the tolerance to cold and ABA stress

WRKY transcription factors constitute one of the largest plant-specific gene families, regulating various aspects of plant growth, development, physiological processes, and responses to abiotic stresses. This study aimed to comprehensively analyze the WRKY gene family of yam (Dioscorea opposita Thunb.), to understand their expression patterns during the growth and development process and their response to different treatments of yam and analyze the function of DoWRKY71 in detail. A total of 25 DoWRKY genes were identified from the transcriptome of yam, which were divided into six clades (I, IIa, IIc, IId, IIe, III) based on phylogenetic analysis. The analysis of conserved motifs revealed 10 motifs, varying in length from 16 to 50 amino acids. Based on real-time quantitative PCR (qRT-PCR) analysis, DoWRKY genes were expressed at different stages of growth and development and responded differentially to various abiotic stresses. The expression level of DoWRKY71 genes was up-regulated in the early stage and then down-regulated in tuber enlargement. This gene showed responsiveness to cold and abiotic stresses, such as abscisic acid (ABA) and methyl jasmonate (MeJA). Therefore, further study was conducted on this gene. Subcellular localization analysis revealed that the DoWRKY71 protein was localized in the nucleus. Moreover, the overexpression of DoWRKY71 enhanced the cold tolerance of transgenic tobacco and promoted ABA mediated stomatal closure. This study presents the first systematic analysis of the WRKY gene family in yam, offering new insights for studying WRKY transcription factors in yam. The functional study of DoWRKY71 lays theoretical foundation for further exploring the regulatory function of the DoWRKY71 gene in the growth and development related signaling pathway of yam.

High resolution maps of chromatin reorganization through mouse meiosis reveal novel features of the 3D meiotic structure

When germ cells transition from the mitotic cycle into meiotic prophase I (MPI), chromosomes condense into an array of chromatin loops that are required to promote homolog pairing and genetic recombination. To identify the changes in chromosomal conformation, we isolated nuclei on a trajectory from spermatogonia to the end of MPI. At each stage along this trajectory, we built genomic interaction maps with the highest temporal and spatial resolution to date. The changes in chromatin folding coincided with a concurrent decline in mitotic cohesion and a rise in meiotic cohesin complexes. We found that the stereotypical large-scale A and B compartmentalization was lost during meiotic prophase I alongside the loss of topological associating domains (TADs). Still, local subcompartments were detected and maintained throughout meiosis. The enhanced Micro-C resolution revealed that, despite the loss of TADs, higher frequency contact sites between two loci were detectable during meiotic prophase I coinciding with CTCF bound sites. The pattern of interactions around these CTCF sites with their neighboring loci showed that CTCF sites were often anchoring the meiotic loops. Additionally, the localization of CTCF to the meiotic axes indicated that these anchors were at the base of loops. Strikingly, even in the face of the dramatic reconfiguration of interphase chromatin into a condensed loop-array, the interactions between regulatory elements remained well preserved. This establishes a potential mechanism for how the meiotic chromatin maintains active transcription within a highly structured genome. In summary, the high temporal and spatial resolution of these data revealed previously unappreciated aspects of mammalian meiotic chromatin organization.

Obesogenic Diet in Mice Leads to Inflammation and Oxidative Stress in the Mother in Association with Sex-Specific Changes in Fetal Development, Inflammatory Markers and Placental Transcriptome

BACKGROUND

Obesity during pregnancy is related to adverse maternal and neonatal outcomes. Factors involved in these outcomes may include increased maternal insulin resistance, inflammation, oxidative stress, and nutrient mishandling. The placenta is the primary determinant of fetal outcomes, and its function can be impacted by maternal obesity. The aim of this study on mice was to determine the effect of obesity on maternal lipid handling, inflammatory and redox state, and placental oxidative stress, inflammatory signaling, and gene expression relative to female and male fetal growth.

METHODS

Female mice were fed control or obesogenic high-fat/high-sugar diet (HFHS) from 9 weeks prior to, and during, pregnancy. On day 18.5 of pregnancy, maternal plasma, and liver, placenta, and fetal serum were collected to examine the immune and redox states. The placental labyrinth zone (Lz) was dissected for RNA-sequencing analysis of gene expression changes.

RESULTS

the HFHS diet induced, in the dams, hepatic steatosis, oxidative stress (reduced catalase, elevated protein oxidation) and the activation of pro-inflammatory pathways (p38-MAPK), along with imbalanced circulating cytokine concentrations (increased IL-6 and decreased IL-5 and IL-17A). HFHS fetuses were asymmetrically growth-restricted, showing sex-specific changes in circulating cytokines (GM-CSF, TNF-α, IL-6 and IFN-γ). The morphology of the placenta Lz was modified by an HFHS diet, in association with sex-specific alterations in the expression of genes and proteins implicated in oxidative stress, inflammation, and stress signaling. Placental gene expression changes were comparable to that seen in models of intrauterine inflammation and were related to a transcriptional network involving transcription factors, LYL1 and PLAG1.

CONCLUSION

This study shows that fetal growth restriction with maternal obesity is related to elevated oxidative stress, inflammatory pathways, and sex-specific placental changes. Our data are important, given the marked consequences and the rising rates of obesity worldwide.

An miR164-resistant mutation in the transcription factor gene CpCUC2B enhances carpel arrest and ectopic boundary specification in Cucurbita pepo flower development

The sex determination process in cucurbits involves the control of stamen or carpel development during the specification of male or female flowers from a bisexual floral meristem, a function coordinated by ethylene. A gain-of-function mutation in the miR164-binding site of CpCUC2B, ortholog of the Arabidopsis transcription factor gene CUC2, not only produced ectopic floral meristems and organs, but also suppressed the development of carpels and promoted the development of stamens. The cuc2b mutation induced the transcription of CpCUC2B in the apical shoots of plants after female flowering but repressed other CUC genes regulated by miR164, suggesting a conserved functional redundancy of these genes in the development of squash flowers. The synergistic androecious phenotype of the double mutant between cuc2b and etr2b, an ethylene-insensitive mutation that enhances the production of male flowers, demonstrated that CpCUC2B arrests the development of carpels independently of ethylene and CpWIP1B. The transcriptional regulation of CpCUC1, CpCUC2, and ethylene genes in cuc2b and ethylene mutants also confirms this conclusion. However, the epistasis of cuc2b over aco1a, a mutation that suppresses stamen arrest in female flowers, and the down-regulation of CpACS27A in cuc2b female apical shoots, indicated that CpCUC2B promotes stamen development by suppressing the late ethylene production.

Genome-wide characterization of heavy metal-associated isoprenylated plant protein gene family from Citrus sinensis in response to huanglongbing

Introduction

Heavy metal-associated isoprenylated plant proteins (HIPPs) play vital roles in maintaining heavy metal balance and responding to both biotic and abiotic stresses in vascular plants. However, the role of HIPPs in the response to Huanglongbing (HLB), a harmful disease of citrus caused by the phloem-colonizing bacterium Candidatus Liberibacter asiaticus (CLas), has not been examined.

Methods and results

In this study, a total of 26 HIPP genes were identified in Citrus sinensis, and they were grouped into 5 clades. The CsHIPP genes are distributed on 8 chromosomes and exhibited considerable synteny with HIPPs found in Arabidopsis thaliana. Additionally, we analyzed the gene structure, conserved motifs and domains of the CsHIPPs. Various cis-acting elements related to plant hormones and stress responses were identified in the promoters of CsHIPPs. Public transcriptome data and RT-qPCR analysis showed that the expression level of CsHIPP03 was significantly reduced in samples infected by CLas and Xanthomonas citri ssp. citri (Xcc). Furthermore, silencing the homologous gene of CsHIPP03 in Nicotiana benthamiana increased the disease resistance of plants to bacteria.

Discussion

Our results provide a basis for functional studies of HIPP gene family in C. sinensis, highlighting their functions in bacterial resistance, and improve our understanding to the susceptibility mechanism of HLB.

The complete mitochondrial genome of Leucoptera coffeella (Lepidoptera: Lyonetiidae) and phylogenetic relationships within the Yponomeutoidea superfamily

The coffee leaf miner (Leucoptera coffeella) is one of the major pests of coffee crops in the neotropical regions, and causes major economic losses. Few molecular data are available to identify this pest and advances in the knowledge of the genome of L. coffeella will contribute to improving pest identification and also clarify taxonomy of this microlepidoptera. L. coffeella DNA was extracted and sequenced using PacBio HiFi technology. Here we report the complete L. coffeella circular mitochondrial genome (16,407 bp) assembled using Aladin software. We found a total of 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and an A + T rich-region and a D-loop. The L. coffeella mitochondrial gene organization is highly conserved with similarities to lepidopteran mitochondrial gene rearrangements (trnM-trnI-trnQ). We concatenated the 13 PCG to construct a phylogenetic tree and inferred the relationship between L. coffeella and other lepidopteran species. L. coffeella is found in the Lyonetiidae clade together with L. malifoliella and Lyonetia clerkella, both leaf miners. Interestingly, this clade is assigned in the Yponomeutoidea superfamily together with Gracillariidae, and both superfamilies displayed species with leaf-mining feeding habits.

Tracking and characterization of a novel conjugative transposon identified by shotgun transposon mutagenesis

The horizontal transfer of mobile genetic elements (MGEs) is an essential process determining the functional and genomic diversity of bacterial populations. MGEs facilitate the exchange of fitness determinant genes like antibiotic resistance and virulence factors. Various computational methods exist to identify potential MGEs, but confirming their ability to transfer requires additional experimental approaches. Here, we apply a transposon (Tn) mutagenesis technique for confirming mobilization without the need for targeted mutations. Using this method, we identified two MGEs, including a previously known conjugative transposon (CTn) called BoCTn found in Bacteroides ovatus and a novel CTn, PvCTn, identified in Phocaeicola vulgatus. In addition, Tn mutagenesis and subsequent genetic deletion enabled our characterization of a helix-turn-helix motif gene, BVU3433 which negatively regulates the conjugation efficiency of PvCTn in vitro. Furthermore, our transcriptomics data revealed that BVU3433 plays a crucial role in the repression of PvCTn genes, including genes involved in forming complete conjugation machinery [Type IV Secretion System (T4SS)]. Finally, analysis of individual strain genomes and community metagenomes identified the widespread prevalence of PvCTn-like elements with putative BVU3433 homologs among human gut-associated bacteria. In summary, this Tn mutagenesis mobilization method (TMMM) enables observation of transfer events in vitro and can ultimately be applied in vivo to identify a broader diversity of functional MGEs that may underly the transfer of important fitness determinants.

Genome-wide analysis of the KNOX gene family in Moso bamboo: insights into their role in promoting the rapid shoot growth

BACKGROUND

KNOTTED1-like homeobox (KNOX) genes, plant-specific homologous box transcription factors (TFs), play a central role in regulating plant growth, development, organ formation, and response to biotic and abiotic stresses. However, a comprehensive genome-wide identification of the KNOX genes in Moso bamboo (Phyllostachys edulis), the fastest growing plant, has not yet been conducted, and the specific biological functions of this family remain unknown.

RESULTS

The expression profiles of 24 KNOX genes, divided into two subfamilies, were determined by integrating Moso bamboo genome and its transcriptional data. The KNOX gene promoters were found to contain several light and stress-related cis-acting elements. Synteny analysis revealed stronger similarity with rice KNOX genes than with Arabidopsis KNOX genes. Additionally, several conserved structural domains and motifs were identified in the KNOX proteins. The expansion of the KNOX gene family was primarily regulated by tandem duplications. Furthermore, the KNOX genes were responsive to naphthaleneacetic acid (NAA) and gibberellin (GA) hormones, exhibiting distinct temporal expression patterns in four different organs of Moso bamboo. Short Time-series Expression Miner (STEM) analysis and quantitative real-time PCR (qRT-PCR) assays demonstrated that PeKNOX genes may play a role in promoting rapid shoot growth. Additionally, Gene Ontology (GO) and Protein-Protein Interaction (PPI) network enrichment analyses revealed several functional annotations for PeKNOXs. By regulating downstream target genes, PeKNOXs are involved in the synthesis of AUX /IAA, ultimately affecting cell division and elongation.

CONCLUSIONS

In the present study, we identified and characterized a total of 24 KNOX genes in Moso bamboo and investigated their physiological properties and conserved structural domains. To understand their functional roles, we conducted an analysis of gene expression profiles using STEM and RNA-seq data. This analysis successfully revealed regulatory networks of the KNOX genes, involving both upstream and downstream genes. Furthermore, the KNOX genes are involved in the AUX/IAA metabolic pathway, which accelerates shoot growth by influencing downstream target genes. These results provide a theoretical foundation for studying the molecular mechanisms underlying the rapid growth and establish the groundwork for future research into the functions and transcriptional regulatory networks of the KNOX gene family.

Identification, Evolutionary Dynamics, and Gene Expression Patterns of the ACP Gene Family in Responding to Salt Stress in Brassica Genus

Acyl carrier proteins (ACPs) have been reported to play a crucial role in responding to biotic and abiotic stresses, regulating growth and development. However, the biological function of the ACP gene family in the Brassica genus has been limited until now. In this study, we conducted a comprehensive analysis and identified a total of 120 ACP genes across six species in the Brassica genus. Among these, there were 27, 26, and 30 ACP genes in the allotetraploid B. napus, B. juncea, and B. carinata, respectively, and 14, 13, and 10 ACP genes in the diploid B. rapa, B. oleracea, and B. nigra, respectively. These ACP genes were further classified into six subclades, each containing conserved motifs and domains. Interestingly, the majority of ACP genes exhibited high conservation among the six species, suggesting that the genome evolution and polyploidization processes had relatively minor effects on the ACP gene family. The duplication modes of the six Brassica species were diverse, and the expansion of most ACPs in Brassica occurred primarily through dispersed duplication (DSD) events. Furthermore, most of the ACP genes were under purifying selection during the process of evolution. Subcellular localization experiments demonstrated that ACP genes in Brassica species are localized in chloroplasts and mitochondria. Cis-acting element analysis revealed that most of the ACP genes were associated with various abiotic stresses. Additionally, RNA-seq data revealed differential expression levels of BnaACP genes across various tissues in B. napus, with particularly high expression in seeds and buds. qRT-PCR analysis further indicated that BnaACP genes play a significant role in salt stress tolerance. These findings provide a comprehensive understanding of ACP genes in Brassica plants and will facilitate further functional analysis of these genes.

BnAP2-12 overexpression delays ramie flowering: evidence from AP2/ERF gene expression

Introduction

The APETALA2/ethylene response factor (AP2/ERF) superfamily plays a significant role in regulating plant gene expression in response to growth and development. To date, there have been no studies into whether the ramie AP2/ERF genes are involved in the regulation of flower development.

Methods

Here, 84 BnAP2/ERF members were identified from the ramie genome database, and various bioinformatics data on the AP2/ERF gene family, structure, replication, promoters and regulatory networks were analysed. BnAP2-12 was transferred into Arabidopsis through the flower-dipping method.

Results

Phylogenetic analysis classified the 84 BnAP2/ERF members into four subfamilies: AP2 (18), RAV (3), ERF (42), and DREB (21). The functional domain analysis of genes revealed 10 conserved motifs. Genetic mapping localised the 84 members on 14 chromosomes, among which chromosomes 1, 3, 5, and 8 had more members. Collinearity analysis revealed that 43.37% possibly resulted from replication events during the evolution of the ramie genome. Promoter sequence analysis identified classified cis-acting elements associated with plant growth and development, and responses to stress, hormones, and light. Transcriptomic comparison identified 3,635 differentially expressed genes (DEGs) between male and female flowers (1,803 and 1,832 upregulated and downregulated genes, respectively). Kyoto Encyclopaedia of Genes and Genomes pathway analysis categorised DEGs involved in metabolic pathways and biosynthesis of secondary metabolites. Gene Ontology enrichment analysis further identified enriched genes associated with pollen and female gamete formations. Of the 84 BnAP2/ERFs genes, 22 and 8 upregulated and downregulated genes, respectively, were present in female flowers. Co-expression network analysis identified AP2/ERF members associated with flower development, including BnAP2-12. Subcellular localisation analysis showed that the BnAP2-12 protein is localised in the nucleus and cell membrane. Overexpression BnAP2-12 delayed the flowering time of Arabidopsis thaliana.

Conclusion

These findings provide insights into the mechanism of ramie flower development.

InPACT: a computational method for accurate characterization of intronic polyadenylation from RNA sequencing data

Alternative polyadenylation can occur in introns, termed intronic polyadenylation (IPA), has been implicated in diverse biological processes and diseases, as it can produce noncoding transcripts or transcripts with truncated coding regions. However, a reliable method is required to accurately characterize IPA. Here, we propose a computational method called InPACT, which allows for the precise characterization of IPA from conventional RNA-seq data. InPACT successfully identifies numerous previously unannotated IPA transcripts in human cells, many of which are translated, as evidenced by ribosome profiling data. We have demonstrated that InPACT outperforms other methods in terms of IPA identification and quantification. Moreover, InPACT applied to monocyte activation reveals temporally coordinated IPA events. Further application on single-cell RNA-seq data of human fetal bone marrow reveals the expression of several IPA isoforms in a context-specific manner. Therefore, InPACT represents a powerful tool for the accurate characterization of IPA from RNA-seq data.

Cancer cell genetics shaping of the tumor microenvironment reveals myeloid cell-centric exploitable vulnerabilities in hepatocellular carcinoma

Myeloid cells are abundant and plastic immune cell subsets in the liver, to which pro-tumorigenic, inflammatory and immunosuppressive roles have been assigned in the course of tumorigenesis. Yet several aspects underlying their dynamic alterations in hepatocellular carcinoma (HCC) progression remain elusive, including the impact of distinct genetic mutations in shaping a cancer-permissive tumor microenvironment (TME). Here, in newly generated, clinically-relevant somatic female HCC mouse models, we identify cancer genetics' specific and stage-dependent alterations of the liver TME associated with distinct histopathological and malignant HCC features. Mitogen-activated protein kinase (MAPK)-activated, NrasG12D-driven tumors exhibit a mixed phenotype of prominent inflammation and immunosuppression in a T cell-excluded TME. Mechanistically, we report a NrasG12D cancer cell-driven, MEK-ERK1/2-SP1-dependent GM-CSF secretion enabling the accumulation of immunosuppressive and proinflammatory monocyte-derived Ly6Clow cells. GM-CSF blockade curbs the accumulation of these cells, reduces inflammation, induces cancer cell death and prolongs animal survival. Furthermore, GM-CSF neutralization synergizes with a vascular endothelial growth factor (VEGF) inhibitor to restrain HCC outgrowth. These findings underscore the profound alterations of the myeloid TME consequential to MAPK pathway activation intensity and the potential of GM-CSF inhibition as a myeloid-centric therapy tailored to subsets of HCC patients.

Comprehensive Analysis of Autophagy-Related Genes in Rice Immunity against Magnaporthe oryzae

Rice blast disease, caused by the fungus Magnaporthe oryzae, is a significant threat to rice production. Resistant cultivars can effectively resist the invasion of M. oryzae. Thus, the identification of disease-resistant genes is of utmost importance for improving rice production. Autophagy, a cellular process that recycles damaged components, plays a vital role in plant growth, development, senescence, stress response, and immunity. To understand the involvement of autophagy-related genes (ATGs) in rice immune response against M. oryzae, we conducted a comprehensive analysis of 37 OsATGs, including bioinformatic analysis, transcriptome analysis, disease resistance analysis, and protein interaction analysis. Bioinformatic analysis revealed that the promoter regions of 33 OsATGs contained cis-acting elements responsive to salicylic acid (SA) or jasmonic acid (JA), two key hormones involved in plant defense responses. Transcriptome data showed that 21 OsATGs were upregulated during M. oryzae infection. Loss-of-function experiments demonstrated that OsATG6c, OsATG8a, OsATG9b, and OsATG13a contribute to rice blast resistance. Additionally, through protein interaction analysis, we identified five proteins that may interact with OsATG13a and potentially contribute to plant immunity. Our study highlights the important role of autophagy in rice immunity and suggests that OsATGs may enhance resistance to rice blast fungus through the involvement of SA, JA, or immune-related proteins. These findings provide valuable insights for future efforts in improving rice production through the identification and utilization of autophagy-related genes.