The MEME Suite

The complete catalog of antimicrobial resistance secondary active transporters in Clostridioides difficile: evolution and drug resistance perspective

Secondary active transporters shuttle substrates across eukaryotic and prokaryotic membranes, utilizing different electrochemical gradients. They are recognized as one of the antimicrobial efflux pumps among pathogens. While primary active transporters within the genome of C. difficile 630 have been completely cataloged, the systematical study of secondary active transporters remains incomplete. Here, we not only identify secondary active transporters but also disclose their evolution and role in drug resistance in C. difficile 630. Our analysis reveals that C. difficile 630 carries 147 secondary active transporters belonging to 27 (super)families. Notably, 50 (34%) of them potentially contribute to antimicrobial resistance (AMR). AMR-secondary active transporters are structurally classified into five (super)families: the p-aminobenzoyl-glutamate transporter (AbgT), drug/metabolite transporter (DMT) superfamily, major facilitator (MFS) superfamily, multidrug and toxic compound extrusion (MATE) family, and resistance-nodulation-division (RND) family. Surprisingly, complete RND genes found in C. difficile 630 are likely an evolutionary leftover from the common ancestor with the diderm. Through protein structure comparisons, we have potentially identified six novel AMR-secondary active transporters from DMT, MATE, and MFS (super)families. Pangenome analysis revealed that half of the AMR-secondary transporters are accessory genes, which indicates an important role in adaptive AMR function rather than innate physiological homeostasis. Gene expression profile firmly supports their ability to respond to a wide spectrum of antibiotics. Our findings highlight the evolution of AMR-secondary active transporters and their integral role in antibiotic responses. This marks AMR-secondary active transporters as interesting therapeutic targets to synergize with other antibiotic activity.

Augusta: From RNA-Seq to gene regulatory networks and Boolean models

Computational models of gene regulations help to understand regulatory mechanisms and are extensively used in a wide range of areas, e.g., biotechnology or medicine, with significant benefits. Unfortunately, there are only a few computational gene regulatory models of whole genomes allowing static and dynamic analysis due to the lack of sophisticated tools for their reconstruction. Here, we describe Augusta, an open-source Python package for Gene Regulatory Network (GRN) and Boolean Network (BN) inference from the high-throughput gene expression data. Augusta can reconstruct genome-wide models suitable for static and dynamic analyses. Augusta uses a unique approach where the first estimation of a GRN inferred from expression data is further refined by predicting transcription factor binding motifs in promoters of regulated genes and by incorporating verified interactions obtained from databases. Moreover, a refined GRN is transformed into a draft BN by searching in the curated model database and setting logical rules to incoming edges of target genes, which can be further manually edited as the model is provided in the SBML file format. The approach is applicable even if information about the organism under study is not available in the databases, which is typically the case for non-model organisms including most microbes. Augusta can be operated from the command line and, thus, is easy to use for automated prediction of models for various genomes. The Augusta package is freely available at github.com/JanaMus/Augusta. Documentation and tutorials are available at augusta.readthedocs.io.

Genome-wide identification of the Hsp70 gene family in Penaeus chinensis and their response to environmental stress

The heat shock protein 70 (HSP70) gene family plays a crucial role in the response of organisms to environmental stress. However, it has not been systematically characterized in shrimp. In this study, we identified 25 PcHsp70 genes in the Penaeus chinensis genome. The encoded proteins were categorized into six subgroups based on phylogenetic relationships. Tandem duplication was the main driver of amplification in the PcHsp70 family, and the genes have experienced strong purifying selection during evolution. Transcriptome data analysis revealed that the 25 PcHsp70 members have different expression patterns in shrimp under conditions of low temperature, low salinity, and white spot syndrome virus infection. Among them, PcHsp70.11 was significantly induced under all three stress conditions, suggesting that this gene plays an important role in response to environmental stress in P. chinensis. To the best of our knowledge, this is the first study to systematically analyze the Hsp70 gene family in shrimp. The results provide important information on shrimp Hsp70s, contributing to a better understanding of the role of these genes in environmental stress and providing a basis for further functional studies.

The expansion of oligopeptide transporters in Melampsora larici-populina may reflect its adaptation to a phytoparasitic lifestyle

The acquisition of nutrients from host plants by phytopathogenic fungi is critically important for their invasion success. Melampsora larici-populina, an obligate biotrophic pathogenic fungus, causes the poplar leaf rust disease and can severely damage host poplar plants. Previously, we found that oligopeptide transporters (OPTs) have undergone a convergent expansion, which might reflect adaptation to a phytoparasitic lifestyle. Here, we used various methods to evaluate this hypothesis, including conserved motif identification, positive selection signal mining, expression pattern clustering analysis, and neutral selection tests. The motif composition of the five clades in the OPT family differed, and positive selection was observed during clade differentiation. This suggests that OPTs in these five clades may be functionally differentiated, which would increase the range of transported substrates and promote the absorption of more types of nitrogen compounds from the hosts. According to clustering analysis of gene expression patterns, the expression of most genes from the two expanded clades (clade 2 and 4) was up-regulated during the infection of poplar trees, indicating that the expansion of OPTs likely occurred to promote the uptake of oligopeptides from host poplar plants. The MellpOPT4g gene was determined to be under significant balancing selection based on the neutral selection tests, suggesting that it plays a role in the pathogenic process. In conclusion, these three observations provide preliminary evidence supporting our hypothesis, as they indicate that the expansion of OPTs in M. larici-populina has aided the ability of this pathogen to acquire nutrients from host plants.

Genome-wide identification and characterization of DIRIGENT gene family (CcDIR) in pigeonpea (Cajanus cajan L.) provide insights on their spatial expression pattern and relevance to stress response

This study is a thorough characterization of pigeonpea dirigent gene (CcDIR) family, an important component of the lignin biosynthesis pathway. Genome-wide analysis identified 25 CcDIR genes followed by a range of analytical approaches employed to unravel their structural and functional characteristics. Structural examination revealed a classic single exon and no intron arrangement in CcDIRs contributing to our understanding on evolutionary dynamics. Phylogenetic analysis elucidated evolutionary relationships among CcDIR genes with six DIR sub-families, while motif distribution analysis displayed and highlighted ten conserved protein motifs in CcDIRs. Promoter analyses of all the dirigent genes detected 18 stress responsive cis-acting elements offering insights into transcriptional regulation. While spatial expression analyses across six plant tissues showed preferential expression of CcDIR genes, exposure to salt (CcDIR2 and CcDIR9) and herbivory (CcDIR1, CcDIR2, CcDIR3 and CcDIR11), demonstrated potential roles of specific DIRs in plant defense. Interestingly, increased gene expression during herbivory, also correlated with increased lignin content authenticating the specific response. Furthermore, exogenous application of stress hormones, SA and MeJA on leaves significantly induced the expression of CcDIRs that responded to herbivory. Taken together, these findings contribute to a comprehensive understanding of CcDIR genes impacting development and stress response in the important legume pigeonpea.

Comprehensive identification and analysis of DUF640 genes associated with rice growth

Protein domains with conserved amino acid sequences and uncharacterized functions are called domains of unknown function (DUF). The DUF640 gene family plays a crucial role in plant growth, particularly in light regulation, floral organ development, and fruit development. However, there exists a lack of systematic understanding of the evolutionary relationships and functional differentiation of DUF640 within the Oryza genus. In this study, 61 DUF640 genes were identified in the Oryza genus. The expression of DUF640s is induced by multiple hormonal stressors including abscisic acid (ABA), cytokinin (CK), ethylene (ETH), and indole-3-acetic acid (IAA). Specifically, OiDUF640-10 expression significantly increased after ETH treatment. Transgenic experiments showed that overexpressing OiDUF640-10 lines were sensitive to ETH, and seedling length was obstructed. Evolutionary analysis revealed differentiation of the OiDUF640-10 gene in O. sativa ssp. indica and japonica varieties, likely driven by natural selection during the domestication of cultivated rice. These results indicate that OiDUF640-10 plays a vital role in the regulation of rice seedling length.

Dynamic DNA methylation modifications in the cold stress response of cassava

