Recruitment of CTCF to the SIRT1 promoter after Oxidative Stress mediates Cardioprotective Transcription.

Because most DNA-binding transcription factors (dbTFs), including the architectural regulator CTCF, bind RNA and exhibit di-/multimerization, a central conundrum is whether these distinct properties are regulated post-transcriptionally to modulate transcriptional programs. Here, investigating stress-dependent activation of SIRT1, encoding an evolutionarily-conserved protein deacetylase, we show that induced phosphorylation of CTCF acts as a rheostat to permit CTCF occupancy of low-affinity promoter DNA sites to precisely the levels necessary. This CTCF recruitment to the SIRT1 promoter is eliciting a cardioprotective cardiomyocyte transcriptional activation program and provides resilience against the stress of the beating heart in vivo . Mice harboring a mutation in the conserved low-affinity CTCF promoter binding site exhibit an altered, cardiomyocyte-specific transcriptional program and a systolic heart failure phenotype. This transcriptional role for CTCF reveals that a covalent dbTF modification regulating signal-dependent transcription serves as a previously unsuspected component of the oxidative stress response.

Reconstructing parent genomes using siblings and other relatives.

Reconstructing the DNA of ancestors from their descendants has the potential to empower phenotypic analyses (including association and genetic nurture studies), improve pedigree reconstruction, and shed light on the ancestral population and phenotypes of ancestors. We developed HAPI-RECAP, a method that reconstructs the DNA of parents from full siblings and their relatives. This tool leverages HAPI2's output, a new phasing approach that applies to siblings (and optionally one or both parents) and reliably infers parent haplotypes but does not link the ungenotyped parents' DNA across chromosomes or between segments flanking ambiguities. By combining IBD between the reconstructed parents and the relatives, HAPI-RECAP resolves the source parent of these segments. Moreover, the method exploits crossovers the children inherited and sex-specific genetic maps to infer the reconstructed parents' sexes. We validated these methods on research participants from both 23andMe, Inc. and the San Antonio Mexican American Family Studies. Given data for one parent, HAPI2 reconstructs large fractions of the missing parent's DNA, between 77.6% and 99.97% among all families, and 90.3% on average in three- and four-child families. When reconstructing both parents, HAPI-RECAP inferred between 33.2% and 96.6% of the parents' genotypes, averaging 70.6% in four-child families. Reconstructed genotypes have average error rates < 10 -3 , or comparable to those from direct genotyping. HAPI-RECAP inferred the parent sexes 100% correctly given IBD-linked segments and can also reconstruct parents without any IBD. As datasets grow in size, more families will be implicitly collected; HAPI-RECAP holds promise to enable high quality parent genotype reconstruction.

Statistical Genetic Approaches to Investigate Genotype-by-Environment Interaction: Review and Novel Extension of Models.

Statistical genetic models of genotype-by-environment (G×E) interaction can be divided into two general classes, one on G×E interaction in response to dichotomous environments (e.g., sex, disease-affection status, or presence/absence of an exposure) and the other in response to continuous environments (e.g., physical activity, nutritional measurements, or continuous socioeconomic measures). Here we develop a novel model to jointly account for dichotomous and continuous environments. We develop the model in terms of a joint genotype-by-sex (for the dichotomous environment) and genotype-by-social determinants of health (SDoH; for the continuous environment). Using this model, we show how a depression variable, as measured by the Beck Depression Inventory-II survey instrument, is not only underlain by genetic effects (as has been reported elsewhere) but is also significantly determined by joint G×Sex and G×SDoH interaction effects. This model has numerous applications leading to potentially transformative research on the genetic and environmental determinants underlying complex diseases.

Role of Neurocellular Endoplasmic Reticulum Stress Response in Alzheimer's Disease and Related Dementias Risk.

