High-throughput identification of functional regulatory SNPs in systemic lupus erythematosus

Genome-wide association studies implicate multiple loci in risk for systemic lupus erythematosus (SLE), but few contain exonic variants, rendering systematic identification of non-coding variants essential to decoding SLE genetics. We utilized SNP-seq and bioinformatic enrichment to interrogate 2180 single-nucleotide polymorphisms (SNPs) from 87 SLE risk loci for potential binding of transcription factors and related proteins from B cells. 52 SNPs that passed initial screening were tested by electrophoretic mobility shift and luciferase reporter assays. To validate the approach, we studied rs2297550 in detail, finding that the risk allele enhanced binding to the transcription factor Ikaros (IKZF1), thereby modulating expression of IKBKE. Correspondingly, primary cells from genotyped healthy donors bearing the risk allele expressed higher levels of the interferon / NF-κB regulator IKKϵ. Together, these findings define a set of likely functional non-coding lupus risk variants and identify a new regulatory pathway involving rs2297550, Ikaros, and IKKϵ implicated by human genetics in risk for SLE.

Deep Resequencing of the 1q22 Locus in Non-Lobar Intracerebral Hemorrhage

OBJECTIVE

Genome-wide association studies have identified 1q22 as a susceptibility locus for cerebral small vessel diseases, including non-lobar intracerebral hemorrhage (ICH) and lacunar stroke. In the present study, we performed targeted high-depth sequencing of 1q22 in ICH cases and controls to further characterize this locus and prioritize potential causal mechanisms, which remain unknown.

METHODS

A total of 95,000 base pairs spanning 1q22, including SEMA4A, SLC25A44, and PMF1/PMF1-BGLAP were sequenced in 1,055 spontaneous ICH cases (534 lobar and 521 non-lobar) and 1,078 controls. Firth regression and Rare Variant Influential Filtering Tool analysis were used to analyze common and rare variants, respectively. Chromatin interaction analyses were performed using Hi-C, chromatin immunoprecipitation followed by sequencing, and chromatin interaction analysis with paired-end tag databases. Multivariable Mendelian randomization assessed whether alterations in gene-specific expression relative to regionally co-expressed genes at 1q22 could be causally related to ICH risk.

RESULTS

Common and rare variant analyses prioritized variants in SEMA4A 5'-UTR and PMF1 intronic regions, overlapping with active promoter and enhancer regions based on ENCODE annotation. Hi-C data analysis determined that 1q22 is spatially organized in a single chromatin loop, and that the genes therein belong to the same topologically associating domain. Chromatin immunoprecipitation followed by sequencing and chromatin interaction analysis with paired-end tag data analysis highlighted the presence of long-range interactions between the SEMA4A-promoter and PMF1-enhancer regions prioritized by association testing. Multivariable Mendelian randomization analyses demonstrated that PMF1 overexpression could be causally related to non-lobar ICH risk.

INTERPRETATION

Altered promoter-enhancer interactions leading to PMF1 overexpression, potentially dysregulating polyamine catabolism, could explain demonstrated associations with non-lobar ICH risk at 1q22, offering a potential new target for prevention of ICH and cerebral small vessel disease. ANN NEUROL 2024;95:325-337.

A cis-regulatory element regulates ERAP2 expression through autoimmune disease risk SNPs

Single-nucleotide polymorphisms (SNPs) near the ERAP2 gene are associated with various autoimmune conditions, as well as protection against lethal infections. Due to high linkage disequilibrium, numerous trait-associated SNPs are correlated with ERAP2 expression; however, their functional mechanisms remain unidentified. We show by reciprocal allelic replacement that ERAP2 expression is directly controlled by the splice region variant rs2248374. However, disease-associated variants in the downstream LNPEP gene promoter are independently associated with ERAP2 expression. Allele-specific conformation capture assays revealed long-range chromatin contacts between the gene promoters of LNPEP and ERAP2 and showed that interactions were stronger in patients carrying the alleles that increase susceptibility to autoimmune diseases. Replacing the SNPs in the LNPEP promoter by reference sequences lowered ERAP2 expression. These findings show that multiple SNPs act in concert to regulate ERAP2 expression and that disease-associated variants can convert a gene promoter region into a potent enhancer of a distal gene.

A Novel Metastatic Estrogen Receptor-Expressing Breast Cancer Model with Antiestrogen Responsiveness

Most women diagnosed with breast cancer (BC) have estrogen receptor alpha-positive (ER+) disease. The current mouse models of ER+ BC often rely on exogenous estrogen to encourage metastasis, which modifies the immune system and the function of some tissues like bone. Other studies use genetically modified or immunocompromised mouse strains, which do not accurately replicate the clinical disease. To create a model of antiestrogen responsive BC with spontaneous metastasis, we developed a mouse model of 4T1.2 triple-negative (TN) breast cancer with virally transduced ER expression that metastasizes spontaneously without exogenous estrogen stimulation and is responsive to antiestrogen drugs. Our mouse model exhibited upregulated ER-responsive genes and multi-organ metastasis without exogenous estrogen administration. Additionally, we developed a second TN BC cell line, E0771/bone, to express ER, and while it expressed ER-responsive genes, it lacked spontaneous metastasis to clinically important tissues. Following antiestrogen treatment (tamoxifen, ICI 182,780, or vehicle control), 4T1.2- and E0771/bone-derived tumor volumes and weights were significantly decreased, exemplifying antiestrogen responsivity in both cell lines. This 4T1.2 tumor model, which expresses the estrogen receptor, metastasizes spontaneously, and responds to antiestrogen treatment, will allow for further investigation into the biology and potential treatment of metastasis.

Identification of a regulatory pathway governing TRAF1 via an arthritis-associated non-coding variant

TRAF1/C5 was among the first loci shown to confer risk for inflammatory arthritis in the absence of an associated coding variant, but its genetic mechanism remains undefined. Using Immunochip data from 3,939 patients with juvenile idiopathic arthritis (JIA) and 14,412 control individuals, we identified 132 plausible common non-coding variants, reduced serially by single-nucleotide polymorphism sequencing (SNP-seq), electrophoretic mobility shift, and luciferase studies to the single variant rs7034653 in the third intron of TRAF1. Genetically manipulated experimental cells and primary monocytes from genotyped donors establish that the risk G allele reduces binding of Fos-related antigen 2 (FRA2), encoded by FOSL2, resulting in reduced TRAF1 expression and enhanced tumor necrosis factor (TNF) production. Conditioning on this JIA variant eliminated attributable risk for rheumatoid arthritis, implicating a mechanism shared across the arthritis spectrum. These findings reveal that rs7034653, FRA2, and TRAF1 mediate a pathway through which a non-coding functional variant drives risk of inflammatory arthritis in children and adults.

Gene-educational attainment interactions in a multi-population genome-wide meta-analysis identify novel lipid loci

Introduction: Educational attainment, widely used in epidemiologic studies as a surrogate for socioeconomic status, is a predictor of cardiovascular health outcomes. Methods: A two-stage genome-wide meta-analysis of low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), and triglyceride (TG) levels was performed while accounting for gene-educational attainment interactions in up to 226,315 individuals from five population groups. We considered two educational attainment variables: "Some College" (yes/no, for any education beyond high school) and "Graduated College" (yes/no, for completing a 4-year college degree). Genome-wide significant (p < 5 × 10-8) and suggestive (p < 1 × 10-6) variants were identified in Stage 1 (in up to 108,784 individuals) through genome-wide analysis, and those variants were followed up in Stage 2 studies (in up to 117,531 individuals). Results: In combined analysis of Stages 1 and 2, we identified 18 novel lipid loci (nine for LDL, seven for HDL, and two for TG) by two degree-of-freedom (2 DF) joint tests of main and interaction effects. Four loci showed significant interaction with educational attainment. Two loci were significant only in cross-population analyses. Several loci include genes with known or suggested roles in adipose (FOXP1, MBOAT4, SKP2, STIM1, STX4), brain (BRI3, FILIP1, FOXP1, LINC00290, LMTK2, MBOAT4, MYO6, SENP6, SRGAP3, STIM1, TMEM167A, TMEM30A), and liver (BRI3, FOXP1) biology, highlighting the potential importance of brain-adipose-liver communication in the regulation of lipid metabolism. An investigation of the potential druggability of genes in identified loci resulted in five gene targets shown to interact with drugs approved by the Food and Drug Administration, including genes with roles in adipose and brain tissue. Discussion: Genome-wide interaction analysis of educational attainment identified novel lipid loci not previously detected by analyses limited to main genetic effects.

Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

Validation of Epigenetic Markers for the Prediction of Response to Topical Corticosteroid Treatment in Eosinophilic Esophagitis

INTRODUCTION

We previously identified 18 CpG methylation biomarkers associated with treatment response to topical corticosteroids (tCS) in eosinophilic esophagitis (EoE). In this study, in an independent cohort, we assessed the validity of these CpG sites as treatment response biomarkers.

METHODS

DNA was extracted from prospectively biobanked esophageal biopsies from patients with newly diagnosed EoE enrolled in a randomized trial of 2 tCS formulations. Histologic response was defined as <15 eosinophils per high-power field. Pretreatment DNA methylation was assayed on the Illumina Human MethylationEPIC BeadChip. Logistic regression and area under the receiver operating characteristic curve analyses, adjusting for chip, position on the chip, age, sex, and baseline eosinophil count, were computed to test for an association between DNA methylation and treatment response at the 18 previously identified CpG sites.

