Whole-genome sequencing in 333,100 individuals reveals rare non-coding single variant and aggregate associations with height.

The role of rare non-coding variation in complex human phenotypes is still largely unknown. To elucidate the impact of rare variants in regulatory elements, we performed a whole-genome sequencing association analysis for height using 333,100 individuals from three datasets: UK Biobank (N = 200,003), TOPMed (N = 87,652) and All of Us (N = 45,445). We performed rare ( < 0.1% minor-allele-frequency) single-variant and aggregate testing of non-coding variants in regulatory regions based on proximal-regulatory, intergenic-regulatory and deep-intronic annotation. We observed 29 independent variants associated with height at P < 6 × 10 - 10 after conditioning on previously reported variants, with effect sizes ranging from -7cm to +4.7 cm. We also identified and replicated non-coding aggregate-based associations proximal to HMGA1 containing variants associated with a 5 cm taller height and of highly-conserved variants in MIR497HG on chromosome 17. We have developed an approach for identifying non-coding rare variants in regulatory regions with large effects from whole-genome sequencing data associated with complex traits.

Evidence of survival bias in the association between APOE-Є4 and age at ischemic stroke onset.

Introduction: Large genome-wide association studies (GWASs) using case-control study designs have now identified tens of loci associated with ischemic stroke (IS). As a complement to these studies, we performed GWAS in a case-only design to identify loci influencing the age at onset (AAO) of ischemic stroke. Methods: Analyses were conducted in a discovery cohort of 10,857 ischemic stroke cases using a linear regression framework. We meta-analyzed all SNPs with p-value <1 x 10-5 in a sexcombined or sex-stratified analysis using summary data from two additional replication cohorts. Results: In the women-only meta-analysis, we detected significant evidence for the association of AAO with rs429358, an exonic variant in apolipoprotein E (APOE) that encodes for the APOE-Є4 allele. Each copy of the rs429358:T>C allele was associated with a 1.29-year earlier stroke AAO (meta p-value = 2.48 x 10-11). This APOE variant has previously been associated with increased mortality and ischemic stroke AAO. We hypothesized that the association with AAO may reflect a survival bias attributable to an age-related decrease in mortality among APOE-Є4 carriers and have no association to stroke AAO per se. A simulation study showed that a variant associated with overall mortality might indeed be detected with an AAO analysis. A variant with a 2-fold increase in mortality risk would lead to an observed effect of AAO that is comparable to what we found. Discussion: In conclusion, we detected a robust association of the APOE locus with stroke AAO and provided simulations to suggest that this association may be unrelated to ischemic stroke per se but related to a general survival bias.

Contrasting association pattern of plasma low-density lipoprotein with white matter integrity in APOE4 carriers versus non-carriers.

Apolipoprotein E ε4 (APOE4) is a strong genetic risk factor of Alzheimer's disease and metabolic dysfunction. However, whether APOE4 and markers of metabolic dysfunction synergistically impact the deterioration of white matter (WM) integrity in older adults remains unknown. In the UK Biobank data, we conducted a multivariate analysis to investigate the interactions between APOE4 and 249 plasma metabolites (measured using nuclear magnetic resonance spectroscopy) with whole-brain WM integrity (measured by diffusion-weighted magnetic resonance imaging) in a cohort of 1917 older adults (aged 65.0-81.0 years; 52.4 % female). Although no main association was observed between either APOE4 or metabolites with WM integrity (adjusted P > 0.05), significant interactions between APOE4 and metabolites with WM integrity were identified. Among the examined metabolites, higher concentrations of low-density lipoprotein and very low-density lipoprotein were associated with a lower level of WM integrity (b=-0.12, CI=-0.14,-0.10) among APOE4 carriers. Conversely, among non-carriers, they were associated with a higher level of WM integrity (b=0.05, CI=0.04,0.07), demonstrating a significant moderation role of APOE4 (b =-0.18, CI=-0.20,-0.15, P<0.00001).

A Preliminary Study of Quantitative MRI Cartilage Loss Fraction and Its Association With Future Arthroplasty Using the Osteoarthritis Initiative Database.

