Characterization and Diversity of 243 Complete Human Papillomavirus Genomes in Cervical Swabs Using Next Generation Sequencing.

In recent years, next generation sequencing (NGS) technology has been widely used for the discovery of novel human papillomavirus (HPV) genotypes, variant characterization and genotyping. Here, we compared the analytical performance of NGS with a commercial PCR-based assay (Anyplex II HPV28) in cervical samples of 744 women. Overall, HPV positivity was 50.2% by the Anyplex and 45.5% by the NGS. With the NGS, we detected 25 genotypes covered by Anyplex and 41 additional genotypes. Agreement between the two methods for HPV positivity was 80.8% (kappa = 0.616) and 84.8% (kappa = 0.652) for 28 HPV genotypes and 14 high-risk genotypes, respectively. We recovered and characterized 243 complete HPV genomes from 153 samples spanning 40 different genotypes. According to phylogenetic analysis and pairwise distance, we identified novel lineages and sublineages of four high-risk and 16 low-risk genotypes. In total, 17 novel lineages and 14 novel sublineages were proposed, including novel lineages of HPV45, HPV52, HPV66 and a novel sublineage of HPV59. Our study provides important genomic insights on HPV types and lineages, where few complete genomes were publicly available.

Deconvolution of transcriptomes and miRNomes by independent component analysis provides insights into biological processes and clinical outcomes of melanoma patients.

BACKGROUND: The amount of publicly available cancer-related "omics" data is constantly growing and can potentially be used to gain insights into the tumour biology of new cancer patients, their diagnosis and suitable treatment options. However, the integration of different datasets is not straightforward and requires specialized approaches to deal with heterogeneity at technical and biological levels. METHODS: Here we present a method that can overcome technical biases, predict clinically relevant outcomes and identify tumour-related biological processes in patients using previously collected large discovery datasets. The approach is based on independent component analysis (ICA) - an unsupervised method of signal deconvolution. We developed parallel consensus ICA that robustly decomposes transcriptomics datasets into expression profiles with minimal mutual dependency. RESULTS: By applying the method to a small cohort of primary melanoma and control samples combined with a large discovery melanoma dataset, we demonstrate that our method distinguishes cell-type specific signals from technical biases and allows to predict clinically relevant patient characteristics. We showed the potential of the method to predict cancer subtypes and estimate the activity of key tumour-related processes such as immune response, angiogenesis and cell proliferation. ICA-based risk score was proposed and its connection to patient survival was validated with an independent cohort of patients. Additionally, through integration of components identified for mRNA and miRNA data, the proposed method helped deducing biological functions of miRNAs, which would otherwise not be possible. CONCLUSIONS: We present a method that can be used to map new transcriptomic data from cancer patient samples onto large discovery datasets. The method corrects technical biases, helps characterizing activity of biological processes or cell types in the new samples and provides the prognosis of patient survival.

A new ALK isoform transported by extracellular vesicles confers drug resistance to melanoma cells.

BACKGROUND: Drug resistance remains an unsolved clinical issue in oncology. Despite promising initial responses obtained with BRAF and MEK kinase inhibitors, resistance to treatment develops within months in virtually all melanoma patients. METHODS: Microarray analyses were performed in BRAF inhibitor-sensitive and resistant cell lines to identify changes in the transcriptome that might play a role in resistance. siRNA approaches and kinase inhibitors were used to assess the involvement of the identified Anaplastic Lymphoma Kinase (ALK) in drug resistance. The capability of extracellular vesicles (EVs) to transfer drug resistant properties was investigated in co-culture assays. RESULTS: Here, we report a new mechanism of acquired drug resistance involving the activation of a novel truncated form of ALK. Knock down or inhibition of ALK re-sensitised resistant cells to BRAF inhibition and induced apoptosis. Interestingly, truncated ALK was also secreted into EVs and we show that EVs were the vehicle for transferring drug resistance. CONCLUSIONS: To our knowledge, this is the first report demonstrating the functional involvement of EVs in melanoma drug resistance by transporting a truncated but functional form of ALK, able to activate the MAPK signalling pathway in target cells. Combined inhibition of ALK and BRAF dramatically reduced tumour growth in vivo. These findings make ALK a promising clinical target in melanoma patients.

