Population frequencies of single nucleotide polymorphisms (SNPs) in immuno-modulatory genes

Inherited polymorphisms in immuno-modulatory genes may contribute to variations in immune function and genetic susceptibility for complex diseases, including cancer. We report results from a comprehensive study to discover novel single nucleotide polymorphisms (SNPs) and to estimate allelic frequency for both novel and known coding and regulatory region SNPs in genes encoding proteins that have been implicated in the immune response to tumors. We identified 12 novel nucleotide substitution variants and one deletion variant in 17 genes analyzed (TGFBETA;R, BETA;2M, IFNGAMMA;, TNFALPHA;, TNFALPHA;R, LTALPHA;, IL-6, IL-12, IL-2, IL-1ALPHA;, IL-1BETA;, IL-1RN, IL-10, CTLA4, CD40L, FAS and FASL). We determined the frequency of these novel polymorphisms, as well as 17 previously identified polymorphisms, in a control sample of 158 individuals, approximately half of which were Caucasian (n = 74) and half of which were African American (n = 84). Significant differences in allele frequencies were observed between the two racial groups for 13/17 genes tested. These allelic variations maybe associated with alterations in immune function and thus susceptibility to a number of complex disease states such as cancer.     

TLR4 single nucleotide polymorphisms (SNPs) associated with Salmonella shedding in pigs

Toll-like receptor 4 (TLR4) is a key factor in the innate immune recognition of lipopolysaccharide (LPS) from Gram-negative bacteria. Previous studies from our group identified differences in the expression profile of TLR4 and genes affected by the TLR4 signaling pathway among pigs that shed varying levels of Salmonella, a Gram-negative bacterium. Therefore, genetic variation in this gene may be involved with the host’s immune response to bacterial infections. The current study screened for single nucleotide polymorphisms (SNPs) in the TLR4 gene and tested their association with Salmonella fecal shedding. Pigs (n?=?117) were intranasally challenged at 7 weeks of age with 1?×?10⁹ CFU of S. Typhimurium χ4232 and were classified as low or persistent Salmonella shedders based on the levels of Salmonella being excreted in fecal material. Salmonella fecal shedding was determined by quantitative bacteriology on days 2, 7, 14, and 20/21 post exposure, and the cumulative levels of Salmonella were calculated to identify the low (n?=?20) and persistent (n?=?20) Salmonella shedder pigs. From those 40 animals, the TLR4 region was sequenced, and 18 single nucleotide polymorphisms (SNPs) in TLR4 were identified. Twelve SNPs have been previously described and six are novel SNPs of which five are in the 5′ untranslated region and one is in intron 2. Single marker association test identified 13 SNPs associated with the qualitative trait of Salmonella fecal shedding, and seven of those SNPs were also associated with a quantitative measurement of fecal shedding (P?<?0.05). Using a stepwise regression process, a haplotype composed of SNPs rs80787918 and rs80907449 (P?≤?4.0?×?10^?3) spanning a region of 4.9 Kb was identified, thereby providing additional information of the influence of those SNPs on Salmonella fecal shedding in pigs.     

Single nucleotide polymorphisms (SNPs) in human lactoferrin gene.

The lactoferrin protein possesses antimicrobial and antiviral activities. It is also involved in the modulation of the immune response. In a normal healthy individual, lactoferrin plays a role in the front-line host defense against infection and in immune and inflammatory responses. Whether genomic variations, such as single nucleotide polymorphisms (SNPs), have an effect on the structure and function of lactoferrin protein and whether these variations contribute to the different susceptibility of individuals in response to environmental insults are interesting health-related issues. In this study, the lactoferrin gene was resequenced as part of the Environmental Genome Project of the National Institute of Environmental Health Sciences, which operates within the National Institutes of Health. Ninety-one healthy donors of different ethnicities were used to establish common SNPs in the exons of the lactoferrin gene in the general population. The data will serve as a basis from which study the association of lactoferrin polymorphism and disease.     

African American-preponderant single nucleotide polymorphisms (SNPs) and risk of breast cancer

