Progranulin AAV gene therapy for frontotemporal dementia: translational studies and phase 1/2 trial interim results.

GRN mutations cause progranulin haploinsufficiency, which eventually leads to frontotemporal dementia (FTD-GRN). PR006 is an investigational gene therapy delivering the granulin gene (GRN) using an adeno-associated virus serotype 9 (AAV9) vector. In non-clinical studies, PR006 transduced neurons derived from induced pluripotent stem cells of patients with FTD-GRN, resulted in progranulin expression and improvement of lipofuscin, lysosomal and neuroinflammation pathologies in Grn-knockout mice, and was well tolerated except for minimal, asymptomatic dorsal root ganglionopathy in non-human primates. We initiated a first-in-human phase 1/2 open-label trial. Here we report results of a pre-specified interim analysis triggered with the last treated patient of the low-dose cohort (n = 6) reaching the 12-month follow-up timepoint. We also include preliminary data from the mid-dose cohort (n = 7). Primary endpoints were safety, immunogenicity and change in progranulin levels in cerebrospinal fluid (CSF) and blood. Secondary endpoints were Clinical Dementia Rating (CDR) plus National Alzheimer's Disease Coordinating Center (NACC) Frontotemporal Lobar Degeneration (FTLD) rating scale and levels of neurofilament light chain (NfL). One-time administration of PR006 into the cisterna magna was generally safe and well tolerated. All patients developed treatment-emergent anti-AAV9 antibodies in the CSF, but none developed anti-progranulin antibodies. CSF pleocytosis was the most common PR006-related adverse event. Twelve serious adverse events occurred, mostly unrelated to PR006. Deep vein thrombosis developed in three patients. There was one death (unrelated) occurring 18 months after treatment. CSF progranulin increased after PR006 treatment in all patients; blood progranulin increased in most patients but only transiently. NfL levels transiently increased after PR006 treatment, likely reflecting dorsal root ganglia toxicity. Progression rates, based on the CDR scale, were within the broad ranges reported for patients with FTD. These data provide preliminary insights into the safety and bioactivity of PR006. Longer follow-up and additional studies are needed to confirm the safety and potential efficacy of PR006. ClinicalTrials.gov identifier: NCT04408625 .

Exome-Based Rare-Variant Analyses in CKD.

BACKGROUND: Studies have identified many common genetic associations that influence renal function and all-cause CKD, but these explain only a small fraction of variance in these traits. The contribution of rare variants has not been systematically examined. METHODS: We performed exome sequencing of 3150 individuals, who collectively encompassed diverse CKD subtypes, and 9563 controls. To detect causal genes and evaluate the contribution of rare variants we used collapsing analysis, in which we compared the proportion of cases and controls carrying rare variants per gene. RESULTS: The analyses captured five established monogenic causes of CKD: variants in PKD1, PKD2, and COL4A5 achieved study-wide significance, and we observed suggestive case enrichment for COL4A4 and COL4A3. Beyond known disease-associated genes, collapsing analyses incorporating regional variant intolerance identified suggestive dominant signals in CPT2 and several other candidate genes. Biallelic mutations in CPT2 cause carnitine palmitoyltransferase II deficiency, sometimes associated with rhabdomyolysis and acute renal injury. Genetic modifier analysis among cases with APOL1 risk genotypes identified a suggestive signal in AHDC1, implicated in Xia-Gibbs syndrome, which involves intellectual disability and other features. On the basis of the observed distribution of rare variants, we estimate that a two- to three-fold larger cohort would provide 80% power to implicate new genes for all-cause CKD. CONCLUSIONS: This study demonstrates that rare-variant collapsing analyses can validate known genes and identify candidate genes and modifiers for kidney disease. In so doing, these findings provide a motivation for larger-scale investigation of rare-variant risk contributions across major clinical CKD categories.

A case-control collapsing analysis identifies retinal dystrophy genes associated with ophthalmic disease in patients with no pathogenic ABCA4 variants.

