Search

Omics data integration is already a reality. However, few omics-based algorithms show enough predictive ability to be implemented into clinics or public health domains. Clinical/epidemiological data tend to explain most of the variation of health-related traits, and its joint modeling with omics data is crucial to increase the algorithm’s predictive ability. Only a small number of published studies performed a “real” integration of omics and non-omics (OnO) data, mainly to predict cancer outcomes. Challenges in OnO data integration regard the nature and heterogeneity of non-omics data, the possibility of integrating large-scale non-omics data with high-throughput omics data, the relationship between OnO data (i.e., ascertainment bias), the presence of interactions, the fairness of the models, and the presence of subphenotypes. These challenges demand the development and application of new analysis strategies to integrate OnO data. In this contribution we discuss different attempts of OnO data integration in clinical and epidemiological studies. Most of the reviewed papers considered only one type of omics data set, mainly RNA expression data. All selected papers incorporated non-omics data in a low-dimensionality fashion. The integrative strategies used in the identified papers adopted three modeling methods: Independent, conditional, and joint modeling. This review presents, discusses, and proposes integrative analytical strategies towards OnO data integration, Fondo de Investigaciones Sanitarias (FIS), España, Instituto de Salud Carlos III, España, European Cooperation in Science and Technology, World Cancer Research, Ministerio de Economía y Competitividad, España, Depto. de Estadística y Ciencia de los Datos, Fac. de Estudios Estadísticos, TRUE, pub

ObjectiveChronic pancreatitis (CP) is a potentially fatal disease of the exocrine pancreas, with no specific or effective approved therapies. Due to difficulty in accessing pancreas tissues, little is known about local immune responses or pathogenesis in human CP. We sought to characterise pancreatic immune responses using tissues derived from patients with different aetiologies of CP and non-CP organ donors in order to identify key signalling molecules associated with human CP.DesignWe performed single-cell level cellular indexing of transcriptomes and epitopes by sequencing and T-cell receptor (TCR) sequencing of pancreatic immune cells isolated from organ donors, hereditary and idiopathic patients with CP who underwent total pancreatectomy. We validated gene expression data by performing flow cytometry and functional assays in a second patient with CP cohort.ResultsDeep single-cell sequencing revealed distinct immune characteristics and significantly enriched CCR6+ CD4+ T cells in hereditary compared with idiopathic CP. In hereditary CP, a reduction in T-cell clonality was observed due to the increased CD4+ T (Th) cells that replaced tissue-resident CD8+ T cells. Shared TCR clonotype analysis among T-cell lineages also unveiled unique interactions between CCR6+ Th and Th1 subsets, and TCR clustering analysis showed unique common antigen binding motifs in hereditary CP. In addition, we observed a significant upregulation of the CCR6 ligand (CCL20) expression among monocytes in hereditary CP as compared with those in idiopathic CP. The functional significance of CCR6 expression in CD4+ T cells was confirmed by flow cytometry and chemotaxis assay.ConclusionSingle-cell sequencing with pancreatic immune cells in human CP highlights pancreas-specific immune crosstalk through the CCR6-CCL20 axis, a signalling pathway that might be leveraged as a potential future target in human hereditary CP.

INTRODUCTION: Genomic alterations in the juxtamembrane exon 14 splice sites in NSCLC lead to increased MET stability and oncogenesis. We present the largest cohort study of MET Exon 14 (METex14) using whole transcriptome sequencing. METHODS: A total of 21,582 NSCLC tumor samples underwent complete genomic profiling with next-generation sequencing of DNA (592 Gene Panel, NextSeq, whole exome sequencing, NovaSeq) and RNA (NovaSeq, whole transcriptome sequencing). Clinicopathologic information including programmed death-ligand 1 and tumor mutational burden were collected and RNA expression for mutation subtypes and MET amplification were quantified. Immunogenic signatures and potential pathways of invasion were characterized using single-sample gene set enrichment analysis and mRNA gene signatures. RESULTS: A total of 533tumors (2.47%) with METex14 were identified. The most common alterations were point mutations (49.5%) at donor splice sites. Most alterations translated to increased MET expression, with MET co-amplification resulting in synergistic increase in expression (q < 0.05). Common coalterations were amplifications of MDM2 (19.0% versus 1.8% wild-type [WT]), HMGA2 (13.2% versus 0.98% WT), and CDK4 (10.0% versus 1.5% WT) (q < 0.05). High programmed death-ligand 1 > 50% (52.5% versus 27.3% WT, q < 0.0001) and lower proportion of high tumor mutational burden (>10 mutations per megabase, 8.3% versus 36.7% WT, p < 0.0001) were associated with METex14, which were also enriched in both immunogenic signatures and immunosuppressive checkpoints. Pathways associated with METex14 included angiogenesis and apical junction pathways (q < 0.05). CONCLUSIONS: METex14 splicing alterations and MET co-amplification translated to higher and synergistic MET expression at the transcriptomic level. High frequencies of MDM2 and CDK4 co-amplifications and association with multiple immunosuppressive checkpoints and angiogenic pathways provide insight into