Background: African American women more often present with more aggressive types of breast cancer than Caucasian women, but little is known whether genetic polymorphisms specific to or disproportionate in African Americans are associated with their risk of breast cancer. Methods: A population-based case-control study was conducted including 194 cases identified through the Metropolitan Detroit Cancer Surveillance System and 189 controls recruited through random digit dialing to examine polymorphisms in genes involved in estrogen metabolism and action. Results: The African American-specific CYP1A1 5639C allele was associated with an increased risk of breast cancer (odds ratio (OR) = 2.34, 95% confidence interval (CI) 1.23–4.44) and this association with the CYP1A1 5639 locus was dependent on another polymorphism in the CYP3A4 gene (P = 0.043 for the interaction). In addition, African American-predominant CYP1B1 432 Val allele was significantly more often found in the cases than in the controls overall and the HSD17B1 312 Gly allele was specifically associated with premenopausal breast cancer risk (OR = 3.00, 95%CI 1.29–6.99). Conclusion: These observations need to be confirmed in larger studies due to the limited statistical power of the study based on a small number of cases.     

Assessment of population structure by single nucleotide polymorphisms (SNPs) in goat breeds

Single nucleotide polymorphisms (SNPs) may be used in biodiversity studies and commercial tasks like traceability, paternity testing and selection for suitable genotypes. Twenty-seven SNPs were characterized and genotyped on 250 individuals belonging to eight Italian goat breeds. Multilocus genotype data were used to infer population structure and assign individuals to populations. To estimate the number of groups (K) to test in population structure analysis we used likelihood values and variance of the bootstrap samples, deriving optimal K from a drop in the likelihood and a rise in the variance plots against K.     

Identification of novel single nucleotide polymorphisms (SNPs) of the lipoprotein lipase (LPL) gene associated with fatty acid composition in Korean cattle

The lipoprotein lipase (LPL) gene can be considered a functional candidate gene that regulates fatty acid composition. In this study, genetic associations between fatty acid composition and exonic single nucleotide polymorphisms (SNPs) in the LPL gene were examined using 612 Korean cattle. We investigated the relationship between unsaturated fatty acids and five novel SNPs (c.322G>A, c.329A>T, c.527T>G, c.988C>T and c.1591G>A), and confirmed that three polymorphic SNPs (c.322G>A, c.329A>T and c.1591G>A) were associated with fatty acid composition. Korean cattle with an AA genotype of c.322G>A, c.329A>T, and GA genotype of c.1591G>A had higher levels of monounsaturated fatty acids and carcass traits (P < 0.05). Our findings confirmed that three novel SNPs we identified in the LPL gene can affect fatty acid composition and carcass traits. Therefore, selection for AA and GA genotypes should be recommended to genetically improve beef quality and flavor.     

Mutational spectrum in the recent human genome inferred by single nucleotide polymorphisms

So far, there is no genome-wide estimation of the mutational spectrum in humans. In this study, we systematically examined the directionality of the point mutations and maintenance of GC content in the human genome using approximately 1.8 million high-quality human single nucleotide polymorphisms and their ancestral sequences in chimpanzees. The frequency of C-->T (G-->A) changes was the highest among all mutation types and the frequency of each type of transition was approximately fourfold that of each type of transversion. In intergenic regions, when the GC content increased, the frequency of changes from G or C increased. In exons, the frequency of G:C-->A:T was the highest among the genomic categories and contributed mainly by the frequent mutations at the CpG sites. In contrast, mutations at the CpG sites, or CpG-->TpG/CpA mutations, occurred less frequently in the CpG islands relative to intergenic regions with similar GC content. Our results suggest that the GC content is overall not in equilibrium in the human genome, with a trend toward shifting the human genome to be AT rich and shifting the GC content of a region to approach the genome average. Our results, which differ from previous estimates based on limited loci or on the rodent lineage, provide the first representative and reliable mutational spectrum in the recent human genome and categorized genomic regions.     

Sequence context analysis in the mouse genome: single nucleotide polymorphisms and CpG island sequences

A genome-wide view of sequence mutability in mice is still limited, although biologists usually assume the same scenario for mice as for humans. In this study, we examined the sequence context in the local environment of 482,528 mouse single nucleotide polymorphisms (SNPs). We found that CpG-containing short sequences, in general, had more representation in the local sequences of SNPs compared to the genome sequences. The extent of this overrepresentation was stronger in mice than in humans, which is inconsistent with previous observations of the weaker neighboring-nucleotide biases on mouse SNPs. To exclude the CpG effect, we compared the distribution patterns of short sequences among the six categories of SNPs. The results revealed an even stronger pattern in the CpG-containing group for C/G substitution compared to for A/G or C/T substitutions. We next performed the first genome-wide sequence context analysis of SNPs in the mouse CpG islands. SNPs occurring at CpG sites were 3.14-fold less prevalent than expected, suggesting the suppression of methylation-dependent deamination in the CpG islands. The extent of this suppression was less in mice than in humans. Finally, compared with humans, the observations of a greater deficit of CpG dinucleotides, a stronger overrepresentation of CpG-containing n-mers surrounding the polymorphic sites, and a higher SNP/genome ratio of CpG dinucleotides in the mouse genome support the "loss of CpG islands" model in the mouse lineage.     

