The CUT&RUN greenlist: genomic regions of consistent noise are effective normalizing factors for quantitative epigenome mapping.

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a recent development for epigenome mapping, but its unique methodology can hamper proper quantitative analyses. As traditional normalization approaches have been shown to be inaccurate, we sought to determine endogenous normalization factors based on the human genome regions of constant nonspecific signal. This constancy was determined by applying Shannon's information entropy, and the set of normalizer regions, which we named the 'Greenlist', was extensively validated using publicly available datasets. We demonstrate here that the greenlist normalization outperforms the current top standards, and remains consistent across different experimental setups, cell lines and antibodies; the approach can even be applied to different species or to CUT&Tag. Requiring no additional experimental steps and no added cost, this approach can be universally applied to CUT&RUN experiments to greatly minimize the interference of technical variation over the biological epigenome changes of interest.

The Human Developing Cerebral Cortex Is Characterized by an Elevated De Novo Expression of Long Noncoding RNAs in Excitatory Neurons.

The outstanding human cognitive capacities are computed in the cerebral cortex, a mammalian-specific brain region and the place of massive biological innovation. Long noncoding RNAs have emerged as gene regulatory elements with higher evolutionary turnover than mRNAs. The many long noncoding RNAs identified in neural tissues make them candidates for molecular sources of cerebral cortex evolution and disease. Here, we characterized the genomic and cellular shifts that occurred during the evolution of the long noncoding RNA repertoire expressed in the developing cerebral cortex and explored putative roles for these long noncoding RNAs in the evolution of the human brain. Using transcriptomics and comparative genomics, we comprehensively annotated the cortical transcriptomes of humans, rhesus macaques, mice, and chickens and classified human cortical long noncoding RNAs into evolutionary groups as a function of their predicted minimal ages. Long noncoding RNA evolutionary groups showed differences in expression levels, splicing efficiencies, transposable element contents, genomic distributions, and transcription factor binding to their promoters. Furthermore, older long noncoding RNAs showed preferential expression in germinative zones, outer radial glial cells, and cortical inhibitory (GABAergic) neurons. In comparison, younger long noncoding RNAs showed preferential expression in cortical excitatory (glutamatergic) neurons, were enriched in primate and human-specific gene co-expression modules, and were dysregulated in neurodevelopmental disorders. These results suggest different evolutionary routes for older and younger cortical long noncoding RNAs, highlighting old long noncoding RNAs as a possible source of molecular evolution of conserved developmental programs; conversely, we propose that the de novo expression of primate- and human-specific young long noncoding RNAs is a putative source of molecular evolution and dysfunction of cortical excitatory neurons, warranting further investigation.

Long non-coding RNAs are essential for Schistosoma mansoni pairing-dependent adult worm homeostasis and fertility.

The trematode parasite Schistosoma mansoni causes schistosomiasis, which affects over 200 million people worldwide. Schistosomes are dioecious, with egg laying depending on the females' obligatory pairing with males. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with low or no protein-coding potential that have been involved in other species with reproduction, stem cell maintenance, and drug resistance. In S. mansoni, we recently showed that the knockdown of one lncRNA affects the pairing status of these parasites. Here, we re-analyzed public RNA-Seq data from paired and unpaired adult male and female worms and their gonads, obtained from mixed-sex or single-sex cercariae infections, and found thousands of differentially expressed pairing-dependent lncRNAs among the 23 biological samples that were compared. The expression levels of selected lncRNAs were validated by RT-qPCR using an in vitro unpairing model. In addition, the in vitro silencing of three selected lncRNAs showed that knockdown of these pairing-dependent lncRNAs reduced cell proliferation in adult worms and their gonads, and are essential for female vitellaria maintenance, reproduction, and/or egg development. Remarkably, in vivo silencing of each of the three selected lncRNAs significantly reduced worm burden in infected mice by 26 to 35%. Whole mount in situ hybridization experiments showed that these pairing-dependent lncRNAs are expressed in reproductive tissues. These results show that lncRNAs are key components intervening in S. mansoni adult worm homeostasis, which affects pairing status and survival in the mammalian host, thus presenting great potential as new therapeutic target candidates.