RESULTS

We analyzed 88 patients (58 histologic responders, 30 nonresponders), with a mean age of 38 ± 16 years, 64% male, 97% White race. Of the 18 CpG sites, 13 met quality control criteria, and 3 were associated with responder status ( P < 0.012), including sites within UNC5B (cg26152017), ITGA6 (cg01044293), and LRRC8A (cg13962589). All 3 showed evidence of reduced methylation in treatment responders, consistent with the original discovery associations. The predictive probability for nonresponse with all 3 CpG sites was strong (area under the receiver operating characteristic curve = 0.79).

DISCUSSION

We validated epigenetic biomarkers (CpG methylation sites) for the prediction of tCS response in patients with EoE in an independent population. While not all previously identified markers replicated, 3 demonstrated a relatively high predictive probability for response to treatment and hold promise for guiding tCS treatment in EoE.

Predicting chronic postsurgical pain: current evidence and a novel program to develop predictive biomarker signatures

ABSTRACT

Chronic pain affects more than 50 million Americans. Treatments remain inadequate, in large part, because the pathophysiological mechanisms underlying the development of chronic pain remain poorly understood. Pain biomarkers could potentially identify and measure biological pathways and phenotypical expressions that are altered by pain, provide insight into biological treatment targets, and help identify at-risk patients who might benefit from early intervention. Biomarkers are used to diagnose, track, and treat other diseases, but no validated clinical biomarkers exist yet for chronic pain. To address this problem, the National Institutes of Health Common Fund launched the Acute to Chronic Pain Signatures (A2CPS) program to evaluate candidate biomarkers, develop them into biosignatures, and discover novel biomarkers for chronification of pain after surgery. This article discusses candidate biomarkers identified by A2CPS for evaluation, including genomic, proteomic, metabolomic, lipidomic, neuroimaging, psychophysical, psychological, and behavioral measures. Acute to Chronic Pain Signatures will provide the most comprehensive investigation of biomarkers for the transition to chronic postsurgical pain undertaken to date. Data and analytic resources generatedby A2CPS will be shared with the scientific community in hopes that other investigators will extract valuable insights beyond A2CPS's initial findings. This article will review the identified biomarkers and rationale for including them, the current state of the science on biomarkers of the transition from acute to chronic pain, gaps in the literature, and how A2CPS will address these gaps.

Maternal health and pregnancy outcomes in autosomal dominant tubulointerstitial kidney disease

Introduction

Autosomal dominant tubulointerstitial kidney disease (ADTKD) is an increasingly recognized cause of chronic kidney disease. ADTKD pregnancy outcomes have not previously been described.

Methods

A cross-sectional survey was sent to women from ADTKD families.

Results

Information was obtained from 85 afffected women (164 term pregnancies) and 23 controls (50 pregnancies). Only 16.5% of genetically affected women knew they had ADTKD during pregnancy. Eighteen percent of ADTKD mothers had hypertension during pregnancy versus 12% in controls (p  =  0.54) and >40% in comparative studies of chronic kidney disease in pregnancy. Eleven percent of births of ADTKD mothers were <37 weeks versus 0 in controls (p < 0.0001). Cesarean section occurred in 19% of pregnancies in affected women versus 38% of unaffected individuals (p  =  0.06). Only 12% of babies required a neonatal intensive care unit stay.

Conclusions

ADTKD pregnancies had lower rates of hypertension during pregnancy versus other forms of chronic kidney disease, which may have contributed to good maternal and fetal outcomes.

Genome-wide association study identifies multiple HLA loci for sarcoidosis susceptibility

Sarcoidosis is a complex systemic disease. Our study aimed to (1) identify novel alleles associated with sarcoidosis susceptibility; (2) provide an in-depth evaluation of HLA alleles and sarcoidosis susceptibility and (3) integrate genetic and transcription data to identify risk loci that may more directly impact disease pathogenesis. We report a genome-wide association study of 1335 sarcoidosis cases and 1264 controls of European descent (EA) and investigate associated alleles in a study of African Americans (AA: 1487 cases and 1504 controls). The EA and AA cohort was recruited from multiple United States sites. HLA alleles were imputed and tested for association with sarcoidosis susceptibility. Expression quantitative locus and colocalization analysis were performed using a subset of subjects with transcriptome data. Forty-nine SNPs in the HLA region in HLA-DRA, -DRB9, -DRB5, -DQA1 and BRD2 genes were significantly associated with sarcoidosis susceptibility in EA, rs3129888 was also a risk variant for sarcoidosis in AA. Classical HLA alleles DRB1*0101, DQA1*0101 and DQB1*0501, which are highly correlated, were also associated with sarcoidosis. rs3135287 near HLA-DRA was associated with HLA-DRA expression in peripheral blood mononuclear cells and bronchoalveolar lavage from subjects and lung tissue and whole blood from GTEx. We identified six novel SNPs (out of the seven SNPs representing the 49 significant SNPs) and nine HLA alleles associated with sarcoidosis susceptibility in the largest EA population. We also replicated our findings in an AA population. Our study reiterates the potential role of antigen recognition and/or presentation HLA class II genes in sarcoidosis pathogenesis.

Social Determinants of Health and Cerebral Small Vessel Disease: Is Epigenetics a Key Mediator?

Cerebral small vessel disease is highly prevalent, particularly in marginalized communities, and its incidence is expected to increase given the aging global population. Cerebral small vessel disease contributes to risk for stroke, vascular cognitive impairment and dementia, late-life depression, and gait disorders. A growing body of evidence suggests that adverse outcomes, including cerebral small vessel disease, caused by traditional cardiovascular risk factors are at least partly mediated by epigenetic changes, some of them already beginning during fetal development. Societal and health care access inequities, summarized under the umbrella term social determinants of health, put a higher burden of cardiovascular risk factors on marginalized populations and expose them to an increased risk for adverse outcomes. Social epigenetics has begun to deliver solid evidence that social determinants of health lead to distinct epigenetic signatures that potentially mediate the biological effect of environmental exposures on cardiovascular risk factors. Here, we provide a review of the most recent advances in the epigenetics of cerebral small vessel disease risk factors and social determinants of health and call for research efforts combining insights from both fields to reach a deeper understanding of the causal pathways, ultimately facilitating discovery of new treatment targets for a disease whose burden is magnified by existing health disparities.

Identification of influential rare variants in aggregate testing using random forest importance measures

Aggregate tests of rare variants are often employed to identify associated regions compared to sequentially testing each individual variant. When an aggregate test is significant, it is of interest to identify which rare variants are "driving" the association. We recently developed the rare variant influential filtering tool (RIFT) to identify influential rare variants and showed RIFT had higher true positive rates compared to other published methods. Here we use importance measures from the standard random forest (RF) and variable importance weighted RF (vi-RF) to identify influential variants. For very rare variants (minor allele frequency [MAF] < 0.001), the vi-RF:Accuracy method had the highest median true positive rate (TPR = 0.24; interquartile range [IQR]: 0.13, 0.42) followed by the RF:Accuracy method (TPR = 0.16; IQR: 0.07, 0.33) and both were superior to RIFT (TPR = 0.05; IQR: 0.02, 0.15). Among uncommon variants (0.001 < MAF < 0.03), the RF methods had higher true positive rates than RIFT while observing comparable false positive rates. Finally, we applied the RF methods to a targeted resequencing study in idiopathic pulmonary fibrosis (IPF), in which the vi-RF approach identified eight and seven variants in TERT and FAM13A, respectively. In summary, the vi-RF provides an improved, objective approach to identifying influential variants following a significant aggregate test. We have expanded our previously developed R package RIFT to include the random forest methods.

Deep resequencing of the 1q22 locus in non-lobar intracerebral hemorrhage

Objective

Genome-wide association studies have identified 1q22 as a susceptibility locus for cerebral small vessel diseases (CSVDs), including non-lobar intracerebral hemorrhage (ICH) and lacunar stroke. In the present study we performed targeted high-depth sequencing of 1q22 in ICH cases and controls to further characterize this locus and prioritize potential causal mechanisms, which remain unknown.

Methods

95,000 base pairs spanning 1q22 , including SEMA4A, SLC25A44 and PMF1 / PMF1-BGLAP were sequenced in 1,055 spontaneous ICH cases (534 lobar and 521 non-lobar) and 1,078 controls. Firth regression and RIFT analysis were used to analyze common and rare variants, respectively. Chromatin interaction analyses were performed using Hi-C, ChIP-Seq and ChIA-PET databases. Multivariable Mendelian randomization (MVMR) assessed whether alterations in gene-specific expression relative to regionally co-expressed genes at 1q22 could be causally related to ICH risk.

Results

Common and rare variant analyses prioritized variants in SEMA4A 5'-UTR and PMF1 intronic regions, overlapping with active promoter and enhancer regions based on ENCODE annotation. Hi-C data analysis determined that 1q22 is spatially organized in a single chromatin loop and that the genes therein belong to the same Topologically Associating Domain. ChIP-Seq and ChIA-PET data analysis highlighted the presence of long-range interactions between the SEMA4A -promoter and PMF1 -enhancer regions prioritized by association testing. MVMR analyses demonstrated that PMF1 overexpression could be causally related to non-lobar ICH risk.

Interpretation

Altered promoter-enhancer interactions leading to PMF1 overexpression, potentially dysregulating polyamine catabolism, could explain demonstrated associations with non-lobar ICH risk at 1q22 , offering a potential new target for prevention of ICH and CSVD.

