The genome-wide patterns of variation expose significant substructure in a founder population

Diverse ancestral representation improves genetic intolerance metrics

The unprecedented scale of genomic databases has revolutionized our ability to identify regions in the human genome intolerant to variation-regions often implicated in disease. However, these datasets remain constrained by limited ancestral diversity. Here, we analyze whole-exome sequencing data from 460,551 UK Biobank and 125,748 Genome Aggregation Database (gnomAD) participants across multiple ancestries to test several key intolerance metrics, including the Residual Variance Intolerance Score (RVIS), Missense Tolerance Ratio (MTR), and Loss-of-Function Observed/Expected ratio (LOF O/E). We demonstrate that increasing ancestral representation, rather than sample size alone, critically drives their performance. Scores trained on variation observed in African and Admixed American ancestral groups show higher resolution in detecting haploinsufficient and neurodevelopmental disease risk genes compared to scores trained on European ancestry groups. Most strikingly, MTR trained on 43,000 multi-ancestry exomes demonstrates greater predictive power than when trained on a nearly 10-fold larger dataset of 440,000 non-Finnish European exomes. We further find that European ancestry group-based scores are likely approaching saturation. These findings highlight the need for enhanced population representation in genomic resources to fully realize the potential of precision medicine and drug discovery. Ancestry group-specific scores are publicly available through an interactive portal: http://intolerance.public.cgr.astrazeneca.com/ .

Identification and functional characterisation of DNA methylation differences between East- and West-originating Finns

Eastern and Western Finns show a striking difference in coronary heart disease-related mortality; genetics is a known contributor for this discrepancy. Here, we discuss the potential role of DNA methylation in mediating the discrepancy in cardiometabolic disease-risk phenotypes between the sub-populations. We used data from the Young Finns Study (n = 969) to compare the genome-wide DNA methylation levels of East- and West-originating Finns. We identified 21 differentially methylated loci (FDR < 0.05; Δβ >2.5%) and 7 regions (smoothed FDR < 0.05; CpGs ≥ 5). Methylation at all loci and regions associates with genetic variants (p < 5 × 10-8). Independently of genetics, methylation at 11 loci and 4 regions associates with transcript expression, including genes encoding zinc finger proteins. Similarly, methylation at 5 loci and 4 regions associates with cardiometabolic disease-risk phenotypes including triglycerides, glucose, cholesterol, as well as insulin treatment. This analysis was also performed in LURIC (n = 2371), a German cardiovascular patient cohort, and results replicated for the association of methylation at cg26740318 and DMR_11p15 with diabetes-related phenotypes and methylation at DMR_22q13 with triglyceride levels. Our results indicate that DNA methylation differences between East and West Finns may have a functional role in mediating the cardiometabolic disease discrepancy between the sub-populations.

An abdominal obesity missense variant in the adipocyte thermogenesis gene TBX15 is implicated in adaptation to cold in Finns

Mechanisms of abdominal obesity GWAS variants have remained largely unknown. To elucidate these mechanisms, we leveraged subcutaneous adipose tissue (SAT) single nucleus RNA-sequencing and genomics data. After discovering that heritability of abdominal obesity is enriched in adipocytes, we focused on a SAT unique adipocyte marker gene, the transcription factor TBX15, and its abdominal obesity-associated deleterious missense variant, rs10494217. The allele frequency of rs10494217 revealed a north-to-south decreasing gradient, with consistent significant FST values observed for 25 different populations when compared to Finns, a population with a history of genetic isolation. Given the role of Tbx15 in mouse thermogenesis, the frequency may have increased as an adaptation to cold in Finns. Our selection analysis provided significant evidence of selection for the abdominal obesity risk allele T of rs10494217 in Finns, with a north-to-south decreasing trend in other populations, and demonstrated that latitude significantly predicts the allele frequency. We also discovered that the risk allele status significantly affects SAT adipocyte expression of multiple adipocyte marker genes in trans in two cohorts. Two of these trans genes have been connected to thermogenesis, supporting the thermogenic effect of the TBX15 missense variant as a possible cause of its selection. Adipose expression of one trans gene, a lncRNA, AC002066.1, was strongly associated with adipocyte size, implicating it in metabolically unhealthy adipocyte hypertrophy. In summary, the abdominal obesity variant rs10494217 was selected in Finns, and individuals with the risk allele have trans effects on adipocyte expression of genes relating to thermogenesis and adipocyte hypertrophy.

Exome sequencing identified rare recurrent copy number variants and hereditary breast cancer susceptibility

Copy number variants (CNVs) are a major source of genetic variation and can disrupt genes or affect gene dosage. They are known to be causal or underlie predisposition to various diseases. However, the role of CNVs in inherited breast cancer susceptibility has not been thoroughly investigated. To address this, we performed whole-exome sequencing based analysis of rare CNVs in 98 high-risk Northern Finnish breast cancer cases. After filtering, selected candidate alleles were validated and characterized with a combination of orthogonal methods, including PCR-based approaches, optical genome mapping and long-read sequencing. This revealed three recurrent alterations: a 31 kb deletion co-occurring with a retrotransposon insertion (delins) in RAD52, a 13.4 kb deletion in HSD17B14 and a 64 kb partial duplication of RAD51C. Notably, all these genes encode proteins involved in pathways previously identified as essential for breast cancer development. Variants were genotyped in geographically matched cases and controls (altogether 278 hereditary and 1983 unselected breast cancer cases, and 1229 controls). The RAD52 delins and HSD17B14 deletion both showed significant enrichment among cases with indications of hereditary disease susceptibility. RAD52 delins was identified in 7/278 cases (2.5%, P = 0.034, OR = 2.86, 95% CI = 1.10-7.45) and HSD17B14 deletion in 8/278 cases (2.9%, P = 0.014, OR = 3.28, 95% CI = 1.31-8.23), the frequency of both variants in the controls being 11/1229 (0.9%). This suggests a role for RAD52 and HSD17B14 in hereditary breast cancer susceptibility. The RAD51C duplication was very rare, identified only in 2/278 of hereditary cases and 2/1229 controls (P = 0.157, OR = 4.45, 95% CI = 0.62-31.70). The identification of recurrent CNVs in these genes, and especially the relatively high frequency of RAD52 and HSD17B14 alterations in the Finnish population, highlights the importance of studying CNVs alongside single nucleotide variants when searching for genetic factors underlying hereditary disease predisposition.

Human genetics uncovers MAP3K15 as an obesity-independent therapeutic target for diabetes

We performed collapsing analyses on 454,796 UK Biobank (UKB) exomes to detect gene-level associations with diabetes. Recessive carriers of nonsynonymous variants in MAP3K15 were 30% less likely to develop diabetes (P = 5.7 × 10-10) and had lower glycosylated hemoglobin (β = -0.14 SD units, P = 1.1 × 10-24). These associations were independent of body mass index, suggesting protection against insulin resistance even in the setting of obesity. We replicated these findings in 96,811 Admixed Americans in the Mexico City Prospective Study (P < 0.05)Moreover, the protective effect of MAP3K15 variants was stronger in individuals who did not carry the Latino-enriched SLC16A11 risk haplotype (P = 6.0 × 10-4). Separately, we identified a Finnish-enriched MAP3K15 protein-truncating variant associated with decreased odds of both type 1 and type 2 diabetes (P < 0.05) in FinnGen. No adverse phenotypes were associated with protein-truncating MAP3K15 variants in the UKB, supporting this gene as a therapeutic target for diabetes.