Functionally Convergent B Cell Receptor Sequences in Transgenic Rats Expressing a Human B Cell Repertoire in Response to Tetanus Toxoid and Measles Antigens.

The identification and tracking of antigen-specific immunoglobulin (Ig) sequences within total Ig repertoires is central to high-throughput sequencing (HTS) studies of infections or vaccinations. In this context, public Ig sequences shared by different individuals exposed to the same antigen could be valuable markers for tracing back infections, measuring vaccine immunogenicity, and perhaps ultimately allow the reconstruction of the immunological history of an individual. Here, we immunized groups of transgenic rats expressing human Ig against tetanus toxoid (TT), Modified Vaccinia virus Ankara (MVA), measles virus hemagglutinin and fusion proteins expressed on MVA, and the environmental carcinogen benzo[a]pyrene, coupled to TT. We showed that these antigens impose a selective pressure causing the Ig heavy chain (IgH) repertoires of the rats to converge toward the expression of antibodies with highly similar IgH CDR3 amino acid sequences. We present a computational approach, similar to differential gene expression analysis, that selects for clusters of CDR3s with 80% similarity, significantly overrepresented within the different groups of immunized rats. These IgH clusters represent antigen-induced IgH signatures exhibiting stereotypic amino acid patterns including previously described TT- and measles-specific IgH sequences. Our data suggest that with the presented methodology, transgenic Ig rats can be utilized as a model to identify antigen-induced, human IgH signatures to a variety of different antigens.

The Effects of Sequence Variation on Genome-wide NRF2 Binding--New Target Genes and Regulatory SNPs.

Transcription factor binding specificity is crucial for proper target gene regulation. Motif discovery algorithms identify the main features of the binding patterns, but the accuracy on the lower affinity sites is often poor. Nuclear factor E2-related factor 2 (NRF2) is a ubiquitous redox-activated transcription factor having a key protective role against endogenous and exogenous oxidant and electrophile stress. Herein, we decipher the effects of sequence variation on the DNA binding sequence of NRF2, in order to identify both genome-wide binding sites for NRF2 and disease-associated regulatory SNPs (rSNPs) with drastic effects on NRF2 binding. Interactions between NRF2 and DNA were studied using molecular modelling, and NRF2 chromatin immunoprecipitation-sequence datasets together with protein binding microarray measurements were utilized to study binding sequence variation in detail. The binding model thus generated was used to identify genome-wide binding sites for NRF2, and genomic binding sites with rSNPs that have strong effects on NRF2 binding and reside on active regulatory elements in human cells. As a proof of concept, miR-126-3p and -5p were identified as NRF2 target microRNAs, and a rSNP (rs113067944) residing on NRF2 target gene (Ferritin, light polypeptide, FTL) promoter was experimentally verified to decrease NRF2 binding and result in decreased transcriptional activity.

Interferon-γ-induced activation of Signal Transducer and Activator of Transcription 1 (STAT1) up-regulates the tumor suppressing microRNA-29 family in melanoma cells.

