Returning actionable genomic results in a research biobank: Analytic validity, clinical implementation, and resource utilization.

Over 100 million research participants around the world have had research array-based genotyping (GT) or genome sequencing (GS), but only a small fraction of these have been offered return of actionable genomic findings (gRoR). Between 2017 and 2021, we analyzed genomic results from 36,417 participants in the Mass General Brigham Biobank and offered to confirm and return pathogenic and likely pathogenic variants (PLPVs) in 59 genes. Variant verification prior to participant recontact revealed that GT falsely identified PLPVs in 44.9% of samples, and GT failed to identify 72.0% of PLPVs detected in a subset of samples that were also sequenced. GT and GS detected verified PLPVs in 1% and 2.5% of the cohort, respectively. Of 256 participants who were alerted that they carried actionable PLPVs, 37.5% actively or passively declined further disclosure. 76.3% of those carrying PLPVs were unaware that they were carrying the variant, and over half of those met published professional criteria for genetic testing but had never been tested. This gRoR protocol cost approximately $129,000 USD per year in laboratory testing and research staff support, representing $14 per participant whose DNA was analyzed or $3,224 per participant in whom a PLPV was confirmed and disclosed. These data provide logistical details around gRoR that could help other investigators planning to return genomic results.

Association of GLP-1 Receptor Agonists with Chronic Obstructive Pulmonary Disease Exacerbations among Patients with Type 2 Diabetes.

Rationale: Patients with chronic obstructive pulmonary disease (COPD) and type 2 diabetes (T2D) have worse clinical outcomes compared with patients without metabolic dysregulation. GLP-1 (glucagon-like peptide 1) receptor agonists (GLP-1RAs) reduce asthma exacerbation risk and improve FVC in patients with COPD. Objectives: To determine whether GLP-1RA use is associated with reduced COPD exacerbation rates, and severe and moderate exacerbation risk, compared with other T2D therapies. Methods: A retrospective, observational, electronic health records-based study was conducted using an active comparator, new-user design of 1,642 patients with COPD in a U.S. health system from 2012 to 2022. The COPD cohort was identified using a previously validated machine learning algorithm that includes a natural language processing tool. Exposures were defined as prescriptions for GLP-1RAs (reference group), DPP-4 (dipeptidyl peptidase 4) inhibitors (DPP-4is), SGLT2 (sodium-glucose cotransporter 2) inhibitors, or sulfonylureas. Measurements and Main Results: Unadjusted COPD exacerbation counts were lower in GLP-1RA users. Adjusted exacerbation rates were significantly higher in DPP-4i (incidence rate ratio, 1.48 [95% confidence interval, 1.08-2.04]; P = 0.02) and sulfonylurea (incidence rate ratio, 2.09 [95% confidence interval, 1.62-2.69]; P < 0.0001) users compared with GLP-1RA users. GLP-1RA use was also associated with significantly reduced risk of severe exacerbations compared with DPP-4i and sulfonylurea use, and of moderate exacerbations compared with sulfonylurea use. After adjustment for clinical covariates, moderate exacerbation risk was also lower in GLP-1RA users compared with DPP-4i users. No statistically significant difference in exacerbation outcomes was seen between GLP-1RA and SGLT2 inhibitor users. Conclusions: Prospective studies of COPD exacerbations in patients with comorbid T2D are warranted. Additional research may elucidate the mechanisms underlying these observed associations with T2D medications.

Returning integrated genomic risk and clinical recommendations: The eMERGE study.

PURPOSE: Assessing the risk of common, complex diseases requires consideration of clinical risk factors as well as monogenic and polygenic risks, which in turn may be reflected in family history. Returning risks to individuals and providers may influence preventive care or use of prophylactic therapies for those individuals at high genetic risk. METHODS: To enable integrated genetic risk assessment, the eMERGE (electronic MEdical Records and GEnomics) network is enrolling 25,000 diverse individuals in a prospective cohort study across 10 sites. The network developed methods to return cross-ancestry polygenic risk scores, monogenic risks, family history, and clinical risk assessments via a genome-informed risk assessment (GIRA) report and will assess uptake of care recommendations after return of results. RESULTS: GIRAs include summary care recommendations for 11 conditions, education pages, and clinical laboratory reports. The return of high-risk GIRA to individuals and providers includes guidelines for care and lifestyle recommendations. Assembling the GIRA required infrastructure and workflows for ingesting and presenting content from multiple sources. Recruitment began in February 2022. CONCLUSION: Return of a novel report for communicating monogenic, polygenic, and family history-based risk factors will inform the benefits of integrated genetic risk assessment for routine health care.

