Assigning roles to DNA regulatory motifs using comparative genomics

Deciphering the mechanisms involved in reduced sensitivity to azoles and fengycin lipopeptide in Venturia inaequalis

Apple scab, caused by the hemibiotrophic fungus Venturia inaequalis, is currently the most common and damaging disease in apple orchards. Two strains of V. inaequalis (S755 and Rs552) with different sensitivities to azole fungicides and the bacterial metabolite fengycin were compared to determine the mechanisms responsible for these differences. Antifungal activity tests showed that Rs552 had reduced sensitivity to tebuconazole and tetraconazole, as well as to fengycin alone or in a binary mixture with other lipopeptides (iturin A, pumilacidin, lichenysin). S755 was highly sensitive to fengycin, whose activity was close to that of tebuconazole. Unlike fengycin, lipopeptides from the iturin family (mycosubtilin, iturin A) had similar activity on both strains, while those from the surfactin family (lichenysin, pumilacidin) were not active, except in binary mixtures with fengycin. The activity of lipopeptides varies according to their family and structure. Analyses to determine the difference in sensitivity to azoles (which target the CYP51 enzyme involved in the ergosterol biosynthesis pathway) showed that the reduced sensitivity in Rs552 is linked to (i) a constitutive increased expression of the Cyp51A gene caused by insertions in the upstream region and (ii) greater efflux by membrane pumps with the involvement of ABC transporters. Microscopic observations revealed that fengycin, known to interact with plasma membranes, induced morphological and cytological changes in cells from both strains. Sterol and phospholipid analyses showed a higher level of ergosta-7,22-dien-3-ol and a lower level of PI(C16:0/C18:1) in Rs552 compared with S755. These differences could therefore influence the composition of the plasma membrane and explain the differential sensitivity of the strains to fengycin. However, the similar antifungal activities of mycosubtilin and iturin A in the two strains indirectly indicate that sterols are probably not involved in the fengycin resistance mechanism. This leads to the conclusion that different mechanisms are responsible for the difference in susceptibility to azoles or fengycin in the strains studied.

Identification of DNA motif pairs on paired sequences based on composite heterogeneous graph

Motivation

The interaction between DNA motifs (DNA motif pairs) influences gene expression through partnership or competition in the process of gene regulation. Potential chromatin interactions between different DNA motifs have been implicated in various diseases. However, current methods for identifying DNA motif pairs rely on the recognition of single DNA motifs or probabilities, which may result in local optimal solutions and can be sensitive to the choice of initial values. A method for precisely identifying DNA motif pairs is still lacking.

Results

Here, we propose a novel computational method for predicting DNA Motif Pairs based on Composite Heterogeneous Graph (MPCHG). This approach leverages a composite heterogeneous graph model to identify DNA motif pairs on paired sequences. Compared with the existing methods, MPCHG has greatly improved the accuracy of motifs prediction. Furthermore, the predicted DNA motifs demonstrate heightened DNase accessibility than the background sequences. Notably, the two DNA motifs forming a pair exhibit functional consistency. Importantly, the interacting TF pairs obtained by predicted DNA motif pairs were significantly enriched with known interacting TF pairs, suggesting their potential contribution to chromatin interactions. Collectively, we believe that these identified DNA motif pairs held substantial implications for revealing gene transcriptional regulation under long-range chromatin interactions.

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.

Interpretable prediction of mRNA abundance from promoter sequence using contextual regression models

While machine learning models have been successfully applied to predicting gene expression from promoter sequences, it remains a great challenge to derive intuitive interpretation of the model and reveal DNA motif grammar such as motif cooperation and distance constraint between motif sites. Previous interpretation approaches are often time-consuming or have difficulty to learn the combinatory rules. In this work, we designed interpretable neural network models to predict the mRNA expression levels from DNA sequences. By applying the Contextual Regression framework we developed, we extracted weighted features to cluster samples into different groups, which have different gene expression levels. We performed motif analysis in each cluster and found motifs with active or repressive regulation on gene expression. By comparing the co-occurrence locations of discovered motifs, we also uncovered multiple grammars of motif combination including communities of cooperative motifs and distance constraints between motif pairs. These results revealed new insights of the regulatory architecture of promoter sequences.

A single dose of cocaine rewires the 3D genome structure of midbrain dopamine neurons

Midbrain dopamine neurons (DNs) respond to a first exposure to addictive drugs and play key roles in chronic drug usage1-3. As the synaptic and transcriptional changes that follow an acute cocaine exposure are mostly resolved within a few days4,5, the molecular changes that encode the long-term cellular memory of the exposure within DNs remain unknown. To investigate whether a single cocaine exposure induces long-term changes in the 3D genome structure of DNs, we applied Genome Architecture Mapping and single nucleus transcriptomic analyses in the mouse midbrain. We found extensive rewiring of 3D genome architecture at 24 hours past exposure which remains or worsens by 14 days, outlasting transcriptional responses. The cocaine-induced chromatin rewiring occurs at all genomic scales and affects genes with major roles in cocaine-induced synaptic changes. A single cocaine exposure triggers extensive long-lasting changes in chromatin condensation in post-synaptic and post-transcriptional regulatory genes, for example the unfolding of Rbfox1 which becomes most prominent 14 days post exposure. Finally, structurally remodeled genes are most expressed in a specific DN sub-type characterized by low expression of the dopamine auto-receptor Drd2, a key feature of highly cocaine-sensitive cells. These results reveal an important role for long-lasting 3D genome remodelling in the cellular memory of a single cocaine exposure, providing new hypotheses for understanding the inception of drug addiction and 3D genome plasticity.

MMGAT: a graph attention network framework for ATAC-seq motifs finding

BACKGROUND

Motif finding in Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data is essential to reveal the intricacies of transcription factor binding sites (TFBSs) and their pivotal roles in gene regulation. Deep learning technologies including convolutional neural networks (CNNs) and graph neural networks (GNNs), have achieved success in finding ATAC-seq motifs. However, CNN-based methods are limited by the fixed width of the convolutional kernel, which makes it difficult to find multiple transcription factor binding sites with different lengths. GNN-based methods has the limitation of using the edge weight information directly, makes it difficult to aggregate the neighboring nodes' information more efficiently when representing node embedding.

RESULTS

To address this challenge, we developed a novel graph attention network framework named MMGAT, which employs an attention mechanism to adjust the attention coefficients among different nodes. And then MMGAT finds multiple ATAC-seq motifs based on the attention coefficients of sequence nodes and k-mer nodes as well as the coexisting probability of k-mers. Our approach achieved better performance on the human ATAC-seq datasets compared to existing tools, as evidenced the highest scores on the precision, recall, F1_score, ACC, AUC, and PRC metrics, as well as finding 389 higher quality motifs. To validate the performance of MMGAT in predicting TFBSs and finding motifs on more datasets, we enlarged the number of the human ATAC-seq datasets to 180 and newly integrated 80 mouse ATAC-seq datasets for multi-species experimental validation. Specifically on the mouse ATAC-seq dataset, MMGAT also achieved the highest scores on six metrics and found 356 higher-quality motifs. To facilitate researchers in utilizing MMGAT, we have also developed a user-friendly web server named MMGAT-S that hosts the MMGAT method and ATAC-seq motif finding results.

CONCLUSIONS

The advanced methodology MMGAT provides a robust tool for finding ATAC-seq motifs, and the comprehensive server MMGAT-S makes a significant contribution to genomics research. The open-source code of MMGAT can be found at https://github.com/xiaotianr/MMGAT , and MMGAT-S is freely available at https://www.mmgraphws.com/MMGAT-S/ .

Employing bimodal representations to predict DNA bendability within a self-supervised pre-trained framework

The bendability of genomic DNA, which measures the DNA looping rate, is crucial for numerous biological processes of DNA. Recently, an advanced high-throughput technique known as 'loop-seq' has made it possible to measure the inherent cyclizability of DNA fragments. However, quantifying the bendability of large-scale DNA is costly, laborious, and time-consuming. To close the gap between rapidly evolving large language models and expanding genomic sequence information, and to elucidate the DNA bendability's impact on critical regulatory sequence motifs such as super-enhancers in the human genome, we introduce an innovative computational model, named MIXBend, to forecast the DNA bendability utilizing both nucleotide sequences and physicochemical properties. In MIXBend, a pre-trained language model DNABERT and convolutional neural network with attention mechanism are utilized to construct both sequence- and physicochemical-based extractors for the sophisticated refinement of DNA sequence representations. These bimodal DNA representations are then fed to a k-mer sequence-physicochemistry matching module to minimize the semantic gap between each modality. Lastly, a self-attention fusion layer is employed for the prediction of DNA bendability. In conclusion, the experimental results validate MIXBend's superior performance relative to other state-of-the-art methods. Additionally, MIXBend reveals both novel and known motifs from the yeast. Moreover, MIXBend discovers significant bendability fluctuations within super-enhancer regions and transcription factors binding sites in the human genome.

