Detailed stratified GWAS analysis for severe COVID-19 in four European populations.

Given the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), a deeper analysis of the host genetic contribution to severe COVID-19 is important to improve our understanding of underlying disease mechanisms. Here, we describe an extended genome-wide association meta-analysis of a well-characterized cohort of 3255 COVID-19 patients with respiratory failure and 12 488 population controls from Italy, Spain, Norway and Germany/Austria, including stratified analyses based on age, sex and disease severity, as well as targeted analyses of chromosome Y haplotypes, the human leukocyte antigen region and the SARS-CoV-2 peptidome. By inversion imputation, we traced a reported association at 17q21.31 to a ~0.9-Mb inversion polymorphism that creates two highly differentiated haplotypes and characterized the potential effects of the inversion in detail. Our data, together with the 5th release of summary statistics from the COVID-19 Host Genetics Initiative including non-Caucasian individuals, also identified a new locus at 19q13.33, including NAPSA, a gene which is expressed primarily in alveolar cells responsible for gas exchange in the lung.

Ancient Yersinia pestis genomes lack the virulence-associated YpfΦ prophage present in modern pandemic strains.

Yersinia pestis is the causative agent of at least three major plague pandemics (Justinianic, Medieval and Modern). Previous studies on ancient Y. pestis genomes revealed that several genomic alterations had occurred approximately 5000-3000 years ago and contributed to the remarkable virulence of this pathogen. How a subset of strains evolved to cause the Modern pandemic is less well-understood. Here, we examined the virulence-associated prophage (YpfΦ), which had been postulated to be exclusively present in the genomes of strains associated with the Modern pandemic. The analysis of two new Y. pestis genomes from medieval/early modern Denmark confirmed that the phage is absent from the genome of strains dating to this time period. An extended comparative genome analysis of over 300 strains spanning more than 5000 years showed that the prophage is found in the genomes of modern strains only and suggests an integration into the genome during recent Y. pestis evolution. The phage-encoded Zot protein showed structural homology to a virulence factor of Vibrio cholerae. Similar to modern Y. pestis, we observed phages with a common origin to YpfΦ in individual strains of other bacterial species. Our findings present an updated view on the prevalence of YpfΦ, which might contribute to our understanding of the host spectrum, geographical spread and virulence of Y. pestis responsible for the Modern pandemic.

DNA methylation QTL analysis identifies new regulators of human longevity.

Human longevity is a complex trait influenced by both genetic and environmental factors, whose interaction is mediated by epigenetic mechanisms like DNA methylation. Here, we generated genome-wide whole-blood methylome data from 267 individuals, of which 71 were long-lived (90-104 years), by applying reduced representation bisulfite sequencing. We followed a stringent two-stage analysis procedure using discovery and replication samples to detect differentially methylated sites (DMSs) between young and long-lived study participants. Additionally, we performed a DNA methylation quantitative trait loci analysis to identify DMSs that underlie the longevity phenotype. We combined the DMSs results with gene expression data as an indicator of functional relevance. This approach yielded 21 new candidate genes, the majority of which are involved in neurophysiological processes or cancer. Notably, two candidates (PVRL2, ERCC1) are located on chromosome 19q, in close proximity to the well-known longevity- and Alzheimer's disease-associated loci APOE and TOMM40. We propose this region as a longevity hub, operating on both a genetic (APOE, TOMM40) and an epigenetic (PVRL2, ERCC1) level. We hypothesize that the heritable methylation and associated gene expression changes reported here are overall advantageous for the LLI and may prevent/postpone age-related diseases and facilitate survival into very old age.

Quantitative comparison of within-sample heterogeneity scores for DNA methylation data.

DNA methylation is an epigenetic mark with important regulatory roles in cellular identity and can be quantified at base resolution using bisulfite sequencing. Most studies are limited to the average DNA methylation levels of individual CpGs and thus neglect heterogeneity within the profiled cell populations. To assess this within-sample heterogeneity (WSH) several window-based scores that quantify variability in DNA methylation in sequencing reads have been proposed. We performed the first systematic comparison of four published WSH scores based on simulated and publicly available datasets. Moreover, we propose two new scores and provide guidelines for selecting appropriate scores to address cell-type heterogeneity, cellular contamination and allele-specific methylation. Most of the measures were sensitive in detecting DNA methylation heterogeneity in these scenarios, while we detected differences in susceptibility to technical bias. Using recently published DNA methylation profiles of Ewing sarcoma samples, we show that DNA methylation heterogeneity provides information complementary to the DNA methylation level. WSH scores are powerful tools for estimating variance in DNA methylation patterns and have the potential for detecting novel disease-associated genomic loci not captured by established statistics. We provide an R-package implementing the WSH scores for integration into analysis workflows.