BACKGROUND: The type-II-cytokine IFN-γ is a pivotal player in innate immune responses but also assumes functions in controlling tumor cell growth by orchestrating cellular responses against neoplastic cells. The role of IFN-γ in melanoma is not fully understood: it is a well-known growth inhibitor of melanoma cells in vitro. On the other hand, IFN-γ may also facilitate melanoma progression. While interferon-regulated genes encoding proteins have been intensively studied since decades, the contribution of miRNAs to effects mediated by interferons is an emerging area of research.We recently described a distinct and dynamic regulation of a whole panel of microRNAs (miRNAs) after IFN-γ-stimulation. The aim of this study was to analyze the transcriptional regulation of miR-29 family members in detail, identify potential interesting target genes and thus further elucidate a potential signaling pathway IFN-γ → Jak→ P-STAT1 → miR-29 → miR-29 target genes and its implication for melanoma growth. RESULTS: Here we show that IFN-γ induces STAT1-dependently a profound up-regulation of the miR-29 primary cluster pri-29a~b-1 in melanoma cell lines. Furthermore, expression levels of pri-29a~b-1 and mature miR-29a and miR-29b were elevated while the pri-29b-2~c cluster was almost undetectable. We observed an inverse correlation between miR-29a/b expression and the proliferation rate of various melanoma cell lines. This finding could be corroborated in cells transfected with either miR-29 mimics or inhibitors. The IFN-γ-induced G1-arrest of melanoma cells involves down-regulation of CDK6, which we proved to be a direct target of miR-29 in these cells. Compared to nevi and normal skin, and metastatic melanoma samples, miR-29a and miR-29b levels were found strikingly elevated in certain patient samples derived from primary melanoma. CONCLUSIONS: Our findings reveal that the miR-29a/b1 cluster is to be included in the group of IFN- and STAT-regulated genes. The up-regulated miR-29 family members may act as effectors of cytokine signalling in melanoma and other cancer cells as well as in the immune system.

Dynamic regulation of microRNA expression following interferon-γ-induced gene transcription.

MicroRNAs are major players in post-transcriptional gene regulation. Even small changes in miRNA levels may have profound consequences for the expression levels of target genes. Hence, miRNAs themselves need to be tightly, albeit dynamically, regulated. Here, we investigated the dynamic behavior of miRNAs over a wide time range following stimulation of melanoma cells with interferon-γ (IFN-γ), which activates the transcription factor STAT1. By applying several bioinformatic and statistical software tools for visualization and identification of differentially expressed miRNAs derived from time-series microarray experiments, 8.9% of 1105 miRNAs appeared to be directly or indirectly regulated by STAT1. Focusing on distinct dynamic expression patterns, we found that the majority of robust miRNA expression changes occurred in the intermediate time range (24-48 h). Three miRNAs (miR-27a, miR-30a, miR-34a) had a delayed regulation occurring at 72 h while none showed significant expression changes at early time points between 30 min and 6 h. Expression patterns of individual miRNAs were altered gradually over time or abruptly increased or decreased between two time points. Furthermore, we observed coordinated dynamic transcription of most miRNA clusters while few were found to be regulated independently of their genetic cluster. Most interestingly, several "star" or passenger strand sequences were specifically regulated over time while their "guide" strands were not.

Genome-wide landscape of liver X receptor chromatin binding and gene regulation in human macrophages.

BACKGROUND: The liver X receptors (LXRs) are oxysterol sensing nuclear receptors with multiple effects on metabolism and immune cells. However, the complete genome-wide cistrome of LXR in cells of human origin has not yet been provided. RESULTS: We performed ChIP-seq in phorbol myristate acetate-differentiated THP-1 cells (macrophage-type) after stimulation with the potent synthetic LXR ligand T0901317 (T09). Microarray gene expression analysis was performed in the same cellular model. We identified 1357 genome-wide LXR locations (FDR < 1%), of which 526 were observed after T09 treatment. De novo analysis of LXR binding sequences identified a DR4-type element as the major motif. On mRNA level T09 up-regulated 1258 genes and repressed 455 genes. Our results show that LXR actions are focused on 112 genomic regions that contain up to 11 T09 target genes per region under the control of highly stringent LXR binding sites with individual constellations for each region. We could confirm that LXR controls lipid metabolism and transport and observed a strong association with apoptosis-related functions. CONCLUSIONS: This first report on genome-wide binding of LXR in a human cell line provides new insights into the transcriptional network of LXR and its target genes with their link to physiological processes, such as apoptosis.The gene expression microarray and sequence data have been submitted collectively to the NCBI Gene Expression Omnibus http://www.ncbi.nlm.nih.gov/geo under accession number GSE28319.