Paternity validation and estimation of genotyping error rate for the BovineSNP50 BeadChip

Incorrect paternity assignment in cattle can have a major effect on rates of genetic gain. Of the 576 Israeli Holstein bulls genotyped by the BovineSNP50 BeadChip, there were 204 bulls for which the father was also genotyped. The results of 38 828 valid single nucleotide polymorphisms (SNPs) were used to validate paternity, determine the genotyping error rates and determine criteria enabling deletion of defective SNPs from further analysis. Based on the criterion of >2% conflicts between the genotype of the putative sire and son, paternity was rejected for seven bulls (3.5%). The remaining bulls had fewer conflicts by one or two orders of magnitude. Excluding these seven bulls, all other discrepancies between sire and son genotypes are assumed to be caused by genotyping mistakes. The frequency of discrepancies was >0.07 for nine SNPs, and >0.025 for 81 SNPs. The overall frequency of discrepancies was reduced from 0.00017 to 0.00010 after deletion of these 81 SNPs, and the total expected fraction of genotyping errors was estimated to be 0.05%. Paternity of bulls that are genotyped for genomic selection may be verified or traced against candidate sires at virtually no additional cost.     

Maternity validation using sire‐only BovineSNP50 BeadChip data

Based on pairwise identity‐by‐state (IBS) distances and whole‐genome SNP data, kinship was investigated in the Israeli Holstein population. A total of 789 bulls, including most of the artificial insemination sires in service since 1987, were genotyped by the BovineSNP50 BeadChip. This sample included up to five generations. For each bull‐by‐bull combination, three states are possible for each marker: no match, a single match and both alleles match. Summing over all markers, the 932 598 IBS scores (three match frequencies*310 866 bull‐by‐bull combinations) were visualized using three‐dimensional coordinates that corresponded to the frequencies of the three possible states. Results were reduced to two dimensions using the transformations x’ = 0.7071(1 + freq1?freq2) and y’ = 1.2247freq0. Bull‐by‐bull pairs were grouped according to their level of kinship, and canonical scores were calculated using discriminant analysis and the x’ and y’ features. Of the 474 pairs of recorded maternal grandsire–grandson with both individuals genotyped, the probability for 28 pairs to belong to this level of kinship was low (P < 0.05), suggesting an error rate of around 3% per generation in pedigree determination.     

Genome-wide mapping of copy number variations in commercial hybrid pigs using a high-density SNP genotyping array

Copy number variations (CNVs) are important forms of structural variation in human and animals and can be considered as a major genetic component of phenotypic diversity. Here we used the Illumina PorcineSNP60 BeadChip V2 and a DLY [Duroc × (Large White × Landrace)] commercial hybrid population to identify 272 CNVs belonging to 165 CNV regions (CNVRs), of which 66 are new. As CNVRs are specific to origin of population, our DLY-specific data is an important complementary to the existing CNV map in the pig genome. Eight CNVRs were selected for validation by quantitative real-time PCR (qRT-PCR) and the accurate rate was high (87.25%). Gene function analysis suggested that a common CNVR may play an important role in multiple traits, including growth rate and carcass quality.     

High-throughput genotyping in citrus accessions using an SNP genotyping array

We developed a 384 multiplexed SNP array, named CitSGA-1, for the genotyping of Citrus cultivars, and evaluated the performance and reliability of the genotyping. SNPs were surveyed by direct sequence comparison of the sequence tagged site (STS) fragment amplified from genomic DNA of cultivars representing the genetic diversity of citrus breeding in Japan. Among 1497 SNPs candidates, 384 SNPs for a high-throughput genotyping array were selected based on physical parameters of Illumina’s bead array criteria. The assay using CitSGA-1 was applied to a hybrid population of 88 progeny and 103 citrus accessions for breeding in Japan, which resulted in 73,726 SNP calls. A total of 351 SNPs (91 %) could call different genotypes among the DNA samples, resulting in a success rate for the assay comparable to previously reported rates for other plant species. To confirm the reliability of SNP genotype calls, parentage analysis was applied, and it indicated that the number of reliable SNPs and corresponding STSs were 276 and 213, respectively. The multiplexed SNP genotyping array reported here will be useful for the efficient construction of linkage map, for the detection of markers for marker-assisted breeding, and for the identification of cultivars.     