Cassava, a crucial tropical crop, faces challenges from cold stress, necessitating an exploration of its molecular response. Here, we investigated the role of DNA methylation in moderating the response to moderate cold stress (10 °C) in cassava. Using whole-genome bisulfite sequencing, we examined DNA methylation patterns in leaf blades and petioles under control conditions, 5 h, and 48 h of cold stress. Tissue-specific responses were observed, with leaf blades exhibiting subtle changes, while petioles displayed a pronounced decrease in methylation levels under cold stress. We identified cold stress-induced differentially methylated regions (DMRs) that demonstrated both tissue and treatment specificity. Importantly, these DMRs were enriched in genes with altered expression, implying functional relevance. The cold-response transcription factor ERF105 associated with DMRs emerged as a significant and conserved regulator across tissues and treatments. Furthermore, we investigated DNA methylation dynamics in transposable elements, emphasizing the sensitivity of MITEs with bHLH binding motifs to cold stress. These findings provide insights into the epigenetic regulation of response to cold stress in cassava, contributing to an understanding of the molecular mechanisms underlying stress adaptation in this tropical plant.

Transcriptional synergy in human aortic endothelial cells is vulnerable to combination p300/CBP and BET bromodomain inhibition

Combinatorial signaling by proinflammatory cytokines synergizes to exacerbate toxicity to cells and tissue injury during acute infections. To explore synergism at the gene-regulatory level, we investigated the dynamics of transcription and chromatin signaling in response to dual cytokines by integrating nascent RNA imaging mass spectrometry, RNA sequencing, amplification-independent mRNA quantification, assay for transposase-accessible chromatin using sequencing (ATAC-seq), and transcription factor profiling. Costimulation with interferon-gamma (IFNγ) and tumor necrosis factor alpha (TNFα) synergistically induced a small subset of genes, including the chemokines CXCL9, -10, and -11. Gene induction coincided with increased chromatin accessibility at non-coding regions enriched for p65 and STAT1 binding sites. To discover coactivator dependencies, we conducted a targeted chemogenomic screen of transcriptional inhibitors followed by modeling of inhibitor dose-response curves. These results identified high efficacy of either p300/CREB-binding protein (CBP) or bromodomain and extra-terminal (BET) bromodomain inhibitors to disrupt induction of synergy genes. Combination p300/CBP and BET bromodomain inhibition at half-maximal inhibitory concentrations (subIC50) synergistically abrogated IFNγ/TNFα-induced chemokine gene and protein levels.

Engineering of a TetR family transcriptional regulator BkdR enhances heterologous spinosad production in Streptomyces albus B4 chassis

Precursor supply plays a significant role in the production of secondary metabolites. In Streptomyces bacteria, propionyl-, malonyl-, and methylmalonyl-CoA are the most common precursors used for polyketide biosynthesis. Although propionyl-CoA synthetases participate in the propionate assimilation pathway and directly convert propionate into propionyl-CoA, malonyl- and methylmalonyl-CoA cannot be formed using common acyl-CoA synthetases. Therefore, both acetyl- and propionyl-CoA carboxylation, catalyzed by acyl-CoA carboxylases, should be considered when engineering a microorganism chassis to increase polyketide production. In this study, we identified a transcriptional regulator of the TetR family, BkdR, in Streptomyces albus B4, which binds directly to the promoter region of the neighboring pccAB operon. This operon encodes acetyl/propionyl-CoA carboxylase and negatively regulates its transcription. In addition to acetate and propionate, the binding of BkdR to pccAB is disrupted by acetyl- and propionyl-CoA ligands. We identified a 16-nucleotide palindromic BkdR-binding motif (GTTAg/CGGTCg/TTAAC) in the intergenic region between pccAB and bkdR. When bkdR was deleted, we found an enhanced supply of malonyl- and methylmalonyl-CoA precursors in S. albus B4. In this study, spinosad production was detected in the recombinant strain after introducing the entire artificial biosynthesized gene cluster into S. albus B4. When supplemented with propionate to provide propionyl-CoA, the novel bkdR-deleted strain produced 29.4% more spinosad than the initial strain in trypticase soy broth (TSB) medium.

IMPORTANCE

In this study, we describe a pccAB operon involved in short-chain acyl-CoA carboxylation in S. albus B4 chassis. The TetR family regulator, BkdR, represses this operon. Our results show that BkdR regulates the precursor supply needed for heterologous spinosad biosynthesis by controlling acetyl- and propionyl-CoA assimilation. The deletion of the BkdR-encoding gene exerts an increase in heterologous spinosad yield. Our research reveals a regulatory mechanism in short-chain acyl-CoA metabolism and suggests new possibilities for S. albus chassis engineering to enhance heterologous polyketide yield.

Genome-wide identification and characterization of SNW/SKIP domain-containing proteins in plants

Sessile organisms, such as plants, developed various ways to sense and respond to external and internal stimuli to maximize their fitness through evolutionary time. Transcripts and protein regulation are, among many, the main mechanisms that plants use to respond to environmental changes. SKIP protein is one such, presenting an SNKW interacting domain, which is highly conserved among eukaryotes, where SKI interacting protein acts in regulating key processes. In the present work, many bioinformatics tools, such as phylogenetic relationships, gene structure, physical-chemical properties, conserved motifs, prediction of regulatory cis-elements, chromosomal localization, and protein-protein interaction network, were used to better understand the genome-wide SNW/SKIP domain-containing proteins. In total, 28 proteins containing the SNW/SKIP domain were identified in different plant species, including plants of agronomic interest. Two main protein clusters were formed in phylogenetic analysis, and gene structure analysis revealed that, in general, the coding region had no introns. Also, expression of these genes is possibly induced by abiotic stress stimuli. Primary structure analysis of the proteins revealed the existence of an evolutionarily conserved functional unit. But physicochemical properties show that proteins containing the SNW/SKIP domain are commonly unstable under in vivo conditions. In addition, the protein network, demonstrated that SKIP homologues could act by modulating plant fitness through gene expression regulation at the transcriptional and post-transcriptional levels. This could be corroborated by the expression number of gene copies of SKIP proteins in many species, highlighting it's crucial role in plant development and tolerance through the course of evolution.

Age-dependent heat shock hormesis to HSF-1 deficiency suggests a compensatory mechanism mediated by the unfolded protein response and innate immunity in young Caenorhabditis elegans

The transcription factor HSF-1 (heat shock factor 1) acts as a master regulator of heat shock response in eukaryotic cells to maintain cellular proteostasis. The protein has a protective role in preventing cells from undergoing ageing, and neurodegeneration, and also mediates tumorigenesis. Thus, modulating HSF-1 activity in humans has a promising therapeutic potential for treating these pathologies. Loss of HSF-1 function is usually associated with impaired stress tolerance. Contrary to this conventional knowledge, we show here that inactivation of HSF-1 in the nematode Caenorhabditis elegans results in increased thermotolerance at young adult stages, whereas HSF-1 deficiency in animals passing early adult stages indeed leads to decreased thermotolerance, as compared to wild-type. Furthermore, a gene expression analysis supports that in young adults, distinct cellular stress response and immunity-related signaling pathways become induced upon HSF-1 deficiency. We also demonstrate that increased tolerance to proteotoxic stress in HSF-1-depleted young worms requires the activity of the unfolded protein response of the endoplasmic reticulum and the SKN-1/Nrf2-mediated oxidative stress response pathway, as well as an innate immunity-related pathway, suggesting a mutual compensatory interaction between HSF-1 and these conserved stress response systems. A similar compensatory molecular network is likely to also operate in higher animal taxa, raising the possibility of an unexpected outcome when HSF-1 activity is manipulated in humans.

Transcriptomic profiling of haloarchaeal denitrification through RNA-Seq analysis

Denitrification, a crucial biochemical pathway prevalent among haloarchaea in hypersaline ecosystems, has garnered considerable attention in recent years due to its ecological implications. Nevertheless, the underlying molecular mechanisms and genetic regulation governing this respiration/detoxification process in haloarchaea remain largely unexplored. In this study, RNA-sequencing was used to compare the transcriptomes of the haloarchaeon Haloferax mediterranei under oxic and denitrifying conditions, shedding light on the intricate metabolic alterations occurring within the cell, such as the accurate control of the metal homeostasis. Furthermore, the investigation identifies several genes encoding transcriptional regulators and potential accessory proteins with putative roles in denitrification. Among these are bacterioopsin-like transcriptional activators, proteins harboring a domain of unknown function (DUF2249), and cyanoglobin. In addition, the study delves into the genetic regulation of denitrification, finding a regulatory motif within promoter regions that activates numerous denitrification-related genes. This research serves as a starting point for future molecular biology studies in haloarchaea, offering a promising avenue to unravel the intricate mechanisms governing haloarchaeal denitrification, a pathway of paramount ecological importance.IMPORTANCEDenitrification, a fundamental process within the nitrogen cycle, has been subject to extensive investigation due to its close association with anthropogenic activities, and its contribution to the global warming issue, mainly through the release of N2O emissions. Although our comprehension of denitrification and its implications is generally well established, most studies have been conducted in non-extreme environments with mesophilic microorganisms. Consequently, there is a significant knowledge gap concerning extremophilic denitrifiers, particularly those inhabiting hypersaline environments. The significance of this research was to delve into the process of haloarchaeal denitrification, utilizing the complete denitrifier haloarchaeon Haloferax mediterranei as a model organism. This research led to the analysis of the metabolic state of this microorganism under denitrifying conditions and the identification of regulatory signals and genes encoding proteins potentially involved in this pathway, serving as a valuable resource for future molecular studies.

