A genome-wide scan study identifies a single nucleotide substitution in ASIP associated with white versus non-white coat-colour variation in sheep (Ovis aries)

In sheep, coat colour (and pattern) is one of the important traits of great biological, economic and social importance. However, the genetics of sheep coat colour has not yet been fully clarified. We conducted a genome-wide association study of sheep coat colours by genotyping 47 303 single-nucleotide polymorphisms (SNPs) in the Finnsheep population in Finland. We identified 35 SNPs associated with all the coat colours studied, which cover genomic regions encompassing three known pigmentation genes (TYRP1, ASIP and MITF) in sheep. Eighteen of these associations were confirmed in further tests between white versus non-white individuals, but none of the 35 associations were significant in the analysis of only non-white colours. Across the tests, the s66432.1 in ASIP showed significant association (P=4.2 × 10⁻¹¹ for all the colours; P=2.3 × 10⁻¹¹ for white versus non-white colours) with the variation in coat colours and strong linkage disequilibrium with other significant variants surrounding the ASIP gene. The signals detected around the ASIP gene were explained by differences in white versus non-white alleles. Further, a genome scan for selection for white coat pigmentation identified a strong and striking selection signal spanning ASIP. Our study identified the main candidate gene for the coat colour variation between white and non-white as ASIP, an autosomal gene that has been directly implicated in the pathway regulating melanogenesis. Together with ASIP, the two other newly identified genes (TYRP1 and MITF) in the Finnsheep, bordering associated SNPs, represent a new resource for enriching sheep coat-colour genetics and breeding.     

Comparative population genetics of a mimicry locus among hybridizing Heliconius butterfly species

The comimetic Heliconius butterfly species pair, H. erato and H. melpomene, appear to use a conserved Mendelian switch locus to generate their matching red wing patterns. Here we investigate whether H. cydno and H. pachinus, species closely related to H. melpomene, use this same switch locus to generate their highly divergent red and brown color pattern elements. Using an F2 intercross between H. cydno and H. pachinus, we first map the genomic positions of two novel red/brown wing pattern elements; the G locus, which controls the presence of red vs brown at the base of the ventral wings, and the Br locus, which controls the presence vs absence of a brown oval pattern on the ventral hind wing. The results reveal that the G locus is tightly linked to markers in the genomic interval that controls red wing pattern elements of H. erato and H. melpomene. Br is on the same linkage group but approximately 26 cM away. Next, we analyze fine-scale patterns of genetic differentiation and linkage disequilibrium throughout the G locus candidate interval in H. cydno, H. pachinus and H. melpomene, and find evidence for elevated differentiation between H. cydno and H. pachinus, but no localized signature of association. Overall, these results indicate that the G locus maps to the same interval as the locus controlling red patterning in H. melpomene and H. erato. This, in turn, suggests that the genes controlling red pattern elements may be homologous across Heliconius, supporting the hypothesis that Heliconius butterflies use a limited suite of conserved genetic switch loci to generate both convergent and divergent wing patterns.     

Genetic association of single nucleotide polymorphisms in dystrobrevin binding protein 1 gene with schizophrenia in a Malaysian population

Dystrobrevin binding protein 1 (DTNBP1) gene is pivotal in regulating the glutamatergic system. Genetic variants of the DTNBP1 affect cognition and thus may be particularly relevant to schizophrenia. We therefore evaluated the association of six single nucleotide polymorphisms (SNPs) with schizophrenia in a Malaysian population (171 cases; 171 controls). Associations between these six SNPs and schizophrenia were tested in two stages. Association signals with p < 0.05 and minor allele frequency > 0.05 in stage 1 were followed by genotyping the SNPs in a replication phase (stage 2). Genotyping was performed with sequenced specific primer (PCR-SSP) and restriction fragment length polymorphism (PCR-RFLP). In our sample, we found significant associations between rs2619522 (allele p = 0.002, OR = 1.902, 95%CI = 1.266 – 2.859; genotype p = 0.002) and rs2619528 (allele p = 0.008, OR = 1.606, 95%CI = 1.130 – 2.281; genotype p = 6.18 × 10⁻⁵) and schizophrenia. Given that these two SNPs may be associated with the pathophysiology of schizophrenia, further studies on the other DTNBP1 variants are warranted.     