Background and objective  Osteoarthritis (OA) is the most common arthritis in the world. Despite the high disease burden, there is no therapy to prevent, halt, or reverse OA, and many clinical trials relied on radiographic biomarkers for therapy response. It is important to identify patients with early OA who will eventually need arthroplasty, the end-stage treatment for osteoarthritis. This pilot study evaluates a novel MRI biomarker, cartilage loss fraction, for association with future arthroplasty and evaluates its feasibility of use and effect size estimates. Materials and methods Publicly available knee MRIs from the Osteoarthritis Initiative were used. A total of 38 participants with Kellgren-Lawrence (K-L) grade >1 and 38 participants with K-L grade ≤ 1 at enrollment were matched in age, sex, race, and BMI, and assessed for the degree of full-thickness cartilage loss, or cartilage loss fraction. Univariate conditional logistic regression analysis was performed for differences in cartilage loss fractions between groups. Receiver operating characteristic (ROC) curve analysis was performed to assess the association of MRI biomarkers and knee arthroplasty during the eight-year follow-up. Results The medial femoral condyle, medial tibial plateau, total, and two-year progression cartilage loss fractions were significantly higher in participants with K-L grade >1 (p < 0.01 for all) and showed high area under the curve (AUC) values on ROC analysis (812, 0.827, 0.917, and 0.933, respectively). These results were comparable or more strongly associated with other OA grading schemes. Conclusion MRI biomarker cartilage loss fractions are significantly higher in subjects with K-L grade >1 and show a strong association with arthroplasty. After further validation, cartilage loss fracture may be used to predict future arthroplasty.

A Guide to Understanding Mendelian Randomization Studies.

Epidemiology provides a powerful framework for characterizing exposure-disease relationships, but its utility for making causal inferences is limited because epidemiologic data are observational in nature and subject to biases stemming from undetected confounding variables and reverse causation. Mendelian randomization (MR) is an increasingly popular method used to circumvent these limitations. MR uses genetic variants, or instruments, as a natural experiment to proxy an exposure, thus allowing estimation of causal effects upon an outcome that are minimally affected by the usual biases present in epidemiologic studies. Notably, MR relies on three core assumptions related to the selection of the genetic instruments, and adherence to these assumptions must be carefully evaluated to assess the validity of the causal estimates. The goal of this review is to provide readers with a basic understanding of MR studies and how to read and evaluate them. Specifically, we outline the basics of how MR analysis is conducted, the assumptions underlying instrument selection, and how to assess the quality of MR studies.

The Impact of Conventional Stroke Risk Factors on Early and Late Onset Ischemic Stroke: a Mendelian Randomization Study.

Objective: Although stroke incidence is decreasing in older ages, it is increasing in young adults. While these divergent trends in stroke incidence are at least partially attributable to diverging prevalence trends in stoke risk factors, age-dependent differences in the impact of stroke risk factors on stroke may also contribute. To address this issue, we utilized Mendelian Randomization (MR) to assess differences in the association of stroke risk factors between early onset ischemic stroke (EOS) and late onset ischemic stroke (LOS). Methods: We employed a two-sample MR design with inverse variance weighting as the primary method of analysis. Using large publicly available genome-wide association summary results, we calculated MR estimates for conventional stroke risk factors (body mass index, total, HDL-and LDL-cholesterol, triglycerides, type 2 diabetes, systolic and diastolic blood pressure, and smoking) in EOS cases (onset 18-59 years, n = 6,728) and controls from the Early Onset Stroke Consortium and in LOS cases (onset ≥ 60 years, n = 9,272) and controls from the Stroke Genetics Network. We then compared odds ratios between EOS and LOS, stratified by TOAST subtypes, to determine if any differences observed between effect sizes could be attributed to differences in the distribution of stroke subtypes. Results: EOS was significantly associated with all risk factors except for total cholesterol levels, and LOS was associated with all risk factors except for triglyceride and total cholesterol levels. The associations of BMI, DBP, SBP, and HDL-cholesterol were significantly stronger in EOS than LOS (all p < 0.004). The differential distribution of stroke subtypes could not explain the difference in effect size observed between EOS and LOS. Conclusion: These results suggest that interventions targeted at lowering body mass index and blood pressure may be particularly important for reducing stroke risk in young adults.

A Linear Relationship between the Number of Cancers among First-degree Relatives and the Risk of Multiple Primary Cancers.

