Genome-wide association study of early-onset and late-onset postpartum depression: the IGEDEPP prospective study

Postpartum depression (PPD) appears at two peak periods: early-onset prior to 2 months after delivery and late-onset (2 months after delivery and beyond). The aim of our study is to evaluate the different genetic factors associated with early- and late-onset PPD. With the French multicenter interaction of gene and environment of depression during postpartum (IGEDEPP) cohort, we conducted a genome-wide association study (GWAS) on 234 women with early-onset PPD and 223 women with late-onset PPD, as well as 1,204 controls with no history of lifetime depression. We performed post-GWAS analyses: functional mapping and annotation of GWAS results using MAGMA thanks to Functional Mapping and Annotation of Genome-Wide Association Studies (FUMA), expression quantitative trait loci (QTL) analyses, mapping using data from the PsychENCODE and GTEx, and polygenic risk score (PRS) analysis based on published GWAS. We found two new significant candidate loci for early-onset PPD, rs6436132 in PTPRN gene on chromosome 2 and rs184644645 in RAD18 on chromosome 14, respectively, and one region of interest with five significant associated SNPs in chromosome 20 for late-onset PPD. Variant rs6436132 is the most significant associated with early-onset PPD, and it is a QTL that significantly modifies the expression and splicing of the PTPRN gene in different brain tissues. We also found an enrichment of uterus tissue in the early expression of PPD genes. PRS analysis showed a genetic overlap between both early and late-onset PPD and major depressive disorder, but only early-onset PPD overlaps with bipolar disorder. Our study presents two GWAS separately, highlighting two candidate loci for early-onset PPD and one different region of interest for late-onset PPD. These results have important consequences in our understanding of these disorders, especially since our data reinforce the hormonal pathophysiological hypotheses for early-onset PPD.

An elevated level of interleukin-17A in a Senegalese malaria cohort is associated with rs8193038 IL-17A genetic variant

Malaria infection is a multifactorial disease partly modulated by host immuno-genetic factors. Recent evidence has demonstrated the importance of Interleukin-17 family proinflammatory cytokines and their genetic variants in host immunity. However, limited knowledge exists about their role in parasitic infections such as malaria. We aimed to investigate IL-17A serum levels in patients with severe and uncomplicated malaria and gene polymorphism's influence on the IL-17A serum levels. In this research, 125 severe (SM) and uncomplicated (UM) malaria patients and 48 free malaria controls were enrolled. IL-17A serum levels were measured with ELISA. PCR and DNA sequencing were used to assess host genetic polymorphisms in IL-17A. We performed a multivariate regression to estimate the impact of human IL-17A variants on IL-17A serum levels and malaria outcomes. Elevated serum IL-17A levels accompanied by increased parasitemia were found in SM patients compared to UM and controls (P < 0.0001). Also, the IL-17A levels were lower in SM patients who were deceased than in those who survived. In addition, the minor allele frequencies (MAF) of two IL-17A polymorphisms (rs3819024 and rs3748067) were more prevalent in SM patients than UM patients, indicating an essential role in SM. Interestingly, the heterozygous rs8193038 AG genotype was significantly associated with higher levels of IL-17A than the homozygous wild type (AA). According to our results, it can be concluded that the IL-17A gene rs8193038 polymorphism significantly affects IL-17A gene expression. Our results fill a gap in the implication of IL-17A gene polymorphisms on the cytokine level in a malaria cohort. IL-17A gene polymorphisms also may influence cytokine production in response to Plasmodium infections and may contribute to the hyperinflammatory responses during severe malaria outcomes.

Highly efficient and sensitive membrane-based concentration process allows quantification, surveillance, and sequencing of viruses in large volumes of wastewater

Wastewater-based epidemiology is experiencing exponential development. Despite undeniable advantages compared to patient-centered approaches (cost, anonymity, survey of large populations without bias, detection of asymptomatic infected peoples…), major technical limitations persist. Among them is the low sensitivity of the current methods used for quantifying and sequencing viral genomes from wastewater. In situations of low viral circulation, during initial stages of viral emergences, or in areas experiencing heavy rains, the extremely low concentrations of viruses in wastewater may fall below the limit of detection of the current methods. The availability during crisis and the cost of the commercial kits, as well as the requirement of expensive materials such as high-speed centrifuge, can also present major blocks to the development of wastewater-based epidemiological survey, specifically in low-income countries. Thereby, highly sensitive, low cost and standardized methods are still needed, to increase the predictability of the viral emergences, to survey low-circulating viruses and to make the results from different labs comparable. Here, we outline and characterize new protocols for concentrating and quantifying SARS-CoV-2 from large volumes (500 mL-1 L) of untreated wastewater. In addition, we report that the methods are applicable for monitoring and sequencing. Our nucleic acid extraction technique (the routine C: 5 mL method) does not require sophisticated equipment such as automatons and is not reliant on commercial kits, making it readily available to a broader range of laboratories for routine epidemiological survey. Furthermore, we demonstrate the efficiency, the repeatability, and the high sensitivity of a new membrane-based concentration method (MBC: 500 mL method) for enveloped (SARS-CoV-2) and non-enveloped (F-specific RNA phages of genogroup II / FRNAPH GGII) viruses. We show that the MBC method allows the quantification and the monitoring of viruses in wastewater with a significantly improved sensitivity compared to the routine C method. In contexts of low viral circulation, we report quantifications of SARS-CoV-2 in wastewater at concentrations as low as 40 genome copies per liter. In highly diluted samples collected in wastewater treatment plants of French Guiana, we confirmed the accuracy of the MBC method compared to the estimations done with the routine C method. Finally, we demonstrate that both the routine C method processing 5 mL and the MBC method processing 500 mL of untreated wastewater are both compatible with SARS-CoV-2 sequencing. We show that the quality of the sequence is correlated with the concentration of the extracted viral genome. Of note, the quality of the sequences obtained with some MBC processed wastewater was improved by dilutions or enzyme substitutions suggesting the presence of specific enzyme inhibitors in some wastewater. To the best of our knowledge, our MBC method is one of the first efficient, sensitive, and repeatable method characterized for SARS-CoV-2 quantification and sequencing from large volumes of wastewater.

Multi-ancestry polygenic risk scores for venous thromboembolism

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium GWAS meta-analyses of European- (71,771 cases and 1,059,740 controls) and African-ancestry samples (7,482 cases and 129,975 controls). We used LDpred2 and PRSCSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6,261 cases and 88,238 controls) and African-ancestry sample (1,385 cases and 12,569 controls). Multi-ancestry PRSs with weights tuned in European- and African-ancestry samples, respectively, outperformed ancestry-specific PRSs in European- (PRSCSXEUR: AUC=0.61 (0.60, 0.61), PRSCSX_combinedEUR: AUC=0.61 (0.60, 0.62)) and African-ancestry test samples (PRSCSXAFR: AUC=0.58 (0.57, 0.6), PRSCSX_combined AFR: AUC=0.59 (0.57, 0.60)). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS may be used to identify individuals at highest risk for VTE and provide guidance for the most effective treatment strategy across diverse populations.

Elevated plasma complement factor H related 5 protein is associated with venous thromboembolism

Venous thromboembolism (VTE) is a common, multi-causal disease with potentially serious short- and long-term complications. In clinical practice, there is a need for improved plasma biomarker-based tools for VTE diagnosis and risk prediction. Here we show, using proteomics profiling to screen plasma from patients with suspected acute VTE, and several case-control studies for VTE, how Complement Factor H Related 5 protein (CFHR5), a regulator of the alternative pathway of complement activation, is a VTE-associated plasma biomarker. In plasma, higher CFHR5 levels are associated with increased thrombin generation potential and recombinant CFHR5 enhanced platelet activation in vitro. GWAS analysis of ~52,000 participants identifies six loci associated with CFHR5 plasma levels, but Mendelian randomization do not demonstrate causality between CFHR5 and VTE. Our results indicate an important role for the regulation of the alternative pathway of complement activation in VTE and that CFHR5 represents a potential diagnostic and/or risk predictive plasma biomarker.

Dissecting the genetic heterogeneity of gastric cancer

BACKGROUND

Gastric cancer (GC) is clinically heterogenous according to location (cardia/non-cardia) and histopathology (diffuse/intestinal). We aimed to characterize the genetic risk architecture of GC according to its subtypes. Another aim was to examine whether cardia GC and oesophageal adenocarcinoma (OAC) and its precursor lesion Barrett's oesophagus (BO), which are all located at the gastro-oesophageal junction (GOJ), share polygenic risk architecture.

