Red blood cell (RBC) traits are important heritable clinical biomarkers and modifiers of disease severity. To identify coding genetic variants associated with these traits, we conducted meta-analyses of seven RBC phenotypes in 130,273 multi-ethnic individuals from studies genotyped on an exome array. After conditional analyses and replication in 27,480 independent individuals, we identified 16 new RBC variants. We found low-frequency missense variants in MAP1A (rs55707100, minor allele frequency [MAF] = 3.3%, p = 2 × 10(-10) for hemoglobin [HGB]) and HNF4A (rs1800961, MAF = 2.4%, p < 3 × 10(-8) for hematocrit [HCT] and HGB). In African Americans, we identified a nonsense variant in CD36 associated with higher RBC distribution width (rs3211938, MAF = 8.7%, p = 7 × 10(-11)) and showed that it is associated with lower CD36 expression and strong allelic imbalance in ex vivo differentiated human erythroblasts. We also identified a rare missense variant in ALAS2 (rs201062903, MAF = 0.2%) associated with lower mean corpuscular volume and mean corpuscular hemoglobin (p < 8 × 10(-9)). Mendelian mutations in ALAS2 are a cause of sideroblastic anemia and erythropoietic protoporphyria. Gene-based testing highlighted three rare missense variants in PKLR, a gene mutated in Mendelian non-spherocytic hemolytic anemia, associated with HGB and HCT (SKAT p < 8 × 10(-7)). These rare, low-frequency, and common RBC variants showed pleiotropy, being also associated with platelet, white blood cell, and lipid traits. Our association results and functional annotation suggest the involvement of new genes in human erythropoiesis. We also confirm that rare and low-frequency variants play a role in the architecture of complex human traits, although their phenotypic effect is generally smaller than originally anticipated.

We conducted a genome-wide association study (GWAS) on multiple sclerosis (MS) susceptibility in German cohorts with 4888 cases and 10,395 controls. In addition to associations within the major histocompatibility complex (MHC) region, 15 non-MHC loci reached genome-wide significance. Four of these loci are novel MS susceptibility loci. They map to the genes L3MBTL3, MAZ, ERG, and SHMT1. The lead variant at SHMT1 was replicated in an independent Sardinian cohort. Products of the genes L3MBTL3, MAZ, and ERG play important roles in immune cell regulation. SHMT1 encodes a serine hydroxymethyltransferase catalyzing the transfer of a carbon unit to the folate cycle. This reaction is required for regulation of methylation homeostasis, which is important for establishment and maintenance of epigenetic signatures. Our GWAS approach in a defined population with limited genetic substructure detected associations not found in larger, more heterogeneous cohorts, thus providing new clues regarding MS pathogenesis.

The human leukocyte antigen (HLA) haplotype DRB1*15:01 is the major risk factor for multiple sclerosis (MS). Here, we find that DRB1*15:01 is hypomethylated and predominantly expressed in monocytes among carriers of DRB1*15:01. A differentially methylated region (DMR) encompassing HLA-DRB1 exon 2 is particularly affected and displays methylation-sensitive regulatory properties in vitro. Causal inference and Mendelian randomization provide evidence that HLA variants mediate risk for MS via changes in the HLA-DRB1 DMR that modify HLA-DRB1 expression. Meta-analysis of 14,259 cases and 171,347 controls confirms that these variants confer risk from DRB1*15:01 and also identifies a protective variant (rs9267649, p < 3.32 × 10-8, odds ratio = 0.86) after conditioning for all MS-associated variants in the region. rs9267649 is associated with increased DNA methylation at the HLA-DRB1 DMR and reduced expression of HLA-DRB1, suggesting a modulation of the DRB1*15:01 effect. Our integrative approach provides insights into the molecular mechanisms of MS susceptibility and suggests putative therapeutic strategies targeting a methylation-mediated regulation of the major risk gene.

In recent years, human microbiota, especially gut microbiota, have emerged as an important yet complex trait influencing human metabolism, immunology, and diseases. Many studies are investigating the forces underlying the observed variation, including the human genetic variants that shape human microbiota. Several preliminary genome-wide association studies (GWAS) have been completed, but more are necessary to achieve a fuller picture.

