Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative.

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

Pulmonary emphysema subtypes defined by unsupervised machine learning on CT scans.

BACKGROUND: Treatment and preventative advances for chronic obstructive pulmonary disease (COPD) have been slow due, in part, to limited subphenotypes. We tested if unsupervised machine learning on CT images would discover CT emphysema subtypes with distinct characteristics, prognoses and genetic associations. METHODS: New CT emphysema subtypes were identified by unsupervised machine learning on only the texture and location of emphysematous regions on CT scans from 2853 participants in the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS), a COPD case-control study, followed by data reduction. Subtypes were compared with symptoms and physiology among 2949 participants in the population-based Multi-Ethnic Study of Atherosclerosis (MESA) Lung Study and with prognosis among 6658 MESA participants. Associations with genome-wide single-nucleotide-polymorphisms were examined. RESULTS: The algorithm discovered six reproducible (interlearner intraclass correlation coefficient, 0.91-1.00) CT emphysema subtypes. The most common subtype in SPIROMICS, the combined bronchitis-apical subtype, was associated with chronic bronchitis, accelerated lung function decline, hospitalisations, deaths, incident airflow limitation and a gene variant near DRD1, which is implicated in mucin hypersecretion (p=1.1 ×10-8). The second, the diffuse subtype was associated with lower weight, respiratory hospitalisations and deaths, and incident airflow limitation. The third was associated with age only. The fourth and fifth visually resembled combined pulmonary fibrosis emphysema and had distinct symptoms, physiology, prognosis and genetic associations. The sixth visually resembled vanishing lung syndrome. CONCLUSION: Large-scale unsupervised machine learning on CT scans defined six reproducible, familiar CT emphysema subtypes that suggest paths to specific diagnosis and personalised therapies in COPD and pre-COPD.

Dissecting the genetics of chronic mucus hypersecretion in smokers with and without COPD.

Smoking is a notorious risk factor for chronic mucus hypersecretion (CMH). CMH frequently occurs in chronic obstructive pulmonary disease (COPD). The question arises whether the same single-nucleotide polymorphisms (SNPs) are related to CMH in smokers with and without COPD. We performed two genome-wide association studies of CMH under an additive genetic model in male heavy smokers (≥20 pack-years) with COPD (n=849, 39.9% CMH) and without COPD (n=1348, 25.4% CMH), followed by replication and meta-analysis in comparable populations, and assessment of the functional relevance of significantly associated SNPs. Genome-wide association analysis of CMH in COPD and non-COPD subjects yielded no genome-wide significance after replication. In COPD, our top SNP (rs10461985, p=5.43×10(-5)) was located in the GDNF-AS1 gene that is functionally associated with the GDNF gene. Expression of GDNF in bronchial biopsies of COPD patients was significantly associated with CMH (p=0.007). In non-COPD subjects, four SNPs had a p-value <10(-5) in the meta-analysis, including a SNP (rs4863687) in the MAML3 gene, the T-allele showing modest association with CMH (p=7.57×10(-6), OR 1.48) and with significantly increased MAML3 expression in lung tissue (p=2.59×10(-12)). Our data suggest the potential for differential genetic backgrounds of CMH in individuals with and without COPD.

APOM and high-density lipoprotein cholesterol are associated with lung function and per cent emphysema.

Chronic obstructive pulmonary disease (COPD) is linked to cardiovascular disease; however, there are few studies on the associations of cardiovascular genes with COPD. We assessed the association of lung function with 2100 genes selected for cardiovascular diseases among 20 077 European-Americans and 6900 African-Americans. We performed replication of significant loci in the other racial group and an independent consortium of Europeans, tested the associations of significant loci with per cent emphysema and examined gene expression in an independent sample. We then tested the association of a related lipid biomarker with forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) ratio and per cent emphysema. We identified one new polymorphism for FEV1/FVC (rs805301) in European-Americans (p=1.3×10(-6)) and a second (rs707974) in the combined European-American and African-American analysis (p=1.38×10(-7)). Both single-nucleotide polymorphisms (SNPs) flank the gene for apolipoprotein M (APOM), a component of high-density lipoprotein (HDL) cholesterol. Both were replicated in an independent cohort. SNPs in a second gene related to apolipoprotein M and HDL, PCSK9, were associated with FEV1/FVC ratio among African-Americans. rs707974 was associated with per cent emphysema among European-Americans and African-Americans and APOM expression was related to FEV1/FVC ratio and per cent emphysema. Higher HDL levels were associated with lower FEV1/FVC ratio and greater per cent emphysema. These findings suggest a novel role for the apolipoprotein M/HDL pathway in the pathogenesis of COPD and emphysema.

