Aire mediates tolerance to insulin through thymic trimming of high-affinity T cell clones.

Insulin is a central autoantigen in the pathogenesis of T1D, and thymic epithelial cell expression of insulin under the control of the Autoimmune Regulator (Aire) is thought to be a key component of maintaining tolerance to insulin. In spite of this general working model, direct detection of this thymic selection on insulin-specific T cells has been somewhat elusive. Here, we used a combination of highly sensitive T cell receptor transgenic models for detecting thymic selection and sorting and sequencing of Insulin-specific CD4+ T cells from Aire-deficient mice as a strategy to further define their selection. This analysis revealed a number of unique t cell receptor (TCR) clones in Aire-deficient hosts with high affinity for insulin/major histocompatibility complex (MHC) ligands. We then modeled the thymic selection of one of these clones in Aire-deficient versus wild-type hosts and found that this model clone could escape thymic negative selection in the absence of thymic Aire. Together, these results suggest that thymic expression of insulin plays a key role in trimming and removing high-affinity insulin-specific T cells from the repertoire to help promote tolerance.

Tracking the role of Aire in immune tolerance to the eye with a TCR transgenic mouse model.

Roughly one-half of mice with partial defects in two immune tolerance pathways (AireGW/+Lyn-/- mice) spontaneously develop severe damage to their retinas due to T cell reactivity to Aire-regulated interphotoreceptor retinoid-binding protein (IRBP). Single-cell T cell receptor (TCR) sequencing of CD4+ T cells specific for a predominate epitope of IRBP showed a remarkable diversity of autoantigen-specific TCRs with greater clonal expansions in mice with disease. TCR transgenic mice made with an expanded IRBP-specific TCR (P2.U2) of intermediate affinity exhibited strong but incomplete negative selection of thymocytes. This negative selection was absent in IRBP-/- mice and greatly defective in AireGW/+ mice. Most P2.U2+/- mice and all P2.U.2+/-AireGW/+ mice rapidly developed inflammation of the retina and adjacent uvea (uveitis). Aire-dependent IRBP expression in the thymus also promoted Treg differentiation, but the niche for this fate determination was small, suggesting differences in antigen presentation leading to negative selection vs. thymic Treg differentiation and a stronger role for negative selection in preventing autoimmune disease in the retina.

RNAi screening comes of age: improved techniques and complementary approaches.

Gene silencing through sequence-specific targeting of mRNAs by RNAi has enabled genome-wide functional screens in cultured cells and in vivo in model organisms. These screens have resulted in the identification of new cellular pathways and potential drug targets. Considerable progress has been made to improve the quality of RNAi screen data through the development of new experimental and bioinformatics approaches. The recent availability of genome-editing strategies, such as the CRISPR (clustered regularly interspaced short palindromic repeats)-Cas9 system, when combined with RNAi, could lead to further improvements in screen data quality and follow-up experiments, thus promoting our understanding of gene function and gene regulatory networks.

Methods for mediation analysis with high-dimensional DNA methylation data: Possible choices and comparisons.

Epigenetic researchers often evaluate DNA methylation as a potential mediator of the effect of social/environmental exposures on a health outcome. Modern statistical methods for jointly evaluating many mediators have not been widely adopted. We compare seven methods for high-dimensional mediation analysis with continuous outcomes through both diverse simulations and analysis of DNAm data from a large multi-ethnic cohort in the United States, while providing an R package for their seamless implementation and adoption. Among the considered choices, the best-performing methods for detecting active mediators in simulations are the Bayesian sparse linear mixed model (BSLMM) and high-dimensional mediation analysis (HDMA); while the preferred methods for estimating the global mediation effect are high-dimensional linear mediation analysis (HILMA) and principal component mediation analysis (PCMA). We provide guidelines for epigenetic researchers on choosing the best method in practice and offer suggestions for future methodological development.

Whole-exome sequencing study identifies four novel gene loci associated with diabetic kidney disease.