Dementia is an age-independent risk factor for severity and death in COVID-19 inpatients.

INTRODUCTION: Dementia has been associated with COVID-19 prevalence, but whether this reflects higher infection, older age of patients, or disease severity remains unclear. METHODS: We investigated a cohort of 12,863 UK Biobank community-dwelling individuals > 65 years old (1814 individuals ≥ 80 years old) tested for COVID-19. Individuals were stratified by age to account for age as a confounder. Risk factors were analyzed for COVID-19-positive diagnosis, hospitalization, and death. RESULTS: All-cause dementia, Alzheimer's disease (AD), and Parkinson's disease (PD) were associated with COVID-19-positive diagnosis, and all-cause dementia and AD remained associated in individuals ≥ 80 years old. All-cause dementia, AD, or PD were not risk factors for overall hospitalization, but increased the risk of hospitalization of COVID-19 patients. All-cause dementia and AD increased the risk of COVID-19-related death, and all-cause dementia was uniquely associated with increased death in ≥ 80-year-old patients. DISCUSSION: All-cause dementia and AD are age-independent risk factors for disease severity and death in COVID-19.

Global analysis of biogenesis, stability and sub-cellular localization of lncRNAs mapping to intragenic regions of the human genome.

Long noncoding RNAs (lncRNAs) that map to intragenic regions of the human genome with the same (intronic lncRNAs) or opposite orientation (antisense lncRNAs) relative to protein-coding mRNAs have been largely dismissed from biochemical and functional characterization due to the belief that they are mRNA precursors, byproducts of RNA splicing or simply transcriptional noise. In this work, we used a custom microarray to investigate aspects of the biogenesis, processing, stability, evolutionary conservation, and cellular localization of ∼ 6,000 intronic lncRNAs and ∼ 10,000 antisense lncRNAs. Most intronic (2,903 of 3,427, 85%) and antisense lncRNAs (4,945 of 5,214, 95%) expressed in HeLa cells showed evidence of 5' cap modification, compatible with their transcription by RNAP II. Antisense lncRNAs (median t1/2 = 3.9 h) were significantly (p < 0.0001) more stable than mRNAs (median t1/2 = 3.2 h), whereas intronic lncRNAs (median t1/2 = 2.1 h) comprised a more heterogeneous class that included both stable (t1/2 > 3 h) and unstable (t1/2 < 1 h) transcripts. Intragenic lncRNAs display evidence of evolutionary conservation, have little/no coding potential and were ubiquitously detected in the cytoplasm. Notably, a fraction of the intronic and antisense lncRNAs (13 and 15%, respectively) were expressed from loci at which the corresponding host mRNA was not detected. The abundances of a subset of intronic/antisense lncRNAs were correlated (r ≥ |0.8|) with those of genes encoding proteins involved in cell division and DNA replication. Taken together, the findings of this study contribute novel biochemical and genomic information regarding intronic and antisense lncRNAs, supporting the notion that these classes include independently transcribed RNAs with potentials for exerting regulatory functions in the cell.

Schistosoma mansoni vaccine candidates identified by unbiased phage display screening in self-cured rhesus macaques.

