Simultaneous detection of DNA variation and methylation at HLA class II locus and immune gene promoters using targeted SureSelect Methyl-Sequencing.

The Human Leukocyte Antigen (HLA) locus associates with a variety of complex diseases, particularly autoimmune and inflammatory conditions. The HLA-DR15 haplotype, for example, confers the major risk for developing Multiple Sclerosis in Caucasians, pinpointing an important role in the etiology of this chronic inflammatory disease of the central nervous system. In addition to the protein-coding variants that shape the functional HLA-antigen-T cell interaction, recent studies suggest that the levels of HLA molecule expression, that are epigenetically controlled, also play a role in disease development. However, deciphering the exact molecular mechanisms of the HLA association has been hampered by the tremendous genetic complexity of the locus and a lack of robust approaches to investigate it. Here, we developed a method to specifically enrich the genomic DNA from the HLA class II locus (chr6:32,426,802-34,167,129) and proximal promoters of 2,157 immune-relevant genes, utilizing the Agilent RNA-based SureSelect Methyl-Seq Capture related method, followed by sequencing to detect genetic and epigenetic variation. We demonstrated successful simultaneous detection of the genetic variation and quantification of DNA methylation levels in HLA locus. Moreover, by the detection of differentially methylated positions in promoters of immune-related genes, we identified relevant pathways following stimulation of cells. Taken together, we present a method that can be utilized to study the interplay between genetic variance and epigenetic regulation in the HLA class II region, potentially, in a wide disease context.

Correction: Self-collected and clinician-collected anal swabs show modest agreement for HPV genotyping.

[This corrects the article DOI: 10.1371/journal.pone.0250426.].

Differences in integration frequencies and APOBEC3 profiles of five high-risk HPV types adheres to phylogeny.

Persistent infection with Human Papillomavirus (HPV) is responsible for almost all cases of cervical cancers, and HPV16 and HPV18 associated with the majority of these. These types differ in the proportion of viral minor nucleotide variants (MNVs) caused by APOBEC3 mutagenesis as well as integration frequencies. Whether these traits extend to other types remains uncertain. This study aimed to investigate and compare genomic variability and chromosomal integration in the two phylogenetically distinct Alpha-7 and Alpha-9 clades of carcinogenic HPV types. The TaME-seq protocol was employed to sequence cervical cell samples positive for HPV31, HPV33 or HPV45 and combine these with data from a previous study on HPV16 and HPV18. APOBEC3 mutation signatures were found in Alpha-9 (HPV16/31/33) but not in Alpha-7 (HPV18/45). HPV45 had significantly more MNVs compared to the other types. Alpha-7 had higher integration frequency compared to Alpha-9. An increase in integration frequency with increased diagnostic severity was found for Alpha-7. The results highlight important differences and broaden our understanding of the molecular mechanisms behind cervical cancer induced by high-risk HPV types from the Alpha-7 and Alpha-9 clades.

HPV16 and HPV18 type-specific APOBEC3 and integration profiles in different diagnostic categories of cervical samples.

Human papillomavirus (HPV) 16 and 18 are the most predominant types in cervical cancer. Only a small fraction of HPV infections progress to cancer, indicating that additional factors and genomic events contribute to the carcinogenesis, such as minor nucleotide variation caused by APOBEC3 and chromosomal integration. We analysed intra-host minor nucleotide variants (MNVs) and integration in HPV16 and HPV18 positive cervical samples with different morphology. Samples were sequenced using an HPV whole genome sequencing protocol TaME-seq. A total of 80 HPV16 and 51 HPV18 positive samples passed the sequencing depth criteria of 300× reads, showing the following distribution: non-progressive disease (HPV16 n = 21, HPV18 n = 12); cervical intraepithelial neoplasia (CIN) grade 2 (HPV16 n = 27, HPV18 n = 9); CIN3/adenocarcinoma in situ (AIS) (HPV16 n = 27, HPV18 n = 30); cervical cancer (HPV16 n = 5). Similar numbers of MNVs in HPV16 and HPV18 samples were observed for most viral genes, with the exception of HPV18 E4 with higher numbers across clinical categories. APOBEC3 signatures were observed in HPV16 lesions, while similar mutation patterns were not detected for HPV18. The proportion of samples with integration was 13% for HPV16 and 59% for HPV18 positive samples, with a noticeable portion located within or close to cancer-related genes.

