Implementation of Nanopore sequencing as a pragmatic workflow for copy number variant confirmation in the clinic.

BACKGROUND: Diagnosis of rare genetic diseases can be a long, expensive and complex process, involving an array of tests in the hope of obtaining an actionable result. Long-read sequencing platforms offer the opportunity to make definitive molecular diagnoses using a single assay capable of detecting variants, characterizing methylation patterns, resolving complex rearrangements, and assigning findings to long-range haplotypes. Here, we demonstrate the clinical utility of Nanopore long-read sequencing by validating a confirmatory test for copy number variants (CNVs) in neurodevelopmental disorders and illustrate the broader applications of this platform to assess genomic features with significant clinical implications. METHODS: We used adaptive sampling on the Oxford Nanopore platform to sequence 25 genomic DNA samples and 5 blood samples collected from patients with known or false-positive copy number changes originally detected using short-read sequencing. Across the 30 samples (a total of 50 with replicates), we assayed 35 known unique CNVs (a total of 55 with replicates) and one false-positive CNV, ranging in size from 40 kb to 155 Mb, and assessed the presence or absence of suspected CNVs using normalized read depth. RESULTS: Across 50 samples (including replicates) sequenced on individual MinION flow cells, we achieved an average on-target mean depth of 9.5X and an average on-target read length of 4805 bp. Using a custom read depth-based analysis, we successfully confirmed the presence of all 55 known CNVs (including replicates) and the absence of one false-positive CNV. Using the same CNV-targeted data, we compared genotypes of single nucleotide variant loci to verify that no sample mix-ups occurred between assays. For one case, we also used methylation detection and phasing to investigate the parental origin of a 15q11.2-q13 duplication with implications for clinical prognosis. CONCLUSIONS: We present an assay that efficiently targets genomic regions to confirm clinically relevant CNVs with a concordance rate of 100%. Furthermore, we demonstrate how integration of genotype, methylation, and phasing data from the Nanopore sequencing platform can potentially simplify and shorten the diagnostic odyssey.

Somatic alterations contributing to metastasis of a castration-resistant prostate cancer.

Metastatic castration-resistant prostate cancer (mCRPC) is a lethal disease, and molecular markers that differentiate indolent from aggressive subtypes are needed. We sequenced the exomes of five metastatic tumors and healthy kidney tissue from an index case with mCRPC to identify lesions associated with disease progression and metastasis. An Ashkenazi Jewish (AJ) germline founder mutation, del185AG in BRCA1, was observed and AJ ancestry was confirmed. Sixty-two somatic variants altered proteins in tumors, including cancer-associated genes, TMPRSS2-ERG, PBRM1, and TET2. The majority (n = 53) of somatic variants were present in all metastases and only a subset (n = 31) was observed in the primary tumor. Integrating tumor next-generation sequencing and DNA copy number showed somatic loss of BRCA1 and TMPRSS2-ERG. We sequenced 19 genes with deleterious mutations in the index case in additional mCRPC samples and detected a frameshift, two somatic missense alterations, tumor loss of heterozygosity, and combinations of germline missense SNPs in TET2. In summary, genetic analysis of metastases from an index case permitted us to infer a chronology for the clonal spread of disease based on sequential accrual of somatic lesions. The role of TET2 in mCRPC deserves additional analysis and may define a subset of metastatic disease.

The new sequencer on the block: comparison of Life Technology's Proton sequencer to an Illumina HiSeq for whole-exome sequencing.

We assessed the performance of the new Life Technologies Proton sequencer by comparing whole-exome sequence data in a Centre d'Etude du Polymorphisme Humain trio (family 1463) to the Illumina HiSeq instrument. To simulate a typical user's results, we utilized the standard capture, alignment and variant calling methods specific to each platform. We restricted data analysis to include the capture region common to both methods. The Proton produced high quality data at a comparable average depth and read length, and the Ion Reporter variant caller identified 96 % of single nucleotide polymorphisms (SNPs) detected by the HiSeq and GATK pipeline. However, only 40 % of small insertion and deletion variants (indels) were identified by both methods. Usage of the trio structure and segregation of platform-specific alleles supported this result. Further comparison of the trio data with Complete Genomics sequence data and Illumina SNP microarray genotypes documented high concordance and accurate SNP genotyping of both Proton and Illumina platforms. However, our study underscored the problem of accurate detection of indels for both the Proton and HiSeq platforms.

