Evaluation of a whole-exome sequencing pipeline and benchmarking of causal germline variant prioritizers.

Most causal variants of Mendelian diseases are exonic. Whole-exome sequencing (WES) has become the diagnostic gold standard, but causative variant prioritization constitutes a bottleneck. Here we assessed an in-house sample-to-sequence pipeline and benchmarked free prioritization tools for germline causal variants from WES data. WES of 61 unselected patients with a known genetic disease cause was obtained. Variant prioritizations were performed by diverse tools and recorded to obtain a diagnostic yield when the causal variant was present in the first, fifth, and 10th top rankings. A fraction of causal variants was not captured by WES (8.2%) or did not pass the quality control criteria (13.1%). Most of the applications inspected were unavailable or had technical limitations, leaving nine tools for complete evaluation. Exomiser performed best in the top first rankings, while LIRICAL led in the top fifth rankings. Based on the more conservative top 10th rankings, Xrare had the highest diagnostic yield, followed by a three-way tie among Exomiser, LIRICAL, and PhenIX, then followed by AMELIE, TAPES, Phen-Gen,  AIVar, and VarNote-PAT. Xrare, Exomiser, LIRICAL, and PhenIX are the most efficient options for variant prioritization in real patient WES data.

A genome-wide association study of survival in patients with sepsis.

BACKGROUND: Sepsis is a severe systemic inflammatory response to infections that is accompanied by organ dysfunction and has a high mortality rate in adult intensive care units. Most genetic studies have identified gene variants associated with development and outcomes of sepsis focusing on biological candidates. We conducted the first genome-wide association study (GWAS) of 28-day survival in adult patients with sepsis. METHODS: This study was conducted in two stages. The first stage was performed on 687 European sepsis patients from the GEN-SEP network and 7.5 million imputed variants. Association testing was conducted with Cox regression models, adjusting by sex, age, and the main principal components of genetic variation. A second stage focusing on the prioritized genetic variants was performed on 2,063 ICU sepsis patients (1362 European Americans and 701 African-Americans) from the MESSI study. A meta-analysis of results from the two stages was conducted and significance was established at p < 5.0 × 10-8. Whole-blood transcriptomic, functional annotations, and sensitivity analyses were evaluated on the identified genes and variants. FINDINGS: We identified three independent low-frequency variants associated with reduced 28-day sepsis survival, including a missense variant in SAMD9 (hazard ratio [95% confidence interval] = 1.64 [1.37-6.78], p = 4.92 × 10-8). SAMD9 encodes a possible mediator of the inflammatory response to tissue injury. INTERPRETATION: We performed the first GWAS of 28-day sepsis survival and identified novel variants associated with reduced survival. Larger sample size studies are needed to better assess the genetic effects in sepsis survival and to validate the findings.

Genome-wide association study of asthma exacerbations despite inhaled corticosteroid use.

RATIONALE: Substantial variability in response to asthma treatment with inhaled corticosteroids (ICS) has been described among individuals and populations, suggesting the contribution of genetic factors. Nonetheless, only a few genes have been identified to date. We aimed to identify genetic variants associated with asthma exacerbations despite ICS use in European children and young adults and to validate the findings in non-Europeans. Moreover, we explored whether a gene-set enrichment analysis could suggest potential novel asthma therapies. METHODS: A genome-wide association study (GWAS) of asthma exacerbations was tested in 2681 children of European descent treated with ICS from eight studies. Suggestive association signals were followed up for replication in 538 European asthma patients. Further evaluation was performed in 1773 non-Europeans. Variants revealed by published GWAS were assessed for replication. Additionally, gene-set enrichment analysis focused on drugs was performed. RESULTS: 10 independent variants were associated with asthma exacerbations despite ICS treatment in the discovery phase (p≤5×10-6). Of those, one variant at the CACNA2D3-WNT5A locus was nominally replicated in Europeans (rs67026078; p=0.010), but this was not validated in non-European populations. Five other genes associated with ICS response in previous studies were replicated. Additionally, an enrichment of associations in genes regulated by trichostatin A treatment was found. CONCLUSIONS: The intergenic region of CACNA2D3 and WNT5A was revealed as a novel locus for asthma exacerbations despite ICS treatment in European populations. Genes associated were related to trichostatin A, suggesting that this drug could regulate the molecular mechanisms involved in treatment response.

