Geographical structure and differential natural selection among North European populations

Complete Mitochondrial DNA Genome Variation in the Swedish Population

The development of complete mitochondrial genome (mitogenome) reference data for inclusion in publicly available population databases is currently underway, and the generation of more high-quality mitogenomes will only enhance the statistical power of this forensically useful locus. To characterize mitogenome variation in Sweden, the mitochondrial DNA (mtDNA) reads from the SweGen whole genome sequencing (WGS) dataset were analyzed. To overcome the interference from low-frequency nuclear mtDNA segments (NUMTs), a 10% variant frequency threshold was applied for the analysis. In total, 934 forensic-quality mitogenome haplotypes were characterized. Almost 45% of the SweGen haplotypes belonged to haplogroup H. Nearly all mitogenome haplotypes (99.1%) were assigned to European haplogroups, which was expected based on previous mtDNA studies of the Swedish population. There were signature northern Swedish and Finnish haplogroups observed in the dataset (e.g., U5b1, W1a), consistent with the nuclear DNA analyses of the SweGen data. The complete mitogenome analysis resulted in high haplotype diversity (0.9996) with a random match probability of 0.15%. Overall, the SweGen mitogenomes provide a large mtDNA reference dataset for the Swedish population and also contribute to the effort to estimate global mitogenome haplotype frequencies.

Molecular Basis for Global Incidence of Pemphigoid Diseases and Differences in Phenotypes

Pemphigoid (Pg) diseases are a group of potentially fatal autoimmune mucocutaneous diseases. They have different clinical phenotypes, involving only the skin or multiple mucous membranes. They occur globally and frequently affect the elderly. The common marker among all variants is the presence of autoantibodies targeting the dermal-epidermal or mucosal-submucosal junctions, or basement membrane zone (BMZ). Four target antigens in the BMZ were studied. These included BPAG1, BPAG2 and subunits of α6 and β4 human integrins. Our objective was to find a molecular basis for the global incidence of Pg diseases and a mechanism that will explain the vast differences in clinical phenotypes and outcomes. All the variants of Pg that were analyzed had a statistically significant association with HLA-DQβ1*03:01 in ten countries on four continents. This explains the reason for global incidence. Prediction models discovered multiple peptides in each of the four antigens that serve as T cell epitopes. These T cell epitopes were shown to bind to HLA-DQβ1*03:01. In addition, structure modelling demonstrated the peptide-HLA complex bound to the T cell receptor. These autoreactive T cells would stimulate B cells to produce specific anti-BMZ autoantibodies. Anti-BMZ autoantibodies with different specificities will produce different phenotypes, which will account for involvement of different tissues and organs in different molecules. The contribution this study makes is that it provides a molecular basis of why a similar disease occurs in different racial groups. Furthermore, it provides the basis for the production of autoantibodies with different specificities, which resultantly produces different phenotypes.

Filaggrin variations are associated with PAH metabolites in urine and DNA alterations in blood

Dermal chemical exposure is common in many professions. The filaggrin protein is important for the skin barrier and variations in the filaggrin gene (FLG) may influence the uptake of chemicals via the skin, and consequently, the degree of systemic effects. The aim of this study was to investigate, in chimney sweeps with occupational exposure to polycyclic aromatic hydrocarbons (PAH) from soot, the influence of variation in FLG on internal PAH dose and DNA alterations, including epigenetic, previously linked to cancer and cardiovascular disease. We used TaqMan PCR to genotype 151 chimney sweeps and 152 controls for four FLG null variants (R501X, R2447X, S3247X and 2282del4) which cause impaired skin barrier, and FLG copy number variation (12th repeat, CNV12) which potentially is beneficial for the skin barrier. The internal dose of PAH was represented by urinary PAH metabolites (e.g. 1-hydroxypyrene and 3-hydroxybenzo[a]pyrene) that we measured by LC-MS/MS. We measured epigenetic alterations (methylation of AHRR and F2RL3) in blood by pyrosequencing; and DNA alterations (telomere length and mitochondrial DNA copy number) by real-time PCR. Hypomethylation of AHRR or F2RL3 is a risk factor for lung cancer and shorter telomere length a risk factor for cardiovascular disease. The frequencies of FLG null were 8.6 and 11.8% (p = 0.35), and CNV12 27.8 and 19.7% (p = 0.09) in chimney sweeps and controls, respectively. We found that among chimney sweeps working predominately with soot sweeping (high PAH exposure), CNV12 carriers had lower concentrations of PAH metabolites in urine compared with non-carriers (median 1-hydroxypyrene = 0.37 vs 0.86 μg/g creatinine respectively; p = 0.025 by linear regression models adjusted for age, BMI and smoking) compared to sweeps not carrying CNV12. Further, FLG null was associated with approximately 2.5% higher methylation of F2RL3 (cg03636183, p = 0.019 after adjustment for exposure group, age, BMI and smoking). FLG null was associated with approximately 7% shorter telomere length (p = 0.015, adjusted model). Our results suggest that FLG variations may influence the dose of PAH in highly exposed workers, possibly via dermal uptake. It also suggests that FLG variation may influence the degree of (epi)genotoxicity in the body. FLG variation is common in the working population and should be considered in risk assessment.

Association of Genetic and Environmental Factors With Autism in a 5-Country Cohort

IMPORTANCE

The origins and development of autism spectrum disorder (ASD) remain unresolved. No individual-level study has provided estimates of additive genetic, maternal, and environmental effects in ASD across several countries.

OBJECTIVE

To estimate the additive genetic, maternal, and environmental effects in ASD.

DESIGN, SETTING, AND PARTICIPANTS

Population-based, multinational cohort study including full birth cohorts of children from Denmark, Finland, Sweden, Israel, and Western Australia born between January 1, 1998, and December 31, 2011, and followed up to age 16 years. Data were analyzed from September 23, 2016 through February 4, 2018.

MAIN OUTCOMES AND MEASURES

Across 5 countries, models were fitted to estimate variance components describing the total variance in risk for ASD occurrence owing to additive genetics, maternal, and shared and nonshared environmental effects.

RESULTS

The analytic sample included 2 001 631 individuals, of whom 1 027 546 (51.3%) were male. Among the entire sample, 22 156 were diagnosed with ASD. The median (95% CI) ASD heritability was 80.8% (73.2%-85.5%) for country-specific point estimates, ranging from 50.9% (25.1%-75.6%) (Finland) to 86.8% (69.8%-100.0%) (Israel). For the Nordic countries combined, heritability estimates ranged from 81.2% (73.9%-85.3%) to 82.7% (79.1%-86.0%). Maternal effect was estimated to range from 0.4% to 1.6%. Estimates of genetic, maternal, and environmental effects for autistic disorder were similar with ASD.

CONCLUSIONS AND RELEVANCE

Based on population data from 5 countries, the heritability of ASD was estimated to be approximately 80%, indicating that the variation in ASD occurrence in the population is mostly owing to inherited genetic influences, with no support for contribution from maternal effects. The results suggest possible modest differences in the sources of ASD risk between countries.

Outcomes and Risk Factors in Patients with Multiple Primary Melanomas

The incidence and patient survival rates of melanoma have increased over the last several decades, with a growing population of patients who develop multiple primary melanomas (MPMs). To determine risk factors for developing MPMs and compare the survival of patients with MPMs to those with single primary melanomas, a prospective, multidisciplinary database of patients with melanoma at a single tertiary care institution was retrospectively reviewed. From 1985 to 2013, 6,963 patients with single primary melanomas and 305 patients with MPMs were identified. Mean follow-up was 8.3 ± 6.3 years for patients with single primary melanomas and 8.8 ± 5.9 years for patients with MPMs. Risk of developing multiple melanomas increased with age at diagnosis of first melanoma (hazard ratio [HR] = 1.20 for a 10-year increase in age, 95% confidence interval [CI] = 1.11-1.29, P < 0.001), male sex (HR = 1.44, 95% CI = 1.12-1.84, P = 0.005), and white race (HR = 3.07, 95% CI = 1.45-6.51). Patients with invasive MPMs had increased risk of melanoma-specific death both before (HR = 1.47, 95% CI = 1.0-2.2) and after adjusting for age, sex, site, race, family history of melanoma, personal history of other cancer, and Surveillance, Epidemiology, and End Results Program (SEER) stage (HR = 1.44, 95% CI = 0.95-2.2); however, this result did not reach statistical significance.