The peripheral immune system response to Intracerebral Hemorrhage (ICH) may differ with ICH in different brain locations. Thus, we investigated peripheral blood mRNA expression of Deep ICH, Lobar ICH, and vascular risk factor-matched control subjects (n = 59). Deep ICH subjects usually had hypertension. Some Lobar ICH subjects had cerebral amyloid angiopathy (CAA). Genes and gene networks in Deep ICH and Lobar ICH were compared to controls. We found 774 differentially expressed genes (DEGs) and 2 co-expressed gene modules associated with Deep ICH, and 441 DEGs and 5 modules associated with Lobar ICH. Pathway enrichment showed some common immune/inflammatory responses between locations including Autophagy, T Cell Receptor, Inflammasome, and Neuroinflammation Signaling. Th2, Interferon, GP6, and BEX2 Signaling were unique to Deep ICH. Necroptosis Signaling, Protein Ubiquitination, Amyloid Processing, and various RNA Processing terms were unique to Lobar ICH. Finding amyloid processing pathways in blood of Lobar ICH patients suggests peripheral immune cells may participate in processes leading to perivascular/vascular amyloid in CAA vessels and/or are involved in its removal. This study identifies distinct peripheral blood transcriptome architectures in Deep and Lobar ICH, emphasizes the need for considering location in ICH studies/clinical trials, and presents potential location-specific treatment targets.

BackgroundLimited studies have functionally evaluated the heterogeneity in early ex vivo immune responses during sepsis. Our aim was to characterize early sepsis ex vivo functional immune response heterogeneity by studying whole blood endotoxin responses and derive a transcriptional metric of ex vivo endotoxin response.MethodsBlood collected within 24 h of hospital presentation from 40 septic patients was divided into two fractions and incubated with media (unstimulated) or endotoxin. Supernatants and cells were isolated, and responses measured using: supernatant cytokines, lung endothelial permeability after supernatant exposure, and RNA expression. A transcriptomic signature was derived in unstimulated cells to predict the ex vivo endotoxin response. The signature was tested in a separate cohort of 191 septic patients to evaluate for association with clinical outcome. Plasma biomarkers were quantified to measure in vivo host inflammation.ResultsEx vivo response to endotoxin varied and was unrelated to immunosuppression, white blood cell count, or the causative pathogen. Thirty-five percent of patients demonstrated a minimal response to endotoxin, suggesting early immunosuppression. High ex vivo cytokine production by stimulated blood cells correlated with increased in vitro pulmonary endothelial cell permeability and was associated with attenuated in vivo host inflammation. A four-gene signature of endotoxin response detectable without the need for a functional assay was identified. When tested in a separate cohort of septic patients, its expression was inversely associated with hospital mortality.ConclusionsAn attenuated ex vivo endotoxin response in early sepsis is associated with greater host in vivo inflammation and a worse clinical outcome.

