Supplementary Data from Association Study between Polymorphisms in DNA Methylation–Related Genes and Testicular Germ Cell Tumor Risk

Supplementary Table 1 from Analysis of Gene Expression Profiles of Microdissected Cell Populations Indicates that Testicular Carcinoma In situ Is an Arrested Gonocyte

Immunoglobulin A (IgA) is acknowledged to play a role in the defence of the mucosal barrier by coating microorganisms. Surprisingly, IgA-deficient humans exhibit few infection-related complications, raising the question if the more specific IgG may help IgM in compensating for the lack of IgA. Here we employ a cohort of IgA-deficient humans, each paired with IgA-sufficient household members, to investigate multi-Ig bacterial coating. In IgA-deficient humans, IgM alone, and together with IgG, recapitulate coating of most bacterial families, despite an overall 3.6-fold lower Ig-coating. Bacterial IgG coating is dominated by IgG1 and IgG4. Single-IgG2 bacterial coating is sparse and linked to enhanced Escherichia coli load and TNF-α. Although single-IgG2 coating is 1.6-fold more prevalent in IgA deficiency than in healthy controls, it is 2-fold less prevalent than in inflammatory bowel disease. Altogether we demonstrate that IgG assists IgM in coating of most bacterial families in the absence of IgA and identify single-IgG2 bacterial coating as an inflammatory marker.

Supplementary Table 2 from Analysis of Gene Expression Profiles of Microdissected Cell Populations Indicates that Testicular Carcinoma In situ Is an Arrested Gonocyte

Supplementary Table 4 from Analysis of Gene Expression Profiles of Microdissected Cell Populations Indicates that Testicular Carcinoma In situ Is an Arrested Gonocyte

Supplementary Figure 3 from Analysis of Gene Expression Profiles of Microdissected Cell Populations Indicates that Testicular Carcinoma In situ Is an Arrested Gonocyte

Testicular cancer is predominantly curable, but the long-term side effects of chemotherapy have a severe impact on life quality. In this research study, we focus on hearing loss as a part of overall chemotherapy-induced ototoxicity. This is a unique approach where we combine clinical data from the acclaimed nationwide Danish Testicular Cancer (DaTeCa)-Late database. Clinical and genetic data on 433 patients were collected from hospital files in October 2014. Hearing loss was classified according to the FACT/GOG-Ntx-11 version 4 self-reported Ntx6. Machine learning models combining a genome-wide association study within a nested cross-validated logistic regression were applied to identify patients at high risk of hearing loss. The model comprising clinical and genetic data identified 67% of the patients with hearing loss; however, this was with a false discovery rate of 49%. For the non-affected patients, the model identified 66% of the patients with a false omission rate of 19%. An area under the receiver operating characteristic (ROC-AUC) curve of 0.73 (95% CI, 0.71-0.74) was obtained, and the model suggests genes SOD2 and MGST3 as important in improving prediction over the clinical-only model with a ROC-AUC of 0.66 (95% CI, 0.65-0.66). Such prediction models may be used to allow earlier detection and prevention of hearing loss. We suggest a possible biological mechanism for cisplatin-induced hearing loss development. On confirmation in larger studies, such models can help balance treatment in clinical practice.

[This corrects the article DOI: 10.3389/fmicb.2017.01934.].

AbstractBackground: Testicular germ cell tumors (TGCT), histologically classified as seminomas and nonseminomas, are believed to arise from primordial gonocytes, with the maturation process blocked when they are subjected to DNA methylation reprogramming. SNPs in DNA methylation machinery and folate-dependent one-carbon metabolism genes have been postulated to influence the proper establishment of DNA methylation. Methods: In this pathway-focused investigation, we evaluated the association between 273 selected tag SNPs from 28 DNA methylation–related genes and TGCT risk. We carried out association analysis at individual SNP and gene-based level using summary statistics from the Genome Wide Association Study meta-analysis recently conducted by the international Testicular Cancer Consortium on 10,156 TGCT cases and 179,683 controls. Results: In individual SNP analyses, seven SNPs, four mapping within MTHFR, were associated with TGCT risk after correction for multiple testing (q ≤ 0.05). Queries of public databases showed that three of these SNPs were associated with MTHFR changes in enzymatic activity (rs1801133) or expression level in testis tissue (rs12121543, rs1476413). Gene-based analyses revealed MTHFR (q = 8.4 × 10–4), methyl-CpG–binding protein 2 (MECP2; q = 2 × 10–3), and ZBTB4 (q = 0.03) as the top TGCT-associated genes. Stratifying by tumor histology, four MTHFR SNPs were associated with seminoma. In gene-based analysis MTHFR was associated with risk of seminoma (q = 2.8 × 10–4), but not with nonseminomatous tumors (q = 0.22). Conclusions: Genetic variants within MTHFR, potentially having an impact on the DNA methylation pattern, are associated with TGCT risk. Impact: This finding suggests that TGCT pathogenesis could be associated with the folate cycle status, and this relation could be partly due to hereditary factors.

