Genetic variation in the Sorbs of eastern Germany in the context of broader European genetic diversity

Association of Three Genetic Loci with Molar Incisor Hypomineralization in Polish Children

BACKGROUND

Molar incisor hypomineralization (MIH) is a qualitative, demarcated enamel defect of hypomineralization affecting one to four first permanent molars, often with incisor involvement. Its etiology is complex. However, evidence suggests the influence of genetic factors, potentially including the single nucleotide polymorphisms (SNPs) rs2889956, rs4811117 and rs13058467, which were previously linked to MIH in a genome-wide association study of German children. The aim was to replicate analyses of possible associations between the SNPs and molar incisor hypomineralization in Polish children.

METHODS

The final study group consisted of 778 children aged 126-168 months old. Saliva samples were taken, and genomic DNA was extracted and genotyped using beadchip microarrays.

RESULTS

Among the 778 subjects, there were 68 (8.7%) subjects with MIH and 710 (91.3%) subjects without MIH. There were no significant differences in distributions in age, sex, or the frequency of caries in permanent dentition between the MIH and non-MIH groups. The rs2889956, rs4811117, and rs13058467 genotype distributions in the studied group conformed to the expected Hardy-Weinberg equilibria, and there were no significant differences in the distributions of their alleles or genotypes between the MIH and non-MIH groups.

CONCLUSION

Our replication study did not confirm highly significant associations between the single nucleotide polymorphisms rs2889956, rs4811117, and rs13058467 with molar incisor hypomineralization in Polish children.

An atlas of genome-wide gene expression and metabolite associations and possible mediation effects towards body mass index

Investigating the cross talk of different omics layers is crucial to understand molecular pathomechanisms of metabolic diseases like obesity. Here, we present a large-scale association meta-analysis of genome-wide whole blood and peripheral blood mononuclear cell (PBMC) gene expressions profiled with Illumina HT12v4 microarrays and metabolite measurements from dried blood spots (DBS) characterized by targeted liquid chromatography tandem mass spectrometry (LC-MS/MS) in three large German cohort studies with up to 7706 samples. We found 37,295 associations comprising 72 amino acids (AA) and acylcarnitine (AC) metabolites (including ratios) and 8579 transcripts. We applied this catalogue of associations to investigate the impact of associating transcript-metabolite pairs on body mass index (BMI) as an example metabolic trait. This is achieved by conducting a comprehensive mediation analysis considering metabolites as mediators of gene expression effects and vice versa. We discovered large mediation networks comprising 27,023 potential mediation effects within 20,507 transcript-metabolite pairs. Resulting networks of highly connected (hub) transcripts and metabolites were leveraged to gain mechanistic insights into metabolic signaling pathways. In conclusion, here, we present the largest available multi-omics integration of genome-wide transcriptome data and metabolite data of amino acid and fatty acid metabolism and further leverage these findings to characterize potential mediation effects towards BMI proposing candidate mechanisms of obesity and related metabolic diseases. KEY MESSAGES: Thousands of associations of 72 amino acid and acylcarnitine metabolites and 8579 genes expand the knowledge of metabolome-transcriptome associations. A mediation analysis of effects on body mass index revealed large mediation networks of thousands of obesity-related gene-metabolite pairs. Highly connected, potentially mediating hub genes and metabolites enabled insight into obesity and related metabolic disease pathomechanisms.

Genome-wide meta-analysis of phytosterols reveals five novel loci and a detrimental effect on coronary atherosclerosis

Phytosterol serum concentrations are under tight genetic control. The relationship between phytosterols and coronary artery disease (CAD) is controversially discussed. We perform a genome-wide meta-analysis of 32 phytosterol traits reflecting resorption, cholesterol synthesis and esterification in six studies with up to 9758 subjects and detect ten independent genome-wide significant SNPs at seven genomic loci. We confirm previously established associations at ABCG5/8 and ABO and demonstrate an extended locus heterogeneity at ABCG5/8 with different functional mechanisms. New loci comprise HMGCR, NPC1L1, PNLIPRP2, SCARB1 and APOE. Based on these results, we perform Mendelian Randomization analyses (MR) revealing a risk-increasing causal relationship of sitosterol serum concentrations and CAD, which is partly mediated by cholesterol. Here we report that phytosterols are polygenic traits. MR add evidence of both, direct and indirect causal effects of sitosterol on CAD.

m6A Regulators in Human Adipose Tissue - Depot-Specificity and Correlation With Obesity

BACKGROUND

N⁶-methyladenosine (m6A) is one of the most abundant post-transcriptional modifications on mRNA influencing mRNA metabolism. There is emerging evidence for its implication in metabolic disease. No comprehensive analyses on gene expression of m6A regulators in human adipose tissue, especially in paired adipose tissue depots, and its correlation with clinical variables were reported so far. We hypothesized that inter-depot specific gene expression of m6A regulators may differentially correlate with clinical variables related to obesity and fat distribution.

METHODS

We extracted intra-individually paired gene expression data (omental visceral adipose tissue (OVAT) N=48; subcutaneous adipose tissue (SAT) N=56) of m6A regulators from an existing microarray dataset. We also measured gene expression in another sample set of paired OVAT and SAT (N=46) using RT-qPCR. Finally, we extracted existing gene expression data from peripheral mononuclear blood cells (PBMCs) and single nucleotide polymorphisms (SNPs) in METTL3 and YTHDF3 from genome wide data from the Sorbs population (N=1049). The data were analysed for differential gene expression between OVAT and SAT; and for association with obesity and clinical variables. We further tested for association of SNP markers with gene expression and clinical traits.

RESULTS

In adipose tissue we observed that several m6A regulators (WTAP, VIRMA, YTHDC1 and ALKBH5) correlate with obesity and clinical variables. Moreover, we found adipose tissue depot specific gene expression for METTL3, WTAP, VIRMA, FTO and YTHDC1. In PBMCs, we identified ALKBH5 and YTHDF3 correlated with obesity. Genetic markers in METTL3 associate with BMI whilst SNPs in YTHDF3 are associated with its gene expression.

CONCLUSIONS

Our data show that expression of m6A regulators correlates with obesity, is adipose tissue depot-specific and related to clinical traits. Genetic variation in m6A regulators adds an additional layer of variability to the functional consequences.

Genetically programmed changes in transcription of the novel progranulin regulator

Abstract

Progranulin is a glycoprotein marking chronic inflammation in obesity and type 2 diabetes. Previous studies suggested PSRC1 (proline and serine rich coiled-coil 1) to be a target of genetic variants associated with serum progranulin levels. We aimed to identify potentially functional variants and characterize their role in regulation of PSRC1. Phylogenetic module complexity analysis (PMCA) prioritized four polymorphisms (rs12740374, rs629301, rs660240, rs7528419) altering transcription factor binding sites with an overall score for potential regulatory function of S_all > 7.0. The effects of these variants on transcriptional activity and binding of transcription factors were tested by luciferase reporter and electrophoretic mobility shift assays (EMSA). In parallel, blood DNA promoter methylation of two regions was tested in subjects with a very high (N = 100) or a very low (N = 100) serum progranulin. Luciferase assays revealed lower activities in vectors carrying the rs629301-A compared with the C allele. Moreover, EMSA indicated a different binding pattern for the two rs629301 alleles, with an additional prominent band for the A allele, which was finally confirmed with the supershift for the Yin Yang 1 transcription factor (YY1). Subjects with high progranulin levels manifested a significantly higher mean DNA methylation (P < 1 × 10⁻⁷) in one promoter region, which was in line with a significantly lower PSRC1 mRNA expression levels in blood (P = 1 × 10⁻³). Consistently, rs629301-A allele was associated with lower PSRC1 mRNA expression (P < 1 × 10⁻⁷). Our data suggest that the progranulin-associated variant rs629301 modifies the transcription of PSRC1 through alteration of YY1 binding capacity. DNA methylation studies further support the role of PSRC1 in regulation of progranulin serum levels.

Key messages

PSRC1 (proline and serine rich coiled-coil 1) SNPs are associated with serum progranulin levels.

rs629301 regulates PSRC1 expression by affecting Yin Yang 1 transcription factor (YY1) binding.

