PCSK6 and Survival in Idiopathic Pulmonary Fibrosis

Rationale: Idiopathic pulmonary fibrosis (IPF) is a devastating disease characterized by limited treatment options and high mortality. A better understanding of the molecular drivers of IPF progression is needed. Objectives: To identify and validate molecular determinants of IPF survival. Methods: A staged genome-wide association study was performed using paired genomic and survival data. Stage I cases were drawn from centers across the United States and Europe and stage II cases from Vanderbilt University. Cox proportional hazards regression was used to identify gene variants associated with differential transplantation-free survival (TFS). Stage I variants with nominal significance (P < 5 × 10-5) were advanced for stage II testing and meta-analyzed to identify those reaching genome-wide significance (P < 5 × 10-8). Downstream analyses were performed for genes and proteins associated with variants reaching genome-wide significance. Measurements and Main Results: After quality controls, 1,481 stage I cases and 397 stage II cases were included in the analysis. After filtering, 9,075,629 variants were tested in stage I, with 158 meeting advancement criteria. Four variants associated with TFS with consistent effect direction were identified in stage II, including one in an intron of PCSK6 (proprotein convertase subtilisin/kexin type 6) reaching genome-wide significance (hazard ratio, 4.11 [95% confidence interval, 2.54-6.67]; P = 9.45 × 10-9). PCSK6 protein was highly expressed in IPF lung parenchyma. PCSK6 lung staining intensity, peripheral blood gene expression, and plasma concentration were associated with reduced TFS. Conclusions: We identified four novel variants associated with IPF survival, including one in PCSK6 that reached genome-wide significance. Downstream analyses suggested that PCSK6 protein plays a potentially important role in IPF progression.

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.

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&gt;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.

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⁻⁴; 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⁻⁴; 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.

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.

Admixture mapping analysis reveals differential genetic ancestry associated with Chagas disease susceptibility in the Colombian population

Chagas disease is an infection caused by the parasite Trypanosoma cruzi, endemic in Latino America. Leveraging the three-way admixture between Native American (AMR), European (EUR) and African (AFR) populations in Latin Americans, we aimed to better understand the genetic basis of Chagas disease by performing an admixture mapping study in a Colombian population. A two-stage study was conducted, and subjects were classified as seropositive and seronegative for T. cruzi. In stage 1, global and local ancestries were estimated using reference data from the 1000 Genomes Project (1KGP) and local ancestry associations were performed by logistic regression models. The AMR ancestry showed a protective association with Chagas disease within the Major Histocompatibility Complex region (OR = 0.74, 95%CI = 0.66-0.83, lowest p-value = 4.53x10-8). The fine mapping assessment on imputed genotypes combining data from stage 1 and 2 from an independent Colombian cohort, revealed nominally associated variants in high linkage disequilibrium with the top signal (rs2032134, OR = 0.93, 95%CI = 0.90-0.97, p-value = 3.54x10-4) in the previously associated locus. To assess ancestry-specific adaptive signals, a selective sweep scan in an AMR reference population from 1KGP together with an in silico functional analysis highlighted the Tripartite Motif family and the Human Leukocyte Antigen (HLA) genes, with crucial role in the immune response against pathogens. Furthermore, these analyses emphasized the macrophages, neutrophils, and eosinophils, as key players in the defense against T. cruzi. This first admixture mapping study in Chagas disease provided novel insights underlying the host immune response in the pathogenesis of this neglected disease.

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 SureSelect^QXT 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. SureSelect^QXT 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, SureSelect^QXT 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 SureSelect^QXT and Illumina DNA Prep with Enrichment to be the best enrichment capture solutions.

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.

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.

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.

Novel idiopathic pulmonary fibrosis susceptibility variants revealed by deep sequencing

Background

Specific common and rare single nucleotide variants (SNVs) increase the likelihood of developing sporadic idiopathic pulmonary fibrosis (IPF). We performed target-enriched sequencing on three loci previously identified by a genome-wide association study to gain a deeper understanding of the full spectrum of IPF genetic risk and performed a two-stage case–control association study.

Methods

A total of 1.7 Mb of DNA from 181 IPF patients was deep sequenced (>100×) across 11p15.5, 14q21.3 and 17q21.31 loci. Comparisons were performed against 501 unrelated controls and replication studies were assessed in 3968 subjects.

Results

36 SNVs were associated with IPF susceptibility in the discovery stage (p<5.0×10⁻⁸). After meta-analysis, the strongest association corresponded to rs35705950 (p=9.27×10⁻⁵⁷) located upstream from the mucin 5B gene (MUC5B). Additionally, a novel association was found for two co-inherited low-frequency SNVs (<5%) in MUC5AC, predicting a missense amino acid change in mucin 5AC (lowest p=2.27×10⁻²²). Conditional and haplotype analyses in 11p15.5 supported the existence of an additional contribution of MUC5AC variants to IPF risk.

Conclusions

This study reinforces the significant IPF associations of these loci and implicates MUC5AC as another key player in IPF susceptibility.

Deep sequencing of genome-wide association study hits identified novel low-frequency variants associated with IPF susceptibility.http://bit.ly/2IF4AT8

Whole-Exome Sequencing Identifies Somatic Mutations Associated With Mortality in Metastatic Clear Cell Kidney Carcinoma

Clear cell renal cell carcinoma (ccRCC) is among the most aggressive histologic subtypes of kidney cancer, representing about 3% of all human cancers. Patients at stage IV have nearly 60% of mortality in 2–3 years after diagnosis. To date, most ccRCC studies have used DNA microarrays and targeted sequencing of a small set of well-established, commonly altered genes. An exception is the large multi-omics study of The Cancer Genome Atlas Kidney Renal Clear Cell Carcinoma (TCGA-KIRC), which identified new ccRCC genes based on whole exome-sequencing (WES) data, and molecular prognostic signatures based on transcriptomics, epigenetics and proteomics data. Applying WES to simultaneously interrogate virtually all exons in the human genome for somatic variation, here we analyzed the burden of coding somatic mutations in metastatic ccRCC primary tumors, and its association with patient mortality from cancer, in patients who received VEGF receptor-targeting drugs as the first-line therapy. To this end, we sequenced the exomes of ten tumor–normal pairs of ccRCC patient tissues from primary biopsies at >100× mean depth and called somatic coding variation. Mutation burden analysis prioritized 138 genes linked to patient mortality. A gene set enrichment analysis evidenced strong statistical support for the abundance of genes involved in the development of kidney cancer (p = 2.31 × 10⁻⁹) and carcinoma (p = 1.22 × 10⁻⁵), with 49 genes having direct links with kidney cancer according to the published records. Two of these genes, SIPA1L2 and EIF3A, demonstrated independent associations with mortality in TCGA-KIRC project data. Besides, three mutational signatures were found to be operative in the tumor exomes, one of which (COSMIC signature 12) has not been previously reported in ccRCC. Selection analysis yielded no detectable evidence of overall positive or negative selection, with the exome-wide number of nonsynonymous substitutions per synonymous site reflecting largely neutral tumor evolution. Despite the limited sample size, our results provide evidence for candidate genes where somatic mutation burden is tentatively associated with patient mortality in metastatic ccRCC, offering new potential pharmacological targets and a basis for further validation studies.