Impact of selection bias on polygenic risk score estimates in healthcare settings.

BACKGROUND: Hospital-based biobanks are being increasingly considered as a resource for translating polygenic risk scores (PRS) into clinical practice. However, since these biobanks originate from patient populations, there is a possibility of bias in polygenic risk estimation due to overrepresentation of patients with higher frequency of healthcare interactions. METHODS: PRS for schizophrenia, bipolar disorder, and depression were calculated using summary statistics from the largest available genomic studies for a sample of 24 153 European ancestry participants in the Mass General Brigham (MGB) Biobank. To correct for selection bias, we fitted logistic regression models with inverse probability (IP) weights, which were estimated using 1839 sociodemographic, clinical, and healthcare utilization features extracted from electronic health records of 1 546 440 non-Hispanic White patients eligible to participate in the Biobank study at their first visit to the MGB-affiliated hospitals. RESULTS: Case prevalence of bipolar disorder among participants in the top decile of bipolar disorder PRS was 10.0% (95% CI 8.8-11.2%) in the unweighted analysis but only 6.2% (5.0-7.5%) when selection bias was accounted for using IP weights. Similarly, case prevalence of depression among those in the top decile of depression PRS was reduced from 33.5% (31.7-35.4%) to 28.9% (25.8-31.9%) after IP weighting. CONCLUSIONS: Non-random selection of participants into volunteer biobanks may induce clinically relevant selection bias that could impact implementation of PRS in research and clinical settings. As efforts to integrate PRS in medical practice expand, recognition and mitigation of these biases should be considered and may need to be optimized in a context-specific manner.

Pathologically expanded peripheral T helper cell subset drives B cells in rheumatoid arthritis.

CD4+ T cells are central mediators of autoimmune pathology; however, defining their key effector functions in specific autoimmune diseases remains challenging. Pathogenic CD4+ T cells within affected tissues may be identified by expression of markers of recent activation. Here we use mass cytometry to analyse activated T cells in joint tissue from patients with rheumatoid arthritis, a chronic immune-mediated arthritis that affects up to 1% of the population. This approach revealed a markedly expanded population of PD-1hiCXCR5-CD4+ T cells in synovium of patients with rheumatoid arthritis. However, these cells are not exhausted, despite high PD-1 expression. Rather, using multidimensional cytometry, transcriptomics, and functional assays, we define a population of PD-1hiCXCR5- 'peripheral helper' T (TPH) cells that express factors enabling B-cell help, including IL-21, CXCL13, ICOS, and MAF. Like PD-1hiCXCR5+ T follicular helper cells, TPH cells induce plasma cell differentiation in vitro through IL-21 secretion and SLAMF5 interaction (refs 3, 4). However, global transcriptomics highlight differences between TPH cells and T follicular helper cells, including altered expression of BCL6 and BLIMP1 and unique expression of chemokine receptors that direct migration to inflamed sites, such as CCR2, CX3CR1, and CCR5, in TPH cells. TPH cells appear to be uniquely poised to promote B-cell responses and antibody production within pathologically inflamed non-lymphoid tissues.

Development of a Definition of Postacute Sequelae of SARS-CoV-2 Infection.

Importance: SARS-CoV-2 infection is associated with persistent, relapsing, or new symptoms or other health effects occurring after acute infection, termed postacute sequelae of SARS-CoV-2 infection (PASC), also known as long COVID. Characterizing PASC requires analysis of prospectively and uniformly collected data from diverse uninfected and infected individuals. Objective: To develop a definition of PASC using self-reported symptoms and describe PASC frequencies across cohorts, vaccination status, and number of infections. Design, Setting, and Participants: Prospective observational cohort study of adults with and without SARS-CoV-2 infection at 85 enrolling sites (hospitals, health centers, community organizations) located in 33 states plus Washington, DC, and Puerto Rico. Participants who were enrolled in the RECOVER adult cohort before April 10, 2023, completed a symptom survey 6 months or more after acute symptom onset or test date. Selection included population-based, volunteer, and convenience sampling. Exposure: SARS-CoV-2 infection. Main Outcomes and Measures: PASC and 44 participant-reported symptoms (with severity thresholds). Results: A total of 9764 participants (89% SARS-CoV-2 infected; 71% female; 16% Hispanic/Latino; 15% non-Hispanic Black; median age, 47 years [IQR, 35-60]) met selection criteria. Adjusted odds ratios were 1.5 or greater (infected vs uninfected participants) for 37 symptoms. Symptoms contributing to PASC score included postexertional malaise, fatigue, brain fog, dizziness, gastrointestinal symptoms, palpitations, changes in sexual desire or capacity, loss of or change in smell or taste, thirst, chronic cough, chest pain, and abnormal movements. Among 2231 participants first infected on or after December 1, 2021, and enrolled within 30 days of infection, 224 (10% [95% CI, 8.8%-11%]) were PASC positive at 6 months. Conclusions and Relevance: A definition of PASC was developed based on symptoms in a prospective cohort study. As a first step to providing a framework for other investigations, iterative refinement that further incorporates other clinical features is needed to support actionable definitions of PASC.