Self-collected and clinician-collected anal swabs show modest agreement for HPV genotyping.

BACKGROUND & AIM: Women with HIV/HPV coinfection and cervical lesions are at increased risk of developing HPV related anal cancer. Self-collection of anal swabs may facilitate HPV molecular testing in anal cancer screening, especially in high-risk groups, and yet it is not adequately studied. We evaluated level of agreement between self-collected anal swabs (SCAS) and clinician-collected anal swabs (CCAS) when used for HPV genotyping. We also described the anal HPV genotype distribution and HIV/HPV coinfection. METHODS: We performed a cross sectional study with participants from a visual-inspection-with-acetic-acid and cervicography (VIAC) clinic, in Harare, Zimbabwe. In a clinic setting, the women aged ≥18 years provided anal swabs in duplicate; first CCAS and then SCAS immediately after. HPV detection and genotyping were performed using next generation amplicon sequencing of a 450bp region of the HPV L1 gene. Level of agreement of HPV genotypes between CCAS and SCAS was calculated using the kappa statistic. McNemar tests were used to evaluate agreement in the proportion of genotypes detected by either method. RESULTS: Three-hundred women provided 600 samples for HPV genotyping. HPV genotypes were detected in 25% of SCAS and in 22% of CCAS. The most common genotypes with CCAS were HPV52, HPV62 and HPV70 and with SCAS were HPV62, HPV44, HPV52, HPV53 and HPV68. Total HPV genotypes detected in CCAS were more than those detected in SCAS, 32 versus 27. The agreement of HPV genotypes between the two methods was 0.55 in kappa value (k). The test of proportions using McNemar gave a Chi-square value of 0.75 (p = 0.39). Multiple HPV infections were detected in 28/75 and 29/67 women for CCAS and SCAS respectively. CONCLUSIONS: SCAS and CCAS anal swabs showed moderate agreement, with no statistically significant difference in the proportion of genotypes detected by either methods. Although the differences between the two methods were not statistically significant, CCAS detected more HPV genotypes than SCAS and more HPV infections were detected in SCAS than in CCAS. Our data suggest that self-collected anal swabs can be used as an alternative to clinician-collected anal swabs for HPV genotyping.

STAT3 Mutation Is Associated with STAT3 Activation in CD30+ ALK- ALCL.

Peripheral T-cell lymphomas (PTCL) are a heterogeneous, and often aggressive group of non-Hodgkin lymphomas. Recent advances in the molecular and genetic characterization of PTCLs have helped to delineate differences and similarities between the various subtypes, and the JAK/STAT pathway has been found to play an important oncogenic role. Here, we aimed to characterize the JAK/STAT pathway in PTCL subtypes and investigate whether the activation of the pathway correlates with the frequency of STAT gene mutations. Patient samples from AITL (n = 30), ALCL (n = 21) and PTCL-NOS (n = 12) cases were sequenced for STAT3, STAT5B, JAK1, JAK3, and RHOA mutations using amplicon sequencing and stained immunohistochemically for pSTAT3, pMAPK, and pAKT. We discovered STAT3 mutations in 13% of AITL, 13% of ALK+ ALCL, 38% of ALK- ALCL and 17% of PTCL-NOS cases. However, no STAT5B mutations were found and JAK mutations were only present in ALK- ALCL (15%). Concurrent mutations were found in all subgroups except ALK+ ALCL where STAT3 mutations were always seen alone. High pY-STAT3 expression was observed especially in AITL and ALCL samples. When studying JAK-STAT pathway mutations, pY-STAT3 expression was highest in PTCLs harboring either JAK1 or STAT3 mutations and CD30+ phenotype representing primarily ALK- ALCLs. Further investigation is needed to elucidate the molecular mechanisms of JAK-STAT pathway activation in PTCL.