Diabetic kidney disease (DKD) is recognized as an important public health challenge. However, its genomic mechanisms are poorly understood. To identify rare variants for DKD, we conducted a whole-exome sequencing (WES) study leveraging large cohorts well-phenotyped for chronic kidney disease and diabetes. Our two-stage WES study included 4372 European and African ancestry participants from the Chronic Renal Insufficiency Cohort and Atherosclerosis Risk in Communities studies (stage 1) and 11 487 multi-ancestry Trans-Omics for Precision Medicine participants (stage 2). Generalized linear mixed models, which accounted for genetic relatedness and adjusted for age, sex and ancestry, were used to test associations between single variants and DKD. Gene-based aggregate rare variant analyses were conducted using an optimized sequence kernel association test implemented within our mixed model framework. We identified four novel exome-wide significant DKD-related loci through initiating diabetes. In single-variant analyses, participants carrying a rare, in-frame insertion in the DIS3L2 gene (rs141560952) exhibited a 193-fold increased odds [95% confidence interval (CI): 33.6, 1105] of DKD compared with noncarriers (P = 3.59 × 10-9). Likewise, each copy of a low-frequency KRT6B splice-site variant (rs425827) conferred a 5.31-fold higher odds (95% CI: 3.06, 9.21) of DKD (P = 2.72 × 10-9). Aggregate gene-based analyses further identified ERAP2 (P = 4.03 × 10-8) and NPEPPS (P = 1.51 × 10-7), which are both expressed in the kidney and implicated in renin-angiotensin-aldosterone system modulated immune response. In the largest WES study of DKD, we identified novel rare variant loci attaining exome-wide significance. These findings provide new insights into the molecular mechanisms underlying DKD.

The construction of cross-population polygenic risk scores using transfer learning.

As most existing genome-wide association studies (GWASs) were conducted in European-ancestry cohorts, and as the existing polygenic risk score (PRS) models have limited transferability across ancestry groups, PRS research on non-European-ancestry groups needs to make efficient use of available data until we attain large sample sizes across all ancestry groups. Here we propose a PRS method using transfer learning techniques. Our approach, TL-PRS, uses gradient descent to fine-tune the baseline PRS model from an ancestry group with large sample GWASs to the dataset of target ancestry. In our application of constructing PRS for seven quantitative and two dichotomous traits for 10,285 individuals of South Asian ancestry and 8,168 individuals of African ancestry in UK Biobank, TL-PRS using PRS-CS as a baseline method obtained 25% average relative improvement for South Asian samples and 29% for African samples compared to the standard PRS-CS method in terms of predicted R2. Our approach increases the transferability of PRSs across ancestries and thereby helps reduce existing inequities in genetics research.

METRO: Multi-ancestry transcriptome-wide association studies for powerful gene-trait association detection.

Integrative analysis of genome-wide association studies (GWASs) and gene expression studies in the form of a transcriptome-wide association study (TWAS) has the potential to better elucidate the molecular mechanisms underlying disease etiology. Here we present a method, METRO, that can leverage gene expression data collected from multiple genetic ancestries to enhance TWASs. METRO incorporates expression prediction models constructed in different genetic ancestries through a likelihood-based inference framework, producing calibrated p values with substantially improved TWAS power. We illustrate the benefits of METRO in both simulations and applications to seven complex traits and diseases obtained from four GWASs. These GWASs include two of primarily European ancestry (n = 188,577 and 339,226) and two of primarily African ancestry (n = 42,752 and 23,827). In the real data applications, we leverage gene expression data measured on 1,032 African Americans and 801 European Americans from the Genetic Epidemiology Network of Arteriopathy (GENOA) study to identify a substantially larger number of gene-trait associations as compared to existing TWAS approaches. The benefits of METRO are most prominent in applications to GWASs of African ancestry where the sample size is much smaller than GWASs of European ancestry and where a more powerful TWAS method is crucial. Among the identified associations are high-density lipoprotein-associated genes including PLTP and PPARG that are critical for maintaining lipid homeostasis and the type II diabetes-associated gene MAPT that supports microtubule-associated protein tau as a key component underlying impaired insulin secretion.

Biological Age, Chronological Age and Survival in Pulmonary Fibrosis: A Causal Mediation Analysis.