RASSF1A Ala133Ser polymorphism is associated with increased susceptibility to hepatocellular carcinoma in a Turkish population

Aim

The tumor suppressor gene Ras association domain family 1 isoform A (RASSF1A) regulates cell cycle regulation, apoptosis and microtubule stability and is inactivated by promoter hypermethylation at a high frequency in hepatocellular carcinoma (HCC). A guanine (G)/thymine (T) common single nucleotide polymorphism (SNP) at first position of codon 133 in RASSF1A gene determines an alanine (Ala) to serine (Ser) (Ala133Ser) amino acidic substitution which may alter cancer risk by influencing the function of RASSF1A protein.

Methods

To determine the association of the RASSF1A Ala133Ser polymorphism with the risk of HCC development in a Turkish population, a hospital-based case–control study was designed consisting of 236 subjects with HCC and 236 cancer-free control subjects matched for age, gender, smoking and alcohol status. The genotype frequency of the RASSF1A Ala133Ser polymorphism was determined by using a polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay.

Results

Allele and genotype associations of RASSF1A Ala133Ser polymorphism with HCC susceptibility were observed in comparisons between the patient and control samples (P < 0.001). Risk of HCC development in this Turkish population was significantly increased in carriers of the Ser133 variant allele of Ala133Ser polymorphism (Ala/Ser and Ser/Ser genotypes) when compared with homozygote Ala/Ala genotype (OR = 5.47, 95% CI = 3.63–8.25, P = 0.001).

Conclusion

Because our results suggest for the first time that the Ser133 allele of RASSF1A Ala133Ser polymorphism may be a genetic susceptibility factor for HCC in the Turkish population, further independent studies are required to validate our findings in a larger series, as well as in patients of different ethnic origins.     

Role of six single nucleotide polymorphisms,risk factors in coronary disease,in OLR1 alternative splicing

The OLR1 gene encodes the oxidized low-density lipoprotein receptor (LOX-1), which is responsible for the cellular uptake of oxidized LDL (Ox-LDL), foam cell formation in atheroma plaques and atherosclerotic plaque rupture. Alternative splicing (AS) of OLR1 exon 5 generates two protein isoforms with antagonistic functions in Ox-LDL uptake. Previous work identified six single nucleotide polymorphisms (SNPs) in linkage disequilibrium that influence the inclusion levels of OLR1 exon 5 and correlate with the risk of cardiovascular disease. Here we use minigenes to recapitulate the effects of two allelic series (Low- and High-Risk) on OLR1 AS and identify one SNP in intron 4 (rs3736234) as the main contributor to the differences in exon 5 inclusion, while the other SNPs in the allelic series attenuate the drastic effects of this key SNP. Bioinformatic, proteomic, mutational and functional high-throughput analyses allowed us to define regulatory sequence motifs and identify SR protein family members (SRSF1, SRSF2) and HMGA1 as factors involved in the regulation of OLR1 AS. Our results suggest that antagonism between SRSF1 and SRSF2/HMGA1, and differential recognition of their regulatory motifs depending on the identity of the rs3736234 polymorphism, influence OLR1 exon 5 inclusion and the efficiency of Ox-LDL uptake, with potential implications for atherosclerosis and coronary disease.     

The maternal folate hydrolase gene polymorphism is associated with neural tube defects in a high-risk Chinese population