METHODS

We did a meta-analysis of ten European genome-wide association studies (GWAS) of GC and its subtypes. All patients had a histopathologically confirmed diagnosis of gastric adenocarcinoma. For the identification of risk genes among GWAS loci we did a transcriptome-wide association study (TWAS) and expression quantitative trait locus (eQTL) study from gastric corpus and antrum mucosa. To test whether cardia GC and OAC/BO share genetic aetiology we also used a European GWAS sample with OAC/BO.

FINDINGS

Our GWAS consisting of 5816 patients and 10,999 controls highlights the genetic heterogeneity of GC according to its subtypes. We newly identified two and replicated five GC risk loci, all of them with subtype-specific association. The gastric transcriptome data consisting of 361 corpus and 342 antrum mucosa samples revealed that an upregulated expression of MUC1, ANKRD50, PTGER4, and PSCA are plausible GC-pathomechanisms at four GWAS loci. At another risk locus, we found that the blood-group 0 exerts protective effects for non-cardia and diffuse GC, while blood-group A increases risk for both GC subtypes. Furthermore, our GWAS on cardia GC and OAC/BO (10,279 patients, 16,527 controls) showed that both cancer entities share genetic aetiology at the polygenic level and identified two new risk loci on the single-marker level.

INTERPRETATION

Our findings show that the pathophysiology of GC is genetically heterogenous according to location and histopathology. Moreover, our findings point to common molecular mechanisms underlying cardia GC and OAC/BO.

FUNDING

German Research Foundation (DFG).

Cost of exome analysis in patients with intellectual disability: a micro-costing study in a French setting

BACKGROUND

With the development of next generation sequencing technologies in France, exome sequencing (ES) has recently emerged as an opportunity to improve the diagnosis rate of patients presenting an intellectual disability (ID). To help French policy makers determine an adequate tariff for ES, we aimed to assess the unit cost per ES diagnostic test for ID from the preparation of the pre-analytical step until the report writing step and to identify its main cost drivers.

METHODS

A micro-costing bottom-up approach was conducted for the year 2018 in a French setting as part of the DISSEQ study, a cost-effectiveness study funded by the Ministry of Health and performed in collaboration with the GAD (Génétique des Anomalies du Développement), a genetic team from the Dijon University Hospital, and a public sequencing platform, the Centre National de Recherche en Génomique Humaine (CNRGH). The analysis was conducted from the point of view of these two ES stakeholders. All of the resources (labor, equipment, disposables and reagents, reusable material) required to analyze blood samples were identified, collected and valued. Several sensitivity analyses were performed.

RESULTS

The unit nominal cost per ES diagnostic test for ID was estimated to be €2,019.39. Labor represented 50.7% of the total cost. The analytical step (from the preparation of libraries to the analysis of sequences) represented 88% of the total cost. Sensitivity analyses suggested that a simultaneous price decrease of 20% for the capture kit and 50% for the sequencing support kit led to an estimation of €1,769 per ES diagnostic test for ID.

CONCLUSION

This is the first estimation of ES cost to be done in the French setting of ID diagnosis. The estimation is especially influenced by the price of equipment kits, but more generally by the organization of the centers involved in the different steps of the analysis and the time period in which the study was conducted. This information can now be used to define an adequate tariff and assess the efficiency of ES.

TRIAL REGISTRATION

ClinicalTrials.gov identifier NCT03287206 on September 19, 2017.

Multiomic analysis of malignant pleural mesothelioma identifies molecular axes and specialized tumor profiles driving intertumor heterogeneity

Malignant pleural mesothelioma (MPM) is an aggressive cancer with rising incidence and challenging clinical management. Through a large series of whole-genome sequencing data, integrated with transcriptomic and epigenomic data using multiomics factor analysis, we demonstrate that the current World Health Organization classification only accounts for up to 10% of interpatient molecular differences. Instead, the MESOMICS project paves the way for a morphomolecular classification of MPM based on four dimensions: ploidy, tumor cell morphology, adaptive immune response and CpG island methylator profile. We show that these four dimensions are complementary, capture major interpatient molecular differences and are delimited by extreme phenotypes that-in the case of the interdependent tumor cell morphology and adapted immune response-reflect tumor specialization. These findings unearth the interplay between MPM functional biology and its genomic history, and provide insights into the variations observed in the clinical behavior of patients with MPM.

Heritable defects in telomere and mitotic function selectively predispose to sarcomas

Cancer genetics has to date focused on epithelial malignancies, identifying multiple histotype-specific pathways underlying cancer susceptibility. Sarcomas are rare malignancies predominantly derived from embryonic mesoderm. To identify pathways specific to mesenchymal cancers, we performed whole-genome germline sequencing on 1644 sporadic cases and 3205 matched healthy elderly controls. Using an extreme phenotype design, a combined rare-variant burden and ontologic analysis identified two sarcoma-specific pathways involved in mitotic and telomere functions. Variants in centrosome genes are linked to malignant peripheral nerve sheath and gastrointestinal stromal tumors, whereas heritable defects in the shelterin complex link susceptibility to sarcoma, melanoma, and thyroid cancers. These studies indicate a specific role for heritable defects in mitotic and telomere biology in risk of sarcomas.

G6PD and HBB polymorphisms in the Senegalese population: prevalence, correlation with clinical malaria

Background

Host genetic factors contribute to the variability of malaria phenotypes and can allow a better understanding of mechanisms involved in susceptibility and/or resistance to Plasmodium falciparum infection outcomes. Several genetic polymorphisms were reported to be prevalent among populations living in tropical malaria-endemic regions and induce protection against malaria. The present study aims to investigate the prevalence of HBB (chr11) and G6PD (chrX) deficiencies polymorphisms among Senegalese populations and their associations with the risk for severe Plasmodium falciparum malaria occurrence.

Methods

We performed a retrospective study with 437 samples, 323 patients recruited in hospitals located in three different endemic areas where malaria episodes were confirmed and 114 free malaria controls. The patients enrolled were classified into two groups: severe malaria (SM) (153 patients) and uncomplicated malaria (UM) (170 patients). PCR and DNA sequencing assessed host genetic polymorphisms in HBB and G6PD. Using a multivariate regression and additive model, estimates of the impact of human HBB and G6PD polymorphisms on malaria incidence were performed.

Results