PURPOSE: Variants in the ABCA4 gene are causal for a variety of retinal dystrophy phenotypes, including Stargardt disease (STGD1). However, 15% of patients who present with symptoms compatible with STGD1/ABCA4 disease do not have identifiable causal ABCA4 variants. We hypothesized that a case-control collapsing analysis in ABCA4-negative patients with compatible symptoms would provide an objective measure to identify additional disease genes. METHODS: We performed a genome-wide enrichment analysis of "qualifying variants"-ultrarare variants predicted to impact protein function-in protein-coding genes in 79 unrelated cases and 9028 unrelated controls. RESULTS: Despite modest sample size, two known retinal dystrophy genes, PRPH2 and CRX, achieved study-wide significance (p < 1.33 × 10-6) under a dominant disease model, and eight additional known retinal dystrophy genes achieved nominal significance (p < 0.05). Across these ten genes, the excess of qualifying variants explained up to 36.8% of affected individuals. Furthermore, under a recessive model, the cone-rod dystrophy gene CERKL approached study-wide significance. CONCLUSION: Our results indicate that case-control collapsing analyses can efficiently identify pathogenic variants in genes in non-ABCA4 retinal dystrophies. The genome-wide collapsing analysis framework is an objective discovery method particularly suitable in settings with overlapping disease phenotypes.

Whole-exome sequencing in the evaluation of fetal structural anomalies: a prospective cohort study.

BACKGROUND: Identification of chromosomal aneuploidies and copy number variants that are associated with fetal structural anomalies has substantial value. Although whole-exome sequencing (WES) has been applied to case series of a few selected prenatal cases, its value in routine clinical settings has not been prospectively assessed in a large unselected cohort of fetuses with structural anomalies. We therefore aimed to determine the incremental diagnostic yield (ie, the added value) of WES following uninformative results of standard investigations with karyotype testing and chromosomal microarray in an unselected cohort of sequential pregnancies showing fetal structural anomalies. METHODS: In this prospective cohort study, the parents of fetuses who were found to have a structural anomaly in a prenatal ultrasound were screened for possible participation in the study. These participants were predominantly identified in or were referred to the Columbia University Carmen and John Thain Center for Prenatal Pediatrics (New York, NY, USA). Fetuses with confirmed aneuploidy or a causal pathogenic copy number variant were excluded from WES analyses. By use of WES of the fetuses and parents (parent-fetus trios), we identified genetic variants that indicated an underlying cause (diagnostic genetic variants) and genetic variants that met the criteria of bioinformatic signatures that had previously been described to be significantly enriched among diagnostic genetic variants. FINDINGS: Between April 24, 2015, and April 19, 2017, 517 sequentially identified pregnant women found to have fetuses with a structural anomaly were screened for their eligibility for inclusion in our study. 71 (14%) couples declined testing, 87 (17%) trios were missing at least one DNA sample (from either parent or the fetus), 69 (13%) trios had a clinically relevant abnormal karyotype or chromosomal microarray finding, 51 (10%) couples did not consent to WES or withdrew consent, and five (1%) samples were not of good enough quality for analysis. DNA samples from 234 (45%) eligible trios were therefore used for analysis of the primary outcome. By use of trio sequence data, we identified diagnostic genetic variants in 24 (10%) families. Mutations with bioinformatic signatures that were indicative of pathogenicity but with insufficient evidence to be considered diagnostic were also evaluated; 46 (20%) of the 234 fetuses assessed were found to have such signatures. INTERPRETATION: Our analysis of WES data in a prospective cohort of unselected fetuses with structural anomalies shows the value added by WES following the use of routine genetic tests. Our findings suggest that, in cases of fetal anomalies in which assessment with karyotype testing and chromosomal microarray fail to determine the underlying cause of a structural anomaly, WES can add clinically relevant information that could assist current management of a pregnancy. The unique challenges of WES-based prenatal diagnostics require analysis by a multidisciplinary team of perinatal practitioners and laboratory specialists. FUNDING: Institute for Genomic Medicine (Columbia University Irving Medical Center).

Improved Pathogenic Variant Localization via a Hierarchical Model of Sub-regional Intolerance.

Different parts of a gene can be of differential importance to development and health. This regional heterogeneity is also apparent in the distribution of disease-associated mutations, which often cluster in particular regions of disease-associated genes. The ability to precisely estimate functionally important sub-regions of genes will be key in correctly deciphering relationships between genetic variation and disease. Previous methods have had some success using standing human variation to characterize this variability in importance by measuring sub-regional intolerance, i.e., the depletion in functional variation from expectation within a given region of a gene. However, the ability to precisely estimate local intolerance was restricted by the fact that only information within a given sub-region is used, leading to instability in local estimates, especially for small regions. We show that borrowing information across regions using a Bayesian hierarchical model stabilizes estimates, leading to lower variability and improved predictive utility. Specifically, our approach more effectively identifies regions enriched for ClinVar pathogenic variants. We also identify significant correlations between sub-region intolerance and the distribution of pathogenic variation in disease-associated genes, with AUCs for classifying de novo missense variants in Online Mendelian Inheritance in Man (OMIM) genes of up to 0.86 using exonic sub-regions and 0.91 using sub-regions defined by protein domains. This result immediately suggests that considering the intolerance of regions in which variants are found may improve diagnostic interpretation. We also illustrate the utility of integrating regional intolerance into gene-level disease association tests with a study of known disease-associated genes for epileptic encephalopathy.

