Genome-wide methylation profiles reveal quantitative views of human aging rates.

The ability to measure human aging from molecular profiles has practical implications in many fields, including disease prevention and treatment, forensics, and extension of life. Although chronological age has been linked to changes in DNA methylation, the methylome has not yet been used to measure and compare human aging rates. Here, we build a quantitative model of aging using measurements at more than 450,000 CpG markers from the whole blood of 656 human individuals, aged 19 to 101. This model measures the rate at which an individual's methylome ages, which we show is impacted by gender and genetic variants. We also show that differences in aging rates help explain epigenetic drift and are reflected in the transcriptome. Moreover, we show how our aging model is upheld in other human tissues and reveals an advanced aging rate in tumor tissue. Our model highlights specific components of the aging process and provides a quantitative readout for studying the role of methylation in age-related disease.

COL18A1 is a candidate eye iridocorneal angle-closure gene in humans.

Primary angle-closure glaucoma (PACG) is a common form of glaucoma in the Far East. Its defining feature is iridocorneal angle closure. In addition to PACG, indications of angle closure are included in the diagnostic criteria of related conditions primary angle-closure suspect (PACS) and primary angle closure (PAC). To the best of our knowledge, a causative gene for iridocorneal angle closure in humans has not been identified. This study aimed to identify the genetic cause of iridocorneal angle closure in a pedigree with at least 10 individuals diagnosed with PACS, PAC or PACG. Results of linkage analysis, segregation analysis of 44 novel variations, whole exome sequencing of 10 individuals, screenings of controls and bioinformatics predictions identified a mutation in COL18A1 that encodes collagen type XVIII as the most likely cause of angle closure in the pedigree. The role of COL18A1 in the etiology of Knobloch syndrome (KS) that is consistently accompanied by optic anomalies, available functional data on the encoded protein and the recognized role of collagens and the extracellular matrix in glaucoma pathogenesis supported the proposed role of the COL18A1 mutation in the pedigree. Subsequent identification of other COL18A1 mutations in PACS affected individuals of two unrelated families further supported that COL18A1 may affect angle closure. These PACS individuals were parents and grandparents of KS-affected children. In conclusion, a gene that affects angle closure in humans, a critical feature of PACG, has been identified. The findings also reinforce the importance of collagens in eye features and functions.

Pisces: an accurate and versatile variant caller for somatic and germline next-generation sequencing data.

MOTIVATION: Next-generation sequencing technology is transitioning quickly from research labs to clinical settings. The diagnosis and treatment selection for many acquired and autosomal conditions necessitate a method for accurately detecting somatic and germline variants. RESULTS: We have developed Pisces, a rapid, versatile and accurate small-variant calling suite designed for somatic and germline amplicon sequencing applications. Accuracy is achieved by four distinct modules, each incorporating a number of novel algorithmic strategies. AVAILABILITY AND IMPLEMENTATION: Pisces is distributed under an open source license and can be downloaded from https://github.com/Illumina/Pisces. Pisces is available on the BaseSpace™ SequenceHub. It is distributed on Illumina sequencing platforms such as the MiSeq™ and is included in the Praxis™ Extended RAS Panel test which was recently approved by the FDA. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DNA methylation profiles in African American prostate cancer patients in relation to disease progression.

This study examined whether differential DNA methylation is associated with clinical features of more aggressive disease at diagnosis and prostate cancer recurrence in African American men, who are more likely to die from prostate cancer than other populations. Tumor tissues from 76 African Americans diagnosed with prostate cancer who had radical prostatectomy as their primary treatment were profiled for epigenome-wide DNA methylation levels. Long-term follow-up identified 19 patients with prostate cancer recurrence. Twenty-three CpGs were differentially methylated (FDR q≤0.25, mean methylation difference≥0.10) in patients with vs. without recurrence, including CpGs in GCK, CDKL2, PRDM13, and ZFR2. Methylation differences were also observed between men with metastatic-lethal prostate cancer vs. no recurrence (five CpGs), regional vs. local pathological stage (two CpGs), and higher vs. lower tumor aggressiveness (one CpG). These results indicate that differentially methylated CpG sites identified in tumor tissues of African American men may contribute to prostate cancer aggressiveness.

