Genome Scan Studies Research Articles

Variation in recombination rate affects detection of outliers in genome scans under neutrality.

Genome scans can potentially identify genetic loci involved in evolutionary processes such as local adaptation and gene flow. Here, we show that recombination rate variation across a neutrally evolving genome gives rise to mixed sampling distributions of mean FST ( ), a common population genetic summary statistic. In particular, we show that in regions of low recombination the distribution of estimates has more variance and a longer tail than in more highly recombining regions. Determining outliers from the genome-wide distribution without taking local recombination rate into consideration may therefore increase the frequency of false positives in low recombination regions and be overly conservative in more highly recombining ones. We perform genome scans on simulated and empirical Drosophila melanogaster data sets and, in both cases, find patterns consistent with this neutral model. Similar patterns are observed for other summary statistics used to capture variation in the coalescent process. Linked selection, particularly background selection, is often invoked to explain heterogeneity in across the genome, but here we point out that even under neutrality, statistical artefacts can arise due to variation in recombination rate. Our results highlight a flaw in the design of genome-scan studies and suggest that without estimates of local recombination rate, interpreting the genomic landscape of any summary statistic that captures variation in the coalescent process will be very difficult.

Signatures of genetic adaptation to extremely varied Australian environments in introduced European house sparrows.

Due to its history of multiple introductions to novel environments worldwide, the house sparrow has been used as a model species to study local adaption in invasive avian species. New genomic resources such as a custom 200K SNP array and a house sparrow reference genome provide great prospects for studying rapid local adaptation in this invasive species. Here, we analyse high-density genomewide genetic data collected across an extensive range of temperate, arid and tropical climates, in Australian populations that were introduced from Europe 150years ago. We used two population differentiation (PD) and two ecological association (EA) methods to identify putative loci subject to selection across these varied climates. A majority of the outlier SNPs were identified through the use of the latent factor mixed models (LFMM) EA method, but the BayeScEnv EA method had the strongest overlap with the outliers from the two PD methods. Out of all the 971 outliers identified across the different methods, 38.3% were physically linked (within 20kbps) to 575 known protein-coding regions in the house sparrow reference genome. Interestingly, some outlier genes had been previously identified in genome scan studies of broadly distributed species or had strong links to traits that are expected to be important to local adaptation, for example, heat-shock proteins, immune response and HOX genes. However, many outliers still have unknown relevance and some outliers can be false positives. Our results identify an opportunity to use the house sparrow model to further study local adaptation in an invasive species.

Assessment of genetic diversity and phylogenetic relationship of Limousin herds in Hungary using microsatellite markers.

ObjectiveThis study was conducted to investigate basic information on genetic structure and characteristics of Limousin population in Hungary. Obtained results will be taken into consideration when adopting the new breeding strategy by the Association of Hungarian Limousin and Blonde d’Aquitaine Breeders (AHLBB).MethodsGenetic diversity and phylogenetic relationship of 3,443 Limousin cattle from 16 different herds were investigated by performing genotyping using 18 microsatellite markers. Amplified DNA was genotyped using an automated genetic analyzer.ResultsMean of effective alleles (ne) of the populations was 3.77. Population C had the lowest number of effective alleles (3.01) and the lowest inbreeding coefficient (FIS) value (−0.15). Principal component analysis of estimated genetic distance (FST) values (p<0.000) revealed two herds (C and E) distinct from the majority of other Limousin herds. The pairwise FST values of population C compared to the others (0.066 to 0.120) fell into the range of moderate genetic distance: 0.050 to 0.150, while population E displayed also moderate genetic distance (FST values in range 0.052 to 0.064) but only to six populations (G, H, J, L, N, and P). FST(C-E) was 0.148, all other pairs -excluding C and E herds- displayed low genetic distance (FST<0.049). Population D, F, I, J, K, L, N, O, and P carried private alleles, which alleles belonged to 1.1% of the individuals. Most probable number of clusters (K) were 2 and 7 determined by Structure and BAPS software.ConclusionThis study showed useful genetic diversity and phylogenetic relationship data that can be utilized for the development of a new breeding strategy by AHLBB. The results presented could also contribute to the proper selection of animals for further whole genome scan studies of Limousins.