Purpose: An altered ocular surface microbiota may contribute to the pathophysiology of dry eye disease. The aim of the study was to explore potential differences in microbiota diversity and composition in aqueous tear deficient dry eye (with and without ocular graft-versus-host disease) compared with controls. Methods: Swab samples from the inferior fornix of the conjunctiva were obtained from patients with aqueous tear-deficient dry eye with and without ocular graft-versus-host disease (n = 18, n = 21, respectively) and controls (n = 28). Isolated bacterial DNA from swabs were analyzed with 16S rRNA gene amplicon sequencing. Results: Decreased microbiota diversity was observed in patients with aqueous tear-deficient dry eye (p <= 0.003) who also showed a difference in microbiota composition compared with controls (p = 0.001). Although several genera were less abundant in aqueous tear-deficient dry eye, a minimal core ocular surface microbiota comprising five genera was shared by >75% of the study participants: Enhydrobacter, Brevibacterium, Staphylococcus, Streptococcus and Cutibacterium. Pseudomonas was identified as a bacterial biomarker for controls and Bacilli for patients with aqueous tear-deficient dry eye. Conclusions: Ocular surface microbiota in patients with aqueous tear-deficient dry eye was characterized by an aberrant microbiota composition in comparison to controls, with decreased diversity and reduced relative abundances of several genera. Additionally, a few genera were present in most of the study population, indicating that a minimal core ocular surface microbiota may exist.

Background: Studies evaluating the association between peripheral blood leukocyte telomere length (LTL) and testicular germ cell tumor (TGCT) risk have produced conflicting results. Methods: Using available genotype data from the Testicular Cancer Consortium (TECAC), polygenic risk score and Mendelian randomization analyses of genetic variants previously associated with LTL were used to assess potential etiologic associations between telomere length and TGCT risk. Results: Genetically inferred telomere length was not associated with TGCT risk among 2,049 cases and 6,921 controls with individual-level genotype data (OR, 1.02; 95% confidence interval, 0.97-1.07). Mendelian randomization analyses using summary statistic data further indicated no evidence for an association between telomere length and TGCT risk among all available TECAC participants (3,558 cases and 13,971 controls). Conclusions: Our analyses in the largest molecular genetic testicular cancer study to date provide no evidence for an association between genetically inferred peripheral blood LTL and TGCT risk. Impact: The lack of evidence for an overall association indicates that peripheral blood LTL is likely not a strong biomarker for TGCT risk.

Our objectives were to investigate whether the conjunctival microbiota is altered by contact lens wear and/or bacterial keratitis and to explore the hypothesis that commensals of conjunctival microbiota contribute to bacterial keratitis. Swab samples from both eyes were collected separately from the inferior fornix of the conjunctiva of non-contact-lens users (nparticipants = 28) and contact lens users (nparticipants = 26) and from patients with contact-lens-associated bacterial keratitis (nparticipants = 9). DNA from conjunctival swab samples was analyzed with 16S rRNA gene amplicon sequencing. Pathogens from the corneal infiltrates were identified by cultivation. In total, we identified 19 phyla and 283 genera; the four most abundant genera were Pseudomonas, Enhydrobacter, Staphylococcus, and Cutibacterium. Several pathogens related to bacterial keratitis were identified in the conjunctival microbiota of the whole study population, and the same bacteria were identified by both methods in the conjunctiva and cornea for four patients with contact-lens-associated bacterial keratitis. The overall conjunctival microbiota profile was not altered by contact lens wear or bacterial keratitis; thus, it does not appear to contribute to the development of bacterial keratitis in contact lens users. However, in some individuals, conjunctival microbiota may harbor opportunistic pathogens causing contact-lens-associated bacterial keratitis.