The latent HIV-1 reservoir is comprised of stably integrated and intact proviruses with limited to no viral transcription. It has been proposed that latent infection may be maintained by methylation of pro-viral DNA. Here, for the first time, we investigate the cytosine methylation of a replication competent provirus (AMBI-1) found in a T cell clone in a donor on antiretroviral therapy (ART). Methylation profiles of the AMBI-1 provirus were compared to other proviruses in the same donor and in samples from three other individuals on ART, including proviruses isolated from lymph node mononuclear cells (LNMCs) and peripheral blood mononuclear cells (PBMCs). We also evaluated the apparent methylation of cytosines outside of CpG (i.e., CpH) motifs. We found no evidence for methylation in AMBI-1 or any other provirus tested within the 5' LTR promoter. In contrast, CpG methylation was observed in the env-tat-rev overlapping reading frame. In addition, we found evidence for differential provirus methylation in cells isolated from LNMCs vs. PBMCs in some individuals, possibly from the expansion of infected cell clones. Finally, we determined that apparent low-level methylation of CpH cytosines is consistent with occasional bisulfite reaction failures. In conclusion, our data do not support the proposition that latent HIV infection is associated with methylation of the HIV 5' LTR promoter.

Background: DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n = 6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled. Results: We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10 kb up- and downstream regions adjacent to the gene body. Conclusion: Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues

Background: DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n=6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled. Results: We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10kb up- and downstream regions adjacent to the gene body. Conclusion: Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues such as the liver in our case.

Background: DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n=6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled. Results: We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10kb up- and downstream regions adjacent to the gene body. Conclusion: Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues such as the liver in our case.

Background: DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n=6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled. Results: We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10kb up- and downstream regions adjacent to the gene body. Conclusion: Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues such as the liver in our case.

The latent HIV-1 reservoir is comprised of stably integrated and intact proviruses with limited to no viral transcription. It has been proposed that latent infection may be maintained by methylation of pro-viral DNA. Here, for the first time, we investigate the cytosine methylation of a replication competent provirus (AMBI-1) found in a T cell clone in a donor on antiretroviral therapy (ART). Methylation profiles of the AMBI-1 provirus were compared to other proviruses in the same donor and in samples from three other individuals on ART, including proviruses isolated from lymph node mononuclear cells (LNMCs) and peripheral blood mononuclear cells (PBMCs). We also evaluated the apparent methylation of cytosines outside of CpG (i.e., CpH) motifs. We found no evidence for methylation in AMBI-1 or any other provirus tested within the 5' LTR promoter. In contrast, CpG methylation was observed in the env-tat-rev overlapping reading frame. In addition, we found evidence for differential provirus methylation in cells isolated from LNMCs vs. PBMCs in some individuals, possibly from the expansion of infected cell clones. Finally, we determined that apparent low-level methylation of CpH cytosines is consistent with occasional bisulfite reaction failures. In conclusion, our data do not support the proposition that latent HIV infection is associated with methylation of the HIV 5' LTR promoter.

Renal oncocytoma (RO) accounts for 5% of renal cancers and generally behaves as a benign tumor with favorable long-term prognosis. It is difficult to confidently distinguish between benign RO and other renal malignancies, particularly chromophobe renal cell carcinoma (chRCC). Therefore, RO is often managed aggressively with surgery. We sought to identify molecular biomarkers to distinguish RO from chRCC and other malignant renal cancer mimics. In a 44-patient discovery cohort, we identified a significant differential abundance of nine genes in RO relative to chRCC. These genes were used to train a classifier to distinguish RO from chRCC in an independent 57-patient cohort. The trained classifier was then validated in five independent cohorts comprising 89 total patients. This nine-gene classifier trained on the basis of differential gene expression showed 93% sensitivity and 98% specificity for distinguishing RO from chRCC across the pooled validation cohorts, with a c-statistic of 0.978. This tool may be a useful adjunct to other diagnostic modalities to decrease the diagnostic and management uncertainty associated with small renal masses and to enable clinicians to recommend more confidently less aggressive management for some tumors. PATIENT SUMMARY: Renal oncocytoma is generally a benign form of kidney cancer that does not necessarily require surgical removal. However, it is difficult to distinguish renal oncocytoma from other more aggressive forms of kidney cancer, so it is treated most commonly with surgery. We built a classification tool based on the RNA levels of nine genes that may help avoid these surgeries by reliably distinguishing renal oncocytoma from other forms of kidney cancer.

Background: DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n=6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled. Results: We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10kb up- and downstream regions adjacent to the gene body. Conclusion: Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues such as the liver in our case.