Schistosomiasis, a challenging neglected tropical disease, affects millions of people worldwide. Developing a prophylactic vaccine against Schistosoma mansoni has been hindered by the parasite's biological complexity. In this study, we utilized the innovative phage-display immunoprecipitation followed by a sequencing approach (PhIP-Seq) to screen the immune response of 10 infected rhesus macaques during self-cure and challenge-resistant phases, identifying vaccine candidates. Our high-throughput S. mansoni synthetic DNA phage-display library encoded 99.6% of 119,747 58-mer peptides, providing comprehensive coverage of the parasite's proteome. Library screening with rhesus macaques' antibodies, from the early phase of establishment of parasite infection, identified significantly enriched epitopes of parasite extracellular proteins known to be expressed in the digestive tract, shifting towards intracellular proteins during the late phase of parasite clearance. Immunization of mice with a selected pool of PhIP-Seq-enriched phage-displayed peptides from MEG proteins, cathepsins B, and asparaginyl endopeptidase significantly reduced worm burden in a vaccination assay. These findings enhance our understanding of parasite-host immune responses and provide promising prospects for developing an effective schistosomiasis vaccine.

Expression analysis and in silico characterization of intronic long noncoding RNAs in renal cell carcinoma: emerging functional associations.

BACKGROUND: Intronic and intergenic long noncoding RNAs (lncRNAs) are emerging gene expression regulators. The molecular pathogenesis of renal cell carcinoma (RCC) is still poorly understood, and in particular, limited studies are available for intronic lncRNAs expressed in RCC. METHODS: Microarray experiments were performed with custom-designed arrays enriched with probes for lncRNAs mapping to intronic genomic regions. Samples from 18 primary RCC tumors and 11 nontumor adjacent matched tissues were analyzed. Meta-analyses were performed with microarray expression data from three additional human tissues (normal liver, prostate tumor and kidney nontumor samples), and with large-scale public data for epigenetic regulatory marks and for evolutionarily conserved sequences. RESULTS: A signature of 29 intronic lncRNAs differentially expressed between RCC and nontumor samples was obtained (false discovery rate (FDR) < 5%). A signature of 26 intronic lncRNAs significantly correlated with the RCC five-year patient survival outcome was identified (FDR < 5%, p-value ≤ 0.01). We identified 4303 intronic antisense lncRNAs expressed in RCC, of which 22% were significantly (p < 0.05) cis correlated with the expression of the mRNA in the same locus across RCC and three other human tissues. Gene Ontology (GO) analysis of those loci pointed to 'regulation of biological processes' as the main enriched category. A module map analysis of the protein-coding genes significantly (p < 0.05) trans correlated with the 20% most abundant lncRNAs, identified 51 enriched GO terms (p < 0.05). We determined that 60% of the expressed lncRNAs are evolutionarily conserved. At the genomic loci containing the intronic RCC-expressed lncRNAs, a strong association (p < 0.001) was found between their transcription start sites and genomic marks such as CpG islands, RNA Pol II binding and histones methylation and acetylation. CONCLUSION: Intronic antisense lncRNAs are widely expressed in RCC tumors. Some of them are significantly altered in RCC in comparison with nontumor samples. The majority of these lncRNAs is evolutionarily conserved and possibly modulated by epigenetic modifications. Our data suggest that these RCC lncRNAs may contribute to the complex network of regulatory RNAs playing a role in renal cell malignant transformation.

Transcriptional profiling of Neurospora crassa Δmak-2 reveals that mitogen-activated protein kinase MAK-2 participates in the phosphate signaling pathway.

The filamentous fungus Neurospora crassa is an excellent model system for examining molecular responses to ambient signals in eukaryotic microorganisms. Inorganic phosphate (Pi) is an essential growth-limiting nutrient in nature and is crucial for the synthesis of nucleic acids and the flow of genetic information. The genetic and molecular mechanisms controlling the response to Pi starvation in N. crassa include at least four genes (nuc-2, preg, pogv, and nuc-1), which are involved in a hierarchical regulatory activation network. In a previous work, we identified a number of genes modulated by NUC-2 protein, including the mak-2 gene, which codes for a mitogen-activated protein kinase (MAPK), suggesting its participation in the phosphate signaling pathway. Thus, to identify other genes involved in metabolic responses to exogenous phosphate sensing and the functioning of the MAPK MAK-2, we performed microarray experiments using a mak-2 knockout strain (Δmak-2) grown under phosphate-shortage conditions by comparing its transcription profile to that of a control strain grown in low- and high-phosphate cultures. These experiments revealed 912 unique differentially expressed genes involved in a number of physiological processes related to phosphate transport, metabolism, and regulation as well as posttranslational modification of proteins, and MAPK signaling pathways. Quantitative Real-time PCR gene expression analysis of 18 selected genes, using independent RNA samples, validated our microarray results. A high Pearson correlation between microarray and quantitative Real-time PCR data was observed. The analysis of these differentially expressed genes in the Δmak-2 strain provide evidence that the mak-2 gene participates in the hierarchical phosphate-signaling pathway in N. crassa in addition to its involvement in other metabolic routes such as the isoprenylation pathway, thus revealing novel aspects of the N. crassa phosphorus-sensing network.