Currently, more than 55 million people around the world suffer from dementia, and Alzheimer's Disease and Related Dementias (ADRD) accounts for nearly 60-70% of all those cases. The spread of Alzheimer's Disease (AD) pathology and progressive neurodegeneration in the hippocampus and cerebral cortex is strongly correlated with cognitive decline in AD patients; however, the molecular underpinning of ADRD's causality is still unclear. Studies of postmortem AD brains and animal models of AD suggest that elevated endoplasmic reticulum (ER) stress may have a role in ADRD pathology through altered neurocellular homeostasis in brain regions associated with learning and memory. To study the ER stress-associated neurocellular response and its effects on neurocellular homeostasis and neurogenesis, we modeled an ER stress challenge using thapsigargin (TG), a specific inhibitor of sarco/endoplasmic reticulum Ca2+ ATPase (SERCA), in the induced pluripotent stem cell (iPSC)-derived neural stem cells (NSCs) of two individuals from our Mexican American Family Study (MAFS). High-content screening and transcriptomic analysis of the control and ER stress-challenged NSCs showed that the NSCs' ER stress response resulted in a significant decline in NSC self-renewal and an increase in apoptosis and cellular oxidative stress. A total of 2300 genes were significantly (moderated t statistics FDR-corrected p-value ≤ 0.05 and fold change absolute ≥ 2.0) differentially expressed (DE). The pathway enrichment and gene network analysis of DE genes suggests that all three unfolded protein response (UPR) pathways, protein kinase RNA-like ER kinase (PERK), activating transcription factor-6 (ATF-6), and inositol-requiring enzyme-1 (IRE1), were significantly activated and cooperatively regulated the NSCs' transcriptional response to ER stress. Our results show that IRE1/X-box binding protein 1 (XBP1) mediated transcriptional regulation of the E2F transcription factor 1 (E2F1) gene, and its downstream targets have a dominant role in inducing G1/S-phase cell cycle arrest in ER stress-challenged NSCs. The ER stress-challenged NSCs also showed the activation of C/EBP homologous protein (CHOP)-mediated apoptosis and the dysregulation of synaptic plasticity and neurotransmitter homeostasis-associated genes. Overall, our results suggest that the ER stress-associated attenuation of NSC self-renewal, increased apoptosis, and dysregulated synaptic plasticity and neurotransmitter homeostasis plausibly play a role in the causation of ADRD.

Validation of human telomere length multi-ancestry meta-analysis association signals identifies POP5 and KBTBD6 as human telomere length regulation genes.

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

Pre-Infection Innate Immunity Attenuates SARS-CoV-2 Infection and Viral Load in iPSC-Derived Alveolar Epithelial Type 2 Cells.

A large portion of the heterogeneity in coronavirus disease 2019 (COVID-19) susceptibility and severity of illness (SOI) remains poorly understood. Recent evidence suggests that SARS-CoV-2 infection-associated damage to alveolar epithelial type 2 cells (AT2s) in the distal lung may directly contribute to disease severity and poor prognosis in COVID-19 patients. Our in vitro modeling of SARS-CoV-2 infection in induced pluripotent stem cell (iPSC)-derived AT2s from 10 different individuals showed interindividual variability in infection susceptibility and the postinfection cellular viral load. To understand the underlying mechanism of the AT2's capacity to regulate SARS-CoV-2 infection and cellular viral load, a genome-wide differential gene expression analysis between the mock and SARS-CoV-2 infection-challenged AT2s was performed. The 1393 genes, which were significantly (one-way ANOVA FDR-corrected p ≤ 0.05; FC abs ≥ 2.0) differentially expressed (DE), suggest significant upregulation of viral infection-related cellular innate immune response pathways (p-value ≤ 0.05; activation z-score ≥ 3.5), and significant downregulation of the cholesterol- and xenobiotic-related metabolic pathways (p-value ≤ 0.05; activation z-score ≤ -3.5). Whilst the effect of post-SARS-CoV-2 infection response on the infection susceptibility and postinfection viral load in AT2s is not clear, interestingly, pre-infection (mock-challenged) expression of 238 DE genes showed a high correlation with the postinfection SARS-CoV-2 viral load (FDR-corrected p-value ≤ 0.05 and r2-absolute ≥ 0.57). The 85 genes whose expression was negatively correlated with the viral load showed significant enrichment in viral recognition and cytokine-mediated innate immune GO biological processes (p-value range: 4.65 × 10-10 to 2.24 × 10-6). The 153 genes whose expression was positively correlated with the viral load showed significant enrichment in cholesterol homeostasis, extracellular matrix, and MAPK/ERK pathway-related GO biological processes (p-value range: 5.06 × 10-5 to 6.53 × 10-4). Overall, our results strongly suggest that AT2s' pre-infection innate immunity and metabolic state affect their susceptibility to SARS-CoV-2 infection and viral load.