Background: DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n = 6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled. Results: We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10 kb up- and downstream regions adjacent to the gene body. Conclusion: Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues

Background: DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n=6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled. Results: We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10kb up- and downstream regions adjacent to the gene body. Conclusion: Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues such as the liver in our case.

BackgroundDNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n = 6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled.ResultsWe simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10 kb up- and downstream regions adjacent to the gene body.ConclusionTemporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues such as the liver in our case. Howe

The level of response (LR) to alcohol as measured with the Self-Report of the Effects of Alcohol Retrospective Questionnaire (SRE) evaluates the number of standard drinks usually required for up to four effects. The need for a higher number of drinks for effects is genetically influenced and predicts higher risks for heavy drinking and alcohol problems. We conducted genome-wide association study (GWAS) in the African-American (COGA-AA, N = 1527 from 309 families) and European-American (COGA-EA, N = 4723 from 956 families) subsamples of the Collaborative Studies on the Genetics of Alcoholism (COGA) for two SRE scores: SRE-T (average of first five times of drinking, the period of heaviest drinking, and the most recent 3 months of consumption) and SRE-5 (the first five times of drinking). We then meta-analyzed the two COGA subsamples (COGA-AA + EA). Both SRE-T and SRE-5 were modestly heritable (h2 : 21%-31%) and genetically correlated with alcohol dependence (AD) and DSM-IV AD criterion count (rg : 0.35-0.76). Genome-wide significant associations were observed (SRE-T: chromosomes 6, rs140154945, COGA-EA P = 3.30E-08 and 11, rs10647170, COGA-AA+EA P = 3.53E-09; SRE-5: chromosome13, rs4770359, COGA-AA P = 2.92E-08). Chromosome 11 was replicated in an EA dataset from the National Institute on Alcohol Abuse and Alcoholism intramural program. In silico functional analyses and RNA expression analyses suggest that the chromosome 6 locus is an eQTL for KIF25. Polygenic risk scores derived using the COGA SRE-T and SRE-5 GWAS predicted 0.47% to 2.48% of variances in AD and DSM-IV AD criterion count in independent datasets. This study highlights the genetic contribution of alcohol response phenotypes to the etiology of alcohol use disorders.

Conventional flow cytometry has been widely used for high-throughput single-cell gene expression analysis using specific antibody staining. However, this is limited by the availability of high-quality antibodies. We developed a novel flow cytometry RNA detection technique termed RCA-Flow for single-cell RNA expression analysis. We showed that it is able to analyze not only mRNAs but also microRNAs and circular RNAs that are otherwise difficult to analyze by other flow cytometry techniques. The versatility for high-throughput analysis of different types of RNA molecules makes our method possess great potential for both biomedical and clinical applications.

The level of response (LR) to alcohol as measured with the Self-Report of the Effects of Alcohol Retrospective Questionnaire (SRE) evaluates the number of standard drinks usually required for up to four effects. The need for a higher number of drinks for effects is genetically influenced and predicts higher risks for heavy drinking and alcohol problems. We conducted genome-wide association study (GWAS) in the African-American (COGA-AA, N = 1527 from 309 families) and European-American (COGA-EA, N = 4723 from 956 families) subsamples of the Collaborative Studies on the Genetics of Alcoholism (COGA) for two SRE scores: SRE-T (average of first five times of drinking, the period of heaviest drinking, and the most recent 3 months of consumption) and SRE-5 (the first five times of drinking). We then meta-analyzed the two COGA subsamples (COGA-AA + EA). Both SRE-T and SRE-5 were modestly heritable (h2 : 21%-31%) and genetically correlated with alcohol dependence (AD) and DSM-IV AD criterion count (rg : 0.35-0.76). Genome-wide significant associations were observed (SRE-T: chromosomes 6, rs140154945, COGA-EA P = 3.30E-08 and 11, rs10647170, COGA-AA+EA P = 3.53E-09; SRE-5: chromosome13, rs4770359, COGA-AA P = 2.92E-08). Chromosome 11 was replicated in an EA dataset from the National Institute on Alcohol Abuse and Alcoholism intramural program. In silico functional analyses and RNA expression analyses suggest that the chromosome 6 locus is an eQTL for KIF25. Polygenic risk scores derived using the COGA SRE-T and SRE-5 GWAS predicted 0.47% to 2.48% of variances in AD and DSM-IV AD criterion count in independent datasets. This study highlights the genetic contribution of alcohol response phenotypes to the etiology of alcohol use disorders.