Synthetic reversed sequences reveal default genomic states

Pervasive transcriptional activity is observed across diverse species. The genomes of extant organisms have undergone billions of years of evolution, making it unclear whether these genomic activities represent effects of selection or 'noise'1-4. Characterizing default genome states could help understand whether pervasive transcriptional activity has biological meaning. Here we addressed this question by introducing a synthetic 101-kb locus into the genomes of Saccharomyces cerevisiae and Mus musculus and characterizing genomic activity. The locus was designed by reversing but not complementing human HPRT1, including its flanking regions, thus retaining basic features of the natural sequence but ablating evolved coding or regulatory information. We observed widespread activity of both reversed and native HPRT1 loci in yeast, despite the lack of evolved yeast promoters. By contrast, the reversed locus displayed no activity at all in mouse embryonic stem cells, and instead exhibited repressive chromatin signatures. The repressive signature was alleviated in a locus variant lacking CpG dinucleotides; nevertheless, this variant was also transcriptionally inactive. These results show that synthetic genomic sequences that lack coding information are active in yeast, but inactive in mouse embryonic stem cells, consistent with a major difference in 'default genomic states' between these two divergent eukaryotic cell types, with implications for understanding pervasive transcription, horizontal transfer of genetic information and the birth of new genes.

Comprehensive transcriptomic meta-analysis unveils new responsive genes to methyl jasmonate and ethylene in Catharanthusroseus

In Catharanthus roseus, vital plant hormones, namely methyl jasmonate (MeJA) and ethylene, serve as abiotic triggers, playing a crucial role in stimulating the production of specific secondary compounds with anticancer properties. Understanding how plants react to various stresses, stimuli, and the pathways involved in biosynthesis holds significant promise. The application of stressors like ethylene and MeJA induces the plant's defense mechanisms, leading to increased secondary metabolite production. To delve into the essential transcriptomic processes linked to hormonal responses, this study employed an integrated approach combining RNA-Seq data meta-analysis and system biology methodologies. Furthermore, the validity of the meta-analysis findings was confirmed using RT-qPCR. Within the meta-analysis, 903 genes exhibited differential expression (DEGs) when comparing normal conditions to those of the treatment. Subsequent analysis, encompassing gene ontology, KEGG, TF, and motifs, revealed that these DEGs were actively engaged in multiple biological processes, particularly in responding to various stresses and stimuli. Additionally, these genes were notably enriched in diverse biosynthetic pathways, including those related to TIAs, housing valuable medicinal compounds found in this plant. Furthermore, by conducting co-expression network analysis, we identified hub genes within modules associated with stress response and the production of TIAs. Most genes linked to the biosynthesis pathway of TIAs clustered within three specific modules. Noteworthy hub genes, including Helicase ATP-binding domain, hbdA, and ALP1 genes within the blue, turquoise, and green module networks, are presumed to play a role in the TIAs pathway. These identified candidate genes hold potential for forthcoming genetic and metabolic engineering initiatives aimed at augmenting the production of secondary metabolites and medicinal compounds within C. roseus.

Cacao pod transcriptome profiling of seven genotypes identifies features associated with post-penetration resistance to Phytophthora palmivora

The oomycete Phytophthora palmivora infects the fruit of cacao trees (Theobroma cacao) causing black pod rot and reducing yields. Cacao genotypes vary in their resistance levels to P. palmivora, yet our understanding of how cacao fruit respond to the pathogen at the molecular level during disease establishment is limited. To address this issue, disease development and RNA-Seq studies were conducted on pods of seven cacao genotypes (ICS1, WFT, Gu133, Spa9, CCN51, Sca6 and Pound7) to better understand their reactions to the post-penetration stage of P. palmivora infection. The pod tissue-P. palmivora pathogen assay resulted in the genotypes being classified as susceptible (ICS1, WFT, Gu133 and Spa9) or resistant (CCN51, Sca6 and Pound7). The number of differentially expressed genes (DEGs) ranged from 1625 to 6957 depending on genotype. A custom gene correlation approach identified 34 correlation groups. De novo motif analysis was conducted on upstream promoter sequences of differentially expressed genes, identifying 76 novel motifs, 31 of which were over-represented in the upstream sequences of correlation groups and associated with gene ontology terms related to oxidative stress response, defense against fungal pathogens, general metabolism and cell function. Genes in one correlation group (Group 6) were strongly induced in all genotypes and enriched in genes annotated with defense-responsive terms. Expression pattern profiling revealed that genes in Group 6 were induced to higher levels in the resistant genotypes. An additional analysis allowed the identification of 17 candidate cis-regulatory modules likely to be involved in cacao defense against P. palmivora. This study is a comprehensive exploration of the cacao pod transcriptional response to P. palmivora spread after infection. We identified cacao genes, promoter motifs, and promoter motif combinations associated with post-penetration resistance to P. palmivora in cacao pods and provide this information as a resource to support future and ongoing efforts to breed P. palmivora-resistant cacao.

Transcriptome analysis in Aegilops tauschii unravels further insights into genetic control of stripe rust resistance

MAIN CONCLUSION

The Aegilops tauschii resistant accession prevented the pathogen colonization by controlling the sugar flow and triggering the hypersensitive reaction. This study suggested that NBS-LRRs probably induce resistance through bHLH by controlling JA- and SA-dependent pathways. Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst) is one of wheat's most destructive fungal diseases that causes a severe yield reduction worldwide. The most effective and economically-friendly strategy to manage this disease is genetic resistance which can be achieved through deploying new and effective resistance genes. Aegilops tauschii, due to its small genome and co-evolution with Pst, can provide detailed information about underlying resistance mechanisms. Hence, we used RNA-sequencing approach to identify the transcriptome variations of two contrasting resistant and susceptible Ae. tauschii accessions in interaction with Pst and differentially expressed genes (DEGs) for resistance to stripe rust. Gene ontology, pathway analysis, and search for functional domains, transcription regulators, resistance genes, and protein-protein interactions were used to interpret the results. The genes encoding NBS-LRR, CC-NBS-kinase, and phenylalanine ammonia-lyase, basic helix-loop-helix (bHLH)-, basic-leucine zipper (bZIP)-, APETALA2 (AP2)-, auxin response factor (ARF)-, GATA-, and LSD-like transcription factors were up-regulated exclusively in the resistant accession. The key genes involved in response to salicylic acid, amino sugar and nucleotide sugar metabolism, and hypersensitive response contributed to plant defense against stripe rust. The activation of jasmonic acid biosynthesis and starch and sucrose metabolism pathways under Pst infection in the susceptible accession explained the colonization of the host. Overall, this study can fill the gaps in the literature on host-pathogen interaction and enrich the Ae. tauschii transcriptome sequence information. It also suggests candidate genes that could guide future breeding programs attempting to develop rust-resistant cultivars.

A comprehensive meta-analysis of transcriptome data to identify signature genes associated with pancreatic ductal adenocarcinoma

PURPOSE

Pancreatic ductal adenocarcinoma (PDAC) has a five-year survival rate of less than 5%. Absence of symptoms at primary tumor stages, as well as high aggressiveness of the tumor can lead to high mortality in cancer patients. Most patients are recognized at the advanced or metastatic stage without surgical symptom, because of the lack of reliable early diagnostic biomarkers. The objective of this work was to identify potential cancer biomarkers by integrating transcriptome data.

METHODS

Several transcriptomic datasets comprising of 11 microarrays were retrieved from the GEO database. After pre-processing, a meta-analysis was applied to identify differentially expressed genes (DEGs) between tumor and nontumor samples for datasets. Next, co-expression analysis, functional enrichment and survival analyses were used to determine the functional properties of DEGs and identify potential prognostic biomarkers. In addition, some regulatory factors involved in PDAC including transcription factors (TFs), protein kinases (PKs), and miRNAs were identified.

RESULTS