Plasma metabolite profiles related to plant-based diets and the risk of type 2 diabetes.

AIMS/HYPOTHESIS: Plant-based diets, especially when rich in healthy plant foods, have been associated with a lower risk of type 2 diabetes. However, whether plasma metabolite profiles related to plant-based diets reflect this association was unknown. The aim of this study was to identify the plasma metabolite profiles related to plant-based diets, and to evaluate the associations between the identified metabolite profiles and the risk of type 2 diabetes. METHODS: Within three prospective cohorts (Nurses' Health Study, Nurses' Health Study II and Health Professionals Follow-up Study), we measured plasma metabolites from 10,684 participants using high-throughput LC MS. Adherence to plant-based diets was assessed by three indices derived from the food frequency questionnaire: an overall Plant-based Diet Index (PDI), a Healthy Plant-based Diet Index (hPDI), and an Unhealthy Plant-based Diet Index (uPDI). Multi-metabolite profiles related to plant-based diet were identified using elastic net regression with a training/testing approach. The prospective associations between metabolite profiles and incident type 2 diabetes were evaluated using multivariable Cox proportional hazards regression. Metabolites potentially mediating the association between plant-based diets and type 2 diabetes risk were further identified. RESULTS: We identified multi-metabolite profiles comprising 55 metabolites for PDI, 93 metabolites for hPDI and 75 metabolites for uPDI. Metabolite profile scores based on the identified metabolite profiles were correlated with the corresponding diet index (Pearson r = 0.33-0.35 for PDI, 0.41-0.45 for hPDI, and 0.37-0.38 for uPDI, all p<0.001). Metabolite profile scores of PDI (HR per 1 SD higher = 0.81 [95% CI 0.75, 0.88]) and hPDI (HR per 1 SD higher = 0.77 [95% CI 0.71, 0.84]) showed an inverse association with incident type 2 diabetes, whereas the metabolite profile score for uPDI was not associated with the risk. Mutual adjustment for metabolites selected in the metabolite profiles, including trigonelline, hippurate, isoleucine and a subset of triacylglycerols, attenuated the associations of diet indices PDI and hPDI with lower type 2 diabetes risk. The explainable proportion of PDI/hPDI-related diabetes risk by these metabolites ranged between 8.5% and 37.2% (all p<0.05). CONCLUSIONS/INTERPRETATION: Plasma metabolite profiles related to plant-based diets, especially a healthy plant-based diet, were associated with a lower risk of type 2 diabetes among a generally healthy population. Our findings support the beneficial role of healthy plant-based diets in diabetes prevention and provide new insights for future investigation.

Antibodies to Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) in All of Us Research Program Participants, 2 January to 18 March 2020.

BACKGROUND: With limited severe acute respiratory syndrome coronavirus (SARS-CoV-2) testing capacity in the United States at the start of the epidemic (January-March 2020), testing was focused on symptomatic patients with a travel history throughout February, obscuring the picture of SARS-CoV-2 seeding and community transmission. We sought to identify individuals with SARS-CoV-2 antibodies in the early weeks of the US epidemic. METHODS: All of Us study participants in all 50 US states provided blood specimens during study visits from 2 January to 18 March 2020. Participants were considered seropositive if they tested positive for SARS-CoV-2 immunoglobulin G (IgG) antibodies with the Abbott Architect SARS-CoV-2 IgG enzyme-linked immunosorbent assay (ELISA) and the EUROIMMUN SARS-CoV-2 ELISA in a sequential testing algorithm. The sensitivity and specificity of these ELISAs and the net sensitivity and specificity of the sequential testing algorithm were estimated, along with 95% confidence intervals (CIs). RESULTS: The estimated sensitivities of the Abbott and EUROIMMUN assays were 100% (107 of 107 [95% CI: 96.6%-100%]) and 90.7% (97 of 107 [83.5%-95.4%]), respectively, and the estimated specificities were 99.5% (995 of 1000 [98.8%-99.8%]) and 99.7% (997 of 1000 [99.1%-99.9%]), respectively. The net sensitivity and specificity of our sequential testing algorithm were 90.7% (97 of 107 [95% CI: 83.5%-95.4%]) and 100.0% (1000 of 1000 [99.6%-100%]), respectively. Of the 24 079 study participants with blood specimens from 2 January to 18 March 2020, 9 were seropositive, 7 before the first confirmed case in the states of Illinois, Massachusetts, Wisconsin, Pennsylvania, and Mississippi. CONCLUSIONS: Our findings identified SARS-CoV-2 infections weeks before the first recognized cases in 5 US states.