Long-Lived Individuals Show a Lower Burden of Variants Predisposing to Age-Related Diseases and a Higher Polygenic Longevity Score.

Longevity is a complex phenotype influenced by both environmental and genetic factors. The genetic contribution is estimated at about 25%. Despite extensive research efforts, only a few longevity genes have been validated across populations. Long-lived individuals (LLI) reach extreme ages with a relative low prevalence of chronic disability and major age-related diseases (ARDs). We tested whether the protection from ARDs in LLI can partly be attributed to genetic factors by calculating polygenic risk scores (PRSs) for seven common late-life diseases (Alzheimer's disease (AD), atrial fibrillation (AF), coronary artery disease (CAD), colorectal cancer (CRC), ischemic stroke (ISS), Parkinson's disease (PD) and type 2 diabetes (T2D)). The examined sample comprised 1351 German LLI (≥94 years, including 643 centenarians) and 4680 German younger controls. For all ARD-PRSs tested, the LLI had significantly lower scores than the younger control individuals (areas under the curve (AUCs): ISS = 0.59, p = 2.84 × 10-35; AD = 0.59, p = 3.16 × 10-25; AF = 0.57, p = 1.07 × 10-16; CAD = 0.56, p = 1.88 × 10-12; CRC = 0.52, p = 5.85 × 10-3; PD = 0.52, p = 1.91 × 10-3; T2D = 0.51, p = 2.61 × 10-3). We combined the individual ARD-PRSs into a meta-PRS (AUC = 0.64, p = 6.45 × 10-15). We also generated two genome-wide polygenic scores for longevity, one with and one without the TOMM40/APOE/APOC1 gene region (AUC (incl. TOMM40/APOE/APOC1) = 0.56, p = 1.45 × 10-5, seven variants; AUC (excl. TOMM40/APOE/APOC1) = 0.55, p = 9.85 × 10-3, 10,361 variants). Furthermore, the inclusion of nine markers from the excluded region (not in LD with each other) plus the APOE haplotype into the model raised the AUC from 0.55 to 0.61. Thus, our results highlight the importance of TOMM40/APOE/APOC1 as a longevity hub.

Ancient genomes reveal a high diversity of Mycobacterium leprae in medieval Europe.

Studying ancient DNA allows us to retrace the evolutionary history of human pathogens, such as Mycobacterium leprae, the main causative agent of leprosy. Leprosy is one of the oldest recorded and most stigmatizing diseases in human history. The disease was prevalent in Europe until the 16th century and is still endemic in many countries with over 200,000 new cases reported annually. Previous worldwide studies on modern and European medieval M. leprae genomes revealed that they cluster into several distinct branches of which two were present in medieval Northwestern Europe. In this study, we analyzed 10 new medieval M. leprae genomes including the so far oldest M. leprae genome from one of the earliest known cases of leprosy in the United Kingdom-a skeleton from the Great Chesterford cemetery with a calibrated age of 415-545 C.E. This dataset provides a genetic time transect of M. leprae diversity in Europe over the past 1500 years. We find M. leprae strains from four distinct branches to be present in the Early Medieval Period, and strains from three different branches were detected within a single cemetery from the High Medieval Period. Altogether these findings suggest a higher genetic diversity of M. leprae strains in medieval Europe at various time points than previously assumed. The resulting more complex picture of the past phylogeography of leprosy in Europe impacts current phylogeographical models of M. leprae dissemination. It suggests alternative models for the past spread of leprosy such as a wide spread prevalence of strains from different branches in Eurasia already in Antiquity or maybe even an origin in Western Eurasia. Furthermore, these results highlight how studying ancient M. leprae strains improves understanding the history of leprosy worldwide.

Early integration of pastoralism and millet cultivation in Bronze Age Eurasia.