Integration of a Cross-Ancestry Polygenic Model With Clinical Risk Factors Improves Breast Cancer Risk Stratification.

PURPOSE: To develop and validate a cross-ancestry integrated risk score (caIRS) that combines a cross-ancestry polygenic risk score (caPRS) with a clinical estimator for breast cancer (BC) risk. We hypothesized that the caIRS is a better predictor of BC risk than clinical risk factors across diverse ancestry groups. METHODS: We used diverse retrospective cohort data with longitudinal follow-up to develop a caPRS and integrate it with the Tyrer-Cuzick (T-C) clinical model. We tested the association between the caIRS and BC risk in two validation cohorts including > 130,000 women. We compared model discrimination for 5-year and remaining lifetime BC risk between the caIRS and T-C and assessed how the caIRS would affect screening in the clinic. RESULTS: The caIRS outperformed T-C alone for all populations tested in both validation cohorts and contributed significantly to risk prediction beyond T-C. The area under the receiver operating characteristic curve improved from 0.57 to 0.65, and the odds ratio per standard deviation increased from 1.35 (95% CI, 1.27 to 1.43) to 1.79 (95% CI, 1.70 to 1.88) in validation cohort 1 with similar improvements observed in validation cohort 2. We observed the largest gain in positive predictive value using the caIRS in Black/African American women across both validation cohorts, with an approximately two-fold increase and an equivalent negative predictive value as the T-C. In a multivariate, age-adjusted logistic regression model including both caIRS and T-C, caIRS remained significant, indicating that caIRS provides information over T-C alone. CONCLUSION: Adding a caPRS to the T-C model improves BC risk stratification for women of multiple ancestries, which could have implications for screening recommendations and prevention.

Promoter variants in the MSMB gene associated with prostate cancer regulate MSMB/NCOA4 fusion transcripts.

Beta-microseminoprotein (MSP)/MSMB is an immunoglobulin superfamily protein synthesized by prostate epithelial cells and secreted into seminal plasma. Variants in the promoter of the MSMB gene have been associated with the risk of prostate cancer (PCa) in several independent genome-wide association studies. Both MSMB and an adjacent gene, NCOA4, are subjected to transcriptional control via androgen response elements. The gene product of NCOA4 interacts directly with the androgen receptor as a co-activator to enhance AR transcriptional activity. Here, we provide evidence for the expression of full-length MSMB-NCOA4 fusion transcripts regulated by the MSMB promoter. The predominant MSMB-NCOA4 transcript arises by fusion of the 5'UTR and exons 1-2 of the MSMB pre-mRNA, with exons 2-10 of the NCOA4 pre-mRNA, producing a stable fusion protein, comprising the essential domains of NCOA4. Analysis of the splice sites of this transcript shows an unusually strong splice acceptor at NCOA4 exon 2 and the presence of Alu repeats flanking the exons potentially involved in the splicing event. Transfection experiments using deletion clones of the promoter coupled with luciferase reporter assays define a core MSMB promoter element located between -27 and -236 of the gene, and a negative regulatory element immediately upstream of the start codon. Computational network analysis reveals that the MSMB gene is functionally connected to NCOA4 and the androgen receptor signaling pathway. The data provide an example of how GWAS-associated variants may have multiple genetic and epigenetic effects.

Whole-genome risk prediction of common diseases in human preimplantation embryos.