Psychiatric manifestations of rare variation in medically actionable genes: a PheWAS approach.

BACKGROUND: As genomic sequencing moves closer to clinical implementation, there has been an increasing acceptance of returning incidental findings to research participants and patients for mutations in highly penetrant, medically actionable genes. A curated list of genes has been recommended by the American College of Medical Genetics and Genomics (ACMG) for return of incidental findings. However, the pleiotropic effects of these genes are not fully known. Such effects could complicate genetic counseling when returning incidental findings. In particular, there has been no systematic evaluation of psychiatric manifestations associated with rare variation in these genes. RESULTS: Here, we leveraged a targeted sequence panel and real-world electronic health records from the eMERGE network to assess the burden of rare variation in the ACMG-56 genes and two psychiatric-associated genes (CACNA1C  and TCF4) across common mental health conditions in 15,181 individuals of European descent. As a positive control, we showed that this approach replicated the established association between rare mutations in LDLR and hypercholesterolemia with no visible inflation from population stratification. However, we did not identify any genes significantly enriched with rare deleterious variants that confer risk for common psychiatric disorders after correction for multiple testing. Suggestive associations were observed between depression and rare coding variation in PTEN (P = 1.5 × 10-4), LDLR (P = 3.6 × 10-4), and CACNA1S (P = 5.8 × 10-4). We also observed nominal associations between rare variants in KCNQ1 and substance use disorders (P = 2.4 × 10-4), and APOB and tobacco use disorder (P = 1.1 × 10-3). CONCLUSIONS: Our results do not support an association between psychiatric disorders and incidental findings in medically actionable gene mutations, but power was limited with the available sample sizes. Given the phenotypic and genetic complexity of psychiatric phenotypes, future work will require a much larger sequencing dataset to determine whether incidental findings in these genes have implications for risk of psychopathology.

Longitudinal SARS-CoV-2 Nucleocapsid Antibody Kinetics, Seroreversion, and Implications for Seroepidemiologic Studies.

Given widespread use of spike antibody in generating coronavirus disease vaccines, SARS-CoV-2 nucleocapsid antibodies are increasingly used to indicate previous infection in serologic surveys. However, longitudinal kinetics and seroreversion are poorly defined. We found substantial seroreversion of nucleocapsid total immunoglobulin, underscoring the need to account for seroreversion in seroepidemiologic studies.

Interaction of obesity polygenic score with lifestyle risk factors in an electronic health record biobank.

BACKGROUND: Genetic and lifestyle factors have considerable effects on obesity and related diseases, yet their effects in a clinical cohort are unknown. This study in a patient biobank examined associations of a BMI polygenic risk score (PRS), and its interactions with lifestyle risk factors, with clinically measured BMI and clinical phenotypes. METHODS: The Mass General Brigham (MGB) Biobank is a hospital-based cohort with electronic health record, genetic, and lifestyle data. A PRS for obesity was generated using 97 genetic variants for BMI. An obesity lifestyle risk index using survey responses to obesogenic lifestyle risk factors (alcohol, education, exercise, sleep, smoking, and shift work) was used to dichotomize the cohort into high and low obesogenic index based on the population median. Height and weight were measured at a clinical visit. Multivariable linear cross-sectional associations of the PRS with BMI and interactions with the obesity lifestyle risk index were conducted. In phenome-wide association analyses (PheWAS), similar logistic models were conducted for 675 disease outcomes derived from billing codes. RESULTS: Thirty-three thousand five hundred eleven patients were analyzed (53.1% female; age 60.0 years; BMI 28.3 kg/m2), of which 17,040 completed the lifestyle survey (57.5% female; age: 60.2; BMI: 28.1 (6.2) kg/m2). Each standard deviation increment in the PRS was associated with 0.83 kg/m2 unit increase in BMI (95% confidence interval (CI) =0.76, 0.90). There was an interaction between the obesity PRS and obesity lifestyle risk index on BMI. The difference in BMI between those with a high and low obesogenic index was 3.18 kg/m2 in patients in the highest decile of PRS, whereas that difference was only 1.55 kg/m2 in patients in the lowest decile of PRS. In PheWAS, the obesity PRS was associated with 40 diseases spanning endocrine/metabolic, circulatory, and 8 other disease groups. No interactions were evident between the PRS and the index on disease outcomes. CONCLUSIONS: In this hospital-based clinical biobank, obesity risk conferred by common genetic variants was associated with elevated BMI and this risk was attenuated by a healthier patient lifestyle. Continued consideration of the role of lifestyle in the context of genetic predisposition in healthcare settings is necessary to quantify the extent to which modifiable lifestyle risk factors may moderate genetic predisposition and inform clinical action to achieve personalized medicine.

