An improved ChIP-seq peak detection system for simultaneously identifying post-translational modified transcription factors by combinatorial fusion, using SUMOylation as an example

Visualizing set relationships: EVenn's comprehensive approach to Venn diagrams

Venn diagrams serve as invaluable tools for visualizing set relationships due to their ease of interpretation. Widely applied across diverse disciplines such as metabolomics, genomics, transcriptomics, and proteomics, their utility is undeniable. However, the operational complexity has been compounded by the absence of standardized data formats and the need to switch between various platforms for generating different Venn diagrams. To address these challenges, we introduce the EVenn platform, a versatile tool offering a unified interface for efficient data exploration and visualization of diverse Venn diagrams. EVenn (http://www.ehbio.com/test/venn) streamlines the data upload process with a standardized format, enhancing the capabilities for multimodule analysis. This comprehensive protocol outlines various applications of EVenn, featuring representative results of multiple Venn diagrams, data uploads in the centralized data center, and step-by-step case demonstrations. Through these functionalities, EVenn emerges as a valuable and user-friendly tool for the in-depth exploration of multiomics data.

De novo transcriptome assembly of Dalbergia sissoo Roxb. (Fabaceae) under Botryodiplodia theobromae-induced dieback disease

Dalbergia sissoo Roxb. (Shisham) is a timber-producing species of economic, cultural, and medicinal importance in the Indian subcontinent. In the past few decades, Shisham's dieback disease caused by the fungus Botryodiplodia theobromae has become an evolving issue in the subcontinent endangering its survival. To gain insights into this issue, a standard transcriptome assembly was deployed to assess the response of D. sissoo at the transcriptomic level under the stress of B. theobromae infection. For RNA isolation, the control and infected leaf tissue samples were taken from 1-year-old greenhouse-grown D. sissoo plants after 20 days of stem-base spore inoculation. cDNA synthesis was performed from these freshly isolated RNA samples that were then sent for sequencing. About 18.14 Gb (Giga base) of data was generated using the BGISEQ-500 sequencing platform. In terms of Unigenes, 513,821 were identified after a combined assembly of all samples and then filtering the abundance. The total length of Unigenes, their average length, N50, and GC-content were 310,523,693 bp, 604 bp, 1,101 bp, and 39.95% respectively. The Unigenes were annotated using 7 functional databases i.e., 200,355 (NR: 38.99%), 164,973 (NT: 32.11%), 123,733 (Swissprot: 24.08%), 142,580 (KOG: 27.75%), 139,588 (KEGG: 27.17%), 99,752 (GO: 19.41%), and 137,281 (InterPro: 26.72%). Furthermore, the Transdecoder detected 115,762 CDS. In terms of SSR (Simple Sequence Repeat) markers, 62,863 of them were distributed on 51,508 Unigenes and on the predicted 4673 TF (Transcription Factor) coding Unigenes. A total of 16,018 up- and 19,530 down-regulated Differentially Expressed Genes (DEGs) were also identified. Moreover, the Plant Resistance Genes (PRGs) had a count of 9230. We are hopeful that in the future, these identified Unigenes, SSR markers, DEGs and PRGs will provide the prerequisites for managing Shisham dieback disease, its breeding, and in tree improvement programs.

Using machine learning to predict protein-protein interactions between a zombie ant fungus and its carpenter ant host

Parasitic fungi produce proteins that modulate virulence, alter host physiology, and trigger host responses. These proteins, classified as a type of "effector," often act via protein-protein interactions (PPIs). The fungal parasite Ophiocordyceps camponoti-floridani (zombie ant fungus) manipulates Camponotus floridanus (carpenter ant) behavior to promote transmission. The most striking aspect of this behavioral change is a summit disease phenotype where infected hosts ascend and attach to an elevated position. Plausibly, interspecific PPIs drive aspects of Ophiocordyceps infection and host manipulation. Machine learning PPI predictions offer high-throughput methods to produce mechanistic hypotheses on how this behavioral manipulation occurs. Using D-SCRIPT to predict host-parasite PPIs, we found ca. 6000 interactions involving 2083 host proteins and 129 parasite proteins, which are encoded by genes upregulated during manipulated behavior. We identified multiple overrepresentations of functional annotations among these proteins. The strongest signals in the host highlighted neuromodulatory G-protein coupled receptors and oxidation-reduction processes. We also detected Camponotus structural and gene-regulatory proteins. In the parasite, we found enrichment of Ophiocordyceps proteases and frequent involvement of novel small secreted proteins with unknown functions. From these results, we provide new hypotheses on potential parasite effectors and host targets underlying zombie ant behavioral manipulation.

Bidirectional promoters: an enigmatic genome architecture and their roles in cancers

Bidirectional promoters are the transcription regulatory regions of genes positioned head-to-head on opposite strands. Specific sequence signals, chromatin modifications and three-dimensional structures of the transcription site facilitate the unconventional yet tightly regulated transcription proceeding in both directions from these promoters. Mutations or aberrant epigenetic changes can lead to abnormal enhanced or reduced expression from either of the bidirectionally transcribed genes resulting in tumorigenesis. Moreover, bidirectionally transcribed genes might also contribute towards the immune regulation in tumor microenvironment. In this review, we aimed to expound the characteristic features of bidirectional promoters alongside their transcriptional regulations, and ultimately, the association of these enigmatic genomic elements in different cancers.

A conserved role for transcription factor sumoylation in binding-site selection

Large numbers of eukaryotic transcription factors (TFs) are modified by SUMO post-translational modifications. Whereas the effect of TF sumoylation on the expression of target genes is largely context-dependent, it is not known whether the modification has a common function in regulating TFs throughout eukaryotic species. Here, I highlight four studies that used genome-wide chromatin-immunoprecipitation analysis (ChIP-seq) to examine whether sumoylation affects the selection of sites on the genome that are bound by human and yeast TFs. The studies found that impairing sumoylation led to deregulated binding-site selection for all four of the examined TFs. Predominantly, compared to wild-type forms, the sumoylation-deficient forms of the TFs bound to numerous additional non-specific sites, pointing to a common role for the modification in restricting TF binding to appropriate sites. Evidence from these studies suggests that TF sumoylation influences binding-site selection by modulating protein-protein interactions with other DNA-binding TFs, or by promoting conformational changes in the TFs that alter their DNA-binding specificity or affinity. I propose a model in which, prior to their sumoylation, TFs initially bind to chromatin with reduced specificity, which leads to spurious binding but also ensures that all functional sites become bound. Once the TFs are bound to DNA, sumoylation then acts to increase specificity and promotes release of the TFs from non-specific sites. The similar observations from these four genome-wide studies across divergent species suggest that binding-site selection is a general and conserved function for TF sumoylation.

DeSUMOylation of MKK7 kinase by the SUMO2/3 protease SENP3 potentiates lipopolysaccharide-induced inflammatory signaling in macrophages

Protein SUMOylation has been reported to play a role in innate immune response, but the enzymes, substrates, and consequences of the specific inflammatory signaling events are largely unknown. Reactive oxygen species (ROS) are abundantly produced during macrophage activation and required for Toll-like receptor 4 (TLR4)-mediated inflammatory signaling. Previously, we demonstrated that SENP3 is a redox-sensitive SUMO2/3 protease. To explore any links between reversible SUMOylation and ROS-related inflammatory signaling in macrophage activation, we generated mice with Senp3 conditional knock-out in myeloid cells. In bacterial lipopolysaccharide (LPS)-induced in vitro and in vivo inflammation models, we found that SENP3 deficiency markedly compromises the activation of TLR4 inflammatory signaling and the production of proinflammatory cytokines in macrophages exposed to LPS. Moreover, Senp3 conditional knock-out mice were significantly less susceptible to septic shock. Of note, SENP3 deficiency was associated with impairment in JNK phosphorylation. We found that MKK7, which selectively phosphorylates JNK, is a SENP3 substrate and that SENP3-mediated deSUMOylation of MKK7 may favor its binding to JNK. Importantly, ROS-dependent SENP3 accumulation and MKK7 deSUMOylation rapidly occurred after LPS stimulation. In conclusion, our findings indicate that SENP3 potentiates LPS-induced TLR4 signaling via deSUMOylation of MKK7 leading to enhancement in JNK phosphorylation and the downstream events. Therefore this work provides novel mechanistic insights into redox regulation of innate immune responses.

Comparative secretome analysis of Colletotrichum falcatum identifies a cerato-platanin protein (EPL1) as a potential pathogen-associated molecular pattern (PAMP) inducing systemic resistance in sugarcane

Colletotrichum falcatum, an intriguing hemibiotrophic fungal pathogen causes red rot, a devastating disease of sugarcane. Repeated in vitro subculturing of C. falcatum under dark condition alters morphology and reduces virulence of the culture. Hitherto, no information is available on this phenomenon at molecular level. In this study, the in vitro secretome of C. falcatum cultured under light and dark conditions was analyzed using 2-DE coupled with MALDI TOF/TOF MS. Comparative analysis identified nine differentially abundant proteins. Among them, seven proteins were less abundant in the dark-cultured C. falcatum, wherein only two protein species of a cerato-platanin protein called EPL1 (eliciting plant response-like protein) were found to be highly abundant. Transcriptional expression of candidate high abundant proteins was profiled during host-pathogen interaction using qRT-PCR. Comprehensively, this comparative secretome analysis identified five putative effectors, two pathogenicity-related proteins and one pathogen-associated molecular pattern (PAMP) of C. falcatum. Functional characterization of three distinct domains of the PAMP (EPL1) showed that the major cerato-platanin domain (EPL1∆N1-92) is exclusively essential for inducing defense and hypersensitive response (HR) in sugarcane and tobacco, respectively. Further, priming with EPL1∆N1-92 protein induced systemic resistance and significantly suppressed the red rot severity in sugarcane.

BIOLOGICAL SIGNIFICANCE

Being the first secretomic investigation of C. falcatum, this study has identified five potential effectors, two pathogenicity-related proteins and a PAMP. Although many reports have highlighted the influence of light on pathogenicity, this study has established a direct link between light and expression of effectors, for the first time. This study has presented the influence of a novel N-terminal domain of EPL1 in physical and biological properties and established the functional role of major cerato-platanin domain of EPL1 as a potential elicitor inducing systemic resistance in sugarcane. Comprehensively, the study has identified proteins that putatively contribute to virulence of C. falcatum and for the first time, demonstrated the potential role of EPL1 in inducing PAMP-triggered immunity (PTI) in sugarcane.

Novel insights into mitochondrial gene rearrangement in thrips (Insecta: Thysanoptera) from the grass thrips, Anaphothrips obscurus

We sequenced the mitochondrial (mt) genome of the grass thrips, Anaphothrips obscurus, which is highly rearranged and differs from the four thrips species reported previously in the arrangement of both tRNA genes and a protein-coding gene, nad3, and in the copy number of the control region (CR). We reconstructed the phylogeny of the thrips with mt genome sequences, and used it as a framework to gain insights into mt genome evolution in thrips. It is evident that A. obscurus is less rearranged in mt genome organization than the other four known thrips. nad3 is in its ancestral location in A. obscurus but was translocated in other four thrips. Also, A. obscurus has one CR, which is ancestral to hexapods whereas other thrips have two or three CRs. All of the five thrips whose mt genomes have been sequenced to date are from the subfamily Thripinae, which represents about a quarter of the species richness in the order Thysanoptera. The high variation in mt genome organization observed in a subfamily challenges our knowledge about animal mt genomes. It remains to be investigated why mt genomes evolved so fast in the subfamily Thripinae and how mt genomes evolved in other lineages of thrips.

Exosomal microRNAs in giant panda (Ailuropoda melanoleuca) breast milk: potential maternal regulators for the development of newborn cubs

The physiological role of miRNAs is widely understood to include fine-tuning the post-transcriptional regulation of a wide array of biological processes. Extensive studies have indicated that exosomal miRNAs in the bodily fluids of various organisms can be transferred between living cells for the delivery of gene silencing signals. Here, we illustrated the expression characteristics of exosomal miRNAs in giant panda breast milk during distinct lactation periods and highlighted the enrichment of immune- and development-related endogenous miRNAs in colostral and mature giant panda milk. These miRNAs are stable, even under certain harsh conditions, via the protection of extracellular vesicles. These findings indicate that breast milk may facilitate the dietary intake of maternal miRNAs by infants for the regulation of postnatal development. We also detected exogenous plant miRNAs from the primary food source of the giant panda (bamboo) in the exosomes of giant panda breast milk that were associated with regulatory roles in basic metabolism and neuron development. This result suggested that dietary plant miRNAs are absorbed by host cells and subsequently secreted into bodily fluids as potential cross-kingdom regulators. In conclusion, exosomal miRNAs in giant panda breast milk may be crucial maternal regulators for the development of intrinsic 'slink' newborn cubs.

SUMO and the robustness of cancer

Post-translational protein modification by small ubiquitin-like modifier (SUMO), termed sumoylation, is an important mechanism in cellular responses to stress and one that appears to be upregulated in many cancers. Here, we examine the role of sumoylation in tumorigenesis as a possibly necessary safeguard that protects the stability and functionality of otherwise easily misregulated gene expression programmes and signalling pathways of cancer cells.

HuMiChip2 for strain level identification and functional profiling of human microbiomes

With the massive data generated by the Human Microbiome Project, how to transform such data into useful information and knowledge remains challenging. Here, with currently available sequencing information (reference genomes and metagenomes), we have developed a comprehensive microarray, HuMiChip2, for strain-level identification and functional characterization of human microbiomes. HuMiChip2 was composed of 29,467 strain-specific probes targeting 2063 microbial strains/species and 133,924 sequence- and group-specific probes targeting 157 key functional gene families involved in various metabolic pathways and host-microbiome interaction processes. Computational evaluation of strain-specific probes suggested that they were not only specific to mock communities of sequenced microorganisms and metagenomes from different human body sites but also to non-sequenced microbial strains. Experimental evaluation of strain-specific probes using single strains/species and mock communities suggested a high specificity of these probes with their corresponding targets. Application of HuMiChip2 to human gut microbiome samples showed the patient microbiomes of alcoholic liver cirrhosis significantly (p < 0.05) shifted their functional structure from the healthy individuals, and the relative abundance of 21 gene families significantly (p < 0.1) differed between the liver cirrhosis patients and healthy individuals. At the strain level, five Bacteroides strains were significantly (p < 0.1) and more frequently detected in liver cirrhosis patients. These results suggest that the developed HuMiChip2 is a useful microbial ecological microarray for both strain-level identification and functional profiling of human microbiomes.

In low protein diets, microRNA-19b regulates urea synthesis by targeting SIRT5

Ammonia detoxification, which takes place via the hepatic urea cycle, is essential for nitrogen homeostasis and physiological well-being. It has been reported that a reduction in dietary protein reduces urea nitrogen. MicroRNAs (miRNAs) are major regulatory non-coding RNAs that have significant effects on several metabolic pathways; however, little is known on whether miRNAs regulate hepatic urea synthesis. The objective of this study was to assess the miRNA expression profile in a low protein diet and identify miRNAs involved in the regulation of the hepatic urea cycle using a porcine model. Weaned 28-days old piglets were fed a corn-soybean normal protein diet (NP) or a corn-soybean low protein diet (LP) for 30 d. Hepatic and blood samples were collected, and the miRNA expression profile was assessed by sequencing and qRT-PCR. Furthermore, we evaluated the possible role of miR-19b in urea synthesis regulation. There were 25 differentially expressed miRNAs between the NP and LP groups. Six of these miRNAs were predicted to be involved in urea cycle metabolism. MiR-19b negatively regulated urea synthesis by targeting SIRT5, which is a positive regulator of CPS1, the rate limiting enzyme in the urea cycle. Our study presented a novel explanation of ureagenesis regulation by miRNAs.

Analysis of tigecycline resistance development in clinical Acinetobacter baumannii isolates through a combined genomic and transcriptomic approach

Tigecycline (Tgc) is considered a last-resort antibiotic for the treatment of multi-drug resistant bacteria. To study Tgc resistance development in the important nosocomial pathogen Acinetobacter baumannii, we adopted six clinical isolates from three patients undergoing antibiotic treatment, and bacterial genomic sequences and seven strand-specific transcriptomes were studied. Interestingly, the Tgc-intermediate 2015ZJAB1 only differed from Tgc-resistant 2015ZJAB2 in an SNP-clustered region including OprD, a sugar-type MFS permease, and a LuxR-type transcriptional regulator. Surprisingly, an almost identical region was found in 2015ZJAB3, which supports the possibility of a homologous recombination event that increased Tgc resistance. Furthermore, comparative transcriptomic analysis identified significantly regulated genes associated with Tgc resistance, which was verified using qRT-PCR. Three enriched COG categories included amino acid transport and metabolism, transcription, and inorganic ion transport and metabolism. KEGG analysis revealed common features under Tgc conditions, including up regulated benzoate degradation and a less active TCA cycle. This may be related to selective antimicrobial pressure in the environment and adaptation by lowering metabolism. This study provides the first report of an in vivo evolutionary process that included a putative homologous recombination event conferring Tgc resistance in clinical A. baumannii isolates in which transcriptome analysis revealed resistance-conferring genes and related metabolism characteristics.

The chromatin modification by SUMO-2/3 but not SUMO-1 prevents the epigenetic activation of key immune-related genes during Kaposi's sarcoma associated herpesvirus reactivation

BACKGROUND

SUMOylation, as part of the epigenetic regulation of transcription, has been intensively studied in lower eukaryotes that contain only a single SUMO protein; however, the functions of SUMOylation during mammalian epigenetic transcriptional regulation are largely uncharacterized. Mammals express three major SUMO paralogues: SUMO-1, SUMO-2, and SUMO-3 (normally referred to as SUMO-1 and SUMO-2/3). Herpesviruses, including Kaposi's sarcoma associated herpesvirus (KSHV), seem to have evolved mechanisms that directly or indirectly modulate the SUMO machinery in order to evade host immune surveillance, thus advancing their survival. Interestingly, KSHV encodes a SUMO E3 ligase, K-bZIP, with specificity toward SUMO-2/3 and is an excellent model for investigating the global functional differences between SUMO paralogues.

RESULTS

We investigated the effect of experimental herpesvirus reactivation in a KSHV infected B lymphoma cell line on genomic SUMO-1 and SUMO-2/3 binding profiles together with the potential role of chromatin SUMOylation in transcription regulation. This was carried out via high-throughput sequencing analysis. Interestingly, chromatin immunoprecipitation sequencing (ChIP-seq) experiments showed that KSHV reactivation is accompanied by a significant increase in SUMO-2/3 modification around promoter regions, but SUMO-1 enrichment was absent. Expression profiling revealed that the SUMO-2/3 targeted genes are primarily highly transcribed genes that show no expression changes during viral reactivation. Gene ontology analysis further showed that these genes are involved in cellular immune responses and cytokine signaling. High-throughput annotation of SUMO occupancy of transcription factor binding sites (TFBS) pinpointed the presence of three master regulators of immune responses, IRF-1, IRF-2, and IRF-7, as potential SUMO-2/3 targeted transcriptional factors after KSHV reactivation.

CONCLUSION

Our study is the first to identify differential genome-wide SUMO modifications between SUMO paralogues during herpesvirus reactivation. Our findings indicate that SUMO-2/3 modification near protein-coding gene promoters occurs in order to maintain host immune-related gene unaltered during viral reactivation.