Cancer-associated Histone H3 N-terminal arginine mutations disrupt PRC2 activity and impair differentiation

Dysregulated epigenetic states are a hallmark of cancer and often arise from genetic alterations in epigenetic regulators. This includes missense mutations in histones, which, together with associated DNA, form nucleosome core particles. However, the oncogenic mechanisms of most histone mutations are unknown. Here, we demonstrate that cancer-associated histone mutations at arginines in the histone H3 N-terminal tail disrupt repressive chromatin domains, alter gene regulation, and dysregulate differentiation. We find that histone H3R2C and R26C mutants reduce transcriptionally repressive H3K27me3. While H3K27me3 depletion in cells expressing these mutants is exclusively observed on the minor fraction of histone tails harboring the mutations, the same mutants recurrently disrupt broad H3K27me3 domains in the chromatin context, including near developmentally regulated promoters. H3K27me3 loss leads to de-repression of differentiation pathways, with concordant effects between H3R2 and H3R26 mutants despite different proximity to the PRC2 substrate, H3K27. Functionally, H3R26C-expressing mesenchymal progenitor cells and murine embryonic stem cell-derived teratomas demonstrate impaired differentiation. Collectively, these data show that cancer-associated H3 N-terminal arginine mutations reduce PRC2 activity and disrupt chromatin-dependent developmental functions, a cancer-relevant phenotype.

Phage-host co-evolution has led to distinct generalized transduction strategies

Generalized transduction is pivotal in bacterial evolution but lacks comprehensive understanding regarding the facilitating features and variations among phages. We addressed this gap by sequencing and comparing the transducing particle content of three different Salmonella Typhimurium phages (i.e. Det7, ES18 and P22) that share a headful packaging mechanism that is typically initiated from a cognate pac site within the phage chromosome. This revealed substantial disparities in both the extent and content of transducing particles among these phages. While Det7 outperformed ES18 in terms of relative number of transducing particles, both phages contrasted with P22 in terms of content. In fact, we found evidence for the presence of conserved P22 pac-like sequences in the host chromosome that direct tremendously increased packaging and transduction frequencies of downstream regions by P22. More specifically, a ca. 561 kb host region between oppositely oriented pac-like sequences in the purF and minE loci was identified as highly packaged and transduced during both P22 prophage induction and lytic infection. Our findings underscore the evolution of phage transducing capacity towards attenuation, promiscuity or directionality, and suggest that pac-like sequences in the host chromosome could become selected as sites directing high frequency of transduction.

Zinc-finger (ZiF) fold secreted effectors form a functionally diverse family across lineages of the blast fungus Magnaporthe oryzae

Filamentous plant pathogens deliver effector proteins into host cells to suppress host defence responses and manipulate metabolic processes to support colonization. Understanding the evolution and molecular function of these effectors provides knowledge about pathogenesis and can suggest novel strategies to reduce damage caused by pathogens. However, effector proteins are highly variable, share weak sequence similarity and, although they can be grouped according to their structure, only a few structurally conserved effector families have been functionally characterized to date. Here, we demonstrate that Zinc-finger fold (ZiF) secreted proteins form a functionally diverse effector family in the blast fungus Magnaporthe oryzae. This family relies on the Zinc-finger motif for protein stability and is ubiquitously present in blast fungus lineages infecting 13 different host species, forming different effector tribes. Homologs of the canonical ZiF effector, AVR-Pii, from rice infecting isolates are present in multiple M. oryzae lineages. Wheat infecting strains of the fungus also possess an AVR-Pii like allele that binds host Exo70 proteins and activates the immune receptor Pii. Furthermore, ZiF tribes may vary in the proteins they bind to, indicating functional diversification and an intricate effector/host interactome. Altogether, we uncovered a new effector family with a common protein fold that has functionally diversified in lineages of M. oryzae. This work expands our understanding of the diversity of M. oryzae effectors, the molecular basis of plant pathogenesis and may ultimately facilitate the development of new sources for pathogen resistance.

Genome-wide identification of starch phosphorylase gene family in Rosa chinensis and expression in response to abiotic stress

Chinese rose (Rosa chinensis) is an important ornamental plant, with economic, cultural, and symbolic significance. During the application of outdoor greening, adverse environments such as high temperature and drought are often encountered, which affect its application scope and ornamental quality. The starch phosphorylase (Pho) gene family participate in the synthesis and decomposition of starch, not only related to plant energy metabolism, but also plays an important role in plant stress resistance. The role of Pho in combating salinity and high temperature stress in R. chinensis remains unknown. In this work, 4 Phos from R. chinensis were detected with Pfam number of Pho (PF00343.23) and predicted by homolog-based prediction (HBP). The Phos are characterized by sequence lengths of 821 to 997 bp, and the proteins are predicted to subcellularly located in the plastid and cytoplasm. The regulatory regions of the Phos contain abundant stress and phytohormone-responsive cis-acting elements. Based on transcriptome analysis, the Phos were found to respond to abiotic stress factors such as drought, salinity, high temperature, and plant phytohormone of jasmonic acid and salicylic acid. The response of Phos to abiotic stress factors such as salinity and high temperature was confirmed by qRT-PCR analysis. To evaluate the genetic characteristics of Phos, a total of 69 Phos from 17 species were analyzed and then classified into 3 groups in phylogenetic tree. The collinearity analysis of Phos in R. chinensis and other species was conducted for the first time. This work provides a view of evolution for the Pho gene family and indicates that Phos play an important role in abiotic stress response of R. chinensis.

Genome-wide analysis of AP2/ERF transcription factors that regulate fruit development of Chinese prickly ash

BACKGROUND

AP2/ERF is a large family of plant transcription factor proteins that play essential roles in signal transduction, plant growth and development, and responses to various stresses. The AP2/ERF family has been identified and verified by functional analysis in various plants, but so far there has been no comprehensive study of these factors in Chinese prickly ash. Phylogenetic, motif, and functional analyses combined with transcriptome analysis of Chinese prickly ash fruits at different developmental stages (30, 60, and 90 days after anthesis) were conducted in this study.

RESULTS

The analysis identified 146 ZbAP2/ERF genes that could be classified into 15 subgroups. The motif analysis revealed the presence of different motifs or elements in each group that may explain the functional differences between the groups. ZbERF13.2, ZbRAP2-12, and ZbERF2.1 showed high levels of expression in the early stages of fruit development. ZbRAP2-4, and ZbERF3.1 were significantly expressed at the fruit coloring stage (R2 and G2). ZbERF16 were significantly expressed at fruit ripening and expression level increased as the fruit continued to develop. Relative gene expression levels of 6 representative ZbAP2/ERFs assessed by RT-qPCR agreed with transcriptome analysis results.

CONCLUSIONS

These genes identified by screening can be used as candidate genes that affect fruit development. The results of the analysis can help guide future genetic improvement of Chinese prickly ash and enrich our understanding of AP2/ERF transcription factors and their regulatory functions in plants.

Genome and transcriptome exploration reveals receptor-like kinases as potential resistance gene analogs against bacterial blight in pomegranate

BACKGROUND

Pomegranate (Punica granatum L.) is a tropical fruit crop of pharma-nutritional importance. However, it faces farming challenges due to pests and diseases, particularly bacterial blight and wilt. Developing resistant cultivars is crucial for sustainable pomegranate cultivation, and understanding resistance's genetic basis is essential.

METHODS AND RESULTS