Reanalysis of eMERGE phase III sequence variants in 10,500 participants and infrastructure to support the automated return of knowledge updates.

PURPOSE: The clinical genomics knowledgebase is dynamic with variant classifications changing as newly identified cases, additional population data, and other evidence become available. This is a challenge for the clinical laboratory because of limited resource availability for variant reassessment. METHODS: Throughout the Electronic Medical Records and Genomics phase III program, clinical sites associated with the Mass General Brigham/Broad sequencing center received automated, real-time notifications when reported variants were reclassified. In this study, we summarized the nature of these reclassifications and described the proactive reassessment framework we used for the Electronic Medical Records and Genomics program data set to identify variants most likely to undergo reclassification. RESULTS: Reanalysis of 1855 variants led to the reclassification of 2% (n = 45) of variants, affecting 0.6% (n = 67) of participants. Of these reclassifications, 78% (n = 35) were high-impact changes affecting reportability, with 8 variants downgraded from likely pathogenic/pathogenic to variants of uncertain significance (VUS) and 27 variants upgraded from VUS to likely pathogenic/pathogenic. Most upgraded variants (67%) were initially classified as VUS-Favor Pathogenic, highlighting the benefit of VUS subcategorization. The most common reason for reclassification was new published case data and/or functional evidence. CONCLUSION: Our results highlight the importance of periodic sequence variant reevaluation and the need for automated approaches to advance routine implementation of variant reevaluations in clinical practice.

Long-term weight changes and risk of rheumatoid arthritis among women in a prospective cohort: a marginal structural model approach.

OBJECTIVE: To examine the association of long-term weight change with RA risk in a large prospective cohort study. METHODS: The Nurses' Health Study II started in 1989 (baseline); after exclusions, we studied 108 505 women 25-42 years old without RA. Incident RA was reported by participants and confirmed by medical record review. Body weight was reported biennially through 2015. We investigated two time-varying exposures: weight changes from baseline and from age 18; change was divided into five categories. We used a marginal structural model approach to account for time-varying weight change and covariates. RESULTS: Over 2 583 266 person-years, with a median follow-up time of 25.3 years, 541 women developed RA. Compared with women with stable weight from baseline, weight change was significantly associated with increased RA risk [weight gain 2-<10 kg: RR = 1.98 (95% CI 1.38, 2.85); 10-<20 kg: RR = 3.28 (95% CI 2.20, 4.89); ≥20 kg: RR = 3.81 (95% CI 2.39, 6.07); and weight loss >2 kg: RR = 2.05 (95% CI 1.28, 3.28)]. Weight gain of 10 kg or more from age 18 compared with stable weight was also associated with increased RA risk [10-< 20 kg: RR = 2.12 (95% CI 1.37, 3.27), ≥20 kg: RR = 2.31 (95% CI 1.50, 3.56)]. Consistent findings were observed for seropositive and seronegative RA. CONCLUSION: Long-term weight gain was strongly associated with increased RA risk in women, with weight gain of ≥20 kg associated with more than a three-fold increased RA risk. Maintenance of healthy weight may be a strategy to prevent or delay RA.

Asthma Exacerbations in Patients with Type 2 Diabetes and Asthma on Glucagon-like Peptide-1 Receptor Agonists.