Folate hydrolase 1 (FOLH1) gene encodes intestinal folate hydrolase, which regulates intestinal absorption of dietary folate. Previous studies on the association between polymorphisms rs202676 and rs61886492 and the risk of neural tube defects (NTDs) were inconclusive. A case–control study of women with NTD-affected pregnancies (n = 160) and controls (n = 320) was conducted in the Chinese population of Lvliang, a high-risk area for NTDs. We genotyped the polymorphic sites rs202676 and rs61886492 and assessed maternal plasma folate and total homocysteine (tHcy). Our results showed that in case group, plasma folate concentrations were 18 % lower compared with those of control group (8.32 vs. 6.79 nmol/L, p = 0.033) and tHcy concentrations were 17 % higher (10.47 vs. 12.65 μmol/L, p = 0.047). Almost all samples had the rs61886492 GG genotype (99.78 %). The result showed that the frequency of GG genotype in rs202676 was significantly higher in group with multiple NTDs than in controls (p = 0.030, OR = 2.157, 95 % CI, 1.06–4.38). The multiple-NTD group showed higher maternal plasma concentrations of tHcy (10.47 vs. 13.96 μmol/L, p = 0.024). The GG genotype of rs202676 had a lower maternal folate and higher tHcy concentrations than other genotypes with no significant differences. The result of structural prediction indicated that this variation might change the spatial structure of the protein. These results suggested that the maternal polymorphism rs202676 was a potential risk factor for multiple NTDs in this Chinese population. The allele G might affect maternal plasma folate and tHcy concentration.     

A missense mutation in AGTPBP1 was identified in sheep with a lower motor neuron disease

A type of lower motor neuron (LMN) disease inherited as autosomal recessive in Romney sheep was characterized with normal appearance at birth, but with progressive weakness and tetraparesis after the first week of life. Here, we carried out genome-wide homozygosity mapping using Illumina Ovine SNP50 BeadChips on lambs descended from one carrier ram, including 19 sheep diagnosed as affected and 11 of their parents that were therefore known carriers. A homozygous region of 136 consecutive single-nucleotide polymorphism (SNP) loci on chromosome 2 was common to all affected sheep and it was the basis for searching for the positional candidate genes. Other homozygous regions shared by all affected sheep spanned eight or fewer SNP loci. The 136-SNP region contained the sheep ATP/GTP-binding protein 1 (AGTPBP1) gene. Mutations in this gene have been shown to be related to Purkinje cell degeneration (pcd) phenotypes including ataxia in mice. One missense mutation c.2909G>C on exon 21 of AGTPBP1 was discovered, which induces an Arg to Pro substitution (p.Arg970Pro) at amino-acid 970, a conserved residue for the catalytic activity of AGTPBP1. Genotyping of this mutation showed 100% concordant rate with the recessive pattern of inheritance in affected, carrier, phenotypically normal and unrelated normal individuals. This is the first report showing a mutant AGTPBP1 is associated with a LMN disease in a large mammal animal model. Our finding raises the possibility of human patients with the same etiology caused by this gene or other genes in the same pathway of neuronal development.     

Further evidence that rat liver microsomal glutathione transferase 1 is not a cellular protein target for S-nitrosylation

By adopting biotin switch method, we recently reported that liver microsomal glutathione transferase 1 (MGST1) might not be a protein target for S-nitrosylation in rat microsomes or in vivo. However, alternative analytic methods are needed to confirm this observation, as a single biotin switch method in judging specific protein S-nitrosylation in biological samples is increasingly recognized as insufficient, or even unreliable. Besides, only MGST1 localized on endoplasmic reticulum (ER), but not mitochondria which favors protein S-nitrosylation was examined in the previous report. Present study was therefore carried out to address these issues. Primary cultured hepatocytes were used. A physiological existing nitric oxide (NO) donor S-nitrosoglutathione (GSNO) was adopted to trigger protein S-nitrosylation. MGST1 was immunoprecipitated and its S-nitrosothiol content was measured by the NO probe 2,3-diaminonaphthalene. In parallel, S-nitrosylated proteins were immunoprecipitated by a monoclonal anti-S-nitrosocysteine antibody and probed with an anti-MGST1 antibody. In hepatocytes, neither ER nor mitochondria were found to contain S-nitrosylated MGST1 after GSNO treatment, showing that differently distributed MGST1 was consistently un-nitrosylable in the cellular environment. But under broken cell conditions, when samples were incubated directly with GSNO, MGST1 S-nitrosylation was indeed detectable in both the microsomal and mitochondrial proteins, indicating that previous failure in detecting MGST1 S-nitrosylation in microsomes is due to the limitations of biotin switch method. These results clearly, if not definitely, demonstrate that MGST1 is not a ready candidate for S-nitrosylation in the cellular content, despite its susceptibility to S-nitrosylation under broken cell conditions.     

