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.

Improved detection of low-frequency within-host variants from deep sequencing: A case study with human papillomavirus.

High-coverage sequencing allows the study of variants occurring at low frequencies within samples, but is susceptible to false-positives caused by sequencing error. Ion Torrent has a very low single nucleotide variant (SNV) error rate and has been employed for the majority of human papillomavirus (HPV) whole genome sequences. However, benchmarking of intrahost SNVs (iSNVs) has been challenging, partly due to limitations imposed by the HPV life cycle. We address this problem by deep sequencing three replicates for each of 31 samples of HPV type 18 (HPV18). Errors, defined as iSNVs observed in only one of three replicates, are dominated by C→T (G→A) changes, independently of trinucleotide context. True iSNVs, defined as those observed in all three replicates, instead show a more diverse SNV type distribution, with particularly elevated C→T rates in CCG context (CCG→CTG; CGG→CAG) and C→A rates in ACG context (ACG→AAG; CGT→CTT). Characterization of true iSNVs allowed us to develop two methods for detecting true variants: (1) VCFgenie, a dynamic binomial filtering tool which uses each variant's allele count and coverage instead of fixed frequency cut-offs; and (2) a machine learning binary classifier which trains eXtreme Gradient Boosting models on variant features such as quality and trinucleotide context. Each approach outperforms fixed-cut-off filtering of iSNVs, and performance is enhanced when both are used together. Our results provide improved methods for identifying true iSNVs in within-host applications across sequencing platforms, specifically using HPV18 as a case study.

Somatic mutations in 3929 HPV positive cervical cells associated with infection outcome and HPV type.

Invasive cervical cancers (ICC), caused by HPV infections, have a heterogeneous molecular landscape. We investigate the detection, timing, and HPV type specificity of somatic mutations in 3929 HPV-positive exfoliated cervical cell samples from individuals undergoing cervical screening in the U.S. using deep targeted sequencing in ICC cases, precancers, and HPV-positive controls. We discover a subset of hotspot mutations rare in controls (2.6%) but significantly more prevalent in precancers, particularly glandular precancer lesions (10.2%), and cancers (25.7%), supporting their involvement in ICC carcinogenesis. Hotspot mutations differ by HPV type, and HPV18/45-positive ICC are more likely to have multiple hotspot mutations compared to HPV16-positive ICC. The proportion of cells containing hotspot mutations is higher (i.e., higher variant allele fraction) in ICC and mutations are detectable up to 6 years prior to cancer diagnosis. Our findings demonstrate the feasibility of using exfoliated cervical cells for detection of somatic mutations as potential diagnostic biomarkers.

Differential long-term bivalent HPV vaccine cross-protection by variants in the Costa Rica HPV vaccine trial.

The AS04-adjuvanted human papillomavirus (HPV)16/18 vaccine, an L1-based vaccine, provides strong vaccine efficacy (VE) against vaccine-targeted type infections, and partial cross-protection to phylogenetically-related types, which may be affected by variant-level heterogeneity. We compared VE against incident HPV31, 33, 35, and 45 detections between lineages and SNPs in the L1 region among 2846 HPV-vaccinated and 5465 HPV-unvaccinated women through 11-years of follow-up in the Costa Rica HPV Vaccine Trial. VE was lower against HPV31-lineage-B (VE=60.7%;95%CI = 23.4%,82.8%) compared to HPV31-lineage-A (VE=94.3%;95%CI = 83.7%,100.0%) (VE-ratio = 0.64;95%CI = 0.25,0.90). Differential VE was observed at several lineage-associated HPV31-L1-SNPs, including a nonsynonymous substitution at position 6372 on the FG-loop, an important neutralization domain. For HPV35, the only SNP-level difference was at position 5939 on the DE-loop, with significant VE against nucleotide-G (VE=65.0%;95%CI = 28.0,87.8) but not for more the common nucleotide-A (VE=7.4%;95%CI = -34.1,36.7). Because of the known heterogeneity in precancer/cancer risk across cross-protected HPV genotype variants by race and region, our results of differential variant-level AS04-adjuvanted HPV16/18 vaccine efficacy has global health implications.

Quantitative proteomic analytic approaches to identify metabolic changes in the medial prefrontal cortex of rats exposed to space radiation.

