Rapid and accurate remethylation of DNA in Dnmt3a-deficient hematopoietic cells with restoration of DNMT3A activity.

Here, we characterize the DNA methylation phenotypes of bone marrow cells from mice with hematopoietic deficiency of Dnmt3a or Dnmt3b (or both enzymes) or expressing the dominant-negative Dnmt3aR878H mutation [R882H in humans; the most common DNMT3A mutation found in acute myeloid leukemia (AML)]. Using these cells as substrates, we defined DNA remethylation after overexpressing wild-type (WT) DNMT3A1, DNMT3B1, DNMT3B3 (an inactive splice isoform of DNMT3B), or DNMT3L (a catalytically inactive "chaperone" for DNMT3A and DNMT3B in early embryogenesis). Overexpression of DNMT3A for 2 weeks reverses the hypomethylation phenotype of Dnmt3a-deficient cells or cells expressing the R878H mutation. Overexpression of DNMT3L (which is minimally expressed in AML cells) also corrects the hypomethylation phenotype of Dnmt3aR878H/+ marrow, probably by augmenting the activity of WT DNMT3A encoded by the residual WT allele. DNMT3L reactivation may represent a previously unidentified approach for restoring DNMT3A activity in hematopoietic cells with reduced DNMT3A function.

Ultra-Deep Sequencing Reveals the Mutational Landscape of Classical Hodgkin Lymphoma.

The malignant Hodgkin and Reed Sternberg (HRS) cells of classical Hodgkin lymphoma (cHL) are scarce in affected lymph nodes, creating a challenge to detect driver somatic mutations. As an alternative to cell purification techniques, we hypothesized that ultra-deep exome sequencing would allow genomic study of HRS cells, thereby streamlining analysis and avoiding technical pitfalls. To test this, 31 cHL tumor/normal pairs were exome sequenced to approximately 1,000× median depth of coverage. An orthogonal error-corrected sequencing approach verified >95% of the discovered mutations. We identified mutations in genes novel to cHL including: CDH5 and PCDH7, novel stop gain mutations in IL4R, and a novel pattern of recurrent mutations in pathways regulating Hippo signaling. As a further application of our exome sequencing, we attempted to identify expressed somatic single-nucleotide variants (SNV) in single-nuclei RNA sequencing (snRNA-seq) data generated from a patient in our cohort. Our snRNA analysis identified a clear cluster of cells containing a somatic SNV identified in our deep exome data. This cluster has differentially expressed genes that are consistent with genes known to be dysregulated in HRS cells (e.g., PIM1 and PIM3). The cluster also contains cells with an expanded B-cell clonotype further supporting a malignant phenotype. This study provides proof-of-principle that ultra-deep exome sequencing can be utilized to identify recurrent mutations in HRS cells and demonstrates the feasibility of snRNA-seq in the context of cHL. These studies provide the foundation for the further analysis of genomic variants in large cohorts of patients with cHL. SIGNIFICANCE: Our data demonstrate the utility of ultra-deep exome sequencing in uncovering somatic variants in Hodgkin lymphoma, creating new opportunities to define the genes that are recurrently mutated in this disease. We also show for the first time the successful application of snRNA-seq in Hodgkin lymphoma and describe the expression profile of a putative cluster of HRS cells in a single patient.

Genomic landscape of TP53-mutated myeloid malignancies.

TP53-mutated myeloid malignancies are associated with complex cytogenetics and extensive structural variants, which complicates detailed genomic analysis by conventional clinical techniques. We performed whole-genome sequencing (WGS) of 42 acute myeloid leukemia (AML)/myelodysplastic syndromes (MDS) cases with paired normal tissue to better characterize the genomic landscape of TP53-mutated AML/MDS. WGS accurately determines TP53 allele status, a key prognostic factor, resulting in the reclassification of 12% of cases from monoallelic to multihit. Although aneuploidy and chromothripsis are shared with most TP53-mutated cancers, the specific chromosome abnormalities are distinct to each cancer type, suggesting a dependence on the tissue of origin. ETV6 expression is reduced in nearly all cases of TP53-mutated AML/MDS, either through gene deletion or presumed epigenetic silencing. Within the AML cohort, mutations of NF1 are highly enriched, with deletions of 1 copy of NF1 present in 45% of cases and biallelic mutations in 17%. Telomere content is increased in TP53-mutated AMLs compared with other AML subtypes, and abnormal telomeric sequences were detected in the interstitial regions of chromosomes. These data highlight the unique features of TP53-mutated myeloid malignancies, including the high frequency of chromothripsis and structural variation, the frequent involvement of unique genes (including NF1 and ETV6) as cooperating events, and evidence for altered telomere maintenance.