After applying meta-analysis, 1074 DEGs including 539 down- and 535 up-regulated genes were identified. Pathway enrichment analyzes using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed that DEGs were significantly enriched in the HIF-1 signaling pathway and focal adhesion. The results also showed that some of the DEGs were assigned to TFs that belonged to 23 conserved families. Sixty-four PKs were identified among the DEGs that showed the CAMK family was the most abundant group. Moreover, investigation of corresponding upstream regions of DEGs identified 11 conserved sequence motifs. Furthermore, weighted gene co-expression network analysis (WGCNA) identified 8 modules, more of them were significantly enriched in Ras signaling, p53 signaling, MAPK signaling pathways. In addition, several hubs in modules were identified, including EMP1, EVL, ELP5, DEF8, MTERF4, GLUP1, CAPN1, IGF1R, HSD17B14, TOM1L2 and RAB11FIP3. According to survival analysis, it was identified that the expression levels of two genes, EMP1 and RAB11FIP3 are related to prognosis.

CONCLUSION

We identified several genes critical for PDAC based on meta-analysis and system biology approach. These genes may serve as potential targets for the treatment and prognosis of PDAC.

Modeling methyl-sensitive transcription factor motifs with an expanded epigenetic alphabet

BACKGROUND

Transcription factors bind DNA in specific sequence contexts. In addition to distinguishing one nucleobase from another, some transcription factors can distinguish between unmodified and modified bases. Current models of transcription factor binding tend not to take DNA modifications into account, while the recent few that do often have limitations. This makes a comprehensive and accurate profiling of transcription factor affinities difficult.

RESULTS

Here, we develop methods to identify transcription factor binding sites in modified DNA. Our models expand the standard A/C/G/T DNA alphabet to include cytosine modifications. We develop Cytomod to create modified genomic sequences and we also enhance the MEME Suite, adding the capacity to handle custom alphabets. We adapt the well-established position weight matrix (PWM) model of transcription factor binding affinity to this expanded DNA alphabet. Using these methods, we identify modification-sensitive transcription factor binding motifs. We confirm established binding preferences, such as the preference of ZFP57 and C/EBPβ for methylated motifs and the preference of c-Myc for unmethylated E-box motifs.

CONCLUSIONS

Using known binding preferences to tune model parameters, we discover novel modified motifs for a wide array of transcription factors. Finally, we validate our binding preference predictions for OCT4 using cleavage under targets and release using nuclease (CUT&RUN) experiments across conventional, methylation-, and hydroxymethylation-enriched sequences. Our approach readily extends to other DNA modifications. As more genome-wide single-base resolution modification data becomes available, we expect that our method will yield insights into altered transcription factor binding affinities across many different modifications.

PAX3-FOXO1 dictates myogenic reprogramming and rhabdomyosarcoma identity in endothelial progenitors

Fusion-positive rhabdomyosarcoma (FP-RMS) driven by the expression of the PAX3-FOXO1 (P3F) fusion oncoprotein is an aggressive subtype of pediatric rhabdomyosarcoma. FP-RMS histologically resembles developing muscle yet occurs throughout the body in areas devoid of skeletal muscle highlighting that FP-RMS is not derived from an exclusively myogenic cell of origin. Here we demonstrate that P3F reprograms mouse and human endothelial progenitors to FP-RMS. We show that P3F expression in aP2-Cre expressing cells reprograms endothelial progenitors to functional myogenic stem cells capable of regenerating injured muscle fibers. Further, we describe a FP-RMS mouse model driven by P3F expression and Cdkn2a loss in endothelial cells. Additionally, we show that P3F expression in TP53-null human iPSCs blocks endothelial-directed differentiation and guides cells to become myogenic cells that form FP-RMS tumors in immunocompromised mice. Together these findings demonstrate that FP-RMS can originate from aberrant development of non-myogenic cells driven by P3F.

Analysis of the chromatin landscape and RNA polymerase II binding at SIN3-regulated genes

The chromatin environment has a significant impact on gene expression. Chromatin structure is highly regulated by histone modifications and RNA polymerase II binding dynamics. The SIN3 histone modifying complex regulates the chromatin environment leading to changes in gene expression. In Drosophila melanogaster, the Sin3A gene is alternatively spliced to produce different protein isoforms, two of which include SIN3 220 and SIN3 187. Both SIN3 isoforms are scaffolding proteins that interact with several other factors to regulate the chromatin landscape. The mechanism through which the SIN3 isoforms regulate chromatin is not well understood. Here, we analyze publicly available data sets to allow us to ask specific questions on how SIN3 isoforms regulate chromatin and gene activity. We determined that genes repressed by the SIN3 isoforms exhibited enrichment in histone H3K4me2, H3K4me3, H3K14ac and H3K27ac near the transcription start site. We observed an increase in the amount of paused RNA polymerase II on the promoter of genes repressed by the isoforms as compared to genes that require SIN3 for maximum activation. Furthermore, we analyzed a subset of genes regulated by SIN3 187 that suggest a mechanism in which SIN3 187 might exhibit hard regulation as well as soft regulation. Data presented here expand our knowledge of how the SIN3 isoforms regulate the chromatin environment and RNA polymerase II binding dynamics.

Karbanowicz T. Gene Networks Analysis of Salmonella Typhimurium Reveals New Insights on Key Genes Involved in Response to Low Water Activity

Background

When Salmonella enterica serovar Typhimurium, a foodborne bacterium, is exposed to osmotic stress, cellular adaptations increase virulence severity and cellular survival.

Objectives

The aim of the gene network analysis of S. Typhimurium was to provide insights into the various interactions between the genes involved in cellular survival under low water activity (aw).

Materials and Methods

We performed a gene network analysis to identify the gene clusters and hub genes of S. Typhimurium using Cytoscape in three food samples subjected to aw stress after 72 hours.

Results

The identified hub genes of S. Typhimurium belonged to down-regulated genes and were related to translation, transcription, and ribosome structure in the food samples. The rpsB and Tig were identified as the most important of the hub genes. Enrichment analysis of the hub genes also revealed the importance of translation and cellular protein metabolic processes. Moreover, the biological process associated with organonitrogen metabolism in milk chocolate was identified. According to the KEGG pathway results of gene cluster analysis, cellular responses to stress were associated with RNA polymerase, ribosome, and oxidative phosphorylation. Genes encoding RNA polymerase activity, including rpoA, rpoB, and rpoZ, were also significantly identified in the KEGG pathways. The identified motifs of hub DEGs included EXPREG_00000850, EXPREG_00000b00, EXPREG_000008e0, and EXPREG_00000850.

Conclusion

Based on the results of the gene network analysis, the identified hub genes may contribute to adaptation to food compositions and be responsible for the development of low water stress tolerance in Salmonella. Among the food samples, the milk chocolate matrix leads to more adaptation pathways for S. Typhimurium survival, as more hub genes were down-regulated and more motifs were detected. The identified motifs were involved in carbohydrate metabolism, carbohydrate transport, electron transfer, and oxygen transfer.

Systematic investigation of the homology sequences around the human fusion gene breakpoints in pan-cancer - bioinformatics study for a potential link to MMEJ

Microhomology-mediated end joining (MMEJ), an error-prone DNA damage repair mechanism, frequently leads to chromosomal rearrangements due to its ability to engage in promiscuous end joining of genomic instability and also leads to increasing mutational load at the sequences flanking the breakpoints (BPs). In this study, we systematically investigated the homology sequences around the genomic breakpoint area of human fusion genes, which were formed by the chromosomal rearrangements initiated by DNA double-strand breakage. Since the RNA-seq data is the typical data set to check the fusion genes, for the known exon junction fusion breakpoints identified from RNA-seq data, we have to infer the high chance of genomic breakpoint regions. For this, we utilized the high feature importance score area calculated from our recently developed fusion BP prediction model, FusionAI and identified 151 K microhomologies among ~24 K fusion BPs in 20 K fusion genes. From our multiple bioinformatics studies, we found a relationship between sequence homologies and the immune system. This in-silico study will provide novel knowledge on the sequence homologies around the coded structural variants.

The effects of maternal care on the developmental transcriptome and metatranscriptome of a wild bee

Maternal care acts as a strong environmental stimulus that can induce phenotypic plasticity in animals and may also alter their microbial communities through development. Here, we characterize the developmental metatranscriptome of the small carpenter bee, Ceratina calcarata, across developmental stages and in the presence or absence of mothers. Maternal care had the most influence during early development, with the greatest number and magnitude of differentially expressed genes between maternal care treatments, and enrichment for transcription factors regulating immune response in motherless early larvae. Metatranscriptomic data revealed fungi to be the most abundant group in the microbiome, with Aspergillus the most abundant in early larvae raised without mothers. Finally, integrative analysis between host transcriptome and metatranscriptome highlights several fungi correlating with developmental and immunity genes. Our results provide characterizations of the influence of maternal care on gene expression and the microbiome through development in a wild bee.