Substantial DNA methylation differences between two major neuronal subtypes in human brain.

The brain is built from a large number of cell types which have been historically classified using location, morphology and molecular markers. Recent research suggests an important role of epigenetics in shaping and maintaining cell identity in the brain. To elucidate the role of DNA methylation in neuronal differentiation, we developed a new protocol for separation of nuclei from the two major populations of human prefrontal cortex neurons--GABAergic interneurons and glutamatergic (GLU) projection neurons. Major differences between the neuronal subtypes were revealed in CpG, non-CpG and hydroxymethylation (hCpG). A dramatically greater number of undermethylated CpG sites in GLU versus GABA neurons were identified. These differences did not directly translate into differences in gene expression and did not stem from the differences in hCpG methylation, as more hCpG methylation was detected in GLU versus GABA neurons. Notably, a comparable number of undermethylated non-CpG sites were identified in GLU and GABA neurons, and non-CpG methylation was a better predictor of subtype-specific gene expression compared to CpG methylation. Regions that are differentially methylated in GABA and GLU neurons were significantly enriched for schizophrenia risk loci. Collectively, our findings suggest that functional differences between neuronal subtypes are linked to their epigenetic specification.

Prostate tumor DNA methylation is associated with cigarette smoking and adverse prostate cancer outcomes.

BACKGROUND: DNA methylation has been hypothesized as a mechanism for explaining the association between smoking and adverse prostate cancer (PCa) outcomes. This study was aimed at assessing whether smoking is associated with prostate tumor DNA methylation and whether these alterations may explain in part the association of smoking with PCa recurrence and mortality. METHODS: A total of 523 men had radical prostatectomy as their primary treatment, detailed smoking history data, long-term follow-up for PCa outcomes, and tumor tissue profiled for DNA methylation. Ninety percent of the men also had matched tumor gene expression data. A methylome-wide analysis was conducted to identify differentially methylated regions (DMRs) by smoking status. To select potential functionally relevant DMRs, their correlation with the messenger RNA (mRNA) expression of corresponding genes was evaluated. Finally, a smoking-related methylation score based on the top-ranked DMRs was created to assess its association with PCa outcomes. RESULTS: Forty DMRs were associated with smoking status, and 10 of these were strongly correlated with mRNA expression (aldehyde oxidase 1 [AOX1], claudin 5 [CLDN5], early B-cell factor 1 [EBF1], homeobox A7 [HOXA7], lectin galactoside-binding soluble 3 [LGALS3], microtubule-associated protein τ [MAPT], protocadherin γ A [PCDHGA]/protocadherin γ B [PCDHGB], paraoxonase 3 [PON3], synaptonemal complex protein 2 like [SYCP2L], and zinc finger and SCAN domain containing 12 [ZSCAN12]). Men who were in the highest tertile for the smoking-methylation score derived from these DMRs had a higher risk of recurrence (odds ratio [OR], 2.29; 95% confidence interval [CI], 1.42-3.72) and lethal disease (OR, 4.21; 95% CI, 1.65-11.78) in comparison with men in the lower 2 tertiles. CONCLUSIONS: This integrative molecular epidemiology study supports the hypothesis that smoking-associated tumor DNA methylation changes may explain at least part of the association between smoking and adverse PCa outcomes. Future studies are warranted to confirm these findings and understand the implications for improving patient outcomes. Cancer 2016;122:2168-77. © 2016 American Cancer Society.

Mutation in ADORA1 identified as likely cause of early-onset parkinsonism and cognitive dysfunction.

