Identification of potential common genetic modifiers of neurofibromas: a genome-wide association study in 1333 patients with neurofibromatosis type 1

BACKGROUND

Neurofibromatosis type 1 (NF1) is characterized by the highly variable and unpredictable development of benign peripheral nerve sheath tumours: cutaneous (cNFs), subcutaneous (scNFs) and plexiform (pNFs) neurofibromas.

OBJECTIVES

To identify neurofibroma modifier genes, in order to develop a database of patients with NF1.

METHODS

All patients were phenotypically evaluated by a medical practitioner using a standardized questionnaire and the causal NF1 variant identified. We enrolled 1333 patients with NF1 who were genotyped for > 7 million common variants.

RESULTS

A genome-wide association case-only study identified a significant association with 9q21.33 in the pNF phenotype in the discovery cohort. Twelve, three and four regions suggestive of association at the P ≤ 1 × 10-6 threshold were identified for pNFs, cNFs and scNFs, respectively. Evidence of replication was observed for 4, 2 and 6 loci, including 168 candidate modifier protein-coding genes. Among the candidate modifier genes, some were implicated in the RAS-mitogen-activated protein kinase pathway, cell-cycle control and myelination. Using an original CRISPR/Cas9-based functional assay, we confirmed GAS1 and SPRED2 as pNF and scNF candidate modifiers, as their inactivation specifically affected NF1-mutant Schwann cell growth.

CONCLUSIONS

Our study may shed new light on the pathogenesis of NF1-associated neurofibromas and will, hopefully, contribute to the development of personalized care for patients with this deleterious and life-threatening condition.

Layer myocardial strain is the most heritable echocardiographic trait

AIMS

Myocardial deformation assessed by strain analysis represents a significant advancement in our assessment of cardiac mechanics. However, whether this variable is genetically heritable or whether all/most of its variability is related to environmental factors is currently unknown. We sought to determine the heritability of echocardiographically determined cardiac mechanics indices in a population setting.

METHODS AND RESULTS

A total of 1357 initially healthy subjects (women 51.6%; 48.2 ± 14.1 years) were included in this study from 20-year follow-up after the fourth visit of the longitudinal familial STANISLAS cohort (Lorraine, France). Data were acquired using state-of-the-art cardiac ultrasound equipment, using acquisition and measurement protocols recommended by the EACVI (European Association of Cardiovascular Imaging)/ASE (American Society of Echocardiography)/Industry Task Force. Layer-specific global longitudinal strain (GLS) and global circumferential strain (full-wall, subendocardial, and subepicardial) and conventional structural and functional cardiac parameters and their potential heritability were assessed using restricted maximum likelihood analysis, with genetic relatedness matrix calculated from genome-wide association data. Indices of longitudinal/circumferential myocardial function and left ventricular (LV) ejection fraction had low heritability (ranging from 10% to 20%). Diastolic and standard LV function parameters had moderate heritability (ranging from 20% to 30%) except for end-systolic and end-diastolic volumes (30% and 45%, respectively). In contrast, global longitudinal subendocardial strain (GLSEndo)/global longitudinal subepicardial strain (GLSEpi) ratio had a high level of heritability (65%). Except for GLSEndo/GLSEpi ratio, a large percentage of variance remained unexplained (>50%).

CONCLUSIONS

In our population cohort, GLSEndo/GLSEpi ratio had a high level of heritability, whereas other classical and mechanical LV function parameters did not. Given the increasing recognition of GLSEndo/GLSEpi ratio as an early/sensitive imaging biomarker of systolic dysfunction, our results suggest the possible existence of individual genetic predispositions to myocardial decline.

Genetic Associations Between Modifiable Risk Factors and Alzheimer Disease

Importance

An estimated 40% of dementia is potentially preventable by modifying 12 risk factors throughout the life course. However, robust evidence for most of these risk factors is lacking. Effective interventions should target risk factors in the causal pathway to dementia.

Objective

To comprehensively disentangle potentially causal aspects of modifiable risk factors for Alzheimer disease (AD) to inspire new drug targeting and improved prevention.

Design, Setting, and Participants