Metabolomic profiling of glucose homeostasis in African Americans: the Insulin Resistance Atherosclerosis Family Study (IRAS-FS)

INTRODUCTION

African Americans are at increased risk for type 2 diabetes.

OBJECTIVES

This work aimed to examine metabolomic signature of glucose homeostasis in African Americans.

METHODS

We used an untargeted liquid chromatography-mass spectrometry metabolomic approach to comprehensively profile 727 plasma metabolites among 571 African Americans from the Insulin Resistance Atherosclerosis Family Study (IRAS-FS) and investigate the associations between these metabolites and both the dynamic (S_I, insulin sensitivity; AIR, acute insulin response; DI, disposition index; and S_G, glucose effectiveness) and basal (HOMA-IR and HOMA-B) measures of glucose homeostasis using univariate and regularized regression models. We also compared the results with our previous findings in the IRAS-FS Mexican Americans.

RESULTS

We confirmed increased plasma metabolite levels of branched-chain amino acids and their metabolic derivatives, 2-aminoadipate, 2-hydroxybutyrate, glutamate, arginine and its metabolic derivatives, carbohydrate metabolites, and medium- and long-chain fatty acids were associated with insulin resistance, while increased plasma metabolite levels in the glycine, serine and threonine metabolic pathway were associated with insulin sensitivity. We also observed a differential ancestral effect of glutamate on glucose homeostasis with significantly stronger effects observed in African Americans than those previously observed in Mexican Americans.

CONCLUSION

We extended the observations that metabolites are useful biomarkers in the identification of prediabetes in individuals at risk of type 2 diabetes in African Americans. We revealed, for the first time, differential ancestral effect of certain metabolites (i.e., glutamate) on glucose homeostasis traits. Our study highlights the need for additional comprehensive metabolomic studies in well-characterized multiethnic cohorts.

A Genomic Risk Score Identifies Individuals at High Risk for Intracerebral Hemorrhage

BACKGROUND

Intracerebral hemorrhage (ICH) has an estimated heritability of 29%. We developed a genomic risk score for ICH and determined its predictive power in comparison to standard clinical risk factors.

METHODS

We combined genome-wide association data from individuals of European ancestry for ICH and related traits in a meta-genomic risk score ([metaGRS]; 2.6 million variants). We tested associations with ICH and its predictive performance in addition to clinical risk factors in a held-out validation dataset (842 cases and 796 controls). We tested associations with risk of incident ICH in the population-based UK Biobank cohort (486 784 individuals, 1526 events, median follow-up 11.3 years).

RESULTS

One SD increment in the metaGRS was significantly associated with 31% higher odds for ICH (95% CI, 1.16-1.48) in age-, sex- and clinical risk factor-adjusted models. The metaGRS identified individuals with almost 5-fold higher odds for ICH in the top score percentile (odds ratio, 4.83 [95% CI, 1.56-21.2]). Predictive models for ICH incorporating the metaGRS in addition to clinical predictors showed superior performance compared to the clinical risk factors alone (c-index, 0.695 versus 0.686). The metaGRS showed similar associations for lobar and nonlobar ICH, independent of the known APOE risk locus for lobar ICH. In the UK Biobank, the metaGRS was associated with higher risk of incident ICH (hazard ratio, 1.15 [95% CI, 1.09-1.21]). The associations were significant within both a relatively high-risk population of antithrombotic medications users, as well as among a relatively low-risk population with a good control of vascular risk factors and no use of anticoagulants.

CONCLUSIONS

We developed and validated a genomic risk score that predicts lifetime risk of ICH beyond established clinical risk factors among individuals of European ancestry. Whether implementation of the score in risk prognostication models for high-risk populations, such as patients under antithrombotic treatment, could improve clinical decision making should be explored in future studies.

Molecular pathways identified from single nucleotide polymorphisms demonstrate mechanistic differences in systemic lupus erythematosus patients of Asian and European ancestry

Systemic lupus erythematosus (SLE) is a multi-organ autoimmune disorder with a prominent genetic component. Individuals of Asian-Ancestry (AsA) disproportionately experience more severe SLE compared to individuals of European-Ancestry (EA), including increased renal involvement and tissue damage. However, the mechanisms underlying elevated severity in the AsA population remain unclear. Here, we utilized available gene expression data and genotype data based on all non-HLA SNP associations in EA and AsA SLE patients detected using the Immunochip genotyping array. We identified 2778 ancestry-specific and 327 trans-ancestry SLE-risk polymorphisms. Genetic associations were examined using connectivity mapping and gene signatures based on predicted biological pathways and were used to interrogate gene expression datasets. SLE-associated pathways in AsA patients included elevated oxidative stress, altered metabolism and mitochondrial dysfunction, whereas SLE-associated pathways in EA patients included a robust interferon response (type I and II) related to enhanced cytosolic nucleic acid sensing and signaling. An independent dataset derived from summary genome-wide association data in an AsA cohort was interrogated and identified similar molecular pathways. Finally, gene expression data from AsA SLE patients corroborated the molecular pathways predicted by SNP associations. Identifying ancestry-related molecular pathways predicted by genetic SLE risk may help to disentangle the population differences in clinical severity that impact AsA and EA individuals with SLE.

Distinct genome-wide DNA methylation and gene expression signatures in classical monocytes from African American patients with systemic sclerosis

BACKGROUND

Systemic sclerosis (SSc) is a multisystem autoimmune disorder that has an unclear etiology and disproportionately affects women and African Americans. Despite this, African Americans are dramatically underrepresented in SSc research. Additionally, monocytes show heightened activation in SSc and in African Americans relative to European Americans. In this study, we sought to investigate DNA methylation and gene expression patterns in classical monocytes in a health disparity population.

METHODS

Classical monocytes (CD14+ + CD16-) were FACS-isolated from 34 self-reported African American women. Samples from 12 SSc patients and 12 healthy controls were hybridized on MethylationEPIC BeadChip array, while RNA-seq was performed on 16 SSc patients and 18 healthy controls. Analyses were computed to identify differentially methylated CpGs (DMCs), differentially expressed genes (DEGs), and CpGs associated with changes in gene expression (eQTM analysis).

RESULTS

We observed modest DNA methylation and gene expression differences between cases and controls. The genes harboring the top DMCs, the top DEGs, as well as the top eQTM loci were enriched for metabolic processes. Genes involved in immune processes and pathways showed a weak upregulation in the transcriptomic analysis. While many genes were newly identified, several other have been previously reported as differentially methylated or expressed in different blood cells from patients with SSc, supporting for their potential dysregulation in SSc.

CONCLUSIONS

While contrasting with results found in other blood cell types in largely European-descent groups, the results of this study support that variation in DNA methylation and gene expression exists among different cell types and individuals of different genetic, clinical, social, and environmental backgrounds. This finding supports the importance of including diverse, well-characterized patients to understand the different roles of DNA methylation and gene expression variability in the dysregulation of classical monocytes in diverse populations, which might help explaining the health disparities.

Multi-omic integration reveals cell-type-specific regulatory networks of insulin resistance in distinct ancestry populations

Our knowledge of the cell-type-specific mechanisms of insulin resistance remains limited. To dissect the cell-type-specific molecular signatures of insulin resistance, we performed a multiscale gene network analysis of adipose and muscle tissues in African and European ancestry populations. In adipose tissues, a comparative analysis revealed ethnically conserved cell-type signatures and two adipocyte subtype-enriched modules with opposite insulin sensitivity responses. The modules enriched for adipose stem and progenitor cells as well as immune cells showed negative correlations with insulin sensitivity. In muscle tissues, the modules enriched for stem cells and fibro-adipogenic progenitors responded to insulin sensitivity oppositely. The adipocyte and muscle fiber-enriched modules shared cellular-respiration-related genes but had tissue-specific rearrangements of gene regulations in response to insulin sensitivity. Integration of the gene co-expression and causal networks further pinpointed key drivers of insulin resistance. Together, this study revealed the cell-type-specific transcriptomic networks and signaling maps underlying insulin resistance in major glucose-responsive tissues. A record of this paper's transparent peer review process is included in the supplemental information.

DNA methylation of the TPMT gene and azathioprine pharmacokinetics in children with very early onset inflammatory bowel disease

BACKGROUND

Thiopurine methyltransferase (TPMT) is a crucial enzyme for azathioprine biotransformation and its activity is higher in very early onset inflammatory bowel disease (VEO-IBD) patients than in adolescents with IBD (aIBD).

AIMS

The aims of this pharmacoepigenetic study were to evaluate differences in peripheral blood DNA methylation of the TPMT gene and in azathioprine pharmacokinetics in patients with VEO-IBD compared to aIBD.

METHODS

The association of age with whole genome DNA methylation profile was evaluated in a pilot group of patients and confirmed by a meta-analysis on 3 cohorts of patients available on the public functional genomics data repository. Effects of candidate CpG sites in the TPMT gene were validated in a larger cohort using pyrosequencing. TPMT activity and azathioprine metabolites (TGN) were measured in patients' erythrocytes by HPLC and associated with patients' age group and TPMT DNA methylation.

RESULTS

Whole genome DNA methylation pilot analysis, combined with the meta-analysis revealed cg22736354, located on TPMT downstream neighboring region, as the only statistically significant CpG whose methylation increases with age, resulting lower in VEO-IBD patients compared to aIBD (median 9.6% vs 12%, p = 0.029). Pyrosequencing confirmed lower cg22736354 methylation in VEO-IBD patients (median 4.0% vs 6.0%, p = 4.6 ×10^-5). No differences in TPMT promoter methylation were found. Reduced cg22736354 methylation was associated with lower TGN concentrations (rho = 0.31, p = 0.01) in patients with VEO-IBD and aIBD.