With advances in the early detection and treatment of cancer, the incidence of multiple primary cancers (MPC), or second primary cancers, has risen over time. Characterization of etiologic risk factors, including family history of cancer, within the general population is critical for assessing MPC risk in patients. We examined the association between family history of cancer among first-degree relatives and MPC risk in a prospective study of 139,958 participants from the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial. Cox proportional hazard models were used to calculate hazard ratios (HR) and 95% confidence intervals (95%CI), adjusting for potential confounders. Over a median follow-up of 16 years (interquartile range: 11-19 years), 6,170 participants were diagnosed with MPC. Having a family history of cancer increased the risk of MPC by 18% (HR=1.18, 95%CI: 1.12-1.24). A positive linear trend was observed between the reported number of cancers in the family history and MPC risk with HRs (95%CI) of 1.13 (1.07-1.20), 1.23 (1.14-1.33), 1.29 (1.15-1.45), and 1.42 (1.20-1.70) for 1, 2, 3, and 4+ cancers among first-degree relatives, respectively (Ptrend=2.36x10-13). No significant differences were observed by cancer histology or specific types of cancer reported in the family history. Our study demonstrates that family history of cancer is an important risk factor for the development of multiple primary cancers and that a comprehensive of assessment of the number of cancers reported among first-degree relatives may identify those at higher risk who may benefit from targeted cancer prevention and screening strategies.

A Linear Relationship between the Number of Cancers among First-Degree Relatives and the Risk of Multiple Primary Cancers.

With advances in the early detection and treatment of cancer, the incidence of multiple primary cancers (MPC) or second primary cancers has increased over time. Characterization of etiologic risk factors, including family history of cancer, within the general population is critical for assessing MPC risk in patients. We examined the association between family history of cancer among first-degree relatives and MPC risk in a prospective study of 139,958 participants from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Cox proportional hazard models were used to calculate HRs and 95% confidence intervals (95% CI), adjusting for potential confounders. Over a median follow-up of 16 years (IQR: 11-19 years), 6,170 participants were diagnosed with MPC. Having a family history of cancer increased the risk of MPC by 18% (HR, 1.18; 95% CI, 1.12-1.24). A positive linear trend was observed between the reported number of cancers in the family history and MPC risk with HRs (95% CI) of 1.13 (1.07-1.20), 1.23 (1.14-1.33), 1.29 (1.15-1.45), and 1.42 (1.20-1.70) for one, two, three, and four or more cancers among first-degree relatives, respectively (Ptrend = 2.36 × 10-13). No significant differences were observed by cancer histology or specific cancer types reported in the family history. Our study demonstrates that the family history of cancer is an important risk factor for the development of MPCs and that a comprehensive assessment of the number of cancers reported among first-degree relatives may identify those at higher risk who may benefit from targeted cancer prevention and screening strategies. Prevention Relevance: Our study makes a substantial contribution to the understanding of risk factors for MPCs in the general population. It demonstrates that individuals with a strong family history of cancer are at higher risk for MPCs and may benefit from more targeted cancer prevention and screening interventions.

Machine learning models for predicting blood pressure phenotypes by combining multiple polygenic risk scores.

We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1 to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8 to 5.1% (SBP) and 4.7 to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs. In summary, non-linear ML models improves BP prediction in models incorporating diverse populations.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

Genomic and phenotypic correlates of mosaic loss of chromosome Y in blood.

Mosaic loss of Y (mLOY) is the most common somatic chromosomal alteration detected in human blood. The presence of mLOY is associated with altered blood cell counts and increased risk of Alzheimer's disease, solid tumors, and other age-related diseases. We sought to gain a better understanding of genetic drivers and associated phenotypes of mLOY through analyses of whole genome sequencing of a large set of genetically diverse males from the Trans-Omics for Precision Medicine (TOPMed) program. This approach enabled us to identify differences in mLOY frequencies across populations defined by genetic similarity, revealing a higher frequency of mLOY in the European American (EA) ancestry group compared to those of Hispanic American (HA), African American (AA), and East Asian (EAS) ancestry. Further, we identified two genes ( CFHR1 and LRP6 ) that harbor multiple rare, putatively deleterious variants associated with mLOY susceptibility, show that subsets of human hematopoietic stem cells are enriched for activity of mLOY susceptibility variants, and that certain alleles on chromosome Y are more likely to be lost than others.

Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Acute pharmacodynamic responses to sitagliptin: Drug-induced increase in early insulin secretion in oral glucose tolerance test.