HLA-DR7 and HLA-DQ2: Transgenic mouse strains tested as a model system for ximelagatran hepatotoxicity.

The oral thrombin inhibitor ximelagatran was withdrawn in the late clinical trial phase because it adversely affected the liver. In approximately 8% of treated patients, drug-induced liver injury (DILI) was expressed as transient alanine transaminase (ALT) elevations. No evidence of DILI had been revealed in the pre-clinical in vivo studies. A whole genome scan study performed on the clinical study material identified a strong genetic association between the major histocompatibility complex alleles for human leucocyte antigens (HLA) (HLA-DR7 and HLA-DQ2) and elevated ALT levels in treated patients. An immune-mediated pathogenesis was suggested. Here, we evaluated whether HLA transgenic mice models could be used to investigate whether the expression of relevant HLA molecules was enough to reproduce the DILI effects in humans. In silico modelling performed in this study revealed association of both ximelagatran (pro-drug) and melagatran (active drug) to the antigen-presenting groove of the homology modelled HLA-DR7 molecule suggesting “altered repertoire” as a key initiating event driving development of DILI in humans. Transgenic mouse strains (tgms) expressing HLA of serotype HLA-DR7 (HLA-DRB1*0701, -DRA*0102), and HLA-DQ2 (HLA-DQB1*0202,–DQA1*0201) were created. These two lines were crossed with a human (h)CD4 transgenic line, generating the two tgms DR7xhCD4 and DQ2xhCD4. To investigate whether the DILI effects observed in humans could be reproduced in tgms, the mice were treated for 28 days with ximelagatran. Results revealed no signs of DILI when biomarkers for liver toxicity were measured and histopathology was evaluated. In the ximelagatran case, presence of relevant HLA-expression in a pre-clinical model did not fulfil the prerequisite for reproducing DILI observed in patients. Nonetheless, for the first time an HLA-transgenic mouse model has been investigated for use in HLA-associated DILI induced by a low molecular weight compound. This study shows that mimicking of genetic susceptibility, expressed as DILI-associated HLA-types in mice, is not sufficient for reproducing the complex pathogenesis leading to DILI in man.

Production of a Novel Amide-Containing Polyene by Activating a Cryptic Biosynthetic Gene Cluster in Streptomyces sp. MSC090213JE08.

Streptomyces sp. MSC090213JE08 seems to have more than 20 cryptic biosynthetic gene clusters for secondary metabolites. We aimed to activate some of them by forced production of Streptomyces antibiotic regulatory protein (SARP) family transcriptional activators. We constructed seven recombinant strains, each of which contained a SARP gene under the control of a constitutive promoter, and subjected them to comparative metabolic profiling analysis. Four of the seven strains produced nine metabolites that were hardly detected in the control strains. We isolated a new metabolite (named ishigamide) from the SARP-7-expressing strain and determined its structure as 3-((2E,4E,6E,8E)-13-hydroxytetradeca-2,4,6,8-tetraenamido)propanoic acid. Genome scanning and gene disruption studies identified the ishigamide biosynthetic gene cluster adjacent to the SARP-7 gene. We think that a new subfamily of type II polyketide synthase is involved in the biosynthesis of the polyene structure of ishigamide.

Genetic structure of colline and montane populations of an endangered plant species.

Due to land-use intensification, lowland and colline populations of many plants of nutrient-poor grasslands have been strongly fragmented in the last decades, with potentially negative consequences for their genetic diversity and persistence. Populations in mountains might represent a genetic reservoir for grassland plants, because they have been less affected by land-use changes. We studied the genetic structure and diversity of colline and montane Vosges populations of the threatened perennial plant Arnica montana in western central Europe using AFLP markers. Our results indicate that in contrast to our expectation even strongly fragmented colline populations of A. montana have conserved a considerable amount of genetic diversity. However, mean seed mass increased with the proportion of polymorphic loci, suggesting inbreeding effects in low diversity populations. At a similar small geographical scale there was a clear IBD pattern for the montane Vosges but not for the colline populations. However, there was a strong IBD-pattern for the colline populations at a large geographical scale suggesting that this pattern is a legacy of historical gene flow, as most of the colline populations are today strongly isolated from each other. Genetic differentiation between colline and montane Vosges populations was strong. Moreover, results of a genome scan study indicated differences in loci under selection, suggesting that plants from montane Vosges populations might be maladapted to conditions at colline sites. Our results suggest caution in using material from montane populations of rare plants for the reinforcement of small genetically depauperate lowland populations.