RATIONALE: Accelerated biological aging has been implicated in the development of interstitial lung disease (ILD) and other diseases of aging but remains poorly understood. OBJECTIVES: To identify plasma proteins that mediate the relationship between chronological age and survival association in patients with ILD. METHODS: Causal mediation analysis was performed to identify plasma proteins that mediated the chronological age-survival relationship in an idiopathic pulmonary fibrosis (IPF) discovery cohort. Proteins mediating this relationship after adjustment for false discovery were advanced for testing in an independent ILD validation cohort and explored in a chronic obstructive pulmonary disease (COPD) cohort. A proteomic-based measure of biological age was constructed and survival analysis performed assessing the impact of biological age and peripheral blood telomere length on the chronological age-survival relationship. RESULTS: Twenty-two proteins mediated the chronological age-survival relationship after adjustment for false discovery in the IPF discovery cohort (n=874), with nineteen remaining significant mediators of this relationship in the ILD validation cohort (n=983) and one mediating this relationship in the COPD cohort. Latent transforming growth factor beta binding protein 2 and ectodysplasin A2 receptor showed the strongest mediation across cohorts. A proteomic measure of biological age completely attenuated the chronological age-survival association and better discriminated survival than chronological age. Results were robust to adjustment for peripheral blood telomere length, which did not mediate the chronological age-survival relationship. CONCLUSIONS: Molecular measures of aging completely mediate the relationship between chronological age and survival, suggesting that chronological age has no direct effect on ILD survival.

Methods for large-scale single mediator hypothesis testing: Possible choices and comparisons.

Mediation hypothesis testing for a large number of mediators is challenging due to the composite structure of the null hypothesis, H 0 : α ß = 0 ${H}_{0}:\alpha \beta =0$ ( α $\alpha $ : effect of the exposure on the mediator after adjusting for confounders; ß $\beta $ : effect of the mediator on the outcome after adjusting for exposure and confounders). In this paper, we reviewed three classes of methods for large-scale one at a time mediation hypothesis testing. These methods are commonly used for continuous outcomes and continuous mediators assuming there is no exposure-mediator interaction so that the product α ß $\alpha \beta $ has a causal interpretation as the indirect effect. The first class of methods ignores the impact of different structures under the composite null hypothesis, namely, (1) α = 0 , ß ≠ 0 $\alpha =0,\beta \ne 0$ ; (2) α ≠ 0 , ß = 0 $\alpha \ne 0,\beta =0$ ; and (3) α = ß = 0 $\alpha =\beta =0$ . The second class of methods weights the reference distribution under each case of the null to form a mixture reference distribution. The third class constructs a composite test statistic using the three p values obtained under each case of the null so that the reference distribution of the composite statistic is approximately U ( 0 , 1 ) $U(0,1)$ . In addition to these existing methods, we developed the Sobel-comp method belonging to the second class, which uses a corrected mixture reference distribution for Sobel's test statistic. We performed extensive simulation studies to compare all six methods belonging to these three classes in terms of the false positive rates (FPRs) under the null hypothesis and the true positive rates under the alternative hypothesis. We found that the second class of methods which uses a mixture reference distribution could best maintain the FPRs at the nominal level under the null hypothesis and had the greatest true positive rates under the alternative hypothesis. We applied all methods to study the mediation mechanism of DNA methylation sites in the pathway from adult socioeconomic status to glycated hemoglobin level using data from the Multi-Ethnic Study of Atherosclerosis (MESA). We provide guidelines for choosing the optimal mediation hypothesis testing method in practice and develop an R package medScan available on the CRAN for implementing all the six methods.

Epigenetic and integrative cross-omics analyses of cerebral white matter hyperintensities on MRI.

Cerebral white matter hyperintensities on MRI are markers of cerebral small vessel disease, a major risk factor for dementia and stroke. Despite the successful identification of multiple genetic variants associated with this highly heritable condition, its genetic architecture remains incompletely understood. More specifically, the role of DNA methylation has received little attention. We investigated the association between white matter hyperintensity burden and DNA methylation in blood at ∼450 000 cytosine-phosphate-guanine (CpG) sites in 9732 middle-aged to older adults from 14 community-based studies. Single CpG and region-based association analyses were carried out. Functional annotation and integrative cross-omics analyses were performed to identify novel genes underlying the relationship between DNA methylation and white matter hyperintensities. We identified 12 single CpG and 46 region-based DNA methylation associations with white matter hyperintensity burden. Our top discovery single CpG, cg24202936 (P = 7.6 × 10-8), was associated with F2 expression in blood (P = 6.4 × 10-5) and co-localized with FOLH1 expression in brain (posterior probability = 0.75). Our top differentially methylated regions were in PRMT1 and in CCDC144NL-AS1, which were also represented in single CpG associations (cg17417856 and cg06809326, respectively). Through Mendelian randomization analyses cg06809326 was putatively associated with white matter hyperintensity burden (P = 0.03) and expression of CCDC144NL-AS1 possibly mediated this association. Differentially methylated region analysis, joint epigenetic association analysis and multi-omics co-localization analysis consistently identified a role of DNA methylation near SH3PXD2A, a locus previously identified in genome-wide association studies of white matter hyperintensities. Gene set enrichment analyses revealed functions of the identified DNA methylation loci in the blood-brain barrier and in the immune response. Integrative cross-omics analysis identified 19 key regulatory genes in two networks related to extracellular matrix organization, and lipid and lipoprotein metabolism. A drug-repositioning analysis indicated antihyperlipidaemic agents, more specifically peroxisome proliferator-activated receptor-alpha, as possible target drugs for white matter hyperintensities. Our epigenome-wide association study and integrative cross-omics analyses implicate novel genes influencing white matter hyperintensity burden, which converged on pathways related to the immune response and to a compromised blood-brain barrier possibly due to disrupted cell-cell and cell-extracellular matrix interactions. The results also suggest that antihyperlipidaemic therapy may contribute to lowering risk for white matter hyperintensities possibly through protection against blood-brain barrier disruption.

