Multiplex SNP genotyping in pooled DNA samples by a four-colour microarray system

We selected 125 candidate single nucleotide polymorphisms (SNPs) in genes belonging to the human type 1 interferon (IFN) gene family and the genes coding for proteins in the main type 1 IFN signalling pathway by screening databases and by in silico comparison of DNA sequences. Using quantitative analysis of pooled DNA samples by solid-phase mini-sequencing, we found that only 20% of the candidate SNPs were polymorphic in the Finnish and Swedish populations. To allow more effective validation of candidate SNPs, we developed a four-colour microarray-based mini-sequencing assay for multiplex, quantitative allele frequency determination in pooled DNA samples. We used cyclic mini-sequencing reactions with primers carrying 5′-tag sequences, followed by capture of the products on microarrays by hybridisation to complementary tag oligonucleotides. Standard curves prepared from mixtures of known amounts of SNP alleles demonstrate the applicability of the system to quantitative analysis, and showed that for about half of the tested SNPs the limit of detection for the minority allele was below 5%. The microarray-based genotyping system established here is universally applicable for genotyping and quantification of any SNP, and the validated system for SNPs in type 1 IFN-related genes should find many applications in genetic studies of this important immunoregulatory pathway.     

Haplotype distribution and association of candidate genes with salt tolerance in Indian wild rice germplasm

Key message

The association of natural genetic variations of salt-responsive candidate genes belonging to different gene families with salt-tolerance phenotype and their haplotype variation in different geographic regions.

Abstract

Soil salinity covers a large part of the arable land of the world and is a major factor for yield losses in salt-sensitive crops, such as rice. Different gene families that respond to salinity have been identified in rice, but limited success has been achieved in developing salt-tolerant cultivars. Therefore, 21 salt stress-responsive candidate genes belonging to different gene families were re-sequenced to analyse their genetic variation and association with salt tolerance. The average single nucleotide polymorphism (SNP) density was 16 SNPs per kbp amongst these genes. The identified nucleotide and haplotype diversity showed comparatively higher genetic variation in the transporter family genes. Linkage disequilibrium (LD) analysis showed significant associations of SNPs in BADH2, HsfC1B, MIPS1, MIPS2, MYB2, NHX1, NHX2, NHX3, P5CS1, P5CS2, PIP1, SIK1, SOS1, and SOS2 genes with the salt-tolerant phenotype. A combined analysis of SNPs in the 21 candidate genes and eight other HKT transporter genes produced two separate clusters of tolerant genotypes, carrying unique SNPs in the ion transporter and osmoticum-related genes. Haplotype network analysis showed all the major and few minor alleles distributed over distant geographic regions. Minor haplotypes may be recently evolved alleles which migrated to distant geographic regions and may represent recent expansion of Indian wild rice. The analysis of genetic variation in different gene families identified the relationship between adaptive variations and functional significance of the genes. Introgression of the identified alleles from wild relatives may enhance the salt tolerance and consequently rice production in the salinity-affected areas.

     

Associations between two novel rSNPs in 5'-flanking region of the bovine casein gene cluster and milk performance traits

Ten evolutionary conservative sequences with high identity level to homological sequences in other mammal species were revealed in 5'-flanking region of casein's genes cluster. Five novel SNPs located inside of the evolutionary conservative regions were identified. The binding sites were revealed to be present in one allelic variant of four detected SNPs. So these SNPs were considered as rSNPs. Significant differences of allelic frequencies were revealed between beef cow's group and dairy cow's group in two rSNPs (NCE4, NCE7, p<0.001). Different alleles of those two rSNPs were shown to be associated with some milk performance traits in Black-and-White Holstein dairy cows. Significant difference of protein percentage has been found between cows with G/G and A/A genotypes (P<0.05) and A/G and A/A genotypes (P<0.05) for NCE4 polymorphism. The groups of animals with genotypes G/G and A/G for NCE7 polymorphism were significantly different in milk yield at the first lactation (kg) (P<0.01), milk fat yield (kg) (P<0.05) and milk protein yield (kg) (P<0.01). For the last trait the difference was significant also between cows with genotypes G/G and A/A for rSNP NCE7 (P<0.05).     

Identification of 96 single nucleotide polymorphisms in eight genes involved in iron metabolism: efficiency of bioinformatic extraction compared with a systematic sequencing approach

