Mid-pass Whole-genome Sequencing in a Malagasy Cohort Uncovers Body Composition Associations.

The majority of human genomic research studies have been conducted in European-ancestry cohorts, reducing the likelihood of detecting potentially novel and globally impactful findings. Here, we present mid-pass whole-genome sequencing data and a genome-wide association study in a cohort of 264 self-reported Malagasy individuals from three locations on the island of Madagascar. We describe genetic variation in this Malagasy cohort, providing insight into the shared and unique patterns of genetic variation across the island. We observe phenotypic variation by location, and find high rates of hypertension particularly in the Southern Highlands sampling site as well as elevated self-reported malaria prevalence in the West Coast site relative to other sites. After filtering to a subset of 214 minimally-related individuals, we find a number of genetic associations with body composition traits, including many variants that are only observed in African populations or populations with admixed African ancestry from the 1000 Genomes Project. This study highlights the importance of including diverse populations in genomic research for the potential to gain novel insights, even with small cohort sizes. This project was conducted in partnership and consultation with local stakeholders in Madagascar and serves as an example of genomic research that prioritizes community engagement and that has potential impacts on our understanding of human health and disease.

The gout epidemic in French Polynesia: a modelling study of data from the Ma'i u'u epidemiological survey.

BACKGROUND: Gout is the most common cause of inflammatory arthritis worldwide, particularly in Pacific regions. We aimed to establish the prevalence of gout and hyperuricaemia in French Polynesia, their associations with dietary habits, their comorbidities, the prevalence of the HLA-B*58:01 allele, and current management of the disease. METHODS: The Ma'i u'u survey was epidemiological, prospective, cross-sectional, and gout-focused and included a random sample of adults from the general adult population of French Polynesia. It was conducted and data were collected between April 13 and Aug 16, 2021. Participants were randomly selected to represent the general adult population of French Polynesia on the basis of housing data collected during the 2017 territorial census. Each selected household was visited by a research nurse from the Ma'i u'u survey who collected data via guided, 1-h interviews with participants. In each household, the participant was the individual older than 18 years with the closest upcoming birthday. To estimate the frequency of HLA-B*58:01, we estimated HLA-B haplotypes on individuals who had whole-genome sequencing to approximately 5× average coverage (mid-pass sequencing). A subset of individuals who self-reported Polynesian ancestry and not European, Chinese, or other ancestry were used to estimate Polynesian-ancestry specific allele frequencies. Bivariate associations were reported for weighted participants; effect sizes were estimated through the odds ratio (OR) of the association calculated on the basis of a logistic model fitted with weighted observations. FINDINGS: Among the random sample of 2000 households, 896 participants were included, 140 individuals declined, and 964 households could not be contacted. 22 participants could not be weighted due to missing data, so the final weighted analysis included 874 participants (449 [51·4%] were female and 425 [48·6%] were male) representing the 196 630 adults living in French Polynesia. The estimated prevalence of gout was 14·5% (95% CI 9·9-19·2), representing 28 561 French Polynesian adults, that is 25·5% (18·2-32·8) of male individuals and 3·5% (1·0-6·0) of female individuals. The prevalence of hyperuricaemia was estimated at 71·6% (66·7-76·6), representing 128 687 French Polynesian adults. In multivariable analysis, age (OR 1·5, 95% CI 1·2-1·8 per year), male sex (10·3, 1·8-60·7), serum urate (1·6, 1·3-2·0 per 1 mg/dL), uraturia (0·8, 0·8-0·8 per 100 mg/L), type 2 diabetes (2·1, 1·4-3·1), BMI more than 30 kg/m2 (1·1, 1·0-1·2 per unit), and percentage of visceral fat (1·7, 1·1-2·7 per 1% increase) were associated with gout. There were seven heterozygous HLA-B*58:01 carriers in the full cohort of 833 individuals (seven [0·4%] of 1666 total alleles) and two heterozygous carriers in a subset of 696 individuals of Polynesian ancestry (two [0·1%]). INTERPRETATION: French Polynesia has an estimated high prevalence of gout and hyperuricaemia, with gout affecting almost 15% of adults. Territorial measures that focus on increasing access to effective urate-lowering therapies are warranted to control this major public health problem. FUNDING: Variant Bio, the French Polynesian Health Administration, Lille Catholic University Hospitals, French Society of Rheumatology, and Novartis.

Somatic whole genome dynamics of precancer in Barrett's esophagus reveals features associated with disease progression.