Novel TMEM173 Mutation and the Role of Disease Modifying Alleles.

Upon binding to pathogen or self-derived cytosolic nucleic acids cyclic GMP-AMP synthase (cGAS) triggers the production of cGAMP that further activates transmembrane protein STING. Upon activation STING translocates from ER via Golgi to vesicles. Monogenic STING gain-of-function mutations cause early-onset type I interferonopathy, with disease presentation ranging from fatal vasculopathy to mild chilblain lupus. Molecular mechanisms underlying the variable phenotype-genotype correlation are presently unclear. Here, we report a novel gain-of-function G207E STING mutation causing a distinct phenotype with alopecia, photosensitivity, thyroid dysfunction, and features of STING-associated vasculopathy with onset in infancy (SAVI), such as livedo reticularis, skin vasculitis, nasal septum perforation, facial erythema, and bacterial infections. Polymorphism in TMEM173 and IFIH1 showed variable penetrance in the affected family, implying contribution to varying phenotype spectrum. The G207E mutation constitutively activates inflammation-related pathways in vitro, and causes aberrant interferon signature and inflammasome activation in patient PBMCs. Treatment with Janus kinase 1 and 2 (JAK1/2) inhibitor baricitinib was beneficiary for a vasculitic ulcer, induced hair regrowth and improved overall well-being in one patient. Protein-protein interactions propose impaired cellular trafficking of G207E mutant. These findings reveal the molecular landscape of STING and propose common polymorphisms in TMEM173 and IFIH1 as likely modifiers of the phenotype.

HPV16 whole genome minority variants in persistent infections from young Dutch women.

BACKGROUND: Chronic infections by one of the oncogenic human papillomaviruses (HPVs) are responsible for near 5% of the global cancer burden and HPV16 is the type most often found in cancers. HPV genomes display unexpected levels of variation when deep-sequenced. Minor nucleotide variations (MNVs) may reveal HPV genomic instability and HPV-related carcinogenic transformation of host cells. OBJECTIVES: The objective of this study was to investigate HPV16 genome variation at the minor variant level on persisting HPV16 cervical infections from a population of young Dutch women. STUDY DESIGN: 15 HPV16 infections were sequenced using a whole-HPV genome deep sequencing protocol (TaME-seq). One infection was followed over a three-year period, eight were followed over a two-year period, three were followed over a one-year period and three infections had a single sampling point. RESULTS AND CONCLUSIONS: Using a 1% variant frequency cutoff, we find on average 48 MNVs per HPV16 genome and 1717 MNVs in total when sequencing coverage was >100 × . We find the transition mutation T > C to be the most common, in contrast to other studies detecting APOBEC-related C > T mutation profiles in pre-cancerous and cancer samples. Our results suggest that the relative mutagenic footprint of HPV16 genomes may differ between the infections in this study and transforming lesions. In addition, we identify a number of MNVs that have previously been associated with higher incidence of high-grade lesions (CIN3+) in a population study. These findings may provide a starting point for future studies exploring causality between emerging HPV minor genomic variants and cancer development.

Gender-Specific Associations Between Saliva Microbiota and Body Size.

OBJECTIVE: The human intestinal microbiota likely play an important role in the development of overweight and obesity. However, the associations between saliva microbiota and body mass index (BMI) have been sparsely studied. The aim of this study was to identify the associations between saliva microbiota and body size in Finnish children. METHODS: The saliva microbiota of 900 Finnish children, aged 11-14 years with measured height and weight, was characterized using 16S rRNA (V3-V4) sequencing. RESULTS: The core saliva microbiota consisted of 14 genera that were present in more than 95% of the Finnish children. The saliva microbiota profiles were gender-specific with higher alpha-diversity in boys than girls and significant differences between the genders in community composition and abundances. Alpha-diversity differed between normal weight and overweight girls and between normal weight and obese boys. The composition was dissimilar between normal weight and obese girls, but not in boys. The relative abundance profiles differed according to body size. Decrease in commensal saliva bacteria were observed in all the body sizes when compared to normal weight children. Notably, the relative abundance of bacteria related to, Veillonella, Prevotella, Selenomonas, and Streptococcus was reduced in obese children. CONCLUSION: Saliva microbiota diversity and composition were significantly associated with body size and gender in Finnish children. Body size-specific saliva microbiota profiles open new avenues for studying the potential roles of microbiota in weight development and management.

