Organ aging signatures in the plasma proteome track health and disease.

Animal studies show aging varies between individuals as well as between organs within an individual1-4, but whether this is true in humans and its effect on age-related diseases is unknown. We utilized levels of human blood plasma proteins originating from specific organs to measure organ-specific aging differences in living individuals. Using machine learning models, we analysed aging in 11 major organs and estimated organ age reproducibly in five independent cohorts encompassing 5,676 adults across the human lifespan. We discovered nearly 20% of the population show strongly accelerated age in one organ and 1.7% are multi-organ agers. Accelerated organ aging confers 20-50% higher mortality risk, and organ-specific diseases relate to faster aging of those organs. We find individuals with accelerated heart aging have a 250% increased heart failure risk and accelerated brain and vascular aging predict Alzheimer's disease (AD) progression independently from and as strongly as plasma pTau-181 (ref. 5), the current best blood-based biomarker for AD. Our models link vascular calcification, extracellular matrix alterations and synaptic protein shedding to early cognitive decline. We introduce a simple and interpretable method to study organ aging using plasma proteomics data, predicting diseases and aging effects.

Blood-Based Proteomics for Adult-Onset Focal Dystonias.

OBJECTIVES: The adult-onset focal dystonias are characterized by over-active muscles leading to abnormal movements. For most cases, the etiology and pathogenesis remain unknown. In the current study, unbiased proteomics methods were used to identify potential changes in blood plasma proteins. METHODS: A large-scale unbiased proteomics screen was used to compare proteins (N = 6,345) in blood plasma of normal healthy controls (N = 49) with adult-onset focal dystonia (N = 143) consisting of specific subpopulations of cervical dystonia (N = 45), laryngeal dystonia (N = 49), and blepharospasm (N = 49). Pathway analyses were conducted to identify relevant biological pathways. Finally, protein changes were used to build a prediction model for dystonia. RESULTS: After correction for multiple comparisons, 15 proteins were associated with adult-onset focal dystonia. Subgroup analyses revealed some proteins were shared across the dystonia subgroups while others were unique to 1 subgroup. The top biological pathways involved changes in the immune system, metal ion transport, and reactive oxygen species. A 4-protein model showed high accuracy in discriminating control individuals from dystonia cases [average area under the curve (AUC) = 0.89]. INTERPRETATION: These studies provide novel insights into the etiopathogenesis of dystonia, as well as novel potential biomarkers. ANN NEUROL 2024;96:110-120.

Harmonization of CSF and imaging biomarkers in Alzheimer's disease: Need and practical applications for genetics studies and preclinical classification.

In Alzheimer's disease (AD) research, cerebrospinal fluid (CSF) Amyloid beta (Aß), Tau and pTau are the most accepted and well validated biomarkers. Several methods and platforms exist to measure those biomarkers, leading to challenges in combining data across studies. Thus, there is a need to identify methods that harmonize and standardize these values. We used a Z-score based approach to harmonize CSF and amyloid imaging data from multiple cohorts and compared GWAS results using this approach with currently accepted methods. We also used a generalized mixture model to calculate the threshold for biomarker-positivity. Based on our findings, our normalization approach performed as well as meta-analysis and did not lead to any spurious results. In terms of dichotomization, cutoffs calculated with this approach were very similar to those reported previously. These findings show that the Z-score based harmonization approach can be applied to heterogeneous platforms and provides biomarker cut-offs consistent with the classical approaches without requiring any additional data.

Alzheimer's polygenic risk scores are associated with cognitive phenotypes in Down syndrome.