Linkage disequilibrium grouping of single nucleotide polymorphisms (SNPs) reflecting haplotype phylogeny for efficient selection of tag SNPs

Single nucleotide polymorphisms (SNPs) have been proposed to be grouped into haplotype blocks harboring a limited number of haplotypes. Within each block, the portion of haplotypes is expected to be tagged by a selected subset of SNPs; however, none of the proposed selection algorithms have been definitive. To address this issue, we developed a tag SNP selection algorithm based on grouping of SNPs by the linkage disequilibrium (LD) coefficient r(2) and examined five genes in three ethnic populations--the Japanese, African Americans, and Caucasians. Additionally, we investigated ethnic diversity by characterizing 979 SNPs distributed throughout the genome. Our algorithm could spare 60% of SNPs required for genotyping and limit the imprecision in allele-frequency estimation of nontag SNPs to 2% on average. We discovered the presence of a mosaic pattern of LD plots within a conventionally inferred haplotype block. This emerged because multiple groups of SNPs with strong intragroup LD were mingled in their physical positions. The pattern of LD plots showed some similarity, but the details of tag SNPs were not entirely concordant among three populations. Consequently, our algorithm utilizing LD grouping allows selection of a more faithful set of tag SNPs than do previous algorithms utilizing haplotype blocks.     

Characterization of 12 single nucleotide polymorphisms (SNPs) in Pacific abalone, Haliotis discus hannai

We report here for the first time 12 polymorphic single nucleotide polymorphisms (SNPs) in a commercially important gastropod, Pacific abalone (Haliotis discus hannai) that were identified by searching expressed sequence tag database. These SNP loci (seven nuclear and five mitochondrial SNPs) were polymorphic among 37 wild abalone individuals, based on a four-primer allele-specific polymerase chain reaction analysis. All loci had two alleles and the minor allele frequency ranged from 0.027 to 0.473. For the seven nuclear SNPs, the expected and observed heterozygosities ranged from 0.053 to 0.499 and from 0.054 to 0.811, respectively.     

Single nucleotide polymorphisms (SNPs) that map to gaps in the human SNP map

An international effort is underway to generate a comprehensive haplotype map (HapMap) of the human genome represented by an estimated 300000 to 1 million ‘tag’ single nucleotide polymorphisms (SNPs). Our analysis indicates that the current human SNP map is not sufficiently dense to support the HapMap project. For example, 24.6% of the genome currently lacks SNPs at the minimal density and spacing that would be required to construct even a conservative tag SNP map containing 300 000 SNPs. In an effort to improve the human SNP map, we identified 140 696 additional SNP candidates using a new bioinformatics pipeline. Over 51 000 of these SNPs mapped to the largest gaps in the human SNP map, leading to significant improvements in these regions. Our SNPs will be immediately useful for the HapMap project, and will allow for the inclusion of many additional genomic intervals in the final HapMap. Nevertheless, our results also indicate that additional SNP discovery projects will be required both to define the haplotype architecture of the human genome and to construct comprehensive tag SNP maps that will be useful for genetic linkage studies in humans.     

Association of single nucleotide polymorphisms (SNPs) of PADI4 gene with rheumatoid arthritis (RA) in Indian population

Rheumatoid arthritis (RA) is a chronic inflammatory autoimmune disease affecting ~ 1% of the population worldwide. The genome wide association studies on RA patients revealed linkage with 1p36 locus containing peptidyl arginine deiminase 4 (PADI4) genes. Case-control association studies and mRNA stability assays reported the association of PADI4 gene with RA in Korean and Japanese populations. However, such association was not found in Spanish population. Differences in the association of PADI4 with RA in different populations prompted the present study in Indian population. Anti-CCP antibodies, RF antibody, disease activity scores at 28 joints (DAS28) and mutations in three exons of PADI4 were investigated in RA patients and control group. Among the patients anti-CCP antibody levels were found to be associated with high DAS28 values (r = 0.4526, P < 0.0001). Polymorphism in exon-4 (padi4_104, [rs1748033]) of PADI4 showed significant association of 'C' allele with RA in the study population (P = 0.0008). Polymorphism in exon-3 (padi4_92, [rs874881]) also exhibited moderate association with the disease (P = 0.075). However, no association of the disease was found with the SNPs padi4_89 [rs11203366] and padi4_90 [rs11203367] in exon-2 of PADI4.     

