DNA-pools targeted-sequencing as a robust cost-effective method to detect rare variants: Application to dilated cardiomyopathy genetic diagnosis.

Dilated cardiomyopathy (DCM) is a heart disease characterized by left ventricular dilatation and systolic dysfunction. In 30% of cases, pathogenic variants, essentially private to each patient, are identified in at least one of almost 50 reported genes. Thus, while costly, exons capture-based Next Generation Sequencing (NGS) of a targeted gene panel appears as the best strategy to genetically diagnose DCM. Here, we report a NGS strategy applied to pools of 8 DNAs from DCM patients and validate its robustness for rare variants detection at 4-fold reduced cost. Our pipeline uses Freebayes to detect variants with the expected 1/16 allele frequency. From the whole set of detected rare variants in 96 pools we set the variants quality parameters optimizing true positives calling. When compared to simplex DNA sequencing in a shared subset of 50 DNAs, 96% of SNVs/InsDel were accurately identified in pools. Extended to the 384 DNAs included in the study, we detected 100 variants (ACMG class 4 and 5), mostly in well-known morbid gene causing DCM such as TTN, MYH7, FLNC, and TNNT2. To conclude, we report an original pool-sequencing NGS method accurately detecting rare variants. This innovative approach is cost-effective for genetic diagnostic in rare diseases.

Changes in circulating miRNA19a-3p precede insulin resistance programmed by intra-uterine growth retardation in mice.

OBJECTIVE: Individuals born with intrauterine growth retardation (IUGR) are more prone to cardio-metabolic diseases as adults, and environmental changes during the perinatal period have been identified as potentially crucial factors. We have studied in a preclinical model early-onset molecular alterations present before the development of a clinical phenotype. METHODS: We used a preclinical mouse model of induced IUGR, in which we modulated the nutrition of the pups during the suckling period, to modify their susceptibility to cardio-metabolic diseases in adulthood. RESULTS: Mice born with IUGR that were overfed (IUGR-O) during lactation rapidly developed obesity, hepatic steatosis and insulin resistance, by three months of age, whereas those subjected to nutrition restriction during lactation (IUGR-R) remained permanently thin and highly sensitive to insulin. Mice born with IUGR and fed normally during lactation (IUGR-N) presented an intermediate phenotype and developed insulin resistance by 12 months of age. Molecular alterations to the insulin signaling pathway with an early onset were observed in the livers of adult IUGR-N mice, nine months before the appearance of insulin resistance. The implication of epigenetic changes was revealed by ChIP sequencing, with both posttranslational H3K4me3 histone modifications and microRNAs involved. CONCLUSIONS: These two changes lead to the coherent regulation of insulin signaling, with a decrease in Akt gene transcription associated with an increase in the translation of its inhibitor, Pten. Moreover, we found that the levels of the implicated miRNA19a-3p also decreased in the blood of young adult IUGR mice nine months before the appearance of insulin resistance, suggesting a possible role for this miRNA as an early circulating biomarker of metabolic fate of potential use for precision medicine.

Anti-integrin αv therapy improves cardiac fibrosis after myocardial infarction by blunting cardiac PW1+ stromal cells.

