Genomic characterization of cervical lymph node metastases in papillary thyroid carcinoma following the Chornobyl accident.

Childhood radioactive iodine exposure from the Chornobyl accident increased papillary thyroid carcinoma (PTC) risk. While cervical lymph node metastases (cLNM) are well-recognized in pediatric PTC, the PTC metastatic process and potential radiation association are poorly understood. Here, we analyze cLNM occurrence among 428 PTC with genomic landscape analyses and known drivers (131I-exposed = 349, unexposed = 79; mean age = 27.9 years). We show that cLNM are more frequent in PTC with fusion (55%) versus mutation (30%) drivers, although the proportion varies by specific driver gene (RET-fusion = 71%, BRAF-mutation = 38%, RAS-mutation = 5%). cLNM frequency is not associated with other characteristics, including radiation dose. cLNM molecular profiling (N = 47) demonstrates 100% driver concordance with matched primary PTCs and highly concordant mutational spectra. Transcriptome analysis reveals 17 differentially expressed genes, particularly in the HOXC cluster and BRINP3; the strongest differentially expressed microRNA also is near HOXC10. Our findings underscore the critical role of driver alterations and provide promising candidates for elucidating the biological underpinnings of PTC cLNM.

MultipleXed Population Selection and Enrichment single nucleus RNA sequencing (XPoSE-seq) enables sample identity retention during transcriptional profiling of rare populations.

Single nucleus RNA-sequencing is critical in deciphering tissue heterogeneity and identifying rare populations. However, current high throughput techniques are not optimized for rare target populations and require tradeoffs in design due to feasibility. We provide a novel snRNA pipeline, MulipleXed Population Selection and Enrichment snRNA-sequencing (XPoSE-seq), to enable targeted snRNA-seq experiments and in-depth transcriptomic characterization of rare target populations while retaining individual sample identity.

Metagenomic analysis to identify novel infectious agents in systemic anaplastic large cell lymphoma.

Systemic anaplastic large cell lymphoma (ALCL) is a rare CD30-expressing T-cell non-Hodgkin lymphoma. Risk of systemic ALCL is highly increased among immunosuppressed individuals. Because risk of cancers associated with viruses is increased with immunosuppression, we conducted a metagenomic analysis of systemic ALCL to determine whether a known or novel pathogen is associated with this malignancy. Total RNA was extracted and sequenced from formalin-fixed paraffin-embedded tumor specimens from 19 systemic ALCL cases (including one case from an immunosuppressed individual with human immunodeficiency virus infection), 3 Epstein-Barr virus positive diffuse large B-cell lymphomas (DLBCLs) occurring in solid organ transplant recipients (positive controls), and 3 breast cancers (negative controls). We used a pipeline based on the Genome Analysis Toolkit (GATK)-PathSeq algorithm to subtract out human RNA reads and map the remaining RNA reads to microbes. No microbial association with ALCL was identified, but we found Epstein-Barr virus in the DLBCL positive controls and determined the breast cancers to be negative. In conclusion, we did not find a pathogen associated with systemic ALCL, but because we analyzed only one ALCL tumor from an immunosuppressed person, we cannot exclude the possibility that a pathogen is associated with some cases that arise in the setting of immunosuppression.

Phylogenomic Analysis of Human Papillomavirus Type 31 and Cervical Carcinogenesis: A Study of 2093 Viral Genomes.

Human papillomavirus (HPV) type 31 (HPV31) is closely related to the most carcinogenic type, HPV16, but only accounts for 4% of cervical cancer cases worldwide. Viral genetic and epigenetic variations have been associated with carcinogenesis for other high-risk HPV types, but little is known about HPV31. We sequenced 2093 HPV31 viral whole genomes from two large studies, one from the U.S. and one international. In addition, we investigated CpG methylation in a subset of 175 samples. We evaluated the association of HPV31 lineages/sublineages, single nucleotide polymorphisms (SNPs) and viral methylation with cervical carcinogenesis. HPV31 A/B clade was >1.8-fold more associated with cervical intraepithelial neoplasia grade 3 and cancer (CIN3+) compared to the most common C lineage. Lineage/sublineage distribution varied by race/ethnicity and geographic region. A viral genome-wide association analysis identified SNPs within the A/B clade associated with CIN3+, including H23Y (C626T) (odds ratio = 1.60, confidence intervals = 1.17-2.19) located in the pRb CR2 binding-site within the E7 oncogene. Viral CpG methylation was higher in lineage B, compared to the other lineages, and was most elevated in CIN3+. In conclusion, these data support the increased oncogenicity of the A/B lineages and suggest variation of E7 as a contributing risk factor.