A distinct strain of Arsenophonus symbiont decreases insecticide resistance in its insect host

Symbiotic bacteria are important drivers of phenotypic diversity in insects. One of the widespread symbionts to have emerged belongs to the genus Arsenophonus, however, its biological functions in most host insects remain entirely unknown. Here we report two distinct Arsenophonus strains in the brown planthopper (BPH), Nilaparvata lugens, a major pest insect in Asian countries that causes significant economic damage through rice crop destruction. Genomic resequencing data suggested that one Arsenophonus strain (S-type) negatively affected the insecticide resistance of the host. Indeed, replacement of the resident Arsenophonus with the S-type Arsenophonus significantly decreased host insecticide resistance. Transcriptome and metabolome analysis revealed down-regulation of xenobiotic metabolism and increased amino acid accumulation in the S-type Arsenophonus infected host. This study demonstrates how a symbiont-mediated phenotypic change can occur. The results of this study will aid in developing strategies that work through imposing an ecological disadvantage on insect pests, which will be of great value for pest control in agricultural industry.

A direct test of the diathesis-stress model for depression

The diathesis-stress theory for depression states that the effects of stress on the depression risk are dependent on the diathesis or vulnerability, implying multiplicative interactive effects on the liability scale. We used polygenic risk scores for major depressive disorder (MDD) calculated from the results of the most recent analysis from the Psychiatric Genomics Consortium as a direct measure of the vulnerability for depression in a sample of 5221 individuals from 3083 families. In the same we also had measures of stressful life events and social support and a depression symptom score, as well as DSM-IV MDD diagnoses for most individuals. In order to estimate the variance in depression explained by the genetic vulnerability, the stressors and their interactions, we fitted linear mixed models controlling for relatedness for the whole sample as well as stratified by sex. We show a significant interaction of the polygenic risk scores with personal life events (0.12% of variance explained, P-value=0.0076) contributing positively to the risk of depression. Additionally, our results suggest possible differences in the aetiology of depression between women and men. In conclusion, our findings point to an extra risk for individuals with combined vulnerability and high number of reported personal life events beyond what would be expected from the additive contributions of these factors to the liability for depression, supporting the multiplicative diathesis-stress model for this disease.

GREB1 genetic variants are associated with bone mineral density in Caucasians

Gaining an understanding of factors contributing to bone quality is key to the development of effective preventative treatments for osteoporosis and reduction in osteoporotic fractures. Oestrogen is a strong regulator of bone remodelling which maintains skeletal structural integrity. The growth regulation by oestrogen in breast cancer 1 (GREB1) gene, with an as yet undefined function, is an early response gene in the oestrogen-regulated pathway. Suggestive evidence of linkage with bone mineral density (BMD) variation has been reported with D2S168, located telomeric of GREB1. The aim of this study was to determine if genetic variation within GREB1 was associated with BMD variation at two sites with high fracture rates-the lumbar spine (LS) and the femoral neck (FN). Informative GREB1 single-nucleotide polymorphisms (SNPs) (n = 12) were selected for genotyping and tested for association in a family-based dataset (n = 508 individuals from 229 families). Significantly associated SNPs were tested further in a postmenopausal dataset from the same geographic region (n = 477 individuals). One intronic SNP, rs5020877, was significantly associated with LS and FN BMD in the family-based dataset (P ≤ 0.005). The association was not observed in the postmenopausal dataset (P > 0.017); however, rs10929757 was significantly associated with FN BMD (P = 0.006). Markers, rs5020877 and rs10929757, were constituent SNPs in one GREB1 linkage disequilibrium block, although not historically correlated (r ² = 0.07). Our findings suggest that GREB1 is a novel gene target for osteoporosis genetics and needs to be investigated further.

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.

PROX1 gene CC genotype as a major determinant of early onset of type 2 diabetes in slavic study participants from Action in Diabetes and Vascular Disease: Preterax and Diamicron MR Controlled Evaluation study

BACKGROUND

The prevalence of diabetic nephropathy varies according to ethnicity. Environmental as well as genetic factors contribute to the heterogeneity in the presentation of diabetic nephropathy. Our objective was to evaluate this heterogeneity within the Caucasian population.

METHODS

The geo-ethnic origin of the 3409 genotyped Caucasian type 2 diabetes (T2D) patients of Action in Diabetes and Vascular Disease: Preterax and Diamicron MR Controlled Evaluation was determined using principal component analysis. Genome-wide association studies analyses of age of onset of T2D were performed for geo-ethnic groups separately and combined.

RESULTS

The first principal component separated the Caucasian study participants into Slavic and Celtic ethnic origins. Age of onset of diabetes was significantly lower in Slavic patients (P = 7.3 × 10), whereas the prevalence of hypertension (P = 4.9 × 10) and albuminuria (5.1 × 10) were significantly higher. Age of onset of T2D and albuminuria appear to have an important genetic component as the values of these traits were also different between Slavic and Celtic individuals living in the same countries. Common and geo-ethnic-specific loci were found to be associated to age of onset of diabetes. Among the latter, the PROX1/PROX1-AS1 genes (rs340841) had the highest impact. Single-nucleotide polymorphism rs340841 CC genotype was associated with a 4.4 year earlier onset of T2D in Slavic patients living or not in countries with predominant Slavic populations.

CONCLUSION

These results reveal the presence of distinct genetic architectures between Caucasian ethnic groups that likely have clinical relevance, among them PROX1 gene is a strong candidate of early onset of diabetes with variations depending on ethnicity.

Genome-wide association study of working memory brain activation

In a population-based genome-wide association (GWA) study of n-back working memory task-related brain activation, we extracted the average percent BOLD signal change (2-back minus 0-back) from 46 regions-of-interest (ROIs) in functional MRI scans from 863 healthy twins and siblings. ROIs were obtained by creating spheres around group random effects analysis local maxima, and by thresholding a voxel-based heritability map of working memory brain activation at 50%. Quality control for test-retest reliability and heritability of ROI measures yielded 20 reliable (r>0.7) and heritable (h2>20%) ROIs. For GWA analysis, the cohort was divided into a discovery (n=679) and replication (n=97) sample. No variants survived the stringent multiple-testing-corrected genome-wide significance threshold (p<4.5×10-9), or were replicated (p<0.0016), but several genes were identified that are worthy of further investigation. A search of 529,379 genomic markers resulted in discovery of 31 independent single nucleotide polymorphisms (SNPs) associated with BOLD signal change at a discovery level of p<1×10-5. Two SNPs (rs7917410 and rs7672408) were associated at a significance level of p<1×10-7. Only one, most strongly affecting BOLD signal change in the left supramarginal gyrus (R2=5.5%), had multiple SNPs associated at p<1×10-5 in linkage disequilibrium with it, all located in and around the BANK1 gene. BANK1 encodes a B-cell-specific scaffold protein and has been shown to negatively regulate CD40-mediated AKT activation. AKT is part of the dopamine-signaling pathway, suggesting a mechanism for the involvement of BANK1 in the BOLD response to working memory. Variants identified here may be relevant to (the susceptibility to) common disorders affecting brain function.

A genome-wide association study of essential hypertension in an Australian population using a DNA pooling approach

Despite the success of genome-wide association studies (GWAS) in detecting genetic loci involved in complex traits, few susceptibility genes have been detected for essential hypertension (EH). We aimed to use pooled DNA GWAS approach to identify and validate novel genomic loci underlying EH susceptibility in an Australian case-control population. Blood samples and questionnaires detailing medical history, blood pressure, and prescribed medications were collected for 409 hypertensives and 409 age-, sex- and ethnicity-matched normotensive controls. Case and control DNA were pooled in quadruplicate and hybridized to Illumina 1 M-Duo arrays. Allele frequencies agreed with those reported in reference data and known EH association signals were represented in the top-ranked SNPs more frequently than expected by chance. Validation showed that pooled DNA GWAS gave reliable estimates of case and control allele frequencies. Although no markers reached Bonferroni-corrected genome-wide significance levels (5.0 × 10^-8), the top marker rs34870220 near ASGR1 approached significance (p = 4.32 × 10^-7), as did several candidate loci (p < 1 × 10^-6) on chromosomes 2, 4, 6, 9, 12, and 17. Four markers (located in or near genes NHSL1, NKFB1, GLI2, and LRRC10) from the top ten ranked SNPs were individually genotyped in pool samples and were tested for association between cases and controls using the χ ² test. Of these, rs1599961 (NFKB1) and rs12711538 (GLI2) showed significant difference between cases and controls (p < 0.01). Additionally, four top-ranking markers within NFKB1 were found to be in LD, suggesting a single strong association signal for this gene.