Diagnostic Utility of Exome Sequencing for Kidney Disease.

BACKGROUND: Exome sequencing is emerging as a first-line diagnostic method in some clinical disciplines, but its usefulness has yet to be examined for most constitutional disorders in adults, including chronic kidney disease, which affects more than 1 in 10 persons globally. METHODS: We conducted exome sequencing and diagnostic analysis in two cohorts totaling 3315 patients with chronic kidney disease. We assessed the diagnostic yield and, among the patients for whom detailed clinical data were available, the clinical implications of diagnostic and other medically relevant findings. RESULTS: In all, 3037 patients (91.6%) were over 21 years of age, and 1179 (35.6%) were of self-identified non-European ancestry. We detected diagnostic variants in 307 of the 3315 patients (9.3%), encompassing 66 different monogenic disorders. Of the disorders detected, 39 (59%) were found in only a single patient. Diagnostic variants were detected across all clinically defined categories, including congenital or cystic renal disease (127 of 531 patients [23.9%]) and nephropathy of unknown origin (48 of 281 patients [17.1%]). Of the 2187 patients assessed, 34 (1.6%) had genetic findings for medically actionable disorders that, although unrelated to their nephropathy, would also lead to subspecialty referral and inform renal management. CONCLUSIONS: Exome sequencing in a combined cohort of more than 3000 patients with chronic kidney disease yielded a genetic diagnosis in just under 10% of cases. (Funded by the National Institutes of Health and others.).

Whole-exome sequencing in 20,197 persons for rare variants in Alzheimer's disease.

OBJECTIVE: The genetic bases of Alzheimer's disease remain uncertain. An international effort to fully articulate genetic risks and protective factors is underway with the hope of identifying potential therapeutic targets and preventive strategies. The goal here was to identify and characterize the frequency and impact of rare and ultra-rare variants in Alzheimer's disease, using whole-exome sequencing in 20,197 individuals. METHODS: We used a gene-based collapsing analysis of loss-of-function ultra-rare variants in a case-control study design with data from the Washington Heights-Inwood Columbia Aging Project, the Alzheimer's Disease Sequencing Project and unrelated individuals from the Institute of Genomic Medicine at Columbia University. RESULTS: We identified 19 cases carrying extremely rare SORL1 loss-of-function variants among a collection of 6,965 cases and a single loss-of-function variant among 13,252 controls (P = 2.17 × 10-8; OR: 36.2 [95% CI: 5.8-1493.0]). Age-at-onset was 7 years earlier for patients with SORL1 qualifying variant compared with noncarriers. No other gene attained a study-wide level of statistical significance, but multiple top-ranked genes, including GRID2IP,WDR76 and GRN, were among candidates for follow-up studies. INTERPRETATION: This study implicates ultra-rare, loss-of-function variants in SORL1 as a significant genetic risk factor for Alzheimer's disease and provides a comprehensive dataset comparing the burden of rare variation in nearly all human genes in Alzheimer's disease cases and controls. This is the first investigation to establish a genome-wide statistically significant association between multiple extremely rare loss-of-function variants in SORL1 and Alzheimer's disease in a large whole-exome study of unrelated cases and controls.

Genome-wide Analyses Identify KIF5A as a Novel ALS Gene.

To identify novel genes associated with ALS, we undertook two lines of investigation. We carried out a genome-wide association study comparing 20,806 ALS cases and 59,804 controls. Independently, we performed a rare variant burden analysis comparing 1,138 index familial ALS cases and 19,494 controls. Through both approaches, we identified kinesin family member 5A (KIF5A) as a novel gene associated with ALS. Interestingly, mutations predominantly in the N-terminal motor domain of KIF5A are causative for two neurodegenerative diseases: hereditary spastic paraplegia (SPG10) and Charcot-Marie-Tooth type 2 (CMT2). In contrast, ALS-associated mutations are primarily located at the C-terminal cargo-binding tail domain and patients harboring loss-of-function mutations displayed an extended survival relative to typical ALS cases. Taken together, these results broaden the phenotype spectrum resulting from mutations in KIF5A and strengthen the role of cytoskeletal defects in the pathogenesis of ALS.