TaME-seq: An efficient sequencing approach for characterisation of HPV genomic variability and chromosomal integration.

HPV genomic variability and chromosomal integration are important in the HPV-induced carcinogenic process. To uncover these genomic events in an HPV infection, we have developed an innovative and cost-effective sequencing approach named TaME-seq (tagmentation-assisted multiplex PCR enrichment sequencing). TaME-seq combines tagmentation and multiplex PCR enrichment for simultaneous analysis of HPV variation and chromosomal integration, and it can also be adapted to other viruses. For method validation, cell lines (n = 4), plasmids (n = 3), and HPV16, 18, 31, 33 and 45 positive clinical samples (n = 21) were analysed. Our results showed deep HPV genome-wide sequencing coverage. Chromosomal integration breakpoints and large deletions were identified in HPV positive cell lines and in one clinical sample. HPV genomic variability was observed in all samples allowing identification of low frequency variants. In contrast to other approaches, TaME-seq proved to be highly efficient in HPV target enrichment, leading to reduced sequencing costs. Comprehensive studies on HPV intra-host variability generated during a persistent infection will improve our understanding of viral carcinogenesis. Efficient identification of both HPV variability and integration sites will be important for the study of HPV evolution and adaptability and may be an important tool for use in cervical cancer diagnostics.

Intra-host sequence variability in human papillomavirus.

Human papillomaviruses (HPVs) co-evolve slowly with the human host and each HPV genotype displays epithelial tropisms. We assessed the evolution of intra HPV genotype variants within samples, and their association to anogenital site, cervical cytology and HIV status. Variability in the L1 gene of 35 HPV genotypes was characterized phylogenetically using maximum likelihood, and portrayed by phenotype. Up to a thousand unique variants were identified within individual samples. In-depth analyses of the most prevalent genotypes, HPV16, HPV18 and HPV52, revealed that the high diversity was dominated by a few abundant variants. This suggests high intra-host mutation rates. Clades of HPV16, HPV18 and HPV52 were associated to anatomical site and HIV co-infection. Particularly, we observed that one HPV16 clade was specific to vaginal cells and one HPV52 clade was specific to anal cells. One major HPV52 clade, present in several samples, was strongly associated with cervical neoplasia. Overall, our data suggest that tissue tropism and HIV immunosuppression are strong shapers of HPV evolution.

Reproducibility and repeatability of six high-throughput 16S rDNA sequencing protocols for microbiota profiling.

Culture-independent molecular techniques and advances in next generation sequencing (NGS) technologies make large-scale epidemiological studies on microbiota feasible. A challenge using NGS is to obtain high reproducibility and repeatability, which is mostly attained through robust amplification. We aimed to assess the reproducibility of saliva microbiota by comparing triplicate samples. The microbiota was produced with simplified in-house 16S amplicon assays taking advantage of large number of barcodes. The assays included primers with Truseq (TS-tailed) or Nextera (NX-tailed) adapters and either with dual index or dual index plus a 6-nt internal index. All amplification protocols produced consistent microbial profiles for the same samples. Although, in our study, reproducibility was highest for the TS-tailed method. Five replicates of a single sample, prepared with the TS-tailed 1-step protocol without internal index sequenced on the HiSeq platform provided high alpha-diversity and low standard deviation (mean Shannon and Inverse Simpson diversity was 3.19 ±â€¯0.097 and 13.56 ±â€¯1.634 respectively). Large-scale profiling of microbiota can consistently be produced by all 16S amplicon assays. The TS-tailed-1S dual index protocol is preferred since it provides repeatable profiles on the HiSeq platform and are less labour intensive.

Driver Gene and Novel Mutations in Asbestos-Exposed Lung Adenocarcinoma and Malignant Mesothelioma Detected by Exome Sequencing.

