Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations.

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci.

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

Electrophysiological Correlates of Dentate Nucleus Deep Brain Stimulation for Poststroke Motor Recovery.

While ipsilesional cortical electroencephalography has been associated with poststroke recovery mechanisms and outcomes, the role of the cerebellum and its interaction with the ipsilesional cortex is still largely unknown. We have previously shown that poststroke motor control relies on increased corticocerebellar coherence (CCC) in the low beta band to maintain motor task accuracy and to compensate for decreased excitability of the ipsilesional cortex. We now extend our work to investigate corticocerebellar network changes associated with chronic stimulation of the dentato-thalamo-cortical pathway aimed at promoting poststroke motor rehabilitation. We investigated the excitability of the ipsilesional cortex, the dentate (DN), and their interaction as a function of treatment outcome measures. Relative to baseline, 10 human participants (two women) at the end of 4-8 months of DN deep brain stimulation (DBS) showed (1) significantly improved motor control indexed by computerized motor tasks; (2) significant increase in ipsilesional premotor cortex event-related desynchronization that correlated with improvements in motor function; and (3) significant decrease in CCC, including causal interactions between the DN and ipsilesional cortex, which also correlated with motor function improvements. Furthermore, we show that the functional state of the DN in the poststroke state and its connectivity with the ipsilesional cortex were predictive of motor outcomes associated with DN-DBS. The findings suggest that as participants recovered, the ipsilesional cortex became more involved in motor control, with less demand on the cerebellum to support task planning and execution. Our data provide unique mechanistic insights into the functional state of corticocerebellar-cortical network after stroke and its modulation by DN-DBS.

Determinants of electrical propagation and propagation block in Arrhythmogenic Cardiomyopathy.

Gap junction and ion channel remodeling occur early in Arrhythmogenic Cardiomyopathy (ACM), but their pathogenic consequences have not been elucidated. Here, we identified the arrhythmogenic substrate, consisting of propagation slowing and conduction block, in ACM models expressing two different desmosomal gene variants. Neonatal rat ventricular myocytes were transduced to express variants in genes encoding desmosomal proteins plakoglobin or plakophilin-2. Studies were performed in engineered cells and anisotropic tissues to quantify changes in conduction velocity, formation of unidirectional propagation, cell-cell electrical coupling, and ion currents. Conduction velocity decreased by 71% and 63% in the two ACM models. SB216763, an inhibitor of glycogen synthase kinase-3 beta, restored conduction velocity to near normal levels. Compared to control, both ACM models showed greater propensity for unidirectional conduction block, which increased further at greater stimulation frequencies. Cell-cell electrical conductance measured in cell pairs was reduced by 86% and 87% in the two ACM models. Computer modeling showed close correspondence between simulated and experimentally determined changes in conduction velocity. The simulation identified that reduced cell-cell electrical coupling was the dominant factor leading to slow conduction, while the combination of reduced cell-cell electrical coupling, reduced sodium current and inward rectifier potassium current explained the development of unidirectional block. Expression of two different ACM variants markedly reduced cell-cell electrical coupling and conduction velocity, and greatly increased the likelihood of developing unidirectional block - both key features of arrhythmogenesis. This study provides the first quantitative analysis of cellular electrophysiological changes leading to the substrate of reentrant arrhythmias in early stage ACM.

Hypoxia-inducible factor-driven glycolytic adaptations in host-microbe interactions.

Mammalian cells utilize glucose as a primary carbon source to produce energy for most cellular functions. However, the bioenergetic homeostasis of cells can be perturbed by environmental alterations, such as changes in oxygen levels which can be associated with bacterial infection. Reduction in oxygen availability leads to a state of hypoxia, inducing numerous cellular responses that aim to combat this stress. Importantly, hypoxia strongly augments cellular glycolysis in most cell types to compensate for the loss of aerobic respiration. Understanding how this host cell metabolic adaptation to hypoxia impacts the course of bacterial infection will identify new anti-microbial targets. This review will highlight developments in our understanding of glycolytic substrate channeling and spatiotemporal enzymatic organization in response to hypoxia, shedding light on the integral role of the hypoxia-inducible factor (HIF) during host-pathogen interactions. Furthermore, the ability of intracellular and extracellular bacteria (pathogens and commensals alike) to modulate host cellular glucose metabolism will be discussed.

A novel restrainer device for acquistion of brain images in awake rats.

Functional neuroimaging methods like fMRI and PET are vital in neuroscience research, but require that subjects remain still throughout the scan. In animal research, anesthetic agents are typically applied to facilitate the acquisition of high-quality data with minimal motion artifact. However, anesthesia can have profound effects on brain metabolism, selectively altering dynamic neural networks and confounding the acquired data. To overcome the challenge, we have developed a novel head fixation device designed to support awake rat brain imaging. A validation experiment demonstrated that the device effectively minimizes animal motion throughout the scan, with mean absolute displacement and mean relative displacement of 0.0256 (SD: 0.001) and 0.009 (SD: 0.002), across eight evaluated subjects throughout fMRI image acquisition (total scanning time per subject: 31 min, 12 s). Furthermore, the awake scans did not induce discernable stress to the animals, with stable physiological parameters throughout the scan (Mean HR: 344, Mean RR: 56, Mean SpO2: 94 %) and unaltered serum corticosterone levels (p = 0.159). In conclusion, the device presented in this paper offers an effective and safe method of acquiring functional brain images in rats, allowing researchers to minimize the confounding effects of anesthetic use.