The impact of damaging epilepsy and cardiac genetic variant burden in sudden death in the young.

BACKGROUND: Sudden unexpected death in children is a tragic event. Understanding the genetics of sudden death in the young (SDY) enables family counseling and cascade screening. The objective of this study was to characterize genetic variation in an SDY cohort using whole genome sequencing. METHODS: The SDY Case Registry is a National Institutes of Health/Centers for Disease Control and Prevention surveillance effort to discern the prevalence, causes, and risk factors for SDY. The SDY Case Registry prospectively collected clinical data and DNA biospecimens from SDY cases < 20 years of age. SDY cases were collected from medical examiner and coroner offices spanning 13 US jurisdictions from 2015 to 2019. The cohort included 211 children (median age 0.33 year; range 0-20 years), determined to have died suddenly and unexpectedly and from whom DNA biospecimens for DNA extractions and next-of-kin consent were ascertained. A control cohort consisted of 211 randomly sampled, sex- and ancestry-matched individuals from the 1000 Genomes Project. Genetic variation was evaluated in epilepsy, cardiomyopathy, and arrhythmia genes in the SDY and control cohorts. American College of Medical Genetics/Genomics guidelines were used to classify variants as pathogenic or likely pathogenic. Additionally, pathogenic and likely pathogenic genetic variation was identified using a Bayesian-based artificial intelligence (AI) tool. RESULTS: The SDY cohort was 43% European, 29% African, 3% Asian, 16% Hispanic, and 9% with mixed ancestries and 39% female. Six percent of the cohort was found to harbor a pathogenic or likely pathogenic genetic variant in an epilepsy, cardiomyopathy, or arrhythmia gene. The genomes of SDY cases, but not controls, were enriched for rare, potentially damaging variants in epilepsy, cardiomyopathy, and arrhythmia-related genes. A greater number of rare epilepsy genetic variants correlated with younger age at death. CONCLUSIONS: While damaging cardiomyopathy and arrhythmia genes are recognized contributors to SDY, we also observed an enrichment in epilepsy-related genes in the SDY cohort and a correlation between rare epilepsy variation and younger age at death. These findings emphasize the importance of considering epilepsy genes when evaluating SDY.

Peripheral Blood Mononuclear Cell Gene Expression Associated with Pulmonary Microvascular Perfusion: The Multi-Ethnic Study of Atherosclerosis Chronic Obstructive Pulmonary Disease.

Rationale: Chronic obstructive pulmonary disease (COPD) and emphysema are associated with endothelial damage and altered pulmonary microvascular perfusion. The molecular mechanisms underlying these changes are poorly understood in patients, in part because of the inaccessibility of the pulmonary vasculature. Peripheral blood mononuclear cells (PBMCs) interact with the pulmonary endothelium. Objectives: To test the association between gene expression in PBMCs and pulmonary microvascular perfusion in COPD. Methods: The Multi-Ethnic Study of Atherosclerosis (MESA) COPD Study recruited two independent samples of COPD cases and controls with ⩾10 pack-years of smoking history. In both samples, pulmonary microvascular blood flow, pulmonary microvascular blood volume, and mean transit time were assessed on contrast-enhanced magnetic resonance imaging, and PBMC gene expression was assessed by microarray. Additional replication was performed in a third sample with pulmonary microvascular blood volume measures on contrast-enhanced dual-energy computed tomography. Differential expression analyses were adjusted for age, gender, race/ethnicity, educational attainment, height, weight, smoking status, and pack-years of smoking. Results: The 79 participants in the discovery sample had a mean age of 69 ± 6 years, 44% were female, 25% were non-White, 34% were current smokers, and 66% had COPD. There were large PBMC gene expression signatures associated with pulmonary microvascular perfusion traits, with several replicated in the replication sets with magnetic resonance imaging (n = 47) or dual-energy contrast-enhanced computed tomography (n = 157) measures. Many of the identified genes are involved in inflammatory processes, including nuclear factor-κB and chemokine signaling pathways. Conclusions: PBMC gene expression in nuclear factor-κB, inflammatory, and chemokine signaling pathways was associated with pulmonary microvascular perfusion in COPD, potentially offering new targetable candidates for novel therapies.