Germline polymorphisms in an enhancer of PSIP1 are associated with progression-free survival in epithelial ovarian cancer

Women with epithelial ovarian cancer (EOC) are usually treated with platinum/taxane therapy after cytoreductive surgery but there is considerable inter-individual variation in response. To identify germline single-nucleotide polymorphisms (SNPs) that contribute to variations in individual responses to chemotherapy, we carried out a multi-phase genome-wide association study (GWAS) in 1,244 women diagnosed with serous EOC who were treated with the same first-line chemotherapy, carboplatin and paclitaxel. We identified two SNPs (rs7874043 and rs72700653) in TTC39B (best P=7x10-5, HR=1.90, for rs7874043) associated with progression-free survival (PFS). Functional analyses show that both SNPs lie in a putative regulatory element (PRE) that physically interacts with the promoters of PSIP1, CCDC171 and an alternative promoter of TTC39B. The C allele of rs7874043 is associated with poor PFS and showed increased binding of the Sp1 transcription factor, which is critical for chromatin interactions with PSIP1. Silencing of PSIP1 significantly impaired DNA damage-induced Rad51 nuclear foci and reduced cell viability in ovarian cancer lines. PSIP1 (PC4 and SFRS1 Interacting Protein 1) is known to protect cells from stress-induced apoptosis, and high expression is associated with poor PFS in EOC patients. We therefore suggest that the minor allele of rs7874043 confers poor PFS by increasing PSIP1 expression.

Characterization of the biological processes shaping the genetic structure of the Italian population

Background

The genetic structure of human populations is the outcome of the combined action of different processes such as demographic dynamics and natural selection. Several efforts toward the characterization of population genetic architectures and the identification of adaptation signatures were recently made. In this study, we provide a genome-wide depiction of the Italian population structure and the analysis of the major determinants of the current existing genetic variation.

Results

We defined and characterized 210 genomic loci associated with the first Principal Component calculated on the Italian genotypic data and correlated to the North–south genetic gradient. Using a gene-enrichment approach we identified the immune function as primarily involved in the Italian population differentiation and we described a locus on chromosome 13 showing combined evidence of North–south diversification in allele frequencies and signs of recent positive selection. In this region our bioinformatics analysis pinpointed an uncharacterized long intergenic non-coding (lincRNA), whose expression appeared specific for immune-related tissues suggesting its relevance for the immune function.

Conclusions

Our study, combining population genetic analyses with biological insights provides a description of the Italian genetic structure that in future could contribute to the evaluation of complex diseases risk in the population context.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0293-x) contains supplementary material, which is available to authorized users.

Individual Aesthetic Preferences for Faces Are Shaped Mostly by Environments, Not Genes

Although certain characteristics of human faces are broadly considered more attractive (e.g., symmetry, averageness), people also routinely disagree with each other on the relative attractiveness of faces. That is, to some significant degree, beauty is in the "eye of the beholder." Here, we investigate the origins of these individual differences in face preferences using a twin design, allowing us to estimate the relative contributions of genetic and environmental variation to individual face attractiveness judgments or face preferences. We first show that individual face preferences (IP) can be reliably measured and are readily dissociable from other types of attractiveness judgments (e.g., judgments of scenes, objects). Next, we show that individual face preferences result primarily from environments that are unique to each individual. This is in striking contrast to individual differences in face identity recognition, which result primarily from variations in genes [1]. We thus complete an etiological double dissociation between two core domains of social perception (judgments of identity versus attractiveness) within the same visual stimulus (the face). At the same time, we provide an example, rare in behavioral genetics, of a reliably and objectively measured behavioral characteristic where variations are shaped mostly by the environment. The large impact of experience on individual face preferences provides a novel window into the evolution and architecture of the social brain, while lending new empirical support to the long-standing claim that environments shape individual notions of what is attractive.

Potential Signals of Natural Selection in the Top Risk Loci for Coronary Artery Disease: 9p21 and 10q11

BACKGROUND

Coronary artery disease (CAD) is a complex disease and the leading cause of death in the world. Populations of different ancestry do not always share the same risk markers. Natural selective processes may be the cause of some of the population differences detected for specific risk mutations.

OBJECTIVE

In this study, 384 single nucleotide polymorphisms (SNPs) located in four genomic regions associated with CAD (1p13, 1q41, 9p21 and 10q11) are analysed in a set of 19 populations from Europe, Middle East and North Africa and also in Asian and African samples from the 1000 Genomes Project. The aim of this survey is to explore for the first time whether the genetic variability in these genomic regions is better explained by demography or by natural selection.

RESULTS

The results indicate significant differences in the structure of genetic variation and in the LD patterns among populations that probably explain the population disparities found in markers of susceptibility to CAD.

CONCLUSIONS

The results are consistent with potential signature of positive selection in the 9p21 region and of balancing selection in the 9p21 and 10q11. Specifically, in Europe three CAD risk markers in the 9p21 region (rs9632884, rs1537371 and rs1333042) show consistent signals of positive selection. The results of this study are consistent with a potential selective role of CAD in the configuration of genetic diversity in current human populations.

The European Genome-phenome Archive of human data consented for biomedical research

The European Genome-phenome Archive (EGA) is a permanent archive that promotes distribution and sharing of genetic and phenotype data consented for specific approved uses, but not fully open public distribution. The EGA follows strict protocols for information management, data storage, security and dissemination. Authorized access to the data is managed in partnership with the data providing organizations. The EGA includes major reference data collections for human genetics research.

Estimation of Genetic Relationships Between Individuals Across Cohorts and Platforms: Application to Childhood Height

Combining genotype data across cohorts increases power to estimate the heritability due to common single nucleotide polymorphisms (SNPs), based on analyzing a Genetic Relationship Matrix (GRM). However, the combination of SNP data across multiple cohorts may lead to stratification, when for example, different genotyping platforms are used. In the current study, we address issues of combining SNP data from different cohorts, the Netherlands Twin Register (NTR) and the Generation R (GENR) study. Both cohorts include children of Northern European Dutch background (N = 3102 + 2826, respectively) who were genotyped on different platforms. We explore imputation and phasing as a tool and compare three GRM-building strategies, when data from two cohorts are (1) just combined, (2) pre-combined and cross-platform imputed and (3) cross-platform imputed and post-combined. We test these three strategies with data on childhood height for unrelated individuals (N = 3124, average age 6.7 years) to explore their effect on SNP-heritability estimates and compare results to those obtained from the independent studies. All combination strategies result in SNP-heritability estimates with a standard error smaller than those of the independent studies. We did not observe significant difference in estimates of SNP-heritability based on various cross-platform imputed GRMs. SNP-heritability of childhood height was on average estimated as 0.50 (SE = 0.10). Introducing cohort as a covariate resulted in ≈2 % drop. Principal components (PCs) adjustment resulted in SNP-heritability estimates of about 0.39 (SE = 0.11). Strikingly, we did not find significant difference between cross-platform imputed and combined GRMs. All estimates were significant regardless the use of PCs adjustment. Based on these analyses we conclude that imputation with a reference set helps to increase power to estimate SNP-heritability by combining cohorts of the same ethnicity genotyped on different platforms. However, important factors should be taken into account such as remaining cohort stratification after imputation and/or phenotypic heterogeneity between and within cohorts. Whether one should use imputation, or just combine the genotype data, depends on the number of overlapping SNPs in relation to the total number of genotyped SNPs for both cohorts, and their ability to tag all the genetic variance related to the specific trait of interest.

Extensive genomic variability of knops blood group polymorphisms is associated with sickle cell disease in Africa