Genomic analysis reveals geography rather than culture as the predominant factor shaping genetic variation in northern Kenyan human populations

OBJECTIVES

The aim of this study was to characterize the genetic relationships within and among four neighboring ethnolinguistic groups in northern Kenya in light of cultural relationships to understand the extent to which geography and culture shape patterns of genetic variation.

MATERIALS AND METHODS

We collected DNA and demographic information pertaining to aspects of social identity and heritage from 572 individuals across the Turkana, Samburu, Waso Borana, and Rendille of northern Kenya. We sampled individuals across a total of nine clans from these four groups and, additionally, three territorial sections within the Turkana and successfully genotyped 376 individuals.

RESULTS

Here we report that geography predominately shapes genetic variation within and among human groups in northern Kenya. We observed a clinal pattern of genetic variation that mirrors the overall geographic distribution of the individuals we sampled. We also found relatively higher rates of intermarriage between the Rendille and Samburu and evidence of gene flow between them that reflect these higher rates of intermarriage. Among the Turkana, we observed strong recent genetic substructuring based on territorial section affiliation. Within ethnolinguistic groups, we found that Y chromosome haplotypes do not consistently cluster by natal clan affiliation. Finally, we found that sampled populations that are geographically closer have lower genetic differentiation, and that cultural similarity does not predict genetic similarity as a whole across these northern Kenyan populations.

DISCUSSION

Overall, the results from this study highlight the importance of geography, even on a local geographic scale, in shaping observed patterns of genetic variation in human populations.

A genealogical estimate of genetic relationships

The application of genetic relationships among individuals, characterized by a genetic relationship matrix (GRM), has far-reaching effects in human genetics. However, the current standard to calculate the GRM treats linked markers as independent and does not explicitly model the underlying genealogical history of the study sample. Here, we propose a coalescent-informed framework, namely the expected GRM (eGRM), to infer the expected relatedness between pairs of individuals given an ancestral recombination graph (ARG) of the sample. Through extensive simulations, we show that the eGRM is an unbiased estimate of latent pairwise genome-wide relatedness and is robust when computed with ARG inferred from incomplete genetic data. As a result, the eGRM better captures the structure of a population than the canonical GRM, even when using the same genetic information. More importantly, our framework allows a principled approach to estimate the eGRM at different time depths of the ARG, thereby revealing the time-varying nature of population structure in a sample. When applied to SNP array genotypes from a population sample from Northern and Eastern Finland, we find that clustering analysis with the eGRM reveals population structure driven by subpopulations that would not be apparent via the canonical GRM and that temporally the population model is consistent with recent divergence and expansion. Taken together, our proposed eGRM provides a robust tree-centric estimate of relatedness with wide application to genetic studies.

Pharmacogenetic distinction of the Croatian population from the European average

AIM

To compare the Croatian and European population in terms of allele frequencies of clinically relevant polymorphisms in drug absorption, distribution, metabolism, and excretion (ADME) genes.

METHODS

In 429 Croatian participants, we genotyped 27 loci in 20 ADME genes. The obtained frequencies were merged with the published frequencies for the Croatian population by sample size weighting. The study sample obtained in this way was compared with the average data for the European population from the gnomAD database.

RESULTS

Variant allele frequencies in the Croatian population were higher in three and lower in two polymorphisms (Benjamini-Hochberg-corrected P values: 0.0027 for CYP2B6*4 rs2279343, CYP2C9*2 rs1799853, and VKORC1 rs9923231; 0.0297 for GSTP1 rs1695; 0.0455 for CYP2A6 rs1801272) compared with the European population. The most marked difference was observed for CYP2B6*4 (9.3% in Europe vs 24.3% in Croatia). The most clinically relevant findings were higher variant allele frequencies in two polymorphisms related to lower warfarin requirements: VKORC1*2 (34.9% in Europe vs 40.1% in Croatia) and CYP2C9*2 (12.3% in Europe vs 14.7% in Croatia). This indicates that three-quarters of Croatian people have at least one variant allele at these loci. Variants in genes GSTP1 and CYP2A6 were significantly less frequently observed in Croatia.

CONCLUSIONS

Croatian population has a higher bleeding and over-anticoagulation risk, which is why we recommend the prescription of lower doses of anticoagulation drugs such as warfarin and acenocoumarol. Lower phenytoin, and higher bupropion and efavirenz doses are also recommended in the Croatian population.

Genome-wide association studies of metabolites in Finnish men identify disease-relevant loci

Few studies have explored the impact of rare variants (minor allele frequency < 1%) on highly heritable plasma metabolites identified in metabolomic screens. The Finnish population provides an ideal opportunity for such explorations, given the multiple bottlenecks and expansions that have shaped its history, and the enrichment for many otherwise rare alleles that has resulted. Here, we report genetic associations for 1391 plasma metabolites in 6136 men from the late-settlement region of Finland. We identify 303 novel association signals, more than one third at variants rare or enriched in Finns. Many of these signals identify genes not previously implicated in metabolite genome-wide association studies and suggest mechanisms for diseases and disease-related traits.

The causal effect of interleukin-17 on the risk of psoriatic arthritis: a Mendelian randomization study

OBJECTIVE

To determine causal associations between genetically predicted TNF-α, IL-12p70 and IL-17 levels and risk of PsA.

METHODS

The publicly available summary-level findings from genome-wide association studies (GWAS) was used to identify loci influencing normal physiological concentrations of TNF-α, IL-12p70 and IL-17 (n = 8293) among healthy individuals as exposure and a GWAS for PsA from the UK Biobank (PsA = 900, control = 462 033) as the outcome. A two-sample Mendelian randomization (MR) analysis was performed using the inverse-variance weighted (IVW), weighted median and MR-Egger regression methods. Sensitivity analysis and MR-Egger regression analysis were performed to evaluate the heterogeneity and pleiotropic effects of each variant.

RESULTS

Single-nucleotide polymorphisms (SNPs) at genome-wide significance from GWASs on TNF-α, IL-12p70 and IL-17 were identified as the instrumental variables. The IVW method indicated a causal association between increased IL-17 level and risk of PsA (β = -0.00186 per allele, s.e. = 0.00043, P = 0.002). Results were consistent in the weighted median method (β = -0.00145 per allele, s.e. = 0.00059, P = 0.014) although the MR-Egger method suggested a non-significant association (β = -0.00133 per allele, s.e. = 0.00087; P = 0.087). Single SNP MR results revealed that the C allele of rs117556572 was robustly associated with risk of PsA (β = 0.00210, s.e. = 0.00069, P = 0.002). However, no evidence for a causal effect was observed between TNF-α, IL-12p70, decreased IL-17 levels and risk of PsA.

CONCLUSION

Our findings provide preliminary evidence that genetic variants predisposing to higher physiological IL-17 level are associated with decreased risk of PsA.

Long runs of homozygosity are associated with Alzheimer's disease