Rationale: GLP-1R (glucagon-like peptide-1 receptor) agonists are approved to treat type 2 diabetes mellitus and obesity. GLP-1R agonists reduce airway inflammation and hyperresponsiveness in preclinical models.Objectives: To compare rates of asthma exacerbations and symptoms between adults with type 2 diabetes and asthma prescribed GLP-1R agonists and those prescribed SGLT-2 (sodium-glucose cotransporter-2) inhibitors, DPP-4 (dipeptidyl peptidase-4) inhibitors, sulfonylureas, or basal insulin for diabetes treatment intensification.Methods: This study was an electronic health records-based new-user, active-comparator, retrospective cohort study of patients with type 2 diabetes and asthma newly prescribed GLP-1R agonists or comparator drugs at an academic healthcare system from January 2000 to March 2018. The primary outcome was asthma exacerbations; the secondary outcome was encounters for asthma symptoms. Propensity scores were calculated for GLP-1R agonist and non-GLP-1R agonist use. Zero-inflated Poisson regression models included adjustment for multiple covariates.Measurements and Main Results: Patients initiating GLP-1R agonists (n = 448), SGLT-2 inhibitors (n = 112), DPP-4 inhibitors (n = 435), sulfonylureas (n = 2,253), or basal insulin (n = 2,692) were identified. At 6 months, asthma exacerbation counts were lower in persons initiating GLP-1R agonists (reference) compared with SGLT-2 inhibitors (incidence rate ratio [IRR], 2.98; 95% confidence interval [CI], 1.30-6.80), DPP-4 inhibitors (IRR, 2.45; 95% CI, 1.54-3.89), sulfonylureas (IRR, 1.83; 95% CI, 1.20-2.77), and basal insulin (IRR, 2.58; 95% CI, 1.72-3.88). Healthcare encounters for asthma symptoms were also lower among GLP-1R agonist users.Conclusions: Adult patients with asthma prescribed GLP-1R agonists for type 2 diabetes had lower counts of asthma exacerbations compared with other drugs initiated for treatment intensification. GLP-1R agonists may represent a novel treatment for asthma associated with metabolic dysfunction.

Medical Records-Based Genetic Studies of the Complement System.

BACKGROUND: Genetic variants in complement genes have been associated with a wide range of human disease states, but well-powered genetic association studies of complement activation have not been performed in large multiethnic cohorts. METHODS: We performed medical records-based genome-wide and phenome-wide association studies for plasma C3 and C4 levels among participants of the Electronic Medical Records and Genomics (eMERGE) network. RESULTS: In a GWAS for C3 levels in 3949 individuals, we detected two genome-wide significant loci: chr.1q31.3 (CFH locus; rs3753396-A; ß=0.20; 95% CI, 0.14 to 0.25; P=1.52x10-11) and chr.19p13.3 (C3 locus; rs11569470-G; ß=0.19; 95% CI, 0.13 to 0.24; P=1.29x10-8). These two loci explained approximately 2% of variance in C3 levels. GWAS for C4 levels involved 3998 individuals and revealed a genome-wide significant locus at chr.6p21.32 (C4 locus; rs3135353-C; ß=0.40; 95% CI, 0.34 to 0.45; P=4.58x10-35). This locus explained approximately 13% of variance in C4 levels. The multiallelic copy number variant analysis defined two structural genomic C4 variants with large effect on blood C4 levels: C4-BS (ß=-0.36; 95% CI, -0.42 to -0.30; P=2.98x10-22) and C4-AL-BS (ß=0.25; 95% CI, 0.21 to 0.29; P=8.11x10-23). Overall, C4 levels were strongly correlated with copy numbers of C4A and C4B genes. In comprehensive phenome-wide association studies involving 102,138 eMERGE participants, we cataloged a full spectrum of autoimmune, cardiometabolic, and kidney diseases genetically related to systemic complement activation. CONCLUSIONS: We discovered genetic determinants of plasma C3 and C4 levels using eMERGE genomic data linked to electronic medical records. Genetic variants regulating C3 and C4 levels have large effects and multiple clinical correlations across the spectrum of complement-related diseases in humans.

Circulating blood metabolite trajectories and risk of rheumatoid arthritis among military personnel in the Department of Defense Biorepository.