BACKGROUND: We aimed to identify the genetic cause of neurological disease in an Iranian family whose manifestations include symptoms of parkinsonism and cognitive dysfunction. METHODS: Clinical data on the patients were gathered by interviews with parents, neurological examinations, and laboratory tests. Genetic analysis was performed by genome-wide single-nucleotide polymorphism homozygosity mapping and exome sequencing. The effect of putative disease-causing mutation was assessed by immunocytochemistry on HEK293 cells and Western blotting on proteins extracted from HEK293 cells transfected with wild-type and mutated genes. RESULTS: Homozygosity mapping and exome sequencing led to identification of a mutation in ADORA1 that causes p.Gly279Ser in the encoded protein, adenosine A1 receptor (A1 R), as the probable cause of disease. The mutation segregated with disease status in the family, affects a highly conserved amino acid, and was absent in 700 controls. CONCLUSIONS: The known biological activities of A1 R in brain functions including its physical interaction with and inhibitory effect on dopamine receptor D1 provide supportive evidence that disruptions of A1 R may result in neurological dysfunction. Also, recent evidence on the related adenosine A2B receptor marks the domain in which the mutation is positioned as important for function. Finally, ADORA1 is located within the Parkinson's disease locus PARK16, which has been identified in several populations. ADORA1 may be the PD susceptibility gene within this locus. The molecular mechanism by which p.Gly279Ser disrupts A1 R function remains unknown, but a quantitative effect on interaction with the dopamine receptor was not shown. © 2016 International Parkinson and Movement Disorder Society.

Differences in DNA methylation between human neuronal and glial cells are concentrated in enhancers and non-CpG sites.

We applied Illumina Human Methylation450K array to perform a genomic-scale single-site resolution DNA methylation analysis in neuronal and nonneuronal (primarily glial) nuclei separated from the orbitofrontal cortex of postmortem human brain. The findings were validated using enhanced reduced representation bisulfite sequencing. We identified thousands of sites differentially methylated (DM) between neuronal and nonneuronal cells. The DM sites were depleted within CpG-island-containing promoters but enriched in predicted enhancers. Classification of the DM sites into those undermethylated in neurons (neuronal type) and those undermethylated in nonneuronal cells (glial type), combined with findings of others that methylation within control elements typically negatively correlates with gene expression, yielded large sets of predicted neuron-specific and non-neuron-specific genes. These sets of predicted genes were in excellent agreement with the available direct measurements of gene expression in human and mouse. We also found a distinct set of DNA methylation patterns that were unique for neuronal cells. In particular, neuronal-type differential methylation was overrepresented in CpG island shores, enriched within gene bodies but not in intergenic regions, and preferentially harbored binding motifs for a distinct set of transcription factors, including neuron-specific activity-dependent factors. Finally, non-CpG methylation was substantially more prevalent in neurons than in nonneuronal cells.

Whole-genome haplotyping by dilution, amplification, and sequencing.

Standard whole-genome genotyping technologies are unable to determine haplotypes. Here we describe a method for rapid and cost-effective long-range haplotyping. Genomic DNA is diluted and distributed into multiple aliquots such that each aliquot receives a fraction of a haploid copy. The DNA template in each aliquot is amplified by multiple displacement amplification, converted into barcoded sequencing libraries using Nextera technology, and sequenced in multiplexed pools. To assess the performance of our method, we combined two male genomic DNA samples at equal ratios, resulting in a sample with diploid X chromosomes with known haplotypes. Pools of the multiplexed sequencing libraries were subjected to targeted pull-down of a 1-Mb contiguous region of the X-chromosome Duchenne muscular dystrophy gene. We were able to phase the Duchenne muscular dystrophy region into two contiguous haplotype blocks with a mean length of 494 kb. The haplotypes showed 99% agreement with the consensus base calls made by sequencing the individual DNAs. We subsequently used the strategy to haplotype two human genomes. Standard genomic sequencing to identify all heterozygous SNPs in the sample was combined with dilution-amplification-based sequencing data to resolve the phase of identified heterozygous SNPs. Using this procedure, we were able to phase >95% of the heterozygous SNPs from the diploid sequence data. The N50 for a Yoruba male DNA was 702 kb whereas the N50 for a European female DNA was 358 kb. Therefore, the strategy described here is suitable for haplotyping of a set of targeted regions as well as of the entire genome.