Genetic correlations between liver fat content, metabolic health, and adiposity distribution in the Fels Longitudinal Study.

BACKGROUND AND AIMS: Hepatic steatosis is known to be heritable, but its genetic basis is mostly uncharacterized. Steatosis is associated with metabolic and adiposity features; recent studies hypothesize that shared genetic effects between these traits could account for some of the unexplained heritability. This study aimed to quantify these genetic associations in a family-based sample of non-Hispanic white adults. METHODS AND RESULTS: 704 participants (18-95 years, 55.8% female) from the Fels Longitudinal Study with an MRI assessment of liver fat were included. Quantitative genetic analyses estimated the age- and sex-adjusted heritability of individual traits and the genetic correlations within trait pairs. Mean liver fat was 5.95% (SE = 0.23) and steatosis (liver fat >5.56%) was present in 29.8% of participants. Heritability (h2± SE) of steatosis was 0.72 ± 0.17 (p = 6.80e-6). All other traits including liver enzymes, fasting glucose, HOMA-IR, visceral and subcutaneous adipose tissue (VAT, SAT), body mass index, body fat percent, waist circumference, lipids and blood pressure were also heritable. Significant genetic correlations were found between liver fat and all traits except aspartate aminotransferase (AST), and among most trait pairs. Highest genetic correlations were between liver fat and HOMA-IR (0.85 ± 0.08, p = 1.73e-8), fasting glucose and ALT (0.89 ± 0.26, p = 6.68e-5), and HOMA-IR with: waist circumference (0.81 ± 0.12, p = 3.76e-6), body fat percent (0.78 ± 0.12 p = 2.42e-5) and VAT (0.73 ± 0.07, p = 6.37e-8). CONCLUSIONS: Common genes may exist between liver fat accumulation, metabolic features and adiposity phenotypes.

Lifespan reference curves for harmonizing multi-site regional brain white matter metrics from diffusion MRI.