PSRC1 is also epigenetically regulated in subjects with high progranulin levels.

Electronic supplementary material

The online version of this article (10.1007/s00109-020-01942-7) contains supplementary material, which is available to authorized users.

Adipocytokines are not associated with gestational diabetes mellitus but with pregnancy status

AIMS

Adipose tissue-secreted proteins, i.e. adipocytokines, have been identified as potential mediators linking fat mass and adipose tissue dysfunction with impaired glucose homeostasis, alterations in the inflammatory status, and risk of diabetes. The aim of this study was to determine whether seven circulating adipocytokines are associated with gestational diabetes mellitus (GDM) or are altered by metabolic and weight changes during pregnancy itself.

METHODS

A panel of seven adipocytokines (i.e. adiponectin, adipocyte fatty acid-binding protein, chemerin, leptin, Pro-Enkephalin, progranulin, and Pro-Neurotensin) was quantified in serum in a cross-sectional cohort of 222 women with the following three groups matched for age and body mass index: (i) 74 pregnant women with GDM; (ii) 74 pregnant women without GDM; and (iii) 74 non-pregnant and healthy women. A stepwise statistical approach was used by performing pairwise comparisons, principal component analysis (PCA), and partial least square discriminant analysis (PLS-DA).

RESULTS

Five out of seven adipocytokines were dysregulated between pregnant and non-pregnant women, i.e. adiponectin, chemerin, leptin, Pro-Enkephalin, and progranulin. None of the adipocytokines significantly differed between GDM and non-GDM status during pregnancy. The same five adipocytokines clustered in a principal component representing pregnancy-induced effects. Fasting insulin was the most relevant parameter in the discrimination of GDM as compared to pregnant women without GDM, whereas chemerin and adiponectin were most relevant factors to discriminate pregnancy status.

CONCLUSIONS

Pregnancy status but not presence of GDM can be distinguished by the seven investigated adipocytokines in discrimination analyses.

Clinical and lifestyle related factors influencing whole blood metabolite levels - A comparative analysis of three large cohorts

Objective

Human blood metabolites are influenced by a number of lifestyle and environmental factors. Identification of these factors and the proper quantification of their relevance provides insights into human biological and metabolic disease processes, is key for standardized translation of metabolite biomarkers into clinical applications, and is a prerequisite for comparability of data between studies. However, so far only limited data exist from large and well-phenotyped human cohorts and current methods for analysis do not fully account for the characteristics of these data. The primary aim of this study was to identify, quantify and compare the impact of a comprehensive set of clinical and lifestyle related factors on metabolite levels in three large human cohorts. To achieve this goal, we improve current methodology by developing a principled analysis approach, which could be translated to other cohorts and metabolite panels.

Methods

63 Metabolites (amino acids, acylcarnitines) were quantified by liquid chromatography tandem mass spectrometry in three cohorts (total N = 16,222). Supported by a simulation study evaluating various analytical approaches, we developed an analysis pipeline including preprocessing, identification, and quantification of factors affecting metabolite levels. We comprehensively identified uni- and multivariable metabolite associations considering 29 environmental and clinical factors and performed metabolic pathway enrichment and network analyses.

Results

Inverse normal transformation of batch corrected and outlier removed metabolite levels accompanied by linear regression analysis proved to be the best suited method to deal with the metabolite data. Association analyses revealed numerous uni- and multivariable significant associations. 15 of the analyzed 29 factors explained >1% of variance for at least one of the metabolites. Strongest factors are application of steroid hormones, reticulocytes, waist-to-hip ratio, sex, haematocrit, and age. Effect sizes of factors are comparable across studies.

Conclusions

We introduced a principled approach for the analysis of MS data allowing identification, and quantification of effects of clinical and lifestyle factors with metabolite levels. We detected a number of known and novel associations broadening our understanding of the regulation of the human metabolome. The large heterogeneity observed between cohorts could almost completely be explained by differences in the distribution of influencing factors emphasizing the necessity of a proper confounder analysis when interpreting metabolite associations.

•

Amino-acids and acylcarnitines analyzed in three studies with >16,000 individuals.

•

Develop a generic and adaptable bioinformatics workflow.

•

Analysis of the impact of 29 clinical and life-style factors on blood metabolites.

•

Analysis of network between factors and metabolites.

•

Comparison of results between studies.

(Epi)genetic regulation of CRTC1 in human eating behaviour and fat distribution

BACKGROUND

In brain, CREB-regulated transcription co-activator 1 (CRTC1) is involved in metabolic dysregulation. In humans a SNP in CRTC1 was associated to body fat percentage and two SNPs affected RNA Pol II binding and chromatin structure, implying epigenetic regulation of CRTC1. We sought to understand the relevance of CRTC1 SNPs, DNA methylation and expression in human eating behaviour and its relationship to clinical variables of obesity in blood and adipose tissue.

METHODS

13 CRTC1 SNPs were included to analyze eating behaviour. For rs7256986, follow up association analyses were applied on DNA methylation, CRTC1 expression and clinical parameters. Linear regression was used throughout the study adjusted for age, sex and BMI. Besides data extraction from previous work, rs7256986 was de-novo genotyped and DNA methylation was evaluated by using pyrosequencing.

FINDINGS

We found several SNPs in the CRTC1 locus nominally associated with human eating behaviour or 2hr postprandial insulin levels and observed a correlation with alcohol and coffee intake (all P < 0.05). G-allele carriers of rs7256986 showed slightly increased hip circumference. We showed that rs7256986 represents a methylation quantitative trait locus (meQTL) in whole blood and adipose tissue. The presence of the SNP and/or DNA methylation correlated with CRTC1 gene expression which in turn, related to BMI and fat distribution.

INTERPRETATION

Our data support the known role of CRCT1 regulating energy metabolism in brain. Here, we highlight relevance of CRTC1 regulation in blood and adipose tissue. FUND: IFB AdiposityDiseases (BMBF); n609020-Scientia Fellows; Helse-SørØst; DFG: CRC 1052/1 and/2; Kompetenznetz Adipositas, German Diabetes Association.

Genome-wide meta-analysis identifies novel determinants of circulating serum progranulin

Progranulin is a secreted protein with important functions in processes including immune and inflammatory response, metabolism and embryonic development. The present study aimed at identification of genetic factors determining progranulin concentrations. We conducted a genome-wide association meta-analysis for serum progranulin in three independent cohorts from Europe: Sorbs (N = 848) and KORA (N = 1628) from Germany and PPP-Botnia (N = 335) from Finland (total N = 2811). Single nucleotide polymorphisms (SNPs) associated with progranulin levels were replicated in two additional German cohorts: LIFE-Heart Study (Leipzig; N = 967) and Metabolic Syndrome Berlin Potsdam (Berlin cohort; N = 833). We measured mRNA expression of genes in peripheral blood mononuclear cells (PBMC) by micro-arrays and performed mRNA expression quantitative trait and expression-progranulin association studies to functionally substantiate identified loci. Finally, we conducted siRNA silencing experiments in vitro to validate potential candidate genes within the associated loci. Heritability of circulating progranulin levels was estimated at 31.8% and 26.1% in the Sorbs and LIFE-Heart cohort, respectively. SNPs at three loci reached study-wide significance (rs660240 in CELSR2-PSRC1-MYBPHL-SORT1, rs4747197 in CDH23-PSAP and rs5848 in GRN) explaining 19.4%/15.0% of the variance and 61%/57% of total heritability in the Sorbs/LIFE-Heart Study. The strongest evidence for association was at rs660240 (P = 5.75 × 10-50), which was also associated with mRNA expression of PSRC1 in PBMC (P = 1.51 × 10-21). Psrc1 knockdown in murine preadipocytes led to a consecutive 30% reduction in progranulin secretion. In conclusion, the present meta-GWAS combined with mRNA expression identified three loci associated with progranulin and supports the role of PSRC1 in the regulation of progranulin secretion.

Circulating Pro-Neurotensin in gestational diabetes mellitus

BACKGROUND AND AIMS

Pro-Neurotensin (NT), a stable surrogate parameter of NT, has recently been introduced as a peptide predicting the development of obesity, diabetes mellitus, cardiovascular diseases, and cardiovascular mortality. However, regulation of Pro-NT in gestational diabetes mellitus (GDM) remains uninvestigated.