The level of response (LR) to alcohol as measured with the Self-Report of the Effects of Alcohol Retrospective Questionnaire (SRE) evaluates the number of standard drinks usually required for up to four effects. The need for a higher number of drinks for effects is genetically influenced and predicts higher risks for heavy drinking and alcohol problems. We conducted genome-wide association study (GWAS) in the African-American (COGA-AA, N = 1527 from 309 families) and European-American (COGA-EA, N = 4723 from 956 families) subsamples of the Collaborative Studies on the Genetics of Alcoholism (COGA) for two SRE scores: SRE-T (average of first five times of drinking, the period of heaviest drinking, and the most recent 3 months of consumption) and SRE-5 (the first five times of drinking). We then meta-analyzed the two COGA subsamples (COGA-AA + EA). Both SRE-T and SRE-5 were modestly heritable (h2 : 21%-31%) and genetically correlated with alcohol dependence (AD) and DSM-IV AD criterion count (rg : 0.35-0.76). Genome-wide significant associations were observed (SRE-T: chromosomes 6, rs140154945, COGA-EA P = 3.30E-08 and 11, rs10647170, COGA-AA+EA P = 3.53E-09; SRE-5: chromosome13, rs4770359, COGA-AA P = 2.92E-08). Chromosome 11 was replicated in an EA dataset from the National Institute on Alcohol Abuse and Alcoholism intramural program. In silico functional analyses and RNA expression analyses suggest that the chromosome 6 locus is an eQTL for KIF25. Polygenic risk scores derived using the COGA SRE-T and SRE-5 GWAS predicted 0.47% to 2.48% of variances in AD and DSM-IV AD criterion count in independent datasets. This study highlights the genetic contribution of alcohol response phenotypes to the etiology of alcohol use disorders.

Major depressive disorder (MDD) is a complex illness caused by both genetic and environmental factors. Antidepressant resistance also has a genetic component. To date, however, very few genes have been identified for major depression or antidepressant resistance. In this study, we investigated whether outbred heterogeneous stock (HS) rats would be a suitable model to uncover the genetics of depression and its connection to antidepressant resistance. The Wistar Kyoto (WKY) rat, one of the eight founders of the HS, is a recognized animal model of juvenile depression and is resistant to fluoxetine antidepressant treatment. We therefore hypothesized that adolescent HS rats would exhibit variation in both despair-like behavior and response to fluoxetine treatment. We assessed heritability of despair-like behavior and response to sub-acute fluoxetine using a modified forced swim test (FST) in 4-week-old HS rats. We also tested whether blood transcript levels previously identified as depression biomarkers in adolescent human subjects are differentially expressed in HS rats with high vs. low FST immobility. We demonstrate heritability of despair-like behavior in 4-week-old HS rats and show that many HS rats are resistant to fluoxetine treatment. In addition, blood transcript levels of Amfr, Cdr2 and Kiaa1539, genes previously identified in human adolescents with MDD, are differentially expressed between HS rats with high vs. low immobility. These data demonstrate that FST despair-like behavior will be amenable to genetic fine-mapping in adolescent HS rats. The overlap between human and HS blood biomarkers suggest that these studies may translate to depression in humans.

Serotonin is a major functional modulator within the brain and affects various behaviours as well as physiological functions. Particularly, throughout development, serotonin interacts with environmental challenges to shape later-life brain function and behaviour. In the thesis “Brain serotonin throughout development-for better and for worse” the role of brain serotonin was interrogated in rodent models, using various molecular tools and behavioural screenings. Results showed that dependent on sex, early-life experiences and experimental context, brain serotonin affected socio-emotional behaviours. Furthermore, alterations in brain serotonin function in interaction with early adversity influenced RNA expression and epigenetic marks in brain regions of emotion and stress regulation.