Conazole fungicides such as epoxiconazole are mostly used on cereals of crops to inhibit fungal growth through direct inhibition of sterol 14α-demethylase (CYP51A1).However, this enzyme is highly conserved and in humans it is part of the steroid hormone biosynthesis pathway.Endocrine disrupting effects of epoxiconazole have been shown in rodents and have been substantiated by in vitro data, however, the underlying molecular mechanisms are not clear.We took advantage of a human stem cell based in vitro model for developmental toxicity to study the molecular effects of epoxiconazole.This model is based on 3D cultures of embryoid bodies and differentiation into cardiomyocytes, which mimics the early stages of embryonic development.We have previously shown that epoxiconazole impairs differentiation of these embryoid bodies and therefore has the potential to affect human embryonic development.We employed global transcriptome analysis using RNA sequencing and found that the steroid biosynthesis pathway including CYP51A1, the human sterol 14α-demethylase, was highly deregulated by epoxiconazole in our model.We confirmed that most genes of the steroid biosynthesis pathway were upregulated, including CYP51A1, suggesting a compensatory mechanism at the gene expression level.Our data suggest that epoxiconazole acts mainly by decreasing cholesterol biosynthesis in the cells.We conclude that epoxiconazole bears the potential to harm human embryonic development through inhibition of the steroid biosynthesis pathway.As this may be a common feature of compounds that target sterol 14α-demethylase, we add evidence to the assumption that conazole fungicides may be human developmental toxicants.

Circulating miRNAs secreted by testicular germ cell tumors (TGCT) show great potential as novel non-invasive biomarkers for diagnosis of TGCT. Seminal plasma (SP) represents a biofluid closer to the primary site. Here, we investigate whether small RNAs in SP can be used to diagnose men with TGCTs or the precursor lesions, germ cell neoplasia in situ (GCNIS). Small RNAs isolated from SP from men with TGCTs (n = 18), GCNIS-only (n = 5), and controls (n = 25) were sequenced. SP from men with TGCT/GCNIS (n = 37) and controls (n = 22) were used for validation by RT-qPCR. In general, piRNAs were found at lower levels in SP from men with TGCTs. Ten small RNAs were found at significantly (q-value < 0.05) different levels in SP from men with TGCT/GCNIS than controls. Random forests classification identified sets of small RNAs that could detect either TGCT/GCNIS or GCNIS-only with an area under the curve of 0.98 and 1 in ROC analyses, respectively. RT-qPCR validated hsa-miR-6782-5p to be present at 2.3-fold lower levels (p = 0.02) in the SP from men with TGCTs compared with controls. Small RNAs in SP show potential as novel biomarkers for diagnosing men with TGCT/GCNIS but validation in larger cohorts is needed.

The analgesic paracetamol (N-acetyl-4-aminophenol, APAP) is commonly used to relieve pain, fever and malaise. While sales have increased worldwide, a growing body of experimental and epidemiological evidence has suggested APAP as a possible risk factor for various health disorders. To perform internal exposure-based risk assessment, the use of accurate and optimized biomonitoring methods is criticical. However, retrospectively assessing pharmaceutical use of APAP in humans is challenging because of its short half-life. The objective of this study was to address the key biomonitoring issues with APAP using current standard analytical methods based on urinary analyses of free APAP and its phase II conjugates. Using non-targeted analyses based on high-resolution mass spectrometry, we identified in a controlled longitudinal exposure study with male volunteers, unrecognized APAP metabolites with delayed formation and excretion rates. We postulate that these metabolites are formed via the thiomethyl shunt after the enterohepatic circulation as already observed in rodents. Importantly, the conjugated thiomethyl metabolites were (i) of comparable diagnostic sensitivity as the free APAP and its phase II conjugates detected by current methods; (ii) had delayed peak levels in blood and urine compared to other APAP metabolites and therefore potentially extend the window of exposure assessment; and (iii) provide relevant information regarding metabolic pathways of interest from a toxicological point of view. Including these metabolites in future APAP biomonitoring methods provide an option to decrease potential underestimation of APAP use and challenges the notion that the standard methods in biomonitoring based exclusively on the parent compound and its phase II metabolites are adequate for human biomonitoring of non-persistant chemical such as APAP.