INTRODUCTION: This study aimed to investigate the influence of the overall Alzheimer's disease (AD) genetic architecture on Down syndrome (DS) status, cognitive measures, and cerebrospinal fluid (CSF) biomarkers. METHODS: AD polygenic risk scores (PRS) were tested for association with DS-related traits. RESULTS: The AD risk PRS was associated with disease status in several cohorts of sporadic late- and early-onset and familial late-onset AD, but not in familial early-onset AD or DS. On the other hand, lower DS Mental Status Examination memory scores were associated with higher PRS, independent of intellectual disability and APOE (PRS including APOE, PRSAPOE , p = 2.84 × 10-4 ; PRS excluding APOE, PRSnonAPOE , p = 1.60 × 10-2 ). PRSAPOE exhibited significant associations with Aß42, tTau, pTau, and Aß42/40 ratio in DS. DISCUSSION: These data indicate that the AD genetic architecture influences cognitive and CSF phenotypes in DS adults, supporting common pathways that influence memory decline in both traits. HIGHLIGHTS: Examination of the polygenic risk of AD in DS presented here is the first of its kind. AD PRS influences memory aspects in DS individuals, independently of APOE genotype. These results point to an overlap between the genes and pathways that leads to AD and those that influence dementia and memory decline in the DS population. APOE ε4 is linked to DS cognitive decline, expanding cognitive insights in adults.

Post-GWAS multiomic functional investigation of the TNIP1 locus in Alzheimer's disease highlights a potential role for GPX3.

INTRODUCTION: Recent genome-wide association studies (GWAS) have reported a genetic association with Alzheimer's disease (AD) at the TNIP1/GPX3 locus, but the mechanism is unclear. METHODS: We used cerebrospinal fluid (CSF) proteomics data to test (n = 137) and replicate (n = 446) the association of glutathione peroxidase 3 (GPX3) with CSF biomarkers (including amyloid and tau) and the GWAS-implicated variants (rs34294852 and rs871269). RESULTS: CSF GPX3 levels decreased with amyloid and tau positivity (analysis of variance P = 1.5 × 10-5) and higher CSF phosphorylated tau (p-tau) levels (P = 9.28 × 10-7). The rs34294852 minor allele was associated with decreased GPX3 (P = 0.041). The replication cohort found associations of GPX3 with amyloid and tau positivity (P = 2.56 × 10-6) and CSF p-tau levels (P = 4.38 × 10-9). DISCUSSION: These results suggest variants in the TNIP1 locus may affect the oxidative stress response in AD via altered GPX3 levels. HIGHLIGHTS: Cerebrospinal fluid (CSF) glutathione peroxidase 3 (GPX3) levels decreased with amyloid and tau positivity and higher CSF phosphorylated tau. The minor allele of rs34294852 was associated with lower CSF GPX3. levels when also controlling for amyloid and tau category. GPX3 transcript levels in the prefrontal cortex were lower in Alzheimer's disease than controls. rs34294852 is an expression quantitative trait locus for GPX3 in blood, neutrophils, and microglia.

Large-scale proteome and metabolome analysis of CSF implicates altered glucose and carbon metabolism and succinylcarnitine in Alzheimer's disease.

INTRODUCTION: A hallmark of Alzheimer's disease (AD) is the aggregation of proteins (amyloid beta [A] and hyperphosphorylated tau [T]) in the brain, making cerebrospinal fluid (CSF) proteins of particular interest. METHODS: We conducted a CSF proteome-wide analysis among participants of varying AT pathology (n = 137 participants; 915 proteins) with nine CSF biomarkers of neurodegeneration and neuroinflammation. RESULTS: We identified 61 proteins significantly associated with the AT category (P < 5.46 × 10-5 ) and 636 significant protein-biomarker associations (P < 6.07 × 10-6 ). Proteins from glucose and carbon metabolism pathways were enriched among amyloid- and tau-associated proteins, including malate dehydrogenase and aldolase A, whose associations with tau were replicated in an independent cohort (n = 717). CSF metabolomics identified and replicated an association of succinylcarnitine with phosphorylated tau and other biomarkers. DISCUSSION: These results implicate glucose and carbon metabolic dysregulation and increased CSF succinylcarnitine levels with amyloid and tau pathology in AD. HIGHLIGHTS: Cerebrospinal fluid (CSF) proteome enriched for extracellular, neuronal, immune, and protein processing. Glucose/carbon metabolic pathways enriched among amyloid/tau-associated proteins. Key glucose/carbon metabolism protein associations independently replicated. CSF proteome outperformed other omics data in predicting amyloid/tau positivity. CSF metabolomics identified and replicated a succinylcarnitine-phosphorylated tau association.

GEM: scalable and flexible gene-environment interaction analysis in millions of samples.