Effect of apolipoprotein E ε4 and its modification by sociodemographic characteristics on cognitive measures in South Asians from LASI-DAD.

BACKGROUND: We investigated the effects of apolipoprotein E (APOE) ε4 and its interactions with sociodemographic characteristics on cognitive measures in South Asians from the Diagnostic Assessment of Dementia for the Longitudinal Aging Study of India (LASI-DAD). METHODS: Linear regression was used to assess the association between APOE ε4 and global- and domain-specific cognitive function in 2563 participants (mean age 69.6 ± 7.3 years; 53% female). Effect modification by age, sex, and education were explored using interaction terms and subgroup analyses. RESULTS: APOE ε4 was inversely associated with most cognitive measures (p < 0.05). This association was stronger with advancing age for the Hindi Mental State Examination (HMSE) score (ßε4×age = -0.44, p = 0.03), orientation (ßε4×age = -0.07, p = 0.01), and language/fluency (ßε4×age = -0.07, p = 0.01), as well as in females for memory (ßε4×male = 0.17, p = 0.02) and language/fluency (ßε4×male = 0.12, p = 0.03). DISCUSSION: APOE Îµ4 is associated with lower cognitive function in South Asians from India, with a more pronounced impact observed in females and older individuals. HIGHLIGHTS: APOE Îµ4 carriers had lower global and domain-specific cognitive performance. Females and older individuals may be more susceptible to ε4 effects. For most cognitive measures, there was no interaction between ε4 and education.

Associations of Early-Life Adversity With Later-Life Epigenetic Aging Profiles in the Multi-Ethnic Study of Atherosclerosis.

Epigenetic biomarkers of accelerated aging have been widely used to predict disease risk and may enhance our understanding of biological mechanisms between early-life adversity and disparities in aging. With respect to childhood adversity, most studies have used parental education or childhood disadvantage and/or have not examined the role played by socioemotional or physical abuse and trauma in epigenetic profiles at older ages. This study leveraged data from the Multi-Ethnic Study of Atherosclerosis (MESA) on experiences of threat and deprivation in participants' early lives (i.e., before the age of 18 years) to examine whether exposure to specific dimensions of early-life adversity is associated with epigenetic profiles at older ages that are indicative of accelerated biological aging. The sample included 842 MESA respondents with DNA methylation data collected between 2010 and 2012 who answered questions on early-life adversities in a 2018-2019 telephone follow-up. We found that experiences of deprivation, but not threat, were associated with later-life GrimAge epigenetic aging signatures that were developed to predict mortality risk. Results indicated that smoking behavior partially mediates this association, which suggests that lifestyle behaviors may act as downstream mechanisms between parental deprivation in early life and accelerated epigenetic aging in later life.