Long runs of homozygosity (ROH) are contiguous stretches of homozygous genotypes, which are a footprint of inbreeding and recessive inheritance. The presence of recessive loci is suggested for Alzheimer's disease (AD); however, their search has been poorly assessed to date. To investigate homozygosity in AD, here we performed a fine-scale ROH analysis using 10 independent cohorts of European ancestry (11,919 AD cases and 9181 controls.) We detected an increase of homozygosity in AD cases compared to controls [βAVROH (CI 95%) = 0.070 (0.037-0.104); P = 3.91 × 10-5; βFROH (CI95%) = 0.043 (0.009-0.076); P = 0.013]. ROHs increasing the risk of AD (OR > 1) were significantly overrepresented compared to ROHs increasing protection (p < 2.20 × 10-16). A significant ROH association with AD risk was detected upstream the HS3ST1 locus (chr4:11,189,482‒11,305,456), (β (CI 95%) = 1.09 (0.48 ‒ 1.48), p value = 9.03 × 10-4), previously related to AD. Next, to search for recessive candidate variants in ROHs, we constructed a homozygosity map of inbred AD cases extracted from an outbred population and explored ROH regions in whole-exome sequencing data (N = 1449). We detected a candidate marker, rs117458494, mapped in the SPON1 locus, which has been previously associated with amyloid metabolism. Here, we provide a research framework to look for recessive variants in AD using outbred populations. Our results showed that AD cases have enriched homozygosity, suggesting that recessive effects may explain a proportion of AD heritability.

Recent Consanguinity and Outbred Autozygosity Are Associated With Increased Risk of Late-Onset Alzheimer's Disease

Prior work in late-onset Alzheimer's disease (LOAD) has resulted in discrepant findings as to whether recent consanguinity and outbred autozygosity are associated with LOAD risk. In the current study, we tested the association between consanguinity and outbred autozygosity with LOAD in the largest such analysis to date, in which 20 LOAD GWAS datasets were retrieved through public databases. Our analyses were restricted to eight distinct ethnic groups: African-Caribbean, Ashkenazi-Jewish European, European-Caribbean, French-Canadian, Finnish European, North-Western European, South-Eastern European, and Yoruba African for a total of 21,492 unrelated subjects (11,196 LOAD and 10,296 controls). Recent consanguinity determination was performed using FSuite v1.0.3, according to subjects' ancestral background. The level of autozygosity in the outbred population was assessed by calculating inbreeding estimates based on the proportion (FROH) and the number (NROH) of runs of homozygosity (ROHs). We analyzed all eight ethnic groups using a fixed-effect meta-analysis, which showed a significant association of recent consanguinity with LOAD (N = 21,481; OR = 1.262, P = 3.6 × 10-4), independently of APOE ∗4 (N = 21,468, OR = 1.237, P = 0.002), and years of education (N = 9,257; OR = 1.274, P = 0.020). Autozygosity in the outbred population was also associated with an increased risk of LOAD, both for F ROH (N = 20,237; OR = 1.204, P = 0.030) and N ROH metrics (N = 20,237; OR = 1.019, P = 0.006), independently of APOE ∗4 [(F ROH, N = 20,225; OR = 1.222, P = 0.029) (N ROH, N = 20,225; OR = 1.019, P = 0.007)]. By leveraging the Alzheimer's Disease Sequencing Project (ADSP) whole-exome sequencing (WES) data, we determined that LOAD subjects do not show an enrichment of rare, risk-enhancing minor homozygote variants compared to the control population. A two-stage recessive GWAS using ADSP data from 201 consanguineous subjects in the discovery phase followed by validation in 10,469 subjects led to the identification of RPH3AL p.A303V (rs117190076) as a rare minor homozygote variant increasing the risk of LOAD [discovery: Genotype Relative Risk (GRR) = 46, P = 2.16 × 10-6; validation: GRR = 1.9, P = 8.0 × 10-4]. These results confirm that recent consanguinity and autozygosity in the outbred population increase risk for LOAD. Subsequent work, with increased samples sizes of consanguineous subjects, should accelerate the discovery of non-additive genetic effects in LOAD.

Benefits and limitations of genome-wide association studies

Genome-wide association studies (GWAS) involve testing genetic variants across the genomes of many individuals to identify genotype-phenotype associations. GWAS have revolutionized the field of complex disease genetics over the past decade, providing numerous compelling associations for human complex traits and diseases. Despite clear successes in identifying novel disease susceptibility genes and biological pathways and in translating these findings into clinical care, GWAS have not been without controversy. Prominent criticisms include concerns that GWAS will eventually implicate the entire genome in disease predisposition and that most association signals reflect variants and genes with no direct biological relevance to disease. In this Review, we comprehensively assess the benefits and limitations of GWAS in human populations and discuss the relevance of performing more GWAS.

Disparities in discovery of pathogenic variants for autosomal recessive non-syndromic hearing impairment by ancestry

Hearing impairment (HI) is characterized by extensive genetic heterogeneity. To determine the population-specific contribution of known autosomal recessive nonsyndromic (ARNS)HI genes and variants to HI etiology; pathogenic and likely pathogenic (PLP) ARNSHI variants were selected from ClinVar and the Deafness Variation Database and their frequencies were obtained from gnomAD for seven populations. ARNSHI prevalence due to PLP variants varies greatly by population ranging from 96.9 affected per 100,000 individuals for Ashkenazi Jews to 5.2 affected per 100,000 individuals for Africans/African Americans. For Europeans, Finns have the lowest prevalence due to ARNSHI PLP variants with 9.5 affected per 100,000 individuals. For East Asians, Latinos, non-Finish Europeans, and South Asians, ARNSHI prevalence due to PLP variants ranges from 17.1 to 33.7 affected per 100,000 individuals. ARNSHI variants that were previously reported in a single ancestry or family were observed in additional populations, e.g., USH1C p.(Q723*) reported in a Chinese family was the most prevalent pathogenic variant observed in gnomAD for African/African Americans. Variability between populations is due to how extensively ARNSHI has been studied, ARNSHI prevalence and ancestry specific ARNSHI variant architecture which is impacted by population history. Our study demonstrates that additional gene and variant discovery studies are necessary for all populations and particularly for individuals of African ancestry.

Exome sequencing of Finnish isolates enhances rare-variant association power

Exome-sequencing studies have generally been underpowered to identify deleterious alleles with a large effect on complex traits as such alleles are mostly rare. Because the population of northern and eastern Finland has expanded considerably and in isolation following a series of bottlenecks, individuals of these populations have numerous deleterious alleles at a relatively high frequency. Here, using exome sequencing of nearly 20,000 individuals from these regions, we investigate the role of rare coding variants in clinically relevant quantitative cardiometabolic traits. Exome-wide association studies for 64 quantitative traits identified 26 newly associated deleterious alleles. Of these 26 alleles, 19 are either unique to or more than 20 times more frequent in Finnish individuals than in other Europeans and show geographical clustering comparable to Mendelian disease mutations that are characteristic of the Finnish population. We estimate that sequencing studies of populations without this unique history would require hundreds of thousands to millions of participants to achieve comparable association power.

Geographic Variation and Bias in the Polygenic Scores of Complex Diseases and Traits in Finland

Polygenic scores (PSs) are becoming a useful tool to identify individuals with high genetic risk for complex diseases, and several projects are currently testing their utility for translational applications. It is also tempting to use PSs to assess whether genetic variation can explain a part of the geographic distribution of a phenotype. However, it is not well known how the population genetic properties of the training and target samples affect the geographic distribution of PSs. Here, we evaluate geographic differences, and related biases, of PSs in Finland in a geographically well-defined sample of 2,376 individuals from the National FINRISK study. First, we detect geographic differences in PSs for coronary artery disease (CAD), rheumatoid arthritis, schizophrenia, waist-hip ratio (WHR), body-mass index (BMI), and height, but not for Crohn disease or ulcerative colitis. Second, we use height as a model trait to thoroughly assess the possible population genetic biases in PSs and apply similar approaches to the other phenotypes. Most importantly, we detect suspiciously large accumulations of geographic differences for CAD, WHR, BMI, and height, suggesting bias arising from the population's genetic structure rather than from a direct genotype-phenotype association. This work demonstrates how sensitive the geographic patterns of current PSs are for small biases even within relatively homogeneous populations and provides simple tools to identify such biases. A thorough understanding of the effects of population genetic structure on PSs is essential for translational applications of PSs.