A case-control collapsing analysis identifies epilepsy genes implicated in trio sequencing studies focused on de novo mutations.

Trio exome sequencing has been successful in identifying genes with de novo mutations (DNMs) causing epileptic encephalopathy (EE) and other neurodevelopmental disorders. Here, we evaluate how well a case-control collapsing analysis recovers genes causing dominant forms of EE originally implicated by DNM analysis. We performed a genome-wide search for an enrichment of "qualifying variants" in protein-coding genes in 488 unrelated cases compared to 12,151 unrelated controls. These "qualifying variants" were selected to be extremely rare variants predicted to functionally impact the protein to enrich for likely pathogenic variants. Despite modest sample size, three known EE genes (KCNT1, SCN2A, and STXBP1) achieved genome-wide significance (p<2.68×10-6). In addition, six of the 10 most significantly associated genes are known EE genes, and the majority of the known EE genes (17 out of 25) originally implicated in trio sequencing are nominally significant (p<0.05), a proportion significantly higher than the expected (Fisher's exact p = 2.33×10-17). Our results indicate that a case-control collapsing analysis can identify several of the EE genes originally implicated in trio sequencing studies, and clearly show that additional genes would be implicated with larger sample sizes. The case-control analysis not only makes discovery easier and more economical in early onset disorders, particularly when large cohorts are available, but also supports the use of this approach to identify genes in diseases that present later in life when parents are not readily available.

Exome-wide Association Study Identifies GREB1L Mutations in Congenital Kidney Malformations.

Renal agenesis and hypodysplasia (RHD) are major causes of pediatric chronic kidney disease and are highly genetically heterogeneous. We conducted whole-exome sequencing in 202 case subjects with RHD and identified diagnostic mutations in genes known to be associated with RHD in 7/202 case subjects. In an additional affected individual with RHD and a congenital heart defect, we found a homozygous loss-of-function (LOF) variant in SLIT3, recapitulating phenotypes reported with Slit3 inactivation in the mouse. To identify genes associated with RHD, we performed an exome-wide association study with 195 unresolved case subjects and 6,905 control subjects. The top signal resided in GREB1L, a gene implicated previously in Hoxb1 and Shha signaling in zebrafish. The significance of the association, which was p = 2.0 × 10-5 for novel LOF, increased to p = 4.1 × 10-6 for LOF and deleterious missense variants combined, and augmented further after accounting for segregation and de novo inheritance of rare variants (joint p = 2.3 × 10-7). Finally, CRISPR/Cas9 disruption or knockdown of greb1l in zebrafish caused specific pronephric defects, which were rescued by wild-type human GREB1L mRNA, but not mRNA containing alleles identified in case subjects. Together, our study provides insight into the genetic landscape of kidney malformations in humans, presents multiple candidates, and identifies SLIT3 and GREB1L as genes implicated in the pathogenesis of RHD.

Characterization of a variant of t(14;18) negative nodal diffuse follicular lymphoma with CD23 expression, 1p36/TNFRSF14 abnormalities, and STAT6 mutations.

A predominantly diffuse growth pattern and CD23 co-expression are uncommon findings in nodal follicular lymphoma and can create diagnostic challenges. A single case series in 2009 (Katzenberger et al) proposed a unique morphologic variant of nodal follicular lymphoma, characterized by a predominantly diffuse architecture, lack of the t(14;18) IGH/BCL2 translocation, presence of 1p36 deletion, frequent inguinal lymph node involvement, CD23 co-expression, and low clinical stage. Other studies on CD23+ follicular lymphoma, while associating inguinal location, have not specifically described this architecture. In addition, no follow-up studies have correlated the histopathologic and cytogenetic/molecular features of these cases, and they remain a diagnostic problem. We identified 11 cases of diffuse, CD23+ follicular lymphoma with histopathologic features similar to those described by Katzenberger et al. Along with pertinent clinical information, we detail their histopathology, IGH/BCL2 translocation status, lymphoma-associated chromosomal gains/losses, and assessment of mutations in 220 lymphoma-associated genes by massively parallel sequencing. All cases showed a diffuse growth pattern around well- to ill-defined residual germinal centers, uniform CD23 expression, mixed centrocytic/centroblastic cytology, and expression of at least one germinal center marker. Ten of 11 involved inguinal lymph nodes, 5 solely. By fluorescence in situ hybridization analysis, the vast majority lacked IGH/BCL2 translocation (9/11). Deletion of 1p36 was observed in five cases and included TNFRSF14. Of the six cases lacking 1p36 deletion, TNFRSF14 mutations were identified in three, highlighting the strong association of 1p36/TNFRSF14 abnormalities with this follicular lymphoma variant. In addition, 9 of the 11 cases tested (82%) had STAT6 mutations and nuclear P-STAT6 expression was detectable in the mutated cases by immunohistochemistry. The proportion of STAT6 mutations is higher than recently reported in conventional follicular lymphoma (11%). These findings lend support for a clinicopathologic variant of t(14;18) negative nodal follicular lymphoma and suggests importance of the interleukin (IL)-4/JAK/STAT6 pathway in this variant.

Brief-exposure to preoperative bevacizumab reveals a TGF-ß signature predictive of response in HER2-negative breast cancers.

To best define biomarkers of response, and to shed insight on mechanism of action of certain clinically important agents for early breast cancer, we used a brief-exposure paradigm in the preoperative setting to study transcriptional changes in patient tumors that occur with one dose of therapy prior to combination chemotherapy. Tumor biopsies from breast cancer patients enrolled in two preoperative clinical trials were obtained at baseline and after one dose of bevacizumab (HER2-negative), trastuzumab (HER2-positive) or nab-paclitaxel, followed by treatment with combination chemo-biologic therapy. RNA-Sequencing based PAM50 subtyping at baseline of 46 HER2-negative patients revealed a strong association between the basal-like subtype and pathologic complete response (pCR) to chemotherapy plus bevacizumab (p ≤ 0.0027), but did not provide sufficient specificity to predict response. However, a single dose of bevacizumab resulted in down-regulation of a well-characterized TGF-ß activity signature in every single breast tumor that achieved pCR (p ≤ 0.004). The TGF-ß signature was confirmed to be a tumor-specific read-out of the canonical TGF-ß pathway using pSMAD2 (p ≤ 0.04), with predictive power unique to brief-exposure to bevacizumab (p ≤ 0.016), but not trastuzumab or nab-paclitaxel. Down-regulation of TGF-ß activity was associated with reduction in tumor hypoxia by transcription and protein levels, suggesting therapy-induced disruption of an autocrine-loop between tumor stroma and malignant cells. Modulation of the TGF-ß pathway upon brief-exposure to bevacizumab may provide an early functional readout of pCR to preoperative anti-angiogenic therapy in HER2-negative breast cancer, thus providing additional avenues for exploration in both preclinical and clinical settings with these agents.

Identification of a novel clone, ST736, among Enterococcus faecium clinical isolates and its association with daptomycin nonsusceptibility.

Resistance to daptomycin in enterococcal clinical isolates remains rare but is being increasingly reported in the United States and worldwide. There are limited data on the genetic relatedness and microbiological and clinical characteristics of daptomycin-nonsusceptible enterococcal clinical isolates. In this study, we assessed the population genetics of daptomycin-nonsusceptible Enterococcus faecium (DNSE) clinical isolates by multilocus sequence typing (MLST) and whole-genome sequencing analysis. Forty-two nonduplicate DNSE isolates and 43 randomly selected daptomycin-susceptible E. faecium isolates were included in the analysis. All E. faecium isolates were recovered from patients at a tertiary care medical center in suburban New York City from May 2009 through December 2013. The daptomycin MICs of the DNSE isolates ranged from 6 to >256 µg/ml. Three major clones of E. faecium (ST18, ST412, and ST736) were identified among these clinical isolates by MLST and whole-genome sequence-based analysis. A newly recognized clone, ST736, was seen in 32 of 42 (76.2%) DNSE isolates and in only 14 of 43 (32.6%) daptomycin-susceptible E. faecium isolates (P < 0.0001). This report provides evidence of the association between E. faecium clone ST736 and daptomycin nonsusceptibility. The identification and potential spread of this novel E. faecium clone and its association with daptomycin nonsusceptibility constitute a challenge for patient management and infection control at our medical center.

Integrative prediction of gene function and platinum-free survival from genomic and epigenetic features in ovarian cancer.

The identification of genetic and epigenetic alterations from primary tumor cells has become a common method to discover genes critical to the development, progression, and therapeutic resistance of cancer. We seek to identify those genetic and epigenetic aberrations that have the most impact on gene function within the tumor. First, we perform a bioinformatics analysis of copy number variation (CNV) and DNA methylation covering the genetic landscape of ovarian cancer tumor cells. We were specifically interested in copy number variation as our base genomic property in the prediction of tumor suppressors and oncogenes in the altered ovarian tumor. We identify changes in DNA methylation and expression specifically for all amplified and deleted genes. We statistically define tumor suppressor and oncogenic gene function from integrative analysis of three modalities: copy number variation, DNA methylation, and gene expression. Our method (1) calculates the extent of genomic and epigenetic alterations of defined tumor suppressor and oncogenic features for the functional prediction of significant ovarian cancer gene candidates and (2) identifies the functional activity or inactivity of known tumor suppressors and oncogenes in ovarian cancer. We applied our protocol on 42 primary serous ovarian cancer samples using MOMA-ROMA representational array assays. Additionally, we provide the basis for incorporating epigenetic profiles of ovarian tumors for the purposes of platinum-free survival prediction in the context of TCGA data.

Translating next generation sequencing to practice: opportunities and necessary steps.

Next-generation sequencing (NGS) approaches for measuring RNA and DNA benefit from greatly increased sensitivity, dynamic range and detection of novel transcripts. These technologies are rapidly becoming the standard for molecular assays and represent huge potential value to the practice of oncology. However, many challenges exist in the transition of these technologies from research application to clinical practice. This review discusses the value of NGS in detecting mutations, copy number changes and RNA quantification and their applications in oncology, the challenges for adoption and the relevant steps that are needed for translating this potential to routine practice.

Loss of DOK2 induces carboplatin resistance in ovarian cancer via suppression of apoptosis.

OBJECTIVE: Ovarian cancers are highly heterogeneous and while chemotherapy is the preferred treatment many patients are intrinsically resistant or quickly develop resistance. Furthermore, all tumors that recur ultimately become resistant. Recent evidence suggests that epigenetic deregulation may be a key factor in the onset and maintenance of chemoresistance. We set out to identify epigenetically silenced genes that affect chemoresistance. METHODS: The epigenomes of a total of 45 ovarian samples were analyzed to identify epigenetically altered genes that segregate with platinum response, and further filtered with expression data to identify genes that were suppressed. A tissue culture carboplatin resistance screen was utilized to functionally validate this set of candidate platinum resistance genes. RESULTS: Our screen correctly identified 19 genes that when suppressed altered the chemoresistance of the cells in culture. Of the genes identified in the screen we further characterized one gene, docking protein 2 (DOK2), an adapter protein downstream of tyrosine kinase, to determine if we could elucidate the mechanism by which it increased resistance. The loss of DOK2 decreased the level of apoptosis in response to carboplatin. Furthermore, in cells with reduced DOK2, the level of anoikis was decreased. CONCLUSIONS: We have developed a screening methodology that analyzes the epigenome and informatically identifies candidate genes followed by in vitro culture screening of the candidate genes. To validate our screening methodology we further characterized one candidate gene, DOK2, and showed that loss of DOK2 induces chemotherapy resistance by decreasing the level of apoptosis in response to treatment.

Effective normalization for copy number variation detection from whole genome sequencing.

BACKGROUND: Whole genome sequencing enables a high resolution view of the human genome and provides unique insights into genome structure at an unprecedented scale. There have been a number of tools to infer copy number variation in the genome. These tools, while validated, also include a number of parameters that are configurable to genome data being analyzed. These algorithms allow for normalization to account for individual and population-specific effects on individual genome CNV estimates but the impact of these changes on the estimated CNVs is not well characterized. We evaluate in detail the effect of normalization methodologies in two CNV algorithms FREEC and CNV-seq using whole genome sequencing data from 8 individuals spanning four populations. METHODS: We apply FREEC and CNV-seq to a sequencing data set consisting of 8 genomes. We use multiple configurations corresponding to different read-count normalization methodologies in FREEC, and statistically characterize the concordance of the CNV calls between FREEC configurations and the analogous output from CNV-seq. The normalization methodologies evaluated in FREEC are: GC content, mappability and control genome. We further stratify the concordance analysis within genic, non-genic, and a collection of validated variant regions. RESULTS: The GC content normalization methodology generates the highest number of altered copy number regions. Both mappability and control genome normalization reduce the total number and length of copy number regions. Mappability normalization yields Jaccard indices in the 0.07 - 0.3 range, whereas using a control genome normalization yields Jaccard index values around 0.4 with normalization based on GC content. The most critical impact of using mappability as a normalization factor is substantial reduction of deletion CNV calls. The output of another method based on control genome normalization, CNV-seq, resulted in comparable CNV call profiles, and substantial agreement in variable gene and CNV region calls. CONCLUSIONS: Choice of read-count normalization methodology has a substantial effect on CNV calls and the use of genomic mappability or an appropriately chosen control genome can optimize the output of CNV analysis.

Major chromosomal breakpoint intervals in breast cancer co-localize with differentially methylated regions.

Solid tumors exhibit chromosomal rearrangements resulting in gain or loss of multiple chromosomal loci (copy number variation, or CNV), and translocations that occasionally result in the creation of novel chimeric genes. In the case of breast cancer, although most individual tumors each have unique CNV landscape, the breakpoints, as measured over large datasets, appear to be non-randomly distributed in the genome. Breakpoints show a significant regional concentration at genomic loci spanning perhaps several megabases. The proximal cause of these breakpoint concentrations is a subject of speculation, but is, as yet, largely unknown. To shed light on this issue, we have performed a bio-statistical analysis on our previously published data for a set of 119 breast tumors and normal controls (Wiedswang et al., 2003), where each sample has both high-resolution CNV and methylation data. The method examined the distribution of closeness of breakpoint regions with differentially methylated regions (DMR), coupled with additional genomic parameters, such as repeat elements and designated "fragile sites" in the reference genome. Through this analysis, we have identified a set of 93 regional loci called breakpoint enriched DMR (BEDMRs) characterized by altered DNA methylation in cancer compared to normal cells that are associated with frequent breakpoint concentrations within a distance of 1 Mb. BEDMR loci are further associated with local hypomethylation (66%), concentrations of the Alu SINE repeats within 3 Mb (35% of the cases), and tend to occur near a number of cancer related genes such as the protocadherins, AKT1, DUB3, GAB2. Furthermore, BEDMRs seem to deregulate members of the histone gene family and chromatin remodeling factors, e.g., JMJD1B, which might affect the chromatin structure and disrupt coordinate signaling and repair. From this analysis we propose that preference for chromosomal breakpoints is related to genome structure coupled with alterations in DNA methylation and hence, chromatin structure, associated with tumorigenesis.

Identification of tumor suppressors and oncogenes from genomic and epigenetic features in ovarian cancer.

The identification of genetic and epigenetic alterations from primary tumor cells has become a common method to identify genes critical to the development and progression of cancer. We seek to identify those genetic and epigenetic aberrations that have the most impact on gene function within the tumor. First, we perform a bioinformatic analysis of copy number variation (CNV) and DNA methylation covering the genetic landscape of ovarian cancer tumor cells. We separately examined CNV and DNA methylation for 42 primary serous ovarian cancer samples using MOMA-ROMA assays and 379 tumor samples analyzed by The Cancer Genome Atlas. We have identified 346 genes with significant deletions or amplifications among the tumor samples. Utilizing associated gene expression data we predict 156 genes with altered copy number and correlated changes in expression. Among these genes CCNE1, POP4, UQCRB, PHF20L1 and C19orf2 were identified within both data sets. We were specifically interested in copy number variation as our base genomic property in the prediction of tumor suppressors and oncogenes in the altered ovarian tumor. We therefore identify changes in DNA methylation and expression for all amplified and deleted genes. We statistically define tumor suppressor and oncogenic features for these modalities and perform a correlation analysis with expression. We predicted 611 potential oncogenes and tumor suppressors candidates by integrating these data types. Genes with a strong correlation for methylation dependent expression changes exhibited at varying copy number aberrations include CDCA8, ATAD2, CDKN2A, RAB25, AURKA, BOP1 and EIF2C3. We provide copy number variation and DNA methylation analysis for over 11,500 individual genes covering the genetic landscape of ovarian cancer tumors. We show the extent of genomic and epigenetic alterations for known tumor suppressors and oncogenes and also use these defined features to identify potential ovarian cancer gene candidates.