Mobile pastoralists are thought to have facilitated the first trans-Eurasian dispersals of domesticated plants during the Early Bronze Age (ca 2500-2300 BC). Problematically, the earliest seeds of wheat, barley and millet in Inner Asia were recovered from human mortuary contexts and do not inform on local cultivation or subsistence use, while contemporaneous evidence for the use and management of domesticated livestock in the region remains ambiguous. We analysed mitochondrial DNA and multi-stable isotopic ratios (δ13C, δ15N and δ18O) of faunal remains from key pastoralist sites in the Dzhungar Mountains of southeastern Kazakhstan. At ca 2700 BC, Near Eastern domesticated sheep and goat were present at the settlement of Dali, which were also winter foddered with the region's earliest cultivated millet spreading from its centre of domestication in northern China. In the following centuries, millet cultivation and caprine management became increasingly intertwined at the nearby site of Begash. Cattle, on the other hand, received low levels of millet fodder at the sites for millennia. By primarily examining livestock dietary intake, this study reveals that the initial transmission of millet across the mountains of Inner Asia coincided with a substantial connection between pastoralism and plant cultivation, suggesting that pastoralist livestock herding was integral for the westward dispersal of millet from farming societies in China.

Atopobium and Fusobacterium as novel candidates for sarcoidosis-associated microbiota.

Sarcoidosis is a granulomatous disease that mainly affects the lung. A role of microbial factors in disease pathogenesis is assumed, but has not been investigated systematically in a large cohort.This cross-sectional study compared the lung microbiota of 71 patients with sarcoidosis, 15 patients with idiopathic pulmonary fibrosis (non-infectious controls) and 10 healthy controls (HCs). Next-generation sequencing of 16S DNA was used on bronchoalveolar lavage samples to characterise the microbial composition, which was analysed for diversity and indicator species. Host genotypes for 13 known sarcoidosis risk variants were determined and correlated with microbial parameters.The microbial composition differed significantly between sarcoidosis and HC samples (redundancy analysis ANOVA, p=0.025) and between radiographic Scadding types. Atopobium spp. was detected in 68% of sarcoidosis samples, but not in HC samples. Fusobacterium spp. was significantly more abundant in sarcoidosis samples compared with those from HCs. Mycobacteria were found in two of 71 sarcoidosis samples. Host-genotype analysis revealed an association of the rs2076530 (BTNL2) risk allele with a decrease in bacterial burden (p=0.002).Our results indicate Scadding type-dependent microbiota in sarcoidosis BAL samples. Atopobium spp. and Fusobacterium spp. were identified as sarcoidosis-associated bacteria, which may enable new insights into the pathogenesis and treatment of the disease.

Genetic Analysis Reveals a Longevity-Associated Protein Modulating Endothelial Function and Angiogenesis.

RATIONALE: Long living individuals show delay of aging, which is characterized by the progressive loss of cardiovascular homeostasis, along with reduced endothelial nitric oxide synthase activity, endothelial dysfunction, and impairment of tissue repair after ischemic injury. OBJECTIVE: Exploit genetic analysis of long living individuals to reveal master molecular regulators of physiological aging and new targets for treatment of cardiovascular disease. METHODS AND RESULTS: We show that the polymorphic variant rs2070325 (Ile229Val) in bactericidal/permeability-increasing fold-containing-family-B-member-4 (BPIFB4) associates with exceptional longevity, under a recessive genetic model, in 3 independent populations. Moreover, the expression of BPIFB4 is instrumental to maintenance of cellular and vascular homeostasis through regulation of protein synthesis. BPIFB4 phosphorylation/activation by protein-kinase-R-like endoplasmic reticulum kinase induces its complexing with 14-3-3 and heat shock protein 90, which is facilitated by the longevity-associated variant. In isolated vessels, BPIFB4 is upregulated by mechanical stress, and its knock-down inhibits endothelium-dependent vasorelaxation. In hypertensive rats and old mice, gene transfer of longevity-associated variant-BPIFB4 restores endothelial nitric oxide synthase signaling, rescues endothelial dysfunction, and reduces blood pressure levels. Furthermore, BPIFB4 is implicated in vascular repair. BPIFB4 is abundantly expressed in circulating CD34(+) cells of long living individuals, and its knock-down in endothelial progenitor cells precludes their capacity to migrate toward the chemoattractant SDF-1. In a murine model of peripheral ischemia, systemic gene therapy with longevity-associated variant-BPIFB4 promotes the recruitment of hematopoietic stem cells, reparative vascularization, and reperfusion of the ischemic muscle. CONCLUSIONS: Longevity-associated variant-BPIFB4 may represent a novel therapeutic tool to fight endothelial dysfunction and promote vascular reparative processes.