Epigenome-wide association study of mitochondrial genome copy number.

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P < 1 × 10-7), with a 0.7-3.0 standard deviation increase (3 CpGs) or decrease (18 CpGs) in mtDNA CN corresponding to a 1% increase in DNA methylation. Several significant CpGs have been reported to be associated with at least two risk factors (e.g. chronological age or smoking) for cardiovascular disease (CVD). Five genes [PR/SET domain 16, nuclear receptor subfamily 1 group H member 3 (NR1H3), DNA repair protein, DNA polymerase kappa and decaprenyl-diphosphate synthase subunit 2], which harbor nine significant CpGs, are known to be involved in mitochondrial biosynthesis and functions. For example, NR1H3 encodes a transcription factor that is differentially expressed during an adipose tissue transition. The methylation level of cg09548275 in NR1H3 was negatively associated with mtDNA CN (effect size = -1.71, P = 4 × 10-8) and was positively associated with the NR1H3 expression level (effect size = 0.43, P = 0.0003), which indicates that the methylation level in NR1H3 may underlie the relationship between mtDNA CN, the NR1H3 transcription factor and energy expenditure. In summary, the study results suggest that mtDNA CN variation in whole blood is associated with DNA methylation levels in genes that are involved in a wide range of mitochondrial activities. These findings will help reveal molecular mechanisms between mtDNA CN and CVD.

Genetic Architecture of Gene Expression in European and African Americans: An eQTL Mapping Study in GENOA.

Most existing expression quantitative trait locus (eQTL) mapping studies have been focused on individuals of European ancestry and are underrepresented in other populations including populations with African ancestry. Lack of large-scale well-powered eQTL mapping studies in populations with African ancestry can both impede the dissemination of eQTL mapping results that would otherwise benefit individuals with African ancestry and hinder the comparable analysis for understanding how gene regulation is shaped through evolution. We fill this critical knowledge gap by performing a large-scale in-depth eQTL mapping study on 1,032 African Americans (AA) and 801 European Americans (EA) in the GENOA cohort. We identified a total of 354,931 eSNPs in AA and 371,309 eSNPs in EA, with 112,316 eSNPs overlapped between the two. We found that eQTL harboring genes (eGenes) are enriched in metabolic pathways and tend to have higher SNP heritability compared to non-eGenes. We found that eGenes that are common in the two populations tend to be less conserved than eGenes that are unique to one population, which are less conserved than non-eGenes. Through conditional analysis, we found that eGenes in AA tend to harbor more independent eQTLs than eGenes in EA, suggesting potentially diverse genetic architecture underlying expression variation in the two populations. Finally, the large sample sizes in GENOA allow us to construct accurate expression prediction models in both AA and EA, facilitating powerful transcriptome-wide association studies. Overall, our results represent an important step toward revealing the genetic architecture underlying expression variation in African Americans.

Circulating metabolites may illustrate relationship of alcohol consumption with cardiovascular disease.

BACKGROUND: Metabolite signatures of long-term alcohol consumption are lacking. To better understand the molecular basis linking alcohol drinking and cardiovascular disease (CVD), we investigated circulating metabolites associated with long-term alcohol consumption and examined whether these metabolites were associated with incident CVD. METHODS: Cumulative average alcohol consumption (g/day) was derived from the total consumption of beer, wine, and liquor on average of 19 years in 2428 Framingham Heart Study Offspring participants (mean age 56 years, 52% women). We used linear mixed models to investigate the associations of alcohol consumption with 211 log-transformed plasma metabolites, adjusting for age, sex, batch, smoking, diet, physical activity, BMI, and familial relationship. Cox models were used to test the association of alcohol-related metabolite scores with fatal and nonfatal incident CVD (myocardial infarction, coronary heart disease, stroke, and heart failure). RESULTS: We identified 60 metabolites associated with cumulative average alcohol consumption (p < 0.05/211 ≈ 0.00024). For example, 1 g/day increase of alcohol consumption was associated with higher levels of cholesteryl esters (e.g., CE 16:1, beta = 0.023 ± 0.002, p = 6.3e - 45) and phosphatidylcholine (e.g., PC 32:1, beta = 0.021 ± 0.002, p = 3.1e - 38). Survival analysis identified that 10 alcohol-associated metabolites were also associated with a differential CVD risk after adjusting for age, sex, and batch. Further, we built two alcohol consumption weighted metabolite scores using these 10 metabolites and showed that, with adjustment age, sex, batch, and common CVD risk factors, the two scores had comparable but opposite associations with incident CVD, hazard ratio 1.11 (95% CI = [1.02, 1.21], p = 0.02) vs 0.88 (95% CI = [0.78, 0.98], p = 0.02). CONCLUSIONS: We identified 60 long-term alcohol consumption-associated metabolites. The association analysis with incident CVD suggests a complex metabolic basis between alcohol consumption and CVD.

Identification of MicroRNAs That Stabilize p53 in Human Papillomavirus-Positive Cancer Cells.