Admixture mapping of asthma in southwestern Europeans with North African ancestry influences.

The prevalence of asthma symptoms in Canary Islanders, a southwestern European population from Spain, is almost three times higher than the country average. Because the genetic risks identified so far explain <5% of asthma heritability, here we aimed to discover new asthma loci by completing the first admixture mapping study in Canary Islanders leveraging their distinctive genetic makeup, where significant northwest African influences coexist in the European genetic diversity landscape. A 2-stage study was conducted in 1,491 unrelated individuals self-declaring having a Canary Islands origin for the 4 grandparents. Local ancestry estimates were obtained for the shared positions with reference data from putative ancestral populations from Europe, North Africa, and sub-Saharan Africa. Case-control deviations in local ancestry were tested for each ancestry separately using logistic regressions adjusted for principal components, followed by fine-mapping analyses based on imputed genotypes and analyses of the likely deleterious exonic variants. The admixture mapping analysis of asthma detected that local North African ancestry in a locus spanning 365 kb of chromosome 16q23.3 was associated with asthma risk at study-wide significance [lowest P = 1.12 × 10-4; odds ratio (OR) = 2.05; 95% confidence interval (CI) = 1.41-3.00]. Fine-mapping studies identified a variant associated with asthma, and results were replicated in independent samples (rs3852738, OR = 1.34; 95% CI = 1.13-1.59, P = 7.58 × 10-4). Whole exome sequencing data from a subset of individuals revealed an enrichment of likely deleterious variants among asthma cases in 16q23.3, particularly in the phospholipase Cγ2 (PLCG2) gene (P = 3.67 × 10-4). By completing the first mapping study of asthma in admixed populations from Europe, our results revealed a new plausible asthma locus.

Interactive Web-Based Resource for Annotation of Genetic Variants Causing Hereditary Angioedema (HADA): Database Development, Implementation, and Validation.

BACKGROUND: Hereditary angioedema is a rare genetic condition caused by C1 esterase inhibitor deficiency, dysfunction, or kinin cascade dysregulation, leading to an increased bradykinin plasma concentration. Hereditary angioedema is a poorly recognized clinical entity and is very often misdiagnosed as a histaminergic angioedema. Despite its genetic nature, first-line genetic screening is not integrated in routine diagnosis. Consequently, a delay in the diagnosis, and inaccurate or incomplete diagnosis and treatment of hereditary angioedema are common. OBJECTIVE: In agreement with recent recommendations from the International Consensus on the Use of Genetics in the Management of Hereditary Angioedema, to facilitate the clinical diagnosis and adapt it to the paradigm of precision medicine and next-generation sequencing-based genetic tests, we aimed to develop a genetic annotation tool, termed Hereditary Angioedema Database Annotation (HADA). METHODS: HADA is built on top of a database of known variants affecting function, including precomputed pathogenic assessment of each variant and a ranked classification according to the current guidelines from the American College of Medical Genetics and Genomics. RESULTS: HADA is provided as a freely accessible, user-friendly web-based interface with versatility for the entry of genetic information. The underlying database can also be incorporated into automated command-line stand-alone annotation tools. CONCLUSIONS: HADA can achieve the rapid detection of variants affecting function for different hereditary angioedema types, and further integrates useful information to reduce the diagnosis odyssey and improve its delay.

Genomic Analyses of Human European Diversity at the Southwestern Edge: Isolation, African Influence and Disease Associations in the Canary Islands.