NASA's planned mission to Mars will result in astronauts being exposed to ∼350 mSv/yr of Galactic Cosmic Radiation (GCR). A growing body of data from ground-based experiments indicates that exposure to space radiation doses (approximating those that astronauts will be exposed to on a mission to Mars) impairs a variety of cognitive processes, including cognitive flexibility tasks. Some studies report that 33% of individuals may experience severe cognitive impairment. Translating the results from ground-based rodent studies into tangible risk estimates for astronauts is an enormous challenge, but it would be germane for NASA to use the vast body of data from the rodent studies to start developing appropriate countermeasures, in the expectation that some level of space radiation (SR) -induced cognitive impairment could occur in astronauts. While some targeted studies have reported radiation-induced changes in the neurotransmission properties and/or increased neuroinflammation within space radiation exposed brains, there remains little information that can be used to start the development of a mechanism-based countermeasure strategy. In this study, we have employed a robust label-free mass spectrometry (MS) -based untargeted quantitative proteomic profiling approach to characterize the composition of the medial prefrontal cortex (mPFC) proteome in rats that have been exposed to 15 cGy of 600 MeV/n28Si ions. A variety of analytical techniques were used to mine the generated expression data, which in such studies is typically hampered by low and variable sample size. We have identified several pathways and proteins whose expression alters as a result of space radiation exposure, including decreased mitochondrial function, and a further subset of proteins differs in rats that have a high level of cognitive performance after SR exposure in comparison with those that have low performance levels. While this study has provided further insight into how SR impacts upon neurophysiology, and what adaptive responses can be invoked to prevent the emergence of SR-induced cognitive impairment, the main objective of this paper is to outline strategies that can be used by others to analyze sub-optimal data sets and to identify new information.

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.

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.

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.

Frequency of Pathogenic Germline Variants in Cancer-Susceptibility Genes in the Childhood Cancer Survivor Study.

Background: Pediatric cancers are the leading cause of death by disease in children despite improved survival rates overall. The contribution of germline genetic susceptibility to pediatric cancer survivors has not been extensively characterized. We assessed the frequency of pathogenic or likely pathogenic (P/LP) variants in 5451 long-term pediatric cancer survivors from the Childhood Cancer Survivor Study. Methods: Exome sequencing was conducted on germline DNA from 5451 pediatric cancer survivors (cases who survived ≥5 years from diagnosis; n = 5105 European) and 597 European cancer-free adults (controls). Analyses focused on comparing the frequency of rare P/LP variants in 237 cancer-susceptibility genes and a subset of 60 autosomal dominant high-to-moderate penetrance genes, for both case-case and case-control comparisons. Results: Of European cases, 4.1% harbored a P/LP variant in high-to-moderate penetrance autosomal dominant genes compared with 1.3% in controls (2-sided P = 3 × 10-4). The highest frequency of P/LP variants was in genes typically associated with adult onset rather than pediatric cancers, including BRCA1/2, FH, PALB2, PMS2, and CDKN2A. A statistically significant excess of P/LP variants, after correction for multiple tests, was detected in patients with central nervous system cancers (NF1, SUFU, TSC1, PTCH2), Wilms tumor (WT1, REST), non-Hodgkin lymphoma (PMS2), and soft tissue sarcomas (SDHB, DICER1, TP53, ERCC4, FGFR3) compared with other pediatric cancers. Conclusion: In long-term pediatric cancer survivors, we identified P/LP variants in cancer-susceptibility genes not previously associated with pediatric cancer as well as confirmed known associations. Further characterization of variants in these genes in pediatric cancer will be important to provide optimal genetic counseling for patients and their families.

Germline Mutation in MUS81 Resulting in Impaired Protein Stability is Associated with Familial Breast and Thyroid Cancer.

Multiple primary thyroid cancer (TC) and breast cancer (BC) are commonly diagnosed, and the lifetime risk for these cancers is increased in patients with a positive family history of both TC and BC. Although this phenotype is partially explained by TP53 or PTEN mutations, a significant number of patients are negative for these alterations. We judiciously recruited patients diagnosed with BC and/or TC having a family history of these tumors and assessed their whole-exome sequencing. After variant prioritization, we selected MUS81 c.1292G>A (p.R431H) for further investigation. This variant was genotyped in a healthy population and sporadic BC/TC tissues and investigated at the protein level and cellular models. MUS81 c.1292G>A was the most frequent variant (25%) and the strongest candidate due to its function of double-strand break repair. This variant was confirmed in four relatives from two families. MUS81 p.R431H protein exhibited lower expression levels in tumors from patients positive for the germline variant, compared with wild-type BC, and normal breast and thyroid tissues. Using cell line models, we showed that c.1292G>A induced protein instability and affected DNA damage response. We suggest that MUS81 is a novel candidate involved in familial BC/TC based on its low frequency in healthy individuals and proven effect in protein stability.

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.