Age-related white matter (WM) microstructure maturation and decline occur throughout the human lifespan, complementing the process of gray matter development and degeneration. Here, we create normative lifespan reference curves for global and regional WM microstructure by harmonizing diffusion MRI (dMRI)-derived data from ten public datasets (N = 40,898 subjects; age: 3-95 years; 47.6% male). We tested three harmonization methods on regional diffusion tensor imaging (DTI) based fractional anisotropy (FA), a metric of WM microstructure, extracted using the ENIGMA-DTI pipeline. ComBat-GAM harmonization provided multi-study trajectories most consistent with known WM maturation peaks. Lifespan FA reference curves were validated with test-retest data and used to assess the effect of the ApoE4 risk factor for dementia in WM across the lifespan. We found significant associations between ApoE4 and FA in WM regions associated with neurodegenerative disease even in healthy individuals across the lifespan, with regional age-by-genotype interactions. Our lifespan reference curves and tools to harmonize new dMRI data to the curves are publicly available as eHarmonize (https://github.com/ahzhu/eharmonize).

Metabolic syndrome traits exhibit genotype-by-environment interaction in relation to socioeconomic status in the Mexican American family heart study.

Background: Socioeconomic Status (SES) is a potent environmental determinant of health. To our knowledge, no assessment of genotype-environment interaction has been conducted to consider the joint effects of socioeconomic status and genetics on risk for metabolic disease. We analyzed data from the Mexican American Family Studies (MAFS) to evaluate the hypothesis that genotype-by-environment interaction (GxE) is an essential determinant of variation in risk factors for metabolic syndrome (MS). Methods: We employed a maximum likelihood estimation of the decomposition of variance components to detect GxE interaction. After excluding individuals with diabetes and individuals on medication for diabetes, hypertension, or dyslipidemia, we analyzed 12 MS risk factors: fasting glucose (FG), fasting insulin (FI), 2-h glucose (2G), 2-h insulin (2I), body mass index (BMI), waist circumference (WC), leptin (LP), high-density lipoprotein-cholesterol (HDL-C), triglycerides (TG), total serum cholesterol (TSC), systolic blood pressure (SBP), and diastolic blood pressure (DBP). Our SES variable used a combined score of Duncan's socioeconomic index and education years. Heterogeneity in the additive genetic variance across the SES continuum and a departure from unity in the genetic correlation coefficient were taken as evidence of GxE interaction. Hypothesis tests were conducted using standard likelihood ratio tests. Results: We found evidence of GxE for fasting glucose, 2-h glucose, 2-h insulin, BMI, and triglycerides. The genetic effects underlying the insulin/glucose metabolism component of MS are upregulated at the lower end of the SES spectrum. We also determined that the household variance for systolic blood pressure decreased with increasing SES. Conclusion: These results show a significant change in the GxE interaction underlying the major components of MS in response to changes in socioeconomic status. Further mRNA sequencing studies will identify genes and canonical gene pathways to support our molecular-level hypotheses.

Imputation of plasma lipid species to facilitate integration of lipidomic datasets.

Recent advancements in plasma lipidomic profiling methodology have significantly increased specificity and accuracy of lipid measurements. This evolution, driven by improved chromatographic and mass spectrometric resolution of newer platforms, has made it challenging to align datasets created at different times, or on different platforms. Here we present a framework for harmonising such plasma lipidomic datasets with different levels of granularity in their lipid measurements. Our method utilises elastic-net prediction models, constructed from high-resolution lipidomics reference datasets, to predict unmeasured lipid species in lower-resolution studies. The approach involves (1) constructing composite lipid measures in the reference dataset that map to less resolved lipids in the target dataset, (2) addressing discrepancies between aligned lipid species, (3) generating prediction models, (4) assessing their transferability into the targe dataset, and (5) evaluating their prediction accuracy. To demonstrate our approach, we used the AusDiab population-based cohort (747 lipid species) as the reference to impute unmeasured lipid species into the LIPID study (342 lipid species). Furthermore, we compared measured and imputed lipids in terms of parameter estimation and predictive performance, and validated imputations in an independent study. Our method for harmonising plasma lipidomic datasets will facilitate model validation and data integration efforts.

Reproducibility and Repeatability of Intravoxel Incoherent Motion MRI Acquisition Methods in Liver.

BACKGROUND: Intravoxel incoherent motion (IVIM) diffusion weighted MRI (DWI) has potential for evaluating hepatic fibrosis but image acquisition technique influence on diffusion parameter estimation bears investigation. PURPOSE: To minimize variability and maximize repeatably in abdominal DWI in terms of IVIM parameter estimates. STUDY TYPE: Prospective test-retest and image quality comparison. SUBJECTS: Healthy volunteers (3F/7M, 29.9 ± 12.9 years) and Family Study subjects (18F/12M, 51.7 ± 16.7 years), without and with liver steatosis. FIELD STRENGTH/SEQUENCE: Abdominal single-shot echo-planar imaging (EPI) and simultaneous multi-slice (SMS) DWI sequences with respiratory triggering (RT), breath-holding (BH), and navigator echo (NE) at 3 Tesla. ASSESSMENT: SMS-BH, EPI-NE, and SMS-RT data from twice-scanned healthy volunteers were analyzed using 6 × b-values (0-800 s⋅mm-2 ) and lower (LO) and higher (HI) b-value ranges. Family Study subjects were scanned using SMS and standard EPI sequences. The biexponential IVIM model was used to estimate fast-diffusion coefficient (Df ), fraction of fast diffusion (f), and slow-diffusion coefficient (Ds ). Scan time, estimated signal-to-noise ratio (eSNR), eSNR per acquisition, and distortion ratio were compared. STATISTICAL TESTS: Coefficients of variation (CoV) and Bland Altman analyses were performed for test-retest repeatability. Interclass correlation coefficient (ICC) assessed interobserver agreement with P < 0.05 deemed significant. RESULTS: Within-subject CoVs among volunteers (N = 10) for f and Ds were lowest in EPI-NE-LO (11.6%) and SMS-RT-HI (11.1%). Inter-observer ICCs for f and Ds were highest for EPI-NE-LO (0.63) and SMS-RT-LO (0.76). Df could not be estimated for most subjects. Estimated eSNR (EPI = 21.9, SMS = 4.7) and eSNR time (EPI = 6.7, SMS = 16.6) were greater for SMS, with less distortion in the liver region (DR-PE: EPI = 23.6, SMS = 13.1). DATA CONCLUSION: Simultaneous multislice acquisitions had significantly less variability and higher ICCs of Ds , higher eSNR, less distortion, and reduced scan time compared to EPI. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 1.

Riemannian frameworks for the harmonization of resting-state functional MRI scans.

Magnetic Resonance Imaging provides unprecedented images of the brain. Unfortunately, scanners and acquisition protocols can significantly impact MRI scans. The development of statistical methods able to reduce this variability without altering the relevant information in the scans, often coined harmonization methods, has been the topic of an increasing research effort supported by the recent growth of publicly available neuroimaging data sets and new possibilities for combining them to achieve greater statistical power. In this work, we focus on the challenges specifically raised by the harmonization of resting-state functional MRI scans. We propose to harmonize resting-state fMRI scans by reducing the impact of covariates such as scanner differences and scanning protocols on their associated functional connectomes and then propagating the changes back to the rs-fMRI time series. We use Riemannian geometric frameworks to preserve the mathematical properties of functional connectomes during their harmonization, and we demonstrate how state-of-the-art harmonization methods can be embedded within these frameworks to reduce covariates effects while preserving the relevant clinical information associated with aging or brain disorders. During our experiments, a large set of synthetic data was generated and processed to compare eighty variants of the proposed approach. The framework achieving the best harmonization was then applied to three low-dimensional data sets made of 712 sets of fMRI time series provided by the ABIDE consortium and two high-dimensional data sets obtained by processing 1527 rs-fMRI scans provided by the Human Connectome Project, the Framingham Heart Study and the Genetics of Brain Structure and Function study. These experiments established that our new framework could successfully harmonize low-dimensional connectomes and voxelwise functional time series and confirmed the need for preserving connectomes properties during their harmonization.

Functional characterization of the disease-associated CCL2 rs1024611G-rs13900T haplotype: The role of the RNA-binding protein HuR.

CC-chemokine ligand 2 (CCL2) is involved in the pathogenesis of several diseases associated with monocyte/macrophage recruitment, such as HIV-associated neurocognitive disorder (HAND), tuberculosis, and atherosclerosis. The rs1024611 (alleles:A>G; G is the risk allele) polymorphism in the CCL2 cis-regulatory region is associated with increased CCL2 expression in vitro and ex vivo, leukocyte mobilization in vivo, and deleterious disease outcomes. However, the molecular basis for the rs1024611-associated differential CCL2 expression remains poorly characterized. It is conceivable that genetic variant(s) in linkage disequilibrium (LD) with rs1024611 could mediate such effects. Previously, we used rs13900 (alleles:_C>T) in the CCL2 3' untranslated region (3' UTR) that is in perfect LD with rs1024611 to demonstrate allelic expression imbalance (AEI) of CCL2 in heterozygous individuals. Here we tested the hypothesis that the rs13900 could modulate CCL2 expression by altering mRNA turnover and/or translatability. The rs13900 T allele conferred greater stability to the CCL2 transcript when compared to the rs13900 C allele. The rs13900 T allele also had increased binding to Human Antigen R (HuR), an RNA-binding protein, in vitro and ex vivo. The rs13900 alleles imparted differential activity to reporter vectors and influenced the translatability of the reporter transcript. We further demonstrated a role for HuR in mediating allele-specific effects on CCL2 expression in overexpression and silencing studies. The presence of the rs1024611G-rs13900T conferred a distinct transcriptomic signature related to inflammation and immunity. Our studies suggest that the differential interactions of HuR with rs13900 could modulate CCL2 expression and explain the interindividual differences in CCL2-mediated disease susceptibility.

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

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

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis.