Despite the genetic resemblance of Canary Islanders to other southern European populations, their geographical isolation and the historical admixture of aborigines (from North Africa) with sub-Saharan Africans and Europeans have shaped a distinctive genetic makeup that likely affects disease susceptibility and health disparities. Based on single nucleotide polymorphism array data and whole genome sequencing (30×), we inferred that the last African admixture took place ∼14 generations ago and estimated that up to 34% of the Canary Islander genome is of recent African descent. The length of regions in homozygosis and the ancestry-related mosaic organization of the Canary Islander genome support the view that isolation has been strongest on the two smallest islands. Furthermore, several genomic regions showed significant and large deviations in African or European ancestry and were significantly enriched in genes involved in prevalent diseases in this community, such as diabetes, asthma, and allergy. The most prominent of these regions were located near LCT and the HLA, two well-known targets of selection, at which 40‒50% of the Canarian genome is of recent African descent according to our estimates. Putative selective signals were also identified in these regions near the SLC6A11-SLC6A1, KCNMB2, and PCDH20-PCDH9 genes. Taken together, our findings provide solid evidence of a significant recent African admixture, population isolation, and adaptation in this part of Europe, with the favoring of African alleles in some chromosome regions. These findings may have medical implications for populations of recent African ancestry.

Common variants of NFE2L2 gene predisposes to acute respiratory distress syndrome in patients with severe sepsis.

INTRODUCTION: The purpose of this study was to investigate whether common variants across the nuclear factor erythroid 2-like 2 (NFE2L2) gene contribute to the development of the acute respiratory distress syndrome (ARDS) in patients with severe sepsis. NFE2L2 is involved in the response to oxidative stress, and it has been shown to be associated with the development of ARDS in trauma patients. METHODS: We performed a case-control study of 321 patients fulfilling international criteria for severe sepsis and ARDS who were admitted to a Spanish network of post-surgical and critical care units, as well as 871 population-based controls. Six tagging single-nucleotide polymorphisms (SNPs) of NFE2L2 were genotyped, and, after further imputation of additional 34 SNPs, association testing with ARDS susceptibility was conducted using logistic regression analysis. RESULTS: After multiple testing adjustments, our analysis revealed 10 non-coding SNPs in tight linkage disequilibrium (0.75 ≤ r (2) ≤ 1) that were associated with ARDS susceptibility as a single association signal. One of those SNPs (rs672961) was previously associated with trauma-induced ARDS and modified the promoter activity of the NFE2L2 gene, showing an odds ratio of 1.93 per T allele (95 % confidence interval, 1.17-3.18; p = 0.0089). CONCLUSIONS: Our findings support the involvement of NFE2L2 gene variants in ARDS susceptibility and reinforce further exploration of the role of oxidant stress response as a risk factor for ARDS in critically ill patients.

No association between genetic ancestry and susceptibility to asthma or atopy in Canary Islanders.

Asthma is a complex respiratory disease characterized by chronic inflammation of airways and frequently associated with atopic symptoms. The population from the Canary Islands, which has resulted from a recent admixture of North African and Iberian populations, shows the highest prevalence of asthma and atopic symptoms among the Spanish populations. Although environmental particularities would account for the majority of such disparity, genetic ancestry might play a role in increasing the susceptibility of asthma or atopy, as have been demonstrated in other recently African-admixed populations. Here, we aimed to explore whether genetic ancestry was associated with asthma or related traits in the Canary Islanders. For that, a total of 734 DNA samples from unrelated individuals of the GOA study, self-reporting at least two generations of ancestors from the Canary Islands (391 asthmatics and 343 controls), were successfully genotyped for 83 ancestry informative markers (AIMs), which allowed to precisely distinguishing between North African and Iberian ancestries. No association was found between genetic ancestry and asthma or related traits after adjusting by demographic variables differing among compared groups. Similarly, none of the individual AIMs was associated with asthma when results were considered in the context of the multiple comparisons performed (0.005 ≤ p value ≤ 0.042; 0.221 ≤ q value ≤ 0.443). Our results suggest that if genetic ancestry were involved in the susceptibility to asthma or related traits among Canary Islanders, its effects would be modest. Larger studies, examining more genetic variants, would be needed to explore such possibility.

Admixture Mapping of Sepsis in European Individuals With African Ancestries.