Association Study of GABRG2 Polymorphisms with Suicidal Behaviour in Schizophrenia Patients with Alcohol Use Disorder

Background: Schizophrenia is a severe neuropsychiatric disorder where the role of γ-aminobutyric acid (GABA), an inhibitory neurotransmitter, has been implicated in its aetiopathophysiology. Several genes coding for GABA_A subunits, including the GABRG2 gene that encodes the γ₂ subunit, are clustered at 5q31-q35, a chromosomal region that is associated with schizophrenia in genome scan studies. We recently reported GABRG2 to be associated with schizophrenia in our case-control and family samples. Methods: We tested eight single-nucleotide polymorphisms spanning the GABRG2 gene for an association with suicidal behaviour in our schizophrenia sample of European ancestry (n = 197), taking into account history of alcohol abuse or dependence. Results: We found the haplotypes of the rs183294 and rs209356 markers to be significantly associated with history of suicide attempt (p < 0.01) as well as suicide specifier scores (p < 0.05). The association appeared to be originating in patients with a history of alcohol dependence or abuse. Conclusions: Taken together, the results of the present study suggest that GABRG2 may be involved in suicidal behaviour in schizophrenia patients with alcohol dependence or abuse, but replications are required. These results may help in the discovery of novel treatments for alcoholism and/or prevention of suicide.

Genome scan study of prostate cancer in Arabs: identification of three genomic regions with multiple prostate cancer susceptibility loci in Tunisians

BackgroundLarge databases focused on genetic susceptibility to prostate cancer have been accumulated from population studies of different ancestries, including Europeans and African-Americans. Arab populations, however, have been only rarely studied.MethodsUsing Affymetrix Genome-Wide Human SNP Array 6, we conducted a genome-wide association study (GWAS) in which 534,781 single nucleotide polymorphisms (SNPs) were genotyped in 221 Tunisians (90 prostate cancer patients and 131 age-matched healthy controls). TaqMan® SNP Genotyping Assays on 11 prostate cancer associated SNPs were performed in a distinct cohort of 337 individuals from Arab ancestry living in Qatar and Saudi Arabia (155 prostate cancer patients and 182 age-matched controls). In-silico expression quantitative trait locus (eQTL) analysis along with mRNA quantification of nearby genes was performed to identify loci potentially cis-regulated by the identified SNPs.ResultsThree chromosomal regions, encompassing 14 SNPs, are significantly associated with prostate cancer risk in the Tunisian population (P = 1 × 10-4 to P = 1 × 10-5). In addition to SNPs located on chromosome 17q21, previously found associated with prostate cancer in Western populations, two novel chromosomal regions are revealed on chromosome 9p24 and 22q13. eQTL analysis and mRNA quantification indicate that the prostate cancer associated SNPs of chromosome 17 could enhance the expression of STAT5B gene.ConclusionOur findings, identifying novel GWAS prostate cancer susceptibility loci, indicate that prostate cancer genetic risk factors could be ethnic specific.

Genetic epidemiology of cleft lip and palate: A review and potentials for Africa