While the genomes of normal tissues undergo dynamic changes over time, little is understood about the temporal-spatial dynamics of genomes in premalignant tissues that progress to cancer compared to those that remain cancer-free. Here we use whole genome sequencing to contrast genomic alterations in 427 longitudinal samples from 40 patients with stable Barrett's esophagus compared to 40 Barrett's patients who progressed to esophageal adenocarcinoma (ESAD). We show the same somatic mutational processes are active in Barrett's tissue regardless of outcome, with high levels of mutation, ESAD gene and focal chromosomal alterations, and similar mutational signatures. The critical distinction between stable Barrett's versus those who progress to cancer is acquisition and expansion of TP53-/- cell populations having complex structural variants and high-level amplifications, which are detectable up to six years prior to a cancer diagnosis. These findings reveal the timing of common somatic genome dynamics in stable Barrett's esophagus and define key genomic features specific to progression to esophageal adenocarcinoma, both of which are critical for cancer prevention and early detection strategies.

Diverse tumorigenic consequences of human papillomavirus integration in primary oropharyngeal cancers.

Human papillomavirus (HPV) causes 5% of all cancers and frequently integrates into host chromosomes. The HPV oncoproteins E6 and E7 are necessary but insufficient for cancer formation, indicating that additional secondary genetic events are required. Here, we investigate potential oncogenic impacts of virus integration. Analysis of 105 HPV-positive oropharyngeal cancers by whole-genome sequencing detects virus integration in 77%, revealing five statistically significant sites of recurrent integration near genes that regulate epithelial stem cell maintenance (i.e., SOX2, TP63, FGFR, MYC) and immune evasion (i.e., CD274). Genomic copy number hyperamplification is enriched 16-fold near HPV integrants, and the extent of focal host genomic instability increases with their local density. The frequency of genes expressed at extreme outlier levels is increased 86-fold within ±150 kb of integrants. Across 95% of tumors with integration, host gene transcription is disrupted via intragenic integrants, chimeric transcription, outlier expression, gene breaking, and/or de novo expression of noncoding or imprinted genes. We conclude that virus integration can contribute to carcinogenesis in a large majority of HPV-positive oropharyngeal cancers by inducing extensive disruption of host genome structure and gene expression.

Mid-pass whole genome sequencing enables biomedical genetic studies of diverse populations.

BACKGROUND: Historically, geneticists have relied on genotyping arrays and imputation to study human genetic variation. However, an underrepresentation of diverse populations has resulted in arrays that poorly capture global genetic variation, and a lack of reference panels. This has contributed to deepening global health disparities. Whole genome sequencing (WGS) better captures genetic variation but remains prohibitively expensive. Thus, we explored WGS at "mid-pass" 1-7x coverage. RESULTS: Here, we developed and benchmarked methods for mid-pass sequencing. When applied to a population without an existing genomic reference panel, 4x mid-pass performed consistently well across ethnicities, with high recall (98%) and precision (97.5%). CONCLUSION: Compared to array data imputed into 1000 Genomes, mid-pass performed better across all metrics and identified novel population-specific variants with potential disease relevance. We hope our work will reduce financial barriers for geneticists from underrepresented populations to characterize their genomes prior to biomedical genetic applications.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

Author Correction: PGBD5 promotes site-specific oncogenic mutations in human tumors.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Correction to: Sequencing and curation strategies for identifying candidate glioblastoma treatments.

Following publication of the original article [1], it was reported that the given name of the fourteenth author was incorrectly published. The incorrect and the correct names are given below.

Sequencing and curation strategies for identifying candidate glioblastoma treatments.

BACKGROUND: Prompted by the revolution in high-throughput sequencing and its potential impact for treating cancer patients, we initiated a clinical research study to compare the ability of different sequencing assays and analysis methods to analyze glioblastoma tumors and generate real-time potential treatment options for physicians. METHODS: A consortium of seven institutions in New York City enrolled 30 patients with glioblastoma and performed tumor whole genome sequencing (WGS) and RNA sequencing (RNA-seq; collectively WGS/RNA-seq); 20 of these patients were also analyzed with independent targeted panel sequencing. We also compared results of expert manual annotations with those from an automated annotation system, Watson Genomic Analysis (WGA), to assess the reliability and time required to identify potentially relevant pharmacologic interventions. RESULTS: WGS/RNAseq identified more potentially actionable clinical results than targeted panels in 90% of cases, with an average of 16-fold more unique potentially actionable variants identified per individual; 84 clinically actionable calls were made using WGS/RNA-seq that were not identified by panels. Expert annotation and WGA had good agreement on identifying variants [mean sensitivity = 0.71, SD = 0.18 and positive predictive value (PPV) = 0.80, SD = 0.20] and drug targets when the same variants were called (mean sensitivity = 0.74, SD = 0.34 and PPV = 0.79, SD = 0.23) across patients. Clinicians used the information to modify their treatment plan 10% of the time. CONCLUSION: These results present the first comprehensive comparison of technical and machine augmented analysis of targeted panel and WGS/RNA-seq to identify potential cancer treatments.