BACKGROUND: Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D. METHODS: We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test. RESULTS: Using a Bonferroni-corrected significance threshold of P<1.6×10-4, we identified 3 genes (ATP1B1, ARVCF, and LIPG) associated with CAC and 2 genes (ABCG8 and EIF2B2) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both ATP1B1 and ARVCF also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis. CONCLUSIONS: These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.

Metabolic phenotyping of BMI to characterize cardiometabolic risk: evidence from large population-based cohorts.

Obesity is a risk factor for type 2 diabetes and cardiovascular disease. However, a substantial proportion of patients with these conditions have a seemingly normal body mass index (BMI). Conversely, not all obese individuals present with metabolic disorders giving rise to the concept of "metabolically healthy obese". We use lipidomic-based models for BMI to calculate a metabolic BMI score (mBMI) as a measure of metabolic dysregulation associated with obesity. Using the difference between mBMI and BMI (mBMIΔ), we identify individuals with a similar BMI but differing in their metabolic health and disease risk profiles. Exercise and diet associate with mBMIΔ suggesting the ability to modify mBMI with lifestyle intervention. Our findings show that, the mBMI score captures information on metabolic dysregulation that is independent of the measured BMI and so provides an opportunity to assess metabolic health to identify "at risk" individuals for targeted intervention and monitoring.