Epigenome-wide ovarian cancer analysis identifies a methylation profile differentiating clear-cell histology with epigenetic silencing of the HERG K+ channel.

Ovarian cancer remains the leading cause of death in women with gynecologic malignancies, despite surgical advances and the development of more effective chemotherapeutics. As increasing evidence indicates that clear-cell ovarian cancer may have unique pathogenesis, further understanding of molecular features may enable us to begin to understand the underlying biology and histology-specific information for improved outcomes. To study epigenetics in clear-cell ovarian cancer, fresh frozen tumor DNA (n = 485) was assayed on Illumina Infinium HumanMethylation450 BeadChips. We identified a clear-cell ovarian cancer tumor methylation profile (n = 163) which we validated in two independent replication sets (set 1, n = 163; set 2, n = 159), highlighting 22 CpG loci associated with nine genes (VWA1, FOXP1, FGFRL1, LINC00340, KCNH2, ANK1, ATXN2, NDRG21 and SLC16A11). Nearly all of the differentially methylated CpGs showed a propensity toward hypermethylation among clear-cell cases. Several loci methylation inversely correlated with tumor gene expression, most notably KCNH2 (HERG, a potassium channel) (P = 9.5 × 10(-7)), indicating epigenetic silencing. In addition, a predicted methylation class mainly represented by the clear-cell cases (20 clear cell out of 23 cases) had improved survival time. Although these analyses included only 30 clear-cell carcinomas, results suggest that loss of expression of KCNH2 (HERG) by methylation could be a good prognostic marker, given that overexpression of the potassium (K(+)) channel Eag family members promotes increased proliferation and results in poor prognosis. Validation in a bigger cohort of clear-cell tumors of the ovary is warranted.

Epigenomic profiling of DNA methylation in paired prostate cancer versus adjacent benign tissue.

BACKGROUND: Aberrant DNA methylation may promote prostate carcinogenesis. We investigated epigenome-wide DNA methylation profiles in prostate cancer (PCa) compared to adjacent benign tissue to identify differentially methylated CpG sites. METHODS: The study included paired PCa and adjacent benign tissue samples from 20 radical prostatectomy patients. Epigenetic profiling was done using the Infinium HumanMethylation450 BeadChip. Linear models that accounted for the paired study design and False Discovery Rate Q-values were used to evaluate differential CpG methylation. mRNA expression levels of the genes with the most differentially methylated CpG sites were analyzed. RESULTS: In total, 2,040 differentially methylated CpG sites were identified in PCa versus adjacent benign tissue (Q-value < 0.001), the majority of which were hypermethylated (n = 1,946; 95%). DNA methylation profiles accurately distinguished between PCa and benign tissue samples. Twenty-seven top-ranked hypermethylated CpGs had a mean methylation difference of at least 40% between tissue types, which included 25 CpGs in 17 genes. Furthermore, for 10 genes over 50% of promoter region CpGs were hypermethylated in PCa versus benign tissue. The top-ranked differentially methylated genes included three genes that were associated with both promoter hypermethylation and reduced gene expression: SCGB3A1, HIF3A, and AOX1. Analysis of The Cancer Genome Atlas (TCGA) data provided confirmatory evidence for our findings. CONCLUSIONS: This study of PCa versus adjacent benign tissue showed many differentially methylated CpGs and regions in and outside gene promoter regions, which may potentially be used for the development of future epigenetic-based diagnostic tests or as therapeutic targets.

Application of a low cost array-based technique - TAB-Array - for quantifying and mapping both 5mC and 5hmC at single base resolution in human pluripotent stem cells.