Predicting the impact of sequence motifs on gene regulation using single-cell data

The binding of transcription factors at proximal promoters and distal enhancers is central to gene regulation. Identifying regulatory motifs and quantifying their impact on expression remains challenging. Using a convolutional neural network trained on single-cell data, we infer putative regulatory motifs and cell type-specific importance. Our model, scover, explains 29% of the variance in gene expression in multiple mouse tissues. Applying scover to distal enhancers identified using scATAC-seq from the developing human brain, we identify cell type-specific motif activities in distal enhancers. Scover can identify regulatory motifs and their importance from single-cell data where all parameters and outputs are easily interpretable.

A systems biology approach to pathogenesis of gastric cancer: gene network modeling and pathway analysis

BACKGROUND

Gastric cancer (GC) ranks among the most common malignancies worldwide. This study aimed to find critical genes/pathways in GC pathogenesis.

METHODS

Gene interactions were analyzed, and the protein-protein interaction network was drawn. Then enrichment analysis of the hub genes was performed and network cluster analysis and promoter analysis of the hub genes were done. Age/sex analysis was done on the identified genes.

RESULTS

Eleven hub genes in GC were identified in the current study (ATP5A1, ATP5B, ATP5D, MT-ATP8, COX7A2, COX6C, ND4, ND6, NDUFS3, RPL8, and RPS16), mostly involved in mitochondrial functions. There was no report on the ATP5D, ND6, NDUFS3, RPL8, and RPS16 in GC. Our results showed that the most affected processes in GC are the metabolic processes, and the oxidative phosphorylation pathway was considerably enriched which showed the significance of mitochondria in GC pathogenesis. Most of the affected pathways in GC were also involved in neurodegenerative diseases. Promoter analysis showed that negative regulation of signal transduction might play an important role in GC pathogenesis. In the analysis of the basal expression pattern of the selected genes whose basal expression presented a change during the age, we found that a change in age may be an indicator of changes in disease insurgence and/or progression at different ages.

CONCLUSIONS

These results might open up new insights into GC pathogenesis. The identified genes might be novel diagnostic/prognostic biomarkers or potential therapeutic targets for GC. This work, being based on bioinformatics analysis act as a hypothesis generator that requires further clinical validation.

Profiling of Chromatin Accessibility in Pigs across Multiple Tissues and Developmental Stages

The study of chromatin accessibility across tissues and developmental stages is essential for elucidating the transcriptional regulation of various phenotypes and biological processes. However, the chromatin accessibility profiles of multiple tissues in newborn pigs and across porcine liver development remain poorly investigated. Here, we used ATAC-seq and rRNA-depleted RNA-seq to profile open chromatin maps and transcriptional features of heart, kidney, liver, lung, skeletal muscle, and spleen in newborn pigs and porcine liver tissue in the suckling and adult stages, respectively. Specifically, by analyzing a union set of protein-coding genes (PCGs) and two types of transcripts (lncRNAs and TUCPs), we obtained a comprehensive annotation of consensus ATAC-seq peaks for each tissue and developmental stage. As expected, the PCGs with tissue-specific accessible promoters had active transcription and were relevant to tissue-specific functions. In addition, other non-coding tissue-specific peaks were involved in both physical activity and the morphogenesis of neonatal tissues. We also characterized stage-specific peaks and observed a close association between dynamic chromatin accessibility and hepatic function transition during liver postnatal development. Overall, this study expands our current understanding of epigenetic regulation in mammalian tissues and organ development, which can benefit both economic trait improvement and improve the biomedical usage of pigs.

FOXA1 Reprogramming Dictates Retinoid X Receptor Response in ESR1-Mutant Breast Cancer

Estrogen receptor alpha (ER/ESR1) mutations occur in 30% to 40% of endocrine resistant ER-positive (ER+) breast cancer. Forkhead box A1 (FOXA1) is a key pioneer factor mediating ER-chromatin interactions and endocrine response in ER+ breast cancer, but its role in ESR1-mutant breast cancer remains unclear. Our previous FOXA1 chromatin immunoprecipitation sequencing (ChIP-seq) identified a large portion of redistributed binding sites in T47D genome-edited Y537S and D538G ESR1-mutant cells. Here, we further integrated FOXA1 genomic binding profile with the isogenic ER cistrome, accessible genome, and transcriptome data of T47D cell model. FOXA1 redistribution was significantly associated with transcriptomic alterations caused by ESR1 mutations. Furthermore, in ESR1-mutant cells, FOXA1-binding sites less frequently overlapped with ER, and differential gene expression was less associated with the canonical FOXA1-ER axis. Motif analysis revealed a unique enrichment of retinoid X receptor (RXR) motifs in FOXA1-binding sites of ESR1-mutant cells. Consistently, ESR1-mutant cells were more sensitive to growth stimulation with the RXR agonist LG268. The mutant-specific response was dependent on two RXR isoforms, RXR-α and RXR-β, with a stronger dependency on the latter. In addition, T3, the agonist of thyroid receptor (TR) also showed a similar growth-promoting effect in ESR1-mutant cells. Importantly, RXR antagonist HX531 blocked growth of ESR1-mutant cells and a patient-derived xenograft (PDX)-derived organoid with an ESR1 D538G mutation. Collectively, our data support the evidence for a stronger RXR response associated with FOXA1 reprograming in ESR1-mutant cells, suggesting development of therapeutic strategies targeting RXR pathways in breast tumors with ESR1 mutation.

IMPLICATIONS

It provides comprehensive characterization of the role of FOXA1 in ESR1-mutant breast cancer and potential therapeutic strategy through blocking RXR activation.

ZBTB7A regulates MDD-specific chromatin signatures and astrocyte-mediated stress vulnerability in orbitofrontal cortex

Hyperexcitability in the orbitofrontal cortex (OFC) is a key clinical feature of anhedonic domains of Major Depressive Disorder (MDD). However, the cellular and molecular substrates underlying this dysfunction remain unknown. Here, cell-population-specific chromatin accessibility profiling in human OFC unexpectedly mapped genetic risk for MDD exclusively to non-neuronal cells, and transcriptomic analyses revealed significant glial dysregulation in this region. Characterization of MDD-specific cis-regulatory elements identified ZBTB7A - a transcriptional regulator of astrocyte reactivity - as an important mediator of MDD-specific chromatin accessibility and gene expression. Genetic manipulations in mouse OFC demonstrated that astrocytic Zbtb7a is both necessary and sufficient to promote behavioral deficits, cell-type-specific transcriptional and chromatin profiles, and OFC neuronal hyperexcitability induced by chronic stress - a major risk factor for MDD. These data thus highlight a critical role for OFC astrocytes in stress vulnerability and pinpoint ZBTB7A as a key dysregulated factor in MDD that mediates maladaptive astrocytic functions driving OFC hyperexcitability.

BIND&MODIFY: a long-range method for single-molecule mapping of chromatin modifications in eukaryotes

Epigenetic modifications of histones are associated with development and pathogenesis of disease. Existing approaches cannot provide insights into long-range interactions and represent the average chromatin state. Here we describe BIND&MODIFY, a method using long-read sequencing for profiling histone modifications and transcription factors on individual DNA fibers. We use recombinant fused protein A-M.EcoGII to tether methyltransferase M.EcoGII to protein binding sites to label neighboring regions by methylation. Aggregated BIND&MODIFY signal matches bulk ChIP-seq and CUT&TAG. BIND&MODIFY can simultaneously measure histone modification status, transcription factor binding, and CpG 5mC methylation at single-molecule resolution and also quantifies correlation between local and distal elements.

Integrated analysis reveals common DNA methylation patterns of alcohol-associated cancers: A pan-cancer analysis