Persistent Molecular Disease in Adult Patients With AML Evaluated With Whole-Exome and Targeted Error-Corrected DNA Sequencing.

PURPOSE: Persistent molecular disease (PMD) after induction chemotherapy predicts relapse in AML. In this study, we used whole-exome sequencing (WES) and targeted error-corrected sequencing to assess the frequency and mutational patterns of PMD in 30 patients with AML. MATERIALS AND METHODS: The study cohort included 30 patients with adult AML younger than 65 years who were uniformly treated with standard induction chemotherapy. Tumor/normal WES was performed for all patients at presentation. PMD analysis was evaluated in bone marrow samples obtained during clinicopathologic remission using repeat WES and analysis of patient-specific mutations and error-corrected sequencing of 40 recurrently mutated AML genes (MyeloSeq). RESULTS: WES for patient-specific mutations detected PMD in 63% of patients (19/30) using a minimum variant allele fraction (VAF) of 2.5%. In comparison, MyeloSeq identified persistent mutations above 0.1% VAF in 77% of patients (23/30). PMD was usually present at relatively high levels (>2.5% VAFs), such that WES and MyeloSeq agreed for 73% of patients despite differences in detection limits. Mutations in DNMT3A, ASXL1, and TET2 (ie, DTA mutations) were persistent in 16 of 17 patients, but WES also detected non-DTA mutations in 14 of these patients, which for some patients distinguished residual AML cells from clonal hematopoiesis. Surprisingly, MyeloSeq detected additional variants not identified at presentation in 73% of patients that were consistent with new clonal cell populations after chemotherapy. CONCLUSION: PMD and clonal hematopoiesis are both common in patients with AML in first remission. These findings demonstrate the importance of baseline testing for accurate interpretation of mutation-based tumor monitoring assays for patients with AML and highlight the need for clinical trials to determine whether these complex mutation patterns correlate with clinical outcomes in AML.

Genomic landscape of TP53 -mutated myeloid malignancies.

TP53 -mutated myeloid malignancies are most frequently associated with complex cytogenetics. The presence of complex and extensive structural variants complicates detailed genomic analysis by conventional clinical techniques. We performed whole genome sequencing of 42 AML/MDS cases with paired normal tissue to characterize the genomic landscape of TP53 -mutated myeloid malignancies. The vast majority of cases had multi-hit involvement at the TP53 genetic locus (94%), as well as aneuploidy and chromothripsis. Chromosomal patterns of aneuploidy differed significantly from TP53 -mutated cancers arising in other tissues. Recurrent structural variants affected regions that include ETV6 on chr12p, RUNX1 on chr21, and NF1 on chr17q. Most notably for ETV6 , transcript expression was low in cases of TP53 -mutated myeloid malignancies both with and without structural rearrangements involving chromosome 12p. Telomeric content is increased in TP53 -mutated AML/MDS compared other AML subtypes, and telomeric content was detected adjacent to interstitial regions of chromosomes. The genomic landscape of TP53 -mutated myeloid malignancies reveals recurrent structural variants affecting key hematopoietic transcription factors and telomeric repeats that are generally not detected by panel sequencing or conventional cytogenetic analyses. Key Points: WGS comprehensively determines TP53 mutation status, resulting in the reclassification of 12% of cases from mono-allelic to multi-hit Chromothripsis is more frequent than previously appreciated, with a preference for specific chromosomes ETV6 is deleted in 45% of cases, with evidence for epigenetic suppression in non-deleted cases NF1 is mutated in 48% of cases, with multi-hit mutations in 17% of these cases TP53 -mutated AML/MDS is associated with altered telomere content compared with other AMLs.

