The power of allele frequency comparisons to detect the footprint of selection in natural and experimental situations

Recently, inter-population comparisons of allele frequencies to detect past selection haven gained popularity. Data from genome-wide scans are used to detect the number and position of genes that have responded to unknown selection pressures in natural populations, or known selection pressures in experimental lines. Yet, the limitations and possibilities of these methods have not been well studied. In this paper, the objectives were (1) to investigate the distance over which a signal of directional selection is detectable under various scenarios, and (2) to study the power of the method depending on the properties of the used markers, for both natural populations and experimental set-ups. A combination of recurrence equations and simulations was used. The results show that intermediate strength selection on new mutations can be detected with a marker spacing of about 0.5 cM in large natural populations, 200 to 400 generations after the divergence of subpopulations. In experimental situations, only strong selection will be detectable, while markers can be spaced a few cM apart. Adaptation from standing variation in the base population will be hard to detect, though some solutions are presented for experimental designs.     

Variants in Neuropeptide Y Receptor 1 and 5 Are Associated with Nutrient-Specific Food Intake and Are Under Recent Selection in Europeans

There is a large variation in caloric intake and macronutrient preference between individuals and between ethnic groups, and these food intake patterns show a strong heritability. The transition to new food sources during the agriculture revolution around 11,000 years ago probably created selective pressure and shaped the genome of modern humans. One major player in energy homeostasis is the appetite-stimulating hormone neuropeptide Y, in which the stimulatory capacity may be mediated by the neuropeptide Y receptors 1, 2 and 5 (NPY1R, NPY2R and NPY5R). We assess association between variants in the NPY1R, NPY2R and NPY5R genes and nutrient intake in a cross-sectional, single-center study of 400 men aged 40 to 80 years, and we examine whether genomic regions containing these genes show signatures of recent selection in 270 HapMap individuals (90 Africans, 90 Asians, and 90 Caucasians) and in 846 Dutch bloodbank controls. Our results show that derived alleles in NPY1R and NPY5R are associated with lower carbohydrate intake, mainly because of a lower consumption of mono- and disaccharides. We also show that carriers of these derived alleles, on average, consume meals with a lower glycemic index and glycemic load and have higher alcohol consumption. One of these variants shows the hallmark of recent selection in Europe. Our data suggest that lower carbohydrate intake, consuming meals with a low glycemic index and glycemic load, and/or higher alcohol consumption, gave a survival advantage in Europeans since the agricultural revolution. This advantage could lie in overall health benefits, because lower carbohydrate intake, consuming meals with a low GI and GL, and/or higher alcohol consumption, are known to be associated with a lower risk of chronic diseases.     

Exome-Wide Association Analysis of Coronary Artery Disease in the Kingdom of Saudi Arabia Population

Coronary Artery Disease (CAD) remains the leading cause of mortality worldwide. Mortality rates associated with CAD have shown an exceptional increase particularly in fast developing economies like the Kingdom of Saudi Arabia (KSA). Over the past twenty years, CAD has become the leading cause of death in KSA and has reached epidemic proportions. This rise is undoubtedly caused by fast urbanization that is associated with a life-style that promotes CAD. However, the question remains whether genetics play a significant role and whether genetic susceptibility is increased in KSA compared to the well-studied Western European populations. Therefore, we performed an Exome-wide association study (EWAS) in 832 patients and 1,076 controls of Saudi Arabian origin to test whether population specific, strong genetic risk factors for CAD exist, or whether the polygenic risk score for known genetic risk factors for CAD, lipids, and Type 2 Diabetes show evidence for an enriched genetic burden. Our results do not show significant associations for a single genetic locus. However, the heritability estimate for CAD for this population was high (h² = 0.53, S.E. = 0.1, p = 4e^-12) and we observed a significant association of the polygenic risk score for CAD that demonstrates that the population of KSA, at least in part, shares the genetic risk associated to CAD in Western populations.     

Carbon and nitrogen in soil and vegetation at sites differing in successional age

We studied vegetation and soil development during primary succession in an inland drift sand area in the Netherlands. We compared five sites at which primary succession had started at different moments in the past, respectively 0, 10, 43 and 121 years ago, and a site at which succession had not yet started. In the three younger sites the vegetation was herbaceous, whereas in the two older sites a pine forest had formed. Forest formation was accompanied by the development of an FH-layer in the soil, an increase in the amount of soil organic matter, and an increase in nitrogen mineralisation rate from 1.9 to 18 g N m^–2 yr^–1. Soil moisture content also increased, whereas pH showed a steady decrease with site age. The vegetation changed from a herbaceous vegetation dominated by mosses and lichens and the grass species Corynephorus canescens and Festuca ovina towards a pine forest with an understorey vegetation dominated by Deschampsia flexuosa and, at the oldest site, with dwarf shrubs Empetrum nigrum and Vaccinium myrtillus. At the same time the total amounts of carbon and nitrogen of the ecosystem increased, with a relatively stronger increase of the carbon pool. The establishment of trees during succession greatly affects the dynamics of the ecosystem, especially its carbon dynamics.     