BACKGROUND: Asbestos is a carcinogen linked to malignant mesothelioma (MM) and lung cancer. Some gene aberrations related to asbestos exposure are recognized, but many associated mutations remain obscure. We performed exome sequencing to determine the association of previously known mutations (driver gene mutations) with asbestos and to identify novel mutations related to asbestos exposure in lung adenocarcinoma (LAC) and MM. METHODS: Exome sequencing was performed on DNA from 47 tumor tissues of MM (21) and LAC (26) patients, 27 of whom had been asbestos-exposed (18 MM, 9 LAC). In addition, 9 normal lung/blood samples of LAC were sequenced. Novel mutations identified from exome data were validated by amplicon-based deep sequencing. Driver gene mutations in BRAF, EGFR, ERBB2, HRAS, KRAS, MET, NRAS, PIK3CA, STK11, and ephrin receptor genes (EPHA1-8, 10 and EPHB1-4, 6) were studied for both LAC and MM, and in BAP1, CUL1, CDKN2A, and NF2 for MM. RESULTS: In asbestos-exposed MM patients, previously non-described NF2 frameshift mutation (one) and BAP1 mutations (four) were detected. Exome data mining revealed some genes potentially associated with asbestos exposure, such as MRPL1 and SDK1. BAP1 and COPG1 mutations were seen exclusively in MM. Pathogenic KRAS mutations were common in LAC patients (42 %), both in non-exposed (n = 5) and exposed patients (n = 6). Pathogenic BRAF mutations were found in two LACs. CONCLUSION: BAP1 mutations occurred in asbestos-exposed MM. MRPL1, SDK1, SEMA5B, and INPP4A could possibly serve as candidate genes for alterations associated with asbestos exposure. KRAS mutations in LAC were not associated with asbestos exposure.

The analysis of clonal diversity and therapy responses using STAT3 mutations as a molecular marker in large granular lymphocytic leukemia.

T-cell large granular lymphocytic leukemia and chronic lymphoproliferative disorder of natural killer cells are intriguing entities between benign and malignant lymphoproliferation. The molecular pathogenesis has partly been uncovered by the recent discovery of somatic activating STAT3 and STAT5b mutations. Here we show that 43% (75/174) of patients with T-cell large granular lymphocytic leukemia and 18% (7/39) with chronic lymphoproliferative disorder of natural killer cells harbor STAT3 mutations when analyzed by quantitative deep amplicon sequencing. Surprisingly, 17% of the STAT3-mutated patients carried multiple STAT3 mutations, which were located in different lymphocyte clones. The size of the mutated clone correlated well with the degree of clonal expansion of the T-cell repertoire analyzed by T-cell receptor beta chain deep sequencing. The analysis of sequential samples suggested that current immunosuppressive therapy is not able to reduce the level of the mutated clone in most cases, thus warranting the search for novel targeted therapies. Our findings imply that the clonal landscape of large granular lymphocytic leukemia is more complex than considered before, and a substantial number of patients have multiple lymphocyte subclones harboring different STAT3 mutations, thus mimicking the situation in acute leukemia.

Mutated ephrin receptor genes in non-small cell lung carcinoma and their occurrence with driver mutations-targeted resequencing study on formalin-fixed, paraffin-embedded tumor material of 81 patients.

Non-small cell lung carcinoma (NSCLC) is the most common subtype of lung cancer. The oncogenic potential of receptor tyrosine kinases (RTKs) is widely known and they are potential targets for tailored therapy. Ephrin receptors (Ephs) form the largest group of RTKs. Nevertheless, Ephs are not widely studied in NSCLC so far. The aim of our study was to investigate novel mutations of Eph genes (EPHA1-8, EPHB1-4, EPHB6) and their association with clinically relevant mutations in BRAF, EML4-ALK, EGFR, INSR, KDR, KRAS, MET, PDGFRA, PDGFRB, PIK3, PTEN, RET, and TP53 in NSCLC patients. Targeted resequencing was conducted on 81 formalin-fixed, paraffin-embedded NSCLC tumor specimens. We analyzed missense and nonsense mutations harbored in the coding regions of the selected genes. We found 18 novel mutations of Ephs in 20% (16 of 81) of the patients. Nearly half of these mutations occurred in the protein kinase domain. The mutations were not mutually exclusive with other clinically relevant mutations. Our study shows that Ephs are frequently mutated in NSCLC patients, and occur together with other known mutations relevant to the pathogenicity of NSCLC.