A Scan of Pleiotropic Immune Mediated Disease Genes Identifies Novel Determinants of Baseline FVIII Inhibitor Status in Hemophilia-A.

Hemophilia-A (HA) is caused by heterogeneous loss-of-function factor (F)VIII gene (F8)-mutations and deficiencies in plasma-FVIII-activity that impair intrinsic-pathway-mediated coagulation-amplification. The standard-of-care for severe-HA-patients is regular infusions of therapeutic-FVIII-proteins (tFVIIIs) but ~30% develop neutralizing-tFVIII-antibodies called "FVIII-inhibitors (FEIs)" and become refractory. We used the PATH study and ImmunoChip to scan immune-mediated-disease (IMD)-genes for novel and/or replicated genomic-sequence-variations associated with baseline-FEI-status while accounting for non-independence of data due to genetic-relatedness and F8-mutational-heterogeneity. The baseline-FEI-status of 450 North American PATH subjects-206 with black-African-ancestry and 244 with white-European-ancestry-was the dependent variable. The F8-mutation-data and a genetic-relatedness matrix were incorporated into a binary linear-mixed model of genetic association with baseline-FEI-status. We adopted a gene-centric-association-strategy to scan, as candidates, pleiotropic-IMD-genes implicated in the development of either ³2 autoimmune-/autoinflammatory-disorders (AADs) or ³1 AAD and FEIs. Baseline-FEI-status was significantly associated with SNPs assigned to NOS2A (rs117382854; p=3.2E-6) and B3GNT2 (rs10176009; p=5.1E-6), which have functions in anti-microbial-/-tumoral-immunity. Among IMD-genes implicated in FEI-risk previously, we identified strong associations with CTLA4 assigned SNPs (p=2.2E-5). The F8-mutation-effect underlies ~15% of the total heritability for baseline-FEI-status. Additive genetic heritability and SNPs in IMD-genes account for >50% of the patient-specific variability in baseline-FEI-status. Race is a significant determinant independent of F8-mutation-effects and non-F8-genetics.

