Metabolome reflects the interplay of genome and exposome at molecular level and thus can provide deep insights into the pathogenesis of a complex disease like major depression. To identify metabolites associated with depression we performed a metabolome-wide association analysis in 13,596 participants from five European population-based cohorts characterized for depression, and circulating metabolites using ultra high-performance liquid chromatography/tandem accurate mass spectrometry (UHPLC/MS/MS) based Metabolon platform. We tested 806 metabolites covering a wide range of biochemical processes including those involved in lipid, amino-acid, energy, carbohydrate, xenobiotic and vitamin metabolism for their association with depression. In a conservative model adjusting for life style factors and cardiovascular and antidepressant medication use we identified 8 metabolites, including 6 novel, significantly associated with depression. In individuals with depression, increased levels of retinol (vitamin A), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1) (lecithin) and mannitol/sorbitol and lower levels of hippurate, 4-hydroxycoumarin, 2-aminooctanoate (alpha-aminocaprylic acid), 10-undecenoate (11:1n1) (undecylenic acid), 1-linoleoyl-GPA (18:2) (lysophosphatidic acid; LPA 18:2) are observed. These metabolites are either directly food derived or are products of host and gut microbial metabolism of food-derived products. Our Mendelian randomization analysis suggests that low hippurate levels may be in the causal pathway leading towards depression. Our findings highlight putative actionable targets for depression prevention that are easily modifiable through diet interventions.

Abstract In the era of personalized medicine with more and more patient-specific targeted therapies being used, we need reliable, dynamic, faster and sensitive biomarkers both to track the causes of disease and to develop and evolve therapies during the course of treatment. Metabolomics recently has shown substantial evidence to support its emerging role in disease diagnosis and prognosis. Aside from biomarkers and development of therapies, it is also an important goal to understand the involvement of mitochondrial DNA (mtDNA) in metabolic regulation, aging and disease development. Somatic mutations of the mitochondrial genome are also heavily implicated in age-related disease and aging. The general hypothesis is that an alteration in the concentration of metabolite profiles (possibly conveyed by lifestyle and environmental factors) influences the increase of mutation rate in the mtDNA and thereby contributes to a range of pathophysiological alterations observed in complex diseases. We performed an inverted mitochondrial genome-wide association analysis between mitochondrial nucleotide variants (mtSNVs) and concentration of metabolites. We used 151 metabolites and the whole sequenced mitochondrial genome from 2718 individuals to identify the genetic variants associated with metabolite profiles. Because of the high coverage, next-generation sequencing-based analysis of the mitochondrial genome allows for an accurate detection of mitochondrial heteroplasmy and for the identification of variants associated with the metabolome. The strongest association was found for mt715G > A located in the MT-12SrRNA with the metabolite ratio of C2/C10:1 (P-value = 6.82*10−09, β = 0.909). The second most significant mtSNV was found for mt3714A > G located in the MT-ND1 with the metabolite ratio of phosphatidylcholine (PC) ae C42:5/PC ae C44:5 (P-value = 1.02*10−08, β = 3.631). A large number of significant metabolite ratios were observed involving PC aa C36:6 and the variant mt10689G > A, located in the MT-ND4L gene. These results show an important interconnection between mitochondria and metabolite concentrations. Considering that some of the significant metabolites found in this study have been previously related to complex diseases, such as neurological disorders and metabolic conditions, these associations found here might play a crucial role for further investigations of such complex diseases. Understanding the mechanisms that control human health and disease, in particular, the role of genetic predispositions and their interaction with environmental factors is a prerequisite for the development of safe and efficient therapies for complex disorders.