Single nucleotide polymorphisms (SNPs) can significantly contribute to the characterization of the genes predisposing to iron overloads or deficiencies. We report an SNP survey of coding and non-coding regions of eight genes involved in iron metabolism, by two successive methods. First, we made use of the public domain sequence data, by using assembled expressed sequence tags, non-redundant sequences, and SNP database screening. We extracted 77 potential SNPs of which only 31 could be further validated by sequencing DNA from 44 unrelated multi-ethnic individuals. Our results indicate that a bioinformatic approach may be effective only in those cases where candidate SNPs are extracted from two different data sources or in cases of experimentally confirmed SNPs. Second, additional systematic sequencing of DNA from 24 unrelated Breton subjects increased the number of SNPs over a total length of 86 kb to 96. The average distance between the SNPs and minor allele frequencies were higher than reported by others authors; this discrepancy may reflect the nature of the genes studied and the ethnic homogeneity of our test population.     

A catalog of polymorphisms falling in microRNA-binding regions of cancer genes

Recent evidence indicates that small, nonprotein-coding RNA molecules, called microRNAs (miRNAs), control cell growth, differentiation, and apoptosis, and are also involved in tumorigenesis. miRNAs can bind to the 3' untranslated regions (3'UTRs) of messenger RNAs and interfere with their translation. We hypothesized that common polymorphisms within their genes or within their targets could have an important impact for an individual's risk to develop complex diseases. In this study, we selected the 3'UTRs of 129 genes involved in pathways commonly acknowledged as important for cancer, and we identified putative miRNA-binding sites by means of specialized algorithms (PicTar, DIANA-MicroT, miRBase, miRanda, TargetScan, and MicroInspector). Then we investigated 79 single-nucleotide polymorphisms (SNPs) within the putative binding sites for their ability to affect or impair the binding with the miRNA by assessing the DeltaDeltaG, the variation of DeltaG (Gibbs free energy), through comparing the wild-type and their corresponding variant alleles. Moreover, we reported seven identified SNPs in seven pre-miRNA hairpin regions and one SNP in the mature sequence of miR-608. Considering the validation status of the SNPs and their frequencies, we found at least 23 candidate polymorphisms of biological relevance that we propose for further investigation in case-control association studies.     

Genomic regions involved in yield potential detected by genome-wide association analysis in Japanese high-yielding rice cultivars

Background

High-yielding cultivars of rice (Oryza sativa L.) have been developed in Japan from crosses between overseas indica and domestic japonica cultivars. Recently, next-generation sequencing technology and high-throughput genotyping systems have shown many single-nucleotide polymorphisms (SNPs) that are proving useful for detailed analysis of genome composition. These SNPs can be used in genome-wide association studies to detect candidate genome regions associated with economically important traits. In this study, we used a custom SNP set to identify introgressed chromosomal regions in a set of high-yielding Japanese rice cultivars, and we performed an association study to identify genome regions associated with yield.

Results

An informative set of 1152 SNPs was established by screening 14 high-yielding or primary ancestral cultivars for 5760 validated SNPs. Analysis of the population structure of high-yielding cultivars showed three genome types: japonica-type, indica-type and a mixture of the two. SNP allele frequencies showed several regions derived predominantly from one of the two parental genome types. Distinct regions skewed for the presence of parental alleles were observed on chromosomes 1, 2, 7, 8, 11 and 12 (indica) and on chromosomes 1, 2 and 6 (japonica). A possible relationship between these introgressed regions and six yield traits (blast susceptibility, heading date, length of unhusked seeds, number of panicles, surface area of unhusked seeds and 1000-grain weight) was detected in eight genome regions dominated by alleles of one parental origin. Two of these regions were near Ghd7, a heading date locus, and Pi-ta, a blast resistance locus. The allele types (i.e., japonica or indica) of significant SNPs coincided with those previously reported for candidate genes Ghd7 and Pi-ta.

Conclusions

Introgression breeding is an established strategy for the accumulation of QTLs and genes controlling high yield. Our custom SNP set is an effective tool for the identification of introgressed genome regions from a particular genetic background. This study demonstrates that changes in genome structure occurred during artificial selection for high yield, and provides information on several genomic regions associated with yield performance.

Electronic supplementary material

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

Patterns of linkage disequilibrium and haplotype distribution in disease candidate genes

Background

The adequacy of association studies for complex diseases depends critically on the existence of linkage disequilibrium (LD) between functional alleles and surrounding SNP markers.