We used an extensive resistance gene analogues (RGA) prediction tool to identify 958 RGAs, classified into Nucleotide Binding Site-leucine-rich repeat (NBS-LRR) proteins, receptor-like kinases (RLKs), receptor-like proteins (RLPs), Transmembrane coiled-coil (TM-CC), and nine non-canonical RGAs. RGAs were distributed across all eight chromosomes, with chromosome 02 containing the most RGAs (161), and chromosome 08 having the highest density (4.42 RGA/Mb). NBS-LRR genes were predominantly present on chromosomes 08 and 02, whereas RLKs and RLPs were primarily located on chromosomes 04 and 07. Gene ontology analysis revealed that 475 RGAs were associated with defence against various biotic stresses. Using RNAseq, we identified 120 differentially expressed RGAs, with RLKs (74) being prominent among the differentially expressed genes.

CONCLUSION

The discovery of these RGAs is a significant step towards breeding pomegranates for pest and disease resistance. The differentially expressed RLKs hold promise for developing resistant cultivars against bacterial blight, thereby contributing to the sustainability of pomegranate cultivation.

Transcriptomic analysis reveals the regulatory mechanisms of messenger RNA (mRNA) and long non-coding RNA (lncRNA) in response to waterlogging stress in rye (Secale cereale L.)

BACKGROUND

Waterlogging stress (WS) negatively impacts crop growth and productivity, making it important to understand crop resistance processes and discover useful WS resistance genes. In this study, rye cultivars and wild rye species were subjected to 12-day WS treatment, and the cultivar Secale cereale L. Imperil showed higher tolerance. Whole transcriptome sequencing was performed on this cultivar to identify differentially expressed (DE) messenger RNAs (DE-mRNAs) and long non-coding RNAs (DE-lncRNAs) involved in WS response.

RESULTS

Among the 6 species, Secale cereale L. Imperil showed higher tolerance than wild rye species against WS. The cultivar effectively mitigated oxidative stress, and regulated hydrogen peroxide and superoxide anion. A total of 728 DE-mRNAs and 60 DE-lncRNAs were discovered. Among these, 318 DE-mRNAs and 32 DE-lncRNAs were upregulated, and 410 DE-mRNAs and 28 DE-lncRNAs were downregulated. GO enrichment analysis discovered metabolic processes, cellular processes, and single-organism processes as enriched biological processes (BP). For cellular components (CC), the enriched terms were membrane, membrane part, cell, and cell part. Enriched molecular functions (MF) terms were catalytic activity, binding, and transporter activity. LncRNA and mRNA regulatory processes were mainly related to MAPK signaling pathway-plant, plant hormone signal transduction, phenylpropanoid biosynthesis, anthocyanin biosynthesis, glutathione metabolism, ubiquitin-mediated proteolysis, ABC transporter, Cytochrome b6/f complex, secondary metabolite biosynthesis, and carotenoid biosynthesis pathways. The signalling of ethylene-related pathways was not mainly dependent on AP2/ERF and WRKY transcription factors (TF), but on other factors. Photosynthetic activity was active, and carotenoid levels increased in rye under WS. Sphingolipids, the cytochrome b6/f complex, and glutamate are involved in rye WS response. Sucrose transportation was not significantly inhibited, and sucrose breakdown occurs in rye under WS.

CONCLUSIONS

This study investigated the expression levels and regulatory functions of mRNAs and lncRNAs in 12-day waterlogged rye seedlings. The findings shed light on the genes that play a significant role in rye ability to withstand WS. The findings from this study will serve as a foundation for further investigations into the mRNA and lncRNA WS responses in rye.

mRNA cap-binding protein eIF4E1 is a novel regulator of Toxoplasma gondii latency

The protozoan parasite Toxoplasma gondii causes serious opportunistic disease due to its ability to persist in patients as latent tissue cysts. The molecular mechanisms coordinating conversion between proliferative parasites (tachyzoites) and latent cysts (bradyzoites) are not fully understood. We previously showed that phosphorylation of eIF2α accompanies bradyzoite formation, suggesting that this clinically relevant process involves regulation of mRNA translation. In this study, we investigated the composition and role of eIF4F multi-subunit complexes in translational control. Using CLIPseq, we find that the cap-binding subunit, eIF4E1, localizes to the 5'-end of all tachyzoite mRNAs, many of which show evidence of stemming from heterogeneous transcriptional start sites. We further show that eIF4E1 operates as the predominant cap-binding protein in two distinct eIF4F complexes. Using genetic and pharmacological approaches, we found that eIF4E1 deficiency triggers efficient spontaneous formation of bradyzoites without stress induction. Consistent with this result, we also show that stress-induced bradyzoites exhibit reduced eIF4E1 expression. Overall, our findings establish a novel role for eIF4F in translational control required for parasite latency and microbial persistence.

IMPORTANCE

Toxoplasma gondii is an opportunistic pathogen important to global human and animal health. There are currently no chemotherapies targeting the encysted form of the parasite. Consequently, a better understanding of the mechanisms controlling encystation is required. Here we show that the mRNA cap-binding protein, eIF4E1, regulates the encystation process. Encysted parasites reduce eIF4E1 levels, and depletion of eIF4E1 decreases the translation of ribosome-associated machinery and drives Toxoplasma encystation. Together, these data reveal a new layer of mRNA translational control that regulates parasite encystation and latency.

Prokaryotic expression of DFP1 and DFP2 in Dermatophagoides farinae and their responses to temperature stress

The functions of highly expressed genes DFP1 and DFP2 in Dermatophagoides farinae remain unknown. DFP1 and DFP2 have been abundantly annotated and were up-regulated under temperature stress at 43 °C and -10 °C in our previous RNA-seq study, indicating that DFP1 and DFP2 may have temperature stress response function. Here, we amplified, cloned, and sequenced to obtain the complete coding sequences of DFP1 and DFP2 and predicted their protein characteristics using bioinformatics analysis. Then, prokaryotic expression systems were constructed and found that DFP1 was expressed in Escherichia coli Rosetta-gami 2 (DE3) but not BL21 (DE3); DFP2 was expressed in both BL21 (DE3) and Rosetta-gami 2 (DE3), with higher expression in BL21 (DE3). Finally, the growth curves of bacteria were drawn and indicated that the DFP1- and DFP2-pET32a carrying recombinant bacteria grew better than the respectiveonly pET32a carrying control bacteria after heat and cold stress. This study confirms for the first time that DFP1 and DFP2 respond to temperature stress at the protein level. The constructed prokaryotic expression systems will provide an experimental foundation for future antibody preparation for western blotting detection to confirm the temperature-stress response functions of DFP1 and DFP2.

Control of endothelial cell function and arteriogenesis by MEG3:EZH2 epigenetic regulation of integrin expression

Epigenetic processes involving long non-coding RNAs regulate endothelial gene expression. However, the underlying regulatory mechanisms causing endothelial dysfunction remain to be elucidated. Enhancer of zeste homolog 2 (EZH2) is an important rheostat of histone H3K27 trimethylation (H3K27me3) that represses endothelial targets, but EZH2 RNA binding capacity and EZH2:RNA functional interactions have not been explored in post-ischemic angiogenesis. We used formaldehyde/UV-assisted crosslinking ligation and sequencing of hybrids and identified a new role for maternally expressed gene 3 (MEG3). MEG3 formed the predominant RNA:RNA hybrid structures in endothelial cells. Moreover, MEG3:EZH2 assists recruitment onto chromatin. By EZH2-chromatin immunoprecipitation, following MEG3 depletion, we demonstrated that MEG3 controls recruitment of EZH2/H3K27me3 onto integrin subunit alpha4 (ITGA4) promoter. Both MEG3 knockdown or EZH2 inhibition (A-395) promoted ITGA4 expression and improved endothelial cell migration and adhesion to fibronectin in vitro. The A-395 inhibitor re-directed MEG3-assisted chromatin remodeling, offering a direct therapeutic benefit by increasing endothelial function and resilience. This approach subsequently increased the expression of ITGA4 in arterioles following ischemic injury in mice, thus promoting arteriogenesis. Our findings show a context-specific role for MEG3 in guiding EZH2 to repress ITGA4. Novel therapeutic strategies could antagonize MEG3:EZH2 interaction for pre-clinical studies.

Massively parallel dissection of RNA in RNA-protein interactions in vivo