Sickle cell disease (SCD) is a multisystem disorder characterized by chronic hemolytic anemia, vaso-occlusive crises, and marked variability in disease severity. Patients require transfusions to manage disease complications, with complements, directed by complement regulatory genes (CR1) and its polymorphisms, implicated in the development of alloantibodies. We hypothesize that CR1 polymorphisms affect complement regulation and function, leading to adverse outcome in SCD. To this end, we determined the genomic diversity of complement regulatory genes by examining single nucleotide polymorphisms associated with Knops blood group antigens. Genomic DNA samples from 130 SCD cases and 356 control Africans, 331 SCD cases and 497 control African Americans, and 254 Caucasians were obtained and analyzed, utilizing a PCR-RFLP (polymerase chain reaction-restriction fragment length polymorphism) assay. Analyzing for ethnic diversity, we found significant differences in the genotypic and allelic frequencies of Sl1/Sl2 (rs17047661) and McCa/b (rs17047660) polymorphisms between Africans, African Americans, and Caucasians (P < 0.05). The homozygote mutant variants had significantly higher frequencies in Africans and African Americans but were insignificant in Caucasians (80.2% and 59.6% vs 5.9% for Sl1/2; and 36% and 24% vs 1.8% for McCa/b). With SCD, we did not detect any difference among cases and controls either in Africa or in the United States. However, we found significant difference in genotypic (P < 0.0001) and allelic frequencies (P < 0.0001) of Sl1/Sl2 (rs17047661) and McCa/b (rs17047660) polymorphisms between SCD groups from Africa and the United States. There was no difference in haplotype frequencies of these polymorphisms among or between groups. The higher frequency of CR1 homozygote mutant variants in Africa but not United States indicates a potential pathogenic role, possibly associated with complicated disease pathophysiology in the former and potentially protective in the latter. The difference between sickle cell groups suggests potential genetic drift or founder effect imposed on the disease in the United States, but not in Africa, and a possible confirmation of the ancestral susceptibility hypothesis. The lower haplotype frequencies among sickle cell and control populations in the United States may be due to the admixture and the dilution of African genetic ancestry in the African American population.

Educational attainment influences levels of homozygosity through migration and assortative mating

Individuals with a higher education are more likely to migrate, increasing the chance of meeting a spouse with a different ancestral background. In this context, the presence of strong educational assortment can result in greater ancestry differences within more educated spouse pairs, while less educated individuals are more likely to mate with someone with whom they share more ancestry. We examined the association between educational attainment and F_roh (= the proportion of the genome consisting of runs of homozygosity [ROHs]) in ~2,000 subjects of Dutch ancestry. The subjects’ own educational attainment showed a nominally significant negative association with F_roh (p = .045), while the contribution of parental education to offspring F_roh was highly significant (father: p < 10^-5; mother: p = 9×10^-5), with more educated parents having offspring with fewer ROHs. This association was significantly and fully mediated by the physical distance between parental birthplaces (paternal education: p_mediation = 2.4 × 10^-4; maternal education: p_mediation = 2.3 × 10^-4), which itself was also significantly associated with F_roh (p = 9 × 10^-5). Ancestry-informative principal components from the offspring showed a significantly decreasing association with geography as parental education increased, consistent with the significantly higher migration rates among more educated parents. Parental education also showed a high spouse correlation (Spearman’s ρ = .66, p = 3 × 10^-262). We show that less educated parents are less likely to mate with the more mobile parents with a higher education, creating systematic differences in homozygosity due to ancestry differences not directly captured by ancestry-informative principal components (PCs). Understanding how behaviors influence the genomic structure of a population is highly valuable for studies on the genetic etiology of behavioral, cognitive, and social traits.

A hidden Markov model for investigating recent positive selection through haplotype structure

Recent positive selection can increase the frequency of an advantageous mutant rapidly enough that a relatively long ancestral haplotype will be remained intact around it. We present a hidden Markov model (HMM) to identify such haplotype structures. With HMM identified haplotype structures, a population genetic model for the extent of ancestral haplotypes is then adopted for parameter inference of the selection intensity and the allele age. Simulations show that this method can detect selection under a wide range of conditions and has higher power than the existing frequency spectrum-based method. In addition, it provides good estimate of the selection coefficients and allele ages for strong selection. The method analyzes large data sets in a reasonable amount of running time. This method is applied to HapMap III data for a genome scan, and identifies a list of candidate regions putatively under recent positive selection. It is also applied to several genes known to be under recent positive selection, including the LCT, KITLG and TYRP1 genes in Northern Europeans, and OCA2 in East Asians, to estimate their allele ages and selection coefficients.

Genome-wide scan of 29,141 African Americans finds no evidence of directional selection since admixture

The extent of recent selection in admixed populations is currently an unresolved question. We scanned the genomes of 29,141 African Americans and failed to find any genome-wide-significant deviations in local ancestry, indicating no evidence of selection influencing ancestry after admixture. A recent analysis of data from 1,890 African Americans reported that there was evidence of selection in African Americans after their ancestors left Africa, both before and after admixture. Selection after admixture was reported on the basis of deviations in local ancestry, and selection before admixture was reported on the basis of allele-frequency differences between African Americans and African populations. The local-ancestry deviations reported by the previous study did not replicate in our very large sample, and we show that such deviations were expected purely by chance, given the number of hypotheses tested. We further show that the previous study's conclusion of selection in African Americans before admixture is also subject to doubt. This is because the FST statistics they used were inflated and because true signals of unusual allele-frequency differences between African Americans and African populations would be best explained by selection that occurred in Africa prior to migration to the Americas.

A powerful tool for genome analysis in maize: development and evaluation of the high density 600 k SNP genotyping array

Background

High density genotyping data are indispensable for genomic analyses of complex traits in animal and crop species. Maize is one of the most important crop plants worldwide, however a high density SNP genotyping array for analysis of its large and highly dynamic genome was not available so far.

Results

We developed a high density maize SNP array composed of 616,201 variants (SNPs and small indels). Initially, 57 M variants were discovered by sequencing 30 representative temperate maize lines and then stringently filtered for sequence quality scores and predicted conversion performance on the array resulting in the selection of 1.2 M polymorphic variants assayed on two screening arrays. To identify high-confidence variants, 285 DNA samples from a broad genetic diversity panel of worldwide maize lines including the samples used for sequencing, important founder lines for European maize breeding, hybrids, and proprietary samples with European, US, semi-tropical, and tropical origin were used for experimental validation. We selected 616 k variants according to their performance during validation, support of genotype calls through sequencing data, and physical distribution for further analysis and for the design of the commercially available Affymetrix® Axiom® Maize Genotyping Array. This array is composed of 609,442 SNPs and 6,759 indels. Among these are 116,224 variants in coding regions and 45,655 SNPs of the Illumina® MaizeSNP50 BeadChip for study comparison. In a subset of 45,974 variants, apart from the target SNP additional off-target variants are detected, which show only a minor bias towards intermediate allele frequencies. We performed principal coordinate and admixture analyses to determine the ability of the array to detect and resolve population structure and investigated the extent of LD within a worldwide validation panel.

Conclusions

The high density Affymetrix® Axiom® Maize Genotyping Array is optimized for European and American temperate maize and was developed based on a diverse sample panel by applying stringent quality filter criteria to ensure its suitability for a broad range of applications. With 600 k variants it is the largest currently publically available genotyping array in crop species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-823) contains supplementary material, which is available to authorized users.

HLA associations in schizophrenia: are we re-discovering the wheel?

Associations between human leukocyte antigen (HLA) polymorphisms on chromosome 6p and schizophrenia (SZ) risk have been evaluated for over five decades. Numerous case-control studies from the candidate gene era analyzed moderately sized samples and reported nominally significant associations with several loci in the HLA region (sample sizes, n = 100-400). The risk conferred by individual alleles was modest (odds ratios < 2.0). The basis for the associations could not be determined, though connections with known immune and auto-immune abnormalities in SZ were postulated. Interest in the HLA associations has re-emerged following several recent genome-wide association studies (GWAS); which utilized 10- to 100-fold larger samples and also identified associations on the short arm of chromosome 6. Unlike the earlier candidate gene studies, the associations are statistically significant following correction for multiple comparisons. Like the earlier studies; they have modest effect sizes, raising questions about their utility in risk prediction or pathogenesis research. In this review, we summarize the GWAS and reflect on possible bases for the associations. Suggestions for future research are discussed. We favor, in particular; efforts to evaluate local population sub-structure as well as further evaluation of immune-related variables in future studies.