5-hydroxymethylcytosine (5hmC), an oxidized derivative of 5-methylcytosine (5mC), has been implicated as an important epigenetic regulator of mammalian development. Current procedures use DNA sequencing methods to discriminate 5hmC from 5mC, limiting their accessibility to the scientific community. Here we report a method that combines TET-assisted bisulfite conversion with Illumina 450K DNA methylation arrays for a low-cost high-throughput approach that distinguishes 5hmC and 5mC signals at base resolution. Implementing this approach, termed "TAB-array", we assessed DNA methylation dynamics in the differentiation of human pluripotent stem cells into cardiovascular progenitors and neural precursor cells. With the ability to discriminate 5mC and 5hmC, we identified a large number of novel dynamically methylated genomic regions that are implicated in the development of these lineages. The increased resolution and accuracy afforded by this approach provides a powerful means to investigate the distinct contributions of 5mC and 5hmC in human development and disease.

Integrative molecular analysis of intrahepatic cholangiocarcinoma reveals 2 classes that have different outcomes.

BACKGROUND & AIMS: Cholangiocarcinoma, the second most common liver cancer, can be classified as intrahepatic cholangiocarcinoma (ICC) or extrahepatic cholangiocarcinoma. We performed an integrative genomic analysis of ICC samples from a large series of patients. METHODS: We performed a gene expression profile, high-density single-nucleotide polymorphism array, and mutation analyses using formalin-fixed ICC samples from 149 patients. Associations with clinicopathologic traits and patient outcomes were examined for 119 cases. Class discovery was based on a non-negative matrix factorization algorithm and significant copy number variations were identified by Genomic Identification of Significant Targets in Cancer (GISTIC) analysis. Gene set enrichment analysis was used to identify signaling pathways activated in specific molecular classes of tumors, and to analyze their genomic overlap with hepatocellular carcinoma (HCC). RESULTS: We identified 2 main biological classes of ICC. The inflammation class (38% of ICCs) is characterized by activation of inflammatory signaling pathways, overexpression of cytokines, and STAT3 activation. The proliferation class (62%) is characterized by activation of oncogenic signaling pathways (including RAS, mitogen-activated protein kinase, and MET), DNA amplifications at 11q13.2, deletions at 14q22.1, mutations in KRAS and BRAF, and gene expression signatures previously associated with poor outcomes for patients with HCC. Copy number variation-based clustering was able to refine these molecular groups further. We identified high-level amplifications in 5 regions, including 1p13 (9%) and 11q13.2 (4%), and several focal deletions, such as 9p21.3 (18%) and 14q22.1 (12% in coding regions for the SAV1 tumor suppressor). In a complementary approach, we identified a gene expression signature that was associated with reduced survival times of patients with ICC; this signature was enriched in the proliferation class (P < .001). CONCLUSIONS: We used an integrative genomic analysis to identify 2 classes of ICC. The proliferation class has specific copy number alterations, activation of oncogenic pathways, and is associated with worse outcome. Different classes of ICC, based on molecular features, therefore might require different treatment approaches.

The novel mutation p.Asp251Asn in the ß-subunit of succinate-CoA ligase causes encephalomyopathy and elevated succinylcarnitine.

SUCLA2 is one of several nuclear-encoded genes that can cause encephalomyopathy accompanied by mitochondrial DNA depletion. The disorder usually manifests in early childhood and leads to early death. The gene encodes one of the subunits of succinyl-CoA synthase, the enzyme that catalyzes the reversible conversion of substrates succinyl-CoA and ADP to products succinate and ATP in the tricarboxylic acid pathway. Thirty-two individuals harboring mutations in SUCLA2 have so far been reported, and five different mutations were observed among these individuals. Here we report identification of a novel mutation in SUCLA2 in two cousins affected with encephalomyopathy. The novel mutation causes p.Asp251Asn; the affected amino acid is likely positioned within the ATP-grasp domain of the encoded protein. As previously reported in other patients, we did not observe elevation of methylmalonic acid, the biochemical hallmark of patients with mutations in SUCLA2. We instead found elevated levels of succinylcarnitine.

PANK2 and C19orf12 mutations are common causes of neurodegeneration with brain iron accumulation.