Coding and non-coding variants in the SHOX2 gene in patients with early-onset atrial fibrillation.

Atrial fibrillation (AF) is the most prevalent cardiac arrhythmia with a strong genetic component. Molecular pathways involving the homeodomain transcription factor Shox2 control the development and function of the cardiac conduction system in mouse and zebrafish. Here we report the analysis of human SHOX2 as a potential susceptibility gene for early-onset AF. To identify causal variants and define the underlying mechanisms, results from 378 patients with early-onset AF before the age of 60 years were analyzed and compared to 1870 controls or reference datasets. We identified two missense mutations (p.G81E, p.H283Q), that were predicted as damaging. Transactivation studies using SHOX2 targets and phenotypic rescue experiments in zebrafish demonstrated that the p.H283Q mutation severely affects SHOX2 pacemaker function. We also demonstrate an association between a 3'UTR variant c.*28T>C of SHOX2 and AF (p = 0.00515). Patients carrying this variant present significantly longer PR intervals. Mechanistically, this variant creates a functional binding site for hsa-miR-92b-5p. Circulating hsa-miR-92b-5p plasma levels were significantly altered in AF patients carrying the 3'UTR variant (p = 0.0095). Finally, we demonstrate significantly reduced SHOX2 expression levels in right atrial appendages of AF patients compared to patients with sinus rhythm. Together, these results suggest a genetic contribution of SHOX2 in early-onset AF.

APOE genotype and stress response - a mini review.

The APOE gene is one of currently only two genes that have consistently been associated with longevity. Apolipoprotein E (APOE) is a plasma protein which plays an important role in lipid and lipoprotein metabolism. In humans, there are three major APOE isoforms, designated APOE2, APOE3, and APOE4. Of these three isoforms, APOE3 is most common while APOE4 was shown to be associated with age-related diseases, including cardiovascular and Alzheimer's disease, and therefore an increased mortality risk with advanced age. Evidence accumulates, showing that oxidative stress and, correspondingly, mitochondrial function is affected in an APOE isoform-dependent manner. Accordingly, several stress response pathways implicated in the aging process, including the endoplasmic reticulum stress response and immune function, appear to be influenced by the APOE genotype. The investigation and development of treatment strategies targeting APOE4 have not resolved any therapeutic yet that could be entirely recommended. This mini-review provides an overview on the state of research concerning the impact of the APOE genotype on stress response-related processes, emphasizing the strong interconnection between mitochondrial function, endoplasmic reticulum stress and the immune response. Furthermore, this review addresses potential treatment strategies and associated pitfalls as well as lifestyle interventions that could benefit people with an at risk APOE4 genotype.

Identification of Immune-Relevant Factors Conferring Sarcoidosis Genetic Risk.

RATIONALE: Genetic variation plays a significant role in the etiology of sarcoidosis. However, only a small fraction of its heritability has been explained so far. OBJECTIVES: To define further genetic risk loci for sarcoidosis, we used the Immunochip for a candidate gene association study of immune-associated loci. METHODS: Altogether the study population comprised over 19,000 individuals. In a two-stage design, 1,726 German sarcoidosis cases and 5,482 control subjects were genotyped for 128,705 single-nucleotide polymorphisms using the Illumina Immunochip for the screening step. The remaining 3,955 cases, 7,514 control subjects, and 684 parents of affected offspring were used for validation and replication of 44 candidate and two established risk single-nucleotide polymorphisms. MEASUREMENTS AND MAIN RESULTS: Four novel susceptibility loci were identified with genome-wide significance in the European case-control populations, located on chromosomes 12q24.12 (rs653178; ATXN2/SH2B3), 5q33.3 (rs4921492; IL12B), 4q24 (rs223498; MANBA/NFKB1), and 2q33.2 (rs6748088; FAM117B). We further defined three independent association signals in the HLA region with genome-wide significance, peaking in the BTNL2 promoter region (rs5007259), at HLA-B (rs4143332/HLA-B*0801) and at HLA-DPB1 (rs9277542), and found another novel independent signal near IL23R (rs12069782) on chromosome 1p31.3. CONCLUSIONS: Functional predictions and protein network analyses suggest a prominent role of the drug-targetable IL23/Th17 signaling pathway in the genetic etiology of sarcoidosis. Our findings reveal a substantial genetic overlap of sarcoidosis with diverse immune-mediated inflammatory disorders, which could be of relevance for the clinical application of modern therapeutics.