Many of the biological functions performed by RNA are mediated by RNA-binding proteins (RBPs), and understanding the molecular basis of these interactions is fundamental to biology. Here, we present massively parallel RNA assay combined with immunoprecipitation (MPRNA-IP) for in vivo high-throughput dissection of RNA-protein interactions and describe statistical models for identifying RNA domains and parsing the structural contributions of RNA. By using custom pools of tens of thousands of RNA sequences containing systematically designed truncations and mutations, MPRNA-IP is able to identify RNA domains, sequences, and secondary structures necessary and sufficient for protein binding in a single experiment. We show that this approach is successful for multiple RNAs of interest, including the long noncoding RNA NORAD, bacteriophage MS2 RNA, and human telomerase RNA, and we use it to interrogate the hitherto unknown sequence or structural RNA-binding preferences of the DNA-looping factor CTCF. By integrating systematic mutation analysis with crosslinking immunoprecipitation, MPRNA-IP provides a novel high-throughput way to elucidate RNA-based mechanisms behind RNA-protein interactions in vivo.

Transcription factor shapes chromosomal conformation and regulates gene expression in bacterial adaptation

Genomic mutations allow bacteria to adapt rapidly to adverse stress environments. The three-dimensional conformation of the genome may also play an important role in transcriptional regulation and environmental adaptation. Here, using chromosome conformation capture, we investigate the high-order architecture of the Zymomonas mobilis chromosome in response to genomic mutation and ambient stimuli (acetic acid and furfural, derived from lignocellulosic hydrolysate). We find that genomic mutation only influences the local chromosome contacts, whereas stress of acetic acid and furfural restrict the long-range contacts and significantly change the chromosome organization at domain scales. Further deciphering the domain feature unveils the important transcription factors, Ferric uptake regulator (Fur) proteins, which act as nucleoid-associated proteins to promote long-range (>200 kb) chromosomal communications and regulate the expression of genes involved in stress response. Our work suggests that ubiquitous transcription factors in prokaryotes mediate chromosome organization and regulate stress-resistance genes in bacterial adaptation.

ProEnd: A Comprehensive Database for Identifying HbYX Motif-Containing Proteins Across the Tree of Life

The proteasome plays a crucial role in cellular homeostasis by degrading misfolded, damaged, or unnecessary proteins. Understanding the regulatory mechanisms of proteasome activity is vital, particularly the interaction with activators containing the hydrophobic-tyrosine-any amino acid (HbYX) motif. Here, we present ProEnd, a comprehensive database designed to identify and catalog HbYX motif-containing proteins across the tree of life. Using a simple bioinformatics pipeline, we analyzed approximately 73 million proteins from 22,000 reference proteomes in the UniProt/SwissProt database. Our findings reveal the widespread presence of HbYX motifs in diverse organisms, highlighting their evolutionary conservation and functional significance. Notably, we observed an interesting prevalence of these motifs in viral proteomes, suggesting strategic interactions with the host proteasome. As validation two novel HbYX proteins found in this database were tested and found to directly interact with the proteasome. ProEnd's extensive dataset and user-friendly interface enable researchers to explore the potential proteasomal regulator landscape, generating new hypotheses to advance proteasome biology. This resource is set to facilitate the discovery of novel therapeutic targets, enhancing our approach to treating diseases such as neurodegenerative disorders and cancer. Link: http://proend.org/.

M. tuberculosis PrrA binds the dosR promoter and regulates mycobacterial adaptation to hypoxia

The PrrAB two-component system (TCS) is essential for Mycobacterium tuberculosis viability. Previously, it was demonstrated that PrrA binds DNA in the absence of PrrB-mediated transphosphorylation and that non-cognate serine/threonine-kinases phosphorylate PrrA threonine-6 (T6). Therefore, we investigated the differential binding affinity and regulatory properties of the M. tuberculosis-derived wild-type PrrA, PrrA phosphomimetic (D58E, T6E), and PrrA phosphoablative (D58A, T6A) proteins with the prrAMtb, dosRMtb, and cydAMtb genes. While we hypothesized greater DNA binding affinity and more pronounced regulation by PrrA phosphomimetic variants, recombinant, wild-type PrrAMtb bound DNA with greatest affinity. Collectively, wild-type PrrAMtb recombinant protein displayed the highest binding affinity to the dosRMtb promoter (KD 3.46 ± 2.09 nM), followed by the prrAMtb promoter (KD 9.00 ± 2.66 nM). To establish PrrAMtb regulatory activity, we constructed M. smegmatis ΔprrABMsmeg::prrAMtb strains with each of the PrrAMtb variants and extrachromosomal prrAMtb, dosRMtb, and cydAMtb promoter-mCherry reporter fusions. Our findings showed that PrrAMtb is autoregulatory and induces dosRMtb expression only during in vitro, hypoxic growth. Combined expression of prrABMtb in M. smegmatis ΔprrAB significantly induced cydAMtb promoter-mCherry expression. Our studies advanced the understanding of PrrA function and PrrAB phosphorylation-mediated regulatory mechanisms and control of mycobacterial dosR and cydA hypoxic and low-oxygen responsive genes.

Genome-wide identification, characterization, molecular evolution and expression profiling analysis of scavenger receptors in black rockfish (Sebastes schlegelii)

The scavenger receptors (SRs) gene family is considered as the membrane-associated pattern recognition receptors that plays important roles in the immune responses of organisms. However, there is currently limited research on the systematic identification of the SRs gene family in teleost and their role in the innate immunity of S. schegelii. In this study, we identified and annotated 15 SRs genes in S. schegelii. Through phylogenetic analysis, analysis of conserved domains, gene structure, and motif composition, we found that SRs gene family within different classes were relatively conserved. Additionally, we used qRT-PCR to analyze the expression patterns of SRs genes in immune-related tissues from healthy and Acinetobacter johnsonii-infected S. schegelii. The results showed that SRs genes exhibited different tissue expression patterns and the expression of SRs genes significantly changed after A. johnsonii infection. These results provided a valuable basis for further understanding of the functions of SRs in the innate immune response of S. schegelii.

Defining the KRAS- and ERK-dependent transcriptome in KRAS-mutant cancers

How the KRAS oncogene drives cancer growth remains poorly understood. Therefore, we established a systemwide portrait of KRAS- and extracellular signal-regulated kinase (ERK)-dependent gene transcription in KRAS-mutant cancer to delineate the molecular mechanisms of growth and of inhibitor resistance. Unexpectedly, our KRAS-dependent gene signature diverges substantially from the frequently cited Hallmark KRAS signaling gene signature, is driven predominantly through the ERK mitogen-activated protein kinase (MAPK) cascade, and accurately reflects KRAS- and ERK-regulated gene transcription in KRAS-mutant cancer patients. Integration with our ERK-regulated phospho- and total proteome highlights ERK deregulation of the anaphase promoting complex/cyclosome (APC/C) and other components of the cell cycle machinery as key processes that drive pancreatic ductal adenocarcinoma (PDAC) growth. Our findings elucidate mechanistically the critical role of ERK in driving KRAS-mutant tumor growth and in resistance to KRAS-ERK MAPK targeted therapies.

Tissue-resident trained immunity in hepatocytes protects against septic liver injury in zebrafish

Trained immunity is classically characterized by long-term functional reprogramming of innate immune cells to combat infectious diseases. Infection-induced organ injury is a common clinical severity phenotype of sepsis. However, whether the induction of trained immunity plays a role in protecting septic organ injury remains largely unknown. Here, through establishing an in vivo β-glucan training and lipopolysaccharide (LPS) challenge model in zebrafish larvae, we observe that induction of trained immunity could inhibit pyroptosis of hepatocytes to alleviate septic liver injury, with an elevated trimethyl-histone H3 lysine 4 (H3K4me3) modification that targets mitophagy-related genes. Moreover, we identify a C-type lectin domain receptor in zebrafish, named DrDectin-1, which is revealed as the orchestrator in gating H3K4me3 rewiring-mediated mitophagy activation and alleviating pyroptosis-engaged septic liver injury in vivo. Taken together, our results uncover tissue-resident trained immunity in maintaining liver homeostasis at the whole-animal level and offer an in vivo model to efficiently integrate trained immunity for immunotherapies.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Paleolithic Gene Duplications Primed Adaptive Evolution of Human Amylase Locus Upon Agriculture