Psychosis is a highly heritable and heterogeneous psychiatric condition. Its genetic architecture is thought to be the result of the joint effect of common and rare variants. Families with high prevalence are an interesting approach to shed light on the rare variant’s contribution without the need of collecting large cohorts. To unravel the genomic architecture of a family enriched for psychosis, with four affected individuals, we applied a system genomic approach based on karyotyping, genotyping by whole-exome sequencing to search for rare single nucleotide variants (SNVs) and SNP array to search for copy-number variants (CNVs). We identified a rare non-synonymous variant, g.39914279 C > G, in the MACF1 gene, segregating with psychosis. Rare variants in the MACF1 gene have been previously detected in SCZ patients. Besides, two rare CNVs, DUP3p26.3 and DUP16q23.3, were also identified in the family affecting relevant genes ( CNTN6 and CDH13 , respectively). We hypothesize that the co-segregation of these duplications with the rare variant g.39914279 C > G of MACF1 gene precipitated with schizophrenia and schizoaffective disorder.

Historical records document medieval immigration from North Africa to Iberia to create Islamic al-Andalus. Here, we present a low-coverage genome of an eleventh century CE man buried in an Islamic necropolis in Segorbe, near Valencia, Spain. Uniparental lineages indicate North African ancestry, but at the autosomal level he displays a mosaic of North African and European-like ancestries, distinct from any present-day population. Altogether, the genome-wide evidence, stable isotope results and the age of the burial indicate that his ancestry was ultimately a result of admixture between recently arrived Amazigh people (Berbers) and the population inhabiting the Peninsula prior to the Islamic conquest. We detect differences between our sample and a previously published group of contemporary individuals from Valencia, exemplifying how detailed, small-scale aDNA studies can illuminate fine-grained regional and temporal differences. His genome demonstrates how ancient DNA studies can capture portraits of past genetic variation that have been erased by later demographic shifts-in this case, most likely the seventeenth century CE expulsion of formerly Islamic communities as tolerance dissipated following the Reconquista by the Catholic kingdoms of the north.

Numerous predictive models for the risk of type 2 diabetes mellitus (T2DM) exist, but a minority of them has implemented nutrition data so far, even though the significant effect of nutrition on the pathogenesis, prevention and management of T2DM has been established. Thus, in the present study, we aimed to build a predictive model for the risk of T2DM that incorporates nutrition data and calculates its predictive performance. We analysed cross-sectional data from 1591 individuals from the population-based Cooperative Health Research in the Region of Augsburg (KORA) FF4 study (2013-14) and used a bootstrap enhanced elastic net penalised multivariate regression method in order to build our predictive model and select among 193 food intake variables. After selecting the significant predictor variables, we built a logistic regression model with these variables as predictors and T2DM status as the outcome. The values of area under the receiver operating characteristic (ROC) curve, sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV) and accuracy of our predictive model were calculated. Eleven out of the 193 food intake variables were selected for inclusion in our model, which yielded a value of area under the ROC curve of 0⋅79 and a maximum PPV, NPV and accuracy of 0⋅37, 0⋅98 and 0⋅91, respectively. The present results suggest that nutrition data should be implemented in predictive models to predict the risk of T2DM, since they improve their performance and they are easy to assess.

We investigated the genetic causes of major mental disorders (MMDs) including schizophrenia, bipolar disorder I, major depressive disorder and attention deficit hyperactive disorder, in a large family pedigree from Alpujarras, South of Spain, a region with high prevalence of psychotic disorders. We applied a systematic genomic approach based on karyotyping (n = 4), genotyping by genome-wide SNP array (n = 34) and whole-genome sequencing (n = 12). We performed genome-wide linkage analysis, family-based association analysis and polygenic risk score estimates. Significant linkage was obtained at chromosome 9 (9q33.1–33.2, LOD score = 4.11), a suggestive region that contains five candidate genes ASTN2 , BRINP1 , C5 , TLR4 and TRIM32, previously associated with MMDs. Comprehensive analysis associated the MMD phenotype with genes of the immune system with dual brain functions. Moreover, the psychotic phenotype was enriched for genes involved in synapsis. These results should be considered once studying the genetics of psychiatric disorders in other families, especially the ones from the same region, since founder effects may be related to the high prevalence.