Somatic Dnmt3a inactivation leads to slow, canonical DNA methylation loss in murine hematopoietic cells.

Mutations in the gene encoding DNA methyltransferase 3A (DNMT3A) are the most common cause of clonal hematopoiesis and are among the most common initiating events of acute myeloid leukemia (AML). Studies in germline and somatic Dnmt3a knockout mice have identified focal, canonical hypomethylation phenotypes in hematopoietic cells; however, the kinetics of methylation loss following acquired DNMT3A inactivation in hematopoietic cells is essentially unknown. Therefore, we evaluated a somatic, inducible model of hematopoietic Dnmt3a loss, and show that inactivation of Dnmt3a in murine hematopoietic cells results in a relatively slow loss of methylation at canonical sites throughout the genome; in contrast, remethylation of Dnmt3a deficient genomes in hematopoietic cells occurs much more quickly. This data suggests that slow methylation loss may contribute, at least in part, to the long latent period that characterizes clonal expansion and leukemia development in individuals with acquired DNMT3A mutations in hematopoietic stem cells.

Functional and epigenetic phenotypes of humans and mice with DNMT3A Overgrowth Syndrome.

Germline pathogenic variants in DNMT3A were recently described in patients with overgrowth, obesity, behavioral, and learning difficulties (DNMT3A Overgrowth Syndrome/DOS). Somatic mutations in the DNMT3A gene are also the most common cause of clonal hematopoiesis, and can initiate acute myeloid leukemia (AML). Using whole genome bisulfite sequencing, we studied DNA methylation in peripheral blood cells of 11 DOS patients and found a focal, canonical hypomethylation phenotype, which is most severe with the dominant negative DNMT3AR882H mutation. A germline mouse model expressing the homologous Dnmt3aR878H mutation phenocopies most aspects of the human DOS syndrome, including the methylation phenotype and an increased incidence of spontaneous hematopoietic malignancies, suggesting that all aspects of this syndrome are caused by this mutation.

Mitochondrial genome copy number measured by DNA sequencing in human blood is strongly associated with metabolic traits via cell-type composition differences.

BACKGROUND: Mitochondrial genome copy number (MT-CN) varies among humans and across tissues and is highly heritable, but its causes and consequences are not well understood. When measured by bulk DNA sequencing in blood, MT-CN may reflect a combination of the number of mitochondria per cell and cell-type composition. Here, we studied MT-CN variation in blood-derived DNA from 19184 Finnish individuals using a combination of genome (N = 4163) and exome sequencing (N = 19034) data as well as imputed genotypes (N = 17718). RESULTS: We identified two loci significantly associated with MT-CN variation: a common variant at the MYB-HBS1L locus (P = 1.6 × 10-8), which has previously been associated with numerous hematological parameters; and a burden of rare variants in the TMBIM1 gene (P = 3.0 × 10-8), which has been reported to protect against non-alcoholic fatty liver disease. We also found that MT-CN is strongly associated with insulin levels (P = 2.0 × 10-21) and other metabolic syndrome (metS)-related traits. Using a Mendelian randomization framework, we show evidence that MT-CN measured in blood is causally related to insulin levels. We then applied an MT-CN polygenic risk score (PRS) derived from Finnish data to the UK Biobank, where the association between the PRS and metS traits was replicated. Adjusting for cell counts largely eliminated these signals, suggesting that MT-CN affects metS via cell-type composition. CONCLUSION: These results suggest that measurements of MT-CN in blood-derived DNA partially reflect differences in cell-type composition and that these differences are causally linked to insulin and related traits.

Mutation Clearance after Transplantation for Myelodysplastic Syndrome.