Orofacial clefts (CLP) are the most common congenital malformation of the craniofacial region with a world-wide prevalence of 1 in 700 births. They are usually classified into syndromic and non syndromic clefts and can be amongst the phenotypes seen in infants with multiple congenital malformations that can not be classified as a syndrome. While the genetic etiologies of syndromic clefts have been well characterized, the genetic etiology of non-syndromic clefts remains a challenge. However, with advances in the molecular genetic technologies and biostatistics tools that includes linkage and association studies, whole genome scanning and genome-wide association studies, candidate genes and loci in chromosomes with significant associations have been identified and validated. These include loci on chromosome 8q.24 and genes such as IRF6, MAFB and VAX1. Interestingly, all these genes were identified using populations of European and Asian ancestry. Considering that African populations are ancestral to European and Asian populations, genetic studies using predominantly African populations in Africa could provide particularly powerful insights into the origin of cleft lip and palate with the prediction that they would exhibit their own unique characteristics. The increasing pace of global genetic knowledge, and its applications for clinical medicine, as well as known geographic and genetic heterogeneity in non- syndromic cleft lip and palate (NSCLP) provides good scientific rational for the identification of new variants in African populations living in Africa.

Malaria parasite genome scan: insights into antimalarial resistance

Genome scan and genotype-phenotype association study offer excellent opportunities to unearth drug/vaccine targets in human pathogens including malaria parasites. A recently conducted such study in worldwide isolates in the most devastating malaria parasite Plasmodium falciparum has reported important genomic information on genetic basis of antimalarial resistance. Several unknown genes were also found to be under strong influence of natural selection. The findings provide important insights into the malaria parasite genome evolution in general and to use this information to develop more focused malaria control strategies, in particular.

Whole‐genome scan for quantitative trait loci associated with birth weight, gestation length and passive immune transfer in a Holstein × Jersey crossbred population

We herein report results from a daughter design genome-scan study aiming to identify quantitative trait loci (QTL) associated with birth weight, direct gestation length and passive immune transfer in a backcross (Holstein x Jersey) x Holstein population. Two-hundred and seventy-six calves, offspring of seven crossbred sires, were genotyped for 161 microsatellite markers distributed along the 29 bovine autosomes. The genome scan was performed through interval mapping using an animal model in order to identify QTL accounting for phenotypic differences between individual animals. Based on significant chi-squared values, we identified putative QTL on BTA7 and BTA14 for gestation length, on BTA2, BTA6 and BTA14 for birth weight and on BTA20 for passive immune transfer. In total, these QTL accounted for 12%, 18% and 1% of the phenotypic variance in gestation length, birth weight and passive immune transfer respectively. We also report results from a supplementary and independent influential grand-daughter Holstein family. In this family, findings on BTA7 and BTA14 for direct gestation length were in agreement with results in the crossbred population. Two other regions on BTA6 and BTA21 putatively underlying QTL for direct gestation length variability were discovered with this analysis.

Association of a polymorphism in the NRXN3 gene with the degree of smoking in schizophrenia: A preliminary study

Whole genome scan studies have recently identified the NRXN1 and NRXN3 genes as potential contributing factors in the risk for nicotine addiction. We have genotyped 15 single nucleotide polymorphisms (SNPs) spanning the NRXN1 and NRXN3 genes in 195 unrelated patients with schizophrenia for whom information about their smoking status and number of cigarettes smoked per day (CPD) was obtained. The NRXN3 marker rs1004212 was significantly associated with quantity of tobacco smoked. Individuals homozygous for the C allele of rs1004212 smoked more cigarettes per day than heterozygous individuals. We found no significant association of markers within the NRXN1 gene with the risk of smoking or the quantity of tobacco smoked. Because of the relatively small sample size, this is a preliminary study. However, this candidate gene study supports the observations of molecular studies implicating the NRXN genes in drug addiction and suggests that variants in the NRXN3 gene could contribute to the degree of nicotine dependence in patients with schizophrenia.

The genetics of intervertebral disc degeneration.Associated genes

Objectives To review current knowledge on genes associated with intervertebral disk degeneration. Methods A literature-based narrative review of the English language medical literature. Results and conclusions There are a number of genes that have been associated with intervertebral disk degeneration in humans, including genes coding for collagen I, collagen IX (COL9A2 and COL9A3), collagen XI (COL11A2), IL-1, aggrecan, vitamin D receptor, MMP-3, and CILP. For specific genes and some environmental factors, gene–gene, gene–environment and gene–age interactions may exist. Candidate–gene association studies have limitations in detecting the genetic basis of the disease because this approach relies on having predicted the correct genes on the basis of a biological hypothesis or the location of known linkage regions. Additional studies, including linkage analyses and whole genome scan studies in different populations and whole range of ages, are required to improve our understanding of the influence of the aforementioned genes on intervertebral disk degeneration and identify novel genes.