OBJECTIVES: We sought to identify metabolic changes potentially related to rheumatoid arthritis (RA) pathogenesis occurring in the blood prior to its diagnosis. METHODS: In a US military biorepository, serum samples collected at two timepoints prior to a diagnosis of RA were identified. These were matched to controls who did not develop RA by subject age, race and time between sample collections and RA diagnosis time to stored serum samples. Relative abundances of 380 metabolites were measured using liquid chromatography-tandem mass spectrometry. We determined whether pre-RA case versus control status predicted metabolite concentration differences and differences over time (trajectories) using linear mixed models, assessing for interactions between time, pre-RA status and metabolite concentrations. We separately examined pre-RA and pre-seropositive RA cases versus matched controls and adjusted for smoking. Multiple comparison adjustment set the false discovery rate to 0.05. RESULTS: 291 pre-RA cases (80.8% pre seropositive RA) were matched to 292 controls, all with two serum samples (2.7±1.6 years; 1.0±0.9 years before RA/matched date). 52.0% were women; 52.8% were White, 26.8% Black and 20.4% other race. Mean age was 31.2 (±8.1) years at earliest blood draw. Fourteen metabolites had statistically significant trajectory differences among pre-RA subjects versus controls, including sex steroids, amino acid/lipid metabolism and xenobiotics. Results were similar when limited to pre seropositive RA and after adjusting for smoking. CONCLUSIONS: In this military case-control study, metabolite concentration trajectory differences in pre-RA cases versus controls implicated steroidogenesis, lipid/amino acid metabolism and xenobiotics in RA pathogenesis. Metabolites may have potential as biomarkers and/or therapeutic targets preceding RA diagnosis.

Improving sensitivity of the connective tissue disease screening questionnaire: A comparative study of various scoring methods.

OBJECTIVES: Early detection of autoimmune rheumatic diseases is crucial given their high morbidity and mortality and short window of opportunity to improve patient outcomes. Self-administered screening questionnaires such as the connective tissue disease screening questionnaire (CSQ) have been shown to promote early detection of autoimmune rheumatic diseases. However, optimal scoring of screening questionnaires may differ with prevalence of clinical features and changes in classification criteria. We compared the performance of 3 scoring methods for the CSQ for early detection of autoimmune rheumatic diseases in a multi-ethnic Asian population. METHODS: Patients who were newly referred for evaluation of possible autoimmune rheumatic diseases were invited to answer the cross-culturally adapted CSQ. Detection of autoimmune rheumatic diseases using 1) the original CSQ scoring, 2) a modified CSQ scoring and 3) a scoring based on current classification criteria, were compared to classification of autoimmune rheumatic diseases by classification criteria. RESULTS: Of 819 participants, 85 were classified as having autoimmune rheumatic diseases screened for by the adapted CSQ. The original CSQ scoring yielded relatively lower sensitivities in detecting both any and individual autoimmune rheumatic diseases (67% and 20-57%, respectively) compared to the modified CSQ scoring (81% and 60-73%, respectively) and the scoring based on current classification criteria (89% and 50-88%, respectively). CONCLUSION: The adapted CSQ with the classification criteria-based scoring achieved relatively high sensitivities in detecting autoimmune rheumatic diseases, suggesting this could be employed as the first step in population screening.

Adverse Effects of Low-Dose Methotrexate: A Randomized Trial.

Background: Low-dose methotrexate (LD-MTX) is the most commonly used drug for systemic rheumatic diseases worldwide and is the recommended first-line agent for rheumatoid arthritis. Despite extensive clinical use for more than 30 years, few data on adverse event (AE) rates derive from randomized, placebo-controlled trials, where both causality and magnitude of risk can be inferred. Objective: To investigate AE rates, risk, and risk differences comparing LD-MTX versus placebo. Design: Prespecified secondary analyses of a double-blind, placebo-controlled, randomized trial. (ClinicalTrials.gov: NCT01594333). Setting: North America. Participants: Adults with known cardiovascular disease and diabetes or metabolic syndrome. Intervention: Random allocation to LD-MTX (≤20 mg/wk) or placebo. All participants received folic acid, 1 mg/d, 6 days per week. Measurements: Risks for specific AEs of interest, as well as for all AEs, were compared across treatment groups after blinded adjudication. Results: After an active run-in period, 6158 patients were enrolled and 4786 randomly assigned to a group; median follow-up was 23 months and median dosage 15 mg/wk. Among the randomly assigned participants, 81.2% were male, median age was 65.7 years, and median body mass index was 31.5 kg/m2. Of 2391 participants assigned to LD-MTX, 2080 (87.0%) had an AE of interest, compared with 1951 of 2395 (81.5%) assigned to placebo (hazard ratio [HR], 1.17 [95% CI, 1.10 to 1.25]). The relative hazards of gastrointestinal (HR, 1.91 [CI, 1.75 to 2.10]), pulmonary (HR, 1.52 [CI, 1.16 to 1.98]), infectious (HR, 1.15 [CI, 1.01 to 1.30]), and hematologic (HR, 1.15 [CI, 1.07 to 1.23]) AEs were elevated for LD-MTX versus placebo. With the exception of increased risk for skin cancer (HR, 2.05 [CI, 1.28 to 3.28]), the treatment groups did not differ in risk for other cancer or mucocutaneous, neuropsychiatric, or musculoskeletal AEs. Renal AEs were reduced in the LD-MTX group (HR, 0.85 [CI, 0.78 to 0.93]). Limitation: The trial was done in patients without rheumatic disease who tolerated LD-MTX during an active run-in period. Conclusion: Use of LD-MTX was associated with small to moderate elevations in risks for skin cancer and gastrointestinal, infectious, pulmonary, and hematologic AEs, whereas renal AEs were decreased. Primary Funding Source: National Institutes of Health.