CONCLUSION

Methylation of cg22736354 in TPMT gene neighborhood is lower in patients with VEO-IBD and is associated with reduced azathioprine inactivation and increased TGN concentrations.

Comparison of different gene addition strategies to modify placental derived-mesenchymal stromal cells to produce FVIII

Introduction

Placenta-derived mesenchymal cells (PLCs) endogenously produce FVIII, which makes them ideally suited for cell-based fVIII gene delivery. We have previously reported that human PLCs can be efficiently modified with a lentiviral vector encoding a bioengineered, expression/secretion-optimized fVIII transgene (ET3) and durably produce clinically relevant levels of functionally active FVIII. The objective of the present study was to investigate whether CRISPR/Cas9 can be used to achieve location-specific insertion of a fVIII transgene into a genomic safe harbor, thereby eliminating the potential risks arising from the semi-random genomic integration inherent to lentiviral vectors. We hypothesized this approach would improve the safety of the PLC-based gene delivery platform and might also enhance the therapeutic effect by eliminating chromatin-related transgene silencing.

Methods

We used CRISPR/Cas9 to attempt to insert the bioengineered fVIII transgene "lcoET3" into the AAVS1 site of PLCs (CRISPR-lcoET3) and determined their subsequent levels of FVIII production, comparing results with this approach to those achieved using lentivector transduction (LV-lcoET3) and plasmid transfection (Plasmid-lcoET3). In addition, since liver-derived sinusoidal endothelial cells (LSECs) are the native site of FVIII production in the body, we also performed parallel studies in human (h)LSECs).

Results

PLCs and hLSECs can both be transduced (LV-lcoET3) with very high efficiency and produce high levels of biologically active FVIII. Surprisingly, both cell types were largely refractory to CRISPR/Cas9-mediated knockin of the lcoET3 fVIII transgene in the AAVS1 genome locus. However, successful insertion of an RFP reporter into this locus using an identical procedure suggests the failure to achieve knockin of the lcoET3 expression cassette at this site is likely a function of its large size. Importantly, using plasmids, alone or to introduce the CRISPR/Cas9 "machinery", resulted in dramatic upregulation of TLR 3, TLR 7, and BiP in PLCs, compromising their unique immune-inertness.

Discussion

Although we did not achieve our primary objective, our results validate the utility of both PLCs and hLSECs as cell-based delivery vehicles for a fVIII transgene, and they highlight the hurdles that remain to be overcome before primary human cells can be gene-edited with sufficient efficiency for use in cell-based gene therapy to treat HA.

Mendelian randomization and pathway analysis demonstrate shared genetic associations between lupus and coronary artery disease

Coronary artery disease (CAD) is a leading cause of death in patients with systemic lupus erythematosus (SLE). Despite clinical evidence supporting an association between SLE and CAD, pleiotropy-adjusted genetic association studies are limited and focus on only a few common risk loci. Here, we identify a net positive causal estimate of SLE-associated non-HLA SNPs on CAD by traditional Mendelian randomization (MR) approaches. Pathway analysis using SNP-to-gene mapping followed by unsupervised clustering based on protein-protein interactions (PPIs) identifies biological networks composed of positive and negative causal sets of genes. In addition, we confirm the casual effects of specific SNP-to-gene modules on CAD using only SNP mapping to each PPI-defined functional gene set as instrumental variables. This PPI-based MR approach elucidates various molecular pathways with causal implications between SLE and CAD and identifies biological pathways likely causative of both pathologies, revealing known and novel therapeutic interventions for managing CAD in SLE.

COVID-19 and systemic lupus erythematosus genetics: A balance between autoimmune disease risk and protection against infection

Genome wide association studies show there is a genetic component to severe COVID-19. We find evidence that the genome-wide genetic association signal with severe COVID-19 is correlated with that of systemic lupus erythematosus (SLE), having formally tested this using genetic correlation analysis by LD score regression. To identify the shared associated loci and gain insight into the shared genetic effects, using summary level data we performed meta-analyses, a local genetic correlation analysis and fine-mapping using stepwise regression and functional annotation. This identified multiple loci shared between the two traits, some of which exert opposing effects. The locus with most evidence of shared association is TYK2, a gene critical to the type I interferon pathway, where the local genetic correlation is negative. Another shared locus is CLEC1A, where the direction of effects is aligned, that encodes a lectin involved in cell signaling, and the anti-fungal immune response. Our analyses suggest that several loci with reciprocal effects between the two traits have a role in the defense response pathway, adding to the evidence that SLE risk alleles are protective against infection.

We observed a correlation between the genetic associations with severe COVID-19 and those with systemic lupus erythematosus (SLE, Lupus), and aimed to discover which genetic loci were shared by these diseases and what biological processes were involved. This resulted in the discovery of several genetic loci, some of which had alleles that were risk for both diseases and some of which were risk for severe COVID-19 yet protective for SLE. The locus with most evidence of shared association (TYK2) is involved in interferon production, a process that is important in response to viral infection and known to be dysregulated in SLE patients. Other shared associated loci contained genes also involved in the defense response and the immune system signaling. These results add to the growing evidence that there are alleles in the human genome that provide protection against viral infection yet are risk for autoimmune disease.

Erratum: Ancestral diversity improves discovery and fine-mapping of genetic loci for anthropometric traits-The Hispanic/Latino Anthropometry Consortium

[This corrects the article DOI: 10.1016/j.xhgg.2022.100099.].

Alzheimer's disease related single nucleotide polymorphisms and correlation with intracerebral hemorrhage incidence

Apolipoprotein E alleles have been associated with both Alzheimer's disease (AD) and intracerebral hemorrhage (ICH). In addition, ICH is associated with a markedly high risk of subsequent dementia compared to other subtypes of stroke. We sought to evaluate if other genetic markers for AD were also associated with ICH. We examined whether published AD risk single nucleotide polymorphisms (SNPs) and haplotypes were associated with ICH utilizing genome-wide association study data from 2 independent studies (genetic and environmental risk factors for hemorrhagic stroke [GERFHS] study and genetics of cerebral hemorrhage with anticoagulation [GOCHA]). Analyses included evaluation by location of ICH. GERFHS and GOCHA cohorts contained 745 ICH cases and 536 controls for analysis. The strongest association was on 1q32 near Complement receptor type 1 (CR1), where rs6701713 was associated with all ICH (P = .0074, odds ratio [OR] = 2.07) and lobar ICH (P = .0073, OR = 2.80). The 51 most significant 2-SNP haplotypes associated with lobar ICH were identified within the Clusterin (CLU) gene. We identified that variation within CR1 and CLU, previously identified risk factors for AD, and are associated with an increased risk for ICH driven primarily by lobar ICH. Previous work implicated CR1 and CLU in cerebral amyloid clearance, the innate immune system, and cellular stress response.

Relationship Between a Vitamin D Genetic Risk Score and Autoantibodies Among First-Degree Relatives of Probands With Rheumatoid Arthritis and Systemic Lupus Erythematosus

Objective

Higher 25-hydroxyvitamin D (25(OH)D) levels have been associated with reduced risk for autoimmune diseases and are influenced by vitamin D metabolism genes. We estimated genetically-determined vitamin D levels by calculating a genetic risk score (GRS) and investigated whether the vitamin D GRS was associated with the presence of autoantibodies related to rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) in those at increased risk for developing RA and SLE, respectively.

Methods

In this cross-sectional study, we selected autoantibody positive (aAb+) and autoantibody negative (aAb-) individuals from the Studies of the Etiologies of Rheumatoid Arthritis (SERA), a cohort study of first-degree relatives (FDRs) of individuals with RA (189 RA aAb+, 181 RA aAb-), and the Lupus Family Registry and Repository (LFRR), a cohort study of FDRs of individuals with SLE (157 SLE aAb+, 185 SLE aAb-). Five SNPs known to be associated with serum 25(OH)D levels were analyzed individually as well as in a GRS: rs4588 (GC), rs12785878 (NADSYN1), rs10741657 (CYP2R1), rs6538691 (AMDHD1), and rs8018720 (SEC23A).

Results

Both cohorts had similar demographic characteristics, with significantly older and a higher proportion of males in the aAb+ FDRs. The vitamin D GRS was inversely associated with RA aAb+ (OR = 0.85, 95% CI = 0.74-0.99), suggesting a possible protective factor for RA aAb positivity in FDRs of RA probands. The vitamin D GRS was not associated with SLE aAb+ in the LFRR (OR = 1.09, 95% CI = 0.94-1.27). The SEC23A SNP was associated with RA aAb+ in SERA (OR = 0.65, 95% CI = 0.43-0.99); this SNP was not associated with SLE aAb+ in LFRR (OR = 1.41, 95% CI = 0.90 - 2.19).

Conclusion

Genes associated with vitamin D levels may play a protective role in the development of RA aAbs in FDRs of RA probands, perhaps through affecting lifelong vitamin D status. The GRS and the SEC23A SNP may be of interest for future investigation in pre-clinical RA. In contrast, these results do not support a similar association in SLE FDRs, suggesting other mechanisms involved in the relationship between vitamin D and SLE aAbs not assessed in this study.

Gut dysbiosis and the clinical spectrum in anti-Ro positive mothers of children with neonatal lupus