The transcriptome of non-coding RNA (ncRNA) species is increasingly focused in Alzheimer's disease (AD) research. NcRNAs comprise, among others, transfer RNAs, long non-coding RNAs and microRNAs (miRs), each with their own specific biological function. We used smallRNASeq to assess miR expression in the hippocampus of young (3 month old) and aged (8 month old) Tg4-42 mice, a model system for sporadic AD, as well as age-matched wildtype controls. Tg4-42 mice express N-truncated Aβ4-42, develop age-related neuron loss, reduced neurogenesis and behavioral deficits. Our results do not only confirm known miR-AD associations in Tg4-42 mice, but more importantly pinpoint 22 additional miRs associated to the disease. Twenty-five miRs were differentially expressed in both aged Tg4-42 and aged wildtype mice while eight miRs were differentially expressed only in aged wildtype mice, and 33 only in aged Tg4-42 mice. No significant alteration in the miRNome was detected in young mice, which indicates that the changes observed in aged mice are down-stream effects of Aβ-induced pathology in the Tg4-42 mouse model for AD. Targets of those miRs were predicted using miRWalk. For miRs that were differentially expressed only in the Tg4-42 model, 128 targets could be identified, whereas 18 genes were targeted by miRs only differentially expressed in wildtype mice and 85 genes were targeted by miRs differentially expressed in both mouse models. Genes targeted by differentially expressed miRs in the Tg4-42 model were enriched for negative regulation of long-term synaptic potentiation, learning or memory, regulation of trans-synaptic signaling and modulation of chemical synaptic transmission obtained. This untargeted miR sequencing approach supports previous reports on the Tg4-42 mice as a valuable model for AD. Furthermore, it revealed miRs involved in AD, which can serve as biomarkers or therapeutic targets.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

A combined OMICS screening approach of human plasma and serum was used to characterize protein and metabolome signatures displaying association to severity of Community-acquired pneumonia (CAP). 240 serum and BD P100 EDTA plasma samples from patients diagnosed with CAP, collected during the day of enrolment to the hospital, were analyzed by a metabolomic and proteomic approach, respectively. Disease severity of CAP patients was stratified using the Sequential Organ Failure Assessment (SOFA) score. Quantitative proteome and metabolome data, derived by LC-MS/MS, were associated to SOFA and specific parameters of SOFA using linear regression models adjusted for age, BMI, sex, smoking and technical variables. Both proteome and metabolome profiling revealed remarkable strong changes in plasma and serum composition in relation to severity of CAP. Proteins and metabolites displaying SOFA associated levels are involved in immune response, particularly in processes of lipid metabolism. Proteins, which show an association to SOFA score, are involved in acute phase response, coagulation, complement activation and inflammation. Many of these metabolites and proteins displayed not only associations to SOFA, but also to parameters of SOFA score, which likely reflect the strong influence of lung-, liver-, kidney- and heart-dysfunction on the metabolome and proteome patterns. SIGNIFICANCE: Community-acquired pneumonia is the most frequent infection disease with high morbidity and mortality. So far, only few studies focused on the identification of proteins or metabolites associated to severity of CAP, often based on smaller sample sets. A screening for new diagnostic markers requires extensive sample collections in combination with high quality clinical data. To characterize the proteomic and metabolomics pattern associated to severity of CAP we performed a combined metabolomics and proteomic approach of serum and plasma sample from a multi-center clinical study focused on patients with CAP, requiring hospitalization. The results of this association study of omics data to the SOFA score enable not only an interpretation of changes in molecular patterns with severity of CAP but also an assignment of altered molecules to dysfunctions of respiratory, renal, coagulation, cardiovascular systems as well as liver.

MicroRNAs (miRNAs) are small non-coding RNA molecules that bind to target sites in different gene regions and regulate post-transcriptional gene expression. Approximately 95% of human multi-exon genes can be spliced alternatively, which enables the production of functionally diverse transcripts and proteins from a single gene. Through alternative splicing, transcripts might lose the exon with the miRNA target site and become unresponsive to miRNA regulation. To check this hypothesis, we studied the role of miRNA target sites in both coding and non-coding regions using six cancer data sets from The Cancer Genome Atlas (TCGA) and Parkinson's disease data from PPMI. First, we predicted miRNA target sites on mRNAs from their sequence using TarPmiR. To check whether alternative splicing interferes with this regulation, we trained linear regression models to predict miRNA expression from transcript expression. Using nested models, we compared the predictive power of transcripts with miRNA target sites in the coding regions to that of transcripts without target sites. Models containing transcripts with target sites perform significantly better. We conclude that alternative splicing does interfere with miRNA regulation by skipping exons with miRNA target sites within the coding region.

