Molecular restoration of archived transcriptional profiles by complementary-template reverse-transcription (CT-RT)

Gene expression profiling of formalin-fixed and paraffin-embedded (FFPE) specimens, banked from completed clinical trials and routine clinical care, has the potential to yield valuable information implicating and linking genes with clinical parameters. In order to prepare high-quality cDNA from highly fragmented FFPE-RNA, previously precluded from high-throughput analyses, we have designed a novel strategy based on the nucleic acid restoration of incomplete cDNA sequences prior to T7 in vitro transcription (IVT) amplification. We describe this strategy as complementary-template reverse-transcription (CT-RT) because short single-stranded T7-oligo-dT24-VN-DNA sequences, obtained from FFPE-RNA, are used as primers for the RT of complementary RNA templates contained in a sense-RNA library. We validated our assay by determining the correlation between expression profiles of a matched 10-year-old frozen and FFPE breast cancer sample. We show that T7 IVT-amplification of cDNA transcripts restored by CT-RT is a specific and reliable process that allows recovery of transcriptional features undetectable by direct T7 IVT-amplification of FFPE-RNA. Furthermore, CT-RT restored 35–41% of the transcripts from archived breast and cervical specimens when compared to matched frozen tissue; and profiles included tissue-specific transcripts. Our results indicate that CT-RT allows microarray profiling of severely degraded RNA that could not be analyzed by previous methods

Nonpromoter methylation of the CDKN2A gene with active transcription is associated with improved locoregional control in laryngeal squamous cell carcinoma

We previously reported a novel association between CDKN2A nonpromoter methylation and transcription (ARF/INK4a) in human papillomavirus associated oropharyngeal tumors. In this study we assessed whether nonpromoter CDKN2A methylation in laryngeal squamous cell carcinomas (LXSCC) conferred a similar association with transcription that predicted patient outcome. We compared DNA methylation and ARF/INK4a RNA expression levels for the CDKN2A locus using the Illumina HumanMethylation27 beadchip and RT-PCR in 43 LXSCC tumor samples collected from a prospective study of head and neck cancer patients treated at Montefiore Medical Center (MMC). Validation was performed using RNAseq data on 111 LXSCC tumor samples from the Cancer Genome Atlas (TCGA). The clinical relevance of combined nonpromoter CDKN2A methylation and transcription was assessed by multivariate Cox regression for locoregional recurrence on a subset of 69 LXSCC patients with complete clinicopathologic data from the MMC and TCGA cohorts. We found evidence of CDKN2A nonpromoter hypermethylation in a third of LXSCC from our MMC cohort, which was significantly associated with increased ARF and INK4a RNA expression (Wilcoxon rank-sum, P = 0.007 and 0.003, respectively). A similar association was confirmed in TCGA samples (Wilcoxon rank-sum test P < 0.0001 for ARF and INK4a). Patients with CDKN2A hypermethylation or high ARF/INK4a expression were significantly less likely to develop a locoregional recurrence compared to those with neither of the features, independent of other clinicopatholgic risk factors (adjusted hazard ratio=0.21, 95% confidence interval:0.05-0.81). These results support the conclusion that CDKN2A nonpromoter methylation is associated with increased ARF and INK4a RNA expression, and improved locoregional control in LXSCC

Genetic variants in glutamate, Aβ and tau related pathways determine polygenic risk for Alzheimer's disease

Synapse loss is an early event in late-onset Alzheimer's disease (LOAD). In this study we have assessed the capacity of a polygenic risk score (PRS) restricted to synapse-encoding loci to predict LOAD. We used summary statistics from the IGAP genome-wide association meta-analysis of 74,046 subjects for model construction and tested the "Synaptic PRS" in two independent datasets of controls and pathologically-confirmed LOAD. The mean Synaptic PRS was 2.3-fold higher in LOAD compared to controls (

Tumor Transcriptome Sequencing Reveals Allelic Expression Imbalances Associated with Copy Number Alterations