Anti-SSA/Ro antibodies, while strongly linked to fetal cardiac injury and neonatal rash, can associate with a spectrum of disease in the mother, ranging from completely asymptomatic to overt Systemic Lupus Erythematosus (SLE) or Sjögren's Syndrome (SS). This study was initiated to test the hypothesis that the microbiome, influenced in part by genetics, contributes to disease state. The stool microbiome of healthy controls (HC) was compared to that of anti-SSA/Ro positive women whose children had neonatal lupus. At the time of sampling, these women were either asymptomatic (Asym), had minor rheumatic symptoms or signs considered as an undifferentiated autoimmune syndrome (UAS), or were diagnosed with SLE or SS. Differences in microbial relative abundances among these three groups were tested assuming an ordering in clinical severity (HC<Asym/UAS<SS/SLE) and then again without the ordinal assumption. Those taxa that showed differential relative abundances were then tested for whether the effect size differed depending on the women's HLA SLE-risk allele genotype (DRB1*03:01, DRB1*15:01, DQB1*02:01 and DQB1*06:02) or anti-SSA/Ro autoantibody levels. Multiple genera within the families Ruminococcaceae and Lachnospiraceae showed evidence of an HLA-by-genus interaction (P < .05). Four genera exhibited evidence of an interaction with anti-Ro52 IgA: Lachnoclostridium, Romboutsia, Bacteroides and Actinomyces (P < .01). In addition to documenting differences in microbial relative abundances across clinical severity of disease, these data provide a first-time demonstration that microbial differences are correlated with HLA SLE-risk alleles. Taken together, these data suggest that the clinical spectrum from benign to overt clinical autoimmunity may partially result from or trigger a complex interplay among specific microbial profiles, anti-Ro autoantibodies, and genetics.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

Social Factors, Epigenomics and Lupus in African American Women (SELA) Study: protocol for an observational mechanistic study examining the interplay of multiple individual and social factors on lupus outcomes in a health disparity population

Introduction

Despite the disproportional impact of SLE on historically marginalised communities, the individual and sociocultural factors underlying these health disparities remain elusive. We report the design and methods for a study aimed at identifying epigenetic biomarkers associated with racism and resiliency that affect gene function and thereby influence SLE in a health disparity population.

Methods and analysis

The Social Factors, Epigenomics and Lupus in African American Women (SELA) Study is a cross-sectional, case–control study. A total of 600 self-reported African American women will be invited to participate. All participants will respond to questionnaires that capture detailed sociodemographic and medical history, validated measures of racial discrimination, social support, as well as disease activity and damage for cases. Participants who wish will receive their genetic ancestry estimates and be involved in research. Blood samples are required to provide peripheral blood mononuclear cell counts, DNA and RNA. The primary goals of SELA are to identify variation in DNA methylation (DNAm) associated with self-reported exposure to racial discrimination and social support, to evaluate whether social DNAm sites affect gene expression, to identify the synergistic effects of social factors on DNAm changes on SLE and to develop a social factors-DNAm predictive model for disease outcomes. This study is conducted in cooperation with the Sea Island Families Project Citizen Advisory Committee.

Discussion and dissemination

SELA will respond to the pressing need to clarify the interplay and regulatory mechanism by which various positive and negative social exposures influence SLE. Results will be published and shared with patients and the community. Knowledge of the biological impact of social exposures on SLE, as informed by the results of this study, can be leveraged by advocacy efforts to develop psychosocial interventions that prevent or mitigate risk exposures, and services or interventions that promote positive exposures. Implementation of such interventions is paramount to the closure of the health disparities gap.

Shared genetic background between SARS-CoV-2 infection and large artery stroke

BACKGROUND AND AIMS

Increased risk of stroke, particularly large artery stroke (LAS), has been observed in patients with COVID-19. The biological processes underlying the observed higher risk are still unknown. We explored the association between stroke subtypes and COVID-19 susceptibility to understand whether biological mechanisms specific to SARS-CoV-2 uptake/infection could be leading to excess stroke risk in this population.

PATIENTS AND METHODS

We constructed a polygenic risk score (PRS) of COVID-19 susceptibility and tested its association with stroke subtypes using individual- and summary-level genetic data (SiGN, MEGASTROKE). We generated co-expression networks of genes involved in SARS-CoV-2 uptake/infection (ACE2, TMPRSS2, BEST3, ISLR2 and ADAM17) based on existing tissue expression libraries. Gene-based association testing was performed using S-PrediXcan and VEGAS2. Permutation independence tests were performed to assess SARS-CoV-2-related gene enrichment in stroke and its subtypes.

RESULTS

Our PRS demonstrated an association between COVID-19 susceptibility and LAS in SiGN (OR = 1.05 per SD increase, 95% CI: (1.00, 1.10), p = 0.04) and MEGASTROKE (β = 0.510, 95% CI: (0.242, 0.779), FDR-p = 0.0019). The SARS-CoV-2 risk-related ISLR2 co-expression gene network was significantly associated with genetic risk of LAS in aorta, tibial arteries, and multiple brain regions (P < 0.05).

CONCLUSION

Presence of genetic correlation and significant pathway enrichment suggest that increases in LAS risk reported in COVID-19 patients may be intrinsic to the viral infection, rather than a more generalized response to severe illness.

Multi-Site Observational Study to Assess Biomarkers for Susceptibility or Resilience to Chronic Pain: The Acute to Chronic Pain Signatures (A2CPS) Study Protocol

Chronic pain has become a global health problem contributing to years lived with disability and reduced quality of life. Advances in the clinical management of chronic pain have been limited due to incomplete understanding of the multiple risk factors and molecular mechanisms that contribute to the development of chronic pain. The Acute to Chronic Pain Signatures (A2CPS) Program aims to characterize the predictive nature of biomarkers (brain imaging, high-throughput molecular screening techniques, or "omics," quantitative sensory testing, patient-reported outcome assessments and functional assessments) to identify individuals who will develop chronic pain following surgical intervention. The A2CPS is a multisite observational study investigating biomarkers and collective biosignatures (a combination of several individual biomarkers) that predict susceptibility or resilience to the development of chronic pain following knee arthroplasty and thoracic surgery. This manuscript provides an overview of data collection methods and procedures designed to standardize data collection across multiple clinical sites and institutions. Pain-related biomarkers are evaluated before surgery and up to 3 months after surgery for use as predictors of patient reported outcomes 6 months after surgery. The dataset from this prospective observational study will be available for researchers internal and external to the A2CPS Consortium to advance understanding of the transition from acute to chronic postsurgical pain.

Initial antihypertensive agent effects on acute blood pressure after intracerebral haemorrhage

Introduction

Current guidelines recommend blood pressure (BP) lowering in patients after acute intracerebral haemorrhage (ICH) without guidance on initial choice of antihypertensive class. This study sought to determine if initial antihypertensive class differentially effects acute BP lowering in a large multiethnic ICH cohort.

Methods

Subjects enrolled in the Ethnic/Racial Variations in ICH study between August 2010 and August 2017 with elevated admission BP and who received labetalol, nicardipine or hydralazine monotherapy as initial antihypertensive were analysed. Primary outcomes were systolic and diastolic BP changes from baseline to first BP measurement after initial antihypertensive treatment. Secondary outcomes included haematoma expansion (HE), hospital length of stay (LOS) and modified Rankin Score (mRS) up to 12 months after ICH. Exploratory outcomes assessed effects of race/ethnicity. Linear and logistic regression analyses, adjusted for relevant covariates, were performed to determine associations of antihypertensive class with outcomes.

Results

In total, 1156 cases were used in analyses. Antihypertensive class was associated with diastolic BP change (p=0.003), but not systolic BP change (p=0.419). Initial dosing with nicardipine lowered acute diastolic BP than labetalol (least square mean difference (labetalol-nicardipine)=5.47 (2.37, 8.57), p<0.001). Initial antihypertensive class was also found to be associated with LOS (p=0.028), but not with HE (p=0.406), mortality (p=0.118), discharge disposition (p=0.083) or mRS score at discharge, 3, 6 and 12 months follow-up (p=0.262, 0.276, 0.152 and 0.36, respectively). Race/ethnicity variably affected multivariable models.

Conclusion

In this large acute ICH cohort, initial antihypertensive class was associated with acute diastolic, but not systolic, BP-lowering suggesting differential effects of antihypertensive agents.

Trial registration number

NCT01202864.

Ancestral diversity improves discovery and fine-mapping of genetic loci for anthropometric traits-The Hispanic/Latino Anthropometry Consortium

Hispanic/Latinos have been underrepresented in genome-wide association studies (GWAS) for anthropometric traits despite their notable anthropometric variability, ancestry proportions, and high burden of growth stunting and overweight/obesity. To address this knowledge gap, we analyzed densely imputed genetic data in a sample of Hispanic/Latino adults to identify and fine-map genetic variants associated with body mass index (BMI), height, and BMI-adjusted waist-to-hip ratio (WHRadjBMI). We conducted a GWAS of 18 studies/consortia as part of the Hispanic/Latino Anthropometry (HISLA) Consortium (stage 1, n = 59,771) and generalized our findings in 9 additional studies (stage 2, n = 10,538). We conducted a trans-ancestral GWAS with summary statistics from HISLA stage 1 and existing consortia of European and African ancestries. In our HISLA stage 1 + 2 analyses, we discovered one BMI locus, as well as two BMI signals and another height signal each within established anthropometric loci. In our trans-ancestral meta-analysis, we discovered three BMI loci, one height locus, and one WHRadjBMI locus. We also identified 3 secondary signals for BMI, 28 for height, and 2 for WHRadjBMI in established loci. We show that 336 known BMI, 1,177 known height, and 143 known WHRadjBMI (combined) SNPs demonstrated suggestive transferability (nominal significance and effect estimate directional consistency) in Hispanic/Latino adults. Of these, 36 BMI, 124 height, and 11 WHRadjBMI SNPs were significant after trait-specific Bonferroni correction. Trans-ancestral meta-analysis of the three ancestries showed a small-to-moderate impact of uncorrected population stratification on the resulting effect size estimates. Our findings demonstrate that future studies may also benefit from leveraging diverse ancestries and differences in linkage disequilibrium patterns to discover novel loci and additional signals with less residual population stratification.