Summary from an international cancer seminar focused on human papillomavirus (HPV)-positive oropharynx cancer, convened by scientists at IARC and NCI.

Cancer of the oropharynx has attracted considerable attention in recent years given: (1) an increasing incidence in selected populations over the past three decades; (2) the discovery of human papillomavirus (HPV) infection as the driver of the increase, as opposed to the traditional risk factors such as tobacco (smoking and chewing) and alcohol; and (3) the promise of new prevention and treatment strategies. As a result of such developments, the International Agency for Research on Cancer (IARC) and the US National Cancer Institute (NCI), convened the fourth Cancer Seminar meeting in November 2018 to focus on this topic. This report summarizes the proceedings: a review of recent science on the descriptive epidemiology, etiology, biology, genetics, early detection, pathology and treatment of HPV-positive oropharyngeal cancer, and the formulation of key research questions to be addressed.

Frequency of Pathogenic Germline Variants in Cancer-Susceptibility Genes in Patients With Osteosarcoma.

Importance: Osteosarcoma, the most common malignant bone tumor in children and adolescents, occurs in a high number of cancer predisposition syndromes that are defined by highly penetrant germline mutations. The germline genetic susceptibility to osteosarcoma outside of familial cancer syndromes remains unclear. Objective: To investigate the germline genetic architecture of 1244 patients with osteosarcoma. Design, Setting, and Participants: Whole-exome sequencing (n = 1104) or targeted sequencing (n = 140) of the DNA of 1244 patients with osteosarcoma from 10 participating international centers or studies was conducted from April 21, 2014, to September 1, 2017. The results were compared with the DNA of 1062 individuals without cancer assembled internally from 4 participating studies who underwent comparable whole-exome sequencing and 27 173 individuals of non-Finnish European ancestry who were identified through the Exome Aggregation Consortium (ExAC) database. In the analysis, 238 high-interest cancer-susceptibility genes were assessed followed by testing of the mutational burden across 736 additional candidate genes. Principal component analyses were used to identify 732 European patients with osteosarcoma and 994 European individuals without cancer, with outliers removed for patient-control group comparisons. Patients were subsequently compared with individuals in the ExAC group. All data were analyzed from June 1, 2017, to July 1, 2019. Main Outcomes and Measures: The frequency of rare pathogenic or likely pathogenic genetic variants. Results: Among 1244 patients with osteosarcoma (mean [SD] age at diagnosis, 16 [8.9] years [range, 2-80 years]; 684 patients [55.0%] were male), an analysis restricted to individuals with European ancestry indicated a significantly higher pathogenic or likely pathogenic variant burden in 238 high-interest cancer-susceptibility genes among patients with osteosarcoma compared with the control group (732 vs 994, respectively; P = 1.3 × 10-18). A pathogenic or likely pathogenic cancer-susceptibility gene variant was identified in 281 of 1004 patients with osteosarcoma (28.0%), of which nearly three-quarters had a variant that mapped to an autosomal-dominant gene or a known osteosarcoma-associated cancer predisposition syndrome gene. The frequency of a pathogenic or likely pathogenic cancer-susceptibility gene variant was 128 of 1062 individuals (12.1%) in the control group and 2527 of 27 173 individuals (9.3%) in the ExAC group. A higher than expected frequency of pathogenic or likely pathogenic variants was observed in genes not previously linked to osteosarcoma (eg, CDKN2A, MEN1, VHL, POT1, APC, MSH2, and ATRX) and in the Li-Fraumeni syndrome-associated gene, TP53. Conclusions and Relevance: In this study, approximately one-fourth of patients with osteosarcoma unselected for family history had a highly penetrant germline mutation requiring additional follow-up analysis and possible genetic counseling with cascade testing.

Kinetic resolution of drug intermediates catalyzed by lipase B from Candida antarctica immobilized on immobead-350.