DPP4 inhibitors are widely prescribed as treatments for type 2 diabetes. Because drug responses vary among individuals, we initiated investigations to identify genetic variants associated with the magnitude of drug responses. Sitagliptin (100 mg) was administered to 47 healthy volunteers. Several endpoints were measured to assess clinically relevant responses - including the effect of sitagliptin on glucose and insulin levels during an oral glucose tolerance test (OGTT). This pilot study confirmed that sitagliptin (100 mg) decreased the area under the curve for glucose during an OGTT (p = 0.0003). Furthermore, sitagliptin promoted insulin secretion during the early portion of the OGTT as reflected by an increase in the ratio of plasma insulin at 30 min divided by plasma insulin at 60 min (T30:T60) from mean ± SEM 0.87 ± 0.05 to 1.62 ± 0.36 mU/L (p = 0.04). The magnitude of sitagliptin's effect on insulin secretion (as judged by the increase in the T30:T60 ratio for insulin) was correlated with the magnitude of sitagliptin-induced increase in the area under the curve for intact plasma GLP1 levels during the first hour of the OGTT. This study confirmed previously reported sex differences in glucose and insulin levels during an OGTT. Specifically, females exhibited higher levels of glucose and insulin at the 90-180 min time points. However, we did not detect significant sex-associated differences in the magnitude of sitagliptin-induced changes in T30:T60 ratios for either glucose or insulin. In conclusion, T30:T60 ratios for insulin and glucose during an OGTT provide useful indices to assess pharmacodynamic responses to DPP4 inhibitors.

Causal Factors for Osteoarthritis: A Scoping Review of Mendelian Randomization Studies.

OBJECTIVE: Mendelian randomization (MR) has increasingly been utilized as a tool for establishing causal relations between modifiable exposures and osteoarthritis (OA). The goal of this review was to summarize available MR studies of OA that evaluate the causal role of modifiable risk factors on OA. METHODS: This review was performed following the Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) Extension for Scoping Reviews model. We performed a literature search for relevant studies published before December 2021 across multiple databases using the search terms "osteoarthritis" and ("Mendelian randomization" or "polygenic risk score"). We reported the MR estimates of causal associations between exposures and OA and then assessed methodologic quality of abstracted studies according to their efforts to validate the three key MR assumptions. RESULTS: Our search identified 45 studies reporting on 141 exposure-association analyses. All studies performed a formal instrumental variable analysis to estimate the causal effect of exposure on OA. Causal associations (P < 0.05) were reported in 60 of these analyses representing 36 unique publications, and MR-Egger sensitivity analyses were performed in 45 of these analyses. MR studies provided support for causal associations of OA with increased levels of adiposity, coffee consumption, bone mineral density, and sleep disturbance, and decreased levels of serum calcium and low-density lipoprotein cholesterol. CONCLUSION: These results highlight the potential benefits of weight reduction and improvement of sleep quality to reduce the risk of OA and call for a better understanding of the relations of coffee consumption and serum calcium to OA risk.

Polygenic Risk Score Associates With Atherosclerotic Plaque Characteristics at Autopsy.

