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

Abstract 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.We would like to thank Helder I Nakaya and Tarik A El-Jundi for help with the 44 k microarray experiments and, together with Rodrigo Louro, for helpful discussions. This work was supported by grants from Fundação de Amparo a Pesquisa do Estado de São Paulo (FAPESP) and from Instituto Nacional de Ciência e Tecnologia em Oncogenômica to SVA and EMR, and by fellowships from FAPESP to AAF, ACT, SAVA and VMC. SVA and EMR received research fellowship awards from Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq), Brazil

Sex differences in DNA methylation of the cord blood are related to sex-bias psychiatric diseases

Sex differences in the prevalence of psychiatric disorders are well documented, with exposure to stress during gestation differentially impacting females and males. We explored sex-specific DNA methylation in the cord blood of 39 females and 32 males born at term and with appropriate weight at birth regarding their potential connection to psychiatric outcomes. Mothers were interviewed to gather information about environmental factors (gestational exposure) that could interfere with the methylation profiles in the newborns. Bisulphite converted DNA was hybridized to Illumina HumanMethylation450 BeadChips. Excluding XYS probes, there were 2,332 differentially methylated CpG sites (DMSs) between sexes, which were enriched within brain modules of co-methylated CpGs during brain development and also differentially methylated in the brains of boys and girls. Genes associated with the DMSs were enriched for neurodevelopmental disorders, particularly for CpG sites found differentially methylated in brain tissue between patients with schizophrenia and controls. Moreover, the DMS had an overlap of 890 (38%) CpG sites with a cohort submitted to toxic exposition during gestation. This study supports the evidences that sex differences in DNA methylation of autosomes act as a primary driver of sex differences that are found in psychiatric outcomes

Global burden of chronic respiratory diseases and risk factors, 1990–2019: an update from the Global Burden of Disease Study 2019

Background: Updated data on chronic respiratory diseases (CRDs) are vital in their prevention, control, and treatment in the path to achieving the third UN Sustainable Development Goals (SDGs), a one-third reduction in premature mortality from non-communicable diseases by 2030. We provided global, regional, and national estimates of the burden of CRDs and their attributable risks from 1990 to 2019. Methods: Using data from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2019, we estimated mortality, years lived with disability, years of life lost, disability-adjusted life years (DALYs), prevalence, and incidence of CRDs, i.e. chronic obstructive pulmonary disease (COPD), asthma, pneumoconiosis, interstitial lung disease and pulmonary sarcoidosis, and other CRDs, from 1990 to 2019 by sex, age, region, and Socio-demographic Index (SDI) in 204 countries and territories. Deaths and DALYs from CRDs attributable to each risk factor were estimated according to relative risks, risk exposure, and the theoretical minimum risk exposure level input. Findings: In 2019, CRDs were the third leading cause of death responsible for 4.0 million deaths (95% uncertainty interval 3.6–4.3) with a prevalence of 454.6 million cases (417.4–499.1) globally. While the total deaths and prevalence of CRDs have increased by 28.5% and 39.8%, the age-standardised rates have dropped by 41.7% and 16.9% from 1990 to 2019, respectively. COPD, with 212.3 million (200.4–225.1) prevalent cases, was the primary cause of deaths from CRDs, accounting for 3.3 million (2.9–3.6) deaths. With 262.4 million (224.1–309.5) prevalent cases, asthma had the highest prevalence among CRDs. The age-standardised rates of all burden measures of COPD, asthma, and pneumoconiosis have reduced globally from 1990 to 2019. Nevertheless, the age-standardised rates of incidence and prevalence of interstitial lung disease and pulmonary sarcoidosis have increased throughout this period. Low- and low-middle SDI countries had the highest age-standardised death and DALYs rates while the high SDI quintile had the highest prevalence rate of CRDs. The highest deaths and DALYs from CRDs were attributed to smoking globally, followed by air pollution and occupational risks. Non-optimal temperature and high body-mass index were additional risk factors for COPD and asthma, respectively. Interpretation: Albeit the age-standardised prevalence, death, and DALYs rates of CRDs have decreased, they still cause a substantial burden and deaths worldwide. The high death and DALYs rates in low and low-middle SDI countries highlights the urgent need for improved preventive, diagnostic, and therapeutic measures. Global strategies for tobacco control, enhancing air quality, reducing occupational hazards, and fostering clean cooking fuels are crucial steps in reducing the burden of CRDs, especially in low- and lower-middle income countries