Background: The role of alcohol in carcinogenesis has received increasing attention in recent years. Evidence shows its impacts on various aspects, including epigenetics alteration. The DNA methylation patterns underlying alcohol-associated cancers are not fully understood. Methods: We investigated the aberrant DNA methylation patterns in four alcohol-associated cancers based on the Illumina HumanMethylation450 BeadChip. Pearson coefficient correlations were identified between differential methylated CpG probes and annotated genes. Transcriptional factor motifs were enriched and clustered using MEME Suite, and a regulatory network was constructed. Results: In each cancer, differential methylated probes (DMPs) were identified, and 172 hypermethylated and 21 hypomethylated pan-cancer DMPs (PDMPs) were examined further. Annotated genes significantly regulated by PDMPs were investigated and enriched in transcriptional misregulation in cancers. The CpG island chr19:58220189-58220517 was hypermethylated in all four cancers and silenced in the transcription factor ZNF154. Various biological effects were exerted by 33 hypermethylated and seven hypomethylated transcriptional factor motifs grouped into five clusters. Eleven pan-cancer DMPs were identified to be associated with clinical outcomes in the four alcohol-associated cancers, which might provide a potential point of view for clinical outcome prediction. Conclusion: This study provides an integrated insight into DNA methylation patterns in alcohol-associated cancers and reveals the corresponding features, influences, and potential mechanisms.

Genome-wide identification and characterization of parthenocarpic fruit set-related gene homologs in cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.), a major horticultural crop, in the family Cucurbitaceae is grown and consumed globally. Parthenocarpy is an ideal trait for many fruit and vegetables which produces seedless fruit desired by consumers. The seedlessness occurs when fruit develops without fertilization which can be either natural or induced. So far, a limited number of genes regulating parthenocarpic fruit set have been reported in several fruit or vegetable crops, most of which are involved in hormone biosynthesis or signalling. Although parthenocarpic cucumber has been widely used in commercial production for a long time; its genetic basis is not well understood. In this study, we retrieved thirty five parthenocarpy fruit-set related genes (PRGs) from bibliomic data in various plants. Thirty-five PRG homologs were identified in the cucumber genome via homology-based search. An in silico analysis was performed on phylogenetic tree, exon-intron structure, cis-regulatory elements in the promoter region, and conserved domains of their deduced proteins, which provided insights into the genetic make-up of parthenocarpy-related genes in cucumber. Simple sequence repeat (SSR) sequences were mined in these PRGs, and 31 SSR markers were designed. SSR genotyping identified three SSRs in two polymorphic genes. Quantitative real-time PCR of selected genes was conducted in five cucumber lines with varying degrees of parthenocarpic fruit set capacities, which revealed possible association of their expression with parthenocarpy. The results revealed that homologs CsWD40 and CsPIN-4 could be considered potential genes for determination of parthenocarpy as these genes showed parental polymorphism and differential gene expression in case of parthenocarpic and non-parthenocarpic parents.

A circadian clock translational control mechanism targets specific mRNAs to cytoplasmic messenger ribonucleoprotein granules

Phosphorylation of Neurospora crassa eukaryotic initiation factor 2 α (eIF2α), a conserved translation initiation factor, is clock controlled. To determine the impact of rhythmic eIF2α phosphorylation on translation, we performed temporal ribosome profiling and RNA sequencing (RNA-seq) in wild-type (WT), clock mutant Δfrq, eIF2α kinase mutant Δcpc-3, and constitutively active cpc-3c cells. About 14% of mRNAs are rhythmically translated in WT cells, and translation rhythms for ∼30% of these mRNAs, which we named circadian translation-initiation-controlled genes (cTICs), are dependent on the clock and CPC-3. Most cTICs are expressed from arrhythmic mRNAs and contain a P-body (PB) localization motif in their 5' leader sequence. Deletion of SNR-1, a component of cytoplasmic messenger ribonucleoprotein granules (cmRNPgs) that include PBs and stress granules (SGs), and the PB motif on one of the cTIC mRNAs, zip-1, significantly alters zip-1 rhythmic translation. These results reveal that the clock regulates rhythmic translation of specific mRNAs through rhythmic eIF2α activity and cmRNPg metabolism.

Transcriptomic profiling implicates PAF1 in both active and repressive immune regulatory networks

BACKGROUND

Sitting at the interface of gene expression and host-pathogen interaction, polymerase associated factor 1 complex (PAF1C) is a rising player in the innate immune response. The complex localizes to the nucleus and associates with chromatin to modulate RNA polymerase II (RNAPII) elongation of gene transcripts. Performing this function at both proximal and distal regulatory elements, PAF1C interacts with many host factors across such sites, along with several microbial proteins during infection. Therefore, translating the ubiquity of PAF1C into specific impacts on immune gene expression remains especially relevant.

RESULTS

Advancing past work, we treat PAF1 knockout cells with a slate of immune stimuli to identify key trends in PAF1-dependent gene expression with broad analytical depth. From our transcriptomic data, we confirm PAF1 is an activator of traditional immune response pathways as well as other cellular pathways correlated with pathogen defense. With this model, we employ computational approaches to refine how PAF1 may contribute to both gene activation and suppression. Specifically focusing on transcriptional motifs and regulons, we predict gene regulatory elements strongly associated with PAF1, including those implicated in an immune response. Overall, our results suggest PAF1 is involved in innate immunity at several distinct axes of regulation.

CONCLUSIONS

By identifying PAF1-dependent gene expression across several pathogenic contexts, we confirm PAF1C to be a key mediator of innate immunity. Combining these transcriptomic profiles with potential regulatory networks corroborates the previously identified functions of PAF1C. With this, we foster new avenues for its study as a regulator of innate immunity, and our results will serve as a basis for targeted study of PAF1C in future validation studies.

Integrator endonuclease drives promoter-proximal termination at all RNA polymerase II-transcribed loci

RNA polymerase II (RNAPII) pausing in early elongation is critical for gene regulation. Paused RNAPII can be released into productive elongation by the kinase P-TEFb or targeted for premature termination by the Integrator complex. Integrator comprises endonuclease and phosphatase activities, driving termination by cleavage of nascent RNA and removal of stimulatory phosphorylation. We generated a degron system for rapid Integrator endonuclease (INTS11) depletion to probe the direct consequences of Integrator-mediated RNA cleavage. Degradation of INTS11 elicits nearly universal increases in active early elongation complexes. However, these RNAPII complexes fail to achieve optimal elongation rates and exhibit persistent Integrator phosphatase activity. Thus, only short transcripts are significantly upregulated following INTS11 loss, including transcription factors, signaling regulators, and non-coding RNAs. We propose a uniform molecular function for INTS11 across all RNAPII-transcribed loci, with differential effects on particular genes, pathways, or RNA biotypes reflective of transcript lengths rather than specificity of Integrator activity.

RTA1 Is Involved in Resistance to 7-Aminocholesterol and Secretion of Fungal Proteins in Cryptococcus neoformans

Cryptococcus neoformans (Cn) is a pathogenic yeast that is the leading cause of fungal meningitis in immunocompromised patients. Various Cn virulence factors, such as the enzyme laccase and its product melanin, phospholipase, and capsular polysaccharide have been identified. During a screen of knockout mutants, the gene resistance to aminocholesterol 1 (RTA1) was identified, the function of which is currently unknown in Cn. Rta1 homologs in S. cerevisiae belong to a lipid-translocating exporter family of fungal proteins with transmembrane regions and confer resistance to the antimicrobial agent 7-aminocholesterol when overexpressed. To determine the role of RTA1 in Cn, the knock-out (rta1Δ) and reconstituted (rta1Δ+RTA1) strains were created and phenotypically tested. RTA1 was involved in resistance to 7-aminocholesterol, and also in exocyst complex component 3 (Sec6)-mediated secretion of urease, laccase, and the major capsule component, glucuronoxylomannan (GXM), which coincided with significantly smaller capsules in the rta1Δ and rta1Δ+RTA1 strains compared to the wild-type H99 strain. Furthermore, RTA1 expression was reduced in a secretory 14 mutant (sec14Δ) and increased in an RNAi Sec6 mutant. Transmission electron microscopy demonstrated vesicle accumulation inside the rta1Δ strain, predominantly near the cell membrane. Given that Rta1 is likely to be a transmembrane protein located at the plasma membrane, these data suggest that Rta1 may be involved in both secretion of various fungal virulence factors and resistance to 7-aminocholesterol in Cn.

Mutant Ras and inflammation-driven skin tumorigenesis is suppressed via a JNK-iASPP-AP1 axis

Concurrent mutation of a RAS oncogene and the tumor suppressor p53 is common in tumorigenesis, and inflammation can promote RAS-driven tumorigenesis without the need to mutate p53. Here, we show, using a well-established mutant RAS and an inflammation-driven mouse skin tumor model, that loss of the p53 inhibitor iASPP facilitates tumorigenesis. Specifically, iASPP regulates expression of a subset of p63 and AP1 targets, including genes involved in skin differentiation and inflammation, suggesting that loss of iASPP in keratinocytes supports a tumor-promoting inflammatory microenvironment. Mechanistically, JNK-mediated phosphorylation regulates iASPP function and inhibits iASPP binding with AP1 components, such as JUND, via PXXP/SH3 domain-mediated interaction. Our results uncover a JNK-iASPP-AP1 regulatory axis that is crucial for tissue homeostasis. We show that iASPP is a tumor suppressor and an AP1 coregulator.