Allelic and copy-number variations of FcγRs affect granulocyte function and susceptibility for autoimmune blistering diseases.

Low-affinity Fcγ receptors (FcγR) bridge innate and adaptive immune responses. In many autoimmune diseases, these receptors act as key mediators of the pathogenic effects of autoantibodies. Genes encoding FcγR exhibit frequent variations in sequence and gene copy number that influence their functional properties. FcγR variations also affect the susceptibility to systemic autoimmunity, e.g. systemic lupus erythematosus and rheumatoid arthritis. This raises the question whether FcγR variations are also associated with organ-specific autoimmunity, particularly autoantibody-mediated diseases, such as subepidermal autoimmune blistering diseases (AIBD). A multitude of evidence suggests a pathogenic role of neutrophil granulocyte interaction with autoantibodies via FcγR. In a two-stage study, we analyzed whether the FcγR genotype affects neutrophil function and mRNA expression, and consequently, bullous pemphigoid (BP) disease risk. We compared this to findings in pemphigus vulgaris/foliaceus (PV/PF), two Fc-independent AIBDs. Our results indicate that both allele and copy number variation of FcγR genes affect FcγR mRNA expression and reactive oxygen species (ROS) release by granulocytes. Susceptibility of BP was associated with FcγR genotypes that led to a decreased ROS release by neutrophils, indicating an unexpected protective role for these cells. BP and PV/PF differed substantially regarding the FcγR genotype association patterns, pointing towards different disease etiologies.

Toward the blood-borne miRNome of human diseases.

In a multicenter study, we determined the expression profiles of 863 microRNAs by array analysis of 454 blood samples from human individuals with different cancers or noncancer diseases, and validated this 'miRNome' by quantitative real-time PCR. We detected consistently deregulated profiles for all tested diseases; pathway analysis confirmed disease association of the respective microRNAs. We observed significant correlations (P = 0.004) between the genomic location of disease-associated genetic variants and deregulated microRNAs.

Genetics of sarcoidosis.

Sarcoidosis is a multifactorial and polygenic disorder. Recently, several novel predisposing genes have been identified by genome-wide association studies, and fast progress in molecular technologies such as systematic and large-scale resequencing will aid the discovery of further risk loci and variants. In this article, the current knowledge of its genetics will be presented, including known and candidate risk variants and loci, with a focus on loci in the human leukocyte antigen region. Some of these factors are shared with other, clinically distinct diseases. This may lead to the development of new hypotheses on pathomechanisms, which associate sarcoidosis with other granulomatous disorders but also with diseases with significantly different phenotypes. In the near future system, biology approaches will help unravel the differing and common features of these disorders and allow the development of new therapeutic strategies and tools to predict the course and response to treatment of individual patients.

SLC23A1 polymorphism rs6596473 in the vitamin C transporter SVCT1 is associated with aggressive periodontitis.