Genes in the ureteric budding pathway: association study on vesico-ureteral reflux patients

Vesico-ureteral reflux (VUR) is the retrograde passage of urine from the bladder to the urinary tract and causes 8.5% of end-stage renal disease in children. It is a complex genetic developmental disorder, in which ectopic embryonal ureteric budding is implicated in the pathogenesis. VUR is part of the spectrum of Congenital Anomalies of the Kidney and Urinary Tract (CAKUT). We performed an extensive association study for primary VUR using a two-stage, case-control design, investigating 44 candidate genes in the ureteric budding pathway in 409 Dutch VUR patients. The 44 genes were selected from the literature and a set of 567 single nucleotide polymorphisms (SNPs) capturing their genetic variation was genotyped in 207 cases and 554 controls. The 14 SNPs with p<0.005 were included in a follow-up study in 202 cases and 892 controls. Of the total cohort, ~50% showed a clear-cut primary VUR phenotype and ~25% had both a duplex collecting system and VUR. We also looked for association in these two extreme phenotype groups. None of the SNPs reached a significant p-value. Common genetic variants in four genes (GREM1, EYA1, ROBO2 and UPK3A) show a trend towards association with the development of primary VUR (GREM1, EYA1, ROBO2) or duplex collecting system (EYA1 and UPK3A). SNPs in three genes (TGFB1, GNB3 and VEGFA) have been shown to be associated with VUR in other populations. Only the result of rs1800469 in TGFB1 hinted at association in our study. This is the first extensive study of common variants in the genes of the ureteric budding pathway and the genetic susceptibility to primary VUR.     

Detection, imputation, and association analysis of small deletions and null alleles on oligonucleotide arrays

Copy-number variation (CNV) is a major contributor to human genetic variation. Recently, CNV associations with human disease have been reported. Many genome-wide association (GWA) studies in complex diseases have been performed with sets of biallelic single-nucleotide polymorphisms (SNPs), but the available CNV methods are still limited. We present a new method (TriTyper) that can infer genotypes in case-control data sets for deletion CNVs, or SNPs with an extra, untyped allele at a high-resolution single SNP level. By accounting for linkage disequilibrium (LD), as well as intensity data, calling accuracy is improved. Analysis of 3102 unrelated individuals with European descent, genotyped with Illumina Infinium BeadChips, resulted in the identification of 1880 SNPs with a common untyped allele, and these SNPs are in strong LD with neighboring biallelic SNPs. Simulations indicate our method has superior power to detect associations compared to biallelic SNPs that are in LD with these SNPs, yet without increasing type I errors, as shown in a GWA analysis in celiac disease. Genotypes for 1204 triallelic SNPs could be fully imputed, with only biallelic-genotype calls, permitting association analysis of these SNPs in many published data sets. We estimate that 682 of the 1655 unique loci reflect deletions; this is on average 99 deletions per individual, four times greater than those detected by other methods. Whereas the identified loci are strongly enriched for known deletions, 61% have not been reported before. Genes overlapping with these loci more often have paralogs (p = 0.006) and biologically interact with fewer genes than expected (p = 0.004).     

Association Analysis of Genetic Variants in the Myosin IXB Gene in Acute Pancreatitis

Introduction

Impairment of the mucosal barrier plays an important role in the pathophysiology of acute pancreatitis. The myosin IXB (MYO9B) gene and the two tight-junction adaptor genes, PARD3 and MAGI2, have been linked to gastrointestinal permeability. Common variants of these genes are associated with celiac disease and inflammatory bowel disease, two other conditions in which intestinal permeability plays a role. We investigated genetic variation in MYO9B, PARD3 and MAGI2 for association with acute pancreatitis.

Methods

Five single nucleotide polymorphisms (SNPs) in MYO9B, two SNPs in PARD3, and three SNPs in MAGI2 were studied in a Dutch cohort of 387 patients with acute pancreatitis and over 800 controls, and in a German cohort of 235 patients and 250 controls.

Results

Association to MYO9B and PARD3 was observed in the Dutch cohort, but only one SNP in MYO9B and one in MAGI2 showed association in the German cohort (p < 0.05). Joint analysis of the combined cohorts showed that, after correcting for multiple testing, only two SNPs in MYO9B remained associated (rs7259292, p = 0.0031, odds ratio (OR) 1.94, 95% confidence interval (95% CI) 1.35-2.78; rs1545620, p = 0.0006, OR 1.33, 95% CI 1.16-1.53). SNP rs1545620 is a non-synonymous SNP previously suspected to impact on ulcerative colitis. None of the SNPs showed association to disease severity or etiology.

Conclusion

Variants in MYO9B may be involved in acute pancreatitis, but we found no evidence for involvement of PARD3 or MAGI2.