METHODS AND RESULTS

Pro-NT was quantified in 74 women with GDM, 74 healthy, gestational age-matched, pregnant controls, as well as in a second cohort comprising of 74 healthy, non-pregnant control women, using a chemiluminometric sandwich immunoassay. Pro-NT was correlated to measures of obesity, hypertension, glucose and lipid metabolism, renal function, and inflammation. Mean ± standard deviation of circulating Pro-NT levels were not significantly different in women with GDM (100.2 ± 75.7 pmol/l) as compared to healthy, pregnant controls (103.2 ± 37.4 pmol/l) and healthy, non-pregnant female controls (105.9 ± 38.9 pmol/l) (p = 0.661). Homeostasis model assessment of insulin resistance (HOMA-IR) and creatinine positively correlated with serum Pro-NT in multivariate regression analysis. In contrast, free fatty acids (FFA) were inversely correlated with circulating Pro-NT. Results sustained adjustment for pregnancy status.

CONCLUSIONS

Circulating Pro-NT is not independently associated with GDM, but is with HOMA-IR, creatinine, and FFA even after adjustment for pregnancy status.

Siberian genetic diversity reveals complex origins of the Samoyedic-speaking populations

OBJECTIVES

We examined autosomal genome-wide SNPs and Y-chromosome data from 15 Siberian and 12 reference populations to study the affinities of Siberian populations, and to address hypotheses about the origin of the Samoyed peoples.

METHODS

Samples were genotyped for 567 096 autosomal SNPs and 147 Y-chromosome polymorphic sites. For several analyses, we used 281 093 SNPs from the intersection of our data with publicly available ancient Siberian samples. To examine genetic relatedness among populations, we applied PCA, F_ST , TreeMix, and ADMIXTURE analyses. To explore the potential effect of demography and evolutionary processes, the distribution of ROH and IBD sharing within population were studied.

RESULTS

Analyses of autosomal and Y-chromosome data reveal high differentiation of the Siberian groups. The Siberian populations have a large proportion of their genome in ROH and IBD segments. Several populations (ie, Nganasans, Evenks, Yukagirs, and Koryaks) do not appear to have experienced admixture with other Siberian populations (ie, producing only positive f3), while for the other tested populations the composition of mixing sources always included Nganasans or Evenks. The Nganasans from the Taymyr Peninsula demonstrate the greatest level of shared shorter ROH and IBD with nearly all other Siberian populations.

CONCLUSIONS

Autosomal SNP and Y-chromosome data demonstrate that Samoyedic populations differ significantly in their genetic composition. Genetic relationship is observed only between Forest and Tundra Nentsi. Selkups are affiliated with the Kets from the Yenisey River, while the Nganasans are separated from their linguistic neighbors, showing closer affinities with the Evenks and Yukagirs.

Metabolic effects of genetic variation in the human REPIN1 gene

BACKGROUND

Replication initiator 1 (Repin1) is a zinc finger protein highly expressed in liver and adipose tissue. The Repin1 resides within a quantitative trait locus (QTL) for body weight and triglyceride levels in the rat, and its hepatic deletion in mice results in improved insulin sensitivity and lower body weight. Here, we analyzed whether genetic variation within the Repin1 affects parameters of glucose and lipid metabolism.

METHODS

We sequenced REPIN1 in 48 non-related Caucasian subjects. We discovered a 12 base pair deletion (12 bp del; rs3832490), which was subsequently genotyped in two well-characterized cohorts (N = 3013) to test for associations with metabolic traits. Functional consequences of the variant were investigated in HepG2 cells in vitro.

RESULTS

In human cohorts, we show that the 12 bp del associates with improved glucose metabolism (lower fasting plasma glucose, fasting plasma insulin, and HOMA IR). Cells transfected with the plasmid carrying the 12 bp del variant are characterized by increased GLUT2 and fatty acid translocase CD36 expression and more lipid droplets.

CONCLUSION

Our data suggest that genetic variation in human REPIN1 plays a role in glucose and lipid metabolism by differentially affecting the expression of REPIN1 target genes including glucose and fatty acid transporters.

Widely Used Commercial ELISA Does Not Detect Precursor of Haptoglobin2, but Recognizes Properdin as a Potential Second Member of the Zonulin Family

BACKGROUND

There is increasing evidence for the role of impaired intestinal permeability in obesity and associated metabolic diseases. Zonulin is an established serum marker for intestinal permeability and identical to pre-haptoglobin2. Here, we aimed to investigate the relationship between circulating zonulin and metabolic traits related to obesity.

METHODS

Serum zonulin was measured by using a widely used commercial ELISA kit in 376 subjects from the metabolically well-characterized cohort of Sorbs from Germany. In addition, haptoglobin genotype was determined in DNA samples from all study subjects.

RESULTS

As zonulin concentrations did not correlate to the haptoglobin genotypes, we investigated the specificity of the zonulin ELISA assay using antibody capture experiments, mass spectrometry, and Western blot analysis. Using serum samples that gave the highest or lowest ELISA signals, we detected several proteins that are likely to be captured by the antibody in the present kit. However, none of these proteins corresponds to pre-haptoglobin2. We used increasing concentrations of recombinant pre-haptoglobin2 and complement C3 as one of the representative captured proteins and the ELISA kit did not detect either. Western blot analysis using both the polyclonal antibodies used in this kit and monoclonal antibodies rose against zonulin showed a similar protein recognition pattern but with different intensity of detection. The protein(s) measured using the ELISA kit was (were) significantly increased in patients with diabetes and obesity and correlated strongly with markers of the lipid and glucose metabolism. Combining mass spectrometry and Western blot analysis using the polyclonal antibodies used in the ELISA kit, we identified properdin as another member of the zonulin family.

CONCLUSION

Our study suggests that the zonulin ELISA does not recognize pre-haptoglobin2, rather structural (and possibly functional) analog proteins belonging to the mannose-associated serine protease family, with properdin being the most likely possible candidate.

On the impact of relatedness on SNP association analysis

BACKGROUND

When testing for SNP (single nucleotide polymorphism) associations in related individuals, observations are not independent. Simple linear regression assuming independent normally distributed residuals results in an increased type I error and the power of the test is also affected in a more complicate manner. Inflation of type I error is often successfully corrected by genomic control. However, this reduces the power of the test when relatedness is of concern. In the present paper, we derive explicit formulae to investigate how heritability and strength of relatedness contribute to variance inflation of the effect estimate of the linear model. Further, we study the consequences of variance inflation on hypothesis testing and compare the results with those of genomic control correction. We apply the developed theory to the publicly available HapMap trio data (N=129), the Sorbs (a self-contained population with N=977 characterised by a cryptic relatedness structure) and synthetic family studies with different sample sizes (ranging from N=129 to N=999) and different degrees of relatedness.

RESULTS

We derive explicit and easily to apply approximation formulae to estimate the impact of relatedness on the variance of the effect estimate of the linear regression model. Variance inflation increases with increasing heritability. Relatedness structure also impacts the degree of variance inflation as shown for example family structures. Variance inflation is smallest for HapMap trios, followed by a synthetic family study corresponding to the trio data but with larger sample size than HapMap. Next strongest inflation is observed for the Sorbs, and finally, for a synthetic family study with a more extreme relatedness structure but with similar sample size as the Sorbs. Type I error increases rapidly with increasing inflation. However, for smaller significance levels, power increases with increasing inflation while the opposite holds for larger significance levels. When genomic control is applied, type I error is preserved while power decreases rapidly with increasing variance inflation.

CONCLUSIONS

Stronger relatedness as well as higher heritability result in increased variance of the effect estimate of simple linear regression analysis. While type I error rates are generally inflated, the behaviour of power is more complex since power can be increased or reduced in dependence on relatedness and the heritability of the phenotype. Genomic control cannot be recommended to deal with inflation due to relatedness. Although it preserves type I error, the loss in power can be considerable. We provide a simple formula for estimating variance inflation given the relatedness structure and the heritability of a trait of interest. As a rule of thumb, variance inflation below 1.05 does not require correction and simple linear regression analysis is still appropriate.

