SUMMARY Although most cervical human papillomavirus type 16 (HPV16) infections become undetectable within 1–2 years, persistent HPV16 causes half of all cervical cancers. We used a novel HPV whole-genome sequencing technique to evaluate an exceptionally large collection of 5,570 HPV16-infected case-control samples to determine whether viral genetic variation influences risk of cervical precancer and cancer. We observed thousands of unique HPV16 genomes; very few women shared the identical HPV16 sequence, which should stimulate a careful re-evaluation of the clinical implications of HPV mutation rates, transmission, clearance, and persistence. In case-control analyses, HPV16 in the controls had significantly more amino acid changing variants throughout the genome. Strikingly, E7 was devoid of variants in precancers/cancers compared to higher levels in the controls; we confirmed this in cancers from around the world. Strict conservation of the 98 amino acids of E7, which disrupts Rb function, is critical for HPV16 carcinogenesis, presenting a highly specific target for etiologic and therapeutic research.

Specific HPV16 variant sublineages strongly influence risk of histologic types of precancer and cancer, and viral genetic variation may help explain its unique carcinogenic properties.

Dyskeratosis congenita (DC) is an inherited bone marrow failure and cancer predisposition syndrome caused by aberrant telomere biology. The classic triad of dysplastic nails, abnormal skin pigmentation, and oral leukoplakia is diagnostic of DC, but substantial clinical heterogeneity exists; the clinically severe variant Hoyeraal Hreidarsson syndrome (HH) also includes cerebellar hypoplasia, severe immunodeficiency, enteropathy, and intrauterine growth retardation. Germline mutations in telomere biology genes account for approximately one-half of known DC families. Using exome sequencing, we identified mutations in RTEL1, a helicase with critical telomeric functions, in two families with HH. In the first family, two siblings with HH and very short telomeres inherited a premature stop codon from their mother who has short telomeres. The proband from the second family has HH and inherited a premature stop codon in RTEL1 from his father and a missense mutation from his mother, who also has short telomeres. Additionally, inheritance of only the missense mutation led to very short telomeres in the proband’s brother. Targeted sequencing identified a different RTEL1 missense mutation in one additional DC proband who has bone marrow failure and short telomeres. Both missense mutations affect the helicase domain of RTEL1, and three in silico prediction algorithms suggest that they are likely deleterious. The nonsense mutations both cause truncation of the RTEL1 protein, resulting in loss of the PIP box; this may abrogate an important protein-protein interaction. These findings implicate a new telomere biology gene, RTEL1, in the etiology of DC.

BackgroundLung adenocarcinoma (LUAD) is the most common histologic subtype of lung cancer and has a high risk of distant metastasis at every disease stage. We aimed to characterize the genomic landscape of LUAD and identify mutation signatures associated with tumor progression.Methods and FindingsWe performed an integrative genomic analysis, incorporating whole exome sequencing (WES), determination of DNA copy number and DNA methylation, and transcriptome sequencing for 101 LUAD samples from the Environment And Genetics in Lung cancer Etiology (EAGLE) study. We detected driver genes by testing whether the nonsynonymous mutation rate was significantly higher than the background mutation rate and replicated our findings in public datasets with 724 samples. We performed subclonality analysis for mutations based on mutant allele data and copy number alteration data. We also tested the association between mutation signatures and clinical outcomes, including distant metastasis, survival, and tumor grade. We identified and replicated two novel candidate driver genes, POU class 4 homeobox 2 (POU4F2) (mutated in 9 [8.9%] samples) and ZKSCAN1 (mutated in 6 [5.9%] samples), and characterized their major deleterious mutations. ZKSCAN1 was part of a mutually exclusive gene set that included the RTK/RAS/RAF pathway genes BRAF, EGFR, KRAS, MET, and NF1, indicating an important driver role for this gene. Moreover, we observed strong associations between methylation in specific genomic regions and somatic mutation patterns. In the tumor evolution analysis, four driver genes had a significantly lower fraction of subclonal mutations (FSM), including TP53 (p = 0.007), KEAP1 (p = 0.012), STK11 (p = 0.0076), and EGFR (p = 0.0078), suggesting a tumor initiation role for these genes. Subclonal mutations were significantly enriched in APOBEC-related signatures (p < 2.5×10−50). The total number of somatic mutations (p = 0.0039) and the fraction of transitions (p = 5.5×10−4) were associated with increased risk of distant metastasis. Our study’s limitations include a small number of LUAD patients for subgroup analyses and a single-sample design for investigation of subclonality.ConclusionsThese data provide a genomic characterization of LUAD pathogenesis and progression. The distinct clonal and subclonal mutation signatures suggest possible diverse carcinogenesis pathways for endogenous and exogenous exposures, and may serve as a foundation for more effective treatments for this lethal disease. LUAD’s high heterogeneity emphasizes the need to further study this tumor type and to associate genomic findings with clinical outcomes.

