Association study of human leukocyte antigen variants and idiopathic pulmonary fibrosis.

Introduction: Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF, and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Owing to the important role that the human leukocyte antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk. Methods: We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case-control studies of IPF (comprising 5159 cases and 27 459 controls, including a prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold of p<4.5×10-4 and a posterior probability of replication >90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report. Results: The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies. Conclusion: Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

Genome-wide SNP-sex interaction analysis of susceptibility to idiopathic pulmonary fibrosis.

Background: Idiopathic pulmonary fibrosis (IPF) is a chronic lung condition that is more prevalent in males than females. The reasons for this are not fully understood, with differing environmental exposures due to historically sex-biased occupations, or diagnostic bias, being possible explanations. To date, over 20 independent genetic variants have been identified to be associated with IPF susceptibility, but these have been discovered when combining males and females. Our aim was to test for the presence of sex-specific associations with IPF susceptibility and assess whether there is a need to consider sex-specific effects when evaluating genetic risk in clinical prediction models for IPF. Methods: We performed genome-wide single nucleotide polymorphism (SNP)-by-sex interaction studies of IPF risk in six independent IPF case-control studies and combined them using inverse-variance weighted fixed effect meta-analysis. In total, 4,561 cases (1,280 females and 2,281 males) and 23,500 controls (8,360 females and 14,528 males) of European genetic ancestry were analysed. We used polygenic risk scores (PRS) to assess differences in genetic risk prediction between males and females. Findings: Three independent genetic association signals were identified. All showed a consistent direction of effect across all individual IPF studies and an opposite direction of effect in IPF susceptibility between females and males. None had been previously identified in IPF susceptibility genome-wide association studies (GWAS). The predictive accuracy of the PRSs were similar between males and females, regardless of whether using combined or sex-specific GWAS results. Interpretation: We prioritised three genetic variants whose effect on IPF risk may be modified by sex, however these require further study. We found no evidence that the predictive accuracy of common SNP-based PRSs varies significantly between males and females.

Shared genetic risk factors: implications for treatment of idiopathic pulmonary fibrosis and systemic hypertension.

The discovery of shared genetic associations, with sometimes different directions of effect, has implications for drug target discovery for idiopathic pulmonary fibrosis and systemic hypertension https://bit.ly/45IqbRv.

Alternative splicing in lung influences COVID-19 severity and respiratory diseases.

Alternative splicing generates functional diversity in isoforms, impacting immune response to infection. Here, we evaluate the causal role of alternative splicing in COVID-19 severity and susceptibility by applying two-sample Mendelian randomization to cis-splicing quantitative trait loci and the results from COVID-19 Host Genetics Initiative. We identify that alternative splicing in lung, rather than total expression of OAS1, ATP11A, DPP9 and NPNT, is associated with COVID-19 severity. MUC1 and PMF1 splicing is associated with COVID-19 susceptibility. Colocalization analyses support a shared genetic mechanism between COVID-19 severity with idiopathic pulmonary fibrosis at the ATP11A and DPP9 loci, and with chronic obstructive lung diseases at the NPNT locus. Last, we show that ATP11A, DPP9, NPNT, and MUC1 are highly expressed in lung alveolar epithelial cells, both in COVID-19 uninfected and infected samples. These findings clarify the importance of alternative splicing in lung for COVID-19 and respiratory diseases, providing isoform-based targets for drug discovery.

Association study of human leukocyte antigen (HLA) variants and idiopathic pulmonary fibrosis.

Introduction: Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Due to the important role that the Human Leukocyte Antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk. Methods: We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case-control studies of IPF (comprising a total of 5,159 cases and 27,459 controls, including the prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold p<4.5×10-4 and a posterior probability of replication >90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report. Results: The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies. Conclusion: Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

Genetic Determinants of the Acute Respiratory Distress Syndrome.

Acute respiratory distress syndrome (ARDS) is a life-threatening lung condition that arises from multiple causes, including sepsis, pneumonia, trauma, and severe coronavirus disease 2019 (COVID-19). Given the heterogeneity of causes and the lack of specific therapeutic options, it is crucial to understand the genetic and molecular mechanisms that underlie this condition. The identification of genetic risks and pharmacogenetic loci, which are involved in determining drug responses, could help enhance early patient diagnosis, assist in risk stratification of patients, and reveal novel targets for pharmacological interventions, including possibilities for drug repositioning. Here, we highlight the basis and importance of the most common genetic approaches to understanding the pathogenesis of ARDS and its critical triggers. We summarize the findings of screening common genetic variation via genome-wide association studies and analyses based on other approaches, such as polygenic risk scores, multi-trait analyses, or Mendelian randomization studies. We also provide an overview of results from rare genetic variation studies using Next-Generation Sequencing techniques and their links with inborn errors of immunity. Lastly, we discuss the genetic overlap between severe COVID-19 and ARDS by other causes.

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.

Longitudinal lung function and gas transfer in individuals with idiopathic pulmonary fibrosis: a genome-wide association study.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is an incurable lung disease characterised by progressive scarring leading to alveolar stiffness, reduced lung capacity, and impeded gas transfer. We aimed to identify genetic variants associated with declining lung capacity or declining gas transfer after diagnosis of IPF. METHODS: We did a genome-wide meta-analysis of longitudinal measures of forced vital capacity (FVC) and diffusing capacity of the lung for carbon monoxide (DLCO) in individuals diagnosed with IPF. Individuals were recruited to three studies between June, 1996, and August, 2017, from across centres in the US, UK, and Spain. Suggestively significant variants were investigated further in an additional independent study (CleanUP-IPF). All four studies diagnosed cases following American Thoracic Society/European Respiratory Society guidelines. Variants were defined as significantly associated if they had a meta-analysis p<5 × 10-8 when meta-analysing across all discovery and follow-up studies, had consistent direction of effects across all four studies, and were nominally significant (p<0·05) in each study. FINDINGS: 1329 individuals with a total of 5216 measures were included in the FVC analysis. 975 individuals with a total of 3361 measures were included in the DLCO analysis. For the discovery genome-wide analyses, 7 611 174 genetic variants were included in the FVC analysis and 7 536 843 in the DLCO analysis. One variant (rs115982800) located in an antisense RNA gene for protein kinase N2 (PKN2) showed a genome-wide significant association with FVC decline (-140 mL/year per risk allele [95% CI -180 to -100]; p=9·14 × 10-12). INTERPRETATION: Our analysis identifies a genetic variant associated with disease progression, which might highlight a new biological mechanism for IPF. We found that PKN2, a Rho and Rac effector protein, is the most likely gene of interest from this analysis. PKN2 inhibitors are currently in development and signify a potential novel therapeutic approach for IPF. FUNDING: Action for Pulmonary Fibrosis, Medical Research Council, Wellcome Trust, and National Institutes of Health National Heart, Lung, and Blood Institute.

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.

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.

Genome-wide association study across five cohorts identifies five novel loci associated with idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a chronic lung condition with poor survival times. We previously published a genome-wide meta-analysis of IPF risk across three studies with independent replication of associated variants in two additional studies. To maximise power and to generate more accurate effect size estimates, we performed a genome-wide meta-analysis across all five studies included in the previous IPF risk genome-wide association studies. We used the distribution of effect sizes across the five studies to assess the replicability of the results and identified five robust novel genetic association signals implicating mTOR (mammalian target of rapamycin) signalling, telomere maintenance and spindle assembly genes in IPF risk.

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.

Leveraging global multi-ancestry meta-analysis in the study of idiopathic pulmonary fibrosis genetics.

The research of rare and devastating orphan diseases, such as idiopathic pulmonary fibrosis (IPF) has been limited by the rarity of the disease itself. The prognosis is poor-the prevalence of IPF is only approximately four times the incidence, limiting the recruitment of patients to trials and studies of the underlying biology. Global biobanking efforts can dramatically alter the future of IPF research. We describe a large-scale meta-analysis of IPF, with 8,492 patients and 1,355,819 population controls from 13 biobanks around the globe. Finally, we combine this meta-analysis with the largest available meta-analysis of IPF, reaching 11,160 patients and 1,364,410 population controls. We identify seven novel genome-wide significant loci, only one of which would have been identified if the analysis had been limited to European ancestry individuals. We observe notable pleiotropy across IPF susceptibility and severe COVID-19 infection and note an unexplained sex-heterogeneity effect at the strongest IPF locus MUC5B.

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.

Genetic Ancestry Inference and Its Application for the Genetic Mapping of Human Diseases.

Admixed populations arise when two or more ancestral populations interbreed. As a result of this admixture, the genome of admixed populations is defined by tracts of variable size inherited from these parental groups and has particular genetic features that provide valuable information about their demographic history. Diverse methods can be used to derive the ancestry apportionment of admixed individuals, and such inferences can be leveraged for the discovery of genetic loci associated with diseases and traits, therefore having important biomedical implications. In this review article, we summarize the most common methods of global and local genetic ancestry estimation and discuss the use of admixture mapping studies in human diseases.

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 [Odds ratio (OR) = 0.74, 95% confidence interval (CI) = 0.66-0.83, lowest P-value = 4.53 × 10-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.54 × 10-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 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 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.

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.