Genome-wide analyses of the Jeju,Thoroughbred, and Jeju crossbred horse populations using the high density SNP array

The Jeju horse is an indigenous Korean horse breed that is currently registered with the Food and Agriculture Organization of the United Nations. However, there is severe lack of genomic studies on Jeju horse. This study was conducted to investigate genetic characteristics of horses including Jeju horse, Thoroughbred and Jeju crossbred (Jeju?×?Thoroughbred) populations. We compared the genomes of three horse populations using the Equine SNP70 Beadchip array. Short-range Linkage disequilibrium was the highest in Thoroughbred, whereas r² values were lowest in Jeju horse. Expected heterozygosity was the highest in Jeju crossbred (0.351), followed by the Thoroughbred (0.337) and Jeju horse (0.311). The level of inbreeding was slightly higher in Thoroughbred (??0.009) than in Jeju crossbred (??0.035) and Jeju horse (??0.038). F_ST value was the highest between Jeju horse and Thoroughbred (0.113), whereas Jeju crossbred and Thoroughbred showed the lowest value (0.031). The genetic relationship was further assessed by principal component analysis, suggesting that Jeju crossbred is more genetically similar to Thoroughbred than Jeju horse population. Additionally, we detected potential selection signatures, for example, in loci located on LCORL/NCAPG and PROP1 genes that are known to influence body. Genome-wide analyses of the three horse populations showed that all the breeds had somewhat a low level of inbreeding within each population. In the population structure analysis, we found that Jeju crossbred was genetically closer to Thoroughbred than Jeju horse. Furthermore, we identified several signatures of selection which might be associated with traits of interest. To our current knowledge, this study is the first genomic research, analyzing genetic relationships of Jeju horse, Thoroughbred and Jeju crossbred.     

Longitudinal analysis and genotyping of infant dominant bifidobacterial populations

Bifidobacterial population dynamics were investigated by the longitudinal analysis of the dominant population isolated from the feces of young infants. After molecular identification and fingerprinting comparison, clone identity of the consecutive strains belonging to the same species for one individual was performed by pulsed-field gel electrophoresis. The results, obtained from 15 individuals sampled four times over a five-week period suggested a turnover of the dominant bifidobacteria in the population not only at the species but also at its species representative levels. This study provides new insights of the in vivo dynamics of commensal bifidobacteria. It highlights the need to take into consideration the fluctuation of bifidobacterial populations that may occur in one individual in order to investigate reliably the impact of dietary components, such as probiotics or prebiotics, on the intestinal ecosystem.     

Phylogenetic relationships among the European and American bison and seven cattle breeds reconstructed using the BovineSNP50 Illumina Genotyping BeadChip

Here we present the first attempt to use the BovineSNP50 Illumina Genotyping BeadChip for genome-wide screening of European bison Bison bonasus bonasus (EB), two subspecies of American bison: the plains bison Bison bison bison (PB), the wood bison Bison bison athabascae (WB) and seven cattle Bos taurus breeds. Our aims were to (1) reconstruct their evolutionary relationships, (2) detect any genetic signature of past bottlenecks and to quantify the consequences of bottlenecks on the genetic distances amongst bison subspecies and cattle, and (3) detect loci under positive or stabilizing selection. A Bayesian clustering procedure (STRUCTURE) detected ten genetically distinct clusters, with separation among all seven cattle breeds and European and American bison, but no separation between plain and wood bison. A linkage disequilibrium based program (LDNE) was used to estimate the effective population size (N _e) for the cattle breeds; N _e was generally low, relative to the census size of the breeds (cattle breeds: mean N _e = 299.5, min N _e = 18.1, max N _e = 755.0). BOTTLENECK 1.2 detected signs of population bottlenecks in EB, PB and WB populations (sign test and standardized sign test: p = 0.0001). Evidence for loci under selection was found in cattle but not in bison. All extant wild populations of bison have shown to have survived severe bottlenecks, which has likely had large effects on genetic diversity within and differentiation among groups.     