Heterogeneity in the extent of linkage disequilibrium among exonic, intronic, non-coding RNA and intergenic chromosome regions

Whole-genome sequence data enable construction of high-resolution linkage disequilibrium (LD) maps revealing the LD structure of functional elements within genic and subgenic sequences. The Malecot-Morton model defines LD map distances in linkage disequilibrium units (LDUs), analogous to the centimorgan scale of linkage maps. For whole-genome sequence-derived LD maps, we introduce the ratio of corresponding map lengths kilobases/LDU to describe the extent of LD within genome components. The extent of LD is highly variable across the genome ranging from ~38 kb for intergenic sequences to ~858 kb for centromeric regions. LD is ~16% more extensive in genic, compared with intergenic sequences, reflecting relatively increased selection and/or reduced recombination in genes. The LD profile across 18,268 autosomal genes reveals reduced extent of LD, consistent with elevated recombination, in exonic regions near the 5' end of genes but more extensive LD, compared with intronic sequences, across more centrally located exons. Genes classified as essential and genes linked to Mendelian phenotypes show more extensive LD compared with genes associated with complex traits, perhaps reflecting differences in selective pressure. Significant differences between exonic, intronic and intergenic components demonstrate that fine-scale LD structure provides important insights into genome function, which cannot be revealed by LD analysis of much lower resolution array-based genotyping and conventional linkage maps.

A genome-wide association study of tramadol metabolism from post-mortem samples

Phase I tramadol metabolism requires cytochrome p450 family 2, subfamily D, polypeptide 6 (CYP2D6) to form O-desmethyltramadol (M1). CYP2D6 genetic variants may infer metabolizer phenotype; however, drug ADME (absorption, distribution, metabolism, and excretion) and response depend on protein pathway(s), not CYP2D6 alone. There is a paucity of data regarding the contribution of trans-acting proteins to idiosyncratic phenotypes following drug exposure. A genome-wide association study identified five markers (rs79983226/kgp11274252, rs9384825, rs62435418/kgp10370907, rs72732317/kgp3743668, and rs184199168/exm1592932) associated with the conversion of tramadol to M1 (M1:T). These SNPs reside within five genes previously implicated with adverse reactions. Analysis of accompanying toxicological meta-data revealed a significant positive linear relationship between M1:T and degree of sample polypharmacy. Taken together, these data identify candidate loci for potential clinical inferences of phenotype following exposure to tramadol and highlight sample polypharmacy as a possible diagnostic covariate in post-mortem genetic studies.

Contribution of rare and common variants to intellectual disability in a sub-isolate of Northern Finland

The contribution of de novo variants in severe intellectual disability (ID) has been extensively studied whereas the genetics of mild ID has been less characterized. To elucidate the genetics of milder ID we studied 442 ID patients enriched for mild ID (>50%) from a population isolate of Finland. Using exome sequencing, we show that rare damaging variants in known ID genes are observed significantly more often in severe (27%) than in mild ID (13%) patients. We further observe a significant enrichment of functional variants in genes not yet associated with ID (OR: 2.1). We show that a common variant polygenic risk significantly contributes to ID. The heritability explained by polygenic risk score is the highest for educational attainment (EDU) in mild ID (2.2%) but lower for more severe ID (0.6%). Finally, we identify a Finland enriched homozygote variant in the CRADD ID associated gene.

Genetic investigation into an increased susceptibility to biliary atresia in an extended New Zealand Māori family

Background

Biliary atresia (BA), a fibrosing disorder of the developing biliary tract leading to liver failure in infancy, has an elevated incidence in indigenous New Zealand (NZ) Māori. We investigated a high rate of BA in a group of children (n = 12) belonging to a single Māori iwi (or ‘tribe’, related through a remote ancestor).

Methods

Population and geographical data was used to estimate the rate of BA in Māori sub-groups, and a pedigree linking most of the affected children was constructed from oral and documented history. Array genotyping was used to examine hypotheses about the inheritance of a possible genetic risk factor, and the history of the affected population, and Exome Sequencing to search for candidate genes.

Results

Most of these affected children (n = 7) link to a self-reported pedigree and carry a 50-fold increase in BA risk over unrelated Māori (χ² = 296P < 0.001, 95% CI 23–111). Genetic analysis using FEstim and SNP array genotypes revealed no evidence for elevated consanguinity between parents of affected children (FEstim: F (2,21) = 0.469, P > 0.63). Genome-wide quantitation of intervals of contiguous, homozygous-by-state markers reached a similar conclusion (F (2,399) = 1.99, P = 0.138). Principal component analysis and investigation with STRUCTURE found no evidence of increased allele frequency of either a recessive variant, or additive, low-risk variants due to reproductive isolation. To identify candidate causal factors, Exome Sequencing datasets were scrutinised for shared rare coding variants across 8 affected individuals. No rare, non-synonymous, phylogenetically conserved variants were common to 6 or more affected children.

Conclusion

The substantially elevated risk for development of BA in this subgroup could be mediated by genetic factors, but the iwi exhibits no properties indicative of recent or remote reproductive isolation. Resolution of any risk loci may rely on extensive genomic sequencing studies in this iwi or investigation of other mechnaisms such as copy number variation.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0440-0) contains supplementary material, which is available to authorized users.

Frequencies of poor metabolizer alleles of 12 pharmacogenomic actionable genes in Punjabi Sikhs of Indian Origin

Diversity in drug response is attributed to both genetic and non-genetic factors. However, there is paucity of pharmacogenetics information across ethnically and genetically diverse populations of India. Here, we have analyzed 21 SNPs from 12 pharmacogenomics genes in Punjabi Sikhs of Indian origin (N = 1,616), as part of the Sikh Diabetes Study (SDS). We compared the allele frequency of poor metabolism (PM) phenotype among Sikhs across other major global populations from the Exome Aggregation Consortium and 1000 Genomes. The PM phenotype of CYP1A2*1 F for slow metabolism of caffeine and carcinogens was significantly higher in Indians (SDS 42%, GIH [Gujarati] 51%, SAS [Pakistani] 45%) compared to Europeans 29% (pgenotype = 5.3E-05). Similarly, South Asians had a significantly higher frequency of CYP2C9*3 (12% SDS, 13% GIH, 11% SAS) vs. 7% in Europeans (pgenotype = <1.0E-05) and 'T' allele of CYP4F2 (36%) SDS, (43%) GIH, 40% (SAS) vs. (29%) in Europeans (pgenotype = <1.0E-05); both associated with a higher risk of bleeding with warfarin. All South Asians -the Sikhs (0.36), GIH (0.34), and SAS (0.36) had a higher frequency of the NAT2*6 allele (linked with slow acetylation of isoniazid) compared to Europeans (0.29). Additionally, the prevalence of the low activity 'C' allele of MTHFR (rs1801131) was highest in Sikhs compared to all other ethnic groups [SDS (44%), GIH (39%), SAS (42%) and European (32%) (pgenotype = <1.0E-05)]. SNPs in MTHFR affect metabolism of statins, 5-fluorouracil and methotrexate-based cancer drugs. These findings underscore the need for evaluation of other endogamous ethnic groups of India and beyond for establishing a global benchmark for pre-emptive genotyping in drug metabolizing genes before beginning therapeutic intervention.