The Origins of Ashkenaz, Ashkenazic Jews, and Yiddish

Recently, the geographical origins of Ashkenazic Jews (AJs) and their native language Yiddish were investigated by applying the Geographic Population Structure (GPS) to a cohort of exclusively Yiddish-speaking and multilingual AJs. GPS localized most AJs along major ancient trade routes in northeastern Turkey adjacent to primeval villages with names that resemble the word "Ashkenaz." These findings were compatible with the hypothesis of an Irano-Turko-Slavic origin for AJs and a Slavic origin for Yiddish and at odds with the Rhineland hypothesis advocating a Levantine origin for AJs and German origins for Yiddish. We discuss how these findings advance three ongoing debates concerning (1) the historical meaning of the term "Ashkenaz;" (2) the genetic structure of AJs and their geographical origins as inferred from multiple studies employing both modern and ancient DNA and original ancient DNA analyses; and (3) the development of Yiddish. We provide additional validation to the non-Levantine origin of AJs using ancient DNA from the Near East and the Levant. Due to the rising popularity of geo-localization tools to address questions of origin, we briefly discuss the advantages and limitations of popular tools with focus on the GPS approach. Our results reinforce the non-Levantine origins of AJs.

Forensic genetic analyses in isolated populations with examples of central European Valachs and Roma

Isolated populations present a constant threat to the correctness of forensic genetic casework. In this review article we present several examples of how analyzing samples from isolated populations can bias the results of the forensic statistics and analyses. We select our examples from isolated populations from central and southeastern Europe, namely the Valachs and the European Roma. We also provide the reader with general strategies and principles to improve the laboratory practice (best practice) and reporting of samples from supposedly isolated populations. These include reporting the precise population data used for computing the forensic statistics, using the appropriate θ correction factor for calculating allele frequencies, typing ancestry informative markers in samples of unknown or uncertain ethnicity and establishing ethnic-specific forensic databases.

Overcoming the dichotomy between open and isolated populations using genomic data from a large European dataset

Human populations are often dichotomized into “isolated” and “open” categories using cultural and/or geographical barriers to gene flow as differential criteria. Although widespread, the use of these alternative categories could obscure further heterogeneity due to inter-population differences in effective size, growth rate, and timing or amount of gene flow. We compared intra and inter-population variation measures combining novel and literature data relative to 87,818 autosomal SNPs in 14 open populations and 10 geographic and/or linguistic European isolates. Patterns of intra-population diversity were found to vary considerably more among isolates, probably due to differential levels of drift and inbreeding. The relatively large effective size estimated for some population isolates challenges the generalized view that they originate from small founding groups. Principal component scores based on measures of intra-population variation of isolated and open populations were found to be distributed along a continuum, with an area of intersection between the two groups. Patterns of inter-population diversity were even closer, as we were able to detect some differences between population groups only for a few multidimensional scaling dimensions. Therefore, different lines of evidence suggest that dichotomizing human populations into open and isolated groups fails to capture the actual relations among their genomic features.

The Slavic NBN Founder Mutation: A Role for Reproductive Fitness?

The vast majority of patients with Nijmegen Breakage Syndrome (NBS) are of Slavic origin and carry a deleterious deletion (c.657del5; rs587776650) in the NBN gene on chromosome 8q21. This mutation is essentially confined to Slavic populations and may thus be considered a Slavic founder mutation. Notably, not a single parenthood of a homozygous c.657del5 carrier has been reported to date, while heterozygous carriers do reproduce but have an increased cancer risk. These observations seem to conflict with the considerable carrier frequency of c.657del5 of 0.5% to 1% as observed in different Slavic populations because deleterious mutations would be eliminated quite rapidly by purifying selection. Therefore, we propose that heterozygous c.657del5 carriers have increased reproductive success, i.e., that the mutation confers heterozygote advantage. In fact, in our cohort study of the reproductive history of 24 NBS pedigrees from the Czech Republic, we observed that female carriers gave birth to more children on average than female non-carriers, while no such reproductive differences were observed for males. We also estimate that c.657del5 likely occurred less than 300 generations ago, thus supporting the view that the original mutation predated the historic split and subsequent spread of the ‘Slavic people’. We surmise that the higher fertility of female c.657del5 carriers reflects a lower miscarriage rate in these women, thereby reflecting the role of the NBN gene product, nibrin, in the repair of DNA double strand breaks and their processing in immune gene rearrangements, telomere maintenance, and meiotic recombination, akin to the previously described role of the DNA repair genes BRCA1 and BRCA2.

Evaluating mtDNA patterns of genetic isolation using a re-sampling procedure: A case study on Italian populations

BACKGROUND

A number of studies which have investigated isolation patterns in human populations rely on the analysis of intra- and inter-population genetic statistics of mtDNA polymorphisms. However, this approach makes it difficult to differentiate between the effects of long-term genetic isolation and the random fluctuations of statistics due to reduced sample size.

AIM

To overcome the confounding effect of sample size when detecting signatures of genetic isolation.

SUBJECTS AND METHODS

A re-sampling based procedure was employed to evaluate reduction in intra-population diversity, departure from surrounding genetic background and demographic stationarity in 34 Italian populations subject to isolation factors.

RESULTS

Signatures of genetic isolation were detected for all three statistics in seven populations: Pusteria valley, Sappada, Sauris, Timau settled in the eastern Italian Alps and Cappadocia, Filettino and Vallepietra settled in the Appenines. However, this study was unable to find signals for any of the statistics analysed in 19 populations. Finally, eight populations showing signals of isolation were found for one or two statistics.

CONCLUSION

The analysis revealed that the use of population genetic statistics combined with re-sampling procedure can help detect signatures of genetic isolation in human populations, even using a single, although highly informative, locus like mtDNA.

Coevolution of genes and languages and high levels of population structure among the highland populations of Daghestan

As a result of the combination of great linguistic and cultural diversity, the highland populations of Daghestan present an excellent opportunity to test the hypothesis of language-gene coevolution at a fine geographic scale. However, previous genetic studies generally have been restricted to uniparental markers and have not included many of the key populations of the region. To improve our understanding of the genetic structure of Daghestani populations and to investigate possible correlations between genetic and linguistic variation, we analyzed ~550,000 autosomal single nucleotide polymorphisms, phylogenetically informative Y chromosome markers and mtDNA haplotypes in 21 ethnic Daghestani groups. We found high levels of population structure in Daghestan consistent with the hypothesis of long-term isolation among populations of the highland Caucasus. Highland Daghestani populations exhibit extremely high levels of between-population diversity for all genetic systems tested, leading to some of the highest FST values observed for any region of the world. In addition, we find a significant positive correlation between gene and language diversity, suggesting that these two aspects of human diversity have coevolved as a result of historical patterns of social interaction among highland farmers at the community level. Finally, our data are consistent with the hypothesis that most Daghestanian-speaking groups descend from a common ancestral population (~6000-6500 years ago) that spread to the Caucasus by demic diffusion followed by population fragmentation and low levels of gene flow.

Integration of Genome-Wide SNP Data and Gene-Expression Profiles Reveals Six Novel Loci and Regulatory Mechanisms for Amino Acids and Acylcarnitines in Whole Blood

Profiling amino acids and acylcarnitines in whole blood spots is a powerful tool in the laboratory diagnosis of several inborn errors of metabolism. Emerging data suggests that altered blood levels of amino acids and acylcarnitines are also associated with common metabolic diseases in adults. Thus, the identification of common genetic determinants for blood metabolites might shed light on pathways contributing to human physiology and common diseases. We applied a targeted mass-spectrometry-based method to analyze whole blood concentrations of 96 amino acids, acylcarnitines and pathway associated metabolite ratios in a Central European cohort of 2,107 adults and performed genome-wide association (GWA) to identify genetic modifiers of metabolite concentrations. We discovered and replicated six novel loci associated with blood levels of total acylcarnitine, arginine (both on chromosome 6; rs12210538, rs17657775), propionylcarnitine (chromosome 10; rs12779637), 2-hydroxyisovalerylcarnitine (chromosome 21; rs1571700), stearoylcarnitine (chromosome 1; rs3811444), and aspartic acid traits (chromosome 8; rs750472). Based on an integrative analysis of expression quantitative trait loci in blood mononuclear cells and correlations between gene expressions and metabolite levels, we provide evidence for putative causative genes: SLC22A16 for total acylcarnitines, ARG1 for arginine, HLCS for 2-hydroxyisovalerylcarnitine, JAM3 for stearoylcarnitine via a trans-effect at chromosome 1, and PPP1R16A for aspartic acid traits. Further, we report replication and provide additional functional evidence for ten loci that have previously been published for metabolites measured in plasma, serum or urine. In conclusion, our integrative analysis of SNP, gene-expression and metabolite data points to novel genetic factors that may be involved in the regulation of human metabolism. At several loci, we provide evidence for metabolite regulation via gene-expression and observed overlaps with GWAS loci for common diseases. These results form a strong rationale for subsequent functional and disease-related studies.