Expression analysis of long noncoding RNAs in pancreatic adecarcinoma

RNAs não codificadores longos (lncRNAs) compõem uma fração significativa do transcriptoma. Alterações na expressão de lncRNAs já foram observadas em vários cânceres humanos, mas ainda não foram exploradas no adenocarcinoma pancreático ductal (PDAC), uma doença devastadora e agressiva para a qual faltam métodos para diagnóstico precoce e tratamentos efetivos. Utilizando uma plataforma de microarranjo de cDNA com sondas para 984 lncRNAs e 2371 mRNAs, o presente estudo identificou conjuntos de lncRNAs expressos em 38 amostras clínicas pancreáticas. O enriquecimento de (i) elementos regulatórios associados às regiões promotoras (H3K4me3); (ii) possíveis inícios de transcrição (CAGE-tags); (iii) presença de elementos conservados sugere que ao menos uma fração desses RNAs seja originada a partir de unidades transcricionais independentes, reguladas e possivelmente funcionais. Foram identificadas assinaturas de expressão gênica compostas por mRNA e lncRNAs associadas ao tumor primário e à metástase pancreática. A assinatura gIenica associada à metástase apresentou enriquecimento RNAs intrônicos de loci gênicos associados à via MAPK quinase. O aumento de expressão dos transcritos intrônicos dos loci PPP3CB, MAP3K14 e DAPK1 foi confirmado por qPCR em metástases. Em conjunto, este trabalho aponta para a importância de lncRNAs intrônicos no PDAC e para a necessidade de estudos mais aprofundados para uma melhor compreensão do papel dessa classe de transcritos na biologia da doença.Long noncoding RNAs (lncRNAs) compose a significant fraction of transcriptome. Altered expression of lncRNAs has been observed in diverse human cancers, but has not being investigated in pancreatic ductal adenocarcinoma (PDAC), a devastating and aggressive disease that lack early diagnosis methods and effective treatments. Using a cDNA microarray platform with probes interrogating 984 lncRNAs and 2371 mRNA, the present study identified subsets of lncRNAs expressed in 38 pancreatic clinical samples. Enrichment of (i) regulatory elements associated to promoter region (H3K4me3); (ii) putative transcription start site (CAGEtags) and (iii) conserved elements, suggest that at least a fraction of these RNAs could be independent transcriptional unit, regulated, an possibly functional. Gene expression signatures comprised of mRNAs and lncRNAs and associated to primary or metastatic tumors were found. A gene signature associated to metastasis was enriched in intronic ncRNAs mapping to gene loci associated to the MAPK pathway. Over expression of intronic RNAs from PPP3CB, MAP3K14 and DAPK1 was confirmed by qPCR in metastatic samples. Taken together, this study points to the importance of intronic lncRNAs in PDAC and for the need to study this class of ncRNAs in greater detail to better understand its role in the biology of PDAC

Assessment of complementarity of WGCNA and NERI results for identification of modules associated to schizophrenia spectrum disorders.

Psychiatric disorders involve both changes in multiple genes as well different types of variations. As such, gene co-expression networks allowed the comparison of different stages and parts of the brain contributing to an integrated view of genetic variation. Two methods based on co-expression networks presents appealing results: Weighted Gene Correlation Network Analysis (WGCNA) and Network-Medicine Relative Importance (NERI). By selecting two different gene expression databases related to schizophrenia, we evaluated the biological modules selected by both WGCNA and NERI along these databases as well combining both WGCNA and NERI results (WGCNA-NERI). Also we conducted a enrichment analysis for the identification of partial biological function of each result (as well a replication analysis). To appraise the accuracy of whether both algorithms (as well our approach, WGCNA-NERI) were pointing to genes related to schizophrenia and its complex genetic architecture we conducted the MSET analysis, based on a reference gene list of schizophrenia database (SZDB) related to DNA Methylation, Exome, GWAS as well as copy number variation mutation studies. The WGCNA results were more associated with inflammatory pathways and immune system response; NERI obtained genes related with cellular regulation, embryological pathways e cellular growth factors. Only NERI were able to provide a statistical meaningful results to the MSET analysis (for Methylation and de novo mutations data). However, combining WGCNA and NERI provided a much more larger overlap in these two categories and additionally on Transcriptome database. Our study suggests that using both methods in combination is better for establishing a group of modules and pathways related to a complex disease than using each method individually. NERI is available at: https://bitbucket.org/sergionery/neri