Two neurostructural subtypes: results of machine learning on brain images from 4,291 individuals with schizophrenia.

Machine learning can be used to define subtypes of psychiatric conditions based on shared clinical and biological foundations, presenting a crucial step toward establishing biologically based subtypes of mental disorders. With the goal of identifying subtypes of disease progression in schizophrenia, here we analyzed cross-sectional brain structural magnetic resonance imaging (MRI) data from 4,291 individuals with schizophrenia (1,709 females, age=32.5 years±11.9) and 7,078 healthy controls (3,461 females, age=33.0 years±12.7) pooled across 41 international cohorts from the ENIGMA Schizophrenia Working Group, non-ENIGMA cohorts and public datasets. Using a machine learning approach known as Subtype and Stage Inference (SuStaIn), we implemented a brain imaging-driven classification that identifies two distinct neurostructural subgroups by mapping the spatial and temporal trajectory of gray matter (GM) loss in schizophrenia. Subgroup 1 (n=2,622) was characterized by an early cortical-predominant loss (ECL) with enlarged striatum, whereas subgroup 2 (n=1,600) displayed an early subcortical-predominant loss (ESL) in the hippocampus, amygdala, thalamus, brain stem and striatum. These reconstructed trajectories suggest that the GM volume reduction originates in the Broca's area/adjacent fronto-insular cortex for ECL and in the hippocampus/adjacent medial temporal structures for ESL. With longer disease duration, the ECL subtype exhibited a gradual worsening of negative symptoms and depression/anxiety, and less of a decline in positive symptoms. We confirmed the reproducibility of these imaging-based subtypes across various sample sites, independent of macroeconomic and ethnic factors that differed across these geographic locations, which include Europe, North America and East Asia. These findings underscore the presence of distinct pathobiological foundations underlying schizophrenia. This new imaging-based taxonomy holds the potential to identify a more homogeneous sub-population of individuals with shared neurobiological attributes, thereby suggesting the viability of redefining existing disorder constructs based on biological factors.

Genotype-by-socioeconomic status interaction influences heart disease risk scores and carotid artery thickness in Mexican Americans: the predominant role of education in comparison to household income and socioeconomic index.

Background: Socioeconomic status (SES) is a potent environmental determinant of health. To our knowledge, no assessment of genotype-environment interaction has been conducted to consider the joint effects of socioeconomic status and genetics on risk for cardiovascular disease (CVD). We analyzed Mexican American Family Studies (MAFS) data to evaluate the hypothesis that genotype-by-environment interaction (GxE) is an important determinant of variation in CVD risk factors. Methods: We employed a linear mixed model to investigate GxE in Mexican American extended families. We studied two proxies for CVD [Pooled Cohort Equation Risk Scores/Framingham Risk Scores (FRS/PCRS) and carotid artery intima-media thickness (CA-IMT)] in relation to socioeconomic status as determined by Duncan's Socioeconomic Index (SEI), years of education, and household income. Results: We calculated heritability for FRS/PCRS and carotid artery intima-media thickness. There was evidence of GxE due to additive genetic variance heterogeneity and genetic correlation for FRS, PCRS, and CA-IMT measures for education (environment) but not for household income or SEI. Conclusion: The genetic effects underlying CVD are dynamically modulated at the lower end of the SES spectrum. There is a significant change in the genetic architecture underlying the major components of CVD in response to changes in education.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.