C-reactive protein (CRP) is a sensitive biomarker of chronic low-grade inflammation and is associated with multiple complex diseases. The genetic determinants of chronic inflammation remain largely unknown, and the causal role of CRP in several clinical outcomes is debated. We performed two genome-wide association studies (GWASs), on HapMap and 1000 Genomes imputed data, of circulating amounts of CRP by using data from 88 studies comprising 204,402 European individuals. Additionally, we performed in silico functional analyses and Mendelian randomization analyses with several clinical outcomes. The GWAS meta-analyses of CRP revealed 58 distinct genetic loci (p < 5 × 10-8). After adjustment for body mass index in the regression analysis, the associations at all except three loci remained. The lead variants at the distinct loci explained up to 7.0% of the variance in circulating amounts of CRP. We identified 66 gene sets that were organized in two substantially correlated clusters, one mainly composed of immune pathways and the other characterized by metabolic pathways in the liver. Mendelian randomization analyses revealed a causal protective effect of CRP on schizophrenia and a risk-increasing effect on bipolar disorder. Our findings provide further insights into the biology of inflammation and could lead to interventions for treating inflammation and its clinical consequences.

Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10(-300), 2.1 × 10(-6), 2.5 × 10(-10) and 1.8 × 10(-10), respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months' less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

Determinations of thyrotropin (TSH) and free thyroxine (FT4) represent the gold standard in evaluation of thyroid function. To screen for novel peripheral biomarkers of thyroid function and to characterize FT4-associated physiological signatures in human plasma we used an untargeted OMICS approach in a thyrotoxicosis model.

Hypertension is the most important cause of death and disability in the elderly. In 9 out of 10 cases, the molecular cause, however, is unknown. One mechanistic hypothesis involves impaired endothelium-dependent vasodilation through reactive oxygen species (ROS) formation. Indeed, ROS forming NADPH oxidase (Nox) genes associate with hypertension, yet target validation has been negative. We re-investigate this association by molecular network analysis and identify NOX5, not present in rodents, as a sole neighbor to human vasodilatory endothelial nitric oxide (NO) signaling. In hypertensive patients, endothelial microparticles indeed contained higher levels of NOX5-but not NOX1, NOX2, or NOX4-with a bimodal distribution correlating with disease severity. Mechanistically, mice expressing human Nox5 in endothelial cells developed-upon aging-severe systolic hypertension and impaired endothelium-dependent vasodilation due to uncoupled NO synthase (NOS). We conclude that NOX5-induced uncoupling of endothelial NOS is a causal mechanism and theragnostic target of an age-related hypertension endotype. Nox5 knock-in (KI) mice represent the first mechanism-based animal model of hypertension.

Traditional drug discovery faces a severe efficacy crisis. Repurposing of registered drugs provides an alternative with lower costs and faster drug development timelines. However, the data necessary for the identification of disease modules, i.e. pathways and sub-networks describing the mechanisms of complex diseases which contain potential drug targets, are scattered across independent databases. Moreover, existing studies are limited to predictions for specific diseases or non-translational algorithmic approaches. There is an unmet need for adaptable tools allowing biomedical researchers to employ network-based drug repurposing approaches for their individual use cases. We close this gap with NeDRex, an integrative and interactive platform for network-based drug repurposing and disease module discovery. NeDRex integrates ten different data sources covering genes, drugs, drug targets, disease annotations, and their relationships. NeDRex allows for constructing heterogeneous biological networks, mining them for disease modules, prioritizing drugs targeting disease mechanisms, and statistical validation. We demonstrate the utility of NeDRex in five specific use-cases.