BACKGROUND: Polygenic risk scores (PRSs) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically. METHODS: From 4327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner for sudden death between 1994 and 2015, 2455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas, and thrombotic CAD. RESULTS: After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age, 48.8±14.7 years; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared with subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% versus 50.4%±38.7%; adjusted P<0.001) and a higher frequency of calcification (69.6% versus 35.8%; adjusted P=0.004) and thin-cap fibroatheroma (26.7% versus 9.5%; adjusted P=0.007). Even after adjustment for traditional CAD risk factors, subjects within the highest PRS quintile had higher odds of severe atherosclerosis (ie, ≥75% stenosis; adjusted odds ratio, 3.77 [95% CI, 2.10-6.78]; P<0.001) and plaque rupture (adjusted odds ratio, 4.05 [95% CI, 2.26-7.24]; P<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged ≤50 years (adjusted odds ratio, 4.08 [95% CI, 2.01-8.30]; P<0.001). No statistically significant associations were observed with plaque erosion after adjusting for covariates. CONCLUSIONS: This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects.

Genomic risk scores and oral contraceptive-associated ischemic stroke risk: a call for collaboration.

Background: Oral contraceptives (OCs) are generally safe but vascular risk factors increase OC-associated ischemic stroke risk. We performed a case-control study to evaluate whether a genomic risk score for ischemic stroke modifies OC-associated ischemic stroke risk. Methods: The Genetics of Early-Onset Stroke study includes 332 premenopausal women (136 arterial ischemic stroke cases and 196 controls) with data on estrogen-containing OC use within 30 days before the index event (for cases) or interview (for controls). Using a previously validated genetic risk score (metaGRS) for ischemic stroke based on 19 polygenic risk scores for stroke and stroke-associated risk factors, we stratified our combined case-control sample into tertiles of genomic risk. We evaluated the association between OC use and ischemic stroke within each tertile. We tested if the association between OC use and ischemic stroke depended on the genomic risk of stroke using logistic regression with an OC use × metaGRS interaction term. These analyses were performed with and without adjustment for smoking, hypertension, diabetes, coronary heart disease, and body mass index. Results: After adjustment for vascular risk factors, the odds ratio of OC use was 3.2 (1.7-6.3) overall and increased from the lower, middle, and upper tertile of genomic risk from 1.6 (0.5-5.4) to 2.5 (0.08-8.2) to 13.7 (3.8-67.3) respectively, and a p-value for interaction of 0.001. Conclusions: Our results suggest that genomic profile may modify the OC-associated ischemic stroke risk. Larger studies are warranted to determine whether a genomic risk score could be clinically useful in reducing OC-associated ischemic stroke.

Evidence of survival bias in the association between APOE-ϵ4 and age of ischemic stroke onset.

Large genome-wide association studies (GWAS) employing case-control study designs have now identified tens of loci associated with ischemic stroke (IS). As a complement to these studies, we performed GWAS in a case-only design to identify loci influencing age at onset (AAO) of ischemic stroke. Analyses were conducted in a Discovery cohort of 10,857 ischemic stroke cases using a linear regression framework. We meta-analyzed all SNPs with p-value < 1×10-5 in a sex-combined or sex-stratified analysis using summary data from two additional replication cohorts. In the women-only meta-analysis, we detected significant evidence for association of AAO with rs429358, an exonic variant in APOE that encodes for the APOE-ϵ4 allele. Each copy of the rs429358:T>C allele was associated with a 1.29 years earlier stroke AOO (meta p-value = 2.48×10-11). This APOE variant has previously been associated with increased mortality and ischemic stroke AAO. We hypothesized that the association with AAO may reflect a survival bias attributable to an age-related decline in mortality among APOE-ϵ4 carriers and have no association to stroke AAO per se. Using a simulation study, we found that a variant associated with overall mortality might indeed be detected with an AAO analysis. A variant with a two-fold increase on mortality risk would lead to an observed effect of AAO that is comparable to what we found. In conclusion, we detected a robust association of the APOE locus with stroke AAO and provided simulations to suggest that this association may be unrelated to ischemic stroke per se but related to a general survival bias.

Missense variants in SORT1 are associated with LDL-C in an Amish population.

Common noncoding variants at the human 1p13.3 locus associated with SORT1 expression are among those most strongly associated with low-density lipoprotein cholesterol (LDL-C) in human genome-wide association studies. However, validation studies in mice and cell lines have produced variable results regarding the directionality of the effect of SORT1 on LDL-C. This, together with the fact that the 1p13.3 variants are associated with expression of several genes, has raised the question of whether SORT1 is the causal gene at this locus. Using whole exome sequencing in members of an Amish population, we identified coding variants in SORT1 that are associated with increased (rs141749679, K302E) and decreased (rs149456022, Q225H) LDL-C. Further, analysis of plasma lipoprotein particle subclasses by ion mobility in a subset of rs141749679 (K302E) carriers revealed higher levels of large LDL particles compared to noncarriers. In contrast to the effect of these variants in the Amish, the sortilin K302E mutation introduced into a C57BL/6J mouse via CRISPR/Cas9 resulted in decreased non-high-density lipoprotein cholesterol, and the sortilin Q225H mutation did not alter cholesterol levels in mice. This is indicative of different effects of these mutations on cholesterol metabolism in the two species. To our knowledge, this is the first evidence that naturally occurring coding variants in SORT1 are associated with LDL-C, thus supporting SORT1 as the gene responsible for the association of the 1p13.3 locus with LDL-C.

APOE4 poses opposite effects of plasma LDL on white matter integrity in older adults.

INTRODUCTION: APOE4 is a strong genetic risk factor of Alzheimer's disease and is associated with changes in metabolism. However, the interactive relationship between APOE4 and plasma metabolites on the brain remains largely unknown. MEHODS: In the UK Biobank, we investigated the moderation effects of APOE4 on the relationship between 249 plasma metabolites derived from nuclear magnetic resonance spectroscopy on whole-brain white matter integrity, measured by fractional anisotropy using diffusion magnetic resonance imaging. RESULTS: The increase in the concentration of metabolites, mainly LDL and VLDL, is associated with a decrease in white matter integrity (b= -0.12, CI= [-0.14, -0.10]) among older APOE4 carriers, whereas an increase (b= 0.05, CI= [0.04, 0.07]) among non-carriers, implying a significant moderation effect of APOE4 (b= -0.18, CI= [-0.20,-0.15]). DISCUSSION: The results suggest that lipid metabolism functions differently in APOE4 carriers compared to non-carriers, which may inform the development of targeted interventions for APOE4 carriers to mitigate cognitive decline.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.