Putative contributions of the sex chromosome proteins SOX3 and SRY to neurodevelopmental disorders

Coordenação de Aperfeiçoamento de Pessoalde Nível Superior, Grant/Award Number: DSProgramDS-1750212PROEX-1669479PROEX-33002010073P7; Fundaçãode Amparo à Pesquisa do Estado de São Paulo,Grant/Award Number: 2011/04956-62011/14658-22014/00041-12014/00591-12014/10488-32015/06281-7; Universidade Federaldo ABC, Grant/Award Number: InstitutionalScholarship; UFABC; CAPES, Grant/AwardNumber: DS-1750212; FAPESP, Grant/AwardNumbers: 2014/10488-3, 2011/04956-6,2014/00591-1, 2014/00041-1,2015/06281-7, 2011/14658-2.Universidade de São Paulo. Faculdade de Medicina. Instituto de Psiquiatria, Hospital das Clinicas HCFMUSP. Laboratório de Psicopatologia e Terapêutica Psiquiátrica . São Paulo, SP, Brazil.Universidade de São Paulo. Faculdade de Medicina. Instituto de Psiquiatria, Hospital das Clinicas HCFMUSP. Laboratório de Psicopatologia e Terapêutica Psiquiátrica . São Paulo, SP, Brazil / University of São Paulo. Inter-institutional Grad Program on Bioinformatics. São Paulo, SP, Brazil.Federal University of ABC. Center of Mathematics, Computing and Cognition. Santo André, SP, Brazil.Universidade de São Paulo. Faculdade de Medicina. Instituto de Psiquiatria, Hospital das Clinicas HCFMUSP. Laboratório de Psicopatologia e Terapêutica Psiquiátrica . São Paulo, SP, Brazil.Universidade de São Paulo. Faculdade de Medicina. Instituto de Psiquiatria, Hospital das Clinicas HCFMUSP. Laboratório de Psicopatologia e Terapêutica Psiquiátrica . São Paulo, SP, Brazil.Brazilian Center for Research in Energy and Materials. Brazilian Biosciences National Laboratory. Campinas, SP, Brazil.Universidade de São Paulo. Faculdade de Medicina. Instituto de Psiquiatria, Hospital das Clinicas HCFMUSP. Laboratório de Psicopatologia e Terapêutica Psiquiátrica . São Paulo, SP, Brazil.Universidade de São Paulo. Faculdade de Medicina. Instituto de Psiquiatria, Hospital das Clinicas HCFMUSP. Laboratório de Psicopatologia e Terapêutica Psiquiátrica . São Paulo, SP, Brazil.Universidade de São Paulo. Faculdade de Medicina. Instituto de Psiquiatria, Hospital das Clinicas HCFMUSP. Laboratório de Psicopatologia e Terapêutica Psiquiátrica . São Paulo, SP, Brazil / Ministério da Saúde. Secretaria de Vigilância em Saúde. Instituto Evandro Chagas. Laboratório de Pesquisas Básicas em Malária - Entomologia. Ananindeua, PA, Brasil.Brazilian Center for Research in Energy and Materials. Brazilian Biosciences National Laboratory. Campinas, SP, Brazil.Universidade de São Paulo. Faculdade de Medicina. Instituto de Psiquiatria, Hospital das Clinicas HCFMUSP. Laboratório de Psicopatologia e Terapêutica Psiquiátrica . São Paulo, SP, Brazil / University of São Paulo. Inter-institutional Grad Program on Bioinformatics. São Paulo, SP, Brazil / Universidade de São Paulo. Faculdade de Medicina. Instituto de Psiquiatria. Hospital das Clinicas HCFMUSP. São Paulo, SP, Brazil / National Institute of Developmental Psychiatry for Children and Adolescents (INPD). São Paulo, SP, Brazil / Universidade de São Paulo. Faculdade de Medicina FMUSP. São Paulo, SP, Brazil.The male-biased prevalence of certain neurodevelopmental disorders and the sex-biased outcomes associated with stress exposure during gestation have been previously described. Here, we hypothesized that genes distinctively targeted by only one or both homologous proteins highly conserved across therian mammals, SOX3 and SRY, could induce sexual adaptive changes that result in a differential risk for neurodevelopmental disorders. ChIP-seq/chip data showed that SOX3/SRY gene targets were expressed in different brain cell types in mice. We used orthologous human genes in rodent genomes to extend the number of SOX3/SRY set (1,721). These genes were later found to be enriched in five modules of coexpressed genes during the early and mid-gestation periods (FDR < 0.05), independent of sexual hormones. Genes with differential expression (24, p < 0.0001) and methylation (40, p < 0.047) between sexes were overrepresented in this set. Exclusive SOX3 or SRY target genes were more associated with the late gestational and postnatal periods. Using autism as a model sex-biased disorder, the SOX3/SRY set was enriched in autism gene databases (FDR ≤ 0.05), and there were more de novo variations from the male autism spectrum disorder (ASD) samples under the SRY peaks compared to the random peaks (p < 0.024). The comparison of coexpressed networks of SOX3/SRY target genes between male autism and control samples revealed low preservation in gene modules related to stress response (99 genes) and neurogenesis (78 genes). This study provides evidence that while SOX3 is a regulatory mechanism for both sexes, the male-exclusive SRY also plays a role in gene regulation, suggesting a potential mechanism for sex bias in ASD