Factors Associated With Participation in Clinical Trials Among Patients With Lupus

BACKGROUND/OBJECTIVE

Participation rates for clinical trials, including lupus trials, in the United States are low, but are even lower for underrepresented minorities. The impact of underrepresentation in trials can be far-reaching and is problematic because female subjects of color with lupus experience greater morbidity and mortality. As such, the overarching goal of this study was to characterize the factors that influence participation in lupus clinical trials.

METHODS

The Lupus and Allied Diseases Association, the Lupus Foundation of America, and the Lupus Research Alliance collected data for their externally led Patient-Focused Drug Development Initiative-for the purpose of understanding and improving the rates of participation in lupus-related clinical trials. Participants completed a 46-question survey (in English or Spanish) electronically or on paper, which was distributed online or at lupus events. Logistic regression was used to test whether demographic and disease characteristics were associated with participation in past lupus trials.

RESULTS

Data were available for 2220 respondents. Black respondents with lupus were more likely, than their White and Hispanic counterparts, to have participated in past clinical trials (p < 0.05). Although not statistically significant, Hispanic respondents were also more likely to have participated than their White counterparts (odds ratio, 1.40; 95% confidence interval, 0.96-2.11). Both demographic (ie, race/ethnicity) and medical (ie, disease severity defined as more organ involvement) factors seem to be important determinants of participation in clinical trials (p < 0.05).

CONCLUSIONS

Combining the results from this study and prior research provides insight into recruitment strategies to increase participation rates of historically underrepresented minorities.

Identification of Novel Loci Shared by Juvenile Idiopathic Arthritis Subtypes through Integrative Genetic Analysis

Objective

Juvenile idiopathic arthritis (JIA) is the most common chronic immune‐mediated joint disease among children and encompasses a heterogeneous group of immune‐mediated joint disorders classified into 7 subtypes according to clinical presentation. However, phenotype overlap and biologic evidence suggest a shared mechanistic basis between subtypes. This study was undertaken to systematically investigate shared genetic underpinnings of JIA subtypes.

Methods

We performed a heterogeneity‐sensitive genome‐wide association study encompassing a total of 1,245 JIA cases (classified into 7 subtypes) and 9,250 controls, followed by fine‐mapping of candidate causal variants at each genome‐wide significant locus, functional annotation, and pathway and network analysis. We further identified candidate drug targets and drug repurposing opportunities by in silico analyses.

Results

In addition to the major histocompatibility complex locus, we identified 15 genome‐wide significant loci shared between at least 2 JIA subtypes, including 10 novel loci. Functional annotation indicated that candidate genes at these loci were expressed in diverse immune cell types.

Conclusion

This study identified novel genetic loci shared by JIA subtypes. Our findings identified candidate mechanisms underlying JIA subtypes and candidate targets with drug repurposing opportunities for JIA treatment.

Risk Factors Associated With Mortality and Neurologic Disability After Intracerebral Hemorrhage in a Racially and Ethnically Diverse Cohort

Introduction

Intracerebral hemorrhage (ICH) is the most severe subtype of stroke. Its mortality rate is high, and most survivors experience significant disability.

Objective

To assess primary patient risk factors associated with mortality and neurologic disability 3 months after ICH in a large, racially and ethnically balanced cohort.

Design, Setting, and Participants

This cohort study included participants from the Ethnic/Racial Variations of Intracerebral Hemorrhage (ERICH) study, which prospectively recruited 1000 non-Hispanic White, 1000 non-Hispanic Black, and 1000 Hispanic patients with spontaneous ICH to study the epidemiological characteristics and genomics associated with ICH. Participants included those with uniform data collection and phenotype definitions, centralized neuroimaging review, and telephone follow-up at 3 months. Analyses were completed in November 2021.

Exposures

Patient demographic and clinical characteristics as well as hospital event and imaging variables were examined, with characteristics meeting P < .20 considered candidates for a multivariate model. Elements included in the ICH score were specifically analyzed.

Main Outcomes and Measures

Individual characteristics were screened for association with 3-month outcome of neurologic disability or mortality, as assessed by a modified Rankin Scale (mRS) score of 4 or greater vs 3 or less under a logistic regression model. A total of 25 characteristics were tested in the final model, which minimized the Akaike information criterion. Analyses were repeated removing individuals who had withdrawal of care.

Results

A total of 2568 patients (mean [SD] age, 62.4 [14.7] years; 1069 [41.6%] women and 1499 [58.4%] men) had a 3-month outcome determination available, including death. The final logistic model had a significantly higher area under the receiver operating characteristics curve (C = 0.88) compared with ICH score alone (C = 0.76; P < .001). Among characteristics associated with neurologic disability and mortality were larger log ICH volume (OR, 2.74; 95% CI, 2.36-3.19; P < .001), older age (OR per 1-year increase, 1.04; 95% CI, 1.02-1.05; P < .001), pre-ICH mRS score (OR, 1.62; 95% CI, 1.41-1.87; P < .001), lobar location (OR, 0.22; 95% CI, 0.16-0.30; P < .001), and presence of infection (OR, 1.85; 95% CI, 1.42-2.41; P < .001).

Conclusions and Relevance

The findings of this cohort study validate ICH score elements and suggest additional baseline and interim patient characteristics were associated with variation in 3-month outcome.

Abstract TP137: Ethnic/racial Variations Of Intracerebral Hemorrhage Genetics (erich-gene) Study Protocol

Introduction: Epidemiologic studies of intracerebral hemorrhage (ICH) have consistently demonstrated variation in incidence, location, comorbidity burden, age-of-onset, and outcome by race and ethnicity, and genetic studies have identified differences in risk mediated by genetic risk factors such as Apolipoprotein E (APOE). We report the design and methods of the largest multi-ethnic genome-wide association study (GWAS) of ICH risk and outcome conducted to date.

Methods: The Ethnic/Racial Variations of ICH Genetics (ERICH-GENE) study is an international, multi-center, genetic case-control study of ICH. Cases are individuals with confirmed primary ICH with biosample availability and consent compatible with shareable genome-wide genotyping or previous genotyping. Central neuroimaging phenotype harmonization of case status, hemorrhage location, and imaging characteristics of cerebral small vessel disease including leukoaraiosis, atrophy, microbleeds, intraventricular hemorrhage severity and volume of ICH measurements will be performed. Controls are ICH-free individuals of compatible age and race/ethnicity from existing studies and biorepositories.

Results: As of August 2021, 2,002 ICH cases have been collected and genotyped, in addition to 2,558 already-genotyped ICH cases from the multi-ethnic ERICH study. 5,600 total new ICH cases will be genotyped under ERICH-GENE, with planned meta-analyses across existing ICH GWAS datasets and international biobanks totaling >20,000 cases. We are on schedule to meet our genotyping goal within the study period. Non-European ancestry cases are being prioritized for genotyping and ~66% of the previously genotyped ERICH cases are black or Hispanic. From available cases, a total of 10,621 neuroimaging studies have been uploaded for central adjudication to date with 6,278 having undergone harmonization.

Conclusions: ERICH-GENE is a large, multi-ethnic, international, centrally harmonized GWAS of ICH risk and outcome that will identify genetic risk factors across diverse populations for biological discovery and population-specific risk stratification.

Abstract WMP82: Diffusion Tensor Imaging Is Independently Associated With Mobility After Intracerebral Hemorrhage

Background: Intracerebral hemorrhage (ICH) patients may appear disabled despite intact cortical spinal tracts (CST). ICH lesions may also vary in volume and location relative to the CST with differential impacts on motor outcome. In this study, we evaluated whether diffusion tensor imaging (DTI) integrity of CST fiber tracking associates with mobility three months post ICH.

Methods: The Recovery and Outcomes from StrokE (ROSE) Study is a prospective, multi-center study of ICH serially assessed with motor, cognitive, and functional outcomes and MRI/DTI. Outcome data were collected by hospital chart review and interview/exam. Acute MRI images were analyzed by a central neuroimaging core for ICH location, size, intraventricular hemorrhage, and fractional anisotropy (FA) derived from DTI of the CST. Mobility was assessed at 3-month follow-up using the Barthel Mobility Index. We computed single variable logistic regression analyses (outcome = Barthel mobility at 3 months). We also computed a stepwise logistic model, with and without age, sex, and race, including variables that met P<0.1 to best fit.

Results: There were 76 ICH patients with DTI FA data available for analysis; 61 with complete data for the multivariate model. Poor mobility after ICH (N = 21; 27.6%) was associated with ICH volume, presence of IVH, Glasgow coma score, ipsilesional CST FA, ipsilesional CP FA and CP FA asymmetry index (Table 1). The multivariate model that adjusted for age, sex, and race identified ipsilesional CP FA, admission GCS and presence of IVH as associated with poor mobility at 3 months; notably, log ICH volume was not associated after controlling for CST injury (Table 1).