BACKGROUND: Neurodegeneration with brain iron accumulation (NBIA) constitutes a group of neurodegenerative disorders with pronounced iron deposition in the basal ganglia. PANK2 mutations are the most common cause of these disorders. C19orf12 was recently reported as another causative gene. We present phenotypic data and results of screening of PANK2 and C19orf12 in 11 unrelated Iranian NBIA patients. METHODS: Phenotypic data were obtained by neurologic examination, magnetic resonance imaging, and interviews. Mutation screening of PANK2 and C19orf12 was performed by sequencing. RESULTS: PANK2 and C19orf12 mutations were found in 7 and 4 patients, respectively. Phenotypic comparisons suggest that C19orf12 mutations as compared with PANK2 mutations result in a milder disease course. CONCLUSIONS: Mutations in both PANK2 and C19orf12 contributed significantly to NBIA in the Iranian patients. To the best of our knowledge, this is the first genetic analysis reported on a cohort of NBIA patients from the Middle East.

LTBP2 mutations cause Weill-Marchesani and Weill-Marchesani-like syndrome and affect disruptions in the extracellular matrix.

Latent transforming growth factor (TGF) beta-binding protein 2 (LTBP2) is an extracellular matrix (ECM) protein that associates with fibrillin-1 containing microfibrils. Various factors prompted considering LTBP2 in the etiology of isolated ectopia lentis and associated conditions such as Weill-Marchesani syndrome (WMS) and Marfan syndrome (MFS). LTBP2 was screened in 30 unrelated Iranian patients. Mutations were found only in one WMS proband and one MFS proband. Homozygous c.3529G>A (p.Val1177Met) was shown to cause autosomal recessive WMS or WM-like syndrome by several approaches, including homozygosity mapping. Light, fluorescent, and electron microscopy evidenced disruptions of the microfibrillar network in the ECM of the proband's skin. In conjunction with recent findings regarding other ECM proteins, the results presented strongly support the contention that anomalies in WMS patients are due to disruptions in the ECM. Heterozygous c.1642C >T (p.Arg548*) possibly contributed to MFS-related phenotypes, including ocular manifestations, mitral valve prolapse, and pectus excavatum, but was not cause of MFS.

Genetic and genomic analyses of RNA polymerase II-pausing factor in regulation of mammalian transcription and cell growth.

Many mammalian genes are occupied by paused RNA polymerase II (pol II) in the promoter-proximal region on both sides of the transcription start site. However, the impact of pol II pausing on gene expression and cell biology is not fully understood. In this study, we used a Cre-Lox system to conditionally knock out the b subunit of mouse negative elongation factor (Nelf-b), a key pol II-pausing factor, in mouse embryonic fibroblasts. We found that Nelf-b was associated with the promoter-proximal region of the majority of expressed genes, yet genetic ablation of Nelf-b only affected the steady-state mRNA levels of a small percentage of the Nelf-b-associated genes. Interestingly, Nelf-b deletion also increased levels of transcription start site upstream transcripts at multiple negative elongation factor-associated genes. The direct target genes of Nelf-b were highly enriched with those involved in the control of cell growth and cell death. Correspondingly, Nelf-b knock-out mouse embryonic fibroblasts exhibited slower progression from quiescence to proliferation, as well as in a cycling cell population. Furthermore, Nelf-b deletion also resulted in increased apoptosis. Thus, the genetic and genomic studies provide new physiological and molecular insight into Nelf-mediated pol II pausing.

High density DNA methylation array with single CpG site resolution.

We have developed a new generation of genome-wide DNA methylation BeadChip which allows high-throughput methylation profiling of the human genome. The new high density BeadChip can assay over 480K CpG sites and analyze twelve samples in parallel. The innovative content includes coverage of 99% of RefSeq genes with multiple probes per gene, 96% of CpG islands from the UCSC database, CpG island shores and additional content selected from whole-genome bisulfite sequencing data and input from DNA methylation experts. The well-characterized Infinium® Assay is used for analysis of CpG methylation using bisulfite-converted genomic DNA. We applied this technology to analyze DNA methylation in normal and tumor DNA samples and compared results with whole-genome bisulfite sequencing (WGBS) data obtained for the same samples. Highly comparable DNA methylation profiles were generated by the array and sequencing methods (average R2 of 0.95). The ability to determine genome-wide methylation patterns will rapidly advance methylation research.