Six frequent SNPs with minor allele frequencies (MAF) > 3% were detected in the HBB gene (rs7946748, rs7480526, rs10768683, rs35209591, HbS (rs334) and rs713040) and two in the G6PD gene (rs762515 and rs1050828 (G6PD-202 G > A). Analysis of selected HbS polymorphism showed significant association with protective effect against severe malaria with a significant p-value = 0.033 (OR 0.38, 95% CI [0.16-0.91]) for SM vs. UM comparison. Surprisingly, our study did not identify the protective effect of variant HbC polymorphism against severe malaria. Finally, we found some of the polymorphisms, like HbS (rs334), are associated with age and biological parameters like eosinophils, basophils, lymphocytes etc.

Conclusion

Our data report HBB and G6PD polymorphisms in the Senegalese population and their correlation with severe/mild malaria and outcome. The G6PD and HBB deficiencies are widespread in West Africa endemic malaria regions such as The Gambia, Mali, and Burkina Faso. The study shows the critical role of genetic factors in malaria outcomes. Indeed, genetic markers could be good tools for malaria endemicity prognosis.

Genome-wide meta-analysis for Alzheimer's disease cerebrospinal fluid biomarkers

Amyloid-beta 42 (Aβ42) and phosphorylated tau (pTau) levels in cerebrospinal fluid (CSF) reflect core features of the pathogenesis of Alzheimer's disease (AD) more directly than clinical diagnosis. Initiated by the European Alzheimer & Dementia Biobank (EADB), the largest collaborative effort on genetics underlying CSF biomarkers was established, including 31 cohorts with a total of 13,116 individuals (discovery n = 8074; replication n = 5042 individuals). Besides the APOE locus, novel associations with two other well-established AD risk loci were observed; CR1 was shown a locus for Aβ42 and BIN1 for pTau. GMNC and C16orf95 were further identified as loci for pTau, of which the latter is novel. Clustering methods exploring the influence of all known AD risk loci on the CSF protein levels, revealed 4 biological categories suggesting multiple Aβ42 and pTau related biological pathways involved in the etiology of AD. In functional follow-up analyses, GMNC and C16orf95 both associated with lateral ventricular volume, implying an overlap in genetic etiology for tau levels and brain ventricular volume.

Integrated genomic analysis identifies driver genes and cisplatin-resistant progenitor phenotype in pediatric liver cancer

Pediatric liver cancers (PLCs) comprise diverse diseases affecting infants, children and adolescents. Despite overall good prognosis, PLCs display heterogeneous response to chemotherapy. Integrated genomic analysis of 126 pediatric liver tumors showed a continuum of driver mechanisms associated with patient age, including new targetable oncogenes. In 10% of hepatoblastoma patients, all before 3 years old, we identified a mosaic premalignant clonal expansion of cells altered at the 11p15.5 locus. Analysis of spatial and longitudinal heterogeneity revealed an important plasticity between 'Hepatocytic', 'Liver Progenitor' and 'Mesenchymal' molecular subgroups of hepatoblastoma. We showed that during chemotherapy, 'Liver Progenitor' cells accumulated massive loads of cisplatin-induced mutations with a specific mutational signature, leading to the development of heavily mutated relapses and metastases. Drug screening in PLC cell lines identified promising targets for cisplatin-resistant progenitor cells, validated in mouse xenograft experiments. These data provide new insights into cisplatin resistance mechanisms in PLC and suggest alternative therapies.

INPP5K and SIL1 associated pathologies with overlapping clinical phenotypes converge through dysregulation of PHGDH

Marinesco-Sjögren syndrome (MSS) is a rare human disorder caused by biallelic mutations in SIL1 characterized by cataracts in infancy, myopathy and ataxia, symptoms that are also associated with a novel disorder caused by mutations in INPP5K. While these phenotypic similarities may suggest commonalties at a molecular level, an overlapping pathomechanism has not been established yet. In this study, we present six new INPP5K patients and expand the current mutational and phenotypical spectrum of the disease showing the clinical overlap between MSS and the INPP5K-phenotype. We applied unbiased proteomic profiling on cells derived from MSS- and INPP5K-patients and identified alterations in D-3-phosphoglycerate dehydrogenase as a common molecular feature. D-3-phosphoglycerate dehydrogenase modulates the production of L-serine and mutations in this enzyme were previously associated with a neurological phenotype, which clinically overlaps with MSS and INPP5K-disease. As, L-serine administration represents a promising therapeutic strategy for D-3-phosphoglycerate dehydrogenase patients, we tested the effect of L-serine in generated sil1, phgdh and inpp5k a + b zebrafish models which showed an improvement in their neuronal phenotype. Thus our study defines a core phenotypical feature underpinning a key common molecular mechanism in three rare diseases and reveals a common and novel therapeutic target for these patients.

Analysis and annotation of DNA methylation in two nonhuman primate species using the Infinium Human Methylation 450K and EPIC BeadChips

Aim: Nonhuman primates are essential for research on many human diseases. The Infinium Human Methylation450/EPIC BeadChips are popular tools for the study of the methylation state across the human genome at affordable cost. Methods: We performed a precise evaluation and re-annotation of the BeadChip probes for the analysis of genome-wide DNA methylation patterns in rhesus macaques and African green monkeys through in silico analyses combined with functional validation by pyrosequencing. Results: Up to 165,847 of the 450K and 261,545 probes of the EPIC BeadChip can be reliably used. The annotation files are provided in a format compatible with a variety of standard bioinformatic pipelines. Conclusion: Our study will facilitate high-throughput DNA methylation analyses in Macaca mulatta and Chlorocebus sabaeus.

A novel pathogenic variant in DYNC1H1 causes various upper and lower motor neuron anomalies

OBJECTIVE

To perform genotype-phenotype, clinical and molecular analysis in a large 3-generation family with autosomal dominant congenital spinal muscular atrophy.

METHODS

Using a combined genetic approach including whole genome scanning, next generation sequencing-based multigene panel, whole genome sequencing, and targeted variant Sanger sequencing, we studied the proband and multiple affected individuals of this family who presented bilateral proximal lower limb muscle weakness and atrophy.

RESULTS

We identified a novel heterozygous variant, c.1826T > C; p.Ile609Thr, in the DYNC1H1 gene localized within the common haplotype in the 14q32.3 chromosomal region which cosegregated with disease in this large family. Within the family, affected individuals were found to have a wide array of clinical variability. Although some individuals presented the typical lower motor neuron phenotype with areflexia and denervation, others presented with muscle weakness and atrophy, hyperreflexia, and absence of denervation suggesting a predominant upper motor neuron disease. In addition, some affected individuals presented with an intermediate phenotype characterized by hyperreflexia and denervation, expressing a combination of lower and upper motor neuron defects.

CONCLUSION

Our study demonstrates the wide clinical variability associated with a single disease causing variant in DYNC1H1 gene and this variant demonstrated a high penetrance within this large family.

PLCG2 protective variant p.P522R modulates tau pathology and disease progression in patients with mild cognitive impairment

A rare coding variant (rs72824905, p.P522R) conferring protection against Alzheimer’s disease (AD) was identified in the gene encoding the enzyme phospholipase-C-γ2 (PLCG2) that is highly expressed in microglia. To explore the protective nature of this variant, we employed latent process linear mixed models to examine the association of p.P522R with longitudinal cognitive decline in 3595 MCI patients, and in 10,097 individuals from population-based studies. Furthermore, association with CSF levels of pTau₁₈₁, total tau, and Aβ_1-42 was assessed in 1261 MCI patients. We found that MCI patients who carried the p.P522R variant showed a slower rate of cognitive decline compared to non-carriers and that this effect was mediated by lower pTau₁₈₁ levels in CSF. The effect size of the association of p.P522R with the cognitive decline and pTau₁₈₁ was similar to that of APOE-ε4, the strongest genetic risk factor for AD. Interestingly, the protective effect of p.P522R was more pronounced in MCI patients with low Aβ_1-42 levels suggesting a role of PLCG2 in the response to amyloid pathology. In line with this hypothesis, we observed no protective effect of the PLCG2 variant on the cognitive decline in population-based studies probably due to the lower prevalence of amyloid positivity in these samples compared to MCI patients. Concerning the potential biological underpinnings, we identified a network of co-expressed proteins connecting PLCG2 to APOE and TREM2 using unsupervised co-regulatory network analysis. The network was highly enriched for the complement cascade and genes differentially expressed in disease-associated microglia. Our data show that p.P522R in PLCG2 reduces AD disease progression by mitigating tau pathology in the presence of amyloid pathology and, as a consequence, maintains cognitive function. Targeting the enzyme PLCG2 might provide a new therapeutic approach for treating AD.

MACARON: a python framework to identify and re-annotate multi-base affected codons in whole genome/exome sequence data

Summary

Predicted deleteriousness of coding variants is a frequently used criterion to filter out variants detected in next-generation sequencing projects and to select candidates impacting on the risk of human diseases. Most available dedicated tools implement a base-to-base annotation approach that could be biased in presence of several variants in the same genetic codon. We here proposed the MACARON program that, from a standard VCF file, identifies, re-annotates and predicts the amino acid change resulting from multiple single nucleotide variants (SNVs) within the same genetic codon. Applied to the whole exome dataset of 573 individuals, MACARON identifies 114 situations where multiple SNVs within a genetic codon induce an amino acid change that is different from those predicted by standard single SNV annotation tool. Such events are not uncommon and deserve to be studied in sequencing projects with inconclusive findings.

Availability and implementation

MACARON is written in python with codes available on the GENMED website (www.genmed.fr).

Supplementary information

Supplementary data are available at Bioinformatics online.

P2864A novel mechanism of sinus node dysfunction: intergenic deletion between PITX2 and ANK2 disrupts chromatin structure in pacemaker cell differentiation

Background

Genetics of sinus node dysfunction (SND) remains poorly understood with few genes identified (HCN4, SCN5A, GJA5, KCNQ1, MYH6 and ANK2) and further genetic heterogeneity.

Purpose

We report two families with SND segregating with an intergenic deletion associated with long-range cis-regulatory elements (CREs) of PITX2.

Methods

We applied 30x whole genome sequencing (WGS) to two French families in which exome sequencing did not allow to detect the causing genes of SND and analyzed them bioinformatically.

Results

We first applied WGS to a 4-generation family presenting 16 patients with SND, atrial fibrillation and/or repolarization abnormalities and identified a 15-Kb deletion within a 16.5-cM linkage interval (Zmax = 5.9, θ = 0) on chromosome 4q25. We also applied WGS on a second family presenting 5 patients with SND and identified a 91-Kb deletion overlapping the initial deletion. This intergenic deletion, which was located between PITX2 and ANK2 and contained a binding motif of CCCTC-binding factor (CTCF), was segregating in all patients and was not found in 855 French control genomes. CTCF functions as a genome organizer mediating genomic interactions between genes and CREs. To clarify a possible regulatory function of the 1.5-Mb intergenic region containing the deletion, we interrogated 3 epigenetic databases. High-throughput chromosome conformation capture (Hi-C) (1) revealed that this intergenic region was placed in a topologically associating domain (TAD) containing PITX2. Chromatin interaction analysis by paired-end tag (ChIA-PET) of CTCF (2) suggested several possible functional chromatin loops in the TAD. Of them, a CTCF-mediated chromatin loop corresponding to a 300-Kb genomic region was clarified as a possible CRE using histone modification data from Roadmap. The region ranging between the intra-deletion CTCF motif and another distal motif was repressed in H3K27me3 profiles of fetal heart, while the region was activated in H3K27ac profiles of H1 ESC-derived mesendoderm.

Conclusion

This CRE may regulate expression of PITX2 and/or ANK2 in pacemaker cell differentiation. Ongoing experiments on patient's iPSC and transgenic mouse will further characterize the disease mechanism.

Sex specific associations in genome wide association analysis of renal cell carcinoma

Renal cell carcinoma (RCC) has an undisputed genetic component and a stable 2:1 male to female sex ratio in its incidence across populations, suggesting possible sexual dimorphism in its genetic susceptibility. We conducted the first sex-specific genome-wide association analysis of RCC for men (3227 cases, 4916 controls) and women (1992 cases, 3095 controls) of European ancestry from two RCC genome-wide scans and replicated the top findings using an additional series of men (2261 cases, 5852 controls) and women (1399 cases, 1575 controls) from two independent cohorts of European origin. Our study confirmed sex-specific associations for two known RCC risk loci at 14q24.2 (DPF3) and 2p21(EPAS1). We also identified two additional suggestive male-specific loci at 6q24.3 (SAMD5, male odds ratio (ORmale) = 0.83 [95% CI = 0.78-0.89], Pmale = 1.71 × 10-8 compared with female odds ratio (ORfemale) = 0.98 [95% CI = 0.90-1.07], Pfemale = 0.68) and 12q23.3 (intergenic, ORmale = 0.75 [95% CI = 0.68-0.83], Pmale = 1.59 × 10-8 compared with ORfemale = 0.93 [95% CI = 0.82-1.06], Pfemale = 0.21) that attained genome-wide significance in the joint meta-analysis. Herein, we provide evidence of sex-specific associations in RCC genetic susceptibility and advocate the necessity of larger genetic and genomic studies to unravel the endogenous causes of sex bias in sexually dimorphic traits and diseases like RCC.

Integrative and comparative genomic analyses identify clinically relevant pulmonary carcinoid groups and unveil the supra-carcinoids

The worldwide incidence of pulmonary carcinoids is increasing, but little is known about their molecular characteristics. Through machine learning and multi-omics factor analysis, we compare and contrast the genomic profiles of 116 pulmonary carcinoids (including 35 atypical), 75 large-cell neuroendocrine carcinomas (LCNEC), and 66 small-cell lung cancers. Here we report that the integrative analyses on 257 lung neuroendocrine neoplasms stratify atypical carcinoids into two prognostic groups with a 10-year overall survival of 88% and 27%, respectively. We identify therapeutically relevant molecular groups of pulmonary carcinoids, suggesting DLL3 and the immune system as candidate therapeutic targets; we confirm the value of OTP expression levels for the prognosis and diagnosis of these diseases, and we unveil the group of supra-carcinoids. This group comprises samples with carcinoid-like morphology yet the molecular and clinical features of the deadly LCNEC, further supporting the previously proposed molecular link between the low- and high-grade lung neuroendocrine neoplasms.

QuantumClone: clonal assessment of functional mutations in cancer based on a genotype-aware method for clonal reconstruction

Motivation

In cancer, clonal evolution is assessed based on information coming from single nucleotide variants and copy number alterations. Nonetheless, existing methods often fail to accurately combine information from both sources to truthfully reconstruct clonal populations in a given tumor sample or in a set of tumor samples coming from the same patient. Moreover, previously published methods detect clones from a single set of variants. As a result, compromises have to be done between stringent variant filtering [reducing dispersion in variant allele frequency estimates (VAFs)] and using all biologically relevant variants.

Results

We present a framework for defining cancer clones using most reliable variants of high depth of coverage and assigning functional mutations to the detected clones. The key element of our framework is QuantumClone, a method for variant clustering into clones based on VAFs, genotypes of corresponding regions and information about tumor purity. We validated QuantumClone and our framework on simulated data. We then applied our framework to whole genome sequencing data for 19 neuroblastoma trios each including constitutional, diagnosis and relapse samples. We confirmed an enrichment of damaging variants within such pathways as MAPK (mitogen-activated protein kinases), neuritogenesis, epithelial-mesenchymal transition, cell survival and DNA repair. Most pathways had more damaging variants in the expanding clones compared to shrinking ones, which can be explained by the increased total number of variants between these two populations. Functional mutational rate varied for ancestral clones and clones shrinking or expanding upon treatment, suggesting changes in clone selection mechanisms at different time points of tumor evolution.

Availability and implementation

Source code and binaries of the QuantumClone R package are freely available for download at https://CRAN.R-project.org/package=QuantumClone.

Contact

gudrun.schleiermacher@curie.fr or valentina.boeva@inserm.fr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Evidence of Austronesian Genetic Lineages in East Africa and South Arabia: Complex Dispersal from Madagascar and Southeast Asia

The Austronesian dispersal across the Indonesian Ocean to Madagascar and the Comoros has been well documented, but in an unexplained anomaly, few to no traces have been found of the Austronesian expansion in East Africa or the Arabian Peninsula. To revisit this peculiarity, we surveyed the Western Indian Ocean rim populations to identify potential Austronesian genetic ancestry. We generated full mitochondrial DNA genomes and genome-wide genotyping data for these individuals and compared them with the Banjar, the Indonesian source population of the westward Austronesian dispersal. We find strong support for Asian genetic contributions to maternal lineages and autosomal variation in modern day Somalia and Yemen. Surprisingly, this input reveals two apparently different geographic origins and timings of admixture for the Austronesian contact; one at a very early phase (likely associated with the early Austronesian dispersals), and a later movement dating to the end of nineteenth century. These Austronesian gene flows come, respectively, from Madagascar and directly from an unidentified location in Island Southeast Asia. This result reveals a far more complex dynamic of Austronesian dispersals through the Western Indian Ocean than has previously been understood and suggests that Austronesian movements within the Indian Ocean may have been part of a lengthy process, probably continuing well into the modern era.

A large-scale exome array analysis of venous thromboembolism

Although recent Genome-Wide Association Studies have identified novel associations for common variants, there has been no comprehensive exome-wide search for low-frequency variants that affect the risk of venous thromboembolism (VTE). We conducted a meta-analysis of 11 studies comprising 8,332 cases and 16,087 controls of European ancestry and 382 cases and 1,476 controls of African American ancestry genotyped with the Illumina HumanExome BeadChip. We used the seqMeta package in R to conduct single variant and gene-based rare variant tests. In the single variant analysis, we limited our analysis to the 64,794 variants with at least 40 minor alleles across studies (minor allele frequency [MAF] ~0.08%). We confirmed associations with previously identified VTE loci, including ABO, F5, F11, and FGA. After adjusting for multiple testing, we observed no novel significant findings in single variant or gene-based analysis. Given our sample size, we had greater than 80% power to detect minimum odds ratios greater than 1.5 and 1.8 for a single variant with MAF of 0.01 and 0.005, respectively. Larger studies and sequence data may be needed to identify novel low-frequency and rare variants associated with VTE risk.

Prediction of treatment (tx)-induced fatigue in breast cancer (BC) patients (pts) using machine learning on genome-wide association (GWAS) data in the prospective CANTO cohort

11515

Background: Many BC survivors report fatigue. The relevant genomic correlates of fatigue after BC are not well understood. We applied a previously validated machine learning methodology (Oh 2017) to GWAS data to identify biological correlates of fatigue induced after tx. Methods: We analyzed 3825 BC pts with GWAS data (Illumina InfiniumExome24 v 1.1) from the CANTO study (NCT01993498). The outcome of this study was post-tx fatigue 1 year after the end of primary chemotherapy/radiotherapy/surgery using the EORTC C30 fatigue subscale (overall fatigue) and the EORTC FA 12 fatigue domains (physical/emotional/cognitive). For each domain, we limited the study group to those with zero baseline fatigue and defined severe fatigue change as score increase above the third quartile. We tested univariate correlations between severe fatigue in each domain and 496539 SNPs as well as relevant clinical variables. The machine learning prediction model based on preconditioning random forest regression (PRFR) (Oh et al., 2017), was then built using the SNPs with ancestry adjusted univariate p-value < 0.001 and clinical variables with Bonferroni adjusted p-value < 0.05. The model was validated in a holdout subset of the cohort. Gene set enrichment analysis (GSEA) was performed using MetaCore to identify key biological correlates relevant to tx-induced fatigue. Results: Distinct results were found by fatigue domain (table). GSEA showed that the cognitive fatigue model SNPs included biomarkers for cognitive disorders (p = 1.6 x 10-12) and glutamatergic synaptic transmission (p = 1.6 x 10-8). Conclusions: A SNP based model had differential performance by fatigue domain, with a potential genetic role on risk and biology for tx induced cognitive fatigue. Further research to explore biomarkers of tx induced fatigue are needed. [Table: see text]

Prognosis value of a genetic score based on germline genetic variants in a prospective cohort of early triple-negative breast cancer patients

1529

Background: Triple-negative breast cancers (TNBC) are a heterogeneous group of tumors with poor outcome. In this study, the association between germline genetic variants and invasive disease-free survival (iDFS) was analyzed in TNBC patients. Methods: A genome wide-association study (GWAS) aimed to identify variants (single nucleotide polymorphisms – SNPs) associated with prognosis in 1121 patients with TNBC in the SIGNAL prospective cohort. Associations between gene variants and iDFS were assessed in univariate Cox regression models. Variants were combined in a score to identify risk categories. A prognostic model based on breast cancer stage and genetic variants was estimated using a multivariate Cox regression. Interaction between stage and genetic score was tested. Discrimination of the model was assessed by the Harrell’s C statistic and internal validity by bootstrap method. Results: The characteristics of the 1121 patients were representative of a population with early TNBC. Four SNPs on chromosomes 9 and 2 were found significantly associated to iDFS in univariate Cox models. Homozygous status for the most frequent allele was associated with poorer iDFS for two SNPs and this status was present in 50% and 57% of the population. For the two other SNPs, the most frequent allele was associated with more favorable iDFS. Three prognostic categories were derived from the genetic score. The following table presents the results from the multivariate Cox model including genetic score and disease stage. Clinical trial information: RECF1098. Conclusions: In a prospective cohort of 1121 patients with early TNBC, 4 genetic variants (SNPs) were associated with iDFS. A score involving SNPs provided similar prognostic indications as breast cancer stages. A search assessing the function and the role of the involved genes is ongoing.[Table: see text]

Genetic meta-analysis of diagnosed Alzheimer's disease identifies new risk loci and implicates Aβ, tau, immunity and lipid processing

Risk for late-onset Alzheimer's disease (LOAD), the most prevalent dementia, is partially driven by genetics. To identify LOAD risk loci, we performed a large genome-wide association meta-analysis of clinically diagnosed LOAD (94,437 individuals). We confirm 20 previous LOAD risk loci and identify five new genome-wide loci (IQCK, ACE, ADAM10, ADAMTS1, and WWOX), two of which (ADAM10, ACE) were identified in a recent genome-wide association (GWAS)-by-familial-proxy of Alzheimer's or dementia. Fine-mapping of the human leukocyte antigen (HLA) region confirms the neurological and immune-mediated disease haplotype HLA-DR15 as a risk factor for LOAD. Pathway analysis implicates immunity, lipid metabolism, tau binding proteins, and amyloid precursor protein (APP) metabolism, showing that genetic variants affecting APP and Aβ processing are associated not only with early-onset autosomal dominant Alzheimer's disease but also with LOAD. Analyses of risk genes and pathways show enrichment for rare variants (P = 1.32 × 10-7), indicating that additional rare variants remain to be identified. We also identify important genetic correlations between LOAD and traits such as family history of dementia and education.

The influence of obesity-related factors in the etiology of renal cell carcinoma-A mendelian randomization study

BACKGROUND

Several obesity-related factors have been associated with renal cell carcinoma (RCC), but it is unclear which individual factors directly influence risk. We addressed this question using genetic markers as proxies for putative risk factors and evaluated their relation to RCC risk in a mendelian randomization (MR) framework. This methodology limits bias due to confounding and is not affected by reverse causation.

METHODS AND FINDINGS

Genetic markers associated with obesity measures, blood pressure, lipids, type 2 diabetes, insulin, and glucose were initially identified as instrumental variables, and their association with RCC risk was subsequently evaluated in a genome-wide association study (GWAS) of 10,784 RCC patients and 20,406 control participants in a 2-sample MR framework. The effect on RCC risk was estimated by calculating odds ratios (ORSD) for a standard deviation (SD) increment in each risk factor. The MR analysis indicated that higher body mass index increases the risk of RCC (ORSD: 1.56, 95% confidence interval [CI] 1.44-1.70), with comparable results for waist-to-hip ratio (ORSD: 1.63, 95% CI 1.40-1.90) and body fat percentage (ORSD: 1.66, 95% CI 1.44-1.90). This analysis further indicated that higher fasting insulin (ORSD: 1.82, 95% CI 1.30-2.55) and diastolic blood pressure (DBP; ORSD: 1.28, 95% CI 1.11-1.47), but not systolic blood pressure (ORSD: 0.98, 95% CI 0.84-1.14), increase the risk for RCC. No association with RCC risk was seen for lipids, overall type 2 diabetes, or fasting glucose.

CONCLUSIONS

This study provides novel evidence for an etiological role of insulin in RCC, as well as confirmatory evidence that obesity and DBP influence RCC risk.

The Transcription Factor STAT6 Mediates Direct Repression of Inflammatory Enhancers and Limits Activation of Alternatively Polarized Macrophages

The molecular basis of signal-dependent transcriptional activation has been extensively studied in macrophage polarization, but our understanding remains limited regarding the molecular determinants of repression. Here we show that IL-4-activated STAT6 transcription factor is required for the direct transcriptional repression of a large number of genes during in vitro and in vivo alternative macrophage polarization. Repression results in decreased lineage-determining transcription factor, p300, and RNA polymerase II binding followed by reduced enhancer RNA expression, H3K27 acetylation, and chromatin accessibility. The repressor function of STAT6 is HDAC3 dependent on a subset of IL-4-repressed genes. In addition, STAT6-repressed enhancers show extensive overlap with the NF-κB p65 cistrome and exhibit decreased responsiveness to lipopolysaccharide after IL-4 stimulus on a subset of genes. As a consequence, macrophages exhibit diminished inflammasome activation, decreased IL-1β production, and pyroptosis. Thus, the IL-4-STAT6 signaling pathway establishes an alternative polarization-specific epigenenomic signature resulting in dampened macrophage responsiveness to inflammatory stimuli.

•

IL-4-activated STAT6 acts as a transcriptional repressor in macrophages

•

IL-4-STAT6-repressed enhancers associate with reduced LDTF and p300 binding

•

Inflammatory responsiveness of the IL-4-repressed enhancers is attenuated

•

IL-4 limits the LPS-induced inflammasome activation, IL-1β production, and pyroptosis

RXR heterodimers orchestrate transcriptional control of neurogenesis and cell fate specification

Retinoid X Receptors (RXRs) are unique and enigmatic members of the nuclear receptor (NR) family with extensive and complex biological functions in cellular differentiation. On the one hand, RXRs through permissive heterodimerization with other NRs are able to integrate multiple lipid signaling pathways and are believed to play a central role to coordinate the development of the central nervous system. On the other hand, RXRs may have heterodimer-independent functions as well. Therefore, a more RXR-centric analysis is warranted to identify its genomic binding sites and regulated gene networks, which are orchestrating the earliest events in neuronal differentiation. Recently developed genome-wide approaches allow systematic analyses of the RXR-driven neural differentiation. Here we applied next generation sequencing-based methodology to track the dynamic redistribution of the RXR cistrome along the path of embryonic stem cell to glutamatergic neuron differentiation. We identified Retinoic Acid Receptor (RAR) and Liver X Receptor (LXR) as dominant heterodimeric partners of RXR in these cellular stages. Our data presented here characterize the RAR:RXR and LXR:RXR-mediated transcriptional program in embryonic stem cells, neural progenitors and terminally differentiated neurons. Considering the growing evidence for dysregulated RXR-mediated signaling in neurodegenerative disorders, such as Alzheimer's Disease or Amyotrophic Lateral Sclerosis, the data presented here will be also a valuable resource for the field of neuro(patho)biology.

Extensive and functional overlap of the STAT6 and RXR cistromes in the active enhancer repertoire of human CD14+ monocyte derived differentiating macrophages

Macrophages are able to differentiate into classically polarized (M1) or alternatively polarized (M2) states upon encountering pro-inflammatory cytokines such as interferon (IFN) γ or anti-inflammatory cytokines such as interleukin (IL) -4/IL-13, respectively. Moreover, macrophages are known to regulate lipid metabolism via multiple members of the nuclear hormone receptor family, including the retinoid X receptors (RXR). It has been also documented that cytokines are able to modulate macrophage responses to lipid signals but the nature of these interactions and the underlying mechanisms of these processes especially at the level of the chromatinized genome are not well understood. Previous work from our laboratory suggested that STAT6 is a facilitator of nuclear receptor mediated transcriptional activity acting at the genome level. This prompted us to investigate genome-wide DNA binding events and the development of cistromes in human CD14+ monocyte-derived macrophages upon exposure to IL-4. We determined the impact of IL-4 on the PU.1, RXR and STAT6 cistromes within the active enhancer regions marked by H3K27-acetylation using chromatin immunoprecipitation followed by deep sequencing and integrated bioinformatics analyses. We found that about 2/3rd of the IL-4 induced STAT6 peaks co-localized with RXR peaks. These STAT6/RXR co-peaks differed at least in part from the non-overlapping RXR peaks regarding the most enriched de novo transcription factor binding motifs. Interestingly, RXR-binding was not regulated at the STAT6/RXR co-bound enhancers following IL-4 stimulation, but differential enhancer interactions were observed between the IL-4/STAT6 and RXR signaling pathways acting in a gene selective manner. Our results suggest that there is a novel, so far uncharacterized cistromic crosstalk between RXR and STAT6 that is likely to contribute to the formation of the active enhancer repertoire, transcriptome and differential signal-specific gene regulation of polarized macrophages.

Strong selection during the last millennium for African ancestry in the admixed population of Madagascar

While admixed populations offer a unique opportunity to detect selection, the admixture in most of the studied populations occurred too recently to produce conclusive signals. By contrast, Malagasy populations originate from admixture between Asian and African populations that occurred ~27 generations ago, providing power to detect selection. We analyze local ancestry across the genomes of 700 Malagasy and identify a strong signal of recent positive selection, with an estimated selection coefficient >0.2. The selection is for African ancestry and affects 25% of chromosome 1, including the Duffy blood group gene. The null allele at this gene provides resistance to Plasmodium vivax malaria, and previous studies have suggested positive selection for this allele in the Malagasy population. This selection event also influences numerous other genes implicated in immunity, cardiovascular diseases, and asthma and decreases the Asian ancestry genome-wide by 10%, illustrating the role played by selection in recent human history.

Genetic variants of RNASE3 (ECP) and susceptibility to severe malaria in Senegalese population

Background

Severe forms of malaria (SM) are an outcome of Plasmodium falciparum infection and can cause death especially in children under 4 years of age. RNASE3 (ECP) has been identified as an inhibitor of Plasmodium parasites growth in vitro, and genetic analysis in hospitalized Ghanaian subjects has revealed the RNASE3 +371G/C (rs2073342) polymorphism as a susceptibility factor for cerebral malaria. The +371 C allele results in an Arg/Thr mutation that abolishes the cytotoxic activity of the ECP protein. The present study aims to investigate RNASE3 gene polymorphisms and their putative link to severe malaria in a malaria cohort from Senegal.

Methods/results

Patients enrolled from hospitals were classified as having either uncomplicated (UM) or severe malaria (SM). The analysis of the RNASE3 gene polymorphisms was performed in 241 subjects: 178 falciparum infected (96 SM, 82 UM) and 63 non-infected subjects as population control group (CTR). Six frequent SNPs (MAF > 3%) were identified, and one SNP was associated with malaria severity by performing a logistic regression analysis SM vs.UM: RNASE3 +499G/C (rs2233860) under age, sex as covariates and HbS/HbC polymorphisms adjustment (p = 0.003, OR 0.43, CI 95% 0.20–0.92). The polymorphisms: +371G/C (rs2073342), +499G/C (rs2233860) and +577A/T (rs8019343) defined a haplotype risk (G-G-T) for malaria severity (Fisher exact test, p = 0.03) (OR 4.1, IC 95% (1.1–14.9).

Conclusion

In addition to the previously described association of +371G/C polymorphism in Ghanaians cohort, the RNASE3 +499G/C polymorphism was associated with susceptibility to SM in a Senegalese population. The haplotype +371G/+499G/+577T defined by RNASE3 polymorphisms was associated with severity. The genetic association identified independently in the Senegalese population provide additional evidence of a role of RNASE3 (ECP) in malaria severity.

The Comoros Show the Earliest Austronesian Gene Flow into the Swahili Corridor

At the dawn of the second millennium, the expansion of the Indian Ocean trading network aligned with the emergence of an outward-oriented community along the East African coast to create a cosmopolitan cultural and trading zone known as the Swahili Corridor. On the basis of analyses of new genome-wide genotyping data and uniparental data in 276 individuals from coastal Kenya and the Comoros islands, along with large-scale genetic datasets from the Indian Ocean rim, we reconstruct historical population dynamics to show that the Swahili Corridor is largely an eastern Bantu genetic continuum. Limited gene flows from the Middle East can be seen in Swahili and Comorian populations at dates corresponding to historically documented contacts. However, the main admixture event in southern insular populations, particularly Comorian and Malagasy groups, occurred with individuals from Island Southeast Asia as early as the 8^th century, reflecting an earlier dispersal from this region. Remarkably, our results support recent archaeological and linguistic evidence-based suggestions that the Comoros archipelago was the earliest location of contact between Austronesian and African populations in the Swahili Corridor.

Genetic Variants Related to Longer Telomere Length are Associated with Increased Risk of Renal Cell Carcinoma

BACKGROUND

Relative telomere length in peripheral blood leukocytes has been evaluated as a potential biomarker for renal cell carcinoma (RCC) risk in several studies, with conflicting findings.

OBJECTIVE

We performed an analysis of genetic variants associated with leukocyte telomere length to assess the relationship between telomere length and RCC risk using Mendelian randomization, an approach unaffected by biases from temporal variability and reverse causation that might have affected earlier investigations.

DESIGN, SETTING, AND PARTICIPANTS

Genotypes from nine telomere length-associated variants for 10 784 cases and 20 406 cancer-free controls from six genome-wide association studies (GWAS) of RCC were aggregated into a weighted genetic risk score (GRS) predictive of leukocyte telomere length.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Odds ratios (ORs) relating the GRS and RCC risk were computed in individual GWAS datasets and combined by meta-analysis.

RESULTS AND LIMITATIONS

Longer genetically inferred telomere length was associated with an increased risk of RCC (OR=2.07 per predicted kilobase increase, 95% confidence interval [CI]:=1.70-2.53, p<0.0001). As a sensitivity analysis, we excluded two telomere length variants in linkage disequilibrium (R2>0.5) with GWAS-identified RCC risk variants (rs10936599 and rs9420907) from the telomere length GRS; despite this exclusion, a statistically significant association between the GRS and RCC risk persisted (OR=1.73, 95% CI=1.36-2.21, p<0.0001). Exploratory analyses for individual histologic subtypes suggested comparable associations with the telomere length GRS for clear cell (N=5573, OR=1.93, 95% CI=1.50-2.49, p<0.0001), papillary (N=573, OR=1.96, 95% CI=1.01-3.81, p=0.046), and chromophobe RCC (N=203, OR=2.37, 95% CI=0.78-7.17, p=0.13).

CONCLUSIONS

Our investigation adds to the growing body of evidence indicating some aspect of longer telomere length is important for RCC risk.

PATIENT SUMMARY

Telomeres are segments of DNA at chromosome ends that maintain chromosomal stability. Our study investigated the relationship between genetic variants associated with telomere length and renal cell carcinoma risk. We found evidence suggesting individuals with inherited predisposition to longer telomere length are at increased risk of developing renal cell carcinoma.

Genome-wide association study identifies multiple risk loci for renal cell carcinoma

Previous genome-wide association studies (GWAS) have identified six risk loci for renal cell carcinoma (RCC). We conducted a meta-analysis of two new scans of 5,198 cases and 7,331 controls together with four existing scans, totalling 10,784 cases and 20,406 controls of European ancestry. Twenty-four loci were tested in an additional 3,182 cases and 6,301 controls. We confirm the six known RCC risk loci and identify seven new loci at 1p32.3 (rs4381241, P=3.1 × 10-10), 3p22.1 (rs67311347, P=2.5 × 10-8), 3q26.2 (rs10936602, P=8.8 × 10-9), 8p21.3 (rs2241261, P=5.8 × 10-9), 10q24.33-q25.1 (rs11813268, P=3.9 × 10-8), 11q22.3 (rs74911261, P=2.1 × 10-10) and 14q24.2 (rs4903064, P=2.2 × 10-24). Expression quantitative trait analyses suggest plausible candidate genes at these regions that may contribute to RCC susceptibility.

Protein S Heerlen mutation heterozygosity is associated with venous thrombosis risk

Hereditary Protein S (PS) deficiency is a rare coagulation disorder associated with an increased risk of venous thrombosis (VT). The PS Heerlen (PSH) mutation is a rare S501P mutation that was initially considered to be a neutral polymorphism. However, it has been later shown that PSH has a reduced half-life in vivo which may explain the association of PSH heterozygosity with mildly reduced levels of plasma free PS (FPS). Whether the risk of VT is increased in PSH carriers remains unknown. We analyzed the association of PSH (rs121918472 A/G) with VT in 4,173 VT patients and 5,970 healthy individuals from four independent case-control studies. Quantitative determination of FPS levels was performed in a subsample of 1257 VT patients. In the investigated populations, the AG genotype was associated with an increased VT risk of 6.57 [4.06–10.64] (p = 1.73 10⁻¹⁴). In VT patients in whom PS deficiency was excluded, plasma FPS levels were significantly lower in individuals with PSH when compared to those without [72 + 13 vs 91 + 21 UI/dL; p = 1.86 10⁻⁶, mean + SD for PSH carriers (n = 21) or controls (n = 1236) respectively]. We provide strong evidence that the rare PSH variant is associated with VT in unselected individuals.

Immortalized human myotonic dystrophy muscle cell lines to assess therapeutic compounds

Myotonic dystrophy type 1 (DM1) and type 2 (DM2) are autosomal dominant neuromuscular diseases caused by microsatellite expansions and belong to the family of RNA-dominant disorders. Availability of cellular models in which the DM mutation is expressed within its natural context is essential to facilitate efforts to identify new therapeutic compounds. Here, we generated immortalized DM1 and DM2 human muscle cell lines that display nuclear RNA aggregates of expanded repeats, a hallmark of myotonic dystrophy. Selected clones of DM1 and DM2 immortalized myoblasts behave as parental primary myoblasts with a reduced fusion capacity of immortalized DM1 myoblasts when compared with control and DM2 cells. Alternative splicing defects were observed in differentiated DM1 muscle cell lines, but not in DM2 lines. Splicing alterations did not result from differentiation delay because similar changes were found in immortalized DM1 transdifferentiated fibroblasts in which myogenic differentiation has been forced by overexpression of MYOD1. As a proof-of-concept, we show that antisense approaches alleviate disease-associated defects, and an RNA-seq analysis confirmed that the vast majority of mis-spliced events in immortalized DM1 muscle cells were affected by antisense treatment, with half of them significantly rescued in treated DM1 cells. Immortalized DM1 muscle cell lines displaying characteristic disease-associated molecular features such as nuclear RNA aggregates and splicing defects can be used as robust readouts for the screening of therapeutic compounds. Therefore, immortalized DM1 and DM2 muscle cell lines represent new models and tools to investigate molecular pathophysiological mechanisms and evaluate the in vitro effects of compounds on RNA toxicity associated with myotonic dystrophy mutations.

Summary: Myotonic dystrophy muscle cell models displaying characteristic disease-associated molecular features can be used to investigate molecular pathophysiological mechanisms and evaluate therapeutic approaches.

OCT4 Acts as an Integrator of Pluripotency and Signal-Induced Differentiation

Cell type specification relies on the capacity of undifferentiated cells to properly respond to specific differentiation-inducing signals. Using genomic approaches along with loss- and gain-of-function genetic models, we identified OCT4-dependent mechanisms that provide embryonic stem cells with the means to customize their response to external cues. OCT4 binds a large set of low-accessible genomic regions. At these sites, OCT4 is required for proper enhancer and gene activation by recruiting co-regulators and RAR:RXR or β-catenin, suggesting an unexpected collaboration between the lineage-determining transcription factor and these differentiation-initiating, signal-dependent transcription factors. As a proof of concept, we demonstrate that overexpression of OCT4 in a kidney cell line is sufficient for signal-dependent activation of otherwise unresponsive genes in these cells. Our results uncover OCT4 as an integral and necessary component of signal-regulated transcriptional processes required for tissue-specific responses.

Deep sequencing of the Trypanosoma cruzi GP63 surface proteases reveals diversity and diversifying selection among chronic and congenital Chagas disease patients

Background

Chagas disease results from infection with the diploid protozoan parasite Trypanosoma cruzi. T. cruzi is highly genetically diverse, and multiclonal infections in individual hosts are common, but little studied. In this study, we explore T. cruzi infection multiclonality in the context of age, sex and clinical profile among a cohort of chronic patients, as well as paired congenital cases from Cochabamba, Bolivia and Goias, Brazil using amplicon deep sequencing technology.

Methodology/ Principal Findings

A 450bp fragment of the trypomastigote TcGP63I surface protease gene was amplified and sequenced across 70 chronic and 22 congenital cases on the Illumina MiSeq platform. In addition, a second, mitochondrial target—ND5—was sequenced across the same cohort of cases. Several million reads were generated, and sequencing read depths were normalized within patient cohorts (Goias chronic, n = 43, Goias congenital n = 2, Bolivia chronic, n = 27; Bolivia congenital, n = 20), Among chronic cases, analyses of variance indicated no clear correlation between intra-host sequence diversity and age, sex or symptoms, while principal coordinate analyses showed no clustering by symptoms between patients. Between congenital pairs, we found evidence for the transmission of multiple sequence types from mother to infant, as well as widespread instances of novel genotypes in infants. Finally, non-synonymous to synonymous (dn:ds) nucleotide substitution ratios among sequences of TcGP63Ia and TcGP63Ib subfamilies within each cohort provided powerful evidence of strong diversifying selection at this locus.

Conclusions/Significance

Our results shed light on the diversity of parasite DTUs within each patient, as well as the extent to which parasite strains pass between mother and foetus in congenital cases. Although we were unable to find any evidence that parasite diversity accumulates with age in our study cohorts, putative diversifying selection within members of the TcGP63I gene family suggests a link between genetic diversity within this gene family and survival in the mammalian host.

Trypanosoma cruzi, the causal agent of Chagas disease in Latin America, infects several million people in some of the most economically deprived regions of Latin America. T. cruzi infection is lifelong and has a variable prognosis: some patients never exhibit symptoms while others experience debilitating and fatal complications. Available data suggest that parasite genetic diversity within and among disease foci can be exceedingly high. However, little is know about the frequency of multiple genotype infections in humans, as well as their distribution among different age classes and possible impact on disease outcome. In this study we develop a next generation amplicon deep sequencing approach to profile parasite diversity within chronic Chagas Disease patients from Bolivia and Brazil. We were also able to compare parasite genetic diversity present in eleven congenitally infants with parasite genetic diversity present in their mothers. We did not detect any specific association between the number and diversity of parasite genotypes in each patient with their age, sex or disease status. We were, however, able to detect the transmission of multiple parasite genotypes between mother and foetus. Furthermore, we also detected powerful evidence for natural selection at the antigenic locus we targeted, suggesting a possible interaction with the host immune system.

Meta-analysis of 65,734 individuals identifies TSPAN15 and SLC44A2 as two susceptibility loci for venous thromboembolism

Venous thromboembolism (VTE), the third leading cause of cardiovascular mortality, is a complex thrombotic disorder with environmental and genetic determinants. Although several genetic variants have been found associated with VTE, they explain a minor proportion of VTE risk in cases. We undertook a meta-analysis of genome-wide association studies (GWASs) to identify additional VTE susceptibility genes. Twelve GWASs totaling 7,507 VTE case subjects and 52,632 control subjects formed our discovery stage where 6,751,884 SNPs were tested for association with VTE. Nine loci reached the genome-wide significance level of 5 × 10(-8) including six already known to associate with VTE (ABO, F2, F5, F11, FGG, and PROCR) and three unsuspected loci. SNPs mapping to these latter were selected for replication in three independent case-control studies totaling 3,009 VTE-affected individuals and 2,586 control subjects. This strategy led to the identification and replication of two VTE-associated loci, TSPAN15 and SLC44A2, with lead risk alleles associated with odds ratio for disease of 1.31 (p = 1.67 × 10(-16)) and 1.21 (p = 2.75 × 10(-15)), respectively. The lead SNP at the TSPAN15 locus is the intronic rs78707713 and the lead SLC44A2 SNP is the non-synonymous rs2288904 previously shown to associate with transfusion-related acute lung injury. We further showed that these two variants did not associate with known hemostatic plasma markers. TSPAN15 and SLC44A2 do not belong to conventional pathways for thrombosis and have not been associated to other cardiovascular diseases nor related quantitative biomarkers. Our findings uncovered unexpected actors of VTE etiology and pave the way for novel mechanistic concepts of VTE pathophysiology.

A small-molecule screen for enhanced homing of systemically infused cells

Poor homing of systemically infused cells to disease sites may limit the success of exogenous cell-based therapy. In this study, we screened 9,000 signal-transduction modulators to identify hits that increase mesenchymal stromal cell (MSC) surface expression of homing ligands that bind to intercellular adhesion molecule 1 (ICAM-1), such as CD11a. Pretreatment of MSCs with Ro-31-8425, an identified hit from this screen, increased MSC firm adhesion to an ICAM-1-coated substrate in vitro and enabled targeted delivery of systemically administered MSCs to inflamed sites in vivo in a CD11a- (and other ICAM-1-binding domains)-dependent manner. This resulted in a heightened anti-inflammatory response. This represents a new strategy for engineering cell homing to enhance therapeutic efficacy and validates CD11a and ICAM-1 as potential targets. Altogether, this multi-step screening process may significantly improve clinical outcomes of cell-based therapies.

Association between the Thr325Ile polymorphism of the thrombin-activatable fibrinolysis inhibitor and stroke in the Ludwigshafen Risk and Cardiovascular Health Study

The thrombin-activatable fibrinolysis inhibitor (TAFI) is a key mediator in the regulation of endogenous fibrinolysis, down-regulating clot lysis by degrading the C-terminal lysine residues from fibrin, which are important for binding and activating plasminogen. Elevated TAFI antigen levels have been suggested to be associated with promoter variants and the Ala147Thr polymorphism; increased TAFI stability and antifibrinolytic potential instead have been associated with the Thr325Ile polymorphism. We investigated the influence of these two polymorphisms on cardiovascular and thrombotic events in patients of the Ludwigshafen Risk and Cardiovascular Health (LURIC) study. The LURIC study is a prospective cohort study comprising more than 3,300 patients aimed at identifying biochemical and genetic markers for metabolic and cardiovascular diseases. We demonstrate that the Ile/Ile genotype at position 325 of TAFI associates with the incidence of stroke and the age at onset of first stroke in patients of the LURIC cohort. Both the incidence of stroke and the risk of a premature event are higher in TAFI Ile325Ile patients with predisposing risk factors for thrombotic events such as diabetes mellitus, myocardial infarction or hypertension, alone or in combination. In contrast, no significant association was identified for the TAFI Ala147Thr polymorphism. The robust association of the TAFI Thr325Ile polymorphism with the incidence and the age at onset of first stroke strongly suggests a key role for TAFI in the pathogenetic mechanism of stroke.

Association between the MnSOD Ala-9Val polymorphism and development of schizophrenia and abnormal involuntary movements in the Xhosa population

Reactive oxygen species (ROS)-mediated damage has been hypothesized to play a role in the development and poor outcome of schizophrenia, as well as the development of neuroleptic-induced abnormal involuntary movements. Recently, the functional polymorphism (Ala-9Val) in the manganese superoxide dismutase (MnSOD) gene (part of the antioxidant defense mechanism) was found to be associated with schizophrenia in a Turkish population. This study was aimed at replicating this finding in a Xhosa population. In addition, the role of Ala-9Val in abnormal involuntary movement and tardive dyskinesia development in the Xhosa population was also investigated. The schizophrenic patient group (n=286) and a healthy control group (n=243) were genotyped for the Ala-9Val polymorphism using heteroduplex-single stranded conformational polymorphism (HEX-SSCP) analysis. No significant difference in genotype or allele frequency could be observed between the schizophrenia and control group (P=0.294 and P=0.528 respectively). In addition no association could be found between the polymorphism and symptom severity (SANS and SAPS). The Xhosa schizophrenia patient group with abnormal involuntary movements (n=54) and a subgroup with tardive dyskinesia (n=30) was found to significantly differ in Ala-9Val genotype frequency (P=0.008 and P=0.011 respectively) compared to the Xhosa schizophrenia patient group without abnormal involuntary movements (n=204). However, no significant difference was found for the allele frequencies (P=0.955 and P=0.161). Further, using ANCOVA no association was found between AIMS score and genotype in the group with abnormal involuntary movements (P=0.1234). However, in the patient group with tardive dyskinesia an association was observed between genotype and AIMS score (P=0.0365). These results do not support a major role of the MnSOD Ala-9Val polymorphism in the development of schizophrenia or symptom severity in the Xhosa population. Yet it seems to be involved in the development of abnormal involuntary movements and tardive dyskinesia and may even modulate the severity of tardive dyskinesia.

Selection of SNP subsets for association studies in candidate genes: comparison of the power of different strategies to detect single disease susceptibility locus effects

BACKGROUND

The recent advances in genotyping and molecular techniques have greatly increased the knowledge of the human genome structure. Millions of polymorphisms are reported and freely available in public databases. As a result, there is now a need to identify among all these data, the relevant markers for genetic association studies. Recently, several methods have been published to select subsets of markers, usually Single Nucleotide Polymorphisms (SNPs), that best represent genetic polymorphisms in the studied candidate gene or region.

RESULTS

In this paper, we compared four of these selection methods, two based on haplotype information and two based on pairwise linkage disequilibrium (LD). The methods were applied to the genotype data on twenty genes with different patterns of LD and different numbers of SNPs. A measure of the efficiency of the different methods to select SNPs was obtained by comparing, for each gene and under several single disease susceptibility models, the power to detect an association that will be achieved with the selected SNP subsets.

CONCLUSION

None of the four selection methods stands out systematically from the others. Methods based on pairwise LD information turn out to be the most interesting methods in a context of association study in candidate gene. In a context where the number of SNPs to be tested in a given region needs to be more limited, as in large-scale studies or wide genome scans, one of the two methods based on haplotype information, would be more suitable.

Positive and negative symptoms in affected sib pairs with schizophrenia: implications for genetic studies in an African Xhosa sample

Careful phenotyping and the identification of subtypes of schizophrenia can contribute significantly to the success of genetic studies in schizophrenia. The phenomenology of schizophrenia in affected sib pairs has been well-described in Caucasian populations, however a paucity of data exists for African populations. This study therefore investigated symptom dimensions in a sizeable group of affected Xhosa sib pairs as a means of evaluating the role of shared familial factors in the psychosis of schizophrenia. Five hundred and thirteen participants were interviewed with the Diagnostic Interview for Genetic Studies (DIGS), which included the Schedules for the Assessment of Negative and Positive symptoms (SANS/SAPS). One hundred and four sib pairs were then extracted (N = 208) for analysis of concordance for lifetime psychotic symptoms and an exploratory factor analysis of the SANS/SAPS. Concordance analysis of life-time symptoms indicated a significant concordance for olfactory hallucinations, persecutory delusions, jealousy, somatic, reference and control delusions as well as thought insertion and withdrawal. The factor analysis of the global scores of the SAPS and SANS revealed a five factor best-fit model and accounted for 92.5% of variance. The factors included a negative symptom factor, a positive symptom factor, a positive thought disorder and a bizarre behaviour component. The core symptomatology of schizophrenia in this sib pair sample was similar to that reported in Caucasian populations with the exception of higher rates of auditory hallucinations and delusions of persecution. In summary therefore; although the factor analysis only supported the concept of the universality of psychotic symptoms in schizophrenia, the concordance analysis of these symptoms did reveal hallucinations as well as delusions of control as possible candidates relevant for future research into genotype-phenotype relationships.