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.
Display omitted
•5,570 cervical HPV16 genomes revealed high inter-individual variability•Higher number of non-silent variants in benign HPV16 infections compared with cases•E7 gene was strikingly devoid of genetic variants in precancer and cancer cases•HPV16 rare variants potentially induced by the antiviral activity of human APOBEC3
A genomic survey of thousands of cases of HPV infection in women around the world identifies a conserved sequence in the viral genome that is critical for carcinogenesis.
Of the ~60 human papillomavirus (HPV) genotypes that infect the cervicovaginal epithelium, only 12-13 "high-risk" types are well-established as causing cervical cancer, with HPV16 accounting for over ...half of all cases worldwide. While HPV16 is the most important carcinogenic type, variants of HPV16 can differ in their carcinogenicity by 10-fold or more in epidemiologic studies. Strong genotype-phenotype associations embedded in the small 8-kb HPV16 genome motivate molecular studies to understand the underlying molecular mechanisms. Understanding the mechanisms of HPV genomic findings is complicated by the linkage of HPV genome variants. A panel of experts in various disciplines gathered on 21 November 2016 to discuss the interdisciplinary science of HPV oncogenesis. Here, we summarize the discussion of the complexity of the viral-host interaction and highlight important next steps for selected applied basic laboratory studies guided by epidemiological genomic findings.
HPV16 is a common sexually transmitted infection although few infections lead to cervical precancer/cancer; we cannot distinguish nor mechanistically explain why only certain infections progress. ...HPV16 can be classified into four main evolutionary-derived variant lineages (A, B, C, D) that have been previously suggested to have varying disease risks.
We used a high-throughput HPV16 whole-genome sequencing assay to investigate variant lineage risk among 3215 HPV16-infected women. Using sublineages A1/A2 as the reference, we assessed all variant lineage associations with infection outcome over three or more years of follow-up: 1107 control subjects (<CIN2), 906 CIN2, 1008 CIN3, 69 squamous cell carcinomas (SCC), 85 adenocarcinomas in situ (AIS), and 40 adenocarcinomas. Logistic regression models were used to estimate odds ratios (ORs) and 95% confidence intervals (CIs). All statistical tests were two-sided.
A4 sublineage was associated with an increased risk of cancer, specifically adenocarcinoma (OR = 9.81, 95% CI = 2.02 to 47.69, P = 4.7x10(-03)). Lineage B had a lower risk of CIN3 (OR = 0.51, 95% CI = 0. 28 to 0.91, P = 02) while lineage C showed increased risk (OR = 2.06, 95% CI = 1.09 to 3.89, P = 03). D2/D3 sublineages were strongly associated with an increased risk of CIN3 and cancer, particularly D2 (OR for cancer = 28.48, 95% CI = 9.27 to 87.55, P = 5.0x10(-09)). D2 had the strongest increased risk of glandular lesions, AIS (OR = 29.22, 95% CI = 8.94 to 95.51, P = 2.3x10(-08)), and adenocarcinomas (OR = 137.34, 95% CI = 37.21 to 506.88, P = 1.5x10(-13)). Moreover, the risk of precancer and cancer for specific variant lineages varied by a women's race/ethnicity; those women whose race/ethnicity matched that of the infecting HPV16 variant had an increased risk of CIN3 + (P < 001).
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.
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
), 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.
Abstract
Background
The authors investigated the durability of vaccine efficacy (VE) against human papillomavirus (HPV)16 or 18 infections and antibody response among nonrandomly assigned women who ...received a single dose of the bivalent HPV vaccine compared with women who received multiple doses and unvaccinated women.
Methods
HPV infections were compared between HPV16 or 18-vaccinated women aged 18 to 25 years who received one (N = 112), two (N = 62), or three (N = 1365) doses, and age- and geography-matched unvaccinated women (N = 1783) in the long-term follow-up of the Costa Rica HPV Vaccine Trial. Cervical HPV infections were measured at two study visits, approximately 9 and 11 years after initial HPV vaccination, using National Cancer Institute next-generation sequencing TypeSeq1 assay. VE and 95% confidence intervals (CIs) were estimated. HPV16 or 18 antibody levels were measured in all one- and two-dose women, and a subset of three-dose women, using a virus-like particle-based enzyme-linked immunosorbent assay (n = 448).
Results
Median follow-up for the HPV-vaccinated group was 11.3 years (interquartile range = 10.9–11.7 years) and did not vary by dose group. VE against prevalent HPV16 or 18 infection was 80.2% (95% CI = 70.7% to 87.0%) among three-dose, 83.8% (95% CI = 19.5% to 99.2%) among two-dose, and 82.1% (95% CI = 40.2% to 97.0%) among single-dose women. HPV16 or 18 antibody levels did not qualitatively decline between years four and 11 regardless of the number of doses given, although one-dose titers continue to be statistically significantly lower compared with two- and three-dose titers.
Conclusion
More than a decade after HPV vaccination, single-dose VE against HPV16 or 18 infection remained high and HPV16 or 18 antibodies remained stable. A single dose of bivalent HPV vaccine may induce sufficiently durable protection that obviates the need for more doses.
Single-cell RNA sequencing (scRNA-seq) is a powerful technique used to explore gene expression at the single cell level. However, appropriate preparation of samples is essential to obtain the most ...information out of this transformative technology. Generating high-quality single-cell suspensions from the retina is critical to preserve the native expression profile that will ensure meaningful transcriptome data analysis.
We modified the conditions for rapid and optimal dissociation of retina sample preparation. We also included additional filtering steps in data analysis for retinal scRNA-seq.
We report a gentle method for dissociation of the mouse retina that minimizes cell death and preserves cell morphology. This protocol also results in detection of higher transcriptional complexity. In addition, the modified computational pipeline leads to better-quality single-cell RNA-sequencing data in retina samples. We also demonstrate the advantages and limitations of using fresh versus frozen retinas to prepare cell or nuclei suspensions for scRNA-seq.
We provide a simple yet robust and reproducible protocol for retinal scRNA-seq analysis, especially for comparative studies.
Multimers of the HPV genome are generated in cervical tumors replicating as extrachromosomal episomes, which is associated with deletion and rearrangement of the HPV genome and provides a mechanism ...for oncogenesis without integration.
Lung 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.
We 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.
These 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 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. In addition, 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.
BackgroundDiamond-Blackfan anaemia (DBA) is an inherited bone marrow failure syndrome (IBMFS) characterised by erythroid hypoplasia. It is associated with congenital anomalies and a high risk of ...developing specific cancers. DBA is caused predominantly by autosomal dominant pathogenic variants in at least 15 genes affecting ribosomal biogenesis and function. Two X-linked recessive genes have been identified.ObjectivesWe aim to identify the genetic aetiology of DBA.MethodsOf 87 families with DBA enrolled in an institutional review board-approved cohort study (ClinicalTrials.gov Identifier:NCT00027274), 61 had genetic testing information available. Thirty-five families did not have a known genetic cause and thus underwent comprehensive genomic evaluation with whole exome sequencing, deletion and CNV analyses to identify their disease-associated pathogenic variant. Controls for functional studies were healthy mutation-negative individuals enrolled in the same study.ResultsOur analyses uncovered heterozygous pathogenic variants in two previously undescribed genes in two families. One family had a non-synonymous variant (p.K77N) in RPL35; the second family had a non-synonymous variant (p. L51S) in RPL18. Both of these variants result in pre-rRNA processing defects. We identified heterozygous pathogenic variants in previously known DBA genes in 16 of 35 families. Seventeen families who underwent genetic analyses are yet to have a genetic cause of disease identified.ConclusionsOverall, heterozygous pathogenic variants in ribosomal genes were identified in 44 of the 61 families (7 2%). De novo pathogenic variants were observed in 57% of patients with DBA. Ongoing studies of DBA genomics will be important to understand this complex disorder.