Haplotype Sharing Provides Insights into Fine-Scale Population History and Disease in Finland

Finland provides unique opportunities to investigate population and medical genomics because of its adoption of unified national electronic health records, detailed historical and birth records, and serial population bottlenecks. We assembled a comprehensive view of recent population history (≤100 generations), the timespan during which most rare-disease-causing alleles arose, by comparing pairwise haplotype sharing from 43,254 Finns to that of 16,060 Swedes, Estonians, Russians, and Hungarians from geographically and linguistically adjacent countries with different population histories. We find much more extensive sharing in Finns, with at least one ≥ 5 cM tract on average between pairs of unrelated individuals. By coupling haplotype sharing with fine-scale birth records from more than 25,000 individuals, we find that although haplotype sharing broadly decays with geographical distance, there are pockets of excess haplotype sharing; individuals from northeast Finland typically share several-fold more of their genome in identity-by-descent segments than individuals from southwest regions. We estimate recent effective population-size changes through time across regions of Finland, and we find that there was more continuous gene flow as Finns migrated from southwest to northeast between the early- and late-settlement regions than was dichotomously described previously. Lastly, we show that haplotype sharing is locally enriched by an order of magnitude among pairs of individuals sharing rare alleles and especially among pairs sharing rare disease-causing variants. Our work provides a general framework for using haplotype sharing to reconstruct an integrative view of recent population history and gain insight into the evolutionary origins of rare variants contributing to disease.

The genetic structure of the Belgian population

BACKGROUND

National and international efforts like the 1000 Genomes Project are leading to increasing insights in the genetic structure of populations worldwide. Variation between different populations necessitates access to population-based genetic reference datasets. These data, which are important not only in clinical settings but also to potentiate future transitions towards a more personalized public health approach, are currently not available for the Belgian population.

RESULTS

To obtain a representative genetic dataset of the Belgian population, participants in the 2013 National Health Interview Survey (NHIS) were invited to donate saliva samples for DNA analysis. DNA was isolated and single nucleotide polymorphisms (SNPs) were determined using a genome-wide SNP array of around 300,000 sites, resulting in a high-quality dataset of 189 samples that was used for further analysis. A principal component analysis demonstrated the typical European genetic constitution of the Belgian population, as compared to other continents. Within Europe, the Belgian population could be clearly distinguished from other European populations. Furthermore, obvious signs from recent migration were found, mainly from Southern Europe and Africa, corresponding with migration trends from the past decades. Within Belgium, a small north-west to south-east gradient in genetic variability was noted, with differences between Flanders and Wallonia.

CONCLUSIONS

This is the first study on the genetic structure of the Belgian population and its regional variation. The Belgian genetic structure mirrors its geographic location in Europe with regional differences and clear signs of recent migration.

Gene-panel testing of breast and ovarian cancer patients identifies a recurrent RAD51C duplication

Gene-panel sequencing allows comprehensive analysis of multiple genes simultaneously and is now routinely used in clinical mutation testing of high-risk breast and ovarian cancer patients. However, only BRCA1 and BRCA2 are often analyzed also for large genomic changes. Here, we have analyzed 10 clinically relevant susceptibility genes in 95 breast or ovarian cancer patients with gene-panel sequencing including also copy number variants (CNV) analysis for genomic changes. We identified 12 different pathogenic BRCA1, BRCA2, TP53, PTEN, CHEK2, or RAD51C mutations in 18 of 95 patients (19%). BRCA1/2 mutations were observed in 8 patients (8.4%) and CHEK2 protein-truncating mutations in 7 patients (7.4%). In addition, we identified a novel duplication encompassing most of the RAD51C gene. We further genotyped the duplication in breast or ovarian cancer families (n = 1149), in unselected breast (n = 1729) and ovarian cancer cohorts (n = 553), and in population controls (n = 1273). Seven additional duplication carries were observed among cases but none among controls. The duplication associated with ovarian cancer risk (3/590 of all ovarian cancer patients, 0.5%, P = .032 compared with controls) and was found to represent a large fraction of all identified RAD51C mutations in the Finnish population. Our data emphasizes the importance of comprehensive mutation analysis including CNV detection in all the relevant genes.

Weighted likelihood inference of genomic autozygosity patterns in dense genotype data

Background

Genomic regions of autozygosity (ROA) arise when an individual is homozygous for haplotypes inherited identical-by-descent from ancestors shared by both parents. Over the past decade, they have gained importance for understanding evolutionary history and the genetic basis of complex diseases and traits. However, methods to infer ROA in dense genotype data have not evolved in step with advances in genome technology that now enable us to rapidly create large high-resolution genotype datasets, limiting our ability to investigate their constituent ROA patterns.

Methods

We report a weighted likelihood approach for inferring ROA in dense genotype data that accounts for autocorrelation among genotyped positions and the possibilities of unobserved mutation and recombination events, and variability in the confidence of individual genotype calls in whole genome sequence (WGS) data.

Results

Forward-time genetic simulations under two demographic scenarios that reflect situations where inbreeding and its effect on fitness are of interest suggest this approach is better powered than existing state-of-the-art methods to infer ROA at marker densities consistent with WGS and popular microarray genotyping platforms used in human and non-human studies. Moreover, we present evidence that suggests this approach is able to distinguish ROA arising via consanguinity from ROA arising via endogamy. Using subsets of The 1000 Genomes Project Phase 3 data we show that, relative to WGS, intermediate and long ROA are captured robustly with popular microarray platforms, while detection of short ROA is more variable and improves with marker density. Worldwide ROA patterns inferred from WGS data are found to accord well with those previously reported on the basis of microarray genotype data. Finally, we highlight the potential of this approach to detect genomic regions enriched for autozygosity signals in one group relative to another based upon comparisons of per-individual autozygosity likelihoods instead of inferred ROA frequencies.

Conclusions

This weighted likelihood ROA inference approach can assist population- and disease-geneticists working with a wide variety of data types and species to explore ROA patterns and to identify genomic regions with differential ROA signals among groups, thereby advancing our understanding of evolutionary history and the role of recessive variation in phenotypic variation and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4312-3) contains supplementary material, which is available to authorized users.

Fine-Scale Genetic Structure in Finland

Coupling dense genotype data with new computational methods offers unprecedented opportunities for individual-level ancestry estimation once geographically precisely defined reference data sets become available. We study such a reference data set for Finland containing 2376 such individuals from the FINRISK Study survey of 1997 both of whose parents were born close to each other. This sampling strategy focuses on the population structure present in Finland before the 1950s. By using the recent haplotype-based methods ChromoPainter (CP) and FineSTRUCTURE (FS) we reveal a highly geographically clustered genetic structure in Finland and report its connections to the settlement history as well as to the current dialectal regions of the Finnish language. The main genetic division within Finland shows striking concordance with the 1323 borderline of the treaty of Nöteborg. In general, we detect genetic substructure throughout the country, which reflects stronger regional genetic differences in Finland compared to, for example, the UK, which in a similar analysis was dominated by a single unstructured population. We expect that similar population genetic reference data sets will become available for many more populations in the near future with important applications, for example, in forensic genetics and in genetic association studies. With this in mind, we report those extensions of the CP + FS approach that we found most useful in our analyses of the Finnish data.