Preimplantation genetic testing (PGT) of in-vitro-fertilized embryos has been proposed as a method to reduce transmission of common disease; however, more comprehensive embryo genetic assessment, combining the effects of common variants and rare variants, remains unavailable. Here, we used a combination of molecular and statistical techniques to reliably infer inherited genome sequence in 110 embryos and model susceptibility across 12 common conditions. We observed a genotype accuracy of 99.0-99.4% at sites relevant to polygenic risk scoring in cases from day-5 embryo biopsies and 97.2-99.1% in cases from day-3 embryo biopsies. Combining rare variants with polygenic risk score (PRS) magnifies predicted differences across sibling embryos. For example, in a couple with a pathogenic BRCA1 variant, we predicted a 15-fold difference in odds ratio (OR) across siblings when combining versus a 4.5-fold or 3-fold difference with BRCA1 or PRS alone. Our findings may inform the discussion of utility and implementation of genome-based PGT in clinical practice.

Haplotype structure in Ashkenazi Jewish BRCA1 and BRCA2 mutation carriers.

Three founder mutations in BRCA1 and BRCA2 contribute to the risk of hereditary breast and ovarian cancer in Ashkenazi Jews (AJ). They are observed at increased frequency in the AJ compared to other BRCA mutations in Caucasian non-Jews (CNJ). Several authors have proposed that elevated allele frequencies in the surrounding genomic regions reflect adaptive or balancing selection. Such proposals predict long-range linkage disequilibrium (LD) resulting from a selective sweep, although genetic drift in a founder population may also act to create long-distance LD. To date, few studies have used the tools of statistical genomics to examine the likelihood of long-range LD at a deleterious locus in a population that faced a genetic bottleneck. We studied the genotypes of hundreds of women from a large international consortium of BRCA1 and BRCA2 mutation carriers and found that AJ women exhibited long-range haplotypes compared to CNJ women. More than 50% of the AJ chromosomes with the BRCA1 185delAG mutation share an identical 2.1 Mb haplotype and nearly 16% of AJ chromosomes carrying the BRCA2 6174delT mutation share a 1.4 Mb haplotype. Simulations based on the best inference of Ashkenazi population demography indicate that long-range haplotypes are expected in the context of a genome-wide survey. Our results are consistent with the hypothesis that a local bottleneck effect from population size constriction events could by chance have resulted in the large haplotype blocks observed at high frequency in the BRCA1 and BRCA2 regions of Ashkenazi Jews.

Association assessment of copy number polymorphism and risk of age-related macular degeneration.

PURPOSE: We previously identified a genetic copy number polymorphism (CNP147) that was statistically associated with age-related macular degeneration (AMD) and that resides downstream of the complement factor H (CFH) gene. Factor H protein is polymorphic at amino acid 402, in which the resulting histidine containing moiety has been established to impart significant risk of AMD. We present a method to precisely determine the exact copy number of CNP147 and examine in more detail the association with AMD. DESIGN: Case-control study. PARTICIPANTS: A total of 421 Age-Related Eye Disease Study (AREDS) subjects, of whom approximately 35% were diagnosed with neovascular disease, 19% were diagnosed with geographic atrophy, 16% were diagnosed with both, 30% were diagnosed with large drusen, and 215 were controls. METHODS: By using copy number assays available from Applied Biosystems Inc. (Carlsbad, CA), we examined 4 loci spanning CNP147 and neighboring CNP148 in an AREDS matched case-control sample set. We analyzed these data by copy number while controlling for 2 high-risk CFH variants, rs1061170 (Y402H) and rs1410996. We phased the high-risk CFH variants with CNP147 and analyzed haplotype frequencies in cases and controls. To further validate copy numbers, 6 Utah Centre D'etude du Polymorphism Humaine (CEPH) families were typed for CNP147, and the segregation was assessed. MAIN OUTCOME MEASURES: Increased or decreased risk of AMD from genetic loci. RESULTS: Having fewer than 2 copies of CNP147 was associated with an estimated 43% reduction in odds of having AMD in this sample set (adjusted odds ratio [OR] = 0.57, P=0.006). CNP148 variation is rare in Caucasians and was not statistically significant. Common haplotypes reveal that the risk alleles for rs1061170 and rs1410996 most frequently segregate with higher copy numbers for CNP147, but not exclusively, and that 1 haplotype that carried a deletion of CNP147 was highly protective (OR = 0.25 P=1.3×10(-13)) when compared with the reference. CONCLUSIONS: In this matched subset of AREDS subjects, after adjusting for 2 known risk variants in CFH, CNP147 deletion statistically associates with diminished risk for AMD. FINANCIAL DISCLOSURE(S): The author(s) have no proprietary or commercial interest in any materials discussed in this article.