Genetic mechanisms of primary chemotherapy resistance in pediatric acute myeloid leukemia.

Acute myeloid leukemias (AML) are characterized by mutations of tumor suppressor and oncogenes, involving distinct genes in adults and children. While certain mutations have been associated with the increased risk of AML relapse, the genomic landscape of primary chemotherapy-resistant AML is not well defined. As part of the TARGET initiative, we performed whole-genome DNA and transcriptome RNA and miRNA sequencing analysis of pediatric AML with failure of induction chemotherapy. We identified at least three genetic groups of patients with induction failure, including those with NUP98 rearrangements, somatic mutations of WT1 in the absence of apparent NUP98 mutations, and additional recurrent variants including those in KMT2C and MLLT10. Comparison of specimens before and after chemotherapy revealed distinct and invariant gene expression programs. While exhibiting overt therapy resistance, these leukemias nonetheless showed diverse forms of clonal evolution upon chemotherapy exposure. This included selection for mutant alleles of FRMD8, DHX32, PIK3R1, SHANK3, MKLN1, as well as persistence of WT1 and TP53 mutant clones, and elimination of FLT3, PTPN11, and NRAS mutant clones. These findings delineate genetic mechanisms of primary chemotherapy resistance in pediatric AML, which should inform improved approaches for its diagnosis and therapy.

Human papillomavirus and the landscape of secondary genetic alterations in oral cancers.

Human papillomavirus (HPV) is a necessary but insufficient cause of a subset of oral squamous cell carcinomas (OSCCs) that is increasing markedly in frequency. To identify contributory, secondary genetic alterations in these cancers, we used comprehensive genomics methods to compare 149 HPV-positive and 335 HPV-negative OSCC tumor/normal pairs. Different behavioral risk factors underlying the two OSCC types were reflected in distinctive genomic mutational signatures. In HPV-positive OSCCs, the signatures of APOBEC cytosine deaminase editing, associated with anti-viral immunity, were strongly linked to overall mutational burden. In contrast, in HPV-negative OSCCs, T>C substitutions in the sequence context 5'-ATN-3' correlated with tobacco exposure. Universal expression of HPV E6*1 and E7 oncogenes was a sine qua non of HPV-positive OSCCs. Significant enrichment of somatic mutations was confirmed or newly identified in PIK3CA, KMT2D, FGFR3, FBXW7, DDX3X, PTEN, TRAF3, RB1, CYLD, RIPK4, ZNF750, EP300, CASZ1, TAF5, RBL1, IFNGR1, and NFKBIA Of these, many affect host pathways already targeted by HPV oncoproteins, including the p53 and pRB pathways, or disrupt host defenses against viral infections, including interferon (IFN) and nuclear factor kappa B signaling. Frequent copy number changes were associated with concordant changes in gene expression. Chr 11q (including CCND1) and 14q (including DICER1 and AKT1) were recurrently lost in HPV-positive OSCCs, in contrast to their gains in HPV-negative OSCCs. High-ranking variant allele fractions implicated ZNF750, PIK3CA, and EP300 mutations as candidate driver events in HPV-positive cancers. We conclude that virus-host interactions cooperatively shape the unique genetic features of these cancers, distinguishing them from their HPV-negative counterparts.

Genome-wide somatic variant calling using localized colored de Bruijn graphs.

Reliable detection of somatic variations is of critical importance in cancer research. Here we present Lancet, an accurate and sensitive somatic variant caller, which detects SNVs and indels by jointly analyzing reads from tumor and matched normal samples using colored de Bruijn graphs. We demonstrate, through extensive experimental comparison on synthetic and real whole-genome sequencing datasets, that Lancet has better accuracy, especially for indel detection, than widely used somatic callers, such as MuTect, MuTect2, LoFreq, Strelka, and Strelka2. Lancet features a reliable variant scoring system, which is essential for variant prioritization, and detects low-frequency mutations without sacrificing the sensitivity to call longer insertions and deletions empowered by the local-assembly engine. In addition to genome-wide analysis, Lancet allows inspection of somatic variants in graph space, which augments the traditional read alignment visualization to help confirm a variant of interest. Lancet is available as an open-source program at https://github.com/nygenome/lancet.