Individualized systems medicine strategy to tailor treatments for patients with chemorefractory acute myeloid leukemia.

UNLABELLED: We present an individualized systems medicine (ISM) approach to optimize cancer drug therapies one patient at a time. ISM is based on (i) molecular profiling and ex vivo drug sensitivity and resistance testing (DSRT) of patients' cancer cells to 187 oncology drugs, (ii) clinical implementation of therapies predicted to be effective, and (iii) studying consecutive samples from the treated patients to understand the basis of resistance. Here, application of ISM to 28 samples from patients with acute myeloid leukemia (AML) uncovered five major taxonomic drug-response subtypes based on DSRT profiles, some with distinct genomic features (e.g., MLL gene fusions in subgroup IV and FLT3-ITD mutations in subgroup V). Therapy based on DSRT resulted in several clinical responses. After progression under DSRT-guided therapies, AML cells displayed significant clonal evolution and novel genomic changes potentially explaining resistance, whereas ex vivo DSRT data showed resistance to the clinically applied drugs and new vulnerabilities to previously ineffective drugs. SIGNIFICANCE: Here, we demonstrate an ISM strategy to optimize safe and effective personalized cancer therapies for individual patients as well as to understand and predict disease evolution and the next line of therapy. This approach could facilitate systematic drug repositioning of approved targeted drugs as well as help to prioritize and de-risk emerging drugs for clinical testing.

STAT3 mutations indicate the presence of subclinical T-cell clones in a subset of aplastic anemia and myelodysplastic syndrome patients.

Large granular lymphocyte leukemia (LGL) is often associated with immune cytopenias and can cooccur in the context of aplastic anemia (AA) and myelodysplastic syndromes (MDS). We took advantage of the recent description of signal transducer and activator of transcription 3 (STAT3) mutations in LGL clonal expansions to test, using sensitive methods, for the presence of these mutations in a large cohort of 367 MDS and 140 AA cases. STAT3 clones can be found not only in known LGL concomitant cases, but in a small proportion of unsuspected ones (7% AA and 2.5% MDS). In STAT3-mutated AA patients, an interesting trend toward better responses of immunosuppressive therapy and an association with the presence of human leukocyte antigen-DR15 were found. MDSs harboring a STAT3 mutant clone showed a lower degree of bone marrow cellularity and a higher frequency of developing chromosome 7 abnormalities. STAT3-mutant LGL clones may facilitate a persistently dysregulated autoimmune activation, responsible for the primary induction of bone marrow failure in a subset of AA and MDS patients.

Discovery of somatic STAT5b mutations in large granular lymphocytic leukemia.

Large granular lymphocytic (LGL) leukemia is characterized by clonal expansion of cytotoxic T cells or natural killer cells. Recently, somatic mutations in the signal transducer and activator of transcription 3 (STAT3) gene were discovered in 28% to 40% of LGL leukemia patients. By exome and transcriptome sequencing of 2 STAT3 mutation-negative LGL leukemia patients, we identified a recurrent, somatic missense mutation (Y665F) in the Src-like homology 2 domain of the STAT5b gene. Targeted amplicon sequencing of 211 LGL leukemia patients revealed 2 additional patients with STAT5b mutations (N642H), resulting in a total frequency of 2% (4 of 211) of STAT5b mutations across all patients. The Y665F and N642H mutant constructs increased the transcriptional activity of STAT5 and tyrosine (Y694) phosphorylation, which was also observed in patient samples. The clinical course of the disease in patients with the N642H mutation was aggressive and fatal, clearly different from typical LGL leukemia with a relatively favorable outcome. This is the first time somatic STAT5 mutations are discovered in human cancer and further emphasizes the role of STAT family genes in the pathogenesis of LGL leukemia.