Results

We examined the patterns of LD and haplotype distribution in eight candidate genes for osteoporosis and/or obesity using 31 SNPs in 1,873 subjects. These eight genes are apolipoprotein E (APOE), type I collagen α1 (COL1A1), estrogen receptor-α (ER-α), leptin receptor (LEPR), parathyroid hormone (PTH)/PTH-related peptide receptor type 1 (PTHR1), transforming growth factor-β1 (TGF-β1), uncoupling protein 3 (UCP3), and vitamin D (1,25-dihydroxyvitamin D₃) receptor (VDR). Yin yang haplotypes, two high-frequency haplotypes composed of completely mismatching SNP alleles, were examined. To quantify LD patterns, two common measures of LD, D' and r², were calculated for the SNPs within the genes. The haplotype distribution varied in the different genes. Yin yang haplotypes were observed only in PTHR1 and UCP3. D' ranged from 0.020 to 1.000 with the average of 0.475, whereas the average r² was 0.158 (ranging from 0.000 to 0.883). A decay of LD was observed as the intermarker distance increased, however, there was a great difference in LD characteristics of different genes or even in different regions within gene.

Conclusion

The differences in haplotype distributions and LD patterns among the genes underscore the importance of characterizing genomic regions of interest prior to association studies.     

Regulatory SNPs in complex diseases: their identification and functional validation

Finding the genetic causes for complex diseases is a challenge. Expression studies have shown that the level of expression of many genes is altered in disease compared with normal conditions, but what lies behind these changes? Linkage studies provide hints as to where in the genome the genetic triggers--the mutations--might be located. Fine-mapping and association studies can give yet more information about which genes, and which changes in the genes, are involved in the disease. Recent examples show that single-nucleotide polymorphisms (SNPs), which are variations at the single-nucleotide level within an individual's DNA, in the regulatory regions of some genes constitute susceptibility factors in many complex diseases. This article discusses the nature of regulatory SNPs (rSNPs) and techniques for their functional validation, and looks towards what rSNPs can tell us about complex diseases.     

Applications of single nucleotide polymorphisms in crop genetics

The discovery of single nucleotide polymorphisms (SNPs) and insertions/deletions, which are the basis of most differences between alleles, has been simplified by recent developments in sequencing technology. SNP discovery in many crop species, such as corn and soybean, is relatively straightforward because of their high level of intraspecific nucleotide diversity, and the availability of many gene and expressed sequence tag (EST) sequences. For these species, direct readout of SNP haplotypes is possible. Haplotype-based analysis is more informative than analysis based on individual SNPs, and has more power in analyzing association with phenotypes. The elite germplasm of some crops may have been subjected to bottlenecks relatively recently, increasing the amount of linkage disequilibrium (LD) present and facilitating the association of SNP haplotypes at candidate gene loci with phenotypes. Whole-genome scans may help identify genome regions that are associated with interesting phenotypes if sufficient LD is present. Technological improvements make the use of SNP and indel markers attractive for high-throughput use in marker-assisted breeding, EST mapping and the integration of genetic and physical maps.     

Next‐generation sequencing for identification of candidate genes for Fusarium wilt and sterility mosaic disease in pigeonpea (Cajanus cajan)

To map resistance genes for Fusarium wilt (FW) and sterility mosaic disease (SMD) in pigeonpea, sequencing‐based bulked segregant analysis (Seq‐BSA) was used. Resistant (R) and susceptible (S) bulks from the extreme recombinant inbred lines of ICPL 20096 × ICPL 332 were sequenced. Subsequently, SNP index was calculated between R‐ and S‐bulks with the help of draft genome sequence and reference‐guided assembly of ICPL 20096 (resistant parent). Seq‐BSA has provided seven candidate SNPs for FW and SMD resistance in pigeonpea. In parallel, four additional genotypes were re‐sequenced and their combined analysis with R‐ and S‐bulks has provided a total of 8362 nonsynonymous (ns) SNPs. Of 8362 nsSNPs, 60 were found within the 2‐Mb flanking regions of seven candidate SNPs identified through Seq‐BSA. Haplotype analysis narrowed down to eight nsSNPs in seven genes. These eight nsSNPs were further validated by re‐sequencing 11 genotypes that are resistant and susceptible to FW and SMD. This analysis revealed association of four candidate nsSNPs in four genes with FW resistance and four candidate nsSNPs in three genes with SMD resistance. Further, In silico protein analysis and expression profiling identified two most promising candidate genes namely C.cajan_01839 for SMD resistance and C.cajan_03203 for FW resistance. Identified candidate genomic regions/SNPs will be useful for genomics‐assisted breeding in pigeonpea.     