Genetics in an isolated population like Finland: a different basis for genomic medicine?

A unique genetic background in an isolated population like that of Finland offers special opportunities for genetic research as well as for applying the genetic developments to the health care. On the other hand, the different genetic background may require local attempts to develop diagnostics and treatment as the selection of diseases and mutations differs from that in the other populations. In this review, we describe the experiences of research and health care in this genetic isolate starting from the identification of specific monogenic diseases enriched in the Finnish population all the way to implementing the knowledge of the unique genetic background to genomic medicine at population level.

Candidate susceptibility variants for esophageal squamous cell carcinoma

Esophageal cancer is common worldwide, and often fatal. The major histological subtype is esophageal squamous cell carcinoma (ESCC). ESCC shows familial aggregation and high heritability. Mutations in RHBDF2 cause tylosis, a very rare disorder characterized by high life-time risk of ESCC, but no other well-established predisposition genes have been identified. To identify candidate susceptibility variants for ESCC we utilized the Population Information System and the Finnish cancer registry to find study materials by clustering ESCC patients by family name at birth and municipality at birth. We collected archival tissue material and exome sequenced a total of 30 ESCC cases. We prioritized shared, deleterious and rare variants that were significantly enriched in our sample set compared to Finnish and population subset specific controls. Six variants passed filtering, the most frequent being a nonsense mutation in DNAH9 (p.Tyr1573Ter) found in four unrelated patients. DNAH9 has been reported to be frequently lost in ESCC tumors. In this study, one patient's tumor showed loss of the wild type allele of DNAH9 suggesting a tumor suppressive function. A missense variant in GKAP1 was shared by three patients, and missense variants in BAG1, NFX1, FUK, and DDOST by two each. EP300 which has previously been implicated in the genesis of ESCC had a missense variant segregating in three affected individuals in a single family. If validated in independent patient sets, these variants could serve as a tool towards prevention and early diagnosis of ESCC.

Whole-genome view of the consequences of a population bottleneck using 2926 genome sequences from Finland and United Kingdom

Isolated populations with enrichment of variants due to recent population bottlenecks provide a powerful resource for identifying disease-associated genetic variants and genes. As a model of an isolate population, we sequenced the genomes of 1463 Finnish individuals as part of the Sequencing Initiative Suomi (SISu) Project. We compared the genomic profiles of the 1463 Finns to a sample of 1463 British individuals that were sequenced in parallel as part of the UK10K Project. Whereas there were no major differences in the allele frequency of common variants, a significant depletion of variants in the rare frequency spectrum was observed in Finns when comparing the two populations. On the other hand, we observed >2.1 million variants that were twice as frequent among Finns compared with Britons and 800 000 variants that were more than 10 times more frequent in Finns. Furthermore, in Finns we observed a relative proportional enrichment of variants in the minor allele frequency range between 2 and 5% (P<2.2 × 10⁻¹⁶). When stratified by their functional annotations, loss-of-function variants showed the highest proportional enrichment in Finns (P=0.0291). In the non-coding part of the genome, variants in conserved regions (P=0.002) and promoters (P=0.01) were also significantly enriched in the Finnish samples. These functional categories represent the highest a priori power for downstream association studies of rare variants using population isolates.

The effects of resource availability and the demographic transition on the genetic correlation between number of children and grandchildren in humans

Studies of evolutionary change require an estimate of fitness, and lifetime reproductive success is widely used for this purpose. However, many species face a trade-off between the number and quality of offspring and in such cases number of grandoffspring may better represent the genetic contribution to future generations. Here, we apply quantitative genetic methods to a genealogical data set on humans from Finland to address how the genetic correlation between number of children and grandchildren is influenced by the severity of the trade-off between offspring quality and quantity, as estimated by different levels of resource access among individuals in the population. Further, we compare the genetic correlation before and after the demographic transition to low mortality and fertility rates. The genetic correlation was consistently high (0.79-0.92) with the strongest correlations occurring in individuals with higher access to resources and before the demographic transition, and a tendency for lower correlations in resource poor individuals and after the transition. These results indicate that number of grandoffspring is a slightly better predictor of long-term genetic fitness than number of offspring in a human population across a range of environmental conditions, and more generally, that patterns of resource availability need to be taken into account when estimating genetic covariances with fitness.

Nationwide Genomic Study in Denmark Reveals Remarkable Population Homogeneity

Denmark has played a substantial role in the history of Northern Europe. Through a nationwide scientific outreach initiative, we collected genetic and anthropometrical data from ∼800 high school students and used them to elucidate the genetic makeup of the Danish population, as well as to assess polygenic predictions of phenotypic traits in adolescents. We observed remarkable homogeneity across different geographic regions, although we could still detect weak signals of genetic structure reflecting the history of the country. Denmark presented genomic affinity with primarily neighboring countries with overall resemblance of decreasing weight from Britain, Sweden, Norway, Germany, and France. A Polish admixture signal was detected in Zealand and Funen, and our date estimates coincided with historical evidence of Wend settlements in the south of Denmark. We also observed considerably diverse demographic histories among Scandinavian countries, with Denmark having the smallest current effective population size compared to Norway and Sweden. Finally, we found that polygenic prediction of self-reported adolescent height in the population was remarkably accurate (R² = 0.639 ± 0.015). The high homogeneity of the Danish population could render population structure a lesser concern for the upcoming large-scale gene-mapping studies in the country.

Exome sequencing in pooled DNA samples to identify maternal pre-eclampsia risk variants

Pre-eclampsia is a common pregnancy disorder that is a major cause for maternal and perinatal mortality and morbidity. Variants predisposing to pre-eclampsia might be under negative evolutionary selection that is likely to keep their population frequencies low. We exome sequenced samples from a hundred Finnish pre-eclamptic women in pools of ten to screen for low-frequency, large-effect risk variants for pre-eclampsia. After filtering and additional genotyping steps, we selected 28 low-frequency missense, nonsense and splice site variants that were enriched in the pre-eclampsia pools compared to reference data, and genotyped the variants in 1353 pre-eclamptic and 699 non-pre-eclamptic women to test the association of them with pre-eclampsia and quantitative traits relevant for the disease. Genotypes from the SISu project (n = 6118 exome sequenced Finnish samples) were included in the binary trait association analysis as a population reference to increase statistical power. In these analyses, none of the variants tested reached genome-wide significance. In conclusion, the genetic risk for pre-eclampsia is likely complex even in a population isolate like Finland, and larger sample sizes will be necessary to detect risk variants.

Genetic structure of the Newfoundland and Labrador population: founder effects modulate variability

The population of the province of Newfoundland and Labrador (NL) has been a resource for genetic studies because of its historical isolation and increased prevalence of several monogenic disorders. Controversy remains regarding the genetic substructure and the extent of genetic homogeneity, which have implications for disease gene mapping. Population substructure has been reported from other isolated populations such as Iceland, Finland and Sardinia. We undertook this study to further our understanding of the genetic architecture of the NL population. We enrolled 494 individuals randomly selected from NL. Genome-wide SNP data were analyzed together with that from 14 other populations including HapMap3, Ireland, Britain and Native American samples from the Human Genome Diversity Project. Using multidimensional scaling and admixture analysis, we observed that the genetic structure of the NL population resembles that of the British population but can be divided into three clusters that correspond to religious/ethnic origins: Protestant English, Roman Catholic Irish and North American aboriginals. We observed reduced heterozygosity and an increased inbreeding coefficient (mean=0.005), which corresponds to that expected in the offspring of third-cousin marriages. We also found that the NL population has a significantly higher number of runs of homozygosity (ROH) and longer lengths of ROH segments. These results are consistent with our understanding of the population history and indicate that the NL population may be ideal for identifying recessive variants for complex diseases that affect populations of European origin.