The GenoChip: a new tool for genetic anthropology

The Genographic Project is an international effort aimed at charting human migratory history. The project is nonprofit and nonmedical, and, through its Legacy Fund, supports locally led efforts to preserve indigenous and traditional cultures. Although the first phase of the project was focused on uniparentally inherited markers on the Y-chromosome and mitochondrial DNA (mtDNA), the current phase focuses on markers from across the entire genome to obtain a more complete understanding of human genetic variation. Although many commercial arrays exist for genome-wide single-nucleotide polymorphism (SNP) genotyping, they were designed for medical genetic studies and contain medically related markers that are inappropriate for global population genetic studies. GenoChip, the Genographic Project’s new genotyping array, was designed to resolve these issues and enable higher resolution research into outstanding questions in genetic anthropology. The GenoChip includes ancestry informative markers obtained for over 450 human populations, an ancient human (Saqqaq), and two archaic hominins (Neanderthal and Denisovan) and was designed to identify all known Y-chromosome and mtDNA haplogroups. The chip was carefully vetted to avoid inclusion of medically relevant markers. To demonstrate its capabilities, we compared the F_ST distributions of GenoChip SNPs to those of two commercial arrays. Although all arrays yielded similarly shaped (inverse J) F_ST distributions, the GenoChip autosomal and X-chromosomal distributions had the highest mean F_ST, attesting to its ability to discern subpopulations. The chip performances are illustrated in a principal component analysis for 14 worldwide populations. In summary, the GenoChip is a dedicated genotyping platform for genetic anthropology. With an unprecedented number of approximately 12,000 Y-chromosomal and approximately 3,300 mtDNA SNPs and over 130,000 autosomal and X-chromosomal SNPs without any known health, medical, or phenotypic relevance, the GenoChip is a useful tool for genetic anthropology and population genetics.

Interrogation of the platelet-derived growth factor receptor alpha locus and corneal astigmatism in Australians of Northern European ancestry: results of a genome-wide association study

PURPOSE

Corneal astigmatism is a common eye disorder characterized by irregularities in corneal curvature. Recently, the rs7677751 single nucleotide polymorphism (SNP) at the platelet-derived growth factor receptor alpha (PDGFRA) locus was found to be associated with corneal astigmatism in people of Asian ancestry. In the present study, we sought to replicate this finding and identify other genetic markers of corneal astigmatism in an Australian population of Northern European ancestry.

METHODS

Data from two cohorts were included in this study. The first cohort consisted of 1,013 individuals who were part of the Western Australian Pregnancy Cohort (Raine) Study: 20-year follow-up Eye Study. The second cohort comprised 1,788 individuals of 857 twin families who were recruited through the Twins Eye Study in Tasmania and the Brisbane Adolescent Twin Study. Corneal astigmatism was calculated as the absolute difference between the keratometry readings in two meridians, and genotype data were extracted from genome-wide arrays. Initially, each cohort was analyzed separately, before being combined for meta- and subsequent genome-wide pathway analysis.

RESULTS

Following meta-analysis, SNP rs7677751 at the PDGFRA locus had a combined p=0.32. No variant was found to be statistically significantly associated with corneal astigmatism at the genome-wide level (p<5.0×10(-8)). The SNP with strongest association was rs1164064 (p=1.86×10(-6)) on chromosome 3q13. Gene-based pathway analysis identified a significant association between the Gene Ontology "segmentation" (GO:0035282) pathway, corrected p=0.009.

CONCLUSIONS

Our data suggest that the PDGFRA locus does not transfer a major risk of corneal astigmatism in people of Northern European ancestry. Better-powered studies are required to validate the novel putative findings of our study.

GWAS of DNA Methylation Variation Within Imprinting Control Regions Suggests Parent-of-Origin Association

Imprinting control regions (ICRs) play a fundamental role in establishing and maintaining the non-random monoallelic expression of certain genes, via common regulatory elements such as non-coding RNAs and differentially methylated regions (DMRs) of DNA. We recently surveyed DNA methylation levels within four ICRs (H19-ICR, IGF2-DMR, KvDMR, and NESPAS-ICR) in whole-blood genomic DNA from 128 monozygotic (MZ) and 128 dizygotic (DZ) human twin pairs. Our analyses revealed high individual variation and intra-domain covariation in methylation levels across CpGs and emphasized the interaction between epigenetic variation and the underlying genetic sequence in a parent-of-origin fashion. Here, we extend our analysis to conduct two genome-wide screenings of single nucleotide polymorphisms (SNPs) underlying either intra-domain covariation or parent-of-origin-dependent association with methylation status at individual CpG sites located within ICRs. Although genome-wide significance was not surpassed due to sample size limitations, the most significantly associated SNPs found through multiple-trait genome-wide association (MQFAM) included the previously described rs10732516, which is located in the vicinity of the H19-ICR. Similarly, we identified an association between rs965808 and methylation status within the NESPAS-ICR. This SNP is positioned within an intronic region of the overlapping genes GNAS and GNAS-AS1, which are imprinted genes regulated by the NESPAS-ICR. Sixteen other SNPs located in regions apart from the analyzed regions displayed suggestive association with intra-domain methylation. Additionally, we identified 13 SNPs displaying parent-of-origin association with individual methylation sites through family-based association testing. In this exploratory study, we show the value and feasibility of using alternative GWAS approaches in the study of the interaction between epigenetic state and genetic sequence within imprinting regulatory domains. Despite the relatively small sample size, we identified a number of SNPs displaying suggestive association either in a domain-wide or in a parent-of-origin fashion. Nevertheless, these associations will require future experimental validation or replication in larger and independent samples.

Genome-wide association study of a quantitative disordered gambling trait

Disordered gambling is a moderately heritable trait, but the underlying genetic basis is largely unknown. We performed a genome-wide association study (GWAS) for disordered gambling using a quantitative factor score in 1312 twins from 894 Australian families. Association was conducted for 2 381 914 single-nucleotide polymorphisms (SNPs) using the family-based association test in Merlin followed by gene and pathway enrichment analyses. Although no SNP reached genome-wide significance, six achieved P-values < 1 × 10(-5) with variants in three genes (MT1X, ATXN1 and VLDLR) implicated in disordered gambling. Secondary case-control analyses found two SNPs on chromosome 9 (rs1106076 and rs12305135 near VLDLR) and rs10812227 near FZD10 on chromosome 12 to be significantly associated with lifetime Diagnostic and Statistical Manual of Mental Disorders, fourth edition pathological gambling and South Oaks Gambling Screen classified probable pathological gambling status. Furthermore, several addiction-related pathways were enriched for SNPs associated with disordered gambling. Finally, gene-based analysis of 24 candidate genes for dopamine agonist-induced gambling in individuals with Parkinson's disease suggested an enrichment of SNPs associated with disordered gambling. We report the first GWAS of disordered gambling. While further replication is required, the identification of susceptibility loci and biological pathways will be important in characterizing the biological mechanisms that underpin disordered gambling.

Population structure, migration, and diversifying selection in the Netherlands

Genetic variation in a population can be summarized through principal component analysis (PCA) on genome-wide data. PCs derived from such analyses are valuable for genetic association studies, where they can correct for population stratification. We investigated how to capture the genetic population structure in a well-characterized sample from the Netherlands and in a worldwide data set and examined whether (1) removing long-range linkage disequilibrium (LD) regions and LD-based SNP pruning significantly improves correlations between PCs and geography and (2) whether genetic differentiation may have been influenced by migration and/or selection. In the Netherlands, three PCs showed significant correlations with geography, distinguishing between: (1) North and South; (2) East and West; and (3) the middle-band and the rest of the country. The third PC only emerged with minimized LD, which also significantly increased correlations with geography for the other two PCs. In addition to geography, the Dutch North-South PC showed correlations with genome-wide homozygosity (r=0.245), which may reflect a serial-founder effect due to northwards migration, and also with height (♂: r=0.142, ♀: r=0.153). The divergence between subpopulations identified by PCs is partly driven by selection pressures. The first three PCs showed significant signals for diversifying selection (545 SNPs - the majority within 184 genes). The strongest signal was observed between North and South for the functional SNP in HERC2 that determines human blue/brown eye color. Thus, this study demonstrates how to increase ancestry signals in a relatively homogeneous population and how those signals can reveal evolutionary history.