Gene expression profiling of human whole blood samples with the Illumina WG-DASL assay.

BACKGROUND: Microarray-based gene expression analysis of peripheral whole blood is a common strategy in the development of clinically relevant biomarker panels for a variety of human diseases. However, the results of such an analysis are often plagued by decreased sensitivity and reliability due to the effects of relatively high levels of globin mRNA in whole blood. Globin reduction assays have been shown to overcome such effects, but they require large amounts of total RNA and may induce distinct gene expression profiles. The Illumina whole genome DASL assay can detect gene expression levels using partially degraded RNA samples and has the potential to detect rare transcripts present in highly heterogeneous whole blood samples without the need for globin reduction. We assessed the utility of the whole genome DASL assay in an analysis of peripheral whole blood gene expression profiles. RESULTS: We find that gene expression detection is significantly increased with the use of whole genome DASL compared to the standard IVT-based direct hybridization. Additionally, globin-probe negative whole genome DASL did not exhibit significant improvements over globin-probe positive whole genome DASL. Globin reduction further increases the detection sensitivity and reliability of both whole genome DASL and IVT-based direct hybridization with little effect on raw intensity correlations. Raw intensity correlations between total RNA and globin reduced RNA were 0.955 for IVT-based direct hybridization and 0.979 for whole genome DASL. CONCLUSIONS: Overall, the detection sensitivity of the whole genome DASL assay is higher than the IVT-based direct hybridization assay, with or without globin reduction, and should be considered in conjunction with globin reduction methods for future blood-based gene expression studies.

Analysis of gene expression in stage I serous tumors identifies critical pathways altered in ovarian cancer.

OBJECTIVE: Despite recent advances in the conceptual understanding of the pathogenesis of ovarian cancer, it remains the foremost cause of death from gynecologic malignancies in developed countries. The main reason for such a high rate of mortality is the lack of sensitive and specific biomarkers and imaging techniques for early detection of ovarian cancer. Additional biological insights into early-stage ovarian carcinogenesis are needed to help speed the development of markers for early detection of ovarian cancer. The objective of this study was to characterize differentially expressed genes in high-grade stage I serous carcinoma of the ovary. METHODS: We analyzed gene expression in macrodissected formalin-fixed, paraffin-embedded samples from 5 high-grade stage I serous carcinomas and 5 stage I borderline tumors of the ovary using the Illumina Whole Genome DASL assay (cDNA-mediated annealing, selection, extension, and ligation) corresponding to 24,000 genes. Significance Analysis of Microarrays was performed to determine differentially expressed genes in stage I serous carcinoma, and class prediction analysis was performed to determine the predictive value of differentially expressed gene sets to correctly classify serous carcinoma from borderline tumors in 3 independent data sets. Altered transcription factor pathways and biological pathways unique to stage I serous carcinoma were identified through class comparison of differentially expressed genes. RESULTS: Unsupervised cluster analysis of gene expression correctly classified stage I serous carcinomas from serous borderline tumors. Supervised analysis identified several known, as well as novel, genes differentially expressed in stage I ovarian cancer. Using a differentially expressed gene set, class comparison prediction analysis correctly identified serous carcinomas from serous borderline tumors in 3 independent data sets at over 80% accuracy, sensitivity, and specificity. Pathway analysis demonstrated the significance of p53 and E2F pathways in serous carcinogenesis and significant involvements of cell cycle and immune response pathways in stage I serous epithelial ovarian cancer. CONCLUSION: We have identified differentially expressed genes associated with the carcinogenesis of high-grade stage I serous EOC. Furthermore, integrative analysis of biological and transcription pathway data contributed to the confirmation of important biological pathways and discovery of additional unique genes and pathways that may have potential importance in ovarian pathogenesis and biomarker development.