Genome-wide copy number variation in Hanwoo, Black Angus, and Holstein cattle

Hanwoo, Korean native cattle, is indigenous to the Korean peninsula. They have been used mainly as draft animals for about 5,000 years; however, in the last 30 years, their main role has been changed to meat production by selective breeding which has led to substantial increases in their productivity. Massively parallel sequencing technology has recently made possible the systematic identification of structural variations in cattle genomes. In particular, copy number variation (CNV) has been recognized as an important genetic variation complementary to single-nucleotide polymorphisms that can be used to account for variations of economically important traits in cattle. Here we report genome-wide copy number variation regions (CNVRs) in Hanwoo cattle obtained by comparing the whole genome sequence of Hanwoo with Black Angus and Holstein sequence datasets. We identified 1,173 and 963 putative CNVRs representing 16.7 and 7.8 Mbp from comparisons between Black Angus and Hanwoo and between Holstein and Hanwoo, respectively. The potential functional roles of the CNVRs were assessed by Gene Ontology enrichment analysis. The results showed that response to stimulus, immune system process, and cellular component organization were highly enriched in the genic-CNVRs that overlapped with annotated cattle genes. Of the 11 CNVRs that were selected for validation by quantitative real-time PCR, 9 exhibited the expected copy number differences. The results reported in this study show that genome-wide CNVs were detected successfully using massively parallel sequencing technology. The CNVs may be a valuable resource for further studies to correlate CNVs and economically important traits in cattle.     

Identification of novel single nucleotide polymorphisms (SNPs) in deer (Odocoileus spp.) using the BovineSNP50 BeadChip

Single nucleotide polymorphisms (SNPs) are growing in popularity as a genetic marker for investigating evolutionary processes. A panel of SNPs is often developed by comparing large quantities of DNA sequence data across multiple individuals to identify polymorphic sites. For non-model species, this is particularly difficult, as performing the necessary large-scale genomic sequencing often exceeds the resources available for the project. In this study, we trial the Bovine SNP50 BeadChip developed in cattle (Bos taurus) for identifying polymorphic SNPs in cervids Odocoileus hemionus (mule deer and black-tailed deer) and O. virginianus (white-tailed deer) in the Pacific Northwest. We found that 38.7% of loci could be genotyped, of which 5% (n = 1068) were polymorphic. Of these 1068 polymorphic SNPs, a mixture of putatively neutral loci (n = 878) and loci under selection (n = 190) were identified with the F(ST)-outlier method. A range of population genetic analyses were implemented using these SNPs and a panel of 10 microsatellite loci. The three types of deer could readily be distinguished with both the SNP and microsatellite datasets. This study demonstrates that commercially developed SNP chips are a viable means of SNP discovery for non-model organisms, even when used between very distantly related species (the Bovidae and Cervidae families diverged some 25.1-30.1 million years before present).     

Determination of genomic breakpoints in an epileptic patient using genotyping array

Recent advances in DNA microarray technology have enabled the identification of small alterations throughout the genome. We used standard karyotype analysis, followed by DNA microarray analysis and PCR to precisely map the chromosomal 4p deletion and determine the deletion breakpoints in the genome of an epileptic patient. The karyotype of the patient was 46,XY,del(4)(p15.2p15.3) as determined by G-banding analysis. We used a high-density oligonucleotide genotyping array to estimate the size of the deletion (4.5 Mb) and to locate the breakpoints within a 9-kb region on one side of the deletion and a 100-kb region on the other side. We amplified by PCR and sequenced the genomic region encompassing the breakpoints, and mapped the deletion to regions extending from 21648457 to 26164287 and from 26164505 to 26167493, respectively (chromosome 4 of NCBI Homo sapiens Genome Build 35.1). The deletion involves 18 genes, one of which (CCKAR) is partially deleted.     

Error rate for imputation from the Illumina BovineSNP50 chip to the Illumina BovineHD chip