BACKGROUND: Allogeneic hematopoietic stem-cell transplantation is the only curative treatment for patients with myelodysplastic syndrome (MDS). The molecular predictors of disease progression after transplantation are unclear. METHODS: We sequenced bone marrow and skin samples from 90 adults with MDS who underwent allogeneic hematopoietic stem-cell transplantation after a myeloablative or reduced-intensity conditioning regimen. We detected mutations before transplantation using enhanced exome sequencing, and we evaluated mutation clearance by using error-corrected sequencing to genotype mutations in bone marrow samples obtained 30 days after transplantation. In this exploratory study, we evaluated the association of a mutation detected after transplantation with disease progression and survival. RESULTS: Sequencing identified at least one validated somatic mutation before transplantation in 86 of 90 patients (96%); 32 of these patients (37%) had at least one mutation with a maximum variant allele frequency of at least 0.5% (equivalent to 1 heterozygous mutant cell in 100 cells) 30 days after transplantation. Patients with disease progression had mutations with a higher maximum variant allele frequency at 30 days than those who did not (median maximum variant allele frequency, 0.9% vs. 0%; P<0.001). The presence of at least one mutation with a variant allele frequency of at least 0.5% at day 30 was associated with a higher risk of progression (53.1% vs. 13.0%; conditioning regimen-adjusted hazard ratio, 3.86; 95% confidence interval [CI], 1.96 to 7.62; P<0.001) and a lower 1-year rate of progression-free survival than the absence of such a mutation (31.3% vs. 59.3%; conditioning regimen-adjusted hazard ratio for progression or death, 2.22; 95% CI, 1.32 to 3.73; P=0.005). The rate of progression-free survival was lower among patients who had received a reduced-intensity conditioning regimen and had at least one persistent mutation with a variant allele frequency of at least 0.5% at day 30 than among patients with other combinations of conditioning regimen and mutation status (P≤0.001). Multivariate analysis confirmed that patients who had a mutation with a variant allele frequency of at least 0.5% detected at day 30 had a higher risk of progression (hazard ratio, 4.48; 95% CI, 2.21 to 9.08; P<0.001) and a lower 1-year rate of progression-free survival than those who did not (hazard ratio for progression or death, 2.39; 95% CI, 1.40 to 4.09; P=0.002). CONCLUSIONS: The risk of disease progression was higher among patients with MDS in whom persistent disease-associated mutations were detected in the bone marrow 30 days after transplantation than among those in whom these mutations were not detected. (Funded by the Leukemia and Lymphoma Society and others.).

Genomic heterogeneity of ALK fusion breakpoints in non-small-cell lung cancer.

In lung adenocarcinoma, canonical EML4-ALK inversion results in a fusion protein with a constitutively active ALK kinase domain. Evidence of ALK rearrangement occurs in a minority (2-7%) of lung adenocarcinoma, and only ~60% of these patients will respond to targeted ALK inhibition by drugs such as crizotinib and ceritinib. Clinically, targeted anti-ALK therapy is often initiated based on evidence of an ALK genomic rearrangement detected by fluorescence in situ hybridization (FISH) of interphase cells in formalin-fixed, paraffin-embedded tissue sections. At the genomic level, however, ALK rearrangements are heterogeneous, with multiple potential breakpoints in EML4, and alternate fusion partners. Using next-generation sequencing of DNA and RNA together with ALK immunohistochemistry, we comprehensively characterized genomic breakpoints in 33 FISH-positive lung adenocarcinomas. Of these 33 cases, 29 (88%) had detectable DNA level ALK rearrangements involving EML4, KIF5B, or non-canonical partners including ASXL2, ATP6V1B1, PRKAR1A, and SPDYA. A subset of 12 cases had material available for RNA-Seq. Of these, eight of eight (100%) cases with DNA rearrangements showed ALK fusion transcripts from RNA-Seq; three of four cases (75%) without detectable DNA rearrangements were similarly negative by RNA-Seq, and one case was positive by RNA-Seq but negative by DNA next-generation sequencing. By immunohistochemistry, 17 of 19 (89%) tested cases were clearly positive for ALK protein expression; the remaining cases had no detectable DNA level rearrangement or had a non-canonical rearrangement not predicted to form a fusion protein. Survival analysis of patients treated with targeted ALK inhibitors demonstrates a significant difference in mean survival between patients with next-generation sequencing confirmed EML4-ALK rearrangements, and those without (20.6 months vs 5.4 months, P<0.01). Together, these data demonstrate abundant genomic heterogeneity among ALK-rearranged lung adenocarcinoma, which may account for differences in treatment response with targeted ALK inhibitors.

Spectrum of mutations in leiomyosarcomas identified by clinical targeted next-generation sequencing.