Major depressive disorder (MDD) is a complex illness caused by both genetic and environmental factors. Antidepressant resistance also has a genetic component. To date, however, very few genes have been identified for major depression or antidepressant resistance. In this study, we investigated whether outbred heterogeneous stock (HS) rats would be a suitable model to uncover the genetics of depression and its connection to antidepressant resistance. The Wistar Kyoto (WKY) rat, one of the eight founders of the HS, is a recognized animal model of juvenile depression and is resistant to fluoxetine antidepressant treatment. We therefore hypothesized that adolescent HS rats would exhibit variation in both despair-like behavior and response to fluoxetine treatment. We assessed heritability of despair-like behavior and response to sub-acute fluoxetine using a modified forced swim test (FST) in 4-week-old HS rats. We also tested whether blood transcript levels previously identified as depression biomarkers in adolescent human subjects are differentially expressed in HS rats with high vs. low FST immobility. We demonstrate heritability of despair-like behavior in 4-week-old HS rats and show that many HS rats are resistant to fluoxetine treatment. In addition, blood transcript levels of Amfr, Cdr2 and Kiaa1539, genes previously identified in human adolescents with MDD, are differentially expressed between HS rats with high vs. low immobility. These data demonstrate that FST despair-like behavior will be amenable to genetic fine-mapping in adolescent HS rats. The overlap between human and HS blood biomarkers suggest that these studies may translate to depression in humans.