Geographic Variation in Obesity at the State Level in the All of Us Research Program.

INTRODUCTION: National obesity prevention strategies may benefit from precision health approaches involving diverse participants in population health studies. We used cohort data from the National Institutes of Health All of Us Research Program (All of Us) Researcher Workbench to estimate population-level obesity prevalence. METHODS: To estimate state-level obesity prevalence we used data from physical measurements made during All of Us enrollment visits and data from participant electronic health records (EHRs) where available. Prevalence estimates were calculated and mapped by state for 2 categories of body mass index (BMI) (kg/m2): obesity (BMI >30) and severe obesity (BMI >35). We calculated and mapped prevalence by state, excluding states with fewer than 100 All of Us participants. RESULTS: Data on height and weight were available for 244,504 All of Us participants from 33 states, and corresponding EHR data were available for 88,840 of these participants. The median and IQR of BMI taken from physical measurements data was 28.4 (24.4- 33.7) and 28.5 (24.5-33.6) from EHR data, where available. Overall obesity prevalence based on physical measurements data was 41.5% (95% CI, 41.3%-41.7%); prevalence of severe obesity was 20.7% (95% CI, 20.6-20.9), with large geographic variations observed across states. Prevalence estimates from states with greater numbers of All of Us participants were more similar to national population-based estimates than states with fewer participants. CONCLUSION: All of Us participants had a high prevalence of obesity, with state-level geographic variation mirroring national trends. The diversity among All of Us participants may support future investigations on obesity prevention and treatment in diverse populations.

The Mediterranean diet, plasma metabolome, and cardiovascular disease risk.

AIMS: To investigate whether metabolic signature composed of multiple plasma metabolites can be used to characterize adherence and metabolic response to the Mediterranean diet and whether such a metabolic signature is associated with cardiovascular disease (CVD) risk. METHODS AND RESULTS: Our primary study cohort included 1859 participants from the Spanish PREDIMED trial, and validation cohorts included 6868 participants from the US Nurses' Health Studies I and II, and Health Professionals Follow-up Study (NHS/HPFS). Adherence to the Mediterranean diet was assessed using a validated Mediterranean Diet Adherence Screener (MEDAS), and plasma metabolome was profiled by liquid chromatography-tandem mass spectrometry. We observed substantial metabolomic variation with respect to Mediterranean diet adherence, with nearly one-third of the assayed metabolites significantly associated with MEDAS (false discovery rate < 0.05). Using elastic net regularized regressions, we identified a metabolic signature, comprised of 67 metabolites, robustly correlated with Mediterranean diet adherence in both PREDIMED and NHS/HPFS (r = 0.28-0.37 between the signature and MEDAS; P = 3 × 10-35 to 4 × 10-118). In multivariable Cox regressions, the metabolic signature showed a significant inverse association with CVD incidence after adjusting for known risk factors (PREDIMED: hazard ratio [HR] per standard deviation increment in the signature = 0.71, P < 0.001; NHS/HPFS: HR = 0.85, P = 0.001), and the association persisted after further adjustment for MEDAS scores (PREDIMED: HR = 0.73, P = 0.004; NHS/HPFS: HR = 0.85, P = 0.004). Further genome-wide association analysis revealed that the metabolic signature was significantly associated with genetic loci involved in fatty acids and amino acids metabolism. Mendelian randomization analyses showed that the genetically inferred metabolic signature was significantly associated with risk of coronary heart disease (CHD) and stroke (odds ratios per SD increment in the genetically inferred metabolic signature = 0.92 for CHD and 0.91 for stroke; P < 0.001). CONCLUSIONS: We identified a metabolic signature that robustly reflects adherence and metabolic response to a Mediterranean diet, and predicts future CVD risk independent of traditional risk factors, in Spanish and US cohorts.