Integrative System Biology Analysis of Transcriptomic Responses to Drought Stress in Soybean (Glycine max L.)

Drought is a major abiotic stressor that causes yield losses and limits the growing area for most crops. Soybeans are an important legume crop that is sensitive to water-deficit conditions and suffers heavy yield losses from drought stress. To improve drought-tolerant soybean cultivars through breeding, it is necessary to understand the mechanisms of drought tolerance in soybeans. In this study, we applied several transcriptome datasets obtained from soybean plants under drought stress in comparison to those grown under normal conditions to identify novel drought-responsive genes and their underlying molecular mechanisms. We found 2168 significant up/downregulated differentially expressed genes (DEGs) and 8 core modules using gene co-expression analysis to predict their biological roles in drought tolerance. Gene Ontology and KEGG analyses revealed key biological processes and metabolic pathways involved in drought tolerance, such as photosynthesis, glyceraldehyde-3-phosphate dehydrogenase and cytokinin dehydrogenase activity, and regulation of systemic acquired resistance. Genome-wide analysis of plants' cis-acting regulatory elements (CREs) and transcription factors (TFs) was performed for all of the identified DEG promoters in soybeans. Furthermore, the PPI network analysis revealed significant hub genes and the main transcription factors regulating the expression of drought-responsive genes in each module. Among the four modules associated with responses to drought stress, the results indicated that GLYMA_04G209700, GLYMA_02G204700, GLYMA_06G030500, GLYMA_01G215400, and GLYMA_09G225400 have high degrees of interconnection and, thus, could be considered as potential candidates for improving drought tolerance in soybeans. Taken together, these findings could lead to a better understanding of the mechanisms underlying drought responses in soybeans, which may useful for engineering drought tolerance in plants.

Transcriptome data reveal gene clusters and key genes in pepper response to heat shock

Climate change and global warming pose a great threat to plant growth and development as well as crop productivity. To better study the genome-wide gene expression under heat, we performed a time-course (0.5 to 24 h) transcriptome analysis in the leaf and root of 40-day-old pepper plants under 40°C as well as in control plants. Clustering analysis (K-means) showed that the expression of 29,249 genes can be grouped into 12 clusters with distinct expression dynamics under stress. Gene ontology (GO) enrichment analysis and transcription factor (TF) identification were performed on the clusters with certain expression patterns. Comparative analysis between the heat-treated and control plants also identified differentially expressed genes (DEGs), which showed the largest degree of change at 24 h. Interestingly, more DEGs were identified in the root than in the leaf. Moreover, we analyzed the gene expression of 25 heat shock factor genes (HSFs) in pepper after heat stress, identified five of these HSFs that responded to heat stress, and characterized the role of these genes in heat-tolerant (17CL30) and heat-susceptible (05S180) pepper lines. The findings of this study improve our understanding of the genome-wide heat stress response in pepper.

Weighted Gene Correlation Network Analysis (WGCNA) of Arabidopsis Somatic Embryogenesis (SE) and Identification of Key Gene Modules to Uncover SE-Associated Hub Genes

Somatic embryogenesis (SE), which occurs naturally in many plant species, serves as a model to elucidate cellular and molecular mechanisms of embryo patterning in plants. Decoding the regulatory landscape of SE is essential for its further application. Hence, the present study was aimed at employing Weighted Gene Correlation Network Analysis (WGCNA) to construct a gene coexpression network (GCN) for Arabidopsis SE and then identifying highly correlated gene modules to uncover the hub genes associated with SE that may serve as potential molecular targets. A total of 17,059 genes were filtered from a microarray dataset comprising four stages of SE, i.e., stage I (zygotic embryos), stage II (proliferating tissues at 7 days of induction), stage III (proliferating tissues at 14 days of induction), and stage IV (mature somatic embryos). This included 1,711 transcription factors and 445 EMBRYO DEFECTIVE genes. GCN analysis identified a total of 26 gene modules with the module size ranging from 35 to 3,418 genes using a dynamic cut tree algorithm. The module-trait analysis revealed that four, four, seven, and four modules were associated with stages I, II, III, and IV, respectively. Further, we identified a total of 260 hub genes based on the degree of intramodular connectivity. Validation of the hub genes using publicly available expression datasets demonstrated that at least 78 hub genes are potentially associated with embryogenesis; of these, many genes remain functionally uncharacterized thus far. In silico promoter analysis of these genes revealed the presence of cis-acting regulatory elements, "soybean embryo factor 4 (SEF4) binding site," and "E-box" of the napA storage-protein gene of Brassica napus; this suggests that these genes may play important roles in plant embryo development. The present study successfully applied WGCNA to construct a GCN for SE in Arabidopsis and identified hub genes involved in the development of somatic embryos. These hub genes could be used as molecular targets to further elucidate the molecular mechanisms underlying SE in plants.

Integrative mRNA and Long Noncoding RNA Analysis Reveals the Regulatory Network of Floral Bud Induction in Longan (Dimocarpus longan Lour.)

Longan (Dimocarpus longan Lour.) is a tropical/subtropical fruit tree of significant economic importance. Floral induction is an essential process for longan flowering and plays decisive effects on the longan yield. Due to the instability of flowering, it is necessary to understand the molecular mechanisms of floral induction in longan. In this study, mRNA and long noncoding RNA (lncRNA) transcriptome sequencing were performed using the apical buds of fruiting branches as materials. A total of 7,221 differential expressions of mRNAs (DEmRNAs) and 3,238 differential expressions of lncRNAs (DElncRNAs) were identified, respectively. KEGG enrichment analysis of DEmRNAs highlighted the importance of starch and sucrose metabolic, circadian rhythms, and plant hormone signal transduction pathways during floral induction. Combining the analysis of weighted gene co-expression network (WGCNA) and expression pattern of DEmRNAs in the three pathways, specific transcriptional characteristics at each stage during floral induction and regulatory network involving co-expressed genes were investigated. The results showed that sucrose metabolism and auxin signal transduction may be crucial for the growth and maturity of autumn shoots in September and October (B1-B2 stage); starch and sucrose metabolic, circadian rhythms, and plant hormone signal transduction pathways participated in the regulation of floral bud physiological differentiation together in November and December (B3-B4 stage) and the crosstalk among three pathways was also found. Hub genes in the co-expression network and key DEmRNAs in three pathways were identified. The circadian rhythm genes FKF1 and GI were found to activate SOC1gene through the photoperiod core factor COL genes, and they were co-expressed with auxin, gibberellin, abscisic acid, ethylene signaling genes, and sucrose biosynthesis genes at B4 stage. A total of 12 hub-DElncRNAs had potential for positively affecting their distant target genes in three putative key pathways, predominantly in a co-transcriptional manner. A hypothetical model of regulatory pathways and key genes and lncRNAs during floral bud induction in longan was proposed finally. Our studies will provide valuable clues and information to help elucidate the potential molecular mechanisms of floral initiation in longan and woody fruit trees.

Comparative Transcriptomic Analyses of Nitrate-Response in Rice Genotypes With Contrasting Nitrogen Use Efficiency Reveals Common and Genotype-Specific Processes, Molecular Targets and Nitrogen Use Efficiency-Candidates

The genetic basis for nitrogen (N)-response and N use efficiency (NUE) must be found in N-responsive gene expression or protein regulation. Our transcriptomic analysis of nitrate response in two contrasting rice genotypes of Oryza sativa ssp. Indica (Nidhi with low NUE and Panvel1 with high NUE) revealed the processes/functions underlying differential N-response/NUE. The microarray analysis of low nitrate response (1.5 mM) relative to normal nitrate control (15 mM) used potted 21-days old whole plants. It revealed 1,327 differentially expressed genes (DEGs) exclusive to Nidhi and 666 exclusive to Panvel1, apart from 70 common DEGs, of which 10 were either oppositely expressed or regulated to different extents. Gene ontology analyses revealed that photosynthetic processes were among the very few processes common to both the genotypes in low N response. Those unique to Nidhi include cell division, nitrogen utilization, cytoskeleton, etc. in low N-response, whereas those unique to Panvel1 include signal transduction, protein import into the nucleus, and mitochondria. This trend of a few common but mostly unique categories was also true for transporters, transcription factors, microRNAs, and post-translational modifications, indicating their differential involvement in Nidhi and Panvel1. Protein-protein interaction networks constructed using DEG-associated experimentally validated interactors revealed subnetworks involved in cytoskeleton organization, cell wall, etc. in Nidhi, whereas in Panvel1, it was chloroplast development. NUE genes were identified by selecting yield-related genes from N-responsive DEGs and their co-localization on NUE-QTLs revealed the differential distribution of NUE-genes between genotypes but on the same chromosomes 1 and 3. Such hotspots are important for NUE breeders.