Genome-wide association study of shank length and diameter at different developmental stages in chicken F2 resource population

In order to find SNPs and genes affecting shank traits, we performed a GWAS in a chicken F2 population of eight half-sib families from five hatches derived from reciprocal crosses between an Arian fast-growing line and an Urmia indigenous slow-growing chicken. A total of 308 birds were genotyped using a 60K chicken SNP chip. Shank traits including shank length and diameter were measured weekly from birth to 12 weeks of age. A generalized linear model and a compressed mixed linear model (CMLM) were applied to achieve the significant regions. The value of the average genomic inflation factor (λ statistic) of the CMLM model (0.99) indicated that the CMLM was more effective than the generalized linear model in controlling the population structure. The genes surrounding significant SNPs and their biological functions were identified from NCBI, Ensembl and UniProt databases. The results indicated that 12 SNPs at 12 different ages passed the LD-adjusted 5% Bonferroni significant threshold. Two SNPs were significant for shank length and nine SNPs were significant for shank diameter. The significant SNPs were located near to or inside 11 candidate genes. The results showed that a number of significant SNPs in the middle ages were higher than the rest. The MXRA8 gene was related to the significant SNP at week 1 that promotes proliferation of growth plate chondrocytes. A unique SNP of Gga_rs16689511 located on chicken Z chromosome within the LOC101747628 gene was related to shank length at three different ages of birds (weeks 8, 9 and 11). The significant SNPs for shank diameter were found at weeks 4 and 7 (four and five SNPs respectively). The identifications of SNPs and genes here could contribute to a better understanding of the genetic control of shank traits in chicken.     

Sequences associated with human iris pigmentation

To determine whether and how common polymorphisms are associated with natural distributions of iris colors, we surveyed 851 individuals of mainly European descent at 335 SNP loci in 13 pigmentation genes and 419 other SNPs distributed throughout the genome and known or thought to be informative for certain elements of population structure. We identified numerous SNPs, haplotypes, and diplotypes (diploid pairs of haplotypes) within the OCA2, MYO5A, TYRP1, AIM, DCT, and TYR genes and the CYP1A2-15q22-ter, CYP1B1-2p21, CYP2C8-10q23, CYP2C9-10q24, and MAOA-Xp11.4 regions as significantly associated with iris colors. Half of the associated SNPs were located on chromosome 15, which corresponds with results that others have previously obtained from linkage analysis. We identified 5 additional genes (ASIP, MC1R, POMC, and SILV) and one additional region (GSTT2-22q11.23) with haplotype and/or diplotypes, but not individual SNP alleles associated with iris colors. For most of the genes, multilocus gene-wise genotype sequences were more strongly associated with iris colors than were haplotypes or SNP alleles. Diplotypes for these genes explain 15% of iris color variation. Apart from representing the first comprehensive candidate gene study for variable iris pigmentation and constituting a first step toward developing a classification model for the inference of iris color from DNA, our results suggest that cryptic population structure might serve as a leverage tool for complex trait gene mapping if genomes are screened with the appropriate ancestry informative markers.     

Allelic variation in gene expression identified through computational analysis of the dbEST database

Differential expression between the two alleles of an individual and between people with different genotypes has been commonly observed. Quantitative differences in gene expression between people may provide the genetic basis for the phenotypic difference between individuals and may be the primary cause of complex diseases. In this paper, we developed a computational method to identify genes that displayed allelic variation in gene expression in human EST libraries. To model allele-specific gene expression, we first identified EST libraries in which both A and B alleles were expressed and then identified allelic variation in gene expression based on the EST counts for each allele using a binomial test. Among 1107 SNPs that had a sufficient number of ESTs for the analysis, 524 (47%) displayed allelic variation in at least one cDNA library. We verified experimentally the allelic variation in gene expression for 6 of these SNPs. The frequency of allelic variation observed in EST libraries was similar to the previous studies using the SNP chip and primer extension method. We found that genes that displayed allelic variation were distributed throughout the human genome and were enriched in certain chromosome regions. The SNPs and genes identified in this study will provide a rich source for evaluating the effects of those SNPs and associated haplotypes in human health and diseases.