Novozyme 435, which is a commercial immobilized lipase B from Candida antarctica (CALB), has been proven to be inadequate for the kinetic resolution of rac-indanyl acetate. As it has been previously described that different immobilization protocols may greatly alter lipase features, in this work, CALB was covalently immobilized on epoxy Immobead-350 (IB-350) and on glyoxyl-agarose to ascertain if better kinetic resolution would result. Afterwards, all CALB biocatalysts were utilized in the hydrolytic resolution of rac-indanyl acetate and rac-(chloromethyl)-2-(o-methoxyphenoxy) ethyl acetate. After optimization of the immobilization protocol on IB-350, its loading capacity was 150 mg protein/g dried support. Furthermore, the CALB-IB-350 thermal and solvent stabilities were higher than that of the soluble enzyme (e.g., by a 14-fold factor at pH 5-70°C and by a 11-fold factor in dioxane 30%-65°C) and that of the glyoxyl-agarose-CALB (e.g., by a 12-fold factor at pH 10-50°C and by a 21-fold factor in dioxane 30%-65°C). The CALB-IB-350 preparation (with 98% immobilization yield and activity versus p-nitrophenyl butyrate of 6.26 ± 0.2 U/g) was used in the hydrolysis of rac-indanyl acetate using a biocatalyst/substrate ratio of 2:1 and a pH value of 7.0 at 30°C for 24 h. The conversion obtained was 48% and the enantiomeric excess of the product (e.e.p ) was 97%. These values were much higher than the ones obtained with Novozyme 435, 13% and 26% of conversion and e.e.p, respectively. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:878-889, 2018.

Prevalence of pathogenic/likely pathogenic variants in the 24 cancer genes of the ACMG Secondary Findings v2.0 list in a large cancer cohort and ethnicity-matched controls.

BACKGROUND: Prior research has established that the prevalence of pathogenic/likely pathogenic (P/LP) variants across all of the American College of Medical Genetics (ACMG) Secondary Findings (SF) genes is approximately 0.8-5%. We investigated the prevalence of P/LP variants in the 24 ACMG SF v2.0 cancer genes in a family-based cancer research cohort (n = 1173) and in cancer-free ethnicity-matched controls (n = 982). METHODS: We used InterVar to classify variants and subsequently conducted a manual review to further examine variants of unknown significance (VUS). RESULTS: In the 24 genes on the ACMG SF v2.0 list associated with a cancer phenotype, we observed 8 P/LP unique variants (8 individuals; 0.8%) in controls and 11 P/LP unique variants (14 individuals; 1.2%) in cases, a non-significant difference. We reviewed 115 VUS. The median estimated per-variant review time required was 30 min; the first variant within a gene took significantly (p = 0.0009) longer to review (median = 60 min) compared with subsequent variants (median = 30 min). The concordance rate was 83.3% for the variants examined by two reviewers. CONCLUSION: The 115 VUS required database and literature review, a time- and labor-intensive process hampered by the difficulty in interpreting conflicting P/LP determinations. By rigorously investigating the 24 ACMG SF v2.0 cancer genes, our work establishes a benchmark P/LP variant prevalence rate in a familial cancer cohort and controls.

HABP2 p.G534E variant in patients with family history of thyroid and breast cancer.

Familial Papillary Thyroid Carcinoma (PTC) has been described as a hereditary predisposition cancer syndrome associated with mutations in candidate genes including HABP2. Two of 20 probands from families with history of PTC and breast carcinoma (BC) were evaluated by whole exome sequencing (WES) revealing HABP2 p.G534E. Sanger sequencing was used to confirm the involvement of this variant in three families (F1: 7 relatives; F2: 3 and F3: 3). The proband and his sister (with no malignant tumor so far) from F1 were homozygous for the variant whereas one relative with PTC from F2 was negative for the variant. Although the proband of the F3 with PTC was HABP2 wild type, three relatives presented the variant. Five of 170 healthy Brazilian individuals with no family history of BC or PTC and three of 50 sporadic PTC presented the p.G534E. These findings suggested no association of this variant with our familial PTC cases. Genes potentially associated with deregulation of the extracellular matrix organization pathway (CTSB, TNXB, COL4A3, COL16A1, COL24A1, COL5A2, NID1, LOXL2, MMP11, TRIM24 and MUSK) and DNA repair function (NBN and MSH2) were detected by WES, suggesting that other cancer-associated genes have pathogenic effects in the risk of familial PTC development.

Rare germline alterations in cancer-related genes associated with the risk of multiple primary tumor development.