Environmental Influences Measured by Epigenetic Clock and Vulnerability Components at Birth Impact Clinical ASD Heterogeneity

Although Autism Spectrum Disorders (ASD) is recognized as being heavily influenced by genetic factors, the role of epigenetic and environmental factors is still being established. This study aimed to identify ASD vulnerability components based on familial history and intrauterine environmental stress exposure, explore possible vulnerability subgroups, access DNA methylation age acceleration (AA) as a proxy of stress exposure during life, and evaluate the association of ASD vulnerability components and AA to phenotypic severity measures. Principal Component Analysis (PCA) was used to search the vulnerability components from 67 mothers of autistic children. We found that PC1 had a higher correlation with psychosocial stress (maternal stress, maternal education, and social class), and PC2 had a higher correlation with biological factors (psychiatric family history and gestational complications). Comparing the methylome between above and below PC1 average subgroups we found 11,879 statistically significant differentially methylated probes (DMPs, p &lt; 0.05). DMPs CpG sites were enriched in variably methylated regions (VMRs), most showing environmental and genetic influences. Hypermethylated probes presented higher rates in different regulatory regions associated with functional SNPs, indicating that the subgroups may have different affected regulatory regions and their liability to disease explained by common variations. Vulnerability components score moderated by epigenetic clock AA was associated with Vineland Total score (p = 0.0036, adjR2 = 0.31), suggesting risk factors with stress burden can influence ASD phenotype

Integrative Variation Analysis Reveals that a Complex Genotype May Specify Phenotype in Siblings with Syndromic Autism Spectrum Disorder

<div><p>It has been proposed that copy number variations (CNVs) are associated with increased risk of autism spectrum disorder (ASD) and, in conjunction with other genetic changes, contribute to the heterogeneity of ASD phenotypes. Array comparative genomic hybridization (aCGH) and exome sequencing, together with systems genetics and network analyses, are being used as tools for the study of complex disorders of unknown etiology, especially those characterized by significant genetic and phenotypic heterogeneity. Therefore, to characterize the complex genotype-phenotype relationship, we performed aCGH and sequenced the exomes of two affected siblings with ASD symptoms, dysmorphic features, and intellectual disability, searching for <i>de novo</i> CNVs, as well as for <i>de novo</i> and rare inherited point variations—single nucleotide variants (SNVs) or small insertions and deletions (indels)—with probable functional impacts. With aCGH, we identified, in both siblings, a duplication in the 4p16.3 region and a deletion at 8p23.3, inherited by a paternal balanced translocation, t(4, 8) (p16; p23). Exome variant analysis found a total of 316 variants, of which 102 were shared by both siblings, 128 were in the male sibling exome data, and 86 were in the female exome data. Our integrative network analysis showed that the siblings’ shared translocation could explain their similar syndromic phenotype, including overgrowth, macrocephaly, and intellectual disability. However, exome data aggregate genes to those already connected from their translocation, which are important to the robustness of the network and contribute to the understanding of the broader spectrum of psychiatric symptoms. This study shows the importance of using an integrative approach to explore genotype-phenotype variability.</p></div

Summary of the characteristics observed in our patients, compared with those of other cases described in the literature.

<p>Summary of the characteristics observed in our patients, compared with those of other cases described in the literature.</p