Recurrent genomic mutations in uterine and non-uterine leiomyosarcomas have not been well established. Using a next generation sequencing (NGS) panel of common cancer-associated genes, 25 leiomyosarcomas arising from multiple sites were examined to explore genetic alterations, including single nucleotide variants (SNV), small insertions/deletions (indels), and copy number alterations (CNA). Sequencing showed 86 non-synonymous, coding region somatic variants within 151 gene targets in 21 cases, with a mean of 4.1 variants per case; 4 cases had no putative mutations in the panel of genes assayed. The most frequently altered genes were TP53 (36%), ATM and ATRX (16%), and EGFR and RB1 (12%). CNA were identified in 85% of cases, with the most frequent copy number losses observed in chromosomes 10 and 13 including PTEN and RB1; the most frequent gains were seen in chromosomes 7 and 17. Our data show that deletions in canonical cancer-related genes are common in leiomyosarcomas. Further, the spectrum of gene mutations observed shows that defects in DNA repair and chromosomal maintenance are central to the biology of leiomyosarcomas, and that activating mutations observed in other common cancer types are rare in leiomyosarcomas.

In Silico Proficiency Testing for Clinical Next-Generation Sequencing.

Quality assurance for clinical next-generation sequencing (NGS)-based assays is difficult given the complex methods and the range of sequence variants such assays can detect. As the number and range of mutations detected by clinical NGS assays has increased, it is difficult to apply standard analyte-specific proficiency testing (PT). Most current proficiency testing challenges for NGS are methods-based PT surveys that use DNA from reference samples engineered to harbor specific mutations that test both sequence generation and bioinformatics analysis. These methods-based PTs are limited by the number and types of mutations that can be physically introduced into a single DNA sample. In silico proficiency testing, which evaluates only the bioinformatics component of NGS assays, is a recently introduced PT method that allows for evaluation of numerous mutations spanning a range of variant classes. In silico PT data sets can be generated from simulated or actual sequencing data and are used to test alignment through variant detection and annotation steps. In silico PT has several advantages over the use of physical samples, including greater flexibility in tested variants, the ability to design laboratory-specific challenges, and lower costs. Herein, we review the use of in silico PT as an alternative to traditional methods-based PT as it is evolving in oncology applications and discuss how the approach is applicable more broadly.

Occult Specimen Contamination in Routine Clinical Next-Generation Sequencing Testing.

OBJECTIVES: To evaluate the extent of human-to-human specimen contamination in clinical next-generation sequencing (NGS) data. METHODS: Using haplotype analysis to detect specimen admixture, with orthogonal validation by short tandem repeat analysis, we determined the rate of clinically significant (>5%) DNA contamination in clinical NGS data from 296 consecutive cases. Haplotype analysis was performed using read haplotypes at common, closely spaced single-nucleotide polymorphisms in low linkage disequilibrium in the population, which were present in regions targeted by the clinical assay. Percent admixture was estimated based on frequencies of the read haplotypes at loci that showed evidence for contamination. RESULTS: We identified nine (3%) cases with at least 5% DNA admixture. Three cases were bone marrow transplant patients known to be chimeric. Six admixed cases were incidents of contamination, and the rate of contamination was strongly correlated with DNA yield from the tissue specimen. CONCLUSIONS: Human-human specimen contamination occurs in clinical NGS testing. Tools for detecting contamination in NGS sequence data should be integrated into clinical bioinformatics pipelines, especially as laboratories trend toward using smaller amounts of input DNA and reporting lower frequency variants. This study provides one estimate of the rate of clinically significant human-human specimen contamination in clinical NGS testing.

T-cell clonality assessment by next-generation sequencing improves detection sensitivity in mycosis fungoides.