Meta-analysis has been established as an effective approach to combining summary statistics of several genome-wide association studies (GWAS). However, the accuracy of meta-analysis can be attenuated in the presence of cross-study heterogeneity. We present sPLINK, a hybrid federated and user-friendly tool, which performs privacy-aware GWAS on distributed datasets while preserving the accuracy of the results. sPLINK is robust against heterogeneous distributions of data across cohorts while meta-analysis considerably loses accuracy in such scenarios. sPLINK achieves practical runtime and acceptable network usage for chi-square and linear/logistic regression tests. sPLINK is available at https://exbio.wzw.tum.de/splink .

Helicobacter (H.) pylori is the most important cause for peptic ulcer disease and a risk factor for gastric carcinoma. How colonization with H. pylori affects the intestinal microbiota composition in humans is unknown. We investigated the association of H. pylori infection with intestinal microbiota composition in the population-based cohort Study-of-Health-in-Pomerania (SHIP)-TREND. Anti-H. pylori serology and H. pylori stool antigen tests were used to determine the H. pylori infection status. The fecal microbiota composition of 212 H. pylori positive subjects and 212 matched negative control individuals was assessed using 16S rRNA gene sequencing. H. pylori infection was found to be significantly associated with fecal microbiota alterations and a general increase in fecal microbial diversity. In infected individuals, the H. pylori stool antigen load determined a larger portion of the microbial variation than age or sex. The highest H. pylori stool antigen loads were associated with a putatively harmful microbiota composition. This study demonstrates profound alterations in human fecal microbiota of H. pylori infected individuals. While the increased microbiota diversity associated with H. pylori infection as well as changes in abundance of specific genera could be considered to be beneficial, others may be associated with adverse health effects, reflecting the complex relationship between H. pylori and its human host.

During disease progression or organism development, alternative splicing may lead to isoform switches that demonstrate similar temporal patterns and reflect the alternative splicing co-regulation of such genes. Tools for dynamic process analysis usually neglect alternative splicing.

Machine learning and artificial intelligence have shown promising results in many areas and are driven by the increasing amount of available data. However, these data are often distributed across different institutions and cannot be easily shared owing to strict privacy regulations. Federated learning (FL) allows the training of distributed machine learning models without sharing sensitive data. In addition, the implementation is time-consuming and requires advanced programming skills and complex technical infrastructures.

A long-term objective of network medicine is to replace our current, mainly phenotype-based disease definitions by subtypes of health conditions corresponding to distinct pathomechanisms. For this, molecular and health data are modeled as networks and are mined for pathomechanisms. However, many such studies rely on large-scale disease association data where diseases are annotated using the very phenotype-based disease definitions the network medicine field aims to overcome. This raises the question to which extent the biases mechanistically inadequate disease annotations introduce in disease association data distort the results of studies which use such data for pathomechanism mining. We address this question using global- and local-scale analyses of networks constructed from disease association data of various types. Our results indicate that large-scale disease association data should be used with care for pathomechanism mining and that analyses of such data should be accompanied by close-up analyses of molecular data for well-characterized patient cohorts.

Ionizing radiation is known to induce genomic lesions, such as DNA double strand breaks, whose repair can lead to mutations that can modulate cellular and organismal fate. Soon after radiation exposure, cells induce transcriptional changes and alterations of cell cycle programs to respond to the received DNA damage. Radiation-induced mutations occur through misrepair in a stochastic manner and increase the risk of developing cancers years after the incident, especially after high dose radiation exposures. Here, the authors analyzed the transcriptomic response of primary human gingival fibroblasts exposed to increasing doses of acute high dose-rate x rays. In the dataset obtained after 0.5 and 5 Gy x-ray exposures and two different repair intervals (0.5 h and 16 h), the authors discovered several radiation-induced fusion transcripts in conjunction with dose-dependent gene expression changes involving a total of 3,383 genes. Principal component analysis of repeated experiments revealed that the duration of the post-exposure repair intervals had a stronger impact than irradiation dose. Subsequent overrepresentation analyses showed a number of KEGG gene sets and WikiPathways, including pathways known to relate to radioresistance in fibroblasts (Wnt, integrin signaling). Moreover, a significant radiation-induced modulation of microRNA targets was detected. The data sets on IR-induced transcriptomic alterations in primary gingival fibroblasts will facilitate genomic comparisons in various genotoxic exposure scenarios.