Sofer, Tamar

ABSTRACT Obstructive sleep apnea (OSA) is a common, heritable disorder with diverse clinical presentations and etiologies. We conducted a genome-wide association meta-analysis of 492,107 individuals, including 46,028 OSA cases, identifying 14 genome-wide significant loci, eight of which are novel. Analyses adjusting for body mass index (BMI) revealed three loci independent of obesity, implicating alternative biological pathways contributing to disease. Integrative functional analyses, including chromatin interaction mapping, fine-mapping, and eQTL colocalization, prioritized candidate effector genes. Notably, implicated genes in chondrocytes were associated with craniofacial morphology in both mouse and human multivariate genotype-phenotype mapping, supporting a role for craniofacial structure in OSA risk. These findings highlight key genetic pathways (obesity-related, neurological, and craniofacial) underlying the development of OSA, offering new insights into its complex etiology.

Abstract Background:Chronic obstructive pulmonary disease (COPD) is a major global health burden, with genome-wide association studies explaining only a fraction of its heritability. Recent advancements, including large-scale, multi-ancestry whole-genome sequencing and a focus on severe, early-onset cases, offer promising avenues to identify additional COPD-related variants. Furthermore, recent studies show that adjusting for polygenic risk scores (PRS) enhances statistical power to detect rare variants, advancing efforts to understand COPD's genetic basis. Methods:We performed a whole genome sequencing association study of severe, early-onset COPD (FEV1 < 50%, age ≤ 55), FEV1 and FEV1/FVC in eight population-based studies (ARIC, CARDIA, CFS, CHS, FHS, HCHS/SOL, JHS, MESA) and seven COPD-enriched studies (COPDGene, ECLIPSE, Boston Early-Onset COPD Study, LTRC, TCGS (Poland), ICGN, SPIROMICS). We performed single-variant analyses using regenie software, applying a genome-wide significance threshold of P < 5.0 x 10-8. For rare variant set (MAF <1%) analysis, we annotated variants with loss of function or deleterious missense using WGS Annotator and prediction tools LOFTEE, SIFT, and PolyPhen2, then analyzed these variants with SKATO-ACAT method in regenie, adjusting for covariates including sex, sequencing center, and principal components of ancestry. We additionally evaluated adjusting the analyses for PRSs for lung function traits. Results: We analyzed 53,028 participants including 34,751 European, 12,508 African, 1,088 Latino American, and 674 East Asian ancestry. Across cohorts, we analyzed 5,063 severe (FEV1 < 50%) COPD cases, among whom 1,103 were early-onset (age ≤ 55). Excluding variants within 1 Mbp of known variants, single-variant analysis identified 68 novel associations (MAC ≥ 30), with 8 common associations (MAF ≥ 1%) in 6 genomic regions. Our most significant novel association for severe early-onset COPD was rs74297390 (intergenic to POM121L12, MAF = 9.6%, Beta = 0.41, P = 3.5 x 10-8). For FEV1 and FEV1/FVC ratio, variants near DKK2 and AKR1E2 were significantly associated and were replicated in an independent dataset (Shrine et al., 2019). DKK2 is an antagonist in the Wnt pathway. In a preliminary variant set analysis, CHRM3-AS1 was associated with severe early-onset COPD, GPAT with FEV1, and COL8A1I and HMCN1 with FEV1/FVC. With PRS adjustment, ST6GALNAC5 was associated with FEV1/FVC. Conclusions: In a whole genome sequencing study of lung function and severe early-onset COPD, we identified new loci for COPD-related phenotypes, including two novel loci that replicated in an independent dataset. Ongoing work involves variant set-based analyses with additional covariate adjustments, as well as PRS adjustments for binary traits.

Background: Risk for venous thromboembolism has a strong genetic component. Whole genome sequencing from the TOPMed program (Trans-Omics for Precision Medicine) allowed us to look for new associations, particularly rare variants missed by standard genome-wide association studies. Methods: The 3793 cases and 7834 controls (11.6% of cases were individuals of African, Hispanic/Latino, or Asian ancestry) were analyzed using a single variant approach and an aggregate gene-based approach using our primary filter (included only loss-of-function and missense variants predicted to be deleterious) and our secondary filter (included all missense variants). Results: Single variant analyses identified associations at 5 known loci. Aggregate gene-based analyses identified only PROC (odds ratio, 6.2 for carriers of rare variants; P =7.4×10 −14 ) when using our primary filter. Employing our secondary variant filter led to a smaller effect size at PROC (odds ratio, 3.8; P =1.6×10 −14 ), while excluding variants found only in rare isoforms led to a larger one (odds ratio, 7.5). Different filtering strategies improved the signal for 2 other known genes: PROS1 became significant (minimum P =1.8×10 −6 with the secondary filter), while SERPINC1 did not (minimum P =4.4×10 −5 with minor allele frequency <0.0005). Results were largely the same when restricting the analyses to include only unprovoked cases; however, one novel gene, MS4A1 , became significant ( P =4.4×10 −7 using all missense variants with minor allele frequency <0.0005). Conclusions: Here, we have demonstrated the importance of using multiple variant filtering strategies, as we detected additional genes when filtering variants based on their predicted deleteriousness, frequency, and presence on the most expressed isoforms. Our primary analyses did not identify new candidate loci; thus larger follow-up studies are needed to replicate the novel MS4A1 locus and to identify additional rare variation associated with venous thromboembolism.

Background: The Trans-Omics for Precision Medicine (TOPMed) program is performing whole-genome sequencing across multiple studies, including longitudinal cohorts, case-control, and case-only sample sets. Detecting the effects of low frequency variants requires large sample sizes, which can only be achieved by combining data across diverse study designs, including matching of case-only sample sets with controls from other studies. Here we present a strategy for combined analysis of venous thromboembolism (VTE) case status in 3 longitudinal cohorts and 3 case-only sample sets, in the context of whole-genome association studies. Methods: For each cohort study, we sampled ‘pure’ controls (without replacement) from the risk set of each incident case, within strata defined by sex, ancestry and birth-cohort. Pure controls had no event throughout their period of observation. A case’s risk set was defined as controls with no prior VTE history and under observation through an age at least as old as the case. For case-only sample sets, controls for each case were sampled from a cohort study, using the same risk set definition. Because of limited overlap in birth years between the cohort studies and the case-only sample sets, this matching was done within strata defined only by sex and ancestry group. Mixed model logistic regression will be used to account for relatedness as a random effect. Although conditional logistic regression is not practical for whole-genome association studies, case-control matching is implicitly recognized by a fixed effect for age-at-event for each matched set (1 case and >=1 matched controls). Additional fixed effects will include sample set and sex. We will also adjust for variations in case-control ratio among the matched sets. Results: In two cohort studies, we matched 1,231 cases to 4,820 controls (overall ratio = 1:3.9); only 500 controls and 5 cases could not be matched. For the case-only sample sets we matched 2,141 cases to 2,780 controls from one cohort study (overall ratio = 1:1.3); zero controls and only 14 cases could not be matched. Conclusions: We were successfully able to match nearly all cases to controls, and more than 90% of controls were also matched. By matching controls to cases based on age at event, we can account for the different risk of VTE by age clusters. Although this strategy does not provide an asymptotically unbiased estimated of the hazard ratio, compared to classical risk set sampling, it uses a large portion of the available data, it provides odds ratio estimates yielding the correct sign of association, and it reduces the potential influence of resampling subjects with rare variants. This strategy also enables combined analysis of multiple studies with different designs.