Anisotropic isolation by distance: the main orientations of human genetic differentiation

Genetic differentiation among human populations is greatly influenced by geography due to the accumulation of local allele frequency differences. However, little is known about the possibly different increment of genetic differentiation along the different geographical axes (north–south, east–west, etc.). Here, we provide new methods to examine the asymmetrical patterns of genetic differentiation. We analyzed genome-wide polymorphism data from populations in Africa (n = 29), Asia (n = 26), America (n = 9), and Europe (n = 38), and we found that the major orientations of genetic differentiation are north–south in Europe and Africa, and east–west in Asia, but no preferential orientation was found in the Americas. Additionally, we showed that the localization of the individual geographic origins based on single nucleotide polymorphism data was not equally precise along all orientations. Confirming our findings, we obtained that, in each continent, the orientation along which the precision is maximal corresponds to the orientation of maximum differentiation. Our results have implications for interpreting human genetic variation in terms of isolation by distance and spatial range expansion processes. In Europe, for instance, the precise northnorthwest–southsoutheast axis of main European differentiation cannot be explained by a simple Neolithic demic diffusion model without admixture with the local populations because in that case the orientation of greatest differentiation should be perpendicular to the direction of expansion. In addition to humans, anisotropic analyses can guide the description of genetic differentiation for other organisms and provide information on expansions of invasive species or the processes of plant dispersal.

Genome-wide scan with nearly 700,000 SNPs in two Sardinian sub-populations suggests some regions as candidate targets for positive selection

This paper explores the genetic structure and signatures of natural selection in different sub-populations from the Island of Sardinia, exploiting information from nearly 700,000 autosomal SNPs genotyped with the Affymetrix Genome-Wide Human SNP 6.0 Array. The genetic structure of the Sardinian population and its position within the context of other Mediterranean and European human groups were investigated in depth by comparing our data with publicly available data sets. Principal components and admixture analyses suggest a clustering of the examined samples in two significantly differentiated sub-populations (Ogliastra and Southern Sardinia), as confirmed by AMOVA (F(ST)=0.011; P<0.001). Differentiation of these sub-populations was still evident when they were pooled together with supplementary Sardinian samples from HGDP and compared with several other European, North-African and Near Eastern populations, confirming the uniqueness of the Sardinian genetic background. Moreover, by applying several statistical approaches aimed at assessing differences at the SNP level, the highest differentiated genomic regions between Ogliastra and Southern Sardinia were thus investigated via an extended haplotype homozygosity (EHH)-based test to point out potential selective sweeps. Using this approach, 40 genomic regions were detected, with significant differences between Ogliastra and Southern Sardinia. These regions were subsequently investigated using a long-range haplotype test, which found significant REHH values for SNPs rs11070188 and rs11070192 in the Ogliastra sub-population. In the light of these results and the overlap of the different computed statistics, the region encompassing these loci can be considered a strong candidate to have undergone selective pressure in Ogliastra.

Genetic variants near PDGFRA are associated with corneal curvature in Australians

PURPOSE

Irregularity in the corneal curvature (CC) is highly associated with various eye disorders such as keratoconus and myopia. The sample had limited power to find genomewide significant (5 × 10(-8)) hits but good power for replication. Thus, an attempt was made to test whether alleles in the FRAP1 and PDGFRA genes, recently found to be associated with CC in Asian populations, also influence CC in Australians of North European ancestry. Results of initial genomewide association studies (GWAS) for CC in Australians were also reported.

METHODS

Two population-based cohorts of 1788 Australian twins and their families, as well as 1013 individuals from a birth cohort from Western Australia, were genotyped using genomewide arrays. Following separate individual analysis and quality control, the results from each cohort underwent meta-analysis.

RESULTS

Meta-analysis revealed significant replication of association between rs2114039 and corneal curvature (P = 0.0045). The SNP rs2114039 near PDGFRA has been previously implicated in Asians. No SNP at the FRAP1 locus was found to be associated in our Australian samples. No SNP surpassed the genomewide significance threshold of 5 × 10(-8). The SNP with strongest association was rs2444240 (P = 3.658 × 10(-7)), which is 31 kb upstream to the TRIM29 gene.

CONCLUSIONS

A significant role of the PDGFRA gene in determining corneal curvature in the Australian population was confirmed in this study, also highlighting the putative association of the TRIM29 locus with CC.

An overview of the genetic structure within the Italian population from genome-wide data

In spite of the common belief of Europe as reasonably homogeneous at genetic level, advances in high-throughput genotyping technology have resolved several gradients which define different geographical areas with good precision. When Northern and Southern European groups were considered separately, there were clear genetic distinctions. Intra-country genetic differences were also evident, especially in Finland and, to a lesser extent, within other European populations. Here, we present the first analysis using the 125,799 genome-wide Single Nucleotide Polymorphisms (SNPs) data of 1,014 Italians with wide geographical coverage. We showed by using Principal Component analysis and model-based individual ancestry analysis, that the current population of Sardinia can be clearly differentiated genetically from mainland Italy and Sicily, and that a certain degree of genetic differentiation is detectable within the current Italian peninsula population. Pair-wise F(ST) statistics Northern and Southern Italy amounts approximately to 0.001 between, and around 0.002 between Northern Italy and Utah residents with Northern and Western European ancestry (CEU). The Italian population also revealed a fine genetic substructure underscoring by the genomic inflation (Sardinia vs. Northern Italy = 3.040 and Northern Italy vs. CEU = 1.427), warning against confounding effects of hidden relatedness and population substructure in association studies.

Epigenetic variation during the adult lifespan: cross-sectional and longitudinal data on monozygotic twin pairs

The accumulation of epigenetic changes was proposed to contribute to the age-related increase in the risk of most common diseases. In this study on 230 monozygotic twin pairs (MZ pairs), aged 18–89 years, we investigated the occurrence of epigenetic changes over the adult lifespan. Using mass spectrometry, we investigated variation in global (LINE1) DNA methylation and in DNA methylation at INS, KCNQ1OT1, IGF2, GNASAS, ABCA1, LEP, and CRH, candidate loci for common diseases. Except for KCNQ1OT1, interindividual variation in locus-specific DNA methylation was larger in old individuals than in young individuals, ranging from 1.2-fold larger at ABCA1 (P = 0.010) to 1.6-fold larger at INS (P = 3.7 × 10⁻⁰⁷). Similarly, there was more within-MZ-pair discordance in old as compared with young MZ pairs, except for GNASAS, ranging from an 8% increase in discordance each decade at CRH (P = 8.9 × 10⁻⁰⁶) to a 16% increase each decade at LEP (P = 2.0 × 10⁻⁰⁸). Still, old MZ pairs with strikingly similar DNA methylation were also observed at these loci. After 10-year follow-up in elderly twins, the variation in DNA methylation showed a similar pattern of change as observed cross-sectionally. The age-related increase in methylation variation was generally attributable to unique environmental factors, except for CRH, for which familial factors may play a more important role. In conclusion, sustained epigenetic differences arise from early adulthood to old age and contribute to an increasing discordance of MZ twins during aging.

Association between in vivo alcohol metabolism and genetic variation in pathways that metabolize the carbon skeleton of ethanol and NADH reoxidation in the alcohol challenge twin study

BACKGROUND

Variation in alcohol metabolism affects the duration of intoxication and alcohol use. While the majority of genetic association studies investigating variation in alcohol metabolism have focused on polymorphisms in alcohol or aldehyde dehydrogenases, we have now tested for association with genes in alternative metabolic pathways that catalyze the carbon skeleton of ethanol (EtOH) and NADH reoxidation.

METHODS

