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.

Purpose Prostate cancer translational research has been hampered by the lack of comprehensive and tractable models that represent the genomic landscape of clinical disease. Metastatic castrate-resistant prostate cancer (mCRPC) patient-derived xenografts (PDXs) recapitulate the genetic and phenotypic diversity of the disease. We sought to establish a representative, preclinical platform of PDX-derived organoids that is experimentally facile for high throughput and mechanistic analysis. Experimental Design Using 20 models from the LuCaP mCRPC PDX cohort, including adenocarcinoma and neuroendocrine lineages, we systematically tested > 20 modifications to prostate organoid conditions. Organoids were evaluated for genomic and phenotypic stability and continued reliance on the AR signaling pathway. The utility of the platform as a genotype-dependent model of drug sensitivity was tested with olaparib and carboplatin. Results All PDX models proliferated as organoids in culture. Greater than fifty percent could be continuously cultured long-term in modified conditions; however, none of the PDXs could be established long-term as organoids under previously-reported conditions. Additionally, the modified conditions improved the establishment of patient biopsies over current methods. The genomic heterogeneity of the PDXs, was conserved in organoids. Lineage markers and transcriptomes were maintained between PDXs and organoids. Dependence on AR signaling, was preserved in adenocarcinoma organoids, replicating a dominant characteristic of CRPC. Finally, we observed maximum cytotoxicity to the PARP inhibitor olaparib in BRCA2−/− organoids, similar to responses observed in patients. Conclusions The LuCaP PDX/organoid models provide an expansive, genetically-characterized platform to investigate mechanisms of pathogenesis as well as therapeutic responses and their molecular correlates in mCRPC.

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.

To investigate large structural clonal mosaicism of chromosome X, we analysed the SNP microarray intensity data of 38,303 women from cancer genome-wide association studies (20,878 cases and 17,425 controls) and detected 124 mosaic X events >2 Mb in 97 (0.25%) women. Here we show rates for X-chromosome mosaicism are four times higher than mean autosomal rates; X mosaic events more often include the entire chromosome and participants with X events more likely harbour autosomal mosaic events. X mosaicism frequency increases with age (0.11% in 50-year olds; 0.45% in 75-year olds), as reported for Y and autosomes. Methylation array analyses of 33 women with X mosaicism indicate events preferentially involve the inactive X chromosome. Our results provide further evidence that the sex chromosomes undergo mosaic events more frequently than autosomes, which could have implications for understanding the underlying mechanisms of mosaic events and their possible contribution to risk for chronic diseases.

Early growth response-1 ( EGR1) is a transcription factor correlated with prostate cancer (PC) progression in a variety of contexts. For example, EGR1 levels increase in response to suppressed androgen receptor signaling or loss of the tumor suppressor, PTEN. EGR1 has been shown to regulate genes influencing proliferation, apoptosis, immune cell activation, and matrix degradation, among others. Despite this, the impact of EGR1 on PC metastatic colonization is unclear. We demonstrate using a PC model (DU145/RasB1) of bone and brain metastasis that EGR1 expression regulates angiogenic and osteoclastogenic properties of metastases. We have shown previously that FN14 (TNFRSF12A) and downstream NF-κB signaling is required for metastasis in this model. Here we demonstrate that FN14 ligation also leads to NF-κB-independent, MEK-dependent EGR1 expression. EGR1 -depletion in DU145/RasB1 cells reduced both the number and size of metastases but did not affect primary tumor growth. Decreased EGR1 expression led to reduced blood vessel density in brain and bone metastases as well as decreased osteolytic bone lesion area and reduced numbers of osteoclasts at the bone-tumor interface. TWEAK (TNFSF12) induced several EGR1 -dependent angiogenic and osteoclastogenic factors (e.g. PDGFA, TGFB1, SPP1, IL6, IL8, and TGFA, among others). Consistent with this, in clinical samples of PC, the level of several genes encoding angiogenic/osteoclastogenic pathway effectors correlated with EGR1 levels. Thus, we show here that EGR1 has a direct effect on prostate cancer metastases. EGR1 regulates angiogenic and osteoclastogenic factors, informing the underlying signaling networks that impact autonomous and microenvironmental mechanisms of cancer metastases.

Whether novel risk variants of Alzheimer’s disease (AD) identified through genome-wide association studies (GWAS) also influence MRI-based intermediate phenotypes of AD in the general population is unclear. We studied association of 24 AD risk loci with intracranial volume (ICV), total brain volume (TBV), hippocampal volume (HV), white matter hyperintensity (WMH) burden, and brain infarcts in a meta-analysis of genetic association studies from large population-based samples (N=8,175–11,550). In single-SNP based tests, AD risk allele of APOE (rs2075650) was associated with smaller HV (p=0.0054) and CD33 (rs3865444) with smaller ICV (p=0.0058) In gene-based tests, there was associations of HLA-DRB1 with TBV (p=0.0006) and BIN1 with HV (p=0.00089). A weighted AD genetic risk score was associated with smaller HV (beta±SE=−0.047±0.013, p=0.00041), even after excluding the APOE locus (p=0.029). However, only association of AD genetic risk score with HV, including APOE, was significant after multiple testing correction (including number of independent phenotypes tested). These results suggest that novel AD genetic risk variants may contribute to structural brain aging in non-demented older community persons.

When analyzing the relationship between allelic variability and traits, a potential source of confounding is population admixture. An approach to adjusting for potential confounding due to population admixture when estimating the influence of allelic variability at a candidate gene is presented. The approach involves augmenting linear regression models with additional regressors. Family genotype data are used to define the regressors, and inclusion of the regressors ensures that, even in the presence of population admixture, the estimates of the regression coefficients that parameterize the influence of allelic variability on the trait are unbiased. The approach is illustrated through an analysis of the influence of apolipoprotein E genotype on plasma low density lipoprotein cholesterol concentrations.