MOTIVATION: Gene-environment interaction (GEI) studies are a general framework that can be used to identify genetic variants that modify the effects of environmental, physiological, lifestyle or treatment effects on complex traits. Moreover, accounting for GEIs can enhance our understanding of the genetic architecture of complex diseases and traits. However, commonly used statistical software programs for GEI studies are either not applicable to testing certain types of GEI hypotheses or have not been optimized for use in large samples. RESULTS: Here, we develop a new software program, GEM (Gene-Environment interaction analysis in Millions of samples), which supports the inclusion of multiple GEI terms, adjustment for GEI covariates and robust inference, while allowing multi-threading to reduce computation time. GEM can conduct GEI tests as well as joint tests of genetic main and interaction effects for both continuous and binary phenotypes. Through simulations, we demonstrate that GEM scales to millions of samples while addressing limitations of existing software programs. We additionally conduct a gene-sex interaction analysis on waist-hip ratio in 352 768 unrelated individuals from the UK Biobank, identifying 24 novel loci in the joint test that have not previously been reported in combined or sex-specific analyses. Our results demonstrate that GEM can facilitate the next generation of large-scale GEI studies and help advance our understanding of the genetic architecture of complex diseases and traits. AVAILABILITY AND IMPLEMENTATION: GEM is freely available as an open source project at https://github.com/large-scale-gxe-methods/GEM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome-wide meta-analysis for Alzheimer's disease cerebrospinal fluid biomarkers.

Amyloid-beta 42 (Aß42) and phosphorylated tau (pTau) levels in cerebrospinal fluid (CSF) reflect core features of the pathogenesis of Alzheimer's disease (AD) more directly than clinical diagnosis. Initiated by the European Alzheimer & Dementia Biobank (EADB), the largest collaborative effort on genetics underlying CSF biomarkers was established, including 31 cohorts with a total of 13,116 individuals (discovery n = 8074; replication n = 5042 individuals). Besides the APOE locus, novel associations with two other well-established AD risk loci were observed; CR1 was shown a locus for Aß42 and BIN1 for pTau. GMNC and C16orf95 were further identified as loci for pTau, of which the latter is novel. Clustering methods exploring the influence of all known AD risk loci on the CSF protein levels, revealed 4 biological categories suggesting multiple Aß42 and pTau related biological pathways involved in the etiology of AD. In functional follow-up analyses, GMNC and C16orf95 both associated with lateral ventricular volume, implying an overlap in genetic etiology for tau levels and brain ventricular volume.

Multi-ancestry genome-wide gene-sleep interactions identify novel loci for blood pressure.

Long and short sleep duration are associated with elevated blood pressure (BP), possibly through effects on molecular pathways that influence neuroendocrine and vascular systems. To gain new insights into the genetic basis of sleep-related BP variation, we performed genome-wide gene by short or long sleep duration interaction analyses on four BP traits (systolic BP, diastolic BP, mean arterial pressure, and pulse pressure) across five ancestry groups in two stages using 2 degree of freedom (df) joint test followed by 1df test of interaction effects. Primary multi-ancestry analysis in 62,969 individuals in stage 1 identified three novel gene by sleep interactions that were replicated in an additional 59,296 individuals in stage 2 (stage 1 + 2 Pjoint < 5 × 10-8), including rs7955964 (FIGNL2/ANKRD33) that increases BP among long sleepers, and rs73493041 (SNORA26/C9orf170) and rs10406644 (KCTD15/LSM14A) that increase BP among short sleepers (Pint < 5 × 10-8). Secondary ancestry-specific analysis identified another novel gene by long sleep interaction at rs111887471 (TRPC3/KIAA1109) in individuals of African ancestry (Pint = 2 × 10-6). Combined stage 1 and 2 analyses additionally identified significant gene by long sleep interactions at 10 loci including MKLN1 and RGL3/ELAVL3 previously associated with BP, and significant gene by short sleep interactions at 10 loci including C2orf43 previously associated with BP (Pint < 10-3). 2df test also identified novel loci for BP after modeling sleep that has known functions in sleep-wake regulation, nervous and cardiometabolic systems. This study indicates that sleep and primary mechanisms regulating BP may interact to elevate BP level, suggesting novel insights into sleep-related BP regulation.

Gene-educational attainment interactions in a multi-ancestry genome-wide meta-analysis identify novel blood pressure loci.