AIM: Identification of variants within genes SLC23A1 and SLC23A2 coding for vitamin C transporter proteins associated with aggressive (AgP) and chronic periodontitis (CP). MATERIAL AND METHODS: Employment of three independent case-control samples of AgP (I. 283 cases, 979 controls; II. 417 cases, 1912 controls; III. 164 cases, 357 controls) and one sample of CP (1359 cases, 1296 controls). RESULTS: Stage 1: Among the tested single-nucleotide polymorphisms (SNPs), the rare allele (RA) of rs6596473 in SLC23A1 showed nominal significant association with AgP (p = 0.026, odds ratio [OR] 1.26, and a highly similar minor allele frequency between different control panels. Stage 2: rs6596473 showed no significant association with AgP in the replication with the German and Dutch case-control samples. After pooling the German AgP populations (674 cases, 2891 controls) to significantly increase the statistical power (SP = 0.81), rs6596473 RA showed significant association with AgP prior to and upon adjustment with the covariates smoking and gender with padj  = 0.005, OR = 1.35. Stage 3: RA of rs6596473 showed no significant association with severe CP. CONCLUSION: SNP rs6596473 of SLC23A1 is suggested to be associated with AgP. These results add to previous reports that vitamin C plays a role in the pathogenesis of periodontitis.

Neolithic Yersinia pestis infections in humans and a dog.

Yersinia pestis has been infecting humans since the Late Neolithic (LN). Whether those early infections were isolated zoonoses or initiators of a pandemic remains unclear. We report Y. pestis infections in two individuals (of 133) from the LN necropolis at Warburg (Germany, 5300-4900 cal BP). Our analyses show that the two genomes belong to distinct strains and reflect independent infection events. All LN genomes known today (n = 4) are basal in the phylogeny and represent separate lineages that probably originated in different animal hosts. In the LN, an opening of the landscape resulted in the introduction of new rodent species, which may have acted as Y. pestis reservoirs. Coincidentally, the number of dogs increased, possibly leading to Y. pestis infections in canines. Indeed, we detect Y. pestis in an LN dog. Collectively, our data suggest that Y. pestis frequently entered human settlements at the time without causing significant outbreaks.

Polymorphisms in the superoxidase dismutase genes reveal no association with human longevity in Germans: a case-control association study.

The role of superoxide dismutases (SODs) in aging and oxidative stress regulation has been widely studied and there is growing evidence that imbalances in these processes influence lifespan in several species. In humans, genetic polymorphisms in SOD genes may play an important role in the development of age-related diseases and genetic variation in SOD2 is thought to be associated with longevity. These observations prompted us to perform a case-control association study using a comprehensive haplotype tagging approach for the three SOD genes (SOD1, SOD2, SOD3) by testing a total of 19 SNPs in our extensive collection of 1,612 long-lived individuals (centenarians and nonagenarians) and 1,104 younger controls. Furthermore, we intended to replicate the previous association of the SOD2 SNP rs4880 with longevity observed in a Danish cohort. In our study, no association was detected between the tested SNPs and the longevity phenotype, neither in the entire long-lived sample set nor in the centenarian subgroup analysis. Our results suggest that there is no considerable influence of sequence variation in the SOD genes on human longevity in Germans.

A novel sarcoidosis risk locus for Europeans on chromosome 11q13.1.

RATIONALE: Sarcoidosis is a complex inflammatory disease with a heterogeneous clinical picture. Among others, an acute and chronic clinical course can be distinguished, for which specific genetic risk factors are known. OBJECTIVES: To identify additional risk loci for sarcoidosis and its acute and chronic subforms, we analyzed imputed data from a genome-wide association scan for these phenotypes. METHODS: After quality control, the genome-wide association scan comprised nearly 1.3 million imputed single-nucleotide polymorphisms based on an Affymetrix 6.0 Gene Chip dataset of 564 German sarcoidosis cases, including 176 acute and 354 chronic cases and 1,575 control subjects. MEASUREMENTS AND MAIN RESULTS: We identified chromosome 11q13.1 (rs479777) as a novel locus influencing susceptibility to sarcoidosis with genome-wide significance. The marker was significantly associated in three distinct German case-control populations and in an additional German family sample with odds ratios ranging from 0.67 to 0.77. This finding was further replicated in two independent European case-control populations from the Czech Republic (odds ratio, 0.75) and from Sweden (odds ratio, 0.79). In a meta-analysis of the included European case-control samples the marker yielded a P value of 2.68 × 10(-18). The locus was previously reported to be associated with Crohn disease, psoriasis, alopecia areata, and leprosy. For sarcoidosis, fine-mapping and expression analysis suggest KCNK4, PRDX5, PCLB3, and most promising CCDC88B as candidates for the underlying risk gene in the associated region. CONCLUSIONS: This study provides striking evidence for association of chromosome 11q13.1 with sarcoidosis in Europeans, and thus identified a further genetic risk locus shared by sarcoidosis, Crohn disease and psoriasis.