Background

Imputation of genotypes from low-density to higher density chips is a cost-effective method to obtain high-density genotypes for many animals, based on genotypes of only a relatively small subset of animals (reference population) on the high-density chip. Several factors influence the accuracy of imputation and our objective was to investigate the effects of the size of the reference population used for imputation and of the imputation method used and its parameters. Imputation of genotypes was carried out from 50 000 (moderate-density) to 777 000 (high-density) SNPs (single nucleotide polymorphisms).

Methods

The effect of reference population size was studied in two datasets: one with 548 and one with 1289 Holstein animals, genotyped with the Illumina BovineHD chip (777 k SNPs). A third dataset included the 548 animals genotyped with the 777 k SNP chip and 2200 animals genotyped with the Illumina BovineSNP50 chip. In each dataset, 60 animals were chosen as validation animals, for which all high-density genotypes were masked, except for the Illumina BovineSNP50 markers. Imputation was studied in a subset of six chromosomes, using the imputation software programs Beagle and DAGPHASE.

Results

Imputation with DAGPHASE and Beagle resulted in 1.91% and 0.87% allelic imputation error rates in the dataset with 548 high-density genotypes, when scale and shift parameters were 2.0 and 0.1, and 1.0 and 0.0, respectively. When Beagle was used alone, the imputation error rate was 0.67%. If the information obtained by Beagle was subsequently used in DAGPHASE, imputation error rates were slightly higher (0.71%). When 2200 moderate-density genotypes were added and Beagle was used alone, imputation error rates were slightly lower (0.64%). The least imputation errors were obtained with Beagle in the reference set with 1289 high-density genotypes (0.41%).

Conclusions

For imputation of genotypes from the 50 k to the 777 k SNP chip, Beagle gave the lowest allelic imputation error rates. Imputation error rates decreased with increasing size of the reference population. For applications for which computing time is limiting, DAGPHASE using information from Beagle can be considered as an alternative, since it reduces computation time and increases imputation error rates only slightly.     

Selective genotyping for QTL detection using sib pair analysis in outbred populations with hierarchical structures

A simulation study illustrates the effects of the inclusion of half-sib pairs as well as the effects of selective genotyping on the power of detection and the parameter estimates in a sib pair analysis of data from an outbred population. The power of QTL detection obtained from samples of sib pairs selected according to their within family variance or according to the mean within family variance within half sib family was compared and contrasted with the power obtained when only full sib pair analysis was used. There was an increase in power (4–16%) and decrease in the bias of parameter estimates with the use of half-sib information. These improvements in power and parameter estimates depended on the number of the half sib pairs (half sib family size). Almost the same power as that obtained using all the available sib pairs could be achieved by selecting only 50–60% the animals. The most effective method was to select both full and half sib pairs on the basis of high within full sib family variance for the trait in question. The QTL position estimates were in general slightly biased towards the center of the chromosome and the QTL variance estimates were biased upwards, there being quite large differences in bias depending on the selection method.     

Development,validation and genetic analysis of a large soybean SNP genotyping array

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under‐use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high‐density genetic mapping and genome‐wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome‐wide association analyses. More than four million high‐quality SNPs identified from high‐depth genome re‐sequencing of 16 soybean accessions and low‐depth genome re‐sequencing of 31 soybean accessions were used to select 180 961 SNPs for creation of the Axiom^® SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170 223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene‐rich chromosomal regions suggest that this array may be suitable for genome‐wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.     

Genome-wide two-locus epistasis scans in prostate cancer using two European populations