Moving out: from sterol transport to drug resistance - the ABCG subfamily of efflux pumps.

The ATP binding cassette (ABC) proteins are typically ATP-driven transmembrane pumps that have been evolutionarily conserved from bacteria to humans. In humans these transporters are subdivided into seven subfamilies, ranging from A to G. The ABCG subfamily of transporters is the primary focus of this review. This subfamily of proteins has been conserved throughout evolution and plays a central role in several cellular processes, such as sterol homeostasis and multidrug resistance. Functional polymorphisms/mutations in some of these G-subfamily transporters have clinical consequences in humans.

Metastatic and recurrent adrenocortical cancer is not defined by its genomic landscape.

BACKGROUND: Adrenocortical carcinoma (ACC) is a rare, often-aggressive neoplasm of the adrenal cortex, with a 14-17 month median overall survival. We asked whether tumors from patients with advanced or metastatic ACC would offer clues as to putative genes that might have critical roles in disease progression or in more aggressive disease biology. METHODS: We conducted comprehensive genomic and expression analyses of ACCs from 43 patients, 30 female, and 42 from metastatic sites, including deep sequencing, copy number analysis, mRNA expression and microRNA arrays. RESULTS: Copy number gains and losses were similar to that previously reported for ACC. We identified a median mutation rate of 3.38 per megabase (Mb). The mutational signature was characterized by a predominance of C > T, C > A and T > C transitions. Only cancer genes TP53 (26%) and beta-catenin (CTNNB1, 14%) were mutated in more than 10% of samples. The TCGA-identified putative cancer genes MEN1 and PRKAR1A were found in low frequency-4.7 and 2.3%, respectively. The majority of the mutations were in genes not implicated in the etiology or maintenance of cancer. Specifically, amongst the 38 genes that were mutated in more than 9% of samples, only four were represented in Tier 1 of the 576 COSMIC Cancer Gene Census (CCGC). Thus, 82% of genes found to have mutations likely have no role in the etiology or biology of ACC; while the role of the other 18%, if any, remains to be proven. Finally, the transcript length for the 38 most frequently mutated genes in ACC is statistically longer than the average of all coding genes, raising the question of whether transcript length in part determined mutation probability. CONCLUSIONS: We conclude that the mutational and expression profiles of advanced and metastatic tumors are very similar to those from newly diagnosed patients-with very little in the way of genomic aberration to explain differences in biology. With relatively low mutation rates, few major oncogenic drivers, and loss of function mutations in several epigenetic regulators, an epigenetic basis for ACC may be postulated and serve as the basis for future studies.

Addressing health disparities in Hispanic breast cancer: accurate and inexpensive sequencing of BRCA1 and BRCA2.

BACKGROUND: Germline mutations in the BRCA1 and BRCA2 genes account for 20-25 % of inherited breast cancers and about 10 % of all breast cancer cases. Detection of BRCA mutation carriers can lead to therapeutic interventions such as mastectomy, oophorectomy, hormonal prevention therapy, improved screening, and targeted therapies such as PARP-inhibition. We estimate that African Americans and Hispanics are 4-5 times less likely to receive BRCA screening, despite having similar mutation frequencies as non-Jewish Caucasians, who have higher breast cancer mortality. To begin addressing this health disparity, we initiated a nationwide trial of BRCA testing of Latin American women with breast cancer. Patients were recruited through community organizations, clinics, public events, and by mail and Internet. Subjects completed the consent process and questionnaire, and provided a saliva sample by mail or in person. DNA from 120 subjects was used to sequence the entirety of BRCA1 and BRCA2 coding regions and splice sites, and validate pathogenic mutations, with a total material cost of $85/subject. Subjects ranged in age from 23 to 81 years (mean age, 51 years), 6 % had bilateral disease, 57 % were ER/PR+, 23 % HER2+, and 17 % had triple-negative disease. RESULTS: A total of seven different predicted deleterious mutations were identified, one newly described and the rest rare. In addition, four variants of unknown effect were found. CONCLUSIONS: Application of this strategy on a larger scale could lead to improved cancer care of minority and underserved populations.