Sepsis is a severe systemic inflammatory response to infections that is accompanied by organ dysfunction. Although the ancestral genetic background is a relevant factor for sepsis susceptibility, there is a lack of studies using the genetic singularities of a recently admixed population to identify loci involved in sepsis susceptibility. Here we aimed to discover new sepsis loci by completing the first admixture mapping study of sepsis in Canary Islanders, leveraging their distinctive genetic makeup as a mixture of Europeans and African ancestries. We used a case-control approach and inferred local ancestry blocks from genome-wide data from 113,414 polymorphisms genotyped in 343 patients with sepsis and 410 unrelated controls, all ascertained for grandparental origin in the Canary Islands (Spain). Deviations in local ancestries between cases and controls were tested using logistic regressions, followed by fine-mapping analyses based on imputed genotypes, in silico functional assessments, and gene expression analysis centered on the region of interest. The admixture mapping analysis detected that local European ancestry in a locus spanning 1.2 megabases of chromosome 8p23.1 was associated with sepsis (lowest p = 1.37 × 10-4; Odds Ratio [OR] = 0.51; 95%CI = 0.40-0.66). Fine-mapping studies prioritized the variant rs13249564 within intron 1 of MFHAS1 gene associated with sepsis (p = 9.94 × 10-4; OR = 0.65; 95%CI = 0.50-0.84). Functional and gene expression analyses focused on 8p23.1 allowed us to identify alternative genes with possible biological plausibility such as defensins, which are well-known effector molecules of innate immunity. By completing the first admixture mapping study of sepsis, our results revealed a new genetic locus (8p23.1) harboring a number of genes with plausible implications in sepsis susceptibility.

Developing CIRdb as a catalog of natural genetic variation in the Canary Islanders.

The current inhabitants of the Canary Islands have a unique genetic makeup in the European diversity landscape due to the existence of African footprints from recent admixture events, especially of North African components (> 20%). The underrepresentation of non-Europeans in genetic studies and the sizable North African ancestry, which is nearly absent from all existing catalogs of worldwide genetic diversity, justify the need to develop CIRdb, a population-specific reference catalog of natural genetic variation in the Canary Islanders. Based on array genotyping of the selected unrelated donors and comparisons against available datasets from European, sub-Saharan, and North African populations, we illustrate the intermediate genetic differentiation of Canary Islanders between Europeans and North Africans and the existence of within-population differences that are likely driven by genetic isolation. Here we describe the overall design and the methods that are being implemented to further develop CIRdb. This resource will help to strengthen the implementation of Precision Medicine in this population by contributing to increase the diversity in genetic studies. Among others, this will translate into improved ability to fine map disease genes and simplify the identification of causal variants and estimate the prevalence of unattended Mendelian diseases.

A catalog of the genetic causes of hereditary angioedema in the Canary Islands (Spain).

Hereditary angioedema (HAE) is a rare disease where known causes involve C1 inhibitor dysfunction or dysregulation of the kinin cascade. The updated HAE management guidelines recommend performing genetic tests to reach a precise diagnosis. Unfortunately, genetic tests are still uncommon in the diagnosis routine. Here, we characterized for the first time the genetic causes of HAE in affected families from the Canary Islands (Spain). Whole-exome sequencing data was obtained from 41 affected patients and unaffected relatives from 29 unrelated families identified in the archipelago. The Hereditary Angioedema Database Annotation (HADA) tool was used for pathogenicity classification and causal variant prioritization among the genes known to cause HAE. Manual reclassification of prioritized variants was used in those families lacking known causal variants. We detected a total of eight different variants causing HAE in this patient series, affecting essentially SERPING1 and F12 genes, one of them being a novel SERPING1 variant (c.686-12A>G) with a predicted splicing effect which was reclassified as likely pathogenic in one family. Altogether, the diagnostic yield by assessing previously reported causal genes and considering variant reclassifications according to the American College of Medical Genetics guidelines reached 66.7% (95% Confidence Interval [CI]: 30.1-91.0) in families with more than one affected member and 10.0% (95% CI: 1.8-33.1) among cases without family information for the disease. Despite the genetic causes of many patients remain to be identified, our results reinforce the need of genetic tests as first-tier diagnostic tool in this disease, as recommended by the international WAO/EAACI guidelines for the management of HAE.