Genetic Heritage of the Balto-Slavic Speaking Populations: A Synthesis of Autosomal, Mitochondrial and Y-Chromosomal Data

The Slavic branch of the Balto-Slavic sub-family of Indo-European languages underwent rapid divergence as a result of the spatial expansion of its speakers from Central-East Europe, in early medieval times. This expansion–mainly to East Europe and the northern Balkans–resulted in the incorporation of genetic components from numerous autochthonous populations into the Slavic gene pools. Here, we characterize genetic variation in all extant ethnic groups speaking Balto-Slavic languages by analyzing mitochondrial DNA (n = 6,876), Y-chromosomes (n = 6,079) and genome-wide SNP profiles (n = 296), within the context of other European populations. We also reassess the phylogeny of Slavic languages within the Balto-Slavic branch of Indo-European. We find that genetic distances among Balto-Slavic populations, based on autosomal and Y-chromosomal loci, show a high correlation (0.9) both with each other and with geography, but a slightly lower correlation (0.7) with mitochondrial DNA and linguistic affiliation. The data suggest that genetic diversity of the present-day Slavs was predominantly shaped in situ, and we detect two different substrata: ‘central-east European’ for West and East Slavs, and ‘south-east European’ for South Slavs. A pattern of distribution of segments identical by descent between groups of East-West and South Slavs suggests shared ancestry or a modest gene flow between those two groups, which might derive from the historic spread of Slavic people.

Genetic variants in AKR1B10 associate with human eating behavior

BACKGROUND

The human Aldoketoreductase 1B10 gene (AKR1B10) encodes one of the enzymes belonging to the family of aldoketoreductases and may be involved in detoxification of nutrients during digestion. Further, AKR1B10 mRNA (messenger ribonucleic acid) expression was diminished in brain regions potentially involved in the regulation of eating behavior in rats which are more sensitive to cocaine and alcohol. We hypothesized that the human AKR1B10 gene may also play a role in the regulation of human eating behavior.

RESULTS

We investigated the effects of 5 genetic variants of AKR1B10 on human eating behavior among 548 subjects from a German self-contained population, the Sorbs, and in 350 subjects from another independent German cohort. Among the Sorbs, we observed nominal associations with disinhibition at the 5' untranslated region (5' UTR) variant rs10232478 and the intragenic variants rs1834150 and rs782881 (all P ≤ 0.05). Further, we detected a relationship of rs1834150 and rs782881 with waist, smoking consumption (rs782881) and coffee consumption (rs1834150) (all P ≤ 0.05). Albeit non-significant, replication analyses revealed similar effect directions for disinhibition at rs1834150 (combined P = 0.0096). Moreover, in the replication cohort we found rs1834150 related to increased restraint scores with a similar direction as in the Sorbs (combined P = 0.0072).

CONCLUSION

Our data suggest that genetic variants in the AKR1B10 locus may influence human eating behavior.

Indications for potential parent-of-origin effects within the FTO gene

Genome-Wide Association Studies (GWAS) were successfully applied to discover associations with obesity. However, the GWAS design is usually based on unrelated individuals and inheritance information on the parental origin of the alleles is missing. Taking into account parent-of-origin may provide further insights into the genetic mechanisms contributing to obesity. We hypothesized that there may be variants within the robustly replicated fat mass and obesity associated (FTO) gene that may confer different risk for obesity depending on transmission from mother or father. Genome-wide genotypes and pedigree information from the Sorbs population were used. Phased genotypes among 525 individuals were generated by AlphaImpute. Subsequently, 22 SNPs within FTO introns 1 to 3 were selected and parent-of-origin specific association analyses were performed using PLINK. Interestingly, we identified several SNPs conferring different genetic effects (P≤0.05) depending on parental origin—among them, rs1861868, rs1121980 and rs9939973 (all in intron 1). To confirm our findings, we investigated the selected variants in 705 German trios comprising an (extremely) obese child or adolescent and both parents. Again, we observed evidence for POE effects in intron 2 and 3 (P≤0.05) as indicated by the parental asymmetry test. Our results suggest that the obesity risk transmitted by several FTO variants may depend on the parental origin of the allele. Larger family-based studies are warranted to replicate our findings.

Extensive genome-wide autozygosity in the population isolates of Daghestan

Isolated populations are valuable resources for mapping disease genes, as inbreeding increases genome-wide homozygosity and enhances the ability to map disease alleles on a genetically uniform background within a relatively homogenous environment. The populations of Daghestan are thought to have resided in the Caucasus Mountains for hundreds of generations and are characterized by a high prevalence of certain complex diseases. To explore the extent to which their unique population history led to increased levels of inbreeding, we genotyped >550 000 autosomal single-nucleotide polymorphisms (SNPs) in a set of 14 population isolates speaking Nakh-Daghestanian (ND) languages. The ND-speaking populations showed greatly elevated coefficients of inbreeding, very high numbers and long lengths of Runs of Homozygosity, and elevated linkage disequilibrium compared with surrounding groups from the Caucasus, the Near East, Europe, Central and South Asia. These results are consistent with the hypothesis that most ND-speaking groups descend from a common ancestral population that fragmented into a series of genetic isolates in the Daghestanian highlands. They have subsequently maintained a long-term small effective population size as a result of constant inbreeding and very low levels of gene flow. Given these findings, Daghestanian population isolates are likely to be useful for mapping genes associated with complex diseases.

Novel heterozygous IGF1R mutation in two brothers with developing impaired glucose tolerance

Infants born small for gestational age (SGA) are at risk to develop metabolic complications. Insulin-like growth factor 1 (IGF-1) resistance due to IGF-1 receptor (IGF1R) mutations is a rare genetic condition that causes proportionate growth retardation. The contribution of an impaired IGF1R function to the development of comorbidities such as disturbed glucose homeostasis is not well understood. Genetic analysis and detailed auxological, endocrine and psychological investigations in two male SGA siblings were performed. The two patients and their father bear a novel heterozygous mutation (p.Cys1248Tyr) in the IGF1R gene. Both brothers displayed very similar growth pattern before and during recombinant human growth hormone treatment, whereas oral glucose tolerance tests showed variable manifestations of progressive impaired glucose tolerance. The father had already developed type 2 diabetes mellitus. Growth retardation in our patients is likely caused by the IGF1R mutation that might predispose to disturbances of carbohydrate homeostasis. Therefore, a close metabolic monitoring of affected patients is indicated, particularly if growth hormone therapy is commenced.

Role of genetic variants in ADIPOQ in human eating behavior

The beneficial effects of adiponectin and its negative correlation with BMI are well described. Adiponectin serum levels are altered in eating disorders such as anorexia nervosa, bulimia nervosa or binge eating. Here, we tested the hypothesis that (1) adiponectin serum levels correlate with human eating behavior factors and (2) that genetic variants of the ADIPOQ locus influence both serum levels and eating behavior. We analyzed 11 SNPs within ADIPOQ and in the 5' UTR and measured serum adiponectin levels in 1,036 individuals from the German Sorbs population. The German version of the three-factor eating questionnaire (FEV) was completed by 548 Sorbs. For replication purposes, we included an independent replication cohort from Germany (N = 350). In the Sorbs, we observed positive correlations of restraint with adiponectin serum levels (P = 0.001; r = 0.148) which, however, did not withstand adjustment for covariates (P = 0.083; r = 0.077). In addition, four SNPs were nominally associated with serum adiponectin levels (all P < 0.05). Of these, two variants (rs3774261; rs1501229, all P < 0.05) were also related to disinhibition. Furthermore, three variants were exclusively associated with hunger (rs2036373, P = 0.049) and disinhibition (rs822396; rs864265, all P < 0.05). However, none of these associations withstood Bonferroni corrections for multiple testing (all P > 9.3 × 10(-4)). In our replication cohort, we observed similar effect directions at rs1501229 for disinhibition and hunger. A meta-analysis resulted in nominal statistical significance P = 0.036 (Z score 2.086) and P = 0.017 (Z score 2.366), respectively. Given the observed relationship of the SNPs with adiponectin levels and eating behavior, our data support a potential role of adiponectin in human eating behavior. Whether the relationship with eating behavior is mediated by the effects of circulating adiponectin warrants further investigations.