Starch digestion is a cornerstone of human nutrition. The amylase genes code for the starch-digesting amylase enzyme. Previous studies suggested that the salivary amylase (AMY1) gene copy number increased in response to agricultural diets. However, the lack of nucleotide resolution of the amylase locus hindered detailed evolutionary analyses. Here, we have resolved this locus at nucleotide resolution in 98 present-day humans and identified 30 distinct haplotypes, revealing that the coding sequences of all amylase gene copies are evolving under negative selection. The phylogenetic reconstruction suggested that haplotypes with three AMY1 gene copies, prevalent across all continents and constituting about 70% of observed haplotypes, originated before the out-of-Africa migrations of ancestral modern humans. Using thousands of unique 25 base pair sequences across the amylase locus, we showed that additional AMY1 gene copies existed in the genomes of four archaic hominin genomes, indicating that the initial duplication of this locus may have occurred as far back 800,000 years ago. We similarly analyzed 73 ancient human genomes dating from 300 - 45,000 years ago and found that the AMY1 copy number variation observed today existed long before the advent of agriculture (~10,000 years ago), predisposing this locus to adaptive increase in the frequency of higher amylase copy number with the spread of agriculture. Mechanistically, the common three-copy haplotypes seeded non-allelic homologous recombination events that appear to be occurring at one of the fastest rates seen for tandem repeats in the human genome. Our study provides a comprehensive population-level understanding of the genomic structure of the amylase locus, identifying the mechanisms and evolutionary history underlying its duplication and copy number variability in relation to the onset of agriculture.

Genome-wide identification and characterization of DTX family genes highlighting their locations, functions, and regulatory factors in banana (Musa acuminata)

The detoxification efflux carriers (DTX) are a significant group of multidrug efflux transporter family members that play diverse functions in all kingdoms of living organisms. However, genome-wide identification and characterization of DTX family transporters have not yet been performed in banana, despite its importance as an economic fruit plant. Therefore, a detailed genome-wide analysis of DTX family transporters in banana (Musa acuminata) was conducted using integrated bioinformatics and systems biology approaches. In this study, a total of 37 DTX transporters were identified in the banana genome and divided into four groups (I, II, III, and IV) based on phylogenetic analysis. The gene structures, as well as their proteins' domains and motifs, were found to be significantly conserved. Gene ontology (GO) annotation revealed that the predicted DTX genes might play a vital role in protecting cells and membrane-bound organelles through detoxification mechanisms and the removal of drug molecules from banana cells. Gene regulatory analyses identified key transcription factors (TFs), cis-acting elements, and post-transcriptional regulators (miRNAs) of DTX genes, suggesting their potential roles in banana. Furthermore, the changes in gene expression levels due to pathogenic infections and non-living factor indicate that banana DTX genes play a role in responses to both biotic and abiotic stresses. The results of this study could serve as valuable tools to improve banana quality by protecting them from a range of environmental stresses.

DNA methylation patterns of transcription factor binding regions characterize their functional and evolutionary contexts

BACKGROUND

DNA methylation is an important epigenetic modification which has numerous roles in modulating genome function. Its levels are spatially correlated across the genome, typically high in repressed regions but low in transcription factor (TF) binding sites and active regulatory regions. However, the mechanisms establishing genome-wide and TF binding site methylation patterns are still unclear.

RESULTS

Here we use a comparative approach to investigate the association of DNA methylation to TF binding evolution in mammals. Specifically, we experimentally profile DNA methylation and combine this with published occupancy profiles of five distinct TFs (CTCF, CEBPA, HNF4A, ONECUT1, FOXA1) in the liver of five mammalian species (human, macaque, mouse, rat, dog). TF binding sites are lowly methylated, but they often also have intermediate methylation levels. Furthermore, biding sites are influenced by the methylation status of CpGs in their wider binding regions even when CpGs are absent from the core binding motif. Employing a classification and clustering approach, we extract distinct and species-conserved patterns of DNA methylation levels at TF binding regions. CEBPA, HNF4A, ONECUT1, and FOXA1 share the same methylation patterns, while CTCF's differ. These patterns characterize alternative functions and chromatin landscapes of TF-bound regions. Leveraging our phylogenetic framework, we find DNA methylation gain upon evolutionary loss of TF occupancy, indicating coordinated evolution. Furthermore, each methylation pattern has its own evolutionary trajectory reflecting its genomic contexts.

CONCLUSIONS

Our epigenomic analyses indicate a role for DNA methylation in TF binding changes across species including that specific DNA methylation profiles characterize TF binding and are associated with their regulatory activity, chromatin contexts, and evolutionary trajectories.

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

Knockdown of NFIC Promotes Bovine Myoblast Proliferation through the CENPF/CDK1 Axis

As cellular transcription factors and DNA replicators, nuclear factor I (NFI) family members play an important role in mammalian development. However, there is still a lack of research on the muscle regeneration of NFI family members in cattle. In this study, the analysis of NFI family factors was conducted on their characterization, phylogenetics, and functional domains. We found that NFI family members were relatively conserved among different species, but there was heterogeneity in amino acid sequences, DNA coding sequences, and functional domain among members. Furthermore, among NFI family factors, we observed that NFIC exhibited highly expression in bovine muscle tissues, particularly influencing the expression of proliferation marker genes in myoblasts. To investigate the influence of NFIC on myoblast proliferation, we knocked down NFIC (si-NFIC) and found that the proliferation of myoblasts was significantly promoted. In terms of regulation mechanism, we identified that si-NFIC could counteract the inhibitory effect of the cell cycle inhibitor RO-3306. Interestingly, CENPF, as the downstream target gene of NFIC, could affect the expression of CDK1, CCNB1, and actively regulate the cell cycle pathway and cell proliferation. In addition, when CENPF was knocked down, the phosphorylation of p53 and the expression of Bax were increased, but the expression of Bcl2 was inhibited. Our findings mainly highlight the mechanism by which NFIC acts on the CENPF/CDK1 axis to regulate the proliferation of bovine myoblasts.

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

Systematic Characterization of p53-Regulated Long Noncoding RNAs across Human Cancers Reveals Remarkable Heterogeneity among Different Tumor Types

The p53 tumor suppressor protein, a sequence-specific DNA binding transcription factor, regulates the expression of a large number of genes, in response to various forms of cellular stress. Although the protein coding target genes of p53 have been well studied, less is known about its role in regulating long noncoding genes and their functional relevance to cancer. Here we report the genome-wide identification of a large set (>1,000) of long noncoding RNAs (lncRNA), which are putative p53 targets in a colon cancer cell line and in human patient datasets from five different common types of cancer. These lncRNAs have not been annotated by other studies of normal unstressed systems. In the colon cancer cell line, a high proportion of these lncRNAs are uniquely induced by different chemotherapeutic agents that activate p53, whereas others are induced by more than one agent tested. Further, subsets of these lncRNAs independently predict overall and disease-free survival of patients across the five different common cancer types. Interestingly, both genetic alterations and patient survival associated with different lncRNAs are unique to each cancer tested, indicating extraordinary tissue-specific variability in the p53 noncoding response. The newly identified noncoding p53 target genes have allowed us to construct a classifier for tumor diagnosis and prognosis.

IMPLICATIONS

Our results not only identify myriad p53-regulated long noncoding (lncRNA), they also reveal marked drug-induced, as well as tissue- and tumor-specific heterogeneity in these putative p53 targets and our findings have enabled the construction of robust classifiers for diagnosis and prognosis.

Phosphorothioate DNA modification by BREX Type 4 systems in the human gut microbiome

Among dozens of microbial DNA modifications regulating gene expression and host defense, phosphorothioation (PT) is the only known backbone modification, with sulfur inserted at a non-bridging oxygen by dnd and ssp gene families. Here we explored the distribution of PT genes in 13,663 human gut microbiome genomes, finding that 6.3% possessed dnd or ssp genes predominantly in Bacillota, Bacteroidota, and Pseudomonadota. This analysis uncovered several putative new PT synthesis systems, including Type 4 Bacteriophage Exclusion (BREX) brx genes, which were genetically validated in Bacteroides salyersiae. Mass spectrometric analysis of DNA from 226 gut microbiome isolates possessing dnd , ssp , and brx genes revealed 8 PT dinucleotide settings confirmed in 6 consensus sequences by PT-specific DNA sequencing. Genomic analysis showed PT enrichment in rRNA genes and depletion at gene boundaries. These results illustrate the power of the microbiome for discovering prokaryotic epigenetics and the widespread distribution of oxidation-sensitive PTs in gut microbes.

One-sentence Summary