Rare variant analyses validate known ALS genes in a multi-ethnic population and identifies ANTXR2 as a candidate in PLS.

Background: Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease affecting over 30,000 people in the United States. It is characterized by the progressive decline of the nervous system that leads to the weakening of muscles which impacts physical function. Approximately, 15% of individuals diagnosed with ALS have a known genetic variant that contributes to their disease. As therapies that slow or prevent symptoms, such as antisense oligonucleotides, continue to develop, it is important to discover novel genes that could be targets for treatment. Additionally, as cohorts continue to grow, performing analyses in ALS subtypes, such as primary lateral sclerosis (PLS), becomes possible due to an increase in power. These analyses could highlight novel pathways in disease manifestation. Methods: Building on our previous discoveries using rare variant association analyses, we conducted rare variant burden testing on a substantially larger cohort of 6,970 ALS patients from a large multi-ethnic cohort as well as 166 PLS patients, and 22,524 controls. We used intolerant domain percentiles based on sub-region Residual Variation Intolerance Score (subRVIS) that have been described previously in conjunction with gene based collapsing approaches to conduct burden testing to identify genes that associate with ALS and PLS. Results: A gene based collapsing model showed significant associations with SOD1, TARDBP, and TBK1 (OR=19.18, p = 3.67 × 10-39; OR=4.73, p = 2 × 10-10; OR=2.3, p = 7.49 × 10-9, respectively). These genes have been previously associated with ALS. Additionally, a significant novel control enriched gene, ALKBH3 (p = 4.88 × 10-7), was protective for ALS in this model. An intolerant domain based collapsing model showed a significant improvement in identifying regions in TARDBP that associated with ALS (OR=10.08, p = 3.62 × 10-16). Our PLS protein truncating variant collapsing analysis demonstrated significant case enrichment in ANTXR2 (p=8.38 × 10-6). Conclusions: In a large multi-ethnic cohort of 6,970 ALS patients, rare variant burden testing validated known ALS genes and identified a novel potentially protective gene, ALKBH3. A first-ever analysis in 166 patients with PLS found a candidate association with loss-of-function mutations in ANTXR2.

Rare variant analyses validate known ALS genes in a multi-ethnic population and identifies ANTXR2 as a candidate in PLS.

Background: Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease affecting over 30,000 people in the United States. It is characterized by the progressive decline of the nervous system that leads to the weakening of muscles which impacts physical function. Approximately, 15% of individuals diagnosed with ALS have a known genetic variant that contributes to their disease. As therapies that slow or prevent symptoms, such as antisense oligonucleotides, continue to develop, it is important to discover novel genes that could be targets for treatment. Additionally, as cohorts continue to grow, performing analyses in ALS subtypes, such as primary lateral sclerosis (PLS), becomes possible due to an increase in power. These analyses could highlight novel pathways in disease manifestation. Methods: Building on our previous discoveries using rare variant association analyses, we conducted rare variant burden testing on a substantially larger cohort of 6,970 ALS patients from a large multi-ethnic cohort as well as 166 PLS patients, and 22,524 controls. We used intolerant domain percentiles based on sub-region Residual Variation Intolerance Score (subRVIS) that have been described previously in conjunction with gene based collapsing approaches to conduct burden testing to identify genes that associate with ALS and PLS. Results: A gene based collapsing model showed significant associations with SOD1, TARDBP, and TBK1 (OR=19.18, p = 3.67 × 10-39; OR=4.73, p = 2 × 10-10; OR=2.3, p = 7.49 × 10-9, respectively). These genes have been previously associated with ALS. Additionally, a significant novel control enriched gene, ALKBH3 (p = 4.88 × 10-7), was protective for ALS in this model. An intolerant domain based collapsing model showed a significant improvement in identifying regions in TARDBP that associated with ALS (OR=10.08, p = 3.62 × 10-16). Our PLS protein truncating variant collapsing analysis demonstrated significant case enrichment in ANTXR2 (p=8.38 × 10-6). Conclusions: In a large multi-ethnic cohort of 6,970 ALS patients, rare variant burden testing validated known ALS genes and identified a novel potentially protective gene, ALKBH3. A first-ever analysis in 166 patients with PLS found a candidate association with loss-of-function mutations in ANTXR2.

The impact of damaging epilepsy and cardiac genetic variant burden in sudden death in the young.