Educational attainment is widely used as a surrogate for socioeconomic status (SES). Low SES is a risk factor for hypertension and high blood pressure (BP). To identify novel BP loci, we performed multi-ancestry meta-analyses accounting for gene-educational attainment interactions using two variables, "Some College" (yes/no) and "Graduated College" (yes/no). Interactions were evaluated using both a 1 degree of freedom (DF) interaction term and a 2DF joint test of genetic and interaction effects. Analyses were performed for systolic BP, diastolic BP, mean arterial pressure, and pulse pressure. We pursued genome-wide interrogation in Stage 1 studies (N = 117 438) and follow-up on promising variants in Stage 2 studies (N = 293 787) in five ancestry groups. Through combined meta-analyses of Stages 1 and 2, we identified 84 known and 18 novel BP loci at genome-wide significance level (P < 5 × 10-8). Two novel loci were identified based on the 1DF test of interaction with educational attainment, while the remaining 16 loci were identified through the 2DF joint test of genetic and interaction effects. Ten novel loci were identified in individuals of African ancestry. Several novel loci show strong biological plausibility since they involve physiologic systems implicated in BP regulation. They include genes involved in the central nervous system-adrenal signaling axis (ZDHHC17, CADPS, PIK3C2G), vascular structure and function (GNB3, CDON), and renal function (HAS2 and HAS2-AS1, SLIT3). Collectively, these findings suggest a role of educational attainment or SES in further dissection of the genetic architecture of BP.

Development and Applications of CRISPR/Cas9-Based Genome Editing in Lactobacillus.

Lactobacillus, a genus of lactic acid bacteria, plays a crucial function in food production preservation, and probiotics. It is particularly important to develop new Lactobacillus strains with superior performance by gene editing. Currently, the identification of its functional genes and the mining of excellent functional genes mainly rely on the traditional gene homologous recombination technology. CRISPR/Cas9-based genome editing is a rapidly developing technology in recent years. It has been widely applied in mammalian cells, plants, yeast, and other eukaryotes, but less in prokaryotes, especially Lactobacillus. Compared with the traditional strain improvement methods, CRISPR/Cas9-based genome editing can greatly improve the accuracy of Lactobacillus target sites and achieve traceless genome modification. The strains obtained by this technology may even be more efficient than the traditional random mutation methods. This review examines the application and current issues of CRISPR/Cas9-based genome editing in Lactobacillus, as well as the development trend of CRISPR/Cas9-based genome editing in Lactobacillus. In addition, the fundamental mechanisms of CRISPR/Cas9-based genome editing are also presented and summarized.

Efficient gene-environment interaction tests for large biobank-scale sequencing studies.

Complex human diseases are affected by genetic and environmental risk factors and their interactions. Gene-environment interaction (GEI) tests for aggregate genetic variant sets have been developed in recent years. However, existing statistical methods become rate limiting for large biobank-scale sequencing studies with correlated samples. We propose efficient Mixed-model Association tests for GEne-Environment interactions (MAGEE), for testing GEI between an aggregate variant set and environmental exposures on quantitative and binary traits in large-scale sequencing studies with related individuals. Joint tests for the aggregate genetic main effects and GEI effects are also developed. A null generalized linear mixed model adjusting for covariates but without any genetic effects is fit only once in a whole genome GEI analysis, thereby vastly reducing the overall computational burden. Score tests for variant sets are performed as a combination of genetic burden and variance component tests by accounting for the genetic main effects using matrix projections. The computational complexity is dramatically reduced in a whole genome GEI analysis, which makes MAGEE scalable to hundreds of thousands of individuals. We applied MAGEE to the exome sequencing data of 41,144 related individuals from the UK Biobank, and the analysis of 18,970 protein coding genes finished within 10.4 CPU hours.

Identifying blood pressure loci whose effects are modulated by multiple lifestyle exposures.