DNA N6-Adenine methylation in HBV-related hepatocellular carcinoma

DNA methylation on N6-adenine (6mA) has recently been found to be a potential epigenetic marker in prokaryotes and eukaryotes. However, its distribution patterns and potential functions in human tumorigenesis remain largely unknown. Here, we reported global profiling of 6mA sites in the genome of hepatocellular carcinoma at single-nucleotide resolution using Nanopore sequencing. 6mA was widely distributed throughout the human genome. The 6mA sites were related to the porphyrin and chlorophyll metabolism in autosomes and were related to oxidative phosphorylation and ATP metabolism in mitochondria. AGG was the most significant motif associated with 6mA modification and the prevalent motifs in tumors were mainly distributed in mitochondria. The density of 6mA was related to the activation of gene transcription and 6mA density in repetitive sequences decreased in hepatocellular carcinoma. DNA 6mA methylation modification may also be a potential biomarker for cancer diagnosis and treatment.

Role of lncRNAs in cis- and trans-regulatory responses to salt in Populus trichocarpa

Long non-coding RNAs (lncRNAs) are emerging as versatile regulators in diverse biological processes. However, little is known about their cis- and trans-regulatory contributions in gene expression under salt stress. Using 27 RNA-seq data sets from Populus trichocarpa leaves, stems and roots, we identified 2988 high-confidence lncRNAs, including 1183 salt-induced differentially expressed lncRNAs. Among them, 301 lncRNAs have potential for positively affecting their neighboring genes, predominantly in a cis-regulatory manner rather than by co-transcription. Additionally, a co-expression network identified six striking salt-associated modules with a total of 5639 genes, including 426 lncRNAs, and in these lncRNA sequences, the DNA/RNA binding motifs are enriched. This suggests that lncRNAs might contribute to distant gene expression of the salt-associated modules in a trans-regulatory manner. Moreover, we found 30 lncRNAs that have potential to simultaneously cis- and trans-regulate salt-responsive homologous genes, and Ptlinc-NAC72, significantly induced under long-term salt stress, was selected for validating its regulation of the expression and functional roles of the homologs PtNAC72.A and PtNAC72.B (PtNAC72.A/B). The transient transformation of Ptlinc-NAC72 and a dual-luciferase assay of Ptlinc-NAC72 and PtNAC72.A/B promoters confirmed that Ptlinc-NAC72 can directly upregulate PtNAC72.A/B expression, and a presence/absence assay was further conducted to show that the regulation is probably mediated by Ptlinc-NAC72 recognizing the tandem elements (GAAAAA) in the PtNAC72.A/B 5' untranslated region (5'-UTR). Finally, the overexpression of Ptlinc-NAC72 produces a hypersensitive phenotype under salt stress. Altogether, our results shed light on the cis- and trans-regulation of gene expression by lncRNAs in Populus and provides an example of long-term salt-induced Ptlinc-NAC72 that could be used to mitigate growth costs by conferring plant resilience to salt stress.

Advanced genomics identifies growth effectors for proteotoxic ER stress recovery in Arabidopsis thaliana

Adverse environmental and pathophysiological situations can overwhelm the biosynthetic capacity of the endoplasmic reticulum (ER), igniting a potentially lethal condition known as ER stress. ER stress hampers growth and triggers a conserved cytoprotective signaling cascade, the unfolded protein response (UPR) for ER homeostasis. As ER stress subsides, growth is resumed. Despite the pivotal role of the UPR in growth restoration, the underlying mechanisms for growth resumption are yet unknown. To discover these, we undertook a genomics approach in the model plant species Arabidopsis thaliana and mined the gene reprogramming roles of the UPR modulators, basic leucine zipper28 (bZIP28) and bZIP60, in ER stress resolution. Through a network modeling and experimental validation, we identified key genes downstream of the UPR bZIP-transcription factors (bZIP-TFs), and demonstrated their functional roles. Our analyses have set up a critical pipeline for functional gene discovery in ER stress resolution with broad applicability across multicellular eukaryotes.

Ko and Brandizzi use Arabidopsis thaliana to investigate the downstream regulators of two major endoplasmic reticulum (ER) stress-related transcription factors, bZIP60 and bZIP28. Their results provide further insight on how two modulators of the unfolded protein response contribute to growth recovery from ER stress.

Heterogeneity in the Epigenetic Landscape of Murine Testis-Specific Histone Variants TH2A and TH2B Sharing the Same Bi-Directional Promoter

Testis-specific histone variants are crucial to promote open chromatin structure to enable nucleosome disassembly in the final stages of spermiogenesis. However, even after histone replacement, mature sperm retain a proportion of these variants, the function of which is unknown. The present study aimed to understand the functional relevance of the retained H2B and H2A variants, TH2B and TH2A. While no literature is available on the phenotype of TH2A knockouts, TH2B/TH2A double knockout male mice are reported to be infertile. In this study, ChIP-seq analysis was done for TH2B and TH2A to understand the epigenomics of the retained TH2B and TH2A, using murine caudal sperm. Distribution across genomic partitions revealed ∼35% of the TH2B peaks within ±5 kb of TSS whereas TH2A peaks distribution was sparse at TSS. Gene Ontology revealed embryo development as the most significant term associated with TH2B. Also, based on genomic regions, TH2B was observed to be associated with spindle assembly and various meiosis-specific genes, which is an important finding as TH2A/TH2B DKO mice have been reported to have defective cohesin release. A comparison of mouse and human TH2B-linked chromatin revealed 26% overlap between murine and human TH2B-associated genes. This overlap included genes crucial for embryogenesis. Most importantly, heterogeneity in the epigenetic landscape of TH2A and TH2B was seen, which is intriguing as TH2B and TH2A are well reported to be present in the same nucleosomes to promote open chromatin. Additionally, unlike TH2B, TH2A was enriched on the mitochondrial chromosome. TH2A was found to be associated with Nuclear insertion of Mitochondrial DNA sequences (NUMTs) in sperm. A comprehensive analysis of these observations indicates novel functions for the sperm-retained TH2B and TH2A.

Comprehensive analysis reveals common DNA methylation patterns of tobacco-associated cancers: A pan-cancer analysis

The role of tobacco in carcinogenesis has received increasing attention across a number of disciplines in recent years. Accumulating evidences reveal that tobacco consumption affects various epigenetic modifications, especially DNA methylation. However, the genetic modifications of methylation patterns involved in tobacco-attributable cancers remain poorly understood. In this manuscript, aberrant DNA methylation patterns were investigated in 9 tobacco-attributable cancers. Differential methylated probes (DMPs) were identified in each cancer type and a total of 2,392 hyper- and 736 hypomethylated pan-cancer DMPs (PDMPs) were screened out for further analysis. PDMP-associated genes were mostly enriched in metabolism-associated pathways, suggesting the potential roles of methylation alternation in reprogramming cancer cell metabolism. Hypomethylated PDMPs cg12422154, cg02772121 and cg06051311 constituted an enhancer region, significantly downregulating TRIM15, TRIM26 and RPP21, which serve as epigenetically therapeutic biomarkers. Forty-three hypermethylated and 13 hypomethylated transcription factor motifs were clustered into 6 groups, and exhibited various biological functions. Forty-nine PDMPs were reported to be associated with prognosis, providing effective tools to predict clinical outcomes. In summary, our studies revealed the characteristics, influences and potential mechanisms of DNA methylation patterns of tobacco-attributable cancer.

Cholecalciferol and metformin protect against lipopolysaccharide-induced endothelial dysfunction and senescence by modulating sirtuin-1 and protein arginine methyltransferase-1