Conclusions: The presence of IVH and the GCS score at baseline were associated with poor mobility. In addition, the CP FA asymmetry index was independently associated with poor mobility even after accounting for log ICH volume, suggesting that DTI is informative of outcome beyond classic ICH outcome predictors.

Antifibrotic factor KLF4 is repressed by the miR-10/TFAP2A/TBX5 axis in dermal fibroblasts: insights from twins discordant for systemic sclerosis

OBJECTIVES

Systemic sclerosis (SSc) is a complex disease of unknown aetiology in which inflammation and fibrosis lead to multiple organ damage. There is currently no effective therapy that can halt the progression of fibrosis or reverse it, thus studies that provide novel insights into disease pathogenesis and identify novel potential therapeutic targets are critically needed.

METHODS

We used global gene expression and genome-wide DNA methylation analyses of dermal fibroblasts (dFBs) from a unique cohort of twins discordant for SSc to identify molecular features of this pathology. We validated the findings using in vitro, ex vivo and in vivo models.

RESULTS

Our results revealed distinct differentially expressed and methylated genes, including several transcription factors involved in stem cell differentiation and developmental programmes (KLF4, TBX5, TFAP2A and homeobox genes) and the microRNAs miR-10a and miR-10b which target several of these deregulated genes. We show that KLF4 expression is reduced in SSc dFBs and its expression is repressed by TBX5 and TFAP2A. We also show that KLF4 is antifibrotic, and its conditional knockout in fibroblasts promotes a fibrotic phenotype.

CONCLUSIONS

Our data support a role for epigenetic dysregulation in mediating SSc susceptibility in dFBs, illustrating the intricate interplay between CpG methylation, miRNAs and transcription factors in SSc pathogenesis, and highlighting the potential for future use of epigenetic modifiers as therapies.

Exploiting three-dimensional human hepatic constructs to investigate the impact of rs174537 on fatty acid metabolism

The Modern Western Diet has been associated with the rise in metabolic and inflammatory diseases, including obesity, diabetes, and cardiovascular disease. This has been attributed, in part, to the increase in dietary omega-6 polyunsaturated fatty acid (PUFA) consumption, specifically linoleic acid (LA), arachidonic acid (ARA), and their subsequent metabolism to pro-inflammatory metabolites which may be driving human disease. Conversion of dietary LA to ARA is regulated by genetic variants near and within the fatty acid desaturase (FADS) haplotype block, most notably single nucleotide polymorphism rs174537 is strongly associated with FADS1 activity and expression. This variant and others within high linkage disequilibrium may potentially explain the diversity in both diet and inflammatory mediators that drive chronic inflammatory disease in human populations. Mechanistic exploration into this phenomenon using human hepatocytes is limited by current two-dimensional culture models that poorly replicate in vivo functionality. Therefore, we aimed to develop and characterize a three-dimensional hepatic construct for the study of human PUFA metabolism. Primary human hepatocytes cultured in 3D hydrogels were characterized for their capacity to represent basic lipid processing functions, including lipid esterification, de novo lipogenesis, and cholesterol efflux. They were then exposed to control and LA-enriched media and reproducibly displayed allele-specific metabolic activity of FADS1, based on genotype at rs174537. Hepatocytes derived from individuals homozygous with the minor allele at rs174537 (i.e., TT) displayed the slowest metabolic conversion of LA to ARA and significantly reduced FADS1 and FADS2 expression. These results support the feasibility of using 3D human hepatic cultures for the study of human PUFA and lipid metabolism and relevant gene-diet interactions, thereby enabling future nutrition targets in humans.

Bulk and Single-Cell Profiling of Breast Tumors Identifies TREM-1 as a Dominant Immune Suppressive Marker Associated With Poor Outcomes

Background

Triggering receptor expressed on myeloid cells (TREM)-1 is a key mediator of innate immunity previously associated with the severity of inflammatory disorders, and more recently, the inferior survival of lung and liver cancer patients. Here, we investigated the prognostic impact and immunological correlates of TREM1 expression in breast tumors.

Methods

Breast tumor microarray and RNAseq expression profiles (n=4,364 tumors) were analyzed for associations between gene expression, tumor immune subtypes, distant metastasis-free survival (DMFS) and clinical response to neoadjuvant chemotherapy (NAC). Single-cell (sc)RNAseq was performed using the 10X Genomics platform. Statistical associations were assessed by logistic regression, Cox regression, Kaplan-Meier analysis, Spearman correlation, Student’s t-test and Chi-square test.

Results

In pre-treatment biopsies, TREM1 and known TREM-1 inducible cytokines (IL1B, IL8) were discovered by a statistical ranking procedure as top genes for which high expression was associated with reduced response to NAC, but only in the context of immunologically “hot” tumors otherwise associated with a high NAC response rate. In surgical specimens, TREM1 expression varied among tumor molecular subtypes, with highest expression in the more aggressive subtypes (Basal-like, HER2-E). High TREM1 significantly and reproducibly associated with inferior distant metastasis-free survival (DMFS), independent of conventional prognostic markers. Notably, the association between high TREM1 and inferior DMFS was most prominent in the subset of immunogenic tumors that exhibited the immunologically hot phenotype and otherwise associated with superior DMFS. Further observations from bulk and single-cell RNAseq analyses indicated that TREM1 expression was significantly enriched in polymorphonuclear myeloid-derived suppressor cells (PMN-MDSCs) and M2-like macrophages, and correlated with downstream transcriptional targets of TREM-1 (IL8, IL-1B, IL6, MCP-1, SPP1, IL1RN, INHBA) which have been previously associated with pro-tumorigenic and immunosuppressive functions.

Conclusions

Together, these findings indicate that increased TREM1 expression is prognostic of inferior breast cancer outcomes and may contribute to myeloid-mediated breast cancer progression and immune suppression.

Single-cell expression quantitative trait loci (eQTL) analysis of SLE-risk loci in lupus patient monocytes

Background

We performed expression quantitative trait locus (eQTL) analysis in single classical (CL) and non-classical (NCL) monocytes from patients with systemic lupus erythematosus (SLE) to quantify the impact of well-established genetic risk alleles on transcription at single-cell resolution.

Methods

Single-cell gene expression was quantified using qPCR in purified monocyte subpopulations (CD14⁺⁺CD16⁻ CL and CD14^dimCD16⁺ NCL) from SLE patients. Novel analysis methods were used to control for the within-person correlations observed, and eQTLs were compared between cell types and risk alleles.

Results

The SLE-risk alleles demonstrated significantly more eQTLs in NCLs as compared to CLs (p = 0.0004). There were 18 eQTLs exclusive to NCL cells, 5 eQTLs exclusive to CL cells, and only one shared eQTL, supporting large differences in the impact of the risk alleles between these monocyte subsets. The SPP1 and TNFAIP3 loci were associated with the greatest number of transcripts. Patterns of shared influence in which different SNPs impacted the same transcript also differed between monocyte subsets, with greater evidence for synergy in NCL cells. IRF1 expression demonstrated an on/off pattern, in which expression was zero in all of the monocytes studied from some individuals, and this pattern was associated with a number of SLE risk alleles. We observed corroborating evidence of this IRF1 expression pattern in public data sets.

Conclusions

We document multiple SLE-risk allele eQTLs in single monocytes which differ greatly between CL and NCL subsets. These data support the importance of the SPP1 and TNFAIP3 risk variants and the IRF1 transcript in SLE patient monocyte function.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13075-021-02660-2.

Nucleic Acid-Sensing and Interferon-Inducible Pathways Show Differential Methylation in MZ Twins Discordant for Lupus and Overexpression in Independent Lupus Samples: Implications for Pathogenic Mechanism and Drug Targeting

Systemic lupus erythematosus (SLE) is a chronic, multisystem, autoimmune inflammatory disease with genomic and non-genomic contributions to risk. We hypothesize that epigenetic factors are a significant contributor to SLE risk and may be informative for identifying pathogenic mechanisms and therapeutic targets. To test this hypothesis while controlling for genetic background, we performed an epigenome-wide analysis of DNA methylation in genomic DNA from whole blood in three pairs of female monozygotic (MZ) twins of European ancestry, discordant for SLE. Results were replicated on the same array in four cell types from a set of four Danish female MZ twin pairs discordant for SLE. Genes implicated by the epigenetic analyses were then evaluated in 10 independent SLE gene expression datasets from the Gene Expression Omnibus (GEO). There were 59 differentially methylated loci between unaffected and affected MZ twins in whole blood, including 11 novel loci. All but two of these loci were hypomethylated in the SLE twins relative to the unaffected twins. The genes harboring these hypomethylated loci exhibited increased expression in multiple independent datasets of SLE patients. This pattern was largely consistent regardless of disease activity, cell type, or renal tissue type. The genes proximal to CpGs exhibiting differential methylation (DM) in the SLE-discordant MZ twins and exhibiting differential expression (DE) in independent SLE GEO cohorts (DM-DE genes) clustered into two pathways: the nucleic acid-sensing pathway and the type I interferon pathway. The DM-DE genes were also informatically queried for potential gene–drug interactions, yielding a list of 41 drugs including a known SLE therapy. The DM-DE genes delineate two important biologic pathways that are not only reflective of the heterogeneity of SLE but may also correlate with distinct IFN responses that depend on the source, type, and location of nucleic acid molecules and the activated receptors in individual patients. Cell- and tissue-specific analyses will be critical to the understanding of genetic factors dysregulating the nucleic acid-sensing and IFN pathways and whether these factors could be appropriate targets for therapeutic intervention.