Background: Sudden unexpected death in children is a tragic event. Understanding the genetics of sudden death in the young (SDY) enables family counseling and cascade screening. The objective of this study was to characterize genetic variation in an SDY cohort using whole genome sequencing. Methods: The SDY Case Registry is a National Institutes of Health/Centers for Disease Control surveillance effort to discern the prevalence, causes, and risk factors for SDY. The SDY Case Registry prospectively collected clinical data and DNA biospecimens from SDY cases <20 years of age. SDY cases were collected from medical examiner and coroner offices spanning 13 US jurisdictions from 2015-2019. The cohort included 211 children (mean age 1 year; range 0-20 years), determined to have died suddenly and unexpectedly and in whom DNA biospecimens and next-of-kin consent were ascertained. A control cohort consisted of 211 randomly sampled, sex-and ancestry-matched individuals from the 1000 Genomes Project. Genetic variation was evaluated in epilepsy, cardiomyopathy and arrhythmia genes in the SDY and control cohorts. American College of Medical Genetics/Genomics guidelines were used to classify variants as pathogenic or likely pathogenic. Additionally, genetic variation predicted to be damaging was identified using a Bayesian-based artificial intelligence (AI) tool. Results: The SDY cohort was 42% European, 30% African, 17% Hispanic, and 11% with mixed ancestries, and 39% female. Six percent of the cohort was found to harbor a pathogenic or likely pathogenic genetic variant in an epilepsy, cardiomyopathy or arrhythmia gene. The genomes of SDY cases, but not controls, were enriched for rare, damaging variants in epilepsy, cardiomyopathy and arrhythmia-related genes. A greater number of rare epilepsy genetic variants correlated with younger age at death. Conclusions: While damaging cardiomyopathy and arrhythmia genes are recognized contributors to SDY, we also observed an enrichment in epilepsy-related genes in the SDY cohort, and a correlation between rare epilepsy variation and younger age at death. These findings emphasize the importance of considering epilepsy genes when evaluating SDY.

Association of cardiovascular health and epigenetic age acceleration.

BACKGROUND: Cardiovascular health (CVH) has been defined by the American Heart Association (AHA) as the presence of the "Life's Simple 7" ideal lifestyle and clinical factors. CVH is known to predict longevity and freedom from cardiovascular disease, the leading cause of death for women in the United States. DNA methylation markers of aging have been aggregated into a composite epigenetic age score, which is associated with cardiovascular morbidity and mortality. However, it is unknown whether poor CVH is associated with acceleration of aging as measured by DNA methylation markers in epigenetic age. METHODS AND RESULTS: We performed a cross-sectional analysis of racially/ethnically diverse post-menopausal women enrolled in the Women's Health Initiative cohort recruited between 1993 and 1998. Epigenetic age acceleration (EAA) was calculated using DNA methylation data on a subset of participants and the published Horvath and Hannum methods for intrinsic and extrinsic EAA. CVH was calculated using the AHA measures of CVH contributing to a 7-point score. We examined the association between CVH score and EAA using linear regression modeling adjusting for self-reported race/ethnicity and education. Among the 2,170 participants analyzed, 50% were white and mean age was 64 (7 SD) years. Higher or more favorable CVH scores were associated with lower extrinsic EAA (~ 6 months younger age per 1 point higher CVH score, p < 0.0001), and lower intrinsic EAA (3 months younger age per 1 point higher CVH score, p < 0.028). CONCLUSIONS: These cross-sectional observations suggest a possible mechanism by which ideal CVH is associated with greater longevity.

Genomic Context Differs Between Human Dilated Cardiomyopathy and Hypertrophic Cardiomyopathy.