Although multiple lifestyle exposures simultaneously impact blood pressure (BP) and cardiovascular health, most analysis so far has considered each single lifestyle exposure (e.g., smoking) at a time. Here, we exploit gene-multiple lifestyle exposure interactions to find novel BP loci. For each of 6,254 Framingham Heart Study participants, we computed lifestyle risk score (LRS) value by aggregating the risk of four lifestyle exposures (smoking, alcohol, education, and physical activity) on BP. Using the LRS, we performed genome-wide gene-environment interaction analysis in systolic and diastolic BP using the joint 2 degree of freedom (DF) and 1 DF interaction tests. We identified one genome-wide significant (p < 5 × 10-8 ) and 11 suggestive (p < 1 × 10-6 ) loci. Gene-environment analysis using single lifestyle exposures identified only one of the 12 loci. Nine of the 12 BP loci detected were novel. Loci detected by the LRS were located within or nearby genes with biologically plausible roles in the pathophysiology of hypertension, including KALRN, VIPR2, SNX1, and DAPK2. Our results suggest that simultaneous consideration of multiple lifestyle exposures in gene-environment interaction analysis can identify additional loci missed by single lifestyle approaches.

A multi-ancestry genome-wide study incorporating gene-smoking interactions identifies multiple new loci for pulse pressure and mean arterial pressure.

Elevated blood pressure (BP), a leading cause of global morbidity and mortality, is influenced by both genetic and lifestyle factors. Cigarette smoking is one such lifestyle factor. Across five ancestries, we performed a genome-wide gene-smoking interaction study of mean arterial pressure (MAP) and pulse pressure (PP) in 129 913 individuals in stage 1 and follow-up analysis in 480 178 additional individuals in stage 2. We report here 136 loci significantly associated with MAP and/or PP. Of these, 61 were previously published through main-effect analysis of BP traits, 37 were recently reported by us for systolic BP and/or diastolic BP through gene-smoking interaction analysis and 38 were newly identified (P < 5 × 10-8, false discovery rate < 0.05). We also identified nine new signals near known loci. Of the 136 loci, 8 showed significant interaction with smoking status. They include CSMD1 previously reported for insulin resistance and BP in the spontaneously hypertensive rats. Many of the 38 new loci show biologic plausibility for a role in BP regulation. SLC26A7 encodes a chloride/bicarbonate exchanger expressed in the renal outer medullary collecting duct. AVPR1A is widely expressed, including in vascular smooth muscle cells, kidney, myocardium and brain. FHAD1 is a long non-coding RNA overexpressed in heart failure. TMEM51 was associated with contractile function in cardiomyocytes. CASP9 plays a central role in cardiomyocyte apoptosis. Identified only in African ancestry were 30 novel loci. Our findings highlight the value of multi-ancestry investigations, particularly in studies of interaction with lifestyle factors, where genomic and lifestyle differences may contribute to novel findings.

Lysobacter profundi sp. nov., isolated from freshwater sediment and reclassification of Lysobacter panaciterrae as Luteimonas panaciterrae comb. nov.

A polyphasic taxonomic study was carried out on strains CHu50b-3-2T and CHu40b-3-1 isolated from a 67 cm-long sediment core collected from the Daechung Reservoir at a water depth of 17 m, Daejeon, Republic of Korea. The cells of the strains were Gram-stain-negative, non-spore-forming, non-motile and rod-shaped. Comparative 16S rRNA gene sequence studies showed a clear affiliation of two strains with γ-Proteobacteria, which showed the highest pairwise sequence similarities to Lysobacter hankyongensis KTce-2T (96.5 %), Lysobacter pocheonensis Gsoil193T (96.3 %), Lysobacter ginsengisoli Gsoil 357T (96.1 %), Lysobacter solanacearum T20R-70T (96.1 %), Lysobacter brunescens KCTC 12130T (95.4 %) and Lysobacter capsici YC5194T (95.3 %). The phylogenetic analysis based on 16S rRNA gene sequences showed that the strains formed a clear phylogenetic lineage with the genus Lysobacter. The major fatty acids were identified as summed feature 9 (iso-C17 : 1 ω9c and/or C18 : 1 10-methyl), iso-C15 : 0, iso-C16 : 0 and iso-C17 : 0. The respiratory quinone was identified as ubiquinone Q-8. The major polar lipids were phosphatidylglycerol, diphosphatidylglycerol, phosphatidylethanolamine and an unidentified phospholipid. The genomic DNA G+C content was determined to be 66.8 mol% (genome) for strain CHu50b-3-2T and 66.4 mol% (HPLC) for strain CHu40b-3-1. Based on the combined genotypic and phenotypic data, we propose that strains CHu50b-3-2T and CHu40b-3-1 represent a novel species of the genus Lysobacter, for which the name Lysobacter profundi sp. nov. is proposed. The type strain is CHu50b-3-2T (=KCTC 72973T=CCTCC AB 2019129T). Besides Lysobacter panaciterrae Gsoil 068T formed a phylogenetic group together with strain Luteimonas aquatica RIB1-20T (EF626688) that is clearly separated from all other known Lysobacter strains. Based on the phylogenetic relationships together with fatty acid compositions, Lysobacter panaciterrae Gsoil 068T should be reclassified as a member of the genus Luteimonas: Luteimonas aquatica comb. nov. (type strain Gsoil 068T=KCTC 12601T=DSM 17927T).