This genetic association study was conducted using 2-sample univariable and multivariable mendelian randomization. Independent genetic variants associated with modifiable risk factors were selected as instrumental variables from genomic consortia. Outcome data for AD were obtained from the European Alzheimer & Dementia Biobank (EADB), generated on August 31, 2021. Main analyses were conducted using the EADB clinically diagnosed end point data. All analyses were performed between April 12 and October 27, 2022.

Exposures

Genetically determined modifiable risk factors.

Main Outcomes and Measures

Odds ratios (ORs) and 95% CIs for AD were calculated per 1-unit change of genetically determined risk factors.

Results

The EADB-diagnosed cohort included 39 106 participants with clinically diagnosed AD and 401 577 control participants without AD. The mean age ranged from 72 to 83 years for participants with AD and 51 to 80 years for control participants. Among participants with AD, 54% to 75% were female, and among control participants, 48% to 60% were female. Genetically determined high-density lipoprotein (HDL) cholesterol concentrations were associated with increased odds of AD (OR per 1-SD increase, 1.10 [95% CI, 1.05-1.16]). Genetically determined high systolic blood pressure was associated with increased risk of AD after adjusting for diastolic blood pressure (OR per 10-mm Hg increase, 1.22 [95% CI, 1.02-1.46]). In a second analysis to minimize bias due to sample overlap, the entire UK Biobank was excluded from the EADB consortium; odds for AD were similar for HDL cholesterol (OR per 1-SD unit increase, 1.08 [95% CI, 1.02-1.15]) and systolic blood pressure after adjusting for diastolic blood pressure (OR per 10-mm Hg increase, 1.23 [95% CI, 1.01-1.50]).

Conclusions and Relevance

This genetic association study found novel genetic associations between high HDL cholesterol concentrations and high systolic blood pressure with higher risk of AD. These findings may inspire new drug targeting and improved prevention implementation.

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.

Risperidone response in patients with schizophrenia drives DNA methylation changes in immune and neuronal systems

Background: The choice of efficient antipsychotic therapy for schizophrenia relies on a time-consuming trial-and-error approach, whereas the social and economic burdens of the disease call for faster alternatives. Material & methods: In a search for predictive biomarkers of antipsychotic response, blood methylomes of 28 patients were analyzed before and 4 weeks into risperidone therapy. Results: Several CpGs exhibiting response-specific temporal dynamics were identified in otherwise temporally stable methylomes and noticeable global response-related differences were observed between good and bad responders. These were associated with genes involved in immunity, neurotransmission and neuronal development. Polymorphisms in many of these genes were previously linked with schizophrenia etiology and antipsychotic response. Conclusion: Antipsychotic response seems to be shaped by both stable and medication-induced methylation differences.

Association and performance of polygenic risk scores for breast cancer among French women presenting or not a familial predisposition to the disease

BACKGROUND

Three partially overlapping breast cancer polygenic risk scores (PRS) comprising 77, 179 and 313 SNPs have been proposed for European-ancestry women by the Breast Cancer Association Consortium (BCAC) for improving risk prediction in the general population. However, the effect of these SNPs may vary from one country to another and within a country because of other factors.

OBJECTIVE

To assess their associated risk and predictive performance in French women from (1) the CECILE population-based case-control study, (2) BRCA1 or BRCA2 (BRCA1/2) pathogenic variant (PV) carriers from the GEMO study, and (3) familial breast cancer cases with no BRCA1/2 PV and unrelated controls from the GENESIS study.

RESULTS

All three PRS were associated with breast cancer in all studies, with odds ratios per standard deviation varying from 1.7 to 2.0 in CECILE and GENESIS, and hazard ratios varying from 1.1 to 1.4 in GEMO. The predictive performance of PRS₃₁₃ in CECILE was similar to that reported in BCAC but lower than that in GENESIS (area under the receiver operating characteristic curve (AUC) = 0.67 and 0.75, respectively). PRS were less performant in BRCA2 and BRCA1 PV carriers (AUC = 0.58 and 0.54 respectively).

CONCLUSION

Our results are in line with previous validation studies in the general population and in BRCA1/2 PV carriers. Additionally, we showed that PRS may be of clinical utility for women with a strong family history of breast cancer and no BRCA1/2 PV, and for those carrying a predicted PV in a moderate-risk gene like ATM, CHEK2 or PALB2.

Weak Association Between Genetic Markers of Hyperuricemia and Cardiorenal Outcomes: Insights From the STANISLAS Study Cohort With a 20-Year Follow-Up