APOBEC Mutagenesis Is Concordant between Tumor and Viral Genomes in HPV-Positive Head and Neck Squamous Cell Carcinoma.

APOBEC is a mutagenic source in human papillomavirus (HPV)-mediated malignancies, including HPV+ oropharyngeal squamous cell carcinoma (HPV + OPSCC), and in HPV genomes. It is unknown why APOBEC mutations predominate in HPV + OPSCC, or if the APOBEC-induced mutations observed in both human cancers and HPV genomes are directly linked. We performed sequencing of host somatic exomes, transcriptomes, and HPV16 genomes from 79 HPV + OPSCC samples, quantifying APOBEC mutational burden and activity in both host and virus. APOBEC was the dominant mutational signature in somatic exomes. In viral genomes, there was a mean of five (range 0-29) mutations per genome. The mean of APOBEC mutations in viral genomes was one (range 0-5). Viral APOBEC mutations, compared to non-APOBEC mutations, were more likely to be low-variant allele fraction mutations, suggesting that APOBEC mutagenesis actively occurrs in viral genomes during infection. HPV16 APOBEC-induced mutation patterns in OPSCC were similar to those previously observed in cervical samples. Paired host and viral analyses revealed that APOBEC-enriched tumor samples had higher viral APOBEC mutation rates (p = 0.028), and APOBEC-associated RNA editing (p = 0.008), supporting the concept that APOBEC mutagenesis in host and viral genomes is directly linked and occurrs during infection. Using paired sequencing of host somatic exomes, transcriptomes, and viral genomes, we demonstrated for the first-time definitive evidence of concordance between tumor and viral APOBEC mutagenesis. This finding provides a missing link connecting APOBEC mutagenesis in host and virus and supports a common mechanism driving APOBEC dysregulation.

Pathogenic germline IKZF1 variant alters hematopoietic gene expression profiles.

IKZF1 encodes Ikaros, a zinc finger-containing transcription factor crucial to the development of the hematopoietic system. Germline pathogenic variants in IKZF1 have been reported in patients with acute lymphocytic leukemia and immunodeficiency syndromes. Diamond-Blackfan anemia (DBA) is a rare inherited bone marrow failure syndrome characterized by erythroid hypoplasia, associated with a spectrum of congenital anomalies and an elevated risk of certain cancers. DBA is usually caused by heterozygous pathogenic variants in genes that function in ribosomal biogenesis; however, in many cases the genetic etiology is unknown. We identified a germline IKZF1 variant, rs757907717 C > T, in identical twins with DBA-like features and autoimmune gastrointestinal disease. rs757907717 C > T results in a p.R381C amino acid change in the IKZF1 Ik-x isoform (p.R423C on isoform Ik-1), which we show is associated with altered global gene expression and perturbation of transcriptional networks involved in hematopoietic system development. These data suggest that this missense substitution caused a DBA-like syndrome in this family because of alterations in hematopoiesis, including dysregulation of networks essential for normal erythropoiesis and the immune system.

Oropharyngeal Squamous Cell Carcinoma Morphology and Subtypes by Human Papillomavirus Type and by 16 Lineages and Sublineages.

Oropharyngeal squamous cell carcinoma (SCC) is increasing in incidence and, in Western countries, strongly associated with transcriptionally-active high-risk human papillomavirus (HPV). Within HPV-positive tumors, there is wide morphologic diversity with numerous histologic subtypes of SCC. There are also variable degrees of keratinization, anaplasia, stromal fibrosis, and maturing squamous differentiation. Unlike in the uterine cervix, where associations between HPV types and lineages/sublineages within types have been investigated with some clear correlations identified, little to no data exists for oropharyngeal SCC. In this study, for a large cohort of oropharyngeal SCC patients, we performed RTPCR for high-risk HPV. For the HPV positive patients, we sequenced the DNA of the entire HPV16 genome and determined lineages and sublineages, correlating HPV status, genotype, and HPV16 lineages/sublineages with SCC subtype and various histologic features. Of the 259 patients, 224 (86.5%) were high-risk HPV positive, of which 210/224 (93.8%) were HPV type 16 and 6/224 (2.7%) HPV type 33. Of the four HPV16 lineages, A was the most frequent (192/214 or 89.8%) and of the HPV16 A sublineages, A1 was the most frequent (112/210 or 53.3%). Patients with HPV negative tumors were more often keratinizing vs other types (23/35 or 65.7%) and thus more likely to have more maturing squamous differentiation and stromal desmoplasia. There was no significant correlation between HPV type (16 versus other), between HPV16 lineage (A versus others), or HPV16 A sublineages (A1 or A2 versus others) and morphologic type of SCC nor the various morphologic features of anaplasia/multinucleation, degree of keratinization, nor amount of stromal desmoplasia. In summary, in our cohort, there was no correlation between the type of HPV, the HPV 16 lineage or sublineage, and any of the histologic features or morphologic SCC subtypes.