Genome wide meta-analysis highlights the role of genetic variation in RARRES2 in the regulation of circulating serum chemerin

Chemerin is an adipokine proposed to link obesity and chronic inflammation of adipose tissue. Genetic factors determining chemerin release from adipose tissue are yet unknown. We conducted a meta-analysis of genome-wide association studies (GWAS) for serum chemerin in three independent cohorts from Europe: Sorbs and KORA from Germany and PPP-Botnia from Finland (total N = 2,791). In addition, we measured mRNA expression of genes within the associated loci in peripheral mononuclear cells by micro-arrays, and within adipose tissue by quantitative RT-PCR and performed mRNA expression quantitative trait and expression-chemerin association studies to functionally substantiate our loci. Heritability estimate of circulating chemerin levels was 16.2% in the Sorbs cohort. Thirty single nucleotide polymorphisms (SNPs) at chromosome 7 within the retinoic acid receptor responder 2 (RARRES2)/Leucine Rich Repeat Containing (LRRC61) locus reached genome-wide significance (p<5.0×10⁻⁸) in the meta-analysis (the strongest evidence for association at rs7806429 with p = 7.8×10⁻¹⁴, beta = −0.067, explained variance 2.0%). All other SNPs within the cluster were in linkage disequilibrium with rs7806429 (minimum r² = 0.43 in the Sorbs cohort). The results of the subgroup analyses of males and females were consistent with the results found in the total cohort. No significant SNP-sex interaction was observed. rs7806429 was associated with mRNA expression of RARRES2 in visceral adipose tissue in women (p<0.05 after adjusting for age and body mass index). In conclusion, the present meta-GWAS combined with mRNA expression studies highlights the role of genetic variation in the RARRES2 locus in the regulation of circulating chemerin concentrations.

Chemerin is an adipokine proposed to link obesity and chronic inflammation of adipose tissue. In the present study we show that circulating chemerin is a heritable trait. In a meta-analysis of genome-wide association studies (GWAS) of 2,791 individuals from Germany and Finland, we identified common genetic variants which associate with serum chemerin levels. The variants map within the retinoic acid receptor responder 2 (RARRES2)/Leucine Rich Repeat Containing (LRRC61) at chromosome 7. To better understand the potential functionality of the identified variants, we also provide insights into the mRNA expression of RARRES2 (encoding chemerin) in blood and adipose tissue. Our results highlight the role and function of genetic variation in the RARRES2 locus in the regulation of circulating chemerin concentrations.

Global DNA methylation levels in human adipose tissue are related to fat distribution and glucose homeostasis

AIMS/HYPOTHESIS

Epigenetic alterations may influence the metabolic pathways involved in human obesity. We hypothesised that global DNA methylation levels in adipose tissue might be associated with obesity and related phenotypes.

METHODS

We measured global DNA methylation levels in paired samples of subcutaneous adipose tissue (SAT) and omental visceral adipose tissue (OVAT) from 51 individuals, and in leucocytes from 559 Sorbs, a population from Germany, using LUminometric Methylation Assay (LUMA). To further investigate the underlying mechanisms of the observed associations, we measured global methylation levels in 3T3-L1 adipocytes exposed to glucose, insulin and lipids.

RESULTS

Global methylation levels (±SD) were significantly higher in OVAT (74.27% ± 2.2%) compared with SAT (71.97% ± 2.4%; paired t test, p < 1 × 10(-9)). Furthermore, global methylation levels in SAT were positive correlates of measures of fat distribution (waist measurement, WHR) and glucose homeostasis (HbA1c) (all p < 0.015 after accounting for multiple testing and covariates). Global methylation levels in the German Sorb cohort were associated with glucose homeostasis, but this association did not withstand adjustment for covariates. Exposure of 3T3-L1 adipocytes to insulin, palmitate and glucose decreased global methylation levels 1 h after treatment relative to controls.

CONCLUSIONS/INTERPRETATION

Our data suggest that the variability in global methylation in adipose tissue might be related to alterations in glucose metabolism.

Demographic events and evolutionary forces shaping European genetic diversity

Europeans have been the focus of some of the largest studies of genetic diversity in any species to date. Recent genome-wide data have reinforced the hypothesis that present-day European genetic diversity is strongly correlated with geography. The remaining challenge now is to understand more precisely how patterns of diversity in Europe reflect ancient demographic events such as postglacial expansions or the spread of farming. It is likely that recent advances in paleogenetics will give us some of these answers. There has also been progress in identifying specific segments of European genomes that reflect adaptations to selective pressures from the physical environment, disease, and dietary shifts. A growing understanding of how modern European genetic diversity has been shaped by demographic and evolutionary forces is not only of basic historical and anthropological interest but also aids genetic studies of disease.

The impact of whole-genome sequencing on the reconstruction of human population history

Examining patterns of molecular genetic variation in both modern-day and ancient humans has proved to be a powerful approach to learn about our origins. Rapid advances in DNA sequencing technology have allowed us to characterize increasing amounts of genomic information. Although this clearly provides unprecedented power for inference, it also introduces more complexity into the way we use and interpret such data. Here, we review ongoing debates that have been influenced by improvements in our ability to sequence DNA and discuss some of the analytical challenges that need to be overcome in order to fully exploit the rich historical information that is contained in the entirety of the human genome.

TAS2R38 and its influence on smoking behavior and glucose homeostasis in the German Sorbs

Background

Genetic variants within the bitter taste receptor gene TAS2R38 are associated with sensitivity to bitter taste and are related to eating behavior in the Amish population. Sensitivity to bitter taste is further related to anthropometric traits in an genetically isolated Italian population. We tested whether the TAS2R38 variants (rs713598; rs1726866 and rs10246939) may be related to eating behavior, anthropometric parameters, metabolic traits and consumer goods intake in the German Sorbs.

Materials and Methods

The three SNPs were genotyped in a total cohort of 1007 individuals (male/female: 405/602). The German version of the three-factor eating questionnaire was completed by 548 individuals. Genetic association analyses for smoking behavior, alcohol and coffee intake, eating behavior factors (restraint, disinhibition and hunger) and other metabolic traits were analyzed. Further, by combining the three SNPs we applied comparative haplotype analyses categorizing PAV (proline-alanine-valine) carriers (tasters) vs. homozygous AVI (alanin-valine-isoleucine) carriers (non-tasters).

Results

Significant associations of genetic variants within TAS2R38 were identified with percentage of body fat, which were driven by associations in women. In men, we observed significant associations with 30 min plasma glucose, and area under the curve for plasma glucose (0–120 min) (all adjusted P≤0.05). Further, we found that carriers of at least one PAV allele show significantly lower cigarette smoking per day (P = 0.002) as well as, albeit non-significant, lower alcohol intake. We did not confirm previously reported associations between genetic variants of TAS2R38 and eating behavior.

Conclusion

Our data suggest that genetic variation in TAS2R38 is related to individual body composition measures and may further influence consumer goods intake in the Sorbs possibly via individual sensitivity to bitter taste.

Genetic substructure of Kuwaiti population reveals migration history