Lacisediminimonas profundi gen. nov., sp. nov., a member of the family Oxalobacteraceae isolated from freshwater sediment.

A novel Gram-stain-negative bacterial strain, CHu64-6-4T, was isolated from a 67-cm-long sediment core collected from the Daechung Reservoir at a water depth of 17 m, Daejeon, Republic of Korea. The cells of strain CHu64-6-4T were aerobic nonmotile and formed colorless colonies on R2A agar. The phylogenetic analysis based on 16S rRNA gene sequencing indicated that the strain formed a separate lineage within the family Oxalobacteraceae, exhibiting 97.2% and 97.1% 16S rRNA gene sequence similarities to Glaciimonas singularis and Paraherbaspirillum soli, respectively. Strain CHu64-6-4T showed less than 74.4% average nucleotide identity compared to the type strains of related genera within the family Oxalobacteraceae. In the UPGMA dendrogram based on the ANI values of genomic sequences, strain CHu64-6-4T formed an evolutionary lineage independent of the genera Glaciimonas and some other taxa. The chemotaxonomic results showed Q-8 as the predominant respiratory ubiquinone, phosphatidylglycerol, diphosphatidylglycerol, and phosphatidylethnolamine as the major polar lipids, Summed Feature 3 (C16:1ω7c and/or iso-C15:0 2-OH), C16:0, and C18:1ω7c as the major fatty acids, and a DNA G+C content of 62.1 mol%. The combined genotypic and phenotypic data showed that strain CHu64-6-4T could be distinguished from all genera within the family Oxalobacteraceae and represents a novel genus, Lacisediminimonas profundi gen. nov., with the name Lacisediminimonas profundi sp. nov., in the family Oxalobacteraceae. The type strain is CHu64-6-4T (=KCTC 62287T=JCM 32676T).

The Promise of Selecting Individuals from the Extremes of Exposure in the Analysis of Gene-Physical Activity Interactions.

BACKGROUND: Dichotomization using the lower quartile as cutoff is commonly used for harmonizing heterogeneous physical activity (PA) measures across studies. However, this may create misclassification and hinder discovery of new loci. OBJECTIVES: This study aimed to evaluate the performance of selecting individuals from the extremes of the exposure (SIEE) as an alternative approach to reduce such misclassification. METHOD: For systolic and diastolic blood pressure in the Framingham Heart Study, we performed a genome-wide association study with gene-PA interaction analysis using three PA variables derived by SIEE and two other dichotomization approaches. We compared number of loci detected and overlap with loci found using a quantitative PA variable. In addition, we performed simulation studies to assess bias, false discovery rates (FDR), and power under synergistic/antagonistic genetic effects in exposure groups and in the presence/absence of measurement error. RESULTS: In the empirical analysis, SIEE's performance was neither the best nor the worst. In most simulation scenarios, SIEE was consistently outperformed in terms of FDR and power. Particularly, in a scenario characterized by antagonistic effects and measurement error, SIEE had the least bias and highest power. CONCLUSION: SIEE's promise appears limited to detecting loci with antagonistic effects. Further studies are needed to evaluate SIEE's full advantage.

Three Approaches to Modeling Gene-Environment Interactions in Longitudinal Family Data: Gene-Smoking Interactions in Blood Pressure.