Approximately 40 single nucleotide polymorphisms (SNPs) that are associated with prostate cancer (PCa) risk have been identified through genome-wide association studies (GWAS). However, these GWAS-identified PCa risk-associated SNPs can explain only a small proportion of heritability (~13%) of PCa risk. Gene-gene interaction is speculated to be one of the major factors contributing to the so-called missing heritability. To evaluate the gene-gene interaction and PCa risk, we performed a two-stage genome-wide gene-gene interaction scan using a novel statistical approach named "Boolean Operation-based Screening and Testing". In the first stage, we exhaustively evaluated all pairs of SNP-SNP interactions for ~500,000 SNPs in 1,176 PCa cases and 1,101 control subjects from the National Cancer Institute Cancer Genetic Markers of Susceptibility (CGEMS) study. No SNP-SNP interaction reached a genome-wide significant level of 4.4E-13. The second stage of the study involved evaluation of the top 1,325 pairs of SNP-SNP interactions (P(interaction) <1.0E-08) implicated in CGEMS in another GWAS population of 1,964 PCa cases from the Johns Hopkins Hospital (JHH) and 3,172 control subjects from the Illumina iControl database. Sixteen pairs of SNP-SNP interactions were significant in the JHH population at a P(interaction) cutoff of 0.01. However, none of the 16 pairs of SNP-SNP interactions were significant after adjusting for multiple tests. The current study represents one of the first attempts to explore the high-dimensional etiology of PCa on a genome-wide scale. Our results suggested a list of SNP-SNP interactions that can be followed in other replication studies.     

A reliable multiplex genotyping assay for HCV using a suspension bead array

The genotyping of the hepatitis C virus (HCV) plays an important role in the treatment of HCV because genotype determination has recently been incorporated into the treatment guidelines for HCV infections. Most current genotyping methods are unable to detect mixed genotypes from two or more HCV infections. We therefore developed a multiplex genotyping assay to determine HCV genotypes using a bead array. Synthetic plasmids, genotype panels and standards were used to verify the target‐specific primer (TSP) design in the assay, and the results indicated that discrimination efforts using 10 TSPs in a single reaction were extremely successful. Thirty‐five specimens were then tested to evaluate the assay performance, and the results were highly consistent with those of direct sequencing, supporting the reliability of the assay. Moreover, the results from samples with mixed HCV genotypes revealed that the method is capable of detecting two different genotypes within a sample. Furthermore, the specificity evaluation results suggested that the assay could correctly identify HCV in HCV/human immunodeficiency virus (HIV) co‐infected patients. This genotyping platform enables the simultaneous detection and identification of more than one genotype in a same sample and is able to test 96 samples simultaneously. It could therefore provide a rapid, efficient and reliable method of determining HCV genotypes in the future.     

Non-invasive transgenic mouse genotyping using stool analysis

Commonly applied genotyping of transgenic mice involves using tail or ear biopsies which may cause discomfort to the animal. We tested the possibility of polymerase chain reaction (PCR)-based mouse genotyping using stool specimens from three transgenic mouse lines that overexpress 10-18 transgene copies of human keratin polypeptide 18, as compared to genotyping using tail biopsies. Stool specimens were obtained with ease and provided easy detection of the human transgene product. The method was also able to detect endogenous mouse actin and keratin genes which presumably are present at two copies each. Nested PCR was not necessary for genotyping using stool-derived genomic material but did increase the relative magnitude of the signal obtained. The non-invasive genotyping method described herein offers a reproducible, sensitive and effective modality that could replace invasive tissue sampling procedures currently used to test thousands of genetically altered mice.     

Identification of BoLA-DRB3.2 alleles in Korean native cattle (Hanwoo) and Holstein populations using a next generation sequencer

Bovine leucocyte antigen (encoded by BoLA) has been widely studied to identify the association with many traits related to immunity. Exon2 of BoLA-DRB3 is extremely polymorphic, and more than 100 alleles have been identified. We investigated polymorphisms of BoLA-DRB3.2 in Korean native cattle and Holstein populations using a next generation sequencer of the GS-FLX Titanium system. We found 38 alleles including 11 new alleles (BoLA-DRB3*1303, *4702, *7101, *7501, *7201, *7301, *7601, *1104, *7701, *7401 and *50021) in Hanwoo, and nine alleles including one new allele (BoLA-DRB3*7601) in Holstein. The 454 sequencing method is a promising alternative technology for high throughput genotyping of BoLA-DRB3.2 because of its technical advantages that allow it to overcome the disadvantages of sequence-based typing methods.