Diet is an important component in weight management strategies, but heterogeneous responses to the same diet make it difficult to foresee individual weight-loss outcomes. Omics-based technologies now allow for analysis of multiple factors for weight loss prediction at the individual level. Here, we classify weight loss responders (N = 106) and non-responders (N = 97) of overweight non-diabetic middle-aged Danes to two earlier reported dietary trials over 8 weeks. Random forest models integrated gut microbiome, host genetics, urine metabolome, measures of physiology and anthropometrics measured prior to any dietary intervention to identify individual predisposing features of weight loss in combination with diet. The most predictive models for weight loss included features of diet, gut bacterial species and urine metabolites (ROC-AUC: 0.84–0.88) compared to a diet-only model (ROC-AUC: 0.62). A model ensemble integrating multi-omics identified 64% of the non-responders with 80% confidence. Such models will be useful to assist in selecting appropriate weight management strategies, as individual predisposition to diet response varies.

Genome-wide association studies (GWAS) have led to the identification of hundreds of susceptibility loci across cancers, but the impact of further studies remains uncertain. Here we analyse summary-level data from GWAS of European ancestry across fourteen cancer sites to estimate the number of common susceptibility variants (polygenicity) and underlying effect-size distribution. All cancers show a high degree of polygenicity, involving at a minimum of thousands of loci. We project that sample sizes required to explain 80% of GWAS heritability vary from 60,000 cases for testicular to over 1,000,000 cases for lung cancer. The maximum relative risk achievable for subjects at the 99th risk percentile of underlying polygenic risk scores (PRS), compared to average risk, ranges from 12 for testicular to 2.5 for ovarian cancer. We show that PRS have potential for risk stratification for cancers of breast, colon and prostate, but less so for others because of modest heritability and lower incidence.

Background Congenital heart disease (CHD) occurs in almost 1% of newborn children and is considered a multifactorial disorder. CHD may segregate in families due to significant contribution of genetic factors in the disease etiology. The aim of the study was to identify pathophysiological mechanisms in families segregating CHD. Methods We used whole exome sequencing to identify rare genetic variants in ninety consenting participants from 32 Danish families with recurrent CHD. We applied a systems biology approach to identify developmental mechanisms influenced by accumulation of rare variants. We used an independent cohort of 714 CHD cases and 4922 controls for replication and performed functional investigations using zebrafish as in vivo model. Results We identified 1785 genes, in which rare alleles were shared between affected individuals within a family. These genes were enriched for known cardiac developmental genes, and 218 of these genes were mutated in more than one family. Our analysis revealed a functional cluster, enriched for proteins with a known participation in calcium signaling. Replication in an independent cohort confirmed increased mutation burden of calcium-signaling genes in CHD patients. Functional investigation of zebrafish orthologues ofITPR1,PLCB2, andADCY2verified a role in cardiac development and suggests a combinatorial effect of inactivation of these genes. Conclusions The study identifies abnormal calcium signaling as a novel pathophysiological mechanism in human CHD and confirms the complex genetic architecture underlying CHD.

Background: Cisplatin-based chemotherapy may induce nephrotoxicity. This study presents a random forest predictive model that identifies testicular cancer patients at risk of nephrotoxicity before treatment. Methods: Clinical data and DNA from saliva samples were collected for 433 patients. These were genotyped on Illumina Human Omni Express Exome-8 v1.2 (964 193 markers). Clinical and genomics-based random forest models generated a risk score for each individual to develop nephrotoxicity defined as a 20% drop in isotopic glomerular filtration rate during chemotherapy. The area under the receiver operating characteristic curve was the primary measure to evaluate models. Sensitivity, specificity, and positive and negative predictive values were used to discuss model clinical utility. Results: Of 433 patients assessed in this study, 26.8% developed nephrotoxicity after bleomycin-etoposide-cisplatin treatment. Genomic markers found to be associated with nephrotoxicity were located at NAT1, NAT2, and the intergenic region of CNTN6 and CNTN4. These, in addition to previously associated markers located at ERCC1, ERCC2, and SLC22A2, were found to improve predictions in a clinical feature-trained random forest model. Using only clinical data for training the model, an area under the receiver operating characteristic curve of 0.635 (95% confidence interval [CI] = 0.629 to 0.640) was obtained. Retraining the classifier by adding genomics markers increased performance to 0.731 (95% CI = 0.726 to 0.736) and 0.692 (95% CI = 0.688 to 0.696) on the holdout set. Conclusions: A clinical and genomics-based machine learning algorithm improved the ability to identify patients at risk of nephrotoxicity compared with using clinical variables alone. Novel genetics associations with cisplatin-induced nephrotoxicity were found for NAT1, NAT2, CNTN6, and CNTN4 that require replication in larger studies before application to clinical practice.