Analytical Validation of Clinical Whole-Genome and Transcriptome Sequencing of Patient-Derived Tumors for Reporting Targetable Variants in Cancer.

We developed and validated a clinical whole-genome and transcriptome sequencing (WGTS) assay that provides a comprehensive genomic profile of a patient's tumor. The ability to fully capture the mappable genome with sufficient sequencing coverage to precisely call DNA somatic single nucleotide variants, insertions/deletions, copy number variants, structural variants, and RNA gene fusions was analyzed. New York State's Department of Health next-generation DNA sequencing guidelines were expanded for establishing performance validation applicable to whole-genome and transcriptome sequencing. Whole-genome sequencing laboratory protocols were validated for the Illumina HiSeq X Ten platform and RNA sequencing for Illumina HiSeq2500 platform for fresh or frozen and formalin-fixed, paraffin-embedded tumor samples. Various bioinformatics tools were also tested, and CIs for sensitivity and specificity thresholds in calling clinically significant somatic aberrations were determined. The validation was performed on a set of 125 tumor normal pairs. RNA sequencing was performed to call fusions and to confirm the DNA variants or exonic alterations. Here, we present our results and WGTS standards for variant allele frequency, reproducibility, analytical sensitivity, and present limit of detection analysis for single nucleotide variant calling, copy number identification, and structural variants. We show that The New York Genome Center WGTS clinical assay can provide a comprehensive patient variant discovery approach suitable for directed oncologic therapeutic applications.

Whole Genome Sequencing-Based Discovery of Structural Variants in Glioblastoma.

Next-generation DNA sequencing (NGS) technologies are currently being applied in both research and clinical settings for the understanding and management of disease. The goal is to use high-throughput sequencing to identify specific variants that drive tumorigenesis within each individual's tumor genomic profile. The significance of copy number and structural variants in glioblastoma makes it essential to broaden the search beyond oncogenic single nucleotide variants toward whole genome profiles of genetic aberrations that may contribute to disease progression. The heterogeneity of glioblastoma and its variability of cancer driver mutations necessitate a more robust examination of a patient's tumor genome. Here, we present patient whole genome sequencing (WGS) information to identify oncogenic structural variants that may contribute to glioblastoma pathogenesis. We provide WGS protocols and bioinformatics approaches to identify copy number and structural variations in 41 glioblastoma patient samples. We present how WGS can identify structural diversity within glioblastoma samples. We specifically show how to apply current bioinformatics tools to detect EGFR variants and other structural aberrations from DNA whole genome sequencing and how to validate those variants within the laboratory. These comprehensive WGS protocols can provide additional information directing more precise therapeutic options in the treatment of glioblastoma.

Comparing sequencing assays and human-machine analyses in actionable genomics for glioblastoma.

OBJECTIVE: To analyze a glioblastoma tumor specimen with 3 different platforms and compare potentially actionable calls from each. METHODS: Tumor DNA was analyzed by a commercial targeted panel. In addition, tumor-normal DNA was analyzed by whole-genome sequencing (WGS) and tumor RNA was analyzed by RNA sequencing (RNA-seq). The WGS and RNA-seq data were analyzed by a team of bioinformaticians and cancer oncologists, and separately by IBM Watson Genomic Analytics (WGA), an automated system for prioritizing somatic variants and identifying drugs. RESULTS: More variants were identified by WGS/RNA analysis than by targeted panels. WGA completed a comparable analysis in a fraction of the time required by the human analysts. CONCLUSIONS: The development of an effective human-machine interface in the analysis of deep cancer genomic datasets may provide potentially clinically actionable calls for individual patients in a more timely and efficient manner than currently possible. CLINICALTRIALSGOV IDENTIFIER: NCT02725684.

PGBD5 promotes site-specific oncogenic mutations in human tumors.

Genomic rearrangements are a hallmark of human cancers. Here, we identify the piggyBac transposable element derived 5 (PGBD5) gene as encoding an active DNA transposase expressed in the majority of childhood solid tumors, including lethal rhabdoid tumors. Using assembly-based whole-genome DNA sequencing, we found previously undefined genomic rearrangements in human rhabdoid tumors. These rearrangements involved PGBD5-specific signal (PSS) sequences at their breakpoints and recurrently inactivated tumor-suppressor genes. PGBD5 was physically associated with genomic PSS sequences that were also sufficient to mediate PGBD5-induced DNA rearrangements in rhabdoid tumor cells. Ectopic expression of PGBD5 in primary immortalized human cells was sufficient to promote cell transformation in vivo. This activity required specific catalytic residues in the PGBD5 transposase domain as well as end-joining DNA repair and induced structural rearrangements with PSS breakpoints. These results define PGBD5 as an oncogenic mutator and provide a plausible mechanism for site-specific DNA rearrangements in childhood and adult solid tumors.