In this thesis, we used support vector machines (SVMs) to build a tissue-of-origin classifier of 17 cancer types. Our classifier, which uses RNA expression data from over 20000 genes, works with high accuracy on primary (97.6%), metastasis (91.9%) and cell-line samples (71.1%). With the goal of enabling cheaper diagnostics for the clinics, we performed feature selection through recursive feature elimination (RFE) and identified a gene signature of just 120 genes that maintains almost all of the predictive power. We explored how our model could achieve such great accuracy and found that it recognises characteristics from healthy tissues rather than cancer. In order to help disseminate our results among clinicians and basic researchers, we released our trained model and its code in the command-line tool TOPOS (Tissue-of-Origin Predictor of OncoSamples). To our knowledge, this is the first time that a metastasis classifier is developed based on RNAseq data, and we hope to pave the way for others to do the same in the future. (I would like to explain that the reason why I cannot upload the thesis to UPCommons is that we are waiting for peer review to publish a paper with the results. Once the paper has been published, I would be happy to upload it to UPCommons as well if the paper's guidelines allow us to do so.)

Serotonin is a major functional modulator within the brain and affects various behaviours as well as physiological functions. Particularly, throughout development, serotonin interacts with environmental challenges to shape later-life brain function and behaviour. In the thesis “Brain serotonin throughout development-for better and for worse” the role of brain serotonin was interrogated in rodent models, using various molecular tools and behavioural screenings. Results showed that dependent on sex, early-life experiences and experimental context, brain serotonin affected socio-emotional behaviours. Furthermore, alterations in brain serotonin function in interaction with early adversity influenced RNA expression and epigenetic marks in brain regions of emotion and stress regulation.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Acute ischemic stroke affects over 800,000 US adults annually, with hundreds of thousands more experiencing a transient ischemic attack. Emergent evaluation, prompt acute treatment, and identification of stroke or TIA (transient ischemic attack) etiology for specific secondary prevention are critical for decreasing further morbidity and mortality of cerebrovascular disease. The Biomarkers of Acute Stroke Etiology (BASE) study is a multicenter observational study to identify serum markers defining the etiology of acute ischemic stroke. Observational trial of patients presenting to the hospital within 24 h of stroke onset. Blood samples are collected at arrival, 24, and 48 h later, and RNA gene expression is utilized to identify stroke etiology marker candidates. The BASE study began January 2014. At the time of writing, there are 22 recruiting sites. Enrollment is ongoing, expected to hit 1000 patients by March 2017. The BASE study could potentially aid in focusing the initial diagnostic evaluation to determine stroke etiology, with more rapidly initiated targeted evaluations and secondary prevention strategies.Clinical Trial Registration Clinicaltrials.gov NCT02014896 https://clinicaltrials.gov/ct2/show/NCT02014896?term=biomarkers+of+acute+stroke+etiology&rank=1.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Acute ischemic stroke affects over 800,000 US adults annually, with hundreds of thousands more experiencing a transient ischemic attack. Emergent evaluation, prompt acute treatment, and identification of stroke or TIA (transient ischemic attack) etiology for specific secondary prevention are critical for decreasing further morbidity and mortality of cerebrovascular disease. The Biomarkers of Acute Stroke Etiology (BASE) study is a multicenter observational study to identify serum markers defining the etiology of acute ischemic stroke. Observational trial of patients presenting to the hospital within 24 h of stroke onset. Blood samples are collected at arrival, 24, and 48 h later, and RNA gene expression is utilized to identify stroke etiology marker candidates. The BASE study began January 2014. At the time of writing, there are 22 recruiting sites. Enrollment is ongoing, expected to hit 1000 patients by March 2017. The BASE study could potentially aid in focusing the initial diagnostic evaluation to determine stroke etiology, with more rapidly initiated targeted evaluations and secondary prevention strategies.Clinical Trial Registration Clinicaltrials.gov NCT02014896 https://clinicaltrials.gov/ct2/show/NCT02014896?term=biomarkers+of+acute+stroke+etiology&rank=1.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

We studied HERV expression in cell lines after hypoxia, mitogenic stimulation, and demethylation, to better understand if hypoxia may play a role in ERV activation also within the nervous system, as represented by neuroblastoma cell lines. The level of RNA of four human ERV groups (HERVs) (HERVE, I/T, H, and W), and three housekeeping genes, of different cell lines including A549, COS-1, Namalwa, RD-L and Vero-E6, as well as human neuroblastoma cell lines SH-SY5Y, SK-N-DZ, and SK-N-AS were studied using reverse transcription and real-time quantitative PCR (QPCR). During the course of recovery from hypoxia a pronounced and selective activation of RNA expression of HERVW-like sequences, but not of HERVE, I/T, H, and three housekeeping genes, was found in the neuroblastoma cell lines, most pronounced in SK-N-DZ. In the SK-N-DZ cell line, we also tested the expression of HERVs after chemical treatments. HERVW-like sequences were selectively upregulated by 5-azacytidine, a demethylating agent. Some HERVW loci seem especially responsive to hypoxia and demethylation. HERV expression in neuroblastoma cells is selectively and profoundly influenced by some physiological and chemical stimuli.

The gastrointestinal tract (GIT) is subdivided into different anatomical organs with many shared functions and characteristics, but also distinct differences. We have combined a genome-wide transcriptomics analysis with immunohistochemistry-based protein profiling to describe the gene and protein expression patterns that define the human GIT. RNA sequencing data derived from stomach, duodenum, jejunum/ileum and colon specimens were compared to gene expression levels in 23 other normal human tissues analysed with the same method. Protein profiling based on immunohistochemistry and tissue microarrays was used to sub-localize the corresponding proteins with GIT-specific expression into sub-cellular compartments and cell types. Approximately 75 % of all human protein-coding genes were expressed in at least one of the GIT tissues. Only 51 genes showed enriched expression in either one of the GIT tissues and an additional 83 genes were enriched in two or more GIT tissues. The list of GIT-enriched genes with validated protein expression patterns included various well-known but also previously uncharacterised or poorly studied genes. For instance, the colon-enriched expression of NXPE family member 1 (NXPE1) was established, while NLR family, pyrin domain-containing 6 (NLRP6) expression was primarily found in the human small intestine. We have applied a genome-wide analysis based on transcriptomics and antibody-based protein profiling to identify genes that are expressed in a specific manner within the human GIT. These genes and proteins constitute important starting points for an improved understanding of the normal function and the different states of disease associated with the GIT., QC 20150324

BACKGROUND: The gastrointestinal tract (GIT) is subdivided into different anatomical organs with many shared functions and characteristics, but also distinct differences. We have combined a genome-wide transcriptomics analysis with immunohistochemistry-based protein profiling to describe the gene and protein expression patterns that define the human GIT. METHODS: RNA sequencing data derived from stomach, duodenum, jejunum/ileum and colon specimens were compared to gene expression levels in 23 other normal human tissues analysed with the same method. Protein profiling based on immunohistochemistry and tissue microarrays was used to sub-localize the corresponding proteins with GIT-specific expression into sub-cellular compartments and cell types. RESULTS: Approximately 75% of all human protein-coding genes were expressed in at least one of the GIT tissues. Only 51 genes showed enriched expression in either one of the GIT tissues and an additional 83 genes were enriched in two or more GIT tissues. The list of GIT-enriched genes with validated protein expression patterns included various well-known but also previously uncharacterised or poorly studied genes. For instance, the colon-enriched expression of NXPE family member 1 (NXPE1) was established, while NLR family, pyrin domain-containing 6 (NLRP6) expression was primarily found in the human small intestine. CONCLUSIONS: We have applied a genome-wide analysis based on transcriptomics and antibody-based protein profiling to identify genes that are expressed in a specific manner within the human GIT. These genes and proteins constitute important starting points for an improved understanding of the normal function and the different states of disease associated with the GIT., First online: 01 May 2014

Current RNA expression profiling methods rely on enrichment steps for specific RNA classes, thereby not detecting all RNA species in an unperturbed manner. We report strand-specific RNAome sequencing that determines expression of small and large RNAs from rRNA-depleted total RNA in a single sequence run. Since current analysis pipelines cannot reliably analyze small and large RNAs simultaneously, we developed TRAP, Total Rna Analysis Pipeline, a robust interface that is also compatible with existing RNA sequencing protocols. RNAome sequencing quantitatively preserved all RNA classes, allowing cross-class comparisons that facilitates the identification of relationships between different RNA classes. We demonstrate the strength of RNAome sequencing in mouse embryonic stem cells treated with cisplatin. MicroRNA and mRNA expression in RNAome sequencing significantly correlated between replicates and was in concordance with both existing RNA sequencing methods and gene expression arrays generated from the same samples. Moreover, RNAome sequencing also detected additional RNA classes such as enhancer RNAs, anti-sense RNAs, novel RNA species and numerous differentially expressed RNAs undetectable by other methods. At the level of complete RNA classes, RNAome sequencing also identified a specific global repression of the microRNA and microRNA isoform classes after cisplatin treatment whereas all other classes such as mRNAs were unchanged. These characteristics of RNAome sequencing will significantly improve expression analysis as well as studies on RNA biology not covered by existing methods.

Standard RNA expression profiling methods rely on enrichment steps for specific RNA classes, thereby not detecting all RNA species. For example, small and large RNAs from the same sample cannot be sequenced in a single sequence run. We designed RNAome sequencing, which is a strand-specific method to determine the expression of small and large RNAs from ribosomal RNA-depleted total RNA in a single sequence run. RNAome sequencing quantitatively preserves all RNA classes. This characteristic allows comparisons between RNA classes, thereby facilitating relationships between different RNA classes. Here, we describe in detail the experimental procedure associated with RNAome sequencing published by Derks and colleagues in RNA Biology (2015) [1]. We also provide the R code for the developed Total Rna Analysis Pipeline (TRAP), an algorithm to analyze RNAome sequencing datasets (deposited at the Gene Expression Omnibus data repository, accession number GSE48084).

Background - Oxytocin is a neuropeptide that is involved in the regulation of mood, anxiety and social biology. Genetic variation in the oxytocin receptor gene (OXTR) has been implicated in anxiety, depression and related stress phenotypes. It is not yet known whether OXTR interacts with other risk factors such as early life trauma to heighten the severity of experienced anxiety and depression.Methods - In this study, we examined genotypes in 653 individuals and tested whether SNP variation in OXTR correlates with severity of features of self-reported experience on the Depression Anxiety and Stress Scale (DASS), and whether this correlation is enhanced when early life trauma is taken into account. We also assessed the effects of OXTR SNPs on RNA expression levels in two separate brain tissue cohorts totaling 365 samples.Results - A significant effect of OXTR genotype on DASS anxiety, stress and depression scores was found and ELS events, in combination with several different OXTR SNPs, were significantly associated with differences in DASS scores with one SNP (rs139832701) showing significant association or a trend towards association for all three measures. Several OXTR SNPs were correlated with alterations in OXTR RNA expression and rs3831817 replicated across both sets of tissues.Conclusions - These results support the hypothesis that the oxytocin system plays a role in the pathophysiology of mood and anxiety disorders.