Application of informatic, mass spectrometric, and sequencing-based mapping tools to human gut bacteria revealed new phosphorothioate epigenetic systems widespread in the gut microbiome.

IGF2BP2 Maintains Retinal Pigment Epithelium Homeostasis by Stabilizing PAX6 and OTX2

Purpose

N6-methyladenosine (m6A) methylation is a chemical modification that occurs on RNA molecules, where the hydrogen atom of adenine (A) nucleotides is replaced by a methyl group, forming N6-methyladenosine. This modification is a dynamic and reversible process that plays a crucial role in regulating various biological processes, including RNA stability, transport, translation, and degradation. Currently, there is a lack of research on the role of m6A modifications in maintaining the characteristics of RPE cells. m6A readers play a crucial role in executing the functions of m6A modifications, which prompted our investigation into their regulatory roles in the RPE.

Methods

Phagocytosis assays, immunofluorescence staining, flow cytometry experiments, β-galactosidase staining, and RNA sequencing (RNA-seq) were conducted to assess the functional and cellular characteristics changes in retinal pigment epithelium (RPE) cells following short-hairpin RNA-mediated knockdown of insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2). RNA-seq and ultraviolet crosslinking immunoprecipitation with high-throughput sequencing (HITS-CLIP) were employed to identify the target genes regulated by IGF2BP2. adeno-associated virus (AAV) subretinal injection was performed in 6- to 8-week-old C57 mice to reduce IGF2BP2 expression in the RPE, and the impact of IGF2BP2 knockdown on mouse visual function was assessed using immunofluorescence, quantitative real-time PCR, optical coherence tomography, and electroretinography.

Results

IGF2BP2 was found to have a pronounced effect on RPE phagocytosis. Subsequent in-depth exploration revealed that IGF2BP2 modulates the mRNA stability of PAX6 and OTX2, and the loss of IGF2BP2 induces inflammatory and aging phenotypes in RPE cells. IGF2BP2 knockdown impaired RPE function, leading to retinal dysfunction in vivo.

Conclusions

Our data suggest a crucial role of IGF2BP2 as an m6A reader in maintaining RPE homeostasis by regulating the stability of PAX6 and OTX2, making it a potential target for preventing the occurrence of retinal diseases related to RPE malfunction.

Multiomic profiling of transcription factor binding and function in human brain

Transcription factors (TFs) orchestrate gene expression programs crucial for brain function, but we lack detailed information about TF binding in human brain tissue. We generated a multiomic resource (ChIP-seq, ATAC-seq, RNA-seq, DNA methylation) on bulk tissues and sorted nuclei from several postmortem brain regions, including binding maps for more than 100 TFs. We demonstrate improved measurements of TF activity, including motif recognition and gene expression modeling, upon identification and removal of high TF occupancy regions. Further, predictive TF binding models demonstrate a bias for these high-occupancy sites. Neuronal TFs SATB2 and TBR1 bind unique regions depleted for such sites and promote neuronal gene expression. Binding sites for TFs, including TBR1 and PKNOX1, are enriched for risk variants associated with neuropsychiatric disorders, predominantly in neurons. This work, titled BrainTF, is a powerful resource for future studies seeking to understand the roles of specific TFs in regulating gene expression in the human brain.

PT-seq: A method for metagenomic analysis of phosphorothioate epigenetics in complex microbial communities

Among dozens of known epigenetic marks, naturally occurring phosphorothioate (PT) DNA modifications are unique in replacing a non-bridging phosphate oxygen with redox-active sulfur and function in prokaryotic restriction-modification and transcriptional regulation. Interest in PTs has grown due to the widespread distribution of the dnd, ssp , and brx genes among bacteria and archaea, as well as the discovery of PTs in 5-10% of gut microbes. Efforts to map PTs in complex microbiomes using existing next-generation and direct sequencing technologies have failed due to poor sensitivity. Here we developed PT-seq as a high-sensitivity method to quantitatively map PTs across genomes and metagenomically identify PT-containing microbes in complex genomic mixtures. Like other methods for mapping PTs in individual genomes, PT-seq exploits targeted DNA strand cleavage at PTs by iodine, followed by sequencing library construction using ligation or template switching approaches. However, PT-specific sequencing reads are dramatically increased by adding steps to heat denature the DNA, block pre-existing 3'-ends, fragment DNA after T-tailing, and enrich iodine-induced breaks using biotin-labeling and streptavidin beads capture. Iterative optimization of the sensitivity and specificity of PT-seq is demonstrated with individual bacteria and human fecal DNA.

A repertoire of intronic lariat RNAs reveals tissue-specific regulation and target mimicry potential in plants

Lariat RNA is concomitantly produced by excised intron during RNA splicing, which is usually debranched by DBR1, an RNA debranching enzyme. However, increasing evidence showed that some lariat RNA could escape debranching. Little is known about how and why these lariat RNAs could be retained. By comparing the atlas of lariat RNAs between the non-dividing cell (mature pollen) and three actively dividing tissues (young shoot apex, young seeds, and young roots), we identified hundreds to thousands of lariat RNA naturally retained in each tissue, and the incidence of lariat RNA retention is much less in shoot apex while much more in pollen. Many lariat RNAs derived from the same intron or different lariat RNAs from the same pre-mRNA could be retained in one tissue while degraded in the other tissues. By deciphering lariat RNA sequences, we identified an AG-rich (RAAAAVAAAR) motif and a UC-rich (UCUCUYUCUC) motif for pollen-specific and the other three tissues-retained lariat RNAs, respectively. Reconstitution of the pollen-specific AG-rich motif indeed enhanced lariat RNA retention in plants. Biologically, hundreds of lariat RNAs harbored miRNA binding sites, and dual-luciferase reporter assay showed that these natural lariat RNAs had the potential to protect expression of miRNA target genes. Collectively, our results uncover that selective retention of lariat RNA is an actively regulatory process, and provide new insights into understanding how lariat RNA metabolism may impact miRNA activity.

De novo structure prediction of meteorin and meteorin-like protein for identification of domains, functional receptor binding regions, and their high-risk missense variants

Meteorin (Metrn) and Meteorin-like (Metrnl) are homologous secreted proteins involved in neural development and metabolic regulation. In this study, we have performed de novo structure prediction and analysis of both Metrn and Metrnl using Alphafold2 (AF2) and RoseTTAfold (RF). Based on the domain and structural homology analysis of the predicted structures, we have identified that these proteins are composed of two functional domains, a CUB domain and an NTR domain, connected by a hinge/loop region. We have identified the receptor binding regions of Metrn and Metrnl using the machine-learning tools ScanNet and Masif. These were further validated by docking Metrnl with its reported KIT receptor, thus establishing the role of each domain in the receptor interaction. Also, we have studied the effect of non-synonymous SNPs on the structure and function of these proteins using an array of bioinformatics tools and selected 16 missense variants in Metrn and 10 in Metrnl that can affect the protein stability. This is the first study to comprehensively characterize the functional domains of Metrn and Metrnl at their structural level and identify the functional domains, and protein binding regions. This study also highlights the interaction mechanism of the KIT receptor and Metrnl. The predicted deleterious SNPs will allow further understanding of the role of these variants in modulating the plasma levels of these proteins in disease conditions such as diabetes.Communicated by Ramaswamy H. Sarma.

Transcriptional control of the Cryptosporidium life cycle

The parasite Cryptosporidium is a leading agent of diarrhoeal disease in young children, and a cause and consequence of chronic malnutrition1,2. There are no vaccines and only limited treatment options3. The parasite infects enterocytes, in which it engages in asexual and sexual replication4, both of which are essential to continued infection and transmission. However, their molecular mechanisms remain largely unclear5. Here we use single-cell RNA sequencing to reveal the gene expression programme of the entire Cryptosporidium parvum life cycle in culture and in infected animals. Diverging from the prevailing model6, we find support for only three intracellular stages: asexual type-I meronts, male gamonts and female gametes. We reveal a highly organized program for the assembly of components at each stage. Dissecting the underlying regulatory network, we identify the transcription factor Myb-M as the earliest determinant of male fate, in an organism that lacks genetic sex determination. Conditional expression of this factor overrides the developmental program and induces widespread maleness, while conditional deletion ablates male development. Both have a profound impact on the infection. A large set of stage-specific genes now provides the opportunity to understand, engineer and disrupt parasite sex and life cycle progression to advance the development of vaccines and treatments.