Constructing gene association networks for rheumatoid arthritis using the backward genotype-trait association (BGTA) algorithm

Rheumatoid arthritis (RA, MIM 180300) is a common and complex inflammatory disorder. The North American Rheumatoid Arthritis Consortium (NARAC) data, as part of the Genetic Analysis Workshop 15 data, consists of both genome scan and candidate gene studies on RA patients. We applied the backward genotype-trait association (BGTA) algorithm to capture marginal and gene x gene interaction effects of multiple susceptibility loci on RA disease status. A two-stage screening approach was used for the genome scan, whereas a comprehensive study of all possible subsets was conducted for the candidate genes. For the genome scan, we constructed an association network among 39 genetic loci that demonstrated strong signals, 19 of which have been reported in the RA literature. For the candidate genes, we found strong signals for PTPN22 and SUMO4. Based on significant association evidence, we built an association network among the loci of PTPN22, PADI4, DLG5, SLC22A4, SUMO4, and CARD15. To control for false positives, we used permutation tests to constrain the family-wise type I error rate to 1%. Using the BGTA algorithm, we identified genetic loci and candidate genes that were associated with RA susceptibility and association networks among them. For the first time, we report possible interactions between single-nucleotide polymorphisms/genes, which may be useful for biological interpretation.

Evidence for induction of the ornithine transcarbamylase expression in Alzheimer's disease

To more rapidly identify candidate genes located within chromosomal regions of interest defined by genome scan studies in Alzheimer's disease (AD), we have developed a customized microarray containing all the ORFs (n=2741) located within nine of these regions. Levels of gene expression were assessed in total RNA from brain tissue of 12 controls and 12 AD patients. Of all genes showing differential expression, we focused on the ornithine transcarbamylase (OTC) gene on Xp21.1., a key enzyme of the urea cycle which we found to be expressed in AD brains but not in controls, as confirmed by RT-PCR. We also detected mRNA expression of all the other urea cycle enzymes in AD brains. Immunochemistry experiments revealed that the OTC expression was strictly restricted to vascular endothelial cells in brain. Furthermore, OTC activity was 880% increased in the CSF of probable AD cases compared with controls. We analysed the association of the OTC -389 G/A and -241 A/G promoter polymorphisms with the risk of developing AD. We observed that rare haplotypes may be associated with the risk of AD through a possible modulation of the methylation of the OTC promoter. In conclusion, our results suggest the involvement of a new pathway in AD brains involving the urea cycle.

Glutamate receptor 6 gene (GluR6 or GRIK2) polymorphisms in the Indian population: a genetic association study on autism spectrum disorder.

Autism is a neurodevelopmental disorder with early manifestation. It is a multifactorial disorder and several susceptible chromosomal regions for autism are identified through genome scan studies. The gene coding for glutamate receptor 6 (GluR6 or GRIK2) has been suggested as a candidate gene for autism based on its localization in the autism specific region on chromosome 6q21 and the involvement of receptor protein in cognitive functions like learning and memory. Despite its importance, so far no studies have been carried out on possible involvement of GluR6 with autism in the Indian population. Therefore in the present study, we have performed genetic analysis of three markers of GluR6 (SNP1: rs2227281, SNP2: rs2227283, SNP3: rs2235076) for possible association with autism through population, and family-based (TDT and HHRR) approaches. DSM-IV criteria and CARS/ADI-R have been utilized for diagnosis. Genotyping analysis for the SNPs has been carried out in 101 probands with autism spectrum disorder, 180 parents and 152 controls from different regions of India. Since the minor allele frequency of SNP3 was too low, the association studies have been carried out only for SNP1 and SNP2. Even though two earlier studies have shown association of these markers with autism, the present case-control and TDT, as well as HHRR analyses have not demonstrated any biased transmission of alleles or haplotypes to the affected offspring. Thus our results suggest that these markers of GluR6 are unlikely to be associated with autism in the Indian population.