Background

Hyperuricemia is associated with poor cardiovascular outcomes, although it is uncertain whether this relationship is causal in nature. This study aimed to: (1) assess the heritability of serum uric acid (SUA) levels, (2) conduct a genome‐wide association study on SUA levels, and (3) investigate the association between certain single‐nucleotide polymorphisms and target organ damage.

Methods and Results

The STANISLAS (Suivi Temporaire Annuel Non‐Invasif de la Santé des Lorrains Assurés Sociaux) study cohort is a single‐center longitudinal cohort recruited between 1993 and 1995 (visit 1), with a last visit (visit 4 [V4]) performed ≈20 years apart. Serum lipid profile, SUA, urinary albumin/creatinine ratio, estimated glomerular filtration rate, 24‐hour ambulatory blood pressure monitoring, transthoracic echocardiography, pulse wave velocity, and genotyping for each participant were assessed at V4. A total of 1573 participants were included at V4, among whom 1417 had available SUA data at visit 1. Genome‐wide association study results highlighted multiple single‐nucleotide polymorphisms on the SLC2A9 gene linked to SUA levels. Carriers of the most associated mutated SLC2A9 allele (rs16890979) had significantly lower SUA levels. Although SUA level at V4 was highly associated with diabetes, prediabetes, higher body mass index, CRP (C‐reactive protein) levels, estimated glomerular filtration rate variation (visit 1–V4), carotid intima‐media thickness, and pulse wave velocity, rs16890979 was only associated with higher carotid intima‐media thickness.

Conclusions

Our findings demonstrate that rs16890979, a genetic determinant of SUA levels located on the SLC2A9 gene, is associated with carotid intima‐media thickness despite significant associations between SUA levels and several clinical outcomes, thereby lending support to the hypothesis of a link between SUA and cardiovascular disease.

Fatty acid desaturase genetic variations and dietary omega-3 fatty acid intake associate with arterial stiffness

Aims

Long-chain polyunsaturated fatty acids (PUFAs) generate diverse bioactive lipid mediators, which tightly regulate vascular inflammation. The effects of omega-3 PUFA supplementation in cardiovascular prevention however remain controversial. In addition to direct dietary intake, fatty acid desaturases (FADS) determine PUFA levels. Increased arterial stiffness represents an independent predictor of mortality and cardiovascular events. The aim of the present study was to determine the association of PUFA intake, FADS1 genotype, and FADS expression with arterial stiffness.

Methods and results

A cross-sectional population-based cohort study of 1464 participants without overt cardiovascular disease was conducted. Dietary intake was assessed using a food frequency questionnaire. Arterial stiffness was assessed by carotid–femoral pulse wave velocity (cfPWV), and the FADS1 locus variant was determined. Blood cell transcriptomics was performed in a subset of 410 individuals. Pulse wave velocity was significantly associated with the FADS1 locus variant. Differential associations between PWV and omega-3 PUFA intake were observed depending on the FADS1 genotype. High omega-3 PUFA intake attenuated the FADS1 genotype-dependent associations. Carriers of the minor FADS1 locus variant exhibited increased expression of FADS2, which is associated with PWV.

Conclusion

Taken together, these findings point to FADS1 genotype-dependent associations of omega-3 PUFA intake on subclinical cardiovascular disease. These findings may have implications for identifying responders and non-responders to omega-3 PUFA supplementation and open up for personalized dietary counselling in cardiovascular prevention.

Accelerated genome sequencing with controlled costs for infants in intensive care units: a feasibility study in a French hospital network

Obtaining a rapid etiological diagnosis for infants with early-onset rare diseases remains a major challenge. These diseases often have a severe presentation and unknown prognosis, and the genetic causes are very heterogeneous. In a French hospital network, we assessed the feasibility of performing accelerated trio-genome sequencing (GS) with limited additional costs by integrating urgent requests into the routine workflow. In addition to evaluating our capacity for such an approach, this prospective multicentre pilot study was designed to identify pitfalls encountered during its implementation. Over 14 months, we included newborns and infants hospitalized in neonatal or paediatric intensive care units with probable genetic disease and in urgent need for etiological diagnosis to guide medical care. The duration of each step and the pitfalls were recorded. We analysed any deviation from the planned schedule and identified obstacles. Trio-GS was performed for 37 individuals, leading to a molecular diagnosis in 18/37 (49%), and 21/37 (57%) after reanalysis. Corrective measures and protocol adaptations resulted in a median duration of 42 days from blood sampling to report. Accelerated trio-GS is undeniably valuable for individuals in an urgent care context. Such a circuit should coexist with a rapid or ultra-rapid circuit, which, although more expensive, can be used in particularly urgent cases. The drop in GS costs should result in its generalized use for diagnostic purposes and lead to a reduction of the costs of rapid GS.