Nine hundred fifty single nucleotide polymorphisms (SNPs) spanning 14 genes (ACN9, ACSS1, ACSS2, ALDH1A1, CAT, CYP2E1, GOT1, GOT2, MDH1, MDH2, SLC25A10, SLC25A11, SLC25A12, SLC25A13) were genotyped in 352 young adults who participated in an alcohol challenge study. Traits tested were blood alcohol concentration (BAC), breath alcohol concentration (BrAC), peak alcohol concentration, and rates of alcohol absorption and elimination. Allelic association was tested using quantitative univariate and multivariate methods.

RESULTS

A CYP2E1 promoter SNP (rs4838767, minor allele frequency 0.008) exceeded the threshold for study-wide significance (4.01 × 10(-5) ) for 2 early BAC, 8 BrAC measures, and the peak BrAC. For each phenotype, the minor C allele was related to a lower alcohol concentration, most strongly for the fourth BrAC (p = 2.07 × 10(-7) ) explaining ~8% of the phenotypic variance. We also observed suggestive patterns of association with variants in ALDH1A1 and on chromosome 17 near SLC25A11 for aspects of blood and breath alcohol metabolism. An SNP upstream of GOT1 (rs2490286) reached study-wide significance for multivariate BAC metabolism (p = 0.000040).

CONCLUSIONS

Overall, we did not find strong evidence that variation in genes coding for proteins that further metabolize the carbon backbone of acetaldehyde, or contribute to mechanisms for regenerating NAD from NADH, affects alcohol metabolism in our European-descent subjects. However, based on the breath alcohol data, variation in the promoter of CYP2E1 may play a role in preabsorptive or early hepatic alcohol metabolism, but more samples are required to validate this finding.

Cryptic distant relatives are common in both isolated and cosmopolitan genetic samples

Although a few hundred single nucleotide polymorphisms (SNPs) suffice to infer close familial relationships, high density genome-wide SNP data make possible the inference of more distant relationships such as 2^nd to 9^th cousinships. In order to characterize the relationship between genetic similarity and degree of kinship given a timeframe of 100–300 years, we analyzed the sharing of DNA inferred to be identical by descent (IBD) in a subset of individuals from the 23andMe customer database (n = 22,757) and from the Human Genome Diversity Panel (HGDP-CEPH, n = 952). With data from 121 populations, we show that the average amount of DNA shared IBD in most ethnolinguistically-defined populations, for example Native American groups, Finns and Ashkenazi Jews, differs from continentally-defined populations by several orders of magnitude. Via extensive pedigree-based simulations, we determined bounds for predicted degrees of relationship given the amount of genomic IBD sharing in both endogamous and ‘unrelated’ population samples. Using these bounds as a guide, we detected tens of thousands of 2^nd to 9^th degree cousin pairs within a heterogenous set of 5,000 Europeans. The ubiquity of distant relatives, detected via IBD segments, in both ethnolinguistic populations and in large ‘unrelated’ populations samples has important implications for genetic genealogy, forensics and genotype/phenotype mapping studies.

Five years of GWAS discovery

The past five years have seen many scientific and biological discoveries made through the experimental design of genome-wide association studies (GWASs). These studies were aimed at detecting variants at genomic loci that are associated with complex traits in the population and, in particular, at detecting associations between common single-nucleotide polymorphisms (SNPs) and common diseases such as heart disease, diabetes, auto-immune diseases, and psychiatric disorders. We start by giving a number of quotes from scientists and journalists about perceived problems with GWASs. We will then briefly give the history of GWASs and focus on the discoveries made through this experimental design, what those discoveries tell us and do not tell us about the genetics and biology of complex traits, and what immediate utility has come out of these studies. Rather than giving an exhaustive review of all reported findings for all diseases and other complex traits, we focus on the results for auto-immune diseases and metabolic diseases. We return to the perceived failure or disappointment about GWASs in the concluding section.

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.

Brazilian urban population genetic structure reveals a high degree of admixture

Advances in genotyping technologies have contributed to a better understanding of human population genetic structure and improved the analysis of association studies. To analyze patterns of human genetic variation in Brazil, we used SNP data from 1129 individuals--138 from the urban population of Sao Paulo, Brazil, and 991 from 11 populations of the HapMap Project. Principal components analysis was performed on the SNPs common to these populations, to identify the composition and the number of SNPs needed to capture the genetic variation of them. Both admixture and local ancestry inference were performed in individuals of the Brazilian sample. Individuals from the Brazilian sample fell between Europeans, Mexicans, and Africans. Brazilians are suggested to have the highest internal genetic variation of sampled populations. Our results indicate, as expected, that the Brazilian sample analyzed descend from Amerindians, African, and/or European ancestors, but intermarriage between individuals of different ethnic origin had an important role in generating the broad genetic variation observed in the present-day population. The data support the notion that the Brazilian population, due to its high degree of admixture, can provide a valuable resource for strategies aiming at using admixture as a tool for mapping complex traits in humans.

GWAS findings for human iris patterns: associations with variants in genes that influence normal neuronal pattern development

Human iris patterns are highly variable. The origins of this variation are of interest in the study of iris-related eye diseases and forensics, as well as from an embryological developmental perspective, with regard to their possible relationship to fundamental processes of neurodevelopment. We have performed genome-wide association scans on four iris characteristics (crypt frequency, furrow contractions, presence of peripupillary pigmented ring, and number of nevi) in three Australian samples of European descent. Both the discovery (n = 2121) and replication (n = 499 and 73) samples showed evidence for association between (1) crypt frequency and variants in the axonal guidance gene SEMA3A (p = 6.6 × 10(-11)), (2) furrow contractions and variants within the cytoskeleton gene TRAF3IP1 (p = 2.3 × 10(-12)), and (3) the pigmented ring and variants in the well-known pigmentation gene SLC24A4 (p = 7.6 × 10(-21)). These replicated findings individually accounted for around 1.5%-3% of the variance for these iris characteristics. Because both SEMA3A and TRAFIP1 are implicated in pathways that control neurogenesis, neural migration, and synaptogenesis, we also examined the evidence of enhancement among such genes, finding enrichment for crypts and furrows. These findings suggest that genes involved in normal neuronal pattern development may also influence tissue structures in the human iris.

The genetic structure of the Swedish population

Patterns of genetic diversity have previously been shown to mirror geography on a global scale and within continents and individual countries. Using genome-wide SNP data on 5174 Swedes with extensive geographical coverage, we analyzed the genetic structure of the Swedish population. We observed strong differences between the far northern counties and the remaining counties. The population of Dalarna county, in north middle Sweden, which borders southern Norway, also appears to differ markedly from other counties, possibly due to this county having more individuals with remote Finnish or Norwegian ancestry than other counties. An analysis of genetic differentiation (based on pairwise F_st) indicated that the population of Sweden's southernmost counties are genetically closer to the HapMap CEU samples of Northern European ancestry than to the populations of Sweden's northernmost counties. In a comparison of extended homozygous segments, we detected a clear divide between southern and northern Sweden with small differences between the southern counties and considerably more segments in northern Sweden. Both the increased degree of homozygosity in the north and the large genetic differences between the south and the north may have arisen due to a small population in the north and the vast geographical distances between towns and villages in the north, in contrast to the more densely settled southern parts of Sweden. Our findings have implications for future genome-wide association studies (GWAS) with respect to the matching of cases and controls and the need for within-county matching. We have shown that genetic differences within a single country may be substantial, even when viewed on a European scale. Thus, population stratification needs to be accounted for, even within a country like Sweden, which is often perceived to be relatively homogenous and a favourable resource for genetic mapping, otherwise inferences based on genetic data may lead to false conclusions.

Ethical concerns related to developing pharmacogenomic treatment strategies for addiction

Pharmacogenomics (PGx) research is poised to enable physicians to identify optimally effective treatments for individual substance abusers based on their genetic profiles. This paper addresses ethical issues related to PGx treatment strategies for addiction, focusing especially on the use of race variables in genomics research and ensuring equitable access to novel PGx treatments. Unless the field addresses the ethical challenges posed by these issues, PGx treatment innovations for addiction threaten to exacerbate already dramatic disparities in the burden of drug dependence for minority and other underserved populations.

Natural selection among Eurasians at genomic regions associated with HIV-1 control

Background