Missense mutation G296S in GATA4 is not responsible for cardiac septal defects

BACKGROUND:

The most common type of congenital heart disease is the cardiac septal defects, which has reported to be caused by a missense mutation (G296S) in exon 3 of the GATA4 gene.

AIMS:

The present study was undertaken to find out whether GATA4 gene is the prime cause of the septal defects in Mysore population.

MATERIALS AND METHODS:

GATA4 gene analyses were undertaken on 21 confirmed CHD cases by PCR and DNA sequencing.

RESULTS AND CONCLUSION:

Analysis of this particular mutation in 21 septal defect patients revealed that none of the patients had the mutation, indicating that this mutation is population specific or septal defect in Mysore population is caused due to mutations in other regions of the GATA4 gene.     

Genomic diversity and affinities in population groups of North West India: An analysis of Alu insertion and a single nucleotide polymorphism

The North West region of India is extremely important to understand the peopling of India, as it acted as a corridor to the foreign invaders from Eurasia and Central Asia. A series of these invasions along with multiple migrations led to intermixture of variable populations, strongly contributing to genetic variations. The present investigation was designed to explore the genetic diversities and affinities among the five major ethnic groups from North West India; Brahmin, Jat Sikh, Bania, Rajput and Gujjar. A total of 327 individuals of the abovementioned ethnic groups were analyzed for 4 Alu insertion marker loci (ACE, PV92, APO and D1) and a Single Nucleotide Polymorphism (SNP) rs2234693 in the intronic region of the ESR1 gene. Statistical analysis was performed to interpret the genetic structure and diversity of the population groups. Genotypes for ACE, APO, ESR1 and PV92 loci were found to be in Hardy–Weinberg equilibrium in all the ethnic groups, while significant departures were observed at the D1 locus in every investigated population after Bonferroni's correction. The average heterozygosity for all the loci in these ethnic groups was fairly substantial ranging from 0.3927 ± 0.1877 to 0.4333 ± 0.1416. Inbreeding coefficient indicated an overall 10% decrease in heterozygosity in these North West Indian populations. The gene differentiation among the populations was observed to be of the order of 0.013. Genetic distance estimates revealed that Gujjars were close to Banias and Jat Sikhs were close to Rajputs. Overall the study favored the recent division of the populations of North West India into largely endogamous groups. It was observed that the populations of North West India represent a more or less homogenous genetic entity, owing to their common ancestral history as well as geographical proximity.     

The NQO1 allelic frequency in hindu population of central India varies from that of other Asian populations

CONTEXT:

The enzymes encoded by the polymorphic genes NAD (P) H: quinone oxidoreductase 1 (NQO1) play an important role in the activation and inactivation of xenobiotics. This enzyme has been associated with xenobiotic related diseases, such as cancer, therapeutic failure and abnormal effects of drugs.

AIM:

The aim of the present study was to determine the allelic and genotypic frequencies of NQO Hinf I polymorphisms in a Hindu population of Central India.

SETTINGS AND DESIGN:

Polymorphisms of NQO1 were determined in 311 unrelated Hindu individuals.

MATERIALS AND METHODS:

Polymerase chain reaction- Restriction Fragment Length Polymorphism (PCR-RFLP) analysis in peripheral blood DNA for NQO1 Hinf I polymorphism was used in 311 unrelated Hindu individuals.

STATISTICAL ANALYSIS:

Allele frequencies were calculated by direct counting. Hardy Weinberg Equilibrium was evaluated using a Chi-square goodness of fit test.

RESULTS:

The observed allelic frequency was 81% for C (wild) and 19% for T (mutant) in the total sample.

CONCLUSIONS:

The allelic frequency of “C” was higher than in other Asians (57%), but similar to Caucasians (81%). The genotype distributions for Hinf I polymorphisms were in Hardy-Weinberg equilibrium.     

A Novel Deletion Mutation in CCM1 Gene (krit1) is Detected in a Chinese Family with Cerebral Cavernous Malformations

Cerebral Cavernous Malformations (CCM) are vascular malformations that are mostly located in the central nervous system (CNS) and occasionally within the skin and retina, which are classified into three types (CCM1, CCM2 and CCM3) by being located at different loci on chromosomes. At present, CCM1 (7q21), CCM2 (7p13-p15) and CCM3 (3q25.2-q27) are respectively linked to krit1 (Krev interaction trapped gene 1), MGC4607 and PDCD10 (programmed cell death 10). In this work, we identified a novel “GTA” deletion mutation at the acceptor splicing site of intron9/exon 10 on krit1. The mutation results in an abnormally spliced protein by creating a premature termination code at the 23rd amino acid downstream from the sequence alteration. Our results are consistent with previous research on krit1 mutations and confirm the conclusion that KR1T1 haploinsufficiency may be the underlying mechanism of CCM1.     

A bi-filtering method for processing single nucleotide polymorphism array data improves the quality of genetic map and accuracy of quantitative trait locus mapping in doubled haploid populations of polyploid Brassica napus

Background

Single nucleotide polymorphism (SNP) markers have a wide range of applications in crop genetics and genomics. Due to their polyploidy nature, many important crops, such as wheat, cotton and rapeseed contain a large amount of repeat and homoeologous sequences in their genomes, which imposes a huge challenge in high-throughput genotyping with sequencing and/or array technologies. Allotetraploid Brassica napus (AACC, 2n = 4x = 38) comprises of two highly homoeologous sub-genomes derived from its progenitor species B. rapa (AA, 2n = 2x = 20) and B. oleracea (CC, 2n = 2x = 18), and is an ideal species to exploit methods for reducing the interference of extensive inter-homoeologue polymorphisms (mHemi-SNPs and Pseudo-simple SNPs) between closely related sub-genomes.

Results

Based on a recent B. napus 6K SNP array, we developed a bi-filtering procedure to identify unauthentic lines in a DH population, and mHemi-SNPs and Pseudo-simple SNPs in an array data matrix. The procedure utilized both monomorphic and polymorphic SNPs in the DH population and could effectively distinguish the mHemi-SNPs and Pseudo-simple SNPs that resulted from superposition of the signals from multiple SNPs. Compared with conventional procedure for array data processing, the bi-filtering method could minimize the pseudo linkage relationship caused by the mHemi-SNPs and Pseudo-simple SNPs, thus improving the quality of SNP genetic map. Furthermore, the improved genetic map could increase the accuracies of mapping of QTLs as demonstrated by the ability to eliminate non-real QTLs in the mapping population.

Conclusions

The bi-filtering analysis of the SNP array data represents a novel approach to effectively assigning the multi-loci SNP genotypes in polyploid B. napus and may find wide applications to SNP analyses in polyploid crops.

Electronic supplementary material

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

Snf1-related protein kinase 1 is needed for growth in a normal day-night light cycle

The yeast Snf1 protein kinase and its animal homologue, the AMP-activated protein kinase, play important roles in metabolic regulation, by serving as energy gauges that turn off energy-consuming processes and mobilize energy reserves during low-energy conditions. The closest homologue of these kinases in plants is Snf1-related protein kinase 1 (SnRK1). We have cloned two SnRK1-encoding genes, PpSNF1a and PpSNF1b, in the moss Physcomitrella patens, where gene function can be studied directly by gene targeting in the haploid gametophyte. A snf1a snf1b double knockout mutant is viable, but lacks all Snf1-like protein kinase activity. The mutant has a complex phenotype that includes developmental abnormalities, premature senescence and altered sensitivities to plant hormones. Remarkably, the double knockout mutant also requires continuous light, and is unable to grow in a normal day-night light cycle. This suggests that SnRK1 is needed for metabolic changes that help the plant cope with the dark hours of the night.     