Endothelial cell activation through nuclear factor-kappa-B (NFkB) and mitogen-activated protein kinases leads to increased biosynthesis of pro-inflammatory mediators, cellular injury and vascular inflammation under lipopolysaccharide (LPS) exposure. Recent studies report that LPS up-regulated global methyltransferase activity. In this study, we observed that a combination treatment with metformin (MET) and cholecalciferol (VD) blocked the LPS-induced S-adenosylmethionine (SAM)-dependent methyltransferase (SDM) activity in Eahy926 cells. We found that LPS challenge (i) increased arginine methylation through up-regulated protein arginine methyltransferase-1 (PRMT1) mRNA, intracellular concentrations of asymmetric dimethylarginine (ADMA) and homocysteine (HCY); (ii) up-regulated cell senescence through mitigated sirtuin-1 (SIRT1) mRNA, nicotinamide adenine dinucleotide (NAD+) concentration, telomerase activity and total antioxidant capacity; and (iii) lead to endothelial dysfunction through compromised nitric oxide (NOx) production. However, these LPS-mediated cellular events in Eahy926 cells were restored by the synergistic effect of MET and VD. Taken together, this study identified that the dual compound effect inhibits LPS-induced protein arginine methylation, endothelial senescence and dysfunction through the components of epigenetic machinery, SIRT1 and PRMT1, which is a previously unidentified function of the test compounds. In silico results identified the presence of vitamin D response element (VDRE) sequence on PRMT1 suggesting that VDR could regulate PRMT1 gene expression. Further characterization of the cellular events associated with the dual compound challenge, using gene silencing approach or adenoviral constructs for SIRT1 and/or PRMT1 under inflammatory stress, could identify therapeutic strategies to address the endothelial consequences in vascular inflammation-mediated atherosclerosis.

A global screening identifies chromatin-enriched RNA-binding proteins and the transcriptional regulatory activity of QKI5 during monocytic differentiation

BACKGROUND

Cellular RNA-binding proteins (RBPs) have multiple roles in post-transcriptional control, and some are shown to bind DNA. However, the global localization and the general chromatin-binding ability of RBPs are not well-characterized and remain undefined in hematopoietic cells.

RESULTS

We first provide a full view of RBPs' distribution pattern in the nucleus and screen for chromatin-enriched RBPs (Che-RBPs) in different human cells. Subsequently, by generating ChIP-seq, CLIP-seq, and RNA-seq datasets and conducting combined analysis, the transcriptional regulatory potentials of certain hematopoietic Che-RBPs are predicted. From this analysis, quaking (QKI5) emerges as a potential transcriptional activator during monocytic differentiation. QKI5 is over-represented in gene promoter regions, independent of RNA or transcription factors. Furthermore, DNA-bound QKI5 activates the transcription of several critical monocytic differentiation-associated genes, including CXCL2, IL16, and PTPN6. Finally, we show that the differentiation-promoting activity of QKI5 is largely dependent on CXCL2, irrespective of its RNA-binding capacity.

CONCLUSIONS

Our study indicates that Che-RBPs are versatile factors that orchestrate gene expression in different cellular contexts, and identifies QKI5, a classic RBP regulating RNA processing, as a novel transcriptional activator during monocytic differentiation.

Memes: A motif analysis environment in R using tools from the MEME Suite

Identification of biopolymer motifs represents a key step in the analysis of biological sequences. The MEME Suite is a widely used toolkit for comprehensive analysis of biopolymer motifs; however, these tools are poorly integrated within popular analysis frameworks like the R/Bioconductor project, creating barriers to their use. Here we present memes, an R package that provides a seamless R interface to a selection of popular MEME Suite tools. memes provides a novel “data aware” interface to these tools, enabling rapid and complex discriminative motif analysis workflows. In addition to interfacing with popular MEME Suite tools, memes leverages existing R/Bioconductor data structures to store the multidimensional data returned by MEME Suite tools for rapid data access and manipulation. Finally, memes provides data visualization capabilities to facilitate communication of results. memes is available as a Bioconductor package at https://bioconductor.org/packages/memes, and the source code can be found at github.com/snystrom/memes.

Biologically active molecules such as DNA, RNA, and proteins are polymers composed of repeated subunits. For example, nucleotides are the subunits of DNA and RNA, and amino acids are the subunits of proteins. Functional properties of biopolymers are determined by short, recurring stretches of subunits known as motifs. Motifs can serve as binding sites between molecules, they can influence the structure of molecules, and they can contribute to enzymatic activities. For these reasons, motif analysis has become a key step in determining the function of biopolymers and in elucidating their roles in biological networks. The MEME Suite is a widely used compilation of tools used for identifying and analyzing motifs found in biological sequences. Here, we describe a new piece of software named “memes” that connects MEME Suite tools to R, the statistical analysis environment. By providing an interface between the MEME Suite and R, memes allows for improved motif analysis workflows and easy access to a wide range of existing data visualization tools, further expanding the utility of MEME Suite tools.

Genome-Wide Transcription Factor DNA Binding Sites and Gene Regulatory Networks in Clostridium thermocellum

Clostridium thermocellum is a thermophilic bacterium recognized for its natural ability to effectively deconstruct cellulosic biomass. While there is a large body of studies on the genetic engineering of this bacterium and its physiology to-date, there is limited knowledge in the transcriptional regulation in this organism and thermophilic bacteria in general. The study herein is the first report of a large-scale application of DNA-affinity purification sequencing (DAP-seq) to transcription factors (TFs) from a bacterium. We applied DAP-seq to > 90 TFs in C. thermocellum and detected genome-wide binding sites for 11 of them. We then compiled and aligned DNA binding sequences from these TFs to deduce the primary DNA-binding sequence motifs for each TF. These binding motifs are further validated with electrophoretic mobility shift assay (EMSA) and are used to identify individual TFs’ regulatory targets in C. thermocellum. Our results led to the discovery of novel, uncharacterized TFs as well as homologues of previously studied TFs including RexA-, LexA-, and LacI-type TFs. We then used these data to reconstruct gene regulatory networks for the 11 TFs individually, which resulted in a global network encompassing the TFs with some interconnections. As gene regulation governs and constrains how bacteria behave, our findings shed light on the roles of TFs delineated by their regulons, and potentially provides a means to enable rational, advanced genetic engineering of C. thermocellum and other organisms alike toward a desired phenotype.

Hyperglycemia Induces Trained Immunity in Macrophages and Their Precursors and Promotes Atherosclerosis

BACKGROUND

Cardiovascular risk in diabetes remains elevated despite glucose-lowering therapies. We hypothesized that hyperglycemia induces trained immunity in macrophages, promoting persistent proatherogenic characteristics.

METHODS

Bone marrow-derived macrophages from control mice and mice with diabetes were grown in physiological glucose (5 mmol/L) and subjected to RNA sequencing (n=6), assay for transposase accessible chromatin sequencing (n=6), and chromatin immunoprecipitation sequencing (n=6) for determination of hyperglycemia-induced trained immunity. Bone marrow transplantation from mice with (n=9) or without (n=6) diabetes into (normoglycemic) Ldlr-/- mice was used to assess its functional significance in vivo. Evidence of hyperglycemia-induced trained immunity was sought in human peripheral blood mononuclear cells from patients with diabetes (n=8) compared with control subjects (n=16) and in human atherosclerotic plaque macrophages excised by laser capture microdissection.

RESULTS

In macrophages, high extracellular glucose promoted proinflammatory gene expression and proatherogenic functional characteristics through glycolysis-dependent mechanisms. Bone marrow-derived macrophages from diabetic mice retained these characteristics, even when cultured in physiological glucose, indicating hyperglycemia-induced trained immunity. Bone marrow transplantation from diabetic mice into (normoglycemic) Ldlr-/- mice increased aortic root atherosclerosis, confirming a disease-relevant and persistent form of trained innate immunity. Integrated assay for transposase accessible chromatin, chromatin immunoprecipitation, and RNA sequencing analyses of hematopoietic stem cells and bone marrow-derived macrophages revealed a proinflammatory priming effect in diabetes. The pattern of open chromatin implicated transcription factor Runt-related transcription factor 1 (Runx1). Similarly, transcriptomes of atherosclerotic plaque macrophages and peripheral leukocytes in patients with type 2 diabetes were enriched for Runx1 targets, consistent with a potential role in human disease. Pharmacological inhibition of Runx1 in vitro inhibited the trained phenotype.

CONCLUSIONS

Hyperglycemia-induced trained immunity may explain why targeting elevated glucose is ineffective in reducing macrovascular risk in diabetes and suggests new targets for disease prevention and therapy.

Evolutionary History of Inversions in Directional Mutational Pressures in Crustacean Mitochondrial Genomes: Implications for Evolutionary Studies

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.