Human remains from the Iron Age in Atlantic Scotland are rare, which makes the assemblage of an adult female and numerous foetal bones at High Pasture Cave, on the Isle of Skye, particularly noteworthy. Archaeological evidence suggests that the female had been deposited as an articulated skeleton when the cave entrance was blocked off, marking the end of use of the site. Particularly intriguing is the deposition of disarticulated remains from a foetus and perinate close to the adult female, which opens the possibility that the female might have been the mother of both of the infants. We used shotgun genome sequencing in order to analyse the mitochondrial genomes of all three individuals and investigate their maternal relationship, and we report here, for the first time, complete ancient mitogenomes from foetal-aged bone fragments. While we could not exclude the possibility that the female was the mother of, or maternally related to, the foetus, we could definitely say that she was not the mother of the perinate buried alongside her. This finding is contrary to the standard archaeological interpretation, that women in such burials most likely died in childbirth and were buried together with their foetuses.

Objectives: Our main goal is to identify mitochondrial susceptibility genes for human complex diseases. Methodology: Using different genotyping platforms and New Generation Sequencing we performed mtGWA analysis in the KORA population (3,080 individuals) with several phenotypes (BMI, cholesterol, post-traumatic stress disorder, thyroid diseases, depression, among others) as well as with life styles (smoking cigarettes, alcohol consumption, and physical activity). Results: we report several mitochondrial genetic variants associated with post-traumatic stress disorders and metabolic syndrome. BMI, HDL cholesterol and TG levels were also associated to mitochondrial variants, indicating that the presence of heteroplasmy in these variants may influence the balance of HDL cholesterol and TG levels. We also observed that no physical activity, alcohol abuse, and smoke cigarettes may increase the heteroplasmy levels of mitochondrial genetic variants. Conclusions: As fluctuations in mtDNA heteroplasmy may signify alterations in cellular activity, our findings highlight the important role of the mtDNA as a biomarker for aging-related diseases and metabolic syndromes.

Despite the fact that mitochondrial dysfunctions are increasingly recognized as key components in stress-related mental disorders, very little is known about the association between posttraumatic stress disorder (PTSD) and mitochondrial variants. To identify susceptibility mitochondrial genes for PTSD, we analyzed a total number of 978 mitochondrial single-nucleotide polymorphisms (mtSNPs) in a sample of 1238 individuals participating in the KORA (Cooperative Health Research in the Region of Augsburg) study. Participants were classified with ‘no PTSD’, ‘partial PTSD’ or ‘full PTSD’ by applying the Posttraumatic Diagnostic Scale and the Impact of Event Scale. To assess PTSD–mtSNP association while taking heteroplasmy into account, we used the raw signal intensity values measured on the microarray and applied linear regression. Significant associations were obtained between full versus no PTSD and two mtSNPs; mt8414C→T (β=−0.954±0.06, Padjusted=0.037) located in adenosine triphosphate (ATP) synthase subunit 8 (MT-ATP8) and mt12501G→A (β=−1.782±0.40, Padjusted=0.015) located in the NADH dehydrogenase subunits 5 (MT-ND5). Heteroplasmy for the two variants towards a larger number of the respective minor alleles increases the risk of having PTSD. NADH dehydrogenase and ATP synthase are both linked to the regulation of reactive oxygen species. Our results highlight the important role of the mitochondrial genome among the factors that contribute to the risk of PTSD. Mitochondrial genetic variants may be more important than has previously been assumed, leading to further insights regarding effects of existing medications, or even to the development of innovative treatments. As this is the first mitochondrial genome-wide association study for PTSDs, further analyses are needed to follow up on the present findings.

It has been suggested that mitochondrial dysfunction has an influence on lipid metabolism. The fact that mitochondrial defects can be accumulated over time as a normal part of aging may explain why cholesterol levels often are altered with age. To test the hypothesis whether mitochondrial variants are associated with lipid profile (total cholesterol, LDL, HDL, and triglycerides) we analyzed a total number of 978 mitochondrial single nucleotide polymorphisms (mtSNPs) in a sample of 2,815 individuals participating in the population-based KORA F4 study. To assess mtSNP association while taking the presence of heteroplasmy into account we used the raw signal intensity values measured on the microarray and applied linear regression. Ten mtSNPs (mt3285, mt3336, mt5285, mt6591, mt6671, mt9163, mt13855, mt13958, mt14000, and mt14580) were significantly associated with HDL cholesterol and one mtSNP (mt15074) with triglycerides levels. These results highlight the importance of the mitochondrial genome among the factors that contribute to the regulation of lipid levels. Focusing on mitochondrial variants may lead to further insights regarding the underlying physiological mechanisms, or even to the development of innovative treatments. Since this is the first mitochondrial genome-wide association analysis (mtGWAS) for lipid profile, further analyses are needed to follow up on the present findings.