BACKGROUND: T-cell receptor (TCR) clonality assessment is a principal diagnostic test in the management of mycosis fungoides (MF). However, current polymerase chain reaction-based methods may produce ambiguous results, often because of low abundance of clonal T lymphocytes, resulting in weak clonal peaks that cannot be size-resolved by contemporary capillary electrophoresis (CE). OBJECTIVE: We sought to determine if next-generation sequencing (NGS)-based detection has increased sensitivity for T-cell clonality over CE-based detection in MF. METHODS: Clonality was determined by an NGS-based method in which the TCR-γ variable region was polymerase chain reaction amplified and the products sequenced to establish the identity of rearranged variable and joining regions. RESULTS: Of the 35 MF cases tested, 29 (85%) showed a clonal T-cell rearrangement by NGS, compared with 15 (44%) by standard CE detection. Three patients with MF had follow-up testing that showed identical, clonal TCR sequences in subsequent skin biopsy specimens. LIMITATIONS: Clonal T-cell populations have been described in benign conditions; evidence of clonality alone, by any method, is not sufficient for diagnosis. CONCLUSION: TCR clonality assessment by NGS has superior sensitivity compared with CE-based detection. Further, NGS enables tracking of specific clones across multiple time points for more accurate identification of recurrent MF.

Whole Exome Sequencing Reveals the Order of Genetic Changes during Malignant Transformation and Metastasis in a Single Patient with NF1-plexiform Neurofibroma.

PURPOSE: Malignant peripheral nerve sheath tumors (MPNST) occur at increased frequency in individuals with neurofibromatosis type 1 (NF1), where they likely arise from benign plexiform neurofibroma precursors. While previous studies have used a variety of discovery approaches to discover genes associated with MPNST pathogenesis, it is currently unclear what molecular events are associated with the evolution of MPNST from plexiform neurofibroma. EXPERIMENTAL DESIGN: Whole-exome sequencing was performed on biopsy materials representing plexiform neurofibroma (n = 3), MPNST, and metastasis from a single individual with NF1 over a 14-year period. Additional validation cases were used to assess candidate genes involved in malignant progression, while a murine MPNST model was used for functional analysis. RESULTS: There was an increasing proportion of cells with a somatic NF1 gene mutation as the tumors progressed from benign to malignant, suggesting a clonal process in MPNST development. Copy number variations, including loss of one copy of the TP53 gene, were identified in the primary tumor and the metastatic lesion, but not in benign precursor lesions. A limited number of genes with nonsynonymous somatic mutations (ßIII-spectrin and ZNF208) were discovered, several of which were validated in additional primary and metastatic MPNST samples. Finally, increased ßIII-spectrin expression was observed in the majority of MPNSTs, and shRNA-mediated knockdown reduced murine MPNST growth in vivo. CONCLUSIONS: Collectively, the ability to track the molecular evolution of MPNST in a single individual with NF1 offers new insights into the sequence of genetic events important for disease pathogenesis and progression for future mechanistic study.

A wide spectrum of EGFR mutations in glioblastoma is detected by a single clinical oncology targeted next-generation sequencing panel.

With the advent of large-scale genomic analysis, the genetic landscape of glioblastoma (GBM) has become more clear, including characteristic genetic alterations in EGFR. In routine clinical practice, genetic alterations in GBMs are identified using several disparate techniques that consume already limited amounts of tissue and add to overall testing costs. In this study, we sought to determine if the full spectrum of EGFR mutations in GBMs could be detected using a single next generation sequencing (NGS) based oncology assay in 34 consecutive cases. Using a battery of informatics tools to identify single nucleotide variants, insertions and deletions, and amplification (including variants EGFRvIII and EGFRvV), twenty-one of the 34 (62%) individuals had at least one alteration in EGFR by sequencing, consistent with published datasets. Mutations detected include several single nucleotide variants, amplification (confirmed by fluorescence in situ hybridization), and the variants EGFRvIII and EGFRvV (confirmed by multiplex ligation-dependent probe amplification). Here we show that a single NGS assay can identify the full spectrum of relevant EGFR mutations. Overall, sequencing based diagnostics have the potential to maximize the amount of genetic information obtained from GBMs and simultaneously reduce the total time, required specimen material, and costs associated with current multimodality studies.

Extended maternal age at birth of last child and women's longevity in the Long Life Family Study.