Metabolomic architecture of obesity implicates metabolonic lactone sulfate in cardiometabolic disease

OBJECTIVE

Identify and characterize circulating metabolite profiles associated with adiposity to inform precision medicine.

METHODS

Untargeted plasma metabolomic profiles in the Insulin Resistance Atherosclerosis Family Study (IRASFS) Mexican American cohort (n = 1108) were analyzed for association with anthropometric (body mass index, BMI; waist circumference, WC; waist-to-hip ratio, WHR) and computed tomography measures (visceral adipose tissue, VAT; subcutaneous adipose tissue, SAT; visceral-to-subcutaneous ratio, VSR) of adiposity. Genetic data, inclusive of genome-wide array-based genotyping, whole exome sequencing (WES) and whole genome sequencing (WGS), were evaluated to identify the genetic contributors. Phenotypic and genetic association signals were replicated across ancestries. Transcriptomic data were analyzed to explore the relationship between genetic and metabolomic data.

RESULTS

A partially characterized metabolite, tentatively named metabolonic lactone sulfate (X-12063), was consistently associated with BMI, WC, WHR, VAT, and SAT in IRASFS Mexican Americans (P_MA <2.02 × 10^-27). Trait associations were replicated in IRASFS African Americans (P_AA < 1.12 × 10^-07). Expanded analyses revealed associations with multiple phenotypic measures of cardiometabolic health, e.g. insulin sensitivity (S_I), triglycerides (TG), diastolic blood pressure (DBP) and plasminogen activator inhibitor-1 (PAI-1) in both ancestries. Metabolonic lactone sulfate levels were heritable (h² > 0.47), and a significant genetic signal at the ZSCAN25/CYP3A5 locus (P_MA = 9.00 × 10^-41, P_AA = 2.31 × 10^-10) was observed, highlighting a putative functional variant (rs776746, CYP3A5∗3). Transcriptomic analysis in the African American Genetics of Metabolism and Expression (AAGMEx) cohort supported the association of CYP3A5 with metabolonic lactone sulfate levels (P_FDR = 6.64 × 10^-07).

CONCLUSIONS

Variant rs776746 is associated with a decrease in the transcript levels of CYP3A5, which in turn is associated with increased metabolonic lactone sulfate levels and poor cardiometabolic health.

Ethnic and Racial Variation in Intracerebral Hemorrhage Risk Factors and Risk Factor Burden

Importance

Black and Hispanic individuals have an increased risk of intracerebral hemorrhage (ICH) compared with their White counterparts, but no large studies of ICH have been conducted in these disproportionately affected populations.

Objective

To examine the prevalence, odds, and population attributable risk (PAR) percentage for established and novel risk factors for ICH, stratified by ICH location and racial/ethnic group.

Design, Setting, and Participants

The Ethnic/Racial Variations of Intracerebral Hemorrhage Study was a case-control study of ICH among 3000 Black, Hispanic, and White individuals who experienced spontaneous ICH (1000 cases in each group). Recruitment was conducted between September 2009 and July 2016 at 19 US sites comprising 42 hospitals. Control participants were identified through random digit dialing and were matched to case participants by age (±5 years), sex, race/ethnicity, and geographic area. Data analyses were conducted from January 2019 to May 2020.

Main Outcomes and Measures

Case and control participants underwent a standardized interview, physical measurement for body mass index, and genotyping for the ɛ2 and ɛ4 alleles of APOE, the gene encoding apolipoprotein E. Prevalence, multivariable adjusted odds ratio (OR), and PAR percentage were calculated for each risk factor in the entire ICH population and stratified by racial/ethnic group and by lobar or nonlobar location.

Results

There were 1000 Black patients (median [interquartile range (IQR)] age, 57 [50-65] years, 425 [42.5%] women), 1000 Hispanic patients (median [IQR] age, 58 [49-69] years; 373 [37.3%] women), and 1000 White patients (median [IQR] age, 71 [59-80] years; 437 [43.7%] women). The mean (SD) age of patients with ICH was significantly lower among Black and Hispanic patients compared with White patients (eg, lobar ICH: Black, 62.2 [15.2] years; Hispanic, 62.5 [15.7] years; White, 71.0 [13.3] years). More than half of all ICH in Black and Hispanic patients was associated with treated or untreated hypertension (PAR for treated hypertension, Black patients: 53.6%; 95% CI, 46.4%-59.8%; Hispanic patients: 46.5%; 95% CI, 40.6%-51.8%; untreated hypertension, Black patients: 45.5%; 95% CI, 39.%-51.1%; Hispanic patients: 42.7%; 95% CI, 37.6%-47.3%). Lack of health insurance also had a disproportionate association with the PAR percentage for ICH in Black and Hispanic patients (Black patients: 21.7%; 95% CI, 17.5%-25.7%; Hispanic patients: 30.2%; 95% CI, 26.1%-34.1%; White patients: 5.8%; 95% CI, 3.3%-8.2%). A high sleep apnea risk score was associated with both lobar (OR, 1.68; 95% CI, 1.36-2.06) and nonlobar (OR, 1.62; 95% CI, 1.37-1.91) ICH, and high cholesterol was inversely associated only with nonlobar ICH (OR, 0.60; 95% CI, 0.52-0.70); both had no interactions with race and ethnicity. In contrast to the association between the ɛ2 and ɛ4 alleles of APOE and ICH in White individuals (eg, presence of APOE ɛ2 allele: OR, 1.84; 95% CI, 1.34-2.52), APOE alleles were not associated with lobar ICH among Black or Hispanic individuals.

Conclusions and Relevance

This study found sleep apnea as a novel risk factor for ICH. The results suggest a strong contribution from inadequately treated hypertension and lack of health insurance to the disproportionate burden and earlier onset of ICH in Black and Hispanic populations. These findings emphasize the importance of addressing modifiable risk factors and the social determinants of health to reduce health disparities.

Multiethnic Genome-Wide Association Study of Subclinical Atherosclerosis in Individuals With Type 2 Diabetes

BACKGROUND

Coronary artery calcification (CAC) and carotid artery intima-media thickness (cIMT) are measures of subclinical atherosclerosis in asymptomatic individuals and strong risk factors for cardiovascular disease. Type 2 diabetes (T2D) is an independent cardiovascular disease risk factor that accelerates atherosclerosis.

METHODS

We performed meta-analyses of genome-wide association studies in up to 2500 T2D individuals of European ancestry (EA) and 1590 T2D individuals of African ancestry with or without exclusion of prevalent cardiovascular disease, for CAC measured by cardiac computed tomography, and 3608 individuals of EA and 838 individuals of African ancestry with T2D for cIMT measured by ultrasonography within the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium.

RESULTS

We replicated 2 loci (rs9369640 and rs9349379 near PHACTR1 and rs10757278 near CDKN2B) for CAC and one locus for cIMT (rs7412 and rs445925 near APOE-APOC1) that were previously reported in the general EA populations. We identified one novel CAC locus (rs8000449 near CSNK1A1L/LINC00547/POSTN at 13q13.3) at P=2.0×10-8 in EA. No additional loci were identified with the meta-analyses of EA and African ancestry. The expression quantitative trait loci analysis with nearby expressed genes derived from arterial wall and metabolic tissues from the Genotype-Tissue Expression project pinpoints POSTN, encoding a matricellular protein involved in bone formation and bone matrix organization, as the potential candidate gene at this locus. In addition, we found significant associations (P<3.1×10-4) for 3 previously reported coronary artery disease loci for these subclinical atherosclerotic phenotypes (rs2891168 near CDKN2B-AS1 and rs11170820 near FLJ12825 for CAC, and rs7412 near APOE for cIMT).

CONCLUSIONS

Our results provide potential biological mechanisms that could link CAC and cIMT to increased cardiovascular disease risk in individuals with T2D.

Genetic landscape of Gullah African Americans

OBJECTIVES

Gullah African Americans are descendants of formerly enslaved Africans living in the Sea Islands along the coast of the southeastern U.S., from North Carolina to Florida. Their relatively high numbers and geographic isolation were conducive to the development and preservation of a unique culture that retains deep African features. Although historical evidence supports a West-Central African ancestry for the Gullah, linguistic and cultural evidence of a connection to Sierra Leone has led to the suggestion of this country/region as their ancestral home. This study sought to elucidate the genetic structure and ancestry of the Gullah.

MATERIALS AND METHODS

We leveraged whole-genome genotype data from Gullah, African Americans from Jackson, Mississippi, African populations from Sierra Leone, and population reference panels from Africa and Europe to infer population structure, ancestry proportions, and global estimates of admixture.

RESULTS

Relative to non-Gullah African Americans from the Southeast US, the Gullah exhibited higher mean African ancestry, lower European admixture, a similarly small Native American contribution, and increased male-biased European admixture. A slightly tighter bottleneck in the Gullah 13 generations ago suggests a largely shared demographic history with non-Gullah African Americans. Despite a slightly higher relatedness to populations from Sierra Leone, our data demonstrate that the Gullah are genetically related to many West African populations.

DISCUSSION

This study confirms that subtle differences in African American population structure exist at finer regional levels. Such observations can help to inform medical genetics research in African Americans, and guide the interpretation of genetic data used by African Americans seeking to explore ancestral identities.