Due to growing throughput and shrinking cost, massively parallel sequencing is rapidly becoming an attractive alternative to microarrays for the genome-wide study of gene expression and copy number alterations in primary tumors. The sequencing of transcripts (RNA-Seq) should offer several advantages over microarray-based methods, including the ability to detect somatic mutations and accurately measure allele-specific expression. To investigate these advantages we have applied a novel, strand-specific RNA-Seq method to tumors and matched normal tissue from three patients with oral squamous cell carcinomas. Additionally, to better understand the genomic determinants of the gene expression changes observed, we have sequenced the tumor and normal genomes of one of these patients. We demonstrate here that our RNA-Seq method accurately measures allelic imbalance and that measurement on the genome-wide scale yields novel insights into cancer etiology. As expected, the set of genes differentially expressed in the tumors is enriched for cell adhesion and differentiation functions, but, unexpectedly, the set of allelically imbalanced genes is also enriched for these same cancer-related functions. By comparing the transcriptomic perturbations observed in one patient to his underlying normal and tumor genomes, we find that allelic imbalance in the tumor is associated with copy number mutations and that copy number mutations are, in turn, strongly associated with changes in transcript abundance. These results support a model in which allele-specific deletions and duplications drive allele-specific changes in gene expression in the developing tumor

Linking protective GAB2 variants, increased cortical GAB2 expression and decreased Alzheimer's Disease pathology

GRB-associated binding protein 2 (GAB2) represents a compelling genome-wide association signal for late-onset Alzheimer’s disease (LOAD) with reported odds ratios (ORs) ranging from 0.75–0.85. We tested eight GAB2 variants in four North American Caucasian case-control series (2,316 LOAD, 2,538 controls) for association with LOAD. Meta-analyses revealed ORs ranging from (0.61–1.20) with no significant association (all p>0.32). Four variants were hetergeneous across the populations (all p<0.02) due to a potentially inflated effect size (OR = 0.61–0.66) only observed in the smallest series (702 LOAD, 209 controls). Despite the lack of association in our series, the previously reported protective association for GAB2 remained after meta-analyses of our data with all available previously published series (11,952-22,253 samples; OR = 0.82–0.88; all p<0.04). Using a freely available database of lymphoblastoid cell lines we found that protective GAB2 variants were associated with increased GAB2 expression (p = 9.5×10−7−9.3×10−6). We next measured GAB2 mRNA levels in 249 brains and found that decreased neurofibrillary tangle (r = −0.34, p = 0.0006) and senile plaque counts (r = −0.32, p = 0.001) were both good predictors of increased GAB2 mRNA levels albeit that sex (r = −0.28, p = 0.005) may have been a contributing factor. In summary, we hypothesise that GAB2 variants that are protective against LOAD in some populations may act functionally to increase GAB2 mRNA levels (in lymphoblastoid cells) and that increased GAB2 mRNA levels are associated with significantly decreased LOAD pathology. These findings support the hypothesis that Gab2 may protect neurons against LOAD but due to significant population heterogeneity, it is still unclear whether this protection is detectable at the genetic level

Whole-genome sequencing reveals host factors underlying critical COVID-19

Critical COVID-19 is caused by immune-mediated inflammatory lung injury. Host genetic variation influences the development of illness requiring critical care1 or hospitalization2–4 after infection with SARS-CoV-2. The GenOMICC (Genetics of Mortality in Critical Care) study enables the comparison of genomes from individuals who are critically ill with those of population controls to find underlying disease mechanisms. Here we use whole-genome sequencing in 7,491 critically ill individuals compared with 48,400 controls to discover and replicate 23 independent variants that significantly predispose to critical COVID-19. We identify 16 new independent associations, including variants within genes that are involved in interferon signalling (IL10RB and PLSCR1), leucocyte differentiation (BCL11A) and blood-type antigen secretor status (FUT2). Using transcriptome-wide association and colocalization to infer the effect of gene expression on disease severity, we find evidence that implicates multiple genes—including reduced expression of a membrane flippase (ATP11A), and increased expression of a mucin (MUC1)—in critical disease. Mendelian randomization provides evidence in support of causal roles for myeloid cell adhesion molecules (SELE, ICAM5 and CD209) and the coagulation factor F8, all of which are potentially druggable targets. Our results are broadly consistent with a multi-component model of COVID-19 pathophysiology, in which at least two distinct mechanisms can predispose to life-threatening disease: failure to control viral replication; or an enhanced tendency towards pulmonary inflammation and intravascular coagulation. We show that comparison between cases of critical illness and population controls is highly efficient for the detection of therapeutically relevant mechanisms of disease