It has been suggested that mitochondrial dysfunction plays a role in metabolic disorders including obesity, diabetes, and hypertension. The fact that mitochondrial defects can be accumulated over time as a normal part of aging may explain why some individuals can eat all sorts of foods and remain at normal weight while they are young. However, around the fourth decade of life there is a trend towards "middle-age spread'' with weight gain and the body's decreasing ability to metabolize calories efficiently. To test the hypothesis that mitochondrial variants are associated with BMI in adults, we analyzed a total number of 984 mitochondrial single nucleotide polymorphisms (mtSNPs) in a sample of 6,528 individuals participating in the KORA studies. To assess mtSNP association while taking heteroplasmy into account we used the raw signal intensity values measured on the microarray and applied linear regression. Significant results were obtained for 2 mtSNPs located in the Cytochrome c oxidase subunit genes (MT-CO1: P-adjusted = 0.0140 and MT-CO3: P-adjusted = 0.0286) and 3 mtSNPs located in the NADH dehydrogenase subunit genes (MT-ND1, MT-ND2 and MT-ND4L: P-adjusted = 0.0286). Polymorphisms located in the MT-CO3 and MT-ND4L genes have not been associated with BMI or related phenotypes in the past. Our results highlight the importance of the mitochondrial genome among the factors that contribute to the risk of high BMI. Focusing on mitochondrial variants may lead to further insights regarding effects of existing medications, or even to the development of innovative treatments.

Background: Congenital heart defects (CHD) is the most common cause of death from a congenital structure abnormality in newborns and is often associated with fetal loss. There are many types of CHD. Human genetic studies have identified genes that are responsible for the inheritance of a particular type of CHD and for some types of CHD previously thought to be sporadic. However, occasionally different members of the same family might have anatomically distinct defects - for instance, one member with atrial septal defect, one with tetralogy of Fallot, and one with ventricular septal defect. Our objective is to identify susceptibility loci for CHD in families affected by distinct defects. The occurrence of these apparently discordant clinical phenotypes within one family might hint at a genetic framework common to most types of CHD.Results: We performed a genome-wide linkage analysis using MOD score analysis in families with diverse CHD. Significant linkage was obtained in two regions, at chromosome 15 (15q26.3, P-empirical = 0.0004) and at chromosome 18 (18q21.2, P-empirical = 0.0005).Conclusions: In these two novel regions four candidate genes are located: SELS, SNRPA1, and PCSK6 on 15q26.3, and TCF4 on 18q21.2. The new loci reported here have not previously been described in connection with CHD. Although further studies in other cohorts are needed to confirm these findings, the results presented here together with recent insight into how the heart normally develops will improve the understanding of CHD.

It has been suggested that mitochondrial dysfunction plays a role in the pathogenesis of asthma. To test whether mitochondrial variants influence the risk of asthma, we analyzed 16,158 mtSNPs in a sample of 372 asthmatic children and 395 healthy children using the DNA pooling technique and genome wide association analysis. Stratified analysis by sex was performed to explain the differences observed between sexes in the etiology of asthma. Different variants were detected to be significant in the sample of girls and boys with the smallest adjusted p values being 1.4 × 10(-09) (mt5295) and 3.6 × 10(-12) (mt16158), respectively. Most of the significant locations found in boys are within the CYB gene and the non-coding region. For girls, most of the significant mtSNPs lie within NADH-dehydrogenase-subunits. The variants reported here have not previously been described in connection with asthma. Although further studies in other cohorts are needed to confirm these findings our study highlights the importance of the mitochondria among the factors that contribute to the risk of asthma.