Single-cell RNA-seq analyses show that long non-coding RNAs are conspicuously expressed in Schistosoma mansoni gamete and tegument progenitor cell populations.

Schistosoma mansoni is a flatworm that causes schistosomiasis, a neglected tropical disease that affects over 200 million people worldwide. New therapeutic targets are needed with only one drug available for treatment and no vaccine. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with low or no protein-coding potential. In other organisms, they have been shown as involved with reproduction, stem cell maintenance and drug resistance, and they tend to exhibit tissue-specific expression patterns. S. mansoni expresses thousands of lncRNA genes; however, the cell type expression patterns of lncRNAs in the parasite remain uncharacterized. Here, we have re-analyzed publicly available single-cell RNA-sequencing (scRNA-seq) data obtained from adult S. mansoni to identify the lncRNAs signature of adult schistosome cell types. A total of 8023 lncRNAs (79% of all lncRNAs) were detected. Analyses of the lncRNAs expression profiles in the cells using statistically stringent criteria were performed to identify 74 lncRNA gene markers of cell clusters. Male gamete and tegument progenitor lineages clusters contained most of the cluster-specific lncRNA markers. We also identified lncRNA markers of specific neural clusters. Whole-mount in situ hybridization (WISH) and double fluorescence in situ hybridization were used to validate the cluster-specific expression of 13 out of 16 selected lncRNA genes (81%) in the male and female adult parasite tissues; for one of these 16 gene loci, probes for two different lncRNA isoforms were used, which showed differential isoform expression in testis and ovary. An atlas of the expression profiles across the cell clusters of all lncRNAs detected in our analysis is available as a public website resource (http://verjolab.usp.br:8081). The results presented here give strong support to a tissue-specific expression and to a regulated expression program of lncRNAs in S. mansoni. This will be the basis for further exploration of lncRNA genes as potential therapeutic targets.

Long noncoding intronic RNAs are differentially expressed in primary and metastatic pancreatic cancer.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is known by its aggressiveness and lack of effective therapeutic options. Thus, improvement in current knowledge of molecular changes associated with pancreatic cancer is urgently needed to explore novel venues of diagnostics and treatment of this dismal disease. While there is mounting evidence that long noncoding RNAs (lncRNAs) transcribed from intronic and intergenic regions of the human genome may play different roles in the regulation of gene expression in normal and cancer cells, their expression pattern and biological relevance in pancreatic cancer is currently unknown. In the present work we investigated the relative abundance of a collection of lncRNAs in patients' pancreatic tissue samples aiming at identifying gene expression profiles correlated to pancreatic cancer and metastasis. METHODS: Custom 3,355-element spotted cDNA microarray interrogating protein-coding genes and putative lncRNA were used to obtain expression profiles from 38 clinical samples of tumor and non-tumor pancreatic tissues. Bioinformatics analyses were performed to characterize structure and conservation of lncRNAs expressed in pancreatic tissues, as well as to identify expression signatures correlated to tissue histology. Strand-specific reverse transcription followed by PCR and qRT-PCR were employed to determine strandedness of lncRNAs and to validate microarray results, respectively. RESULTS: We show that subsets of intronic/intergenic lncRNAs are expressed across tumor and non-tumor pancreatic tissue samples. Enrichment of promoter-associated chromatin marks and over-representation of conserved DNA elements and stable secondary structure predictions suggest that these transcripts are generated from independent transcriptional units and that at least a fraction is under evolutionary selection, and thus potentially functional.Statistically significant expression signatures comprising protein-coding mRNAs and lncRNAs that correlate to PDAC or to pancreatic cancer metastasis were identified. Interestingly, loci harboring intronic lncRNAs differentially expressed in PDAC metastases were enriched in genes associated to the MAPK pathway. Orientation-specific RT-PCR documented that intronic transcripts are expressed in sense, antisense or both orientations relative to protein-coding mRNAs. Differential expression of a subset of intronic lncRNAs (PPP3CB, MAP3K14 and DAPK1 loci) in metastatic samples was confirmed by Real-Time PCR. CONCLUSION: Our findings reveal sets of intronic lncRNAs expressed in pancreatic tissues whose abundance is correlated to PDAC or metastasis, thus pointing to the potential relevance of this class of transcripts in biological processes related to malignant transformation and metastasis in pancreatic cancer.

The evolution of TNF signaling in platyhelminths suggests the cooptation of TNF receptor in the host-parasite interplay.

BACKGROUND: The TNF signaling pathway is involved in the regulation of many cellular processes (such as apoptosis and cell proliferation). Previous reports indicated the effect of human TNF-α on metabolism, physiology, gene expression and protein phosphorylation of the human parasite Schistosoma mansoni and suggested that its TNF receptor was responsible for this response. The lack of an endogenous TNF ligand reinforced the idea of the use of an exogenous ligand, but also opens the possibility that the receptor actually binds a non-canonical ligand, as observed for NGFRs. METHODS: To obtain a more comprehensive view, we analyzed platyhelminth genomes deposited in the Wormbase ParaSite database to investigate the presence of TNF receptors and their respective ligands. Using different bioinformatics approaches, such as HMMer and BLAST search tools we identified and characterized the sequence of TNF receptors and ligand homologs. We also used bioinformatics resources for the identification of conserved protein domains and Bayesian inference for phylogenetic analysis. RESULTS: Our analyses indicate the presence of 31 TNF receptors in 30 platyhelminth species. All platyhelminths display a single TNF receptor, and all are structurally remarkably similar to NGFR. It suggests no events of duplication and diversification occurred in this phylum, with the exception of a single species-specific duplication. Interestingly, we also identified TNF ligand homologs in five species of free-living platyhelminths. CONCLUSIONS: These results suggest that the TNF receptor from platyhelminths may be able to bind canonical TNF ligands, thus strengthening the idea that these receptors are able to bind human TNF-α. This also raises the hypothesis that an endogenous ligand was substituted by the host ligand in parasitic platyhelminths. Moreover, our analysis indicates that death domains (DD) may be present in the intracellular region of most platyhelminth TNF receptors, thus pointing to a previously unreported apoptotic action of such receptors in platyhelminths. Our data highlight the idea that host-parasite crosstalk using the TNF pathway may be widespread in parasitic platyhelminths to mediate apoptotic responses. This opens up a new hypothesis to uncover what might be an important component to understand platyhelminth infections.

Genome-wide DNA methylation profile in the peripheral blood of cocaine and crack dependents.

OBJECTIVE: Cocaine use disorders (CUDs) represent a major public health problem in many countries. To better understand the interaction between the environmental modulations and phenotype, the aim of the present study was to investigate the DNA methylation pattern of CUD patients, who had concomitant cocaine and crack dependence, and healthy controls. METHODS: We studied DNA methylation profiles in the peripheral blood of 23 CUD patients and 24 healthy control subjects using the Illumina Infinium HumanMethylation450 BeadChip arrays. RESULTS: Comparison between CUD patients and controls revealed 186 differentially methylated positions (DMPs; adjusted p-value [adjP] < 10-5) related to 152 genes, with a subset of CpGs confirmed by pyrosequencing. DNA methylation patterns discriminated CUD patients and control groups. A gene network approach showed that the EHMT1, EHMT2, MAPK1, MAPK3, MAP2K1, and HDAC5 genes, which are involved in transcription and chromatin regulation cellular signaling pathways, were also associated with cocaine dependence. CONCLUSION: The investigation of DNA methylation patterns may contribute to a better understanding of the biological mechanisms involved in CUD.

Co-expression network of neural-differentiation genes shows specific pattern in schizophrenia.

BACKGROUND: Schizophrenia is a neurodevelopmental disorder with genetic and environmental factors contributing to its pathogenesis, although the mechanism is unknown due to the difficulties in accessing diseased tissue during human neurodevelopment. The aim of this study was to find neuronal differentiation genes disrupted in schizophrenia and to evaluate those genes in post-mortem brain tissues from schizophrenia cases and controls. METHODS: We analyzed differentially expressed genes (DEG), copy number variation (CNV) and differential methylation in human induced pluripotent stem cells (hiPSC) derived from fibroblasts from one control and one schizophrenia patient and further differentiated into neuron (NPC). Expression of the DEG were analyzed with microarrays of post-mortem brain tissue (frontal cortex) cohort of 29 schizophrenia cases and 30 controls. A Weighted Gene Co-expression Network Analysis (WGCNA) using the DEG was used to detect clusters of co-expressed genes that were non-conserved between adult cases and controls brain samples. RESULTS: We identified methylation alterations potentially involved with neuronal differentiation in schizophrenia, which displayed an over-representation of genes related to chromatin remodeling complex (adjP = 0.04). We found 228 DEG associated with neuronal differentiation. These genes were involved with metabolic processes, signal transduction, nervous system development, regulation of neurogenesis and neuronal differentiation. Between adult brain samples from cases and controls there were 233 DEG, with only four genes overlapping with the 228 DEG, probably because we compared single cell to tissue bulks and more importantly, the cells were at different stages of development. The comparison of the co-expressed network of the 228 genes in adult brain samples between cases and controls revealed a less conserved module enriched for genes associated with oxidative stress and negative regulation of cell differentiation. CONCLUSION: This study supports the relevance of using cellular approaches to dissect molecular aspects of neurogenesis with impact in the schizophrenic brain. We showed that, although generated by different approaches, both sets of DEG associated to schizophrenia were involved with neocortical development. The results add to the hypothesis that critical metabolic changes may be occurring during early neurodevelopment influencing faulty development of the brain and potentially contributing to further vulnerability to the illness.

Genome-wide DNA methylation profile in the peripheral blood of cocaine and crack dependents

Objective: Cocaine use disorders (CUDs) represent a major public health problem in many countries. To better understand the interaction between the environmental modulations and phenotype, the aim of the present study was to investigate the DNA methylation pattern of CUD patients, who had concomitant cocaine and crack dependence, and healthy controls. Methods: We studied DNA methylation profiles in the peripheral blood of 23 CUD patients and 24 healthy control subjects using the Illumina Infinium HumanMethylation450 BeadChip arrays. Results: Comparison between CUD patients and controls revealed 186 differentially methylated positions (DMPs; adjusted p-value [adjP] < 10-5) related to 152 genes, with a subset of CpGs confirmed by pyrosequencing. DNA methylation patterns discriminated CUD patients and control groups. A gene network approach showed that the EHMT1, EHMT2, MAPK1, MAPK3, MAP2K1, and HDAC5 genes, which are involved in transcription and chromatin regulation cellular signaling pathways, were also associated with cocaine dependence. Conclusion: The investigation of DNA methylation patterns may contribute to a better understanding of the biological mechanisms involved in CUD.