The State of Kuwait is characterized by settlers from Saudi Arabia, Iran, and other regions of the Arabian Peninsula. The settlements and subsequent admixtures have shaped the genetics of Kuwait. High prevalence of recessive disorders and metabolic syndromes (that increase risk of diabetes) is seen in the peninsula. Understanding the genetic structure of its population will aid studies designed to decipher the underlying causes of these disorders. In this study, we analyzed 572,366 SNP markers from 273 Kuwaiti natives genotyped using the illumina HumanOmniExpress BeadChip. Model-based clustering identified three genetic subgroups with different levels of admixture. A high level of concordance (Mantel test, p=0.0001 for 9999 repeats) was observed between the derived genetic clusters and the surname-based ancestries. Use of Human Genome Diversity Project (HGDP) data to understand admixtures in each group reveals the following: the first group (Kuwait P) is largely of West Asian ancestry, representing Persians with European admixture; the second group (Kuwait S) is predominantly of city-dwelling Saudi Arabian tribe ancestry, and the third group (Kuwait B) includes most of the tent-dwelling Bedouin surnames and is characterized by the presence of 17% African ancestry. Identity by Descent and Homozygosity analyses find Kuwait's population to be heterogeneous (placed between populations that have large amount of ROH and the ones with low ROH) with Kuwait S as highly endogamous, and Kuwait B as diverse. Population differentiation FST estimates place Kuwait P near Asian populations, Kuwait S near Negev Bedouin tribes, and Kuwait B near the Mozabite population. FST distances between the groups are in the range of 0.005 to 0.008; distances of this magnitude are known to cause false positives in disease association studies. Results of analysis for genetic features such as linkage disequilibrium decay patterns conform to Kuwait's geographical location at the nexus of Africa, Europe, and Asia.

THOC5: a novel gene involved in HDL-cholesterol metabolism

Although numerous genes are known to regulate serum lipid traits, identified variants explain only a small proportion of the expected heritability. We intended to identify further genetic variants associated with lipid phenotypes in a self-contained population of Sorbs in Germany. We performed a genome-wide association study (GWAS) on LDL-cholesterol, HDL-cholesterol (HDL-C), and triglyceride (TG) levels in 839 Sorbs. All single-nucleotide polymorphisms with a P value <0.01 were subjected to a meta-analysis, including an independent Swedish cohort (Diabetes Genetics Initiative; n = ∼3,100). Novel association signals with the strongest effects were subjected to replication studies in an additional German cohort (Berlin, n = 2,031). In the initial GWAS in the Sorbs, we identified 14 loci associated with lipid phenotypes reaching P values <10⁻⁵ and confirmed significant effects for 18 previously reported loci. The combined meta-analysis of the three study cohorts (n(HDL) = 6041; n(LDL) = 5,995; n(TG) = 6,087) revealed a novel association for a variant in THOC5 (rs8135828) with serum HDL-C levels (P = 1.78 × 10⁻⁷; Z-score = -5.221). Consistently, the variant was also associated with circulating APOA1 levels in Sorbs. The small interfering RNA-mediated mRNA silencing of THOC5 in HepG2 cells resulted in lower mRNA levels of APOA1, SCARB1, and ABCG8 (all P < 0.05). We propose THOC5 to be a novel gene involved in the regulation of serum HDL-C levels.

The role of rs2237781 within GRM8 in eating behavior

INTRODUCTION

The glutamate receptor, metabotropic 8 gene (GRM8) encodes a G-protein-coupled glutamate receptor and has been associated with smoking behavior and liability to alcoholism implying a role in addiction vulnerability. Data from animal studies suggest that GRM8 may be involved in the regulation of the neuropeptide Y and melanocortin pathways and might influence food intake and metabolism. This study aimed to investigate the effects of the genetic variant rs2237781 within GRM8 on human eating behavior.

METHODS

The initial analysis included 548 Sorbs from Germany who have been extensively phenotyped for metabolic traits and who completed the German version of the three-factor eating questionnaire. In addition, we analyzed two independent sample sets comprising 293 subjects from another German cohort and 430 Old Order Amish individuals. Genetic associations with restraint, disinhibition, and hunger were assessed in an additive linear regression model.

RESULTS

Among the Sorbs the major G allele of rs2237781 was significantly associated with increased restraint scores in eating behavior (P = 1.9 × 10(-4); β = +1.936). The German cohort and the Old Order Amish population revealed a trend in the same direction for restraint (P = 0.242; β = +0.874; P = 0.908; β = +0.096; respectively). A meta-analysis resulted in a combined P = 3.1 × 10(-3) (Z-score 2.948).

CONCLUSION

Our data suggest that rs2237781 within GRM8 may influence human eating behavior factors probably via pathways involved in addictive behavior.

Analysis of a rare functional truncating mutation rs61757459 in vaspin (SERPINA12) on circulating vaspin levels

A recent genome-wide association study suggests that genetic variation within the vaspin gene might contribute to the variability in circulating serum visceral adipose tissue-derived serine protease inhibitor (vaspin) concentrations. Here, we analyzed the functional consequences of the rare variant rs61757459 predicting a premature stop codon and its impact on circulating serum vaspin concentrations. In order to identify genetic variation, we sequenced the vaspin gene in 48 nonrelated Caucasian subjects. Rs61757459 was subsequently genotyped in three metabolically well-characterized German cohorts (N = 4,019). We addressed the impact of rs61757459 on the crystal structure of vaspin and investigated its effects on vaspin expression in vivo as well as in vitro using various cell lines (Escherichia coli, HEK293). Along with previously reported common genetic variants, sequencing of vaspin revealed a rare variant (rs61757459; minor allele frequency: 1 %) which predicts a premature stop codon p.R211X. Heterozygous carriers of this mutation had lower circulating vaspin levels when compared with noncarriers. In silico structure analysis of the truncated vaspin, which was estimated to be 24.5 kDa, suggested misfolding and potential instability due to the absence of core structural domains. Indeed, the truncated protein was detected after recombinant expression in E. coli and in lysate, but not in supernatant of HEK293 cells. We conclude that rs61757459 is a functional mutation that results in a truncated protein whose instability likely results in reduced serum vaspin levels.

KEY MESSAGE

A rare variant (rs61757459) in vaspin coding for the stop codon p.R211X is related to lower circulating vaspin concentrations. Structure analysis suggests misfolding and instability due to the absence of core structural domains. The truncated protein is detectable after recombinant expression in E. coli and in lysate, but not in supernatant of HEK293 cells.

Role of vaspin in human eating behaviour

Objective

The adipokine vaspin (visceral adipose tissue derived serine protease inhibitor, serpinA12) follows a meal-related diurnal variation in humans and intracerebroventricular vaspin administration leads to acutely reduced food intake in db/db mice. We therefore hypothesized that vaspin may play a role in human eating behaviour.

Materials and Methods

We measured serum vaspin concentrations in 548 subjects from a self-contained population of Sorbs (Germany) who underwent detailed metabolic testing including eating behaviour assessments using the three-factor eating questionnaire. In addition, genetic variation within vaspin was assessed by genotyping 28 single nucleotide polymorphisms (SNPs) in all study subjects.

Results

Serum vaspin concentrations correlated positively with restraint, disinhibition and hunger (all P<0.05), although the correlations did not withstand further adjustments for age, gender and BMI (all P>0.05). Independent of observed correlations, genetic variants in vaspin were associated with serum vaspin levels but showed no significant association with any of the eating behaviour phenotypes after accounting for multiple testing (P≥0.05 after adjusting for age, gender and BMI).

Conclusion

Our data suggest that serum vaspin concentrations might modulate human eating behaviour, which does not seem to be affected by common genetic variation in vaspin.

Genetic characterization of northeastern Italian population isolates in the context of broader European genetic diversity

Population genetic studies on European populations have highlighted Italy as one of genetically most diverse regions. This is possibly due to the country's complex demographic history and large variability in terrain throughout the territory. This is the reason why Italy is enriched for population isolates, Sardinia being the best-known example. As the population isolates have a great potential in disease-causing genetic variants identification, we aimed to genetically characterize a region from northeastern Italy, which is known for isolated communities. Total of 1310 samples, collected from six geographically isolated villages, were genotyped at >145 000 single-nucleotide polymorphism positions. Newly genotyped data were analyzed jointly with the available genome-wide data sets of individuals of European descent, including several population isolates. Despite the linguistic differences and geographical isolation the village populations still show the greatest genetic similarity to other Italian samples. The genetic isolation and small effective population size of the village populations is manifested by higher levels of genomic homozygosity and elevated linkage disequilibrium. These estimates become even more striking when the detected substructure is taken into account. The observed level of genetic isolation in Friuli-Venezia Giulia region is more extreme according to several measures of isolation compared with Sardinians, French Basques and northern Finns, thus proving the status of an isolate.