Genome-wide identification and expression profiles of sex-related gene families in the Pacific abalone Haliotis discus hannai

In recent years, members of the Dmrt family, TGF-β superfamily and Sox family have been recognized as crucial genes for sex determination/differentiation across diverse animal species. Nevertheless, knowledge regarding the abundance and potential functions of these genes in abalone remains limited. In this study, a total of 5, 10, and 7 members of the Dmrt family, the TGF-β superfamily and the Sox family, respectively, were identified in the Pacific abalone Haliotis discus hannai. Sequence characteristics, phylogenetic relationships and spatiotemporal expression profiles of these genes were investigated. Notably, HdDmrt-04 (Dmrt1/1L-like) emerged as a potential mollusc-specific gene with a preponderance for expression in the testis. Interestingly, none of the TGF-β superfamily members exhibited specific or elevated expression in the gonads, highlighting the need for further investigation into their role in abalone sex differentiation. The Sox proteins in H. discus hannai were categorized into 7 subfamilies: B1, B2, C, D, E, F, and H. Among them, HdSox-07 (SoxH-like) was observed to play a crucial role in testis development, while HdSox-03 (SoxB1-like) and HdSox-04 (SoxC-like) probably cooperate in abalone ovary development. Taken together, the results of the present study suggested that HdDmrt-04 and HdSox-07 can be used as male-specific markers for gonad differentiation in H. discus hannai and imply conservation of their functions across invertebrates and vertebrates. Our findings provide new insights into the evolution and genetic structure of the Dmrt family, the TGF-β superfamily and the Sox family in abalone and pave the way for a deeper understanding of sex differentiation in gastropods.

Manganese-dependent transcription regulation by MntR and PerR in Thermus thermophilus HB8

Bacteria contain conserved mechanisms to control the intracellular levels of metal ions. Metalloregulatory transcription factors bind metal cations and play a central role in regulating gene expression of metal transporters. Often, these transcription factors regulate transcription by binding to a specific DNA sequence in the promoter region of target genes. Understanding the preferred DNA-binding sequence for transcriptional regulators can help uncover novel gene targets and provide insight into the biological role of the transcription factor in the host organism. Here, we identify consensus DNA-binding sequences and subsequent transcription regulatory networks for two metalloregulators from the ferric uptake regulator (FUR) and diphtheria toxin repressor (DtxR) superfamilies in Thermus thermophilus HB8. By homology search, we classify the DtxR homolog as a manganese-specific, MntR (TtMntR), and the FUR homolog as a peroxide-sensing, PerR (TtPerR). Both transcription factors repress separate ZIP transporter genes in vivo, and TtPerR acts as a bifunctional transcription regulator by activating the expression of ferric and hemin transport systems. We show TtPerR and TtMntR bind DNA in the presence of manganese in vitro and in vivo; however, TtPerR is unable to bind DNA in the presence of iron, likely due to iron-mediated histidine oxidation. Unlike canonical PerR homologs, TtPerR does not appear to contribute to peroxide detoxification. Instead, the TtPerR regulon and DNA binding sequence are more reminiscent of Fur or Mur homologs. Collectively, these results highlight the similarities and differences between two metalloregulatory superfamilies and underscore the interplay of manganese and iron in transcription factor regulation.

Structure-based learning to predict and model protein-DNA interactions and transcription-factor co-operativity in cis-regulatory elements

Transcription factor (TF) binding is a key component of genomic regulation. There are numerous high-throughput experimental methods to characterize TF-DNA binding specificities. Their application, however, is both laborious and expensive, which makes profiling all TFs challenging. For instance, the binding preferences of ∼25% human TFs remain unknown; they neither have been determined experimentally nor inferred computationally. We introduce a structure-based learning approach to predict the binding preferences of TFs and the automated modelling of TF regulatory complexes. We show the advantage of using our approach over the classical nearest-neighbor prediction in the limits of remote homology. Starting from a TF sequence or structure, we predict binding preferences in the form of motifs that are then used to scan a DNA sequence for occurrences. The best matches are either profiled with a binding score or collected for their subsequent modeling into a higher-order regulatory complex with DNA. Co-operativity is modelled by: (i) the co-localization of TFs and (ii) the structural modeling of protein-protein interactions between TFs and with co-factors. We have applied our approach to automatically model the interferon-β enhanceosome and the pioneering complexes of OCT4, SOX2 (or SOX11) and KLF4 with a nucleosome, which are compared with the experimentally known structures.

A comprehensive analysis of pneumococcal two-component system regulatory networks

Two-component systems are key signal-transduction systems that enable bacteria to respond to a wide variety of environmental stimuli. The human pathogen, Streptococcus pneumoniae (pneumococcus) encodes 13 two-component systems and a single orphan response regulator, most of which are significant for pneumococcal pathogenicity. Mapping the regulatory networks governed by these systems is key to understand pneumococcal host adaptation. Here we employ a novel bioinformatic approach to predict the regulons of each two-component system based on publicly available whole-genome sequencing data. By employing pangenome-wide association studies (panGWAS) to predict genotype-genotype associations for each two-component system, we predicted regulon genes of 11 of the pneumococcal two-component systems. Through validation via next-generation RNA-sequencing on response regulator overexpression mutants, several top candidate genes predicted by the panGWAS analysis were confirmed as regulon genes. The present study presents novel details on multiple pneumococcal two-component systems, including an expansion of regulons, identification of candidate response regulator binding motifs, and identification of candidate response regulator-regulated small non-coding RNAs. We also demonstrate a use for panGWAS as a complementary tool in target gene identification via identification of genotype-to-genotype links. Expanding our knowledge on two-component systems in pathogens is crucial to understanding how these bacteria sense and respond to their host environment, which could prove useful in future drug development.

Rice small secreted peptide, OsRALF26, recognized by FERONIA-like receptor 1 induces immunity in rice and Arabidopsis

Rapid alkalinization factors (RALFs), belonging to a family of small secreted peptides, have been considered as important signaling molecules in diverse biological processes, including immunity. Current studies on RALF-modulated immunity mainly focus on Arabidopsis, but little is reported in crop plants. The rice immune receptor XA21 confers immunity to the bacterial blight pathogen, Xanthomonas oryzae pv. oryzae (Xoo). Here, we pursued functional characterization of rice RALF26 (OsRALF26) up-regulated by Xoo during XA21-mediated immune response. When applied exogenously as a recombinant peptide, OsRALF26 induced a series of immune responses, including pathogenesis-related genes (PRs) induction, reactive oxygen species (ROS) production, and callose deposition in rice and/or Arabidopsis. Transgenic rice and Arabidopsis overexpressing OsRALF26 exhibited significantly enhanced resistance to Xoo and Pseudomonas syringae pv. tomato DC3000 (Pst DC3000), respectively. In yeast two-hybrid, pull-down assays, and co-immunoprecipitation analyses, rice FER-like receptor 1 (OsFLR1) was identified as a receptor of OsRALF26. Transient expression of OsFLR1 in Nicotiana benthamiana leaves displayed significantly increased ROS production and callose deposition after OsRALF26 treatment. Together, we propose that OsRALF26 induced by Xoo in an XA21-dependent manner is perceived by OsFLR1 and may play a novel role in the enforcement of XA21-mediated immunity.

Comprehensive structural overview of the C-terminal ligand-binding domains of the TetR family regulators

TetR family regulators (TFRs) represent a large group of one-component bacterial signal transduction systems which recognize environmental signals, like the presence of antibiotics or other bactericidal compounds, and trigger the cell response by regulating the expression of genes that secure bacterial survival in harsh environmental conditions. TFRs act as homodimers, each protomer is composed of a conserved DNA-binding N-terminal domain (NTD) and a variable ligand-binding C-terminal domain (CTD). Currently, there are about 500 structures of TFRs available in the Protein Data Bank and one-fourth of them represent the structures of TFR-ligand complexes. In this review, we summarized information on the ligands interacting with TFRs and based on structural data, we compared the CTDs of the TFR family members, as well as their ligand-binding cavities. Additionally, we divided the whole TFR family, including more than half of a million sequences, into subfamilies according to calculated multiple sequence alignment and phylogenetic tree. We also highlighted structural elements characteristic of some of the subfamilies. The presented comprehensive overview of the TFR CTDs provides good bases and future directions for further studies on TFRs that are not only important targets for battling multidrug resistance but also good candidates for many biotechnological approaches, like TFR-based biosensors.