There is currently no therapy to limit the development of cardiac fibrosis and consequent heart failure. We have recently shown that cardiac fibrosis post-myocardial infarction (MI) can be regulated by resident cardiac cells with a fibrogenic signature and identified by the expression of PW1 (Peg3). Here we identify αV-integrin (CD51) as an essential regulator of cardiac PW1+ cells fibrogenic behavior. We used transcriptomic and proteomic approaches to identify specific cell-surface markers for cardiac PW1+ cells and found that αV-integrin (CD51) was expressed in almost all cardiac PW1+ cells (93% ± 1%), predominantly as the αVß1 complex. αV-integrin is a subunit member of the integrin family of cell adhesion receptors and was found to activate complex of latent transforming growth factor beta (TGFß at the surface of cardiac PW1+ cells. Pharmacological inhibition of αV-integrin reduced the profibrotic action of cardiac PW1+CD51+ cells and was associated with improved cardiac function and animal survival following MI coupled with a reduced infarct size and fibrotic lesion. These data identify a targetable pathway that regulates cardiac fibrosis in response to an ischemic injury and demonstrate that pharmacological inhibition of αV-integrin could reduce pathological outcomes following cardiac ischemia.

Bayesian network analysis of plasma microRNA sequencing data in patients with venous thrombosis.

MicroRNAs (miRNAs) are small regulatory RNAs participating to several biological processes and known to be involved in various pathologies. Measurable in body fluids, miRNAs have been proposed to serve as efficient biomarkers for diseases and/or associated traits. Here, we performed a next-generation-sequencing based profiling of plasma miRNAs in 344 patients with venous thrombosis (VT) and assessed the association of plasma miRNA levels with several haemostatic traits and the risk of VT recurrence. Among the most significant findings, we detected an association between hsa-miR-199b-3p and haematocrit levels (P = 0.0016), these two markers having both been independently reported to associate with VT risk. We also observed suggestive evidence for association of hsa-miR-370-3p (P = 0.019), hsa-miR-27b-3p (P = 0.016) and hsa-miR-222-3p (P = 0.049) with VT recurrence, the observations at the latter two miRNAs confirming the recent findings of Wang et al. Besides, by conducting Genome-Wide Association Studies on miRNA levels and meta-analyzing our results with some publicly available, we identified 21 new associations of single nucleotide polymorphisms with plasma miRNA levels at the statistical significance threshold of P < 5 × 10-8, some of these associations pertaining to thrombosis associated mechanisms. In conclusion, this study provides novel data about the impact of miRNAs' variability in haemostasis and new arguments supporting the association of few miRNAs with the risk of recurrence in patients with venous thrombosis.

Los micro-ARN (miARN) son pequeñas moléculas de ARN reguladoras que participan en varios procesos biológicos y están implicados en diversas patologías. Mensurables en los líquidos corporales, se ha planteado que los miARN pueden ser biomarcadores eficaces para el diagnóstico de enfermedades y/o características asociadas. Aquí hemos llevado a cabo un análisis de miARN plasmático con tecnología de secuenciación de última generación en 344 pacientes con trombosis venosa (TV) y hemos evaluado la asociación de los niveles de miARN con distintas características hemostáticas y el riesgo de recidiva de TV. Entre los hallazgos más significativos, hemos detectado una asociación entre hsa-miR-199b-3p y los niveles de hematocritos (p = 0,0016); dos marcadores que se habían asociado de forma independiente con el riesgo de sufrir TV. Asimismo, hemos observado una evidencia indicativa de asociación entre hsa-miR-370-3p (p = 0,019), hsa-miR-27b-3p (p = 0,016) y hsa-miR-222-3p (p = 0,049) y la recidiva de TV; los resultados los dos últimos miARN confirman los hallazgos recientes de Wang et al. (Clin Epigenetics, 2019). Además, al efectuar estudios de asociación del genoma completo sobre los niveles de miARN y al metaanalizar nuestros resultados con otros disponibles públicamente, hemos identificado 21 asociaciones nuevas de polimorfismos de un solo nucleótido (PSN) con niveles de miARN plasmático con un umbral de significación estadística de p < 5 × 10−8; algunas de estas asociaciones pertenecen a los mecanismos patogénicos de la trombosis.Como conclusión, en este estudio se proporcionan nuevos datos sobre el impacto de la variabilidad de miARN en la hemostasia y nuevos argumentos que apoyan la asociación de algunas secuencias de miARN con el riesgo de recidiva en pacientes con trombosis venosa.

Correction: Exome-wide association study reveals novel susceptibility genes to sporadic dilated cardiomyopathy.

[This corrects the article DOI: 10.1371/journal.pone.0172995.].

Whole blood levels of S1PR4 mRNA associated with cerebral vasospasm after aneurysmal subarachnoid hemorrhage.

OBJECTIVE: The authors sought to identify mRNA biomarkers of cerebral vasospasm in whole blood of patients suffering from aneurysmal subarachnoid hemorrhage (aSAH). METHODS: A prospective transcriptomic study for vasospasm was conducted in whole blood samples of 44 aSAH patients who developed (VSP+ group, n = 22) or did not develop (VSP- group, n = 22) vasospasm. Samples from all patients were profiled for 21,460 mRNA probes using the Illumina Human HT12v4.0 array. Differential statistical analysis was performed using a linear mixed model. RESULTS: Levels of sphingosine-1-phosphate receptor 4 (S1PR4) mRNA were significantly higher (p = 8.03 × 10-6) at presentation in patients who developed vasospasm after aSAH than in patients who did not. CONCLUSIONS: The results, which are consistent with findings of previous experimental investigations conducted in animal models, support the role of S1PR4 and its ligand, sphingosine-1-phosphate (S1P), in arterial-associated vasoconstriction, which suggests that S1PR4 could be used as a biomarker for cerebral vasospasm in aSAH patients.

FLNC pathogenic variants in patients with cardiomyopathies: Prevalence and genotype-phenotype correlations.

Pathogenic variants in FLNC encoding filamin C have been firstly reported to cause myopathies, and were recently linked to isolated cardiac phenotypes. Our aim was to estimate the prevalence of FLNC pathogenic variants in subtypes of cardiomyopathies and to study the relations between phenotype and genotype. DNAs from a cohort of 1150 unrelated index-patients with isolated cardiomyopathy (700 hypertrophic, 300 dilated, 50 restrictive cardiomyopathies, and 100 left ventricle non-compactions) have been sequenced on a custom panel of 51 cardiomyopathy disease-causing genes. An FLNC pathogenic variant was identified in 28 patients corresponding to a prevalence ranging from 1% to 8% depending on the cardiomyopathy subtype. Truncating variants were always identified in patients with dilated cardiomyopathy, while missense or in-frame indel variants were found in other phenotypes. A personal or family history of sudden cardiac death (SCD) was significantly higher in patients with truncating variants than in patients carrying missense variants (P = .01). This work reported the first observation of a left ventricular non-compaction associated with a unique probably causal variant in FLNC which highlights the role of FLNC in cardiomyopathies. A correlation between the nature of the variant and the cardiomyopathy subtype was observed as well as with SCD risk.

OPTIMIR, a novel algorithm for integrating available genome-wide genotype data into miRNA sequence alignment analysis.

Next-generation sequencing is an increasingly popular and efficient approach to characterize the full set of microRNAs (miRNAs) present in human biosamples. MiRNAs' detection and quantification still remain a challenge as they can undergo different posttranscriptional modifications and might harbor genetic variations (polymiRs) that may impact on the alignment step. We present a novel algorithm, OPTIMIR, that incorporates biological knowledge on miRNA editing and genome-wide genotype data available in the processed samples to improve alignment accuracy. OPTIMIR was applied to 391 human plasma samples that had been typed with genome-wide genotyping arrays. OPTIMIR was able to detect genotyping errors, suggested the existence of novel miRNAs and highlighted the allelic imbalance expression of polymiRs in heterozygous carriers. OPTIMIR is written in python, and freely available on the GENMED website (http://www.genmed.fr/index.php/fr/) and on Github (github.com/FlorianThibord/OptimiR).

Whole-Blood miRNA Sequencing Profiling for Vasospasm in Patients With Aneurysmal Subarachnoid Hemorrhage.

Background and Purpose- Arterial vasospasm is a well-known delayed complication of aneurysmal subarachnoid hemorrhage (aSAH). However, no validated biomarker exists to help clinicians discriminating patients with aSAH who will develop vasospasm (VSP+) and identifying those who then deserve aggressive preventive therapy. We hypothesized that whole-blood miRNAs could be a source of candidate biomarkers for vasospasm. Methods- Using a next-generation sequencing approach, we performed whole-blood miRNA profiling between VSP+patients with aSAH and patients who did not develop vasospasm (VSP-) in a prospective cohort of 32 patients. Profiling was performed on the admission day and 3 days before vasospasm. Results- Four hundred forty-two miRNAs were highly expressed in whole blood of patients with aSAH. Among them, hsa-miR-3177-3p demonstrated significant ( P=5.9×10-5; PBonferronicorrected=0.03) lower levels in VSP- compared with VSP+ patients. Looking for whole-blood mRNA correlates of hsa-miR-3177-3p, we observed some evidence that the decrease in hsa-miR-3177-3p levels after aSAH was associated with an increase in LDHA mRNA levels in VSP- ( P<10-3) but not in VSP+ ( P=0.66) patients. Conclusions- Whole-blood miRNA levels of hsa-miR-3177-3p could serve as a biomarker for vasospasm. Clinical Trial Registration- URL: https://www.clinicaltrials.gov . Unique identifier: NCT01779713.

Plasma levels of hsa-miR-152-3p are associated with diabetic nephropathy in patients with type 2 diabetes.

Background: MicroRNAs (miRNAs) are small non-coding RNAs participating in post-transcriptional regulation of genes. Their key role in modulating the susceptibility to human diseases is now widely recognized, in particular in the context of cardiometabolic disorders. The aim of the present study was to identify miRNAs associated with diabetic nephropathy (DN) in patients with type 2 diabetes (T2D). Methods: A next-generation sequencing-based miRNA profiling was performed in a case-control study for DN in plasma samples of 23 T2D patients with DN (cases) and 23 T2D without (controls). The main associations were confirmed using quantitative reverse transcription-polymerase chain reaction and tested for replication in an independent case-control collection of 100 T2D patients, 50 with DN and 50 without. Results: From the 381 known mature miRNAs that were found highly expressed in the discovery samples, we observed and replicated an association between increased plasma levels of hsa-miR-152-3p and DN (P = 4.03 × 10-4 in the combined samples). Hsa-miR-152-3p plasma levels were further found to be positively correlated (P = 0.003) to plasma osmolarity, a surrogate marker for solute carrier net activity, whose regulation is controlled by several genes including SLC5A3, one of the predicted targets of hsa-miR-152-3p. Conclusions: We observed strong evidence for the association of hsa-miR-152-3p plasma levels and DN in patients with T2D, confirming an association previously observed in patients with type 1 diabetes.

Blood triglyceride levels are associated with DNA methylation at the serine metabolism gene PHGDH.

Efficient interventions to reduce blood triglycerides are few; newer and more tolerable intervention targets are needed. Understanding the molecular mechanisms underlying blood triglyceride levels variation is key to identifying new therapies. To explore the role of epigenetic mechanisms on triglyceride levels, a blood methylome scan was conducted in 199 individuals from 5 French-Canadian families ascertained on venous thromboembolism, and findings were replicated in 324 French unrelated patients with venous thromboembolism. Genetic context and functional relevance were investigated. Two DNA methylation sites associated with triglyceride levels were identified. The first one, located in the ABCG1 gene, was recently reported, whereas the second one, located in the promoter of the PHGDH gene, is novel. The PHGDH methylation site, cg14476101, was found to be associated with variation in triglyceride levels in a threshold manner: cg14476101 was inversely associated with triglyceride levels only when triglyceride levels were above 1.12 mmol/L (discovery P-value = 8.4 × 10-6; replication P-value = 0.0091). Public databases findings supported a functional role of cg14476101 on PHGDH expression. PHGDH catalyses the first step in the serine biosynthesis pathway. These findings highlight the role of epigenetic regulation of the PHGDH gene in triglyceride metabolism, providing novel insights on putative intervention targets.

Fibrogenic Potential of PW1/Peg3 Expressing Cardiac Stem Cells.

BACKGROUND: Pw1 gene expression is a marker of adult stem cells in a wide range of tissues. PW1-expressing cells are detected in the heart but are not well characterized. OBJECTIVES: The authors characterized cardiac PW1-expressing cells and their cell fate potentials in normal hearts and during cardiac remodeling following myocardial infarction (MI). METHODS: A human cardiac sample was obtained from a patient presenting with reduced left ventricular (LV) function following a recent MI. The authors used the PW1nLacZ+/- reporter mouse to identify, track, isolate, and characterize PW1-expressing cells in the LV myocardium in normal and ischemic conditions 7 days after complete ligature of the left anterior descending coronary artery. RESULTS: In both human and mouse ischemic hearts, PW1 expression was found in cells that were mainly located in the infarct and border zones. Isolated cardiac resident PW1+ cells form colonies and have the potential to differentiate into multiple cardiac and mesenchymal lineages, with preferential differentiation into fibroblast-like cells but not into cardiomyocytes. Lineage-tracing experiments revealed that PW1+ cells differentiated into fibroblasts post-MI. Although the expression of c-Kit and PW1 showed little overlap in normal hearts, a marked increase in cells coexpressing both markers was observed in ischemic hearts (0.1 ± 0.0% in control vs. 5.7 ± 1.2% in MI; p < 0.001). In contrast to the small proportion of c-Kit+/PW1- cells that showed cardiogenic potential, c-Kit+/PW1+ cells were fibrogenic. CONCLUSIONS: This study demonstrated the existence of a novel population of resident adult cardiac stem cells expressing PW1+ and their involvement in fibrotic remodeling after MI.

Exome-wide association study reveals novel susceptibility genes to sporadic dilated cardiomyopathy.

AIMS: Dilated cardiomyopathy (DCM) is an important cause of heart failure with a strong familial component. We performed an exome-wide array-based association study (EWAS) to assess the contribution of missense variants to sporadic DCM. METHODS AND RESULTS: 116,855 single nucleotide variants (SNVs) were analyzed in 2796 DCM patients and 6877 control subjects from 6 populations of European ancestry. We confirmed two previously identified associations with SNVs in BAG3 and ZBTB17 and discovered six novel DCM-associated loci (Q-value<0.01). The lead-SNVs at novel loci are common and located in TTN, SLC39A8, MLIP, FLNC, ALPK3 and FHOD3. In silico fine mapping identified HSPB7 as the most likely candidate at the ZBTB17 locus. Rare variant analysis (MAF<0.01) demonstrated significant association for TTN variants only (P = 0.0085). All candidate genes but one (SLC39A8) exhibit preferential expression in striated muscle tissues and mutations in TTN, BAG3, FLNC and FHOD3 are known to cause familial cardiomyopathy. We also investigated a panel of 48 known cardiomyopathy genes. Collectively, rare (n = 228, P = 0.0033) or common (n = 36, P = 0.019) variants with elevated in silico severity scores were associated with DCM, indicating that the spectrum of genes contributing to sporadic DCM extends beyond those identified here. CONCLUSION: We identified eight loci independently associated with sporadic DCM. The functions of the best candidate genes at these loci suggest that proteostasis regulation might play a role in DCM pathophysiology.

Bariatric Surgery Induces Disruption in Inflammatory Signaling Pathways Mediated by Immune Cells in Adipose Tissue: A RNA-Seq Study.

BACKGROUND: Bariatric surgery is associated to improvements in obesity-associated comorbidities thought to be mediated by a decrease of adipose inflammation. However, the molecular mechanisms behind these beneficial effects are poorly understood. METHODOLOGY/PRINCIPAL FINDINGS: We analyzed RNA-seq expression profiles in adipose tissue from 22 obese women before and 3 months after surgery. Of 15,972 detected genes, 1214 were differentially expressed after surgery at a 5% false discovery rate. Upregulated genes were mostly involved in the basal cellular machinery. Downregulated genes were enriched in metabolic functions of adipose tissue. At baseline, 26 modules of coexpressed genes were identified. The four most stable modules reflected the innate and adaptive immune responses of adipose tissue. A first module reflecting a non-specific signature of innate immune cells, mainly macrophages, was highly conserved after surgery with the exception of DUSP2 and CD300C. A second module reflected the adaptive immune response elicited by T lymphocytes; after surgery, a disconnection was observed between genes involved in T-cell signaling and mediators of the signal transduction such as CXCR1, CXCR2, GPR97, CCR7 and IL7R. A third module reflected neutrophil-mediated inflammation; after surgery, several genes were dissociated from the module, including S100A8, S100A12, CD300E, VNN2, TUBB1 and FAM65B. We also identified a dense network of 19 genes involved in the interferon-signaling pathway which was strongly preserved after surgery, with the exception of DDX60, an antiviral factor involved in RIG-I-mediated interferon signaling. A similar loss of connection was observed in lean mice compared to their obese counterparts. CONCLUSIONS/SIGNIFICANCE: These results suggest that improvements of the inflammatory state following surgery might be explained by a disruption of immuno-inflammatory cascades involving a few crucial molecules which could serve as potential therapeutic targets.

PLA2G10 Gene Variants, sPLA2 Activity, and Coronary Heart Disease Risk.

BACKGROUND: Observational studies report that secretory phospholipase A2 (sPLA2) activity is a marker for coronary heart disease (CHD) risk, and activity measures are thought to represent the composite activity of sPLA2-IIA, -V, and -X. The aim of this study was to use genetic variants of PLA2G10, encoding sPLA2-X, to investigate the contribution of sPLA2-X to the measure of sPLA2 activity and coronary heart disease (CHD) risk traits and outcome. METHODS AND RESULTS: Three PLA2G10 tagging single-nucleotide polymorphisms (rs72546339, rs72546340, and rs4003232) and a previously studied PLA2G10 coding single-nucleotide polymorphism rs4003228, R38C, were genotyped in a nested case: control cohort drawn from the prospective EPIC-Norfolk Study (2175 cases and 2175 controls). Meta-analysis of rs4003228 (R38C) and CHD was performed using data from the Northwick Park Heart Study II and 2 published cohorts AtheroGene and SIPLAC, providing in total an additional 1884 cases and 3119 controls. EPIC-Norfolk subjects in the highest tertile of sPLA2 activity were older and had higher inflammatory markers compared with those in the lowest tertile for sPLA2 activity. None of the PLA2G10 tagging single-nucleotide polymorphism nor R38C, a functional variant, were significantly associated with sPLA2 activity, intermediate CHD risk traits, or CHD risk. In meta-analysis, the summary odds ratio for R38C was odds ratio=0.97 (95% confidence interval, 0.77-1.22). CONCLUSIONS: PLA2G10 variants are not significantly associated with plasma sPLA2 activity or with CHD risk.

Human CalDAG-GEFI gene (RASGRP2) mutation affects platelet function and causes severe bleeding.

The nature of an inherited platelet disorder was investigated in three siblings affected by severe bleeding. Using whole-exome sequencing, we identified the culprit mutation (cG742T) in the RAS guanyl-releasing protein-2 (RASGRP2) gene coding for calcium- and DAG-regulated guanine exchange factor-1 (CalDAG-GEFI). Platelets from individuals carrying the mutation present a reduced ability to activate Rap1 and to perform proper αIIbß3 integrin inside-out signaling. Expression of CalDAG-GEFI mutant in HEK293T cells abolished Rap1 activation upon stimulation. Nevertheless, the PKC- and ADP-dependent pathways allow residual platelet activation in the absence of functional CalDAG-GEFI. The mutation impairs the platelet's ability to form thrombi under flow and spread normally as a consequence of reduced Rac1 GTP-binding. Functional deficiencies were confined to platelets and megakaryocytes with no leukocyte alteration. This contrasts with the phenotype seen in type III leukocyte adhesion deficiency caused by the absence of kindlin-3. Heterozygous did not suffer from bleeding and have normal platelet aggregation; however, their platelets mimicked homozygous ones by failing to undergo normal adhesion under flow and spreading. Rescue experiments on cultured patient megakaryocytes corrected the functional deficiency after transfection with wild-type RASGRP2. Remarkably, the presence of a single normal allele is sufficient to prevent bleeding, making CalDAG-GEFI a novel and potentially safe therapeutic target to prevent thrombosis.

EIF2AK4 mutations cause pulmonary veno-occlusive disease, a recessive form of pulmonary hypertension.

Pulmonary veno-occlusive disease (PVOD) is a rare and devastating cause of pulmonary hypertension that is characterized histologically by widespread fibrous intimal proliferation of septal veins and preseptal venules and is frequently associated with pulmonary capillary dilatation and proliferation. PVOD is categorized into a separate pulmonary arterial hypertension-related group in the current classification of pulmonary hypertension. PVOD presents either sporadically or as familial cases with a seemingly recessive mode of transmission. Using whole-exome sequencing, we detected recessive mutations in EIF2AK4 (also called GCN2) that cosegregated with PVOD in all 13 families studied. We also found biallelic EIF2AK4 mutations in 5 of 20 histologically confirmed sporadic cases of PVOD. All mutations, either in a homozygous or compound-heterozygous state, disrupted the function of the gene. These findings point to EIF2AK4 as the major gene that is linked to PVOD development and contribute toward an understanding of the complex genetic architecture of pulmonary hypertension.

Whole-Genome Sequence of Mycobacterium abscessus Clinical Strain V06705.

Infection caused by Mycobacterium abscessus strains is a growing cause of concern in both community-acquired and health care-associated diseases, as these organisms naturally display multiple drug resistances. We report an annotated draft genome sequence of M. abscessus strain V06705 obtained from a patient in France.

Integrating genome-wide genetic variations and monocyte expression data reveals trans-regulated gene modules in humans.

One major expectation from the transcriptome in humans is to characterize the biological basis of associations identified by genome-wide association studies. So far, few cis expression quantitative trait loci (eQTLs) have been reliably related to disease susceptibility. Trans-regulating mechanisms may play a more prominent role in disease susceptibility. We analyzed 12,808 genes detected in at least 5% of circulating monocyte samples from a population-based sample of 1,490 European unrelated subjects. We applied a method of extraction of expression patterns-independent component analysis-to identify sets of co-regulated genes. These patterns were then related to 675,350 SNPs to identify major trans-acting regulators. We detected three genomic regions significantly associated with co-regulated gene modules. Association of these loci with multiple expression traits was replicated in Cardiogenics, an independent study in which expression profiles of monocytes were available in 758 subjects. The locus 12q13 (lead SNP rs11171739), previously identified as a type 1 diabetes locus, was associated with a pattern including two cis eQTLs, RPS26 and SUOX, and 5 trans eQTLs, one of which (MADCAM1) is a potential candidate for mediating T1D susceptibility. The locus 12q24 (lead SNP rs653178), which has demonstrated extensive disease pleiotropy, including type 1 diabetes, hypertension, and celiac disease, was associated to a pattern strongly correlating to blood pressure level. The strongest trans eQTL in this pattern was CRIP1, a known marker of cellular proliferation in cancer. The locus 12q15 (lead SNP rs11177644) was associated with a pattern driven by two cis eQTLs, LYZ and YEATS4, and including 34 trans eQTLs, several of them tumor-related genes. This study shows that a method exploiting the structure of co-expressions among genes can help identify genomic regions involved in trans regulation of sets of genes and can provide clues for understanding the mechanisms linking genome-wide association loci to disease.

Genetics of venous thrombosis: insights from a new genome wide association study.

BACKGROUND: Venous Thrombosis (VT) is a common multifactorial disease associated with a major public health burden. Genetics factors are known to contribute to the susceptibility of the disease but how many genes are involved and their contribution to VT risk still remain obscure. We aimed to identify genetic variants associated with VT risk. METHODOLOGY/PRINCIPAL FINDINGS: We conducted a genome-wide association study (GWAS) based on 551,141 SNPs genotyped in 1,542 cases and 1,110 controls. Twelve SNPs reached the genome-wide significance level of 2.0×10(-8) and encompassed four known VT-associated loci, ABO, F5, F11 and FGG. By means of haplotype analyses, we also provided novel arguments in favor of a role of HIVEP1, PROCR and STAB2, three loci recently hypothesized to participate in the susceptibility to VT. However, no novel VT-associated loci came out of our GWAS. Using a recently proposed statistical methodology, we also showed that common variants could explain about 35% of the genetic variance underlying VT susceptibility among which 3% could be attributable to the main identified VT loci. This analysis additionally suggested that the common variants left to be identified are not uniformly distributed across the genome and that chromosome 20, itself, could contribute to ∼7% of the total genetic variance. CONCLUSIONS/SIGNIFICANCE: This study might also provide a valuable source of information to expand our understanding of biological mechanisms regulating quantitative biomarkers for VT.