Interleukin-1 receptor-associated kinase 3 gene associates with susceptibility to acute lung injury.

Sepsis is the most common cause of acute lung injury (ALI), leading to organ dysfunction and death in critically ill patients. Previous studies associated variants of interleukin-1 receptor-associated kinase genes (IRAKs) with differential immune responses to pathogens and with outcomes during sepsis, and revealed that increased expression levels of the IRAK3 gene were correlated with poor outcomes during sepsis. Here we explored whether common variants of the IRAK3 gene were associated with susceptibility to, and outcomes of, severe sepsis. After our discovery of polymorphism, we genotyped a subset of seven single-nucleotide polymorphisms (SNPs) in 336 population-based control subjects and 214 patients with severe sepsis, collected as part of a prospective study of adults from a Spanish network of intensive care units. Whereas IRAK3 SNPs were not associated with susceptibility to severe sepsis, rs10506481 showed a significant association with the development of ALI among patients with sepsis (P = 0.007). The association remained significant after adjusting for multiple comparisons, population stratification, and clinical variables (odds ratio, 2.50; 95% confidence interval, 1.15-5.47; P = 0.021). By imputation, we revealed three additional SNPs independently associated with ALI (P < 0.01). One of these (rs1732887) predicted the disruption of a putative human-mouse conserved transcription factor binding site, and demonstrated functional effects in vitro (P = 0.017). Despite the need for replication in independent studies, our data suggest that common SNPs in the IRAK3 gene may be determinants of sepsis-induced ALI.

Targeted analysis of genomic regions enriched in African ancestry reveals novel classical HLA alleles associated with asthma in Southwestern Europeans.

Despite asthma has a considerable genetic component, an important proportion of genetic risks remain unknown, especially for non-European populations. Canary Islanders have the largest African genetic ancestry observed among Southwestern Europeans and the highest asthma prevalence in Spain. Here we examined broad chromosomal regions previously associated with an excess of African genetic ancestry in Canary Islanders, with the aim of identifying novel risk variants associated with asthma susceptibility. In a two-stage cases-control study, we revealed a variant within HLA-DQB1 significantly associated with asthma risk (rs1049213, meta-analysis p = 1.30 × 10-7, OR [95% CI] = 1.74 [1.41-2.13]) previously associated with asthma and broad allergic phenotype. Subsequent fine-mapping analyses of classical HLA alleles revealed a novel allele significantly associated with asthma protection (HLA-DQA1*01:02, meta-analysis p = 3.98 × 10-4, OR [95% CI] = 0.64 [0.50-0.82]) that had been linked to infectious and autoimmune diseases, and peanut allergy. HLA haplotype analyses revealed a novel haplotype DQA1*01:02-DQB1*06:04 conferring asthma protection (meta-analysis p = 4.71 × 10-4, OR [95% CI] = 0.47 [0.29- 0.73]).

Whole-Blood Mitochondrial DNA Copies Are Associated With the Prognosis of Acute Respiratory Distress Syndrome After Sepsis.

Acute respiratory distress syndrome (ARDS) is an inflammatory process of the lungs that develops primarily in response to pulmonary or systemic sepsis, resulting in a disproportionate death toll in intensive care units (ICUs). Given its role as a critical activator of the inflammatory and innate immune responses, previous studies have reported that an increase of circulating cell-free mitochondrial DNA (mtDNA) is a biomarker for fatal outcome in the ICU. Here we analyzed the association of whole-blood mtDNA (wb-mtDNA) copies with 28-day survival from sepsis and sepsis-associated ARDS. We analyzed mtDNA data from 687 peripheral whole-blood samples within 24 h of sepsis diagnosis from unrelated Spanish patients with sepsis (264 with ARDS) included in the GEN-SEP study. The wb-mtDNA copies were obtained from the array intensities of selected probes, with 100% identity with mtDNA and with the largest number of mismatches with the nuclear sequences, and normalized across the individual-probe intensities. We used Cox regression models for testing the association with 28-day survival. We observed that wb-mtDNA copies were significantly associated with 28-day survival in ARDS patients (hazard ratio = 3.65, 95% confidence interval = 1.39-9.59, p = 0.009) but not in non-ARDS patients. Our findings support that wb-mtDNA copies at sepsis diagnosis could be considered an early prognostic biomarker in sepsis-associated ARDS patients. Future studies will be needed to evaluate the mechanistic links of this observation with the pathogenesis of ARDS.