Contemporary paternal genetic landscape of Polish and German populations: from early medieval Slavic expansion to post-World War II resettlements

Homogeneous Proto-Slavic genetic substrate and/or extensive mixing after World War II were suggested to explain homogeneity of contemporary Polish paternal lineages. Alternatively, Polish local populations might have displayed pre-war genetic heterogeneity owing to genetic drift and/or gene flow with neighbouring populations. Although sharp genetic discontinuity along the political border between Poland and Germany indisputably results from war-mediated resettlements and homogenisation, it remained unknown whether Y-chromosomal diversity in ethnically/linguistically defined populations was clinal or discontinuous before the war. In order to answer these questions and elucidate early Slavic migrations, 1156 individuals from several Slavic and German populations were analysed, including Polish pre-war regional populations and an autochthonous Slavic population from Germany. Y chromosomes were assigned to 39 haplogroups and genotyped for 19 STRs. Genetic distances revealed similar degree of differentiation of Slavic-speaking pre-war populations from German populations irrespective of duration and intensity of contacts with German speakers. Admixture estimates showed minor Slavic paternal ancestry (~20%) in modern eastern Germans and hardly detectable German paternal ancestry in Slavs neighbouring German populations for centuries. BATWING analysis of isolated Slavic populations revealed that their divergence was preceded by rapid demographic growth, undermining theory that Slavic expansion was primarily linguistic rather than population spread. Polish pre-war regional populations showed within-group heterogeneity and lower STR variation within R-M17 subclades compared with modern populations, which might have been homogenised by war resettlements. Our results suggest that genetic studies on early human history in the Vistula and Oder basins should rely on reconstructed pre-war rather than modern populations.

The population genetics of Quechuas, the largest native South American group: autosomal sequences, SNPs, and microsatellites evidence high level of diversity

Elucidating the pattern of genetic diversity for non-European populations is necessary to make the benefits of human genetics research available to individuals from these groups. In the era of large human genomic initiatives, Native American populations have been neglected, in particular, the Quechua, the largest South Amerindian group settled along the Andes. We characterized the genetic diversity of a Quechua population in a global setting, using autosomal noncoding sequences (nine unlinked loci for a total of 16 kb), 351 unlinked SNPs and 678 microsatellites and tested predictions of the model of the evolution of Native Americans proposed by (Tarazona-Santos et al.: Am J Hum Genet 68 (2001) 1485-1496). European admixture is <5% and African ancestry is barely detectable in the studied population. The largest genetic distances were between African versus Quechua or Melanesian populations, which is concordant with the African origin of modern humans and the fact that South America was the last part of the world to be peopled. The diversity in the Quechua population is comparable with that of Eurasian populations, and the allele frequency spectrum based on resequencing data does not reflect a reduction in the proportion of rare alleles. Thus, the Quechua population is a large reservoir of common and rare genetic variants of South Amerindians. These results are consistent with and complement our evolutionary model of South Amerindians (Tarazona-Santos et al.: Am J Hum Genet 68 (2001) 1485-1496), proposed based on Y-chromosome data, which predicts high genomic diversity due to the high level of gene flow between Andean populations and their long-term effective population size.

Association and evolutionary studies of the melatonin receptor 1B gene (MTNR1B) in the self-contained population of Sorbs from Germany

AIMS

Several polymorphisms of the melatonin receptor 1B gene (MTNR1B) have been shown to be associated with elevated fasting plasma glucose and impaired early insulin release. The aim of this study was to assess the effects of MTNR1B variants on traits related to the metabolic syndrome in the self-contained population of Sorbs from Germany. As comprehensive studies concerning the conservation of MTNR1B are lacking, we also evaluated natural selection in vertebrates and human populations at this locus.

METHODS

Five single nucleotide polymorphisms representing all blocks of linkage disequilibrium within and surrounding the MTNR1B locus were genotyped in 937 Sorbs for association analyses on metabolic traits related to Type 2 diabetes. The associations were assessed by regression analyses, the conservation between species was investigated with phylogenetic analysis by maximum likelihood (PAML). In addition, various tests of population genetic measures (e.g. fixation index, Tajima's D) were performed.

RESULTS

Previously reported association between MTNR1B variants (rs10830963, rs4753426) and oral glucose tolerance test-derived indices of β-cell function (homeostasis model assessment-B, P = 3.7 × 10⁻⁶ and P = 0.004, respectively), as well as insulin (fasting insulin: P=2×10⁻³ and P=0.02; 30-min insulin: P = 2.1 × 10⁻⁴ and P=0.03, respectively) and fasting glucose (rs10830963, P=1.2×10⁻⁶) parameters could be replicated in the present study. Phylogenetic analysis by maximum likelihood analyses showed that the gene was strongly conserved between species (ω=0.2583). Structures important for the receptor function are also conserved. On the lineage leading to human adaptive selection was present (ω=1.1030). Population genetic measures further indicated natural selection.

CONCLUSIONS

Our data support the physiologic importance of MTNR1B in the context of glucose homeostasis and suggest evidence of selection at this locus.

Perspectives on human population structure at the cusp of the sequencing era

Human groups show structured levels of genetic similarity as a consequence of factors such as geographical subdivision and genetic drift. Surveying this structure gives us a scientific perspective on human origins, sheds light on evolutionary processes that shape both human adaptation and disease, and is integral to effectively carrying out the mission of global medical genetics and personalized medicine. Surveys of population structure have been ongoing for decades, but in the past three years, single-nucleotide-polymorphism (SNP) array technology has provided unprecedented detail on human population structure at global and regional scales. These studies have confirmed well-known relationships between distantly related populations and uncovered previously unresolvable relationships among closely related human groups. SNPs represent the first dense genome-wide markers, and as such, their analysis has raised many challenges and insights relevant to the study of population genetics with whole-genome sequences. Here we draw on the lessons from these studies to anticipate the directions that will be most fruitful to pursue during the emerging whole-genome sequencing era.

Population-genetic comparison of the Sorbian isolate population in Germany with the German KORA population using genome-wide SNP arrays

Background

The Sorbs are an ethnic minority in Germany with putative genetic isolation, making the population interesting for disease mapping. A sample of N = 977 Sorbs is currently analysed in several genome-wide meta-analyses. Since genetic differences between populations are a major confounding factor in genetic meta-analyses, we compare the Sorbs with the German outbred population of the KORA F3 study (N = 1644) and other publically available European HapMap populations by population genetic means. We also aim to separate effects of over-sampling of families in the Sorbs sample from effects of genetic isolation and compare the power of genetic association studies between the samples.

Results

The degree of relatedness was significantly higher in the Sorbs. Principal components analysis revealed a west to east clustering of KORA individuals born in Germany, KORA individuals born in Poland or Czech Republic, Half-Sorbs (less than four Sorbian grandparents) and Full-Sorbs. The Sorbs cluster is nearest to the cluster of KORA individuals born in Poland. The number of rare SNPs is significantly higher in the Sorbs sample. FST between KORA and Sorbs is an order of magnitude higher than between different regions in Germany. Compared to the other populations, Sorbs show a higher proportion of individuals with runs of homozygosity between 2.5 Mb and 5 Mb. Linkage disequilibrium (LD) at longer range is also slightly increased but this has no effect on the power of association studies.

Oversampling of families in the Sorbs sample causes detectable bias regarding higher FST values and higher LD but the effect is an order of magnitude smaller than the observed differences between KORA and Sorbs. Relatedness in the Sorbs also influenced the power of uncorrected association analyses.

Conclusions

Sorbs show signs of genetic isolation which cannot be explained by over-sampling of relatives, but the effects are moderate in size. The Slavonic origin of the Sorbs is still genetically detectable.

Regarding LD structure, a clear advantage for genome-wide association studies cannot be deduced. The significant amount of cryptic relatedness in the Sorbs sample results in inflated variances of Beta-estimators which should be considered in genetic association analyses.