Recovering Power in Association Mapping Panels with Variable Levels of Linkage Disequilibrium

Association mapping has permitted the discovery of major QTL in many species. It can be applied to existing populations and, as a consequence, it is generally necessary to take into account structure and relatedness among individuals in the statistical model to control false positives. We analytically studied power in association studies by computing noncentrality parameter of the tests and its relationship with parameters characterizing diversity (genetic differentiation between groups and allele frequencies) and kinship between individuals. Investigation of three different maize diversity panels genotyped with the 50k SNPs array highlighted contrasted average power among panels and revealed gaps of power of classical mixed models in regions with high linkage disequilibrium (LD). These gaps could be related to the fact that markers are used for both testing association and estimating relatedness. We thus considered two alternative approaches to estimating the kinship matrix to recover power in regions of high LD. In the first one, we estimated the kinship with all the markers that are not located on the same chromosome than the tested SNP. In the second one, correlation between markers was taken into account to weight the contribution of each marker to the kinship. Simulations revealed that these two approaches were efficient to control false positives and were more powerful than classical models.     

Genome-wide association study reveals the genetic architecture of flowering time in rapeseed (Brassica napus L.)

Flowering time adaptation is a major breeding goal in the allopolyploid species Brassica napus. To investigate the genetic architecture of flowering time, a genome-wide association study (GWAS) of flowering time was conducted with a diversity panel comprising 523 B. napus cultivars and inbred lines grown in eight different environments. Genotyping was performed with a Brassica 60K Illumina Infinium SNP array. A total of 41 single-nucleotide polymorphisms (SNPs) distributed on 14 chromosomes were found to be associated with flowering time, and 12 SNPs located in the confidence intervals of quantitative trait loci (QTL) identified in previous researches based on linkage analyses. Twenty-five candidate genes were orthologous to Arabidopsis thaliana flowering genes. To further our understanding of the genetic factors influencing flowering time in different environments, GWAS was performed on two derived traits, environment sensitivity and temperature sensitivity. The most significant SNPs were found near Bn-scaff_16362_1-p380982, just 13 kb away from BnaC09g41990D, which is orthologous to A. thaliana CONSTANS (CO), an important gene in the photoperiod flowering pathway. These results provide new insights into the genetic control of flowering time in B. napus and indicate that GWAS is an effective method by which to reveal natural variations of complex traits in B. napus.     

APOA1 gene polymorphisms in the South Asian immigrant population in the United States

BACKGROUND:

Coronary artery disease (CAD) is a leading cause of death in the United States. South Asian immigrants (SAIs) from the Indian subcontinent living in the US are disproportionately at higher risk of CAD than other immigrant populations. Unique genetic factors may predispose SAIs to increased risk of developing CAD when adopting a Western lifestyle including a higher-fat diet, more sedentary behavior and additional gene-environment interactions. SAIs are known to have low levels of the protective high density lipoprotein (HDL) and an altered function for Apo-lipoprotein A-1 (ApoA1), the main protein component of HDL cholesterol. One gene that may be genetically distinctive in this population is APOA1 which codes for ApoA-1 protein, a potentially important contributing factor in the development of CAD.

MATERIALS AND METHODS:

DNA sequencing was performed to determine the status of the seven single-nucleotide polymorphisms (SNPs) in the APOA1 gene from 94 unrelated SAI adults. Genotypes, allelic frequencies, and intragenic linkage disequilibrium of the APOA1 SNPs were calculated.

RESULTS:

Several polymorphisms and patterns were common among persons of south Asian ethnicity. Frequencies for SNPs T655C, T756C and T1001C were found to be different than those reported in European Caucasian individuals. Linkage disequilibrium was found to be present between most (13 of 15) SNP pairings indicating common inheritance patterns.

CONCLUSIONS:

SAIs showed variability in the sequence of the APOA1 gene and linkage disequilibrium for most SNPS. This pattern of APOA1 SNPs may contribute to decreased levels of HDL cholesterol reported in SAIs, leading to an increased risk for developing CAD in this population.