Tailed primer base excision sequence scanning (TP-BESS) for detection of single nucleotide polymorphisms (SNPs)

Base excision sequence scanning (BESS) is used to detect and localise point mutations in mammalian genes. The cost of BESS can be significant in large-scale projects. however, due to the requirement for end-labelling of one of the two PCR primers. This is especially true when a fluorescent label is used for detection. If a universal label could be used for all of the PCR reactions, it would reduce the cost of this technique significantly. Here we describe a TP-BESS procedure where one fluorescence-labelled primer is used as a universal label for all the PCR reactions in BESS analysis. The universal label makes BESS financially more feasible for large-scale detection of single nucleotide polymorphisms (SNPs).     

The association of two single nucleotide polymorphisms (SNPs) in growth hormone (GH) gene with litter size and superovulation response in goat-breeds

Two active mutations (A 781 G and A 1575 G) in growth hormone (GH) gene, and their associations with litter size (LS), were investigated in both a high prolificacy (Matou, n = 182) and a low prolificacy breed (Boer, n = 352) by using the PCR-RFLP method. Superovulation experiments were designed in 57 dams, in order to evaluate the effect of different genotypes of the GH gene on superovulation response. Two genotypes (AA and AB, CC and CD) in each mutation were detected in these two goat breeds. Neither BB nor DD homozygous genotypes were observed. The genotypic frequencies of AB and CC were significantly higher than those of AA and CD. In the third parity, Matou dams with AB or CC genotypes had significantly larger litter sizes than those with AA and CD (p < 0.05). On combining the two loci, both Matou and Boer dams with ABCD genotype had the largest litter sizes when compared to the other genotypes (p < 0.05). When undergoing like superovulation treatments, a significantly higher number of corpora lutea and ova, with a lower incidence of ovarian cysts, were harvested in the AB and CC genotypes than in AA and CD. These results show that the two loci of GH gene are highly associated with abundant prolificacy and superovulation response in goat breeds.     

Amino acid substitutions in the human genome: evolutionary implications of single nucleotide polymorphisms

Functional differences between amino acids have long been of interest in understanding protein evolution. Several indices exist for comparing residues on the basis of their physicochemical properties and frequencies of occurrence in conserved protein alignments. Here we present a residue dissimilarity index based on coding single nucleotide polymorphisms (SNPs) in the human genome. The index represents an average, organism-wide set of differences between residues and provides important insight into evolutionary restraints on residue substitutions in the human genome. Unlike previous models, it is not restricted to highly conserved protein structures, nor confounded by evolutionary differences between species. Our results confirm earlier observations regarding residue mutabilities but also suggest that in addition to the established key properties, such as size and polarity, charge conservation may be an important and currently underestimated factor in protein evolution. We also estimate that less than 51% of amino acid substitutions occurring in the human genome are evolutionarily neutral.     

EFIN: predicting the functional impact of nonsynonymous single nucleotide polymorphisms in human genome

Background

Predicting the functional impact of amino acid substitutions (AAS) caused by nonsynonymous single nucleotide polymorphisms (nsSNPs) is becoming increasingly important as more and more novel variants are being discovered. Bioinformatics analysis is essential to predict potentially causal or contributing AAS to human diseases for further analysis, as for each genome, thousands of rare or private AAS exist and only a very small number of which are related to an underlying disease. Existing algorithms in this field still have high false prediction rate and novel development is needed to take full advantage of vast amount of genomic data.

Results

Here we report a novel algorithm that features two innovative changes: 1. making better use of sequence conservation information by grouping the homologous protein sequences into six blocks according to evolutionary distances to human and evaluating sequence conservation in each block independently, and 2. including as many such homologous sequences as possible in analyses. Random forests are used to evaluate sequence conservation in each block and to predict potential impact of an AAS on protein function. Testing of this algorithm on a comprehensive dataset showed significant improvement on prediction accuracy upon currently widely-used programs. The algorithm and a web-based application tool implementing it, EFIN (Evaluation of Functional Impact of Nonsynonymous SNPs) were made freely available (http://paed.hku.hk/efin/) to the public.

Conclusions

Grouping homologous sequences into different blocks according to the evolutionary distance of the species to human and evaluating sequence conservation in each group independently significantly improved prediction accuracy. This approach may help us better understand the roles of genetic variants in human disease and health.

Electronic supplementary material

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