Transmission disequilibrium study of an oligodendrocyte and myelin glycoprotein gene allele in 431 families with an autistic proband

Autistic disorder is a neurodevelopmental disorder where genetic factors play an important role. We previously described an association between a subgroup of French autistic patients and an allele of a non-synonymous single nucleotide polymorphism (nsSNP: OMGP62 G > A or rs11080149) in the gene coding for the oligodendrocyte and myelin glycoprotein (OMG), located at 7 Mb from the marker D17S250, linked to autism in two independent genome scan studies. We report a study on 431 families with 1 affected child from different origins: French Canada ( n = 262), Italy ( n = 123) and United States ( n = 46). We analyzed the transmission of the rs11080149 alleles from parents to their affected children. There was a preferential transmission of the G allele from parents to affected children ( p = 0.0017) in the overall sample. Paternal and maternal transmission rates were both skewed. Taking into account our previous results obtained in a French group of patients, where we observed an association with allele A, a direct role of this polymorphism is improbable in autism. The associations observed in Japanese and French patients, the linkage studies and the present work speak in favor of the existence of a susceptibility gene for autism in the NF1 locus.

Inroads to Mechanisms of Disease in Child Psychiatric Disorders

Inroads to Mechanisms of Disease in Child Psychiatric Disorders

A functional Fas promoter polymorphism is associated with a severe phenotype in type 1 autoimmune hepatitis characterized by early development of cirrhosis

Genome scanning studies suggest an important role for genes outside the major histocompatibility complex in autoimmunity. Key candidates are those genes involved in immune regulation and preservation of immune homeostasis, including the genes involved in apoptosis. Our aim was to determine the association between the Fas gene polymorphism at position -670 and susceptibility, clinical expression, and outcome in type 1 autoimmune hepatitis (AIH). An adenosine to guanine single nucleotide polymorphism in the Fas gene (TNFRSF6) promoter was assessed in 149 well-characterized Caucasoid patients and 172 matched controls. Patients and normal subjects had the similar TNFRSF6-670 allele and genotype frequencies. Serum aspartate aminotransferase (510 +/- 77 vs 283 +/- 53 U/l), gamma-globulin (3.3 +/- 0.2 vs 2.6 +/- 0.2 g/dl), and immunoglobulin G (2976 +/- 223 vs 2324 +/- 203 mg/dl) levels were higher in patients with the guanine/guanine genotype than in those with the adenosine/adenosine genotype. Cirrhosis at presentation was more common in patients with the adenosine/adenosine or adenosine/guanine genotypes than in those with the guanine/guanine genotype (29% vs 6%). Polymorphism of the Fas gene at position -670 does not influence susceptibility to AIH, but may affect the early development of cirrhosis.

Controlling the false-positive rate in multilocus genome scans for selection.

Rapid typing of genetic variation at many regions of the genome is an efficient way to survey variability in natural populations in an effort to identify segments of the genome that have experienced recent natural selection. Following such a genome scan, individual regions may be chosen for further sequencing and a more detailed analysis of patterns of variability, often to perform a parametric test for selection and to estimate the strength of a recent selective sweep. We show here that not accounting for the ascertainment of loci in such analyses leads to false inference of natural selection when the true model is selective neutrality, because the procedure of choosing unusual loci (in comparison to the rest of the genome-scan data) selects regions of the genome with genealogies similar to those expected under models of recent directional selection. We describe a simple and efficient correction for this ascertainment bias, which restores the false-positive rate to near-nominal levels. For the parameters considered here, we find that obtaining a test with the expected distribution of P-values depends on accurately accounting both for ascertainment of regions and for demography. Finally, we use simulations to explore the utility of relying on outlier loci to detect recent selective sweeps. We find that measures of diversity and of population differentiation are more effective than summaries of the site-frequency spectrum and that sequencing larger regions (2.5 kbp) in genome-scan studies leads to more power to detect recent selective sweeps.