Evaluation of Whole-Exome Enrichment Solutions: Lessons from the High-End of the Short-Read Sequencing Scale.

Whole-exome sequencing has become a popular technique in research and clinical settings, assisting in disease diagnosis and increasing the understanding of disease pathogenesis. In this study, we aimed to compare common enrichment capture solutions available in the market. Peripheral blood-purified DNA samples were enriched with SureSelectQXT V6 (Agilent) and various Illumina solutions: TruSeq DNA Nano, TruSeq DNA Exome, Nextera DNA Exome, and Illumina DNA Prep with Enrichment, and sequenced on a HiSeq 4000. We found that their percentage of duplicate reads was as much as 2 times higher than previously reported values for the previous HiSeq series. SureSelectQXT and Illumina DNA Prep with Enrichment showed the best average on-target coverage, which improved when off-target regions were included. At high coverage levels and in shared bases, these two solutions and TruSeq DNA Exome provided three of the best performances. With respect to the number of small variants detected, SureSelectQXT presented the lowest number of detected variants in target regions. When off-target regions were considered, its ability equalized to other solutions. Our results show SureSelectQXT and Illumina DNA Prep with Enrichment to be the best enrichment capture solutions.

Sepsis-associated acute respiratory distress syndrome in individuals of European ancestry: a genome-wide association study.

BACKGROUND: Acute respiratory distress syndrome (ARDS) is a lung inflammatory process caused mainly by sepsis. Most previous studies that identified genetic risks for ARDS focused on candidates with biological relevance. We aimed to identify novel genetic variants associated with ARDS susceptibility and to provide complementary functional evidence of their effect in gene regulation. METHODS: We did a case-control genome-wide association study (GWAS) of 1935 European individuals, using patients with sepsis-associated ARDS as cases and patients with sepsis without ARDS as controls. The discovery stage included 672 patients admitted into a network of Spanish intensive care units between January, 2002, and January, 2017. The replication stage comprised 1345 individuals from two independent datasets from the MESSI cohort study (Sep 22, 2008-Nov 30, 2017; USA) and the VISEP (April 1, 2003-June 30, 2005) and MAXSEP (Oct 1, 2007-March 31, 2010) trials of the SepNet study (Germany). Results from discovery and replication stages were meta-analysed to identify association signals. We then used RNA sequencing data from lung biopsies, in-silico analyses, and luciferase reporter assays to assess the functionallity of associated variants. FINDINGS: We identified a novel genome-wide significant association with sepsis-associated ARDS susceptibility (rs9508032, odds ratio [OR] 0·61, 95% CI 0·41-0·91, p=5·18 × 10-8) located within the Fms-related tyrosine kinase 1 (FLT1) gene, which encodes vascular endothelial growth factor receptor 1 (VEGFR-1). The region containing the sentinel variant and its best proxies acted as a silencer for the FLT1 promoter, and alleles with protective effects in ARDS further reduced promoter activity (p=0·0047). A literature mining of all previously described ARDS genes validated the association of vascular endothelial growth factor A (VEGFA; OR 0·55, 95% CI 0·41-0·73; p=4·69 × 10-5). INTERPRETATION: A common variant within the FLT1 gene is associated with sepsis-associated ARDS. Our findings support a role for the vascular endothelial growth factor signalling pathway in ARDS pathogenesis and identify VEGFR-1 as a potential therapeutic target. FUNDING: Instituto de Salud Carlos III, European Regional Development Funds, Instituto Tecnológico y de Energías Renovables.