Multiple primary tumors (MPT) have been described in carriers of inherited cancer predisposition genes. However, the genetic etiology of a large proportion of MPT cases remains unclear. We reviewed 267 patients with hereditary cancer predisposition syndromes (HCPS) that underwent genetic counseling and selected 22 patients with MPT to perform genomic analysis (CytoScan HD Array, Affymetrix) aiming to identify new alterations related to a high risk of developing MPT. Twenty patients had a positive family history of cancer and 11 met phenotypic criteria for HCPS. Genetic testing for each of the genes associated with these syndromes revealed negative results for pathogenic mutations. Seventeen rare germline copy number variations (CNVs) covering 40 genes were identified in 11 patients, including an EPCAM/MSH2 deletion in one Lynch syndrome patient. An enrichment analysis revealed a significant number of genes (where the CNVs are mapped) associated with carcinogenesis and/or related to functions implicated with tumor development, such as proliferation and cell survival. An interaction network analysis highlighted the importance of TP53 pathway in cancer emergence. A high number of germline copy-neutral loss of heterozygosity (cnLOH) was identified in nine cases, particularly in two patients. Eighteen genes were covered by both rare CNVs and cnLOH, including 14 related to tumorigenesis and seven genes (ABCC1, KDM4C, KIAA0430, MYH11, NDE1, PIWIL2, and ULK2) specifically associated with cellular growth and proliferation. Overall, we identified 14 cases with rare CNVs and/or cnLOH that may contribute to the risk of MPT development. KEY MESSAGE: CNVs may explain the risk of hereditary cancer syndromes in MPT patients. CNVs affecting genes related to cancer are candidates to be involved in MPT risk. EPCAM/MSH2 deletions should be investigated in patients suspected to have LS. Gene enrichment related to the TP53 network is associated with MPT development. cnLOH and CNVs contribute to the risk of MPT development.

Integrative miRNA and mRNA analysis in penile carcinomas reveals markers and pathways with potential clinical impact.

Penile carcinoma (PeCa) is an important public health issue in poor and developing countries, and has only recently been explored in terms of genetic and epigenetic studies. Integrative data analysis is a powerful method for the identification of molecular drivers involved in cancer development and progression. miRNA and mRNA expression profiles followed by integrative analysis were investigated in 23 PeCa and 12 non-neoplastic penile tissues (NPT). Expression levels of eight miRNAs and 10 mRNAs were evaluated in the same set of samples used for microarray and in a validation set of cases (PeCa = 36; NPT = 27). Eighty-one miRNAs and 2,697 mRNAs were identified as differentially expressed in PeCa. Integrative data analysis revealed 255 mRNAs potentially regulated by 68 miRNAs. Using RT-qPCR, eight miRNAs and nine transcripts were confirmed as altered in PeCa. We identified that MMP1, MMP12 and PPARG and hsa-miR-31-5p, hsa-miR-224-5p, and hsa-miR-223-3p were able to distinguish tumors from NPT with high sensitivity and specificity. Higher MMP1 expression was detected as a better predictor of lymph node metastasis than the clinical-pathological data. In addition, PPARG and EGFR were highlighted as potential pathways for targeted therapy in PeCa. The analysis based on HPV positivity (7 of 23 cases) revealed five miRNA and 13 mRNA differentially expressed. Although in a limited number of cases, HPV positive PeCa presented less aggressive phenotype in comparison with negative cases. Overall, an integrative analysis using mRNA and miRNA profiles revealed markers related with tumor development and progression. Furthermore, MMP1 expression level was a predictive marker for lymph node metastasis in patients with PeCa.

Integrative transcriptome analysis identifies deregulated microRNA-transcription factor networks in lung adenocarcinoma.

Herein, we aimed at identifying global transcriptome microRNA (miRNA) changes and miRNA target genes in lung adenocarcinoma. Samples were selected as training (N = 24) and independent validation (N = 34) sets. Tissues were microdissected to obtain >90% tumor or normal lung cells, subjected to miRNA transcriptome sequencing and TaqMan quantitative PCR validation. We further integrated our data with published miRNA and mRNA expression datasets across 1,491 lung adenocarcinoma and 455 normal lung samples. We identified known and novel, significantly over- and under-expressed (p ≤ 0.01 and FDR≤0.1) miRNAs in lung adenocarcinoma compared to normal lung tissue: let-7a, miR-10a, miR-15b, miR-23b, miR-26a, miR-26b, miR-29a, miR-30e, miR-99a, miR-146b, miR-181b, miR-181c, miR-421, miR-181a, miR-574 and miR-1247. Validated miRNAs included let-7a-2, let-7a-3, miR-15b, miR-21, miR-155 and miR-200b; higher levels of miR-21 expression were associated with lower patient survival (p = 0.042). We identified a regulatory network including miR-15b and miR-155, and transcription factors with prognostic value in lung cancer. Our findings may contribute to the development of treatment strategies in lung adenocarcinoma.