HIV susceptibility and pathogenicity exhibit both interindividual and intergroup variability. The etiology of intergroup variability is still poorly understood, and could be partly linked to genetic differences among racial/ethnic groups. These genetic differences may be traceable to different regimes of natural selection in the 60,000 years since the human radiation out of Africa. Here, we examine population differentiation and haplotype patterns at several loci identified through genome-wide association studies on HIV-1 control, as determined by viral-load setpoint, in European and African-American populations. We use genome-wide data from the Human Genome Diversity Project, consisting of 53 world-wide populations, to compare measures of F_ST and relative extended haplotype homozygosity (REHH) at these candidate loci to the rest of the respective chromosome.

Results

We find that the Europe-Middle East and Europe-South Asia pairwise F_ST in the most strongly associated region are elevated compared to most pairwise comparisons with the sub-Saharan African group, which exhibit very low F_ST. We also find genetic signatures of recent positive selection (higher REHH) at these associated regions among all groups except for sub-Saharan Africans and Native Americans. This pattern is consistent with one in which genetic differentiation, possibly due to diversifying/positive selection, occurred at these loci among Eurasians.

Conclusions

These findings are concordant with those from earlier studies suggesting recent evolutionary change at immunity-related genomic regions among Europeans, and shed light on the potential genetic and evolutionary origin of population differences in HIV-1 control.

Human population dispersal "Out of Africa" estimated from linkage disequilibrium and allele frequencies of SNPs

Genetic and fossil evidence supports a single, recent (<200,000 yr) origin of modern Homo sapiens in Africa, followed by later population divergence and dispersal across the globe (the "Out of Africa" model). However, there is less agreement on the exact nature of this migration event and dispersal of populations relative to one another. We use the empirically observed genetic correlation structure (or linkage disequilibrium) between 242,000 genome-wide single nucleotide polymorphisms (SNPs) in 17 global populations to reconstruct two key parameters of human evolution: effective population size (N(e)) and population divergence times (T). A linkage disequilibrium (LD)-based approach allows changes in human population size to be traced over time and reveals a substantial reduction in N(e) accompanying the "Out of Africa" exodus as well as the dramatic re-expansion of non-Africans as they spread across the globe. Secondly, two parallel estimates of population divergence times provide clear evidence of population dispersal patterns "Out of Africa" and subsequent dispersal of proto-European and proto-East Asian populations. Estimates of divergence times between European-African and East Asian-African populations are inconsistent with its simplest manifestation: a single dispersal from the continent followed by a split into Western and Eastern Eurasian branches. Rather, population divergence times are consistent with substantial ancient gene flow to the proto-European population after its divergence with proto-East Asians, suggesting distinct, early dispersals of modern H. sapiens from Africa. We use simulated genetic polymorphism data to demonstrate the validity of our conclusions against alternative population demographic scenarios.

Swedish population substructure revealed by genome-wide single nucleotide polymorphism data

The use of genome-wide single nucleotide polymorphism (SNP) data has recently proven useful in the study of human population structure. We have studied the internal genetic structure of the Swedish population using more than 350,000 SNPs from 1525 Swedes from all over the country genotyped on the Illumina HumanHap550 array. We have also compared them to 3212 worldwide reference samples, including Finns, northern Germans, British and Russians, based on the more than 29,000 SNPs that overlap between the Illumina and Affymetrix 250K Sty arrays. The Swedes--especially southern Swedes--were genetically close to the Germans and British, while their genetic distance to Finns was substantially longer. The overall structure within Sweden appeared clinal, and the substructure in the southern and middle parts was subtle. In contrast, the northern part of Sweden, Norrland, exhibited pronounced genetic differences both within the area and relative to the rest of the country. These distinctive genetic features of Norrland probably result mainly from isolation by distance and genetic drift caused by low population density. The internal structure within Sweden (F(ST) = 0.0005 between provinces) was stronger than that in many Central European populations, although smaller than what has been observed for instance in Finland; importantly, it is of the magnitude that may hamper association studies with a moderate number of markers if cases and controls are not properly matched geographically. Overall, our results underline the potential of genome-wide data in analyzing substructure in populations that might otherwise appear relatively homogeneous, such as the Swedes.

Consanguineous Marriage and Human Evolution

Mate choice among early human groups and in many historical populations was subject to both demographic and social constraints, ensuring that most unions were between couples who had coinherited substantial proportions of their genomes from common ancestors. Even in populations in which close consanguineous marriage was proscribed, community endogamy would have been sufficient to ensure high levels of homozygosity. Consanguineous marriage remains the choice of an estimated 10.4% of the global population, although there has been an overall decline in its popularity, especially in developed countries. Recent studies have indicated that the shift from consanguineous marriage to panmixia has been accompanied by a reduction in homozygosity. The concomitant predicted decrease in incidence of both recessive single-gene disorders and more common adult-onset diseases will have a significant impact on the health of future generations.

Fine-scale population structure and the era of next-generation sequencing

Fine-scale population structure characterizes most continents and is especially pronounced in non-cosmopolitan populations. Roughly half of the world's population remains non-cosmopolitan and even populations within cities often assort along ethnic and linguistic categories. Barriers to random mating can be ecologically extreme, such as the Sahara Desert, or cultural, such as the Indian caste system. In either case, subpopulations accumulate genetic differences if the barrier is maintained over multiple generations. Genome-wide polymorphism data, initially with only a few hundred autosomal microsatellites, have clearly established differences in allele frequency not only among continental regions, but also within continents and within countries. We review recent evidence from the analysis of genome-wide polymorphism data for genetic boundaries delineating human population structure and the main demographic and genomic processes shaping variation, and discuss the implications of population structure for the distribution and discovery of disease-causing genetic variants, in the light of the imminent availability of sequencing data for a multitude of diverse human genomes.

Sequencing and analysis of an Irish human genome

Background

Recent studies generating complete human sequences from Asian, African and European subgroups have revealed population-specific variation and disease susceptibility loci. Here, choosing a DNA sample from a population of interest due to its relative geographical isolation and genetic impact on further populations, we extend the above studies through the generation of 11-fold coverage of the first Irish human genome sequence.

Results

Using sequence data from a branch of the European ancestral tree as yet unsequenced, we identify variants that may be specific to this population. Through comparisons with HapMap and previous genetic association studies, we identified novel disease-associated variants, including a novel nonsense variant putatively associated with inflammatory bowel disease. We describe a novel method for improving SNP calling accuracy at low genome coverage using haplotype information. This analysis has implications for future re-sequencing studies and validates the imputation of Irish haplotypes using data from the current Human Genome Diversity Cell Line Panel (HGDP-CEPH). Finally, we identify gene duplication events as constituting significant targets of recent positive selection in the human lineage.

Conclusions

Our findings show that there remains utility in generating whole genome sequences to illustrate both general principles and reveal specific instances of human biology. With increasing access to low cost sequencing we would predict that even armed with the resources of a small research group a number of similar initiatives geared towards answering specific biological questions will emerge.

Inferring geographic coordinates of origin for Europeans using small panels of ancestry informative markers

Recent large-scale studies of European populations have demonstrated the existence of population genetic structure within Europe and the potential to accurately infer individual ancestry when information from hundreds of thousands of genetic markers is used. In fact, when genomewide genetic variation of European populations is projected down to a two-dimensional Principal Components Analysis plot, a surprising correlation with actual geographic coordinates of self-reported ancestry has been reported. This substructure can hamper the search of susceptibility genes for common complex disorders leading to spurious correlations. The identification of genetic markers that can correct for population stratification becomes therefore of paramount importance. Analyzing 1,200 individuals from 11 populations genotyped for more than 500,000 SNPs (Population Reference Sample), we present a systematic exploration of the extent to which geographic coordinates of origin within Europe can be predicted, with small panels of SNPs. Markers are selected to correlate with the top principal components of the dataset, as we have previously demonstrated. Performing thorough cross-validation experiments we show that it is indeed possible to predict individual ancestry within Europe down to a few hundred kilometers from actual individual origin, using information from carefully selected panels of 500 or 1,000 SNPs. Furthermore, we show that these panels can be used to correctly assign the HapMap Phase 3 European populations to their geographic origin. The SNPs that we propose can prove extremely useful in a variety of different settings, such as stratification correction or genetic ancestry testing, and the study of the history of European populations.