Blood pressure (BP) has been shown to be substantially heritable, yet identified genetic variants explain only a small fraction of the heritability. Gene-smoking interactions have detected novel BP loci in cross-sectional family data. Longitudinal family data are available and have additional promise to identify BP loci. However, this type of data presents unique analysis challenges. Although several methods for analyzing longitudinal family data are available, which method is the most appropriate and under what conditions has not been fully studied. Using data from three clinic visits from the Framingham Heart Study, we performed association analysis accounting for gene-smoking interactions in BP at 31,203 markers on chromosome 22. We evaluated three different modeling frameworks: generalized estimating equations (GEE), hierarchical linear modeling, and pedigree-based mixed modeling. The three models performed somewhat comparably, with multiple overlaps in the most strongly associated loci from each model. Loci with the greatest significance were more strongly supported in the longitudinal analyses than in any of the component single-visit analyses. The pedigree-based mixed model was more conservative, with less inflation in the variant main effect and greater deflation in the gene-smoking interactions. The GEE, but not the other two models, resulted in substantial inflation in the tail of the distribution when variants with minor allele frequency <1% were included in the analysis. The choice of analysis method should depend on the model and the structure and complexity of the familial and longitudinal data.

An Empirical Comparison of Joint and Stratified Frameworks for Studying G × E Interactions: Systolic Blood Pressure and Smoking in the CHARGE Gene-Lifestyle Interactions Working Group.

Studying gene-environment (G × E) interactions is important, as they extend our knowledge of the genetic architecture of complex traits and may help to identify novel variants not detected via analysis of main effects alone. The main statistical framework for studying G × E interactions uses a single regression model that includes both the genetic main and G × E interaction effects (the "joint" framework). The alternative "stratified" framework combines results from genetic main-effect analyses carried out separately within the exposed and unexposed groups. Although there have been several investigations using theory and simulation, an empirical comparison of the two frameworks is lacking. Here, we compare the two frameworks using results from genome-wide association studies of systolic blood pressure for 3.2 million low frequency and 6.5 million common variants across 20 cohorts of European ancestry, comprising 79,731 individuals. Our cohorts have sample sizes ranging from 456 to 22,983 and include both family-based and population-based samples. In cohort-specific analyses, the two frameworks provided similar inference for population-based cohorts. The agreement was reduced for family-based cohorts. In meta-analyses, agreement between the two frameworks was less than that observed in cohort-specific analyses, despite the increased sample size. In meta-analyses, agreement depended on (1) the minor allele frequency, (2) inclusion of family-based cohorts in meta-analysis, and (3) filtering scheme. The stratified framework appears to approximate the joint framework well only for common variants in population-based cohorts. We conclude that the joint framework is the preferred approach and should be used to control false positives when dealing with low-frequency variants and/or family-based cohorts.

The association of the vanin-1 N131S variant with blood pressure is mediated by endoplasmic reticulum-associated degradation and loss of function.

High blood pressure (BP) is the most common cardiovascular risk factor worldwide and a major contributor to heart disease and stroke. We previously discovered a BP-associated missense SNP (single nucleotide polymorphism)-rs2272996-in the gene encoding vanin-1, a glycosylphosphatidylinositol (GPI)-anchored membrane pantetheinase. In the present study, we first replicated the association of rs2272996 and BP traits with a total sample size of nearly 30,000 individuals from the Continental Origins and Genetic Epidemiology Network (COGENT) of African Americans (P=0.01). This association was further validated using patient plasma samples; we observed that the N131S mutation is associated with significantly lower plasma vanin-1 protein levels. We observed that the N131S vanin-1 is subjected to rapid endoplasmic reticulum-associated degradation (ERAD) as the underlying mechanism for its reduction. Using HEK293 cells stably expressing vanin-1 variants, we showed that N131S vanin-1 was degraded significantly faster than wild type (WT) vanin-1. Consequently, there were only minimal quantities of variant vanin-1 present on the plasma membrane and greatly reduced pantetheinase activity. Application of MG-132, a proteasome inhibitor, resulted in accumulation of ubiquitinated variant protein. A further experiment demonstrated that atenolol and diltiazem, two current drugs for treating hypertension, reduce the vanin-1 protein level. Our study provides strong biological evidence for the association of the identified SNP with BP and suggests that vanin-1 misfolding and degradation are the underlying molecular mechanism.