Impact of natriuretic peptide polymorphisms on diastolic and metabolic function in a populational cohort: insights from the STANISLAS cohort

AIMS

Elevated brain natriuretic peptide (BNP) and the N-terminal fragment of its pro-hormone (NT-proBNP) have become established biomarkers for heart failure and are associated with cardiovascular morbidity and mortality. Investigating sources of inter-individual heterogeneity, particularly genetic factors, could help better identify patients at risk of future cardiovascular disease. The aim of this study was to estimate the heritability of circulating NT-proBNP levels, to perform a genome-wide association study (GWAS) and gene-candidate analysis focused on NPPB-NPPA genes on these levels, and to examine their association with cardiovascular or metabolic outcomes.

METHODS AND RESULTS

A total of 1555 individuals from the STANISLAS study were included. The heritability of circulating NT-proBNP levels was estimated at 15%, with seven single nucleotide polymorphisms (SNPs) reaching the significant threshold in the GWAS. All above SNPs were located on the same gene cluster constituted of MTHFR, CLCN6, NPPA, NPPB, and C1orf167. NPPA gene expression was also associated with NT-proBNP levels. Moreover, six other SNPs from NPPA-NPPB genes were associated with diastolic function (lateral e' on echocardiography) and metabolic features (glycated haemoglobin).

CONCLUSIONS

The heritability of natriuretic peptides appears relatively low (15%) and mainly based on the same gene cluster constituted of MTHFR, CLCN6, NPPA, NPPB, and C1orf167. Natriuretic peptide polymorphisms are associated with natriuretic peptide levels and diastolic function. These results suggest that natriuretic peptide polymorphisms may have an impact in the early stages of cardiovascular and metabolic disease.

Heritability of a resting heart rate in a 20-year follow-up family cohort with GWAS data: Insights from the STANISLAS cohort

BACKGROUND

The association between resting heart rate (HR) and cardiovascular outcomes, especially heart failure, is now well established. However, whether HR is mainly an integrated marker of risk associated with other features, or rather a genetic origin risk marker, is still a matter for debate. Previous studies reported a heritability ranging from 14% to 65%.

DESIGN

We assessed HR heritability in the STANISLAS family-study, based on the data of four visits performed over a 20-year period, and adjusted for most known confounding effects.

METHODS

These analyses were conducted using a linear mixed model, adjusted on age, sex, tea or coffee consumption, beta-blocker use, physical activity, tobacco use, and alcohol consumption to estimate the variance captured by additive genetic effects, via average information restricted maximum likelihood analysis, with both self-reported pedigree and genetic relatedness matrix (GRM) calculated from genome-wide association study data.

RESULTS

Based on the data of all visits, the HR heritability (h2) estimate was 23.2% with GRM and 24.5% with pedigree. However, we found a large heterogeneity of HR heritability estimations when restricting the analysis to each of the four visits (h2 from 19% to 39% using pedigree, and from 14% to 32% using GRM). Moreover, only a little part of variance was explained by the common household effect (<5%), and half of the variance remained unexplained.

CONCLUSION

Using a comprehensive analysis based on a family cohort, including the data of multiple visits and GRM, we found that HR variability is about 25% from genetic origin, 25% from repeated measures and 50% remains unexplained.

Multiethnic genome-wide association study of differentiated thyroid cancer in the EPITHYR consortium

Incidence of differentiated thyroid carcinoma (DTC) varies considerably between ethnic groups, with particularly high incidence rates in Pacific Islanders. DTC is one of the cancers with the highest familial risk suggesting a major role of genetic risk factors, but only few susceptibility loci were identified so far. In order to assess the contribution of known DTC susceptibility loci and to identify new ones, we conducted a multiethnic genome-wide association study (GWAS) in individuals of European ancestry and of Oceanian ancestry from Pacific Islands. Our study included 1554 cases/1973 controls of European ancestry and 301 cases/348 controls of Oceanian ancestry from seven population-based case-control studies participating to the EPITHYR consortium. All participants were genotyped using the OncoArray-500K Beadchip (Illumina). We confirmed the association with the known DTC susceptibility loci at 2q35, 8p12, 9q22.33 and 14q13.3 in the European ancestry population and suggested two novel signals at 1p31.3 and 16q23.2, which were associated with thyroid-stimulating hormone levels in previous GWAS. We additionally replicated an association with 5p15.33 reported previously in Chinese and European populations. Except at 1p31.3, all associations were in the same direction in the population of Oceanian ancestry. We also observed that the frequencies of risk alleles at 2q35, 5p15.33 and 16q23.2 were significantly higher in Oceanians than in Europeans. However, additional GWAS and epidemiological studies in Oceanian populations are needed to fully understand the highest incidence observed in these populations.

Genome-wide association analysis in dilated cardiomyopathy reveals two new players in systolic heart failure on chromosomes 3p25.1 and 22q11.23

AIMS

Our objective was to better understand the genetic bases of dilated cardiomyopathy (DCM), a leading cause of systolic heart failure.

METHODS AND RESULTS

We conducted the largest genome-wide association study performed so far in DCM, with 2719 cases and 4440 controls in the discovery population. We identified and replicated two new DCM-associated loci on chromosome 3p25.1 [lead single-nucleotide polymorphism (SNP) rs62232870, P = 8.7 × 10-11 and 7.7 × 10-4 in the discovery and replication steps, respectively] and chromosome 22q11.23 (lead SNP rs7284877, P = 3.3 × 10-8 and 1.4 × 10-3 in the discovery and replication steps, respectively), while confirming two previously identified DCM loci on chromosomes 10 and 1, BAG3 and HSPB7. A genetic risk score constructed from the number of risk alleles at these four DCM loci revealed a 3-fold increased risk of DCM for individuals with 8 risk alleles compared to individuals with 5 risk alleles (median of the referral population). In silico annotation and functional 4C-sequencing analyses on iPSC-derived cardiomyocytes identify SLC6A6 as the most likely DCM gene at the 3p25.1 locus. This gene encodes a taurine transporter whose involvement in myocardial dysfunction and DCM is supported by numerous observations in humans and animals. At the 22q11.23 locus, in silico and data mining annotations, and to a lesser extent functional analysis, strongly suggest SMARCB1 as the candidate culprit gene.

CONCLUSION

This study provides a better understanding of the genetic architecture of DCM and sheds light on novel biological pathways underlying heart failure.

Gene network and biological pathways associated with susceptibility to differentiated thyroid carcinoma

Variants identified in earlier genome-wide association studies (GWAS) on differentiated thyroid carcinoma (DTC) explain about 10% of the overall estimated genetic contribution and could not provide complete insights into biological mechanisms involved in DTC susceptibility. Integrating systems biology information from model organisms, genome-wide expression data from tumor and matched normal tissue and GWAS data could help identifying DTC-associated genes, and pathways or functional networks in which they are involved. We performed data mining of GWAS data of the EPITHYR consortium (1551 cases and 1957 controls) using various pathways and protein-protein interaction (PPI) annotation databases and gene expression data from The Cancer Genome Atlas. We identified eight DTC-associated genes at known loci 2q35 (DIRC3), 8p12 (NRG1), 9q22 (FOXE1, TRMO, HEMGN, ANP32B, NANS) and 14q13 (MBIP). Using the EW_dmGWAS approach we found that gene networks related to glycogenolysis, glycogen metabolism, insulin metabolism and signal transduction pathways associated with muscle contraction were overrepresented with association signals (false discovery rate adjusted p-value < 0.05). Additionally, suggestive association of 21 KEGG and 75 REACTOME pathways with DTC indicate a link between DTC susceptibility and functions related to metabolism of cholesterol, amino sugar and nucleotide sugar metabolism, steroid biosynthesis, and downregulation of ERBB2 signaling pathways. Together, our results provide novel insights into biological mechanisms contributing to DTC risk.

Fine-mapping of two differentiated thyroid carcinoma susceptibility loci at 2q35 and 8p12 in Europeans, Melanesians and Polynesians

Differentiated thyroid carcinoma (DTC) incidence is characterized by wide ethnic and geographic variations, with high incidence rates observed in Oceanian populations. Genome-wide association studies (GWAS) identified mainly four DTC susceptibility loci at 9q22.33, 14q13.3, 2q35 and 8p12. Here we performed fine-mapping of the 2q35 and 8p12 loci in the population of the EPITHYR consortium that includes Europeans, Melanesians and Polynesians to identify likely causal variants for DTC risk. We conducted a colocalization analysis using eQTLs data to determine the SNPs with the highest probability of causality. At 2q35, we highlighted rs16857609 located in DIRC3. This SNP has a high probability of causality in the three populations, and a significant association in Europeans (OR = 1.4, p = 1.9 x 10-10). It is also associated with expression of DIRC3 and of the nearby gene IGFBP5 in thyroid tumour cells. At 8p12, we identified rs7844425 which was significantly associated with DTC in Europeans (OR = 1.32, p = 7.6 x 10-8) and rs2439304, which was highlighted by the colocalization analysis but only moderately associated with DTC in our dataset (OR = 1.2, p = 0.001). These SNPs are linked to the expression of NRG1 in thyroid tissue. Hence, our study identified novel variants at 2q35 and 8p12 to be prioritized for further functional studies.

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.

Performance comparison of three DNA extraction kits on human whole-exome data from formalin-fixed paraffin-embedded normal and tumor samples

Next-generation sequencing (NGS) studies are becoming routinely used for the detection of novel and clinically actionable DNA variants at a pangenomic scale. Such analyses are now used in the clinical practice to enable precision medicine. Formalin-fixed paraffin-embedded (FFPE) tissues are still one of the most abundant source of cancer clinical specimen, unfortunately this method of preparation is known to degrade DNA and therefore compromise subsequent analysis. Some studies have reported that variant detection can be performed on FFPE samples sequenced with NGS techniques, but few or none have done an in-depth coverage analysis and compared the influence of different state-of-the-art FFPE DNA extraction kits on the quality of the variant calling. Here, we generated 42 human whole-exome sequencing data sets from fresh-frozen (FF) and FFPE samples. These samples include normal and tumor tissues from two different organs (liver and colon), that we extracted with three different FFPE extraction kits (QIAamp DNA FFPE Tissue kit and GeneRead DNA FFPE kit from Qiagen, Maxwell™ RSC DNA FFPE Kit from Promega). We determined the rate of concordance of called variants between matched FF and FFPE samples on all common variants (representing at least 86% of the total number of variants for SNVs). The concordance rate is very high between all matched FF / FFPE pairs, with equivalent values for the three kits we analyzed. On the other hand, when looking at the difference between the total number of variants in FF and FFPE, we find a significant variation for the three different FFPE DNA extraction kits. Coverage analysis shows that FFPE samples have less good indicators than FF samples, yet the coverage quality remains above accepted thresholds. We detect limited but statistically significant variations in coverage indicator values between the three FFPE extraction kits. Globally, the GeneRead and QIAamp kits have better variant calling and coverage indicators than the Maxwell kit on the samples used in this study, although this kit performs better on some indicators and has advantages in terms of practical usage. Taken together, our results confirm the potential of FFPE samples analysis for clinical genomic studies, but also indicate that the choice of a FFPE DNA extraction kit should be done with careful testing and analysis beforehand in order to maximize the accuracy of the results.

New chimeric RNAs in acute myeloid leukemia

Background: High-throughput next generation sequencing (NGS) technologies enable the detection of biomarkers used for tumor classification, disease monitoring and cancer therapy. Whole-transcriptome analysis using RNA-seq is important, not only as a means of understanding the mechanisms responsible for complex diseases but also to efficiently identify novel genes/exons, splice isoforms, RNA editing, allele-specific mutations, differential gene expression and fusion-transcripts or chimeric RNA (chRNA).

Methods: We used Crac, a tool that uses genomic locations and local coverage to classify biological events and directly infer splice and chimeric junctions within a single read. Crac’s algorithm extracts transcriptional chimeric events irrespective of annotation with a high sensitivity, and CracTools was used to aggregate, annotate and filter the chRNA reads. The selected chRNA candidates were validated by real time PCR and sequencing.  In order to check the tumor specific expression of chRNA, we analyzed a publicly available dataset using a new tag search approach.

Results:  We present data related to acute myeloid leukemia (AML) RNA-seq analysis. We highlight novel biological cases of chRNA, in addition to previously well characterized leukemia chRNA. We have identified and validated 17 chRNAs among 3 AML patients: 10 from an AML patient with a translocation between chromosomes 15 and 17 (AML-t(15;17), 4  from patient with normal karyotype (AML-NK) 3 from a patient with chromosomal 16 inversion (AML-inv16). The new fusion transcripts can be classified into four groups according to the exon organization.

Conclusions:  All groups suggest complex but distinct synthesis mechanisms involving either collinear exons of different genes, non-collinear exons, or exons of different chromosomes. Finally, we check tumor-specific expression in a larger RNA-seq AML cohort and identify new AML biomarkers that could improve diagnosis and prognosis of AML.

A major locus on chromosome 3p22 conferring predisposition to human herpesvirus 8 infection

Infection with human herpesvirus 8 (HHV-8), the etiological agent of Kaposi's sarcoma, has been shown to display strong familial aggregation, in countries in which HHV-8 infection is endemic. We investigated 40 large families (608 subjects aged one to 88 years) living in an isolated area of Cameroon in which HHV-8 is highly endemic. We performed a two-step genetic analysis for HHV-8 infection status (HHV-8+/HHV-8- determined by immunofluorescence) consisting of an initial segregation analysis followed by a model-based genome-wide linkage analysis. Overall HHV-8 seroprevalence was 60%, increasing with age. Segregation analysis provided strong evidence for a recessive major gene conferring predisposition to HHV-8 infection. This gene is predicted to have a major effect during childhood, with almost all homozygous predisposed subjects (∼7% of the population) becoming infected by the age of 10. Linkage analysis was carried out on the 15 most informative families, corresponding to 205 genotyped subjects. A single region on chromosome 3p22 was significantly linked to HHV-8 infection (LOD score=3.83, P=2.0 × 10(-5)). This study provides the first evidence that HHV-8 infection in children in endemic areas has a strong genetic basis involving at least one recessive major locus on chromosome 3p22.

Risk HLA-DQA1 and PLA(2)R1 alleles in idiopathic membranous nephropathy

BACKGROUND

Idiopathic membranous nephropathy is a major cause of the nephrotic syndrome in adults, but its etiologic basis is not fully understood. We investigated the genetic basis of biopsy-proven cases of idiopathic membranous nephropathy in a white population.

METHODS

We performed independent genomewide association studies of single-nucleotide polymorphisms (SNPs) in patients with idiopathic membranous nephropathy from three populations of white ancestry (75 French, 146 Dutch, and 335 British patients). The patients were compared with racially matched control subjects; population stratification and quality controls were carried out according to standard criteria. Associations were calculated by means of a chi-square basic allele test; the threshold for significance was adjusted for multiple comparisons (with the Bonferroni method).

RESULTS

In a joint analysis of data from the 556 patients studied (398 men), we identified significant alleles at two genomic loci associated with idiopathic membranous nephropathy. Chromosome 2q24 contains the gene encoding M-type phospholipase A(2) receptor (PLA(2)R1) (SNP rs4664308, P=8.6×10(-29)), previously shown to be the target of an autoimmune response. Chromosome 6p21 contains the gene encoding HLA complex class II HLA-DQ alpha chain 1 (HLA-DQA1) (SNP rs2187668, P=8.0×10(-93)). The association with HLA-DQA1 was significant in all three populations (P=1.8×10(-9), P=5.6×10(-27), and P=5.2×10(-36) in the French, Dutch, and British groups, respectively). The odds ratio for idiopathic membranous nephropathy with homozygosity for both risk alleles was 78.5 (95% confidence interval, 34.6 to 178.2).

CONCLUSIONS

An HLA-DQA1 allele on chromosome 6p21 is most closely associated with idiopathic membranous nephropathy in persons of white ancestry. This allele may facilitate an autoimmune response against targets such as variants of PLA2R1. Our findings suggest a basis for understanding this disease and illuminate how adaptive immunity is regulated by HLA.