Radiation-related genomic profile of papillary thyroid carcinoma after the Chernobyl accident.

The 1986 Chernobyl nuclear power plant accident increased papillary thyroid carcinoma (PTC) incidence in surrounding regions, particularly for radioactive iodine (131I)-exposed children. We analyzed genomic, transcriptomic, and epigenomic characteristics of 440 PTCs from Ukraine (from 359 individuals with estimated childhood 131I exposure and 81 unexposed children born after 1986). PTCs displayed radiation dose-dependent enrichment of fusion drivers, nearly all in the mitogen-activated protein kinase pathway, and increases in small deletions and simple/balanced structural variants that were clonal and bore hallmarks of nonhomologous end-joining repair. Radiation-related genomic alterations were more pronounced for individuals who were younger at exposure. Transcriptomic and epigenomic features were strongly associated with driver events but not radiation dose. Our results point to DNA double-strand breaks as early carcinogenic events that subsequently enable PTC growth after environmental radiation exposure.

HPV+ oropharyngeal squamous cell carcinomas from patients with two tumors display synchrony of viral genomes yet discordant mutational profiles and signatures.

Human papillomavirus (HPV) positive oropharyngeal squamous cell carcinoma (HPV + OPSCC) is increasing in prevalence in the USA, as are cases of patients with multiple HPV + OPSCCs (mHPV + OPSCC). mHPV + OPSCCs present a unique opportunity to examine HPV + OPSCC mutation acquisition and evolution. We performed sequencing of the viral genome, somatic exome and somatic transcriptome from 8 patients each with 2 spatially distinct HPV + OPSCCs, and 37 'traditional' HPV + OPSCCs to first address if paired tumors are caused by the same viral isolate and next, if acquired alterations, and the underlying processes driving mutagenesis, are shared within pairs. All tumor pairs contained viral genomes from the same HPV type 16 sublineage and differed by 0-2 clonal single nucleotide polymorphisms (SNPs), suggesting infection with the same viral isolate. Despite this, there was significant discordance in expression profiles, mutational burden and mutational profiles between tumors in a pair, with only two pairs sharing any overlapping mutations (3/3343 variants). Within tumor pairs there was a striking discrepancy of mutational signatures, exemplified by no paired tumors sharing high APOBEC mutational burden. Here, leveraging mHPV + OPSCCs as a model system to study mutation acquisition in virally mediated tumors, in which the germline, environmental exposures, immune surveillance and tissue/organ type were internally controlled, we demonstrate that despite infection by the same viral isolate, paired mHPV + OPSCCs develop drastically different somatic alterations and even more strikingly, appear to be driven by disparate underlying mutational processes. Thus, despite a common starting point, HPV + OPSCCs evolve through variable mutational processes with resultant stochastic mutational profiles.

The D2 and D3 Sublineages of Human Papilloma Virus 16-Positive Cervical Cancer in Guatemala Differ in Integration Rate and Age of Diagnosis.

Human papillomavirus (HPV) 16 displays substantial sequence variation; four HPV16 lineages (A, B, C, and D) have been described as well as multiple sublineages. To identify molecular events associated with HPV16 carcinogenesis, we evaluated viral variation, the integration of HPV16, and somatic mutation in 96 cervical cancer samples from Guatemala. A total of 65% (62/96) of the samples had integrated HPV16 sequences and integration was associated with an earlier age of diagnosis and premenopausal disease. HPV16 integration sites were broadly distributed in the genome, but in one tumor, HPV16 integrated into the promoter of the IFN regulatory factor 4 (IRF4) gene, which plays an important role in the regulation of the IFN response to viral infection. The HPV16 D2 and D3 sublineages were found in 23% and 30% of the tumors, respectively, and were significantly associated with adenocarcinoma. D2-positive tumors had a higher rate of integration, earlier age of diagnosis, and a lower rate of somatic mutation, whereas D3-positive tumors were less likely to integrate, had later age of diagnosis, and exhibited a higher rate of somatic mutation. In conclusion, Guatemalan cervical tumors have a high frequency of very high-risk HPV16 D2 and D3 sublineages harboring distinct histology, which may help guide future therapeutic strategies to target the tumor and reduce recurrence. SIGNIFICANCE: This study details the biological and molecular properties of the most pathogenic forms of HPV16, the cause of the majority of cervical cancers.