Etiologically, 5% of all cancers worldwide are caused by the high-risk human papillomaviruses (hrHPVs). These viruses encode two oncoproteins (E6 and E7) whose expression is required for cancer initiation and maintenance. Among their cellular targets are the p53 and the retinoblastoma tumor suppressor proteins. Inhibition of the hrHPV E6-mediated ubiquitylation of p53 through the E6AP ubiquitin ligase results in the stabilization of p53, leading to cellular apoptosis. We utilized a live cell high-throughput screen to determine whether exogenous microRNA (miRNA) transfection had the ability to stabilize p53 in hrHPV-positive cervical cancer cells expressing a p53-fluorescent protein as an in vivo reporter of p53 stability. Among the miRNAs whose transfection resulted in the greatest p53 stabilization was 375-3p, which has previously been reported to stabilize p53 in HeLa cells, providing validation of the screen. The top 32 miRNAs, in addition to 375-3p, were further assessed using a second cell-based p53 stability reporter system, as well as in nonreporter HeLa cells to examine their effects on endogenous p53 protein levels, resulting in the identification of 23 miRNAs whose transfection increased p53 levels in HeLa cells. While a few miRNAs that stabilized p53 led to decreases in E6AP protein levels, all targeted HPV oncoprotein expression. We further examined subsets of these miRNAs for their abilities to induce apoptosis and determined whether it was p53-mediated. The introduction of specific miRNAs revealed surprisingly heterogeneous responses in different cell lines. Nonetheless, some of the miRNAs described here have potential as therapeutics for treating HPV-positive cancers. IMPORTANCE Human papillomaviruses cause approximately 5% of all cancers worldwide and encode genes that contribute to both the initiation and maintenance of these cancers. The viral oncoprotein E6 is expressed in all HPV-positive cancers and functions by targeting the degradation of p53 through the engagement of the cellular ubiquitin ligase E6AP. Inhibiting the degradation of p53 leads to apoptosis in HPV-positive cancer cells. Using a high-throughput live cell assay, we identified several miRNAs whose transfection stabilize p53 in HPV-positive cells. These miRNAs have the potential to be used in the treatment of HPV-positive cancers.

Gene-mapping study of extremes of cerebral small vessel disease reveals TRIM47 as a strong candidate.

Cerebral small vessel disease is a leading cause of stroke and a major contributor to cognitive decline and dementia, but our understanding of specific genes underlying the cause of sporadic cerebral small vessel disease is limited. We report a genome-wide association study and a whole-exome association study on a composite extreme phenotype of cerebral small vessel disease derived from its most common MRI features: white matter hyperintensities and lacunes. Seventeen population-based cohorts of older persons with MRI measurements and genome-wide genotyping (n = 41 326), whole-exome sequencing (n = 15 965), or exome chip (n = 5249) data contributed 13 776 and 7079 extreme small vessel disease samples for the genome-wide association study and whole-exome association study, respectively. The genome-wide association study identified significant association of common variants in 11 loci with extreme small vessel disease, of which the chr12q24.11 locus was not previously reported to be associated with any MRI marker of cerebral small vessel disease. The whole-exome association study identified significant associations of extreme small vessel disease with common variants in the 5' UTR region of EFEMP1 (chr2p16.1) and one probably damaging common missense variant in TRIM47 (chr17q25.1). Mendelian randomization supports the causal association of extensive small vessel disease severity with increased risk of stroke and Alzheimer's disease. Combined evidence from summary-based Mendelian randomization studies and profiling of human loss-of-function allele carriers showed an inverse relation between TRIM47 expression in the brain and blood vessels and extensive small vessel disease severity. We observed significant enrichment of Trim47 in isolated brain vessel preparations compared to total brain fraction in mice, in line with the literature showing Trim47 enrichment in brain endothelial cells at single cell level. Functional evaluation of TRIM47 by small interfering RNAs-mediated knockdown in human brain endothelial cells showed increased endothelial permeability, an important hallmark of cerebral small vessel disease pathology. Overall, our comprehensive gene-mapping study and preliminary functional evaluation suggests a putative role of TRIM47 in the pathophysiology of cerebral small vessel disease, making it an important candidate for extensive in vivo explorations and future translational work.

Cardiometabolic disease and obesity patterns differentially predict acute kidney injury after total joint replacement: a retrospective analysis.