Genome Analysis of Latin American Cervical Cancer: Frequent Activation of the PIK3CA Pathway.

PURPOSE: Cervical cancer is one of the most common causes of cancer mortality for women living in poverty, causing more than 28,000 deaths annually in Latin America and 266,000 worldwide. To better understand the molecular basis of the disease, we ascertained blood and tumor samples from Guatemala and Venezuela and performed genomic characterization. EXPERIMENTAL DESIGN: We performed human papillomavirus (HPV) typing and identified somatically mutated genes using exome and ultra-deep targeted sequencing with confirmation in samples from Mexico. Copy number changes were also assessed in the exome sequence. RESULTS: Cervical cancer cases in Guatemala and Venezuela have an average age of diagnosis of 50 years and 5.6 children. Analysis of 675 tumors revealed activation of PIK3CA and other PI3K/AKT pathway genes in 31% of squamous carcinomas and 24% of adeno- and adenosquamous tumors, predominantly at two sites (E542K, E545K) in the helical domain of the PIK3CA gene. This distribution of PIK3CA mutations is distinct from most other cancer types and does not result in the in vitro phosphorylation of AKT. Somatic mutations were more frequent in squamous carcinomas diagnosed after the age of 50 years. Frequent gain of chromosome 3q was found, and low PIK3CA mutation fractions in many tumors suggest that PI3K mutation can be a late event in tumor progression. CONCLUSIONS: PI3K pathway mutation is important to cervical carcinogenesis in Latin America. Therapeutic agents that directly target PI3K could play a role in the therapy of this common malignancy.

Concurrent alterations in TERT, KDM6A, and the BRCA pathway in bladder cancer.

PURPOSE: Genetic analysis of bladder cancer has revealed a number of frequently altered genes, including frequent alterations of the telomerase (TERT) gene promoter, although few altered genes have been functionally evaluated. Our objective is to characterize alterations observed by exome sequencing and sequencing of the TERT promoter, and to examine the functional relevance of histone lysine (K)-specific demethylase 6A (KDM6A/UTX), a frequently mutated histone demethylase, in bladder cancer. EXPERIMENTAL DESIGN: We analyzed bladder cancer samples from 54 U.S. patients by exome and targeted sequencing and confirmed somatic variants using normal tissue from the same patient. We examined the biologic function of KDM6A using in vivo and in vitro assays. RESULTS: We observed frequent somatic alterations in BRCA1 associated protein-1 (BAP1) in 15% of tumors, including deleterious alterations to the deubiquitinase active site and the nuclear localization signal. BAP1 mutations contribute to a high frequency of tumors with breast cancer (BRCA) DNA repair pathway alterations and were significantly associated with papillary histologic features in tumors. BAP1 and KDM6A mutations significantly co-occurred in tumors. Somatic variants altering the TERT promoter were found in 69% of tumors but were not correlated with alterations in other bladder cancer genes. We examined the function of KDM6A, altered in 24% of tumors, and show depletion in human bladder cancer cells, enhanced in vitro proliferation, in vivo tumor growth, and cell migration. CONCLUSIONS: This study is the first to identify frequent BAP1 and BRCA pathway alterations in bladder cancer, show TERT promoter alterations are independent of other bladder cancer gene alterations, and show KDM6A loss is a driver of the bladder cancer phenotype.

Differential gene and microRNA expression between etoposide resistant and etoposide sensitive MCF7 breast cancer cell lines.

In order to develop targeted strategies for combating drug resistance it is essential to understand it's basic molecular mechanisms. In an exploratory study we have found several possible indicators of etoposide resistance operating in MCF7VP cells, including up-regulation of ABC transporter genes, modulation of miRNA, and alteration in copy numbers of genes.

Single-cell sequencing analysis characterizes common and cell-lineage-specific mutations in a muscle-invasive bladder cancer.