Association of HPV35 with cervical carcinogenesis among women of African ancestry: Evidence of viral-host interaction with implications for disease intervention.

HPV35 has been found in only ∼2% of invasive cervical cancers (ICC) worldwide but up to 10% in Sub-Saharan Africa, warranting further investigation and consideration of impact on preventive strategies. We studied HPV35 and ethnicity, in relation to the known steps in cervical carcinogenesis, using multiple large epidemiologic studies in the U.S. and internationally. Combining five U.S. studies, we measured HPV35 positivity and, in Northern California, observed HPV35 type-specific population prevalence and estimated 5-year risk of developing precancer when HPV35-positive. HPV35 genetic variation was examined for differences in carcinogenicity in 1053 HPV35+ cervical specimens from a U.S. cohort and an international collection. African-American women had more HPV35 (12.1% vs 5.1%, P < .001) and more HPV35-associated precancers (7.4% vs 2.1%, P < .001) compared to other ethnicities. Precancer risks after HPV35 infection did not vary by ethnicity (global P = .52). The HPV35 A2 sublineage showed an increased association with precancer/cancer in African-Americans (OR = 5.6 vs A1, 95% CI = 1.3-24.8) and A2 was more prevalent among ICC in Africa than other world regions (41.9% vs 10.4%, P < .01). Our analyses support a strong link between HPV35 and cervical carcinogenesis in women of African ancestry. Current HPV vaccines cover the majority of cervical precancer/cancer across all ethnic groups; additional analyses are required to determine whether the addition of HPV35 to the already highly effective nine-valent HPV vaccine would provide better protection for women in Africa or of African ancestry.

Mutations in the HPV16 genome induced by APOBEC3 are associated with viral clearance.

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.

1q21.1 deletion and a rare functional polymorphism in siblings with thrombocytopenia-absent radius-like phenotypes.

Thrombocytopenia-absent radii (TAR) syndrome, characterized by neonatal thrombocytopenia and bilateral radial aplasia with thumbs present, is typically caused by the inheritance of a 1q21.1 deletion and a single-nucelotide polymorphism in RBM8A on the nondeleted allele. We evaluated two siblings with TAR-like dysmorphology but lacking thrombocytopenia in infancy. Family NCI-107 participated in an IRB-approved cohort study and underwent comprehensive clinical and genomic evaluations, including aCGH, whole-exome, whole-genome, and targeted sequencing. Gene expression assays and electromobility shift assays (EMSAs) were performed to evaluate the variant of interest. The previously identified TAR-associated 1q21.1 deletion was present in the affected siblings and one healthy parent. Multiple sequencing approaches did not identify previously described TAR-associated SNPs or mutations in relevant genes. We discovered rs61746197 A > G heterozygosity in the parent without the deletion and apparent hemizygosity in both siblings. rs61746197 A > G overlaps a RelA-p65 binding motif, and EMSAs indicate the A allele has higher transcription factor binding efficiency than the G allele. Stimulation of K562 cells to induce megakaryocyte differentiation abrogated the shift of both reference and alternative probes. The 1q21.1 TAR-associated deletion in combination with the G variant of rs61746197 on the nondeleted allele is associated with a TAR-like phenotype. rs61746197 G could be a functional enhancer/repressor element, but more studies are required to identify the specific factor(s) responsible. Overall, our findings suggest a role of rs61746197 A > G and human disease in the setting of a 1q21.1 deletion on the other chromosome.

Human papillomavirus 16 sub-lineage dispersal and cervical cancer risk worldwide: Whole viral genome sequences from 7116 HPV16-positive women.