OBJECTIVE: This study investigated the association between maternal age at birth of last child and likelihood of survival to advanced age. METHODS: This was a nested case-control study using Long Life Family Study data. Three hundred eleven women who survived past the oldest 5th percentile of survival (according to birth cohort-matched life tables) were identified as cases, and 151 women who died at ages younger than the top 5th percentile of survival were identified as controls. A Bayesian mixed-effect logistic regression model was used to estimate the association between maternal age at birth of last child and exceptional longevity among these 462 women. RESULTS: We found a significant association for older maternal age, whereby women who had their last child beyond age 33 years had twice the odds for survival to the top 5th percentile of survival for their birth cohorts compared with women who had their last child by age 29 years (age between 33 and 37 y: odds ratio, 2.08; 95% CI, 1.13 to 3.92; older age: odds ratio, 1.92; 95% CI, 1.03 to 3.68). CONCLUSIONS: This study supports findings from other studies demonstrating a positive association between older maternal age and greater odds for surviving to an unusually old age.

Admixture mapping identifies a quantitative trait locus associated with FEV1/FVC in the COPDGene Study.

African Americans are admixed with genetic contributions from European and African ancestral populations. Admixture mapping leverages this information to map genes influencing differential disease risk across populations. We performed admixture and association mapping in 3,300 African American current or former smokers from the COPDGene Study. We analyzed estimated local ancestry and SNP genotype information to identify regions associated with FEV1 /FVC, the ratio of forced expiratory volume in one second to forced vital capacity, measured by spirometry performed after bronchodilator administration. Global African ancestry inversely associated with FEV1 /FVC (P = 0.035). Genome-wide admixture analysis, controlling for age, gender, body mass index, current smoking status, pack-years smoked, and four principal components summarizing the genetic background of African Americans in the COPDGene Study, identified a region on chromosome 12q14.1 associated with FEV1 /FVC (P = 2.1 × 10(-6) ) when regressed on local ancestry. Allelic association in this region of chromosome 12 identified an intronic variant in FAM19A2 (rs348644) as associated with FEV1 /FVC (P = 1.76 × 10(-6) ). By combining admixture and association mapping, a marker on chromosome 12q14.1 was identified as being associated with reduced FEV1 /FVC ratio among African Americans in the COPDGene Study.

Copy number variants in clinical next-generation sequencing data can define the relationship between simultaneous tumors in an individual patient.

Targeted next-generation sequencing (NGS) cancer panels have become a popular method for the identification of clinically predictive mutations in cancer. Such methods typically detect single nucleotide variants (SNVs) and small insertions/deletions (indels) in known cancer genes and can provide further information regarding diagnosis in challenging surgical pathology cases, as well as identify therapeutic targets and prognostically significant mutations. However, in addition to SNVs and indels, other mutation classes, including copy number variants (CNVs) and translocations, can be simultaneously detected from targeted NGS data. Here, as proof of methods, we present clinical data which demonstrate that targeted NGS panels can separate synchronous liver tumors based on CNV status, in the absence of distinct SNVs and indels. Such CNV-based analysis can be performed without additional cost using existing targeted cancer panel data and publically available software.

Detection of gene rearrangements in targeted clinical next-generation sequencing.

The identification of recurrent gene rearrangements in the clinical laboratory is the cornerstone for risk stratification and treatment decisions in many malignant tumors. Studies have reported that targeted next-generation sequencing assays have the potential to identify such rearrangements; however, their utility in the clinical laboratory is unknown. We examine the sensitivity and specificity of ALK and KMT2A (MLL) rearrangement detection by next-generation sequencing in the clinical laboratory. We analyzed a series of seven ALK rearranged cancers, six KMT2A rearranged leukemias, and 77 ALK/KMT2A rearrangement-negative cancers, previously tested by fluorescence in situ hybridization (FISH). Rearrangement detection was tested using publicly available software tools, including Breakdancer, ClusterFAST, CREST, and Hydra. Using Breakdancer and ClusterFAST, we detected ALK rearrangements in seven of seven FISH-positive cases and KMT2A rearrangements in six of six FISH-positive cases. Among the 77 ALK/KMT2A FISH-negative cases, no false-positive identifications were made by Breakdancer or ClusterFAST. Further, we identified one ALK rearranged case with a noncanonical intron 16 breakpoint, which is likely to affect its response to targeted inhibitors. We report that clinically relevant chromosomal rearrangements can be detected from targeted gene panel-based next-generation sequencing with sensitivity and specificity equivalent to that of FISH while providing finer-scale information and increased efficiency for molecular oncology testing.