Dyskeratosis congenita (DC) is a heterogeneous inherited bone marrow failure and cancer predisposition syndrome in which germline mutations in telomere biology genes account for approximately one-half of known families. Hoyeraal Hreidarsson syndrome (HH) is a clinically severe variant of DC in which patients also have cerebellar hypoplasia and may present with severe immunodeficiency and enteropathy. We discovered a germline autosomal recessive mutation in RTEL1, a helicase with critical telomeric functions, in two unrelated families of Ashkenazi Jewish (AJ) ancestry. The affected individuals in these families are homozygous for the same mutation, R1264H, which affects three isoforms of RTEL1. Each parent was a heterozygous carrier of one mutant allele. Patient-derived cell lines revealed evidence of telomere dysfunction, including significantly decreased telomere length, telomere length heterogeneity, and the presence of extra-chromosomal circular telomeric DNA. In addition, RTEL1 mutant cells exhibited enhanced sensitivity to the interstrand cross-linking agent mitomycin C. The molecular data and the patterns of inheritance are consistent with a hypomorphic mutation in RTEL1 as the underlying basis of the clinical and cellular phenotypes. This study further implicates RTEL1 in the etiology of DC/HH and immunodeficiency, and identifies the first known homozygous autosomal recessive disease-associated mutation in RTEL1.

For unknown reasons, there is huge variability in risk conferred by different HPV types and, remarkably, strong differences even between closely related variant lineages within each type. HPV16 is a uniquely powerful carcinogenic type, causing approximately half of cervical cancer and most other HPV-related cancers. To permit the large-scale study of HPV genome variability and precancer/cancer, starting with HPV16 and cervical cancer, we developed a high-throughput next-generation sequencing (NGS) whole-genome method. We designed a custom HPV16 AmpliSeq™ panel that generated 47 overlapping amplicons covering 99% of the genome sequenced on the Ion Torrent Proton platform. After validating with Sanger, the current “gold standard” of sequencing, in 89 specimens with concordance of 99.9%, we used our NGS method and custom annotation pipeline to sequence 796 HPV16-positive exfoliated cervical cell specimens. The median completion rate per sample was 98.0%. Our method enabled us to discover novel SNPs, large contiguous deletions suggestive of viral integration (OR of 27.3, 95% CI 3.3–222, P=0.002), and the sensitive detection of variant lineage coinfections. This method represents an innovative high-throughput, ultra-deep coverage technique for HPV genomic sequencing, which, in turn, enables the investigation of the role of genetic variation in HPV epidemiology and carcinogenesis.

HPV16 causes half of cervical cancers worldwide; for unknown reasons, most infections resolve within two years. Here, we analyze the viral genomes of 5,328 HPV16-positive case-control samples to investigate mutational signatures and the role of human APOBEC3-induced mutations in viral clearance and cervical carcinogenesis. We identify four de novo mutational signatures, one of which matches the COSMIC APOBEC-associated signature 2. The viral genomes of the precancer/cancer cases are less likely to contain within-host somatic HPV16 APOBEC3-induced mutations (Fisher’s exact test, P = 6.2 x 10−14), and have a 30% lower nonsynonymous APOBEC3 mutation burden compared to controls. We replicate the low prevalence of HPV16 APOBEC3-induced mutations in 1,749 additional cases. APOBEC3 mutations also historically contribute to the evolution of HPV16 lineages. We demonstrate that cervical infections with a greater burden of somatic HPV16 APOBEC3-induced mutations are more likely to be benign or subsequently clear, suggesting they may reduce persistence, and thus progression, within the host.

Metastasis is the leading cause of death in osteosarcoma patients, the most common pediatric bone malignancy. We conducted a multi-stage genome-wide association study of osteosarcoma metastasis at diagnosis in 935 osteosarcoma patients to determine whether germline genetic variation contributes to risk of metastasis. We identified a SNP, rs7034162, in NFIB significantly associated with metastasis in European osteosarcoma cases, as well as in cases of African and Brazilian ancestry (meta-analysis of all cases: P=1.2×10−9, OR 2.43, 95% CI 1.83–3.24). The risk allele was significantly associated with lowered NFIB expression, which led to increased osteosarcoma cell migration, proliferation, and colony formation. Additionally, a transposon screen in mice identified a significant proportion of osteosarcomas harboring inactivating insertions in Nfib, and had lowered Nfib expression. These data suggest that germline genetic variation at rs7034162 is important in osteosarcoma metastasis, and that NFIB is an osteosarcoma metastasis susceptibility gene.