BACKGROUND: Human papillomavirus (HPV)16 can be separated into genetic sub-lineages (A1-4, B1-4, C1-4, D1-4) which may have differential cervical cancer risk. METHODS: A next-generation sequencing assay was used to whole-genome sequence 7116 HPV16-positive cervical samples from well-characterised international epidemiological studies, including 2076 controls, 1878 squamous cell carcinoma (SCC) and 186 adenocarcinoma/adenosquamous cell carcinoma (ADC), and to assign HPV16 sub-lineage. Logistic regression was used to estimate region-stratified country-adjusted odds ratios (OR) and 95%CI. RESULTS: A1 was the most globally widespread sub-lineage, with others showing stronger regional specificity (A3 and A4 for East Asia, B1-4 and C1-4 for Africa, D2 for the Americas, B4, C4 and D4 for North Africa). Increased cancer risks versus A1 were seen for A3, A4 and D (sub)lineages in regions where they were common: A3 in East Asia (OR=2.2, 95%CI:1.0-4.7); A4 in East Asia (6.6, 3.1-14.1) and North America (3.8, 1.7-8.3); and D in North (6.2, 4.1-9.3) and South/Central America (2.2, 0.8-5.7), where D lineages were also more frequent in ADC than SCC (3.2, 1.5-6.5; 12.1, 5.7-25.6, respectively). CONCLUSIONS: HPV16 genetic variation can strongly influence cervical cancer risk. However, burden of cervical cancer attributable to different sub-lineages worldwide is largely driven by historical HPV16 sub-lineage dispersal.

HPV16 E7 Genetic Conservation Is Critical to Carcinogenesis.

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.

A Goal Congruity Model of Role Entry, Engagement, and Exit: Understanding Communal Goal Processes in STEM Gender Gaps.

The goal congruity perspective provides a theoretical framework to understand how motivational processes influence and are influenced by social roles. In particular, we invoke this framework to understand communal goal processes as proximal motivators of decisions to engage in science, technology, engineering, and mathematics (STEM). STEM fields are not perceived as affording communal opportunities to work with or help others, and understanding these perceived goal affordances can inform knowledge about differences between (a) STEM and other career pathways and (b) women's and men's choices. We review the patterning of gender disparities in STEM that leads to a focus on communal goal congruity (Part I), provide evidence for the foundational logic of the perspective (Part II), and explore the implications for research and policy (Part III). Understanding and transmitting the opportunities for communal goal pursuit within STEM can reap widespread benefits for broadening and deepening participation.

Molecular methods resolve the bacterial composition of natural marine biofilms on galvanically coupled stainless steel cathodes.

Navy vessels consist of various metal alloys and biofilm accumulation at the metal surface is thought to play a role in influencing metal deterioration. To develop better strategies to monitor and control metallic biofilms, it is necessary to resolve the bacterial composition within the biofilm. This study aimed to determine if differences in electrochemical current could influence the composition of dominant bacteria in a metallic biofilm, and if so, determine the level of resolution using metagenomic amplicon sequencing. Current was generated by creating galvanic couples between cathodes made from stainless steel and anodes made from carbon steel, aluminum, or copper nickel and exposing them in the Delaware Bay. Stainless steel cathodes (SSCs) coupled to aluminum or carbon steel generated a higher mean current (0.39 mA) than that coupled to copper nickel (0.17 mA). Following 3 months of exposure, the bacterial composition of biofilms collected from the SSCs was determined and compared. Dominant bacterial taxa from the two higher current SSCs were different from that of the low-current SSC as determined by DGGE and verified by Illumina DNA-seq analysis. These results demonstrate that electrochemical current could influence the composition of dominant bacteria in metallic biofilms and that amplicon sequencing is sufficient to complement current methods used to study metallic biofilms in marine environments.

Malleability in communal goals and beliefs influences attraction to stem careers: evidence for a goal congruity perspective.

The goal congruity perspective posits that 2 distinct social cognitions predict attraction to science, technology, engineering, or mathematics (STEM) fields. First, individuals may particularly value communal goals (e.g., working with or helping others), due to either chronic individual differences or the salience of these goals in particular contexts. Second, individuals hold beliefs about the activities that facilitate or impede these goals, or goal affordance stereotypes. Women's tendency to endorse communal goals more highly than do men, along with consensual stereotypes that STEM careers impede communal goals, intersect to produce disinterest in STEM careers. We provide evidence for the foundational predictions that gender differences emerge primarily on communal rather than agentic goals (Studies 1a and 3) and that goal affordance stereotypes reflect beliefs that STEM careers are relatively dissociated from communal goals (Studies 1b and 1c). Most critically, we provide causal evidence that activated communal goals decrease interest in STEM fields (Study 2) and that the potential for a STEM career to afford communal goals elicits greater positivity (Study 3). These studies thus provide a novel demonstration that understanding communal goals and goal affordance stereotypes can lend insight into attitudes toward STEM pursuits.