BACKGROUND: Cancers arise through an evolutionary process in which cell populations are subjected to selection; however, to date, the process of bladder cancer, which is one of the most common cancers in the world, remains unknown at a single-cell level. RESULTS: We carried out single-cell exome sequencing of 66 individual tumor cells from a muscle-invasive bladder transitional cell carcinoma (TCC). Analyses of the somatic mutant allele frequency spectrum and clonal structure revealed that the tumor cells were derived from a single ancestral cell, but that subsequent evolution occurred, leading to two distinct tumor cell subpopulations. By analyzing recurrently mutant genes in an additional cohort of 99 TCC tumors, we identified genes that might play roles in the maintenance of the ancestral clone and in the muscle-invasive capability of subclones of this bladder cancer, respectively. CONCLUSIONS: This work provides a new approach of investigating the genetic details of bladder tumoral changes at the single-cell level and a new method for assessing bladder cancer evolution at a cell-population level.

Adeno-associated viral vector-mediated transgene expression is independent of DNA methylation in primate liver and skeletal muscle.

Recombinant adeno-associated viral (rAAV) vectors can support long-term transgene expression in quiescent tissues. Intramuscular (i.m.) administration of a single-stranded AAV vector (ssAAV) in the nonhuman primate (NHP) results in a peak protein level at 2-3 months, followed by a decrease over several months before reaching a steady-state. To investigate transgene expression and vector genome persistence, we previously demonstrated that rAAV vector genomes associate with histones and form a chromatin structure in NHP skeletal muscle more than one year after injection. In the mammalian nucleus, chromatin remodeling via epigenetic modifications plays key role in transcriptional regulation. Among those, CpG hyper-methylation of promoters is a known hallmark of gene silencing. To assess the involvement of DNA methylation on the transgene expression, we injected NHP via the i.m. or the intravenous (i.v.) route with a recombinant ssAAV2/1 vector. The expression cassette contains the transgene under the transcriptional control of the constitutive Rous Sarcoma Virus promoter (RSVp). Total DNA isolated from NHP muscle and liver biopsies from 1 to 37 months post-injection was treated with sodium bisulfite and subsequently analyzed by pyrosequencing. No significant CpG methylation of the RSVp was found in rAAV virions or in vector DNA isolated from NHP transduced tissues. Direct de novo DNA methylation appears not to be involved in repressing transgene expression in NHP after gene transfer mediated by ssAAV vectors. The study presented here examines host/vector interactions and the impact on transgene expression in a clinically relevant model.

Analysis of the ABCA4 gene by next-generation sequencing.

PURPOSE: To find all possible disease-associated variants in coding sequences of the ABCA4 gene in a large cohort of patients diagnosed with ABCA4-associated diseases. METHODS: One hundred sixty-eight patients who had been clinically diagnosed with Stargardt disease, cone-rod dystrophy, and other ABCA4-associated phenotypes were prescreened for mutations in ABCA4 with the ABCA4 microarray, resulting in finding 1 of 2 expected mutations in 111 patients and 0 of 2 mutations in 57 patients. The next-generation sequencing (NGS) strategy was applied to these patients to sequence the entire coding region and the splice sites of the ABCA4 gene. Identified new variants were confirmed or rejected by Sanger sequencing and analyzed for possible pathogenicity by in silico programs and, where possible, by segregation analyses. RESULTS: Sequencing was successful in 159 of 168 patients and identified the second disease-associated allele in 49 of 103 (~48%) of patients with one previously identified mutation. Among those with no mutations, both disease-associated alleles were detected in 4 of 56 patients, and one mutation was detected in 10 of 56 patients. The authors detected a total of 57 previously unknown, possibly pathogenic, variants: 29 missense, 4 nonsense, 9 small deletions and 15 splice-site-altering variants. Of these, 55 variants were deemed pathogenic by a combination of predictive methods and segregation analyses. CONCLUSIONS: Many mutations in the coding sequences of the ABCA4 gene are still unknown, and many possibly reside in noncoding regions of the ABCA4 locus. Although the ABCA4 array remains a good first-pass screening option, the NGS platform is a time- and cost-efficient tool for screening large cohorts.