Whole-genome sequencing reveals host factors underlying critical COVID-19

Critical COVID-19 is caused by immune-mediated inflammatory lung injury. Host genetic variation influences the development of illness requiring critical care1 or hospitalization2,3,4 after infection with SARS-CoV-2. The GenOMICC (Genetics of Mortality in Critical Care) study enables the comparison of genomes from individuals who are critically ill with those of population controls to find underlying disease mechanisms. Here we use whole-genome sequencing in 7,491 critically ill individuals compared with 48,400 controls to discover and replicate 23 independent variants that significantly predispose to critical COVID-19. We identify 16 new independent associations, including variants within genes that are involved in interferon signalling (IL10RB and PLSCR1), leucocyte differentiation (BCL11A) and blood-type antigen secretor status (FUT2). Using transcriptome-wide association and colocalization to infer the effect of gene expression on disease severity, we find evidence that implicates multiple genes—including reduced expression of a membrane flippase (ATP11A), and increased expression of a mucin (MUC1)—in critical disease. Mendelian randomization provides evidence in support of causal roles for myeloid cell adhesion molecules (SELE, ICAM5 and CD209) and the coagulation factor F8, all of which are potentially druggable targets. Our results are broadly consistent with a multi-component model of COVID-19 pathophysiology, in which at least two distinct mechanisms can predispose to life-threatening disease: failure to control viral replication; or an enhanced tendency towards pulmonary inflammation and intravascular coagulation. We show that comparison between cases of critical illness and population controls is highly efficient for the detection of therapeutically relevant mechanisms of disease

Kruppel-family zinc finger proteins as emerging epigenetic biomarkers in head and neck squamous cell carcinoma

Abstract Background Krüppel-type zinc finger protein genes located on chromosome 19q13 are aberrantly hypermethylated with high frequency in all anatomic sub-sites of head and neck cancers as well as other epithelial tumours resulting in decreased expression. Methods We examined prognostic significance of ZNF154 and ZNF132 expression and DNA methylation in independent patient cohort of about 500 head and neck cancer patients in the Cancer Genome Atlas (TCGA). We also overexpressed these genes in HEK-293 cells, as well as the oral cancer cell line UM-SCC-1. Results In 20 patients from the TCGA cohort of HNSCC patients where ZNF154 and ZNF132 DNA methylation and RNA expression could be compared in tumor and adjacent normal tissue, there was increased DNA methylation and decreased expression of both ZNF154 and ZNF132 in primary tumours. Low ZNF154 and low ZNF132 expression were associated with shorter overall survival in both head and neck squamous cell carcinoma (HNSCC) and lung adenocarcinoma (LUAC patients). While expression of these proteins in HEK-293 cells produced full-length protein, only truncated copies could be expressed in head and neck cancer cells (UM-SCC-1). The truncated version of ZNF154 protein increased doubling time and reduced cell migration in UM-SCC-1 cancer cells. Conclusions Both ZNF132 and ZNF154 represent novel clinically significant biomarkers in head and neck cancer with potential tumour suppressive properties. Future studies will address the underlying molecular mechanisms by which ZNF154 expression in HNSCC contributes to the control of cell growth and migration. Graphical abstrac