Fat metabolism is associated with telomere length in six population-based studies.

Telomeres are repetitive DNA sequences located at the end of chromosomes, which are associated to biological aging, cardiovascular disease, cancer and mortality. Lipid and fatty acid metabolism have been associated with telomere shortening. We have conducted an in-depth study investigating the association of metabolic biomarkers with telomere length (LTL). We performed an association analysis of 226 metabolic biomarkers with LTL using data from 11 775 individuals from six independent population-based cohorts (BBMRI-NL consortium). Metabolic biomarkers include lipoprotein lipids and subclasses, fatty acids, amino acids, glycolysis measures and ketone bodies. LTL was measured by quantitative polymerase chain reaction or FlowFISH. Linear regression analysis was performed adjusting for age, sex, lipid-lowering medication and cohort-specific covariates (model 1) and additionally for body mass index (BMI) and smoking (model 2), followed by inverse variance-weighted meta-analyses (significance threshold Pmeta = 6.5 × 10-4). We identified four metabolic biomarkers positively associated with LTL, including two cholesterol to lipid ratios in small VLDL (S-VLDL-C % and S-VLDL-CE %) and two omega-6 fatty acid ratios (FAw6/FA and LA/FA). After additionally adjusting for BMI and smoking, these metabolic biomarkers remained associated with LTL with similar effect estimates. In addition, cholesterol esters in very small VLDL (XS-VLDL-CE) became significantly associated with LTL (P = 3.6 × 10-4). We replicated the association of FAw6/FA with LTL in an independent dataset of 7845 individuals (P = 1.9 × 10-4). To conclude, we identified multiple metabolic biomarkers involved in lipid and fatty acid metabolism that may be involved in LTL biology. Longitudinal studies are needed to exclude reversed causation.

Whole-exome sequencing of 14 389 individuals from the ESP and CHARGE consortia identifies novel rare variation associated with hemostatic factors.

Plasma levels of fibrinogen, coagulation factors VII and VIII and von Willebrand factor (vWF) are four intermediate phenotypes that are heritable and have been associated with the risk of clinical thrombotic events. To identify rare and low-frequency variants associated with these hemostatic factors, we conducted whole-exome sequencing in 10 860 individuals of European ancestry (EA) and 3529 African Americans (AAs) from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and the National Heart, Lung and Blood Institute's Exome Sequencing Project. Gene-based tests demonstrated significant associations with rare variation (minor allele frequency < 5%) in fibrinogen gamma chain (FGG) (with fibrinogen, P = 9.1 × 10-13), coagulation factor VII (F7) (with factor VII, P = 1.3 × 10-72; seven novel variants) and VWF (with factor VIII and vWF; P = 3.2 × 10-14; one novel variant). These eight novel rare variant associations were independent of the known common variants at these loci and tended to have much larger effect sizes. In addition, one of the rare novel variants in F7 was significantly associated with an increased risk of venous thromboembolism in AAs (Ile200Ser; rs141219108; P = 4.2 × 10-5). After restricting gene-based analyses to only loss-of-function variants, a novel significant association was detected and replicated between factor VIII levels and a stop-gain mutation exclusive to AAs (rs3211938) in CD36 molecule (CD36). This variant has previously been linked to dyslipidemia but not with the levels of a hemostatic factor. These efforts represent the largest integration of whole-exome sequence data from two national projects to identify genetic variation associated with plasma hemostatic factors.

Genetic and clinical determinants of abdominal aortic diameter: genome-wide association studies, exome array data and Mendelian randomization study.

Progressive dilation of the infrarenal aortic diameter is a consequence of the ageing process and is considered the main determinant of abdominal aortic aneurysm (AAA). We aimed to investigate the genetic and clinical determinants of abdominal aortic diameter (AAD). We conducted a meta-analysis of genome-wide association studies in 10 cohorts (n = 13 542) imputed to the 1000 Genome Project reference panel including 12 815 subjects in the discovery phase and 727 subjects [Partners Biobank cohort 1 (PBIO)] as replication. Maximum anterior-posterior diameter of the infrarenal aorta was used as AAD. We also included exome array data (n = 14 480) from seven epidemiologic studies. Single-variant and gene-based associations were done using SeqMeta package. A Mendelian randomization analysis was applied to investigate the causal effect of a number of clinical risk factors on AAD. In genome-wide association study (GWAS) on AAD, rs74448815 in the intronic region of LDLRAD4 reached genome-wide significance (beta = -0.02, SE = 0.004, P-value = 2.10 × 10-8). The association replicated in the PBIO1 cohort (P-value = 8.19 × 10-4). In exome-array single-variant analysis (P-value threshold = 9 × 10-7), the lowest P-value was found for rs239259 located in SLC22A20 (beta = 0.007, P-value = 1.2 × 10-5). In the gene-based analysis (P-value threshold = 1.85 × 10-6), PCSK5 showed an association with AAD (P-value = 8.03 × 10-7). Furthermore, in Mendelian randomization analyses, we found evidence for genetic association of pulse pressure (beta = -0.003, P-value = 0.02), triglycerides (beta = -0.16, P-value = 0.008) and height (beta = 0.03, P-value < 0.0001), known risk factors for AAA, consistent with a causal association with AAD. Our findings point to new biology as well as highlighting gene regions in mechanisms that have previously been implicated in the genetics of other vascular diseases.