Background Inherited cardiomyopathies display variable penetrance and expression, and a component of phenotypic variation is genetically determined. To evaluate the genetic contribution to this variable expression, we compared protein coding variation in the genomes of those with hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM). Methods and Results Nonsynonymous single-nucleotide variants (nsSNVs) were ascertained using whole genome sequencing from familial cases of HCM (n=56) or DCM (n=70) and correlated with echocardiographic information. Focusing on nsSNVs in 102 genes linked to inherited cardiomyopathies, we correlated the number of nsSNVs per person with left ventricular measurements. Principal component analysis and generalized linear models were applied to identify the probability of cardiomyopathy type as it related to the number of nsSNVs in cardiomyopathy genes. The probability of having DCM significantly increased as the number of cardiomyopathy gene nsSNVs per person increased. The increase in nsSNVs in cardiomyopathy genes significantly associated with reduced left ventricular ejection fraction and increased left ventricular diameter for individuals carrying a DCM diagnosis, but not for those with HCM. Resampling was used to identify genes with aberrant cumulative allele frequencies, identifying potential modifier genes for cardiomyopathy. Conclusions Participants with DCM had more nsSNVs per person in cardiomyopathy genes than participants with HCM. The nsSNV burden in cardiomyopathy genes did not correlate with the probability or manifestation of left ventricular measures in HCM. These findings support the concept that increased variation in cardiomyopathy genes creates a genetic background that predisposes to DCM and increased disease severity.

Exome sequencing in obsessive-compulsive disorder reveals a burden of rare damaging coding variants.

Obsessive-compulsive disorder (OCD) affects 1-2% of the population, and, as with other complex neuropsychiatric disorders, it is thought that rare variation contributes to its genetic risk. In this study, we performed exome sequencing in the largest OCD cohort to date (1,313 total cases, consisting of 587 trios, 41 quartets and 644 singletons of affected individuals) and describe contributions to disease risk from rare damaging coding variants. In case-control analyses (n = 1,263/11,580), the most significant single-gene result was observed in SLITRK5 (odds ratio (OR) = 8.8, 95% confidence interval 3.4-22.5, P = 2.3 × 10-6). Across the exome, there was an excess of loss of function (LoF) variation specifically within genes that are LoF-intolerant (OR = 1.33, P = 0.01). In an analysis of trios, we observed an excess of de novo missense predicted damaging variants relative to controls (OR = 1.22, P = 0.02), alongside an excess of de novo LoF mutations in LoF-intolerant genes (OR = 2.55, P = 7.33 × 10-3). These data support a contribution of rare coding variants to OCD genetic risk.

Pathogenic and Uncertain Genetic Variants Have Clinical Cardiac Correlates in Diverse Biobank Participants.

Background Genome sequencing coupled with electronic heath record data can uncover medically important genetic variation. Interpretation of rare genetic variation and its role in mediating cardiovascular phenotypes is confounded by variants of uncertain significance. Methods and Results We analyzed the whole genome sequence of 900 racially and ethnically diverse biobank participants selected from a single US center. Participants were equally divided among European, African, Hispanic, and mixed races/ethnicities. We evaluated the American College of Medical Genetics and Genomics medically actionable list of 59 genes, focusing on the cardiac genes. Variation was interpreted using the most recent reports in ClinVar, a database of medically relevant human variation. We identified 19 individuals with pathogenic or likely pathogenic variants in cardiac actionable genes (2%) and found evidence of related clinical correlates in the electronic health record. Participants of African ancestry, compared with those of European ancestry, had more variants of uncertain significance in the medically actionable genes including the 30 cardiac actionable genes, even when normalized to total variant count per person. Longitudinal measures of left ventricle size from ≈400 biobank participants (1723 patient-years) were correlated with genetic findings. The presence of ≥1 uncertain variant in the actionable cardiac genes and a cardiomyopathy diagnosis correlated with increased left ventricular internal diameter in diastole and in systole. In particular, MYBPC3 was identified as a gene with excess variants of uncertain significance. Conclusions These data indicate that a subset of uncertain genetic variants may confer risk and should not be considered benign.

Meta-analysis of exome array data identifies six novel genetic loci for lung function.

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Genome-wide association analysis identifies six new loci associated with forced vital capacity.

Forced vital capacity (FVC), a spirometric measure of pulmonary function, reflects lung volume and is used to diagnose and monitor lung diseases. We performed genome-wide association study meta-analysis of FVC in 52,253 individuals from 26 studies and followed up the top associations in 32,917 additional individuals of European ancestry. We found six new regions associated at genome-wide significance (P < 5 × 10(-8)) with FVC in or near EFEMP1, BMP6, MIR129-2-HSD17B12, PRDM11, WWOX and KCNJ2. Two loci previously associated with spirometric measures (GSTCD and PTCH1) were related to FVC. Newly implicated regions were followed up in samples from African-American, Korean, Chinese and Hispanic individuals. We detected transcripts for all six newly implicated genes in human lung tissue. The new loci may inform mechanisms involved in lung development and the pathogenesis of restrictive lung disease.