BACKGROUND: Acute kidney injury (AKI) is a frequent yet understudied postoperative total joint arthroplasty complication. This study aimed to describe cardiometabolic disease co-occurrence using latent class analysis, and associated postoperative AKI risk. METHODS: This retrospective analysis examined patients ≥18 years old undergoing primary total knee or hip arthroplasties within the US Multicenter Perioperative Outcomes Group of hospitals from 2008 to 2019. AKI was defined using modified Kidney Disease: Improving Global Outcomes (KDIGO) criteria. Latent classes were constructed from eight cardiometabolic diseases including hypertension, diabetes, and coronary artery disease, excluding obesity. A mixed-effects logistic regression model was constructed for the outcome of any AKI and the exposure of interaction between latent class and obesity status adjusting for preoperative and intraoperative covariates. RESULTS: Of 81 639 cases, 4007 (4.9%) developed AKI. Patients with AKI were more commonly older and non-Hispanic Black, with more significant comorbidity. A latent class model selected three groups of cardiometabolic patterning, labelled 'hypertension only' (n=37 223), 'metabolic syndrome (MetS)' (n=36 503), and 'MetS+cardiovascular disease (CVD)' (n=7913). After adjustment, latent class/obesity interaction groups had differential risk of AKI compared with those in 'hypertension only'/non-obese. Those 'hypertension only'/obese had 1.7-fold increased odds of AKI (95% confidence interval [CI]: 1.5-2.0). Compared with 'hypertension only'/non-obese, those 'MetS+CVD'/obese had the highest odds of AKI (odds ratio 3.1, 95% CI: 2.6-3.7), whereas 'MetS+CVD'/non-obese had 2.2 times the odds of AKI (95% CI: 1.8-2.7; model area under the curve 0.76). CONCLUSIONS: The risk of postoperative AKI varies widely between patients. The current study suggests that the co-occurrence of metabolic conditions (diabetes mellitus, hypertension), with or without obesity, is a more important risk factor for acute kidney injury than individual comorbid diseases.

Leveraging gene co-expression patterns to infer trait-relevant tissues in genome-wide association studies.

Genome-wide association studies (GWASs) have identified many SNPs associated with various common diseases. Understanding the biological functions of these identified SNP associations requires identifying disease/trait relevant tissues or cell types. Here, we develop a network method, CoCoNet, to facilitate the identification of trait-relevant tissues or cell types. Different from existing approaches, CoCoNet incorporates tissue-specific gene co-expression networks constructed from either bulk or single cell RNA sequencing (RNAseq) studies with GWAS data for trait-tissue inference. In particular, CoCoNet relies on a covariance regression network model to express gene-level effect measurements for the given GWAS trait as a function of the tissue-specific co-expression adjacency matrix. With a composite likelihood-based inference algorithm, CoCoNet is scalable to tens of thousands of genes. We validate the performance of CoCoNet through extensive simulations. We apply CoCoNet for an in-depth analysis of four neurological disorders and four autoimmune diseases, where we integrate the corresponding GWASs with bulk RNAseq data from 38 tissues and single cell RNAseq data from 10 cell types. In the real data applications, we show how CoCoNet can help identify specific glial cell types relevant for neurological disorders and identify disease-targeted colon tissues as relevant for autoimmune diseases.

Rare coding variants in RCN3 are associated with blood pressure.

BACKGROUND: While large genome-wide association studies have identified nearly one thousand loci associated with variation in blood pressure, rare variant identification is still a challenge. In family-based cohorts, genome-wide linkage scans have been successful in identifying rare genetic variants for blood pressure. This study aims to identify low frequency and rare genetic variants within previously reported linkage regions on chromosomes 1 and 19 in African American families from the Trans-Omics for Precision Medicine (TOPMed) program. Genetic association analyses weighted by linkage evidence were completed with whole genome sequencing data within and across TOPMed ancestral groups consisting of 60,388 individuals of European, African, East Asian, Hispanic, and Samoan ancestries. RESULTS: Associations of low frequency and rare variants in RCN3 and multiple other genes were observed for blood pressure traits in TOPMed samples. The association of low frequency and rare coding variants in RCN3 was further replicated in UK Biobank samples (N = 403,522), and reached genome-wide significance for diastolic blood pressure (p = 2.01 × 10- 7). CONCLUSIONS: Low frequency and rare variants in RCN3 contributes blood pressure variation. This study demonstrates that focusing association analyses in linkage regions greatly reduces multiple-testing burden and improves power to identify novel rare variants associated with blood pressure traits.