Next-Generation Rapid Autopsies Enable Tumor Evolution Tracking and Generation of Preclinical Models.

PURPOSE: Patients with cancer who graciously consent for autopsy represent an invaluable resource for the study of cancer biology. To advance the study of tumor evolution, metastases, and resistance to treatment, we developed a next-generation rapid autopsy program integrated within a broader precision medicine clinical trial that interrogates pre- and postmortem tissue samples for patients of all ages and cancer types. MATERIALS AND METHODS: One hundred twenty-three (22%) of 554 patients who consented to the clinical trial also consented for rapid autopsy. This report comprises the first 15 autopsies, including patients with metastatic carcinoma (n = 10), melanoma (n = 1), and glioma (n = 4). Whole-exome sequencing (WES) was performed on frozen autopsy tumor samples from multiple anatomic sites and on non-neoplastic tissue. RNA sequencing (RNA-Seq) was performed on a subset of frozen samples. Tissue was also used for the development of preclinical models, including tumor organoids and patient-derived xenografts. RESULTS: Three hundred forty-six frozen samples were procured in total. WES was performed on 113 samples and RNA-Seq on 72 samples. Successful cell strain, tumor organoid, and/or patient-derived xenograft development was achieved in four samples, including an inoperable pediatric glioma. WES data were used to assess clonal evolution and molecular heterogeneity of tumors in individual patients. Mutational profiles of primary tumors and metastases yielded candidate mediators of metastatic spread and organotropism including CUL9 and PIGM in metastatic ependymoma and ANKRD52 in metastatic melanoma to the lung. RNA-Seq data identified novel gene fusion candidates. CONCLUSION: A next-generation sequencing-based autopsy program in conjunction with a pre-mortem precision medicine pipeline for diverse tumors affords a valuable window into clonal evolution, metastasis, and alterations underlying treatment. Moreover, such an autopsy program yields robust preclinical models of disease.

Whole-Exome Sequencing of Metastatic Cancer and Biomarkers of Treatment Response.

IMPORTANCE: Understanding molecular mechanisms of response and resistance to anticancer therapies requires prospective patient follow-up and clinical and functional validation of both common and low-frequency mutations. We describe a whole-exome sequencing (WES) precision medicine trial focused on patients with advanced cancer. OBJECTIVE: To understand how WES data affect therapeutic decision making in patients with advanced cancer and to identify novel biomarkers of response. DESIGN, SETTING, AND PATIENTS: Patients with metastatic and treatment-resistant cancer were prospectively enrolled at a single academic center for paired metastatic tumor and normal tissue WES during a 19-month period (February 2013 through September 2014). A comprehensive computational pipeline was used to detect point mutations, indels, and copy number alterations. Mutations were categorized as category 1, 2, or 3 on the basis of actionability; clinical reports were generated and discussed in precision tumor board. Patients were observed for 7 to 25 months for correlation of molecular information with clinical response. MAIN OUTCOMES AND MEASURES: Feasibility, use of WES for decision making, and identification of novel biomarkers. RESULTS: A total of 154 tumor-normal pairs from 97 patients with a range of metastatic cancers were sequenced, with a mean coverage of 95X and 16 somatic alterations detected per patient. In total, 16 mutations were category 1 (targeted therapy available), 98 were category 2 (biologically relevant), and 1474 were category 3 (unknown significance). Overall, WES provided informative results in 91 cases (94%), including alterations for which there is an approved drug, there are therapies in clinical or preclinical development, or they are considered drivers and potentially actionable (category 1-2); however, treatment was guided in only 5 patients (5%) on the basis of these recommendations because of access to clinical trials and/or off-label use of drugs. Among unexpected findings, a patient with prostate cancer with exceptional response to treatment was identified who harbored a somatic hemizygous deletion of the DNA repair gene FANCA and putative partial loss of function of the second allele through germline missense variant. Follow-up experiments established that loss of FANCA function was associated with platinum hypersensitivity both in vitro and in patient-derived xenografts, thus providing biologic rationale and functional evidence for his extreme clinical response. CONCLUSIONS AND RELEVANCE: The majority of advanced, treatment-resistant tumors across tumor types harbor biologically informative alterations. The establishment of a clinical trial for WES of metastatic tumors with prospective follow-up of patients can help identify candidate predictive biomarkers of response.