Estimation of Recent and Ancient Inbreeding in a Small Endogamous Tunisian Community Through Genomic Runs of Homozygosity

Runs of homozygosity (ROHs) are extended genomic regions of homozygous genotypes that record populations' mating patterns in the past. We performed microarray genotyping on 15 individuals from a small isolated Tunisian community. We estimated the individual and population genome-wide level of homozygosity from data on ROH above 0.5 Mb in length. We found a high average number of ROH per individual (48.2). The smallest ROH category (0.5-1.49 Mb) represents 0.93% of the whole genome, while medium-size (1.5-4.99 Mb) and long-size ROH (≥5 Mb) cover 1.18% and 0.95%, respectively. We found that genealogical individual inbreeding coefficients (Fped ) based on three- to four-generation pedigrees are not reliable indicators of the current proportion of genome-wide homozygosity inferred from ROH (FROH ) either for 0.5 or 1.5 Mb ROH length thresholds, while identity-by-descent sharing is a function of shared coancestry. This study emphasizes the effect of reproductive isolation and a prolonged practice of consanguinity that limits the genetic heterogeneity. It also provides evidence of both recent and ancient parental relatedness contribution to the current level of genome-wide homozygosity in the studied population. These findings may be useful for evaluation of long-term effects of inbreeding on human health and for future applications of ROHs in identifying recessive susceptibility genes.

CXCR3 Polymorphism and Expression Associate with Spontaneous Preterm Birth

Spontaneous preterm birth (SPTB) is a major factor associating with deaths and with lowered quality of life in humans. Environmental and genetic factors influence the susceptibility. Previously, by analyzing families with recurrent SPTB in linkage analysis, we identified a linkage peak close to the gene encoding CXCR3. Present objectives were to investigate the association of CXCR3 with SPTB in Finnish mothers (n = 443) and infants (n = 747), to analyze CXCR3 expression levels in human placenta and levels of its ligands in umbilical cord blood, and to verify the influence of Cxcr3 on SPTB-associating cytokines in mice. We detected an association between an intronic CXCR3 polymorphism, rs2280964, and SPTB in infants from families with recurrent preterm births (p = 0.009 versus term controls, odds ratio 0.52, 95% confidence interval 0.32-0.86). The minor allele was protective and undertransmitted to SPTB infants (p = 0.007). In the placenta and fetal membranes, the rs2280964 major allele homozygotes had higher expression levels than minor allele homozygotes; decidual trophoblasts showed strong CXCR3 immunoreactivity. Expression was higher in SPTB placentas compared with those from elective deliveries. Concentration of a CXCR3 ligand, CXCL9, was increased in cord blood from SPTB, and the protective rs2280964 allele was associated with low CXCL9. In CXCR3-deficient mice (Mus musculus), SPTB-associating cytokines were not acutely increased in amniotic fluid after preterm birth-inducing dose of maternal LPS. Our results indicate that CXCR3 contributes to SPTB. Activation of CXCR3 signaling may disturb the maternal-fetal tolerance, and this may promote labor.

Effects of Self-Rated Health and Self-Rated Economic Situation on Depressed Mood Via Life Satisfaction Among Older Adults in Costa Rica

Objective: The study examined the relationship of self-rated health and self-rated economic situation with depressed mood, and life satisfaction as mediator of this relationship among older adults in Costa Rica. Method: A longitudinal study was conducted with a subsample (N = 1,618) from the Costa Rican Longevity and Healthy Aging Study (CRELES). Self-rated health, self-rated economic situation, depressed mood, and life satisfaction were measured at baseline, and depressed mood was reassessed 18 months later. Putative mechanisms for changes in depressed mood were examined by means of conditional process analysis. Results: Self-rated health was negatively associated to depressed mood. This effect took place via life satisfaction. An interaction showed that better economic situation compensated the effect of a low self-rated health on life satisfaction. Discussion: This study suggests that subjective variables such as self-rated health, economic situation, and life satisfaction should be considered when addressing the onset of depressed mood.

Targeted resequencing of the pericentromere of chromosome 2 linked to constitutional delay of growth and puberty

Constitutional delay of growth and puberty (CDGP) is the most common cause of pubertal delay. CDGP is defined as the proportion of the normal population who experience pubertal onset at least 2 SD later than the population mean, representing 2.3% of all adolescents. While adolescents with CDGP spontaneously enter puberty, they are at risk for short stature, decreased bone mineral density, and psychosocial problems. Genetic factors contribute heavily to the timing of puberty, but the vast majority of CDGP cases remain biologically unexplained, and there is no definitive test to distinguish CDGP from pathological absence of puberty during adolescence. Recently, we published a study identifying significant linkage between a locus at the pericentromeric region of chromosome 2 (chr 2) and CDGP in Finnish families. To investigate this region for causal variation, we sequenced chr 2 between the genomic coordinates of 79-124 Mb (genome build GRCh37) in the proband and affected parent of the 13 families contributing most to this linkage signal. One gene, DNAH6, harbored 6 protein-altering low-frequency variants (< 6% in the Finnish population) in 10 of the CDGP probands. We sequenced an additional 135 unrelated Finnish CDGP subjects and utilized the unique Sequencing Initiative Suomi (SISu) population reference exome set to show that while 5 of these variants were present in the CDGP set, they were also present in the Finnish population at similar frequencies. Additional variants in the targeted region could not be prioritized for follow-up, possibly due to gaps in sequencing coverage or lack of functional knowledge of non-genic genomic regions. Thus, despite having a well-characterized sample collection from a genetically homogeneous population with a large population-based reference sequence dataset, we were unable to pinpoint variation in the linked region predisposing delayed puberty. This study highlights the difficulties of detecting genetic variants under linkage regions for complex traits and suggests that advancements in annotation of gene function and regulatory regions of the genome will be critical for solving the genetic background of complex phenotypes like CDGP.

Genome-wide insights into the genetic history of human populations

Although mtDNA and the non-recombining Y chromosome (NRY) studies continue to provide valuable insights into the genetic history of human populations, recent technical, methodological and computational advances and the increasing availability of large-scale, genome-wide data from contemporary human populations around the world promise to reveal new aspects, resolve finer points, and provide a more detailed look at our past demographic history. Genome-wide data are particularly useful for inferring migrations, admixture, and fine structure, as well as for estimating population divergence and admixture times and fluctuations in effective population sizes. In this review, we highlight some of the stories that have emerged from the analyses of genome-wide SNP genotyping data concerning the human history of Southern Africa, India, Oceania, Island South East Asia, Europe and the Americas and comment on possible future study directions. We also discuss advantages and drawbacks of using SNP-arrays, with a particular focus on the ascertainment bias, and ways to circumvent it.

Reference values of spirometry for Finnish adults

BACKGROUND

Diagnostic assessment of lung function necessitates up-to-date reference values. The aim of this study was to estimate reference values for spirometry for the Finnish population between 18 and 80 years and to compare them with the existing Finnish, European and the recently published global GLI2012 reference values.

METHODS

Spirometry was performed for 1380 adults in the population-based FinEsS studies and for 662 healthy non-smoking volunteer adults. Detailed predefined questionnaire screening of diseases and symptoms, and quality control of spirometry yielded a sample of 1000 native Finns (387 men) healthy non-smokers aged 18-83 years. Sex-specific reference values, which are estimated using the GAMLSS method and adjusted for age and height, are provided.

RESULTS

The predicted values for lung volumes are larger than those obtained by GLI2012 prediction for the Caucasian subgroup for forced vital capacity (FVC) by an average 6·2% and 5·1% and forced expiratory volume in 1 s (FEV1) by an average 4·2% and 3·0% in men and women, respectively. GLI2012 slightly overestimated the ratio FEV1/FVC with an age-dependent trend. Most reference equations from other European countries, with the exception of the Swiss SAPALDIA study, showed an underestimation of FVC and FEV1 to varying degrees, and a slight overestimation of FEV1/FVC.

CONCLUSION

This study offers up-to-date reference values of spirometry for native Finns with a wide age range. The GLI2012 predictions seem not to be suitable for clinical use for native Finns due to underestimation of lung volumes.

Polymorphisms of the gene encoding Kit ligand are associated with bronchopulmonary dysplasia

Bronchopulmonary dysplasia (BPD) is a chronic inflammatory lung disease that affects infants born preterm. Family studies indicate that BPD has a significant genetic component.

RATIONALE

We assessed the gene encoding Kit ligand (KITLG) as a candidate for genetic predisposition to moderate-to-severe BPD (controls were infants with no or mild BPD).

STUDY DESIGN

Eight KITLG-tagging single nucleotide polymorphisms (SNPs) were analyzed in cohorts of very preterm infants originating from northern Finland (56 cases and 197 controls), southern Finland (n = 59 + 52), and Canada (n = 58 + 68). Additional replication populations included infants born in Finland (n = 41 + 241) and Hungary (n = 29 + 40). All infants were of European origin. Results were controlled for risk factors of BPD. Kit ligand concentration in umbilical cord blood, collected from very preterm infants (n = 120), was studied.

RESULTS

Six SNPs of KITLG and a haplotype including all eight genotyped SNPs were associated with moderate-to-severe BPD in the northern Finnish population. When all the populations were combined, SNP rs11104948 was significantly associated with BPD. Kit ligand concentration in umbilical cord blood of infants born very preterm was an independent risk factor of BPD.

CONCLUSIONS

We show that KITLG polymorphisms are associated with susceptibility to moderate-to-severe BPD. In addition, higher Kit ligand concentrations were observed in infants that subsequently developed BPD. These results support the possibility that KITLG gene is involved in predisposition to BPD. Pediatr Pulmonol. 2015; 50:260-270. © 2014 Wiley Periodicals, Inc.

LD Score regression distinguishes confounding from polygenicity in genome-wide association studies

Both polygenicity (many small genetic effects) and confounding biases, such as cryptic relatedness and population stratification, can yield an inflated distribution of test statistics in genome-wide association studies (GWAS). However, current methods cannot distinguish between inflation from a true polygenic signal and bias. We have developed an approach, LD Score regression, that quantifies the contribution of each by examining the relationship between test statistics and linkage disequilibrium (LD). The LD Score regression intercept can be used to estimate a more powerful and accurate correction factor than genomic control. We find strong evidence that polygenicity accounts for the majority of the inflation in test statistics in many GWAS of large sample size.

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.

Genes Encoding Vascular Endothelial Growth Factor A (VEGF-A) and VEGF Receptor 2 (VEGFR-2) and Risk for Bronchopulmonary Dysplasia

BACKGROUND

Bronchopulmonary dysplasia (BPD) is one of the main consequences of prematurity, with notably high heritability. Vascular endothelial growth factor A (VEGF-A) and its main receptor, vascular endothelial growth factor receptor 2 (VEGFR-2), have been implicated in the pathogenesis of BPD.

OBJECTIVE

To study whether common polymorphisms of the genes encoding VEGF-A and VEGFR-2 are associated with BPD.

METHODS

In this association study, six tagging single nucleotide polymorphism (tSNPs) for VEGFA and 25 tSNPs for VEGFR2 were genotyped in a prospectively collected, genetically homogeneous discovery population of 160 infants (44 infants with grade 2-3 BPD) born before 30 completed gestational weeks. The replication population of 328 infants included 120 cases of BPD.

RESULTS

VEGFR2 SNP rs4576072 was associated with BPD grade 2-3 with a minor allele frequency in 23.9% of the cases compared to 9.1% in controls (p = 0.0005, odds ratio 3.15, 95% CI: 1.62-6.12) in the discovery population. This association was not observed in the more heterogeneous replication population.

CONCLUSIONS

In line with the results of recent large-scale genetic studies, our findings indicate that common polymorphisms of the genes encoding VEGF-A and VEGFR-2 are not consistently associated with BPD. This finding does not rule out the involvement of VEGFA and VEGFR2 in BPD pathogenesis since, in addition to common variations within the gene region, other mechanisms also play important roles in the regulation of gene function.

A study of genes encoding cytokines (IL6, IL10, TNF), cytokine receptors (IL6R, IL6ST), and glucocorticoid receptor (NR3C1) and susceptibility to bronchopulmonary dysplasia

Background

Bronchopulmonary dysplasia (BPD) is a common chronic lung disease associated with very preterm birth. The major risk factors include lung inflammation and lung immaturity. In addition, genetic factors play an important role in susceptibility to moderate-to-severe BPD. In this study, the aim was to investigate whether common polymorphisms of specific genes that are involved in inflammation or differentiation of the lung have influence on BPD susceptibility.

Methods

Genes encoding interleukin-6 (IL6) and its receptors (IL6R and IL6ST), IL-10 (IL10), tumor necrosis factor (TNF), and glucocorticoid receptor (NR3C1) were assessed for associations with moderate-to-severe BPD susceptibility. Five IL6, nine IL6R, four IL6ST, one IL10, two TNF, and 23 NR3C1 single nucleotide polymorphisms (SNPs) were analyzed in very preterm infants born in northern Finland (56 cases and 197 controls) and Canada (58 cases and 68 controls). IL-6, TNF and gp130 contents in umbilical cord blood, collected from very preterm infants, were studied for associations with the polymorphisms. Epistasis (i.e., interactions between SNPs in BPD susceptibility) was also examined. SNPs showing suggestive associations were analyzed in additional replication populations from Finland (39 cases and 188 controls) and Hungary (29 cases and 40 controls).

Results

None of the studied SNPs were associated with BPD nor were the IL6, TNF or IL6ST SNPs associated with cord blood IL-6, TNF and gp130, respectively. However, epistasis analysis suggested that SNPs in IL6ST and IL10 were associated interactively with risk of BPD in the northern Finnish population; however, this finding did not remain significant after correction for multiple testing and the finding was not replicated in the other populations.

Conclusions

We conclude that the analyzed SNPs within IL6, IL6R, IL6ST, IL10, TNF, and NR3C1 were not associated with BPD. Furthermore, there was no evidence that the studied SNPs directly contribute to the cord blood protein contents.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-014-0120-7) contains supplementary material, which is available to authorized users.