Identification and analysis of error types in high-throughput genotyping

Although it is clear that errors in genotyping data can lead to severe errors in linkage analysis, there is as yet no consensus strategy for identification of genotyping errors. Strategies include comparison of duplicate samples, independent calling of alleles, and Mendelian-inheritance-error checking. This study aimed to develop a better understanding of error types associated with microsatellite genotyping, as a first step toward development of a rational error-detection strategy. Two microsatellite marker sets (a commercial genomewide set and a custom-designed fine-resolution mapping set) were used to generate 118,420 and 22,500 initial genotypes and 10,088 and 8,328 duplicates, respectively. Mendelian-inheritance errors were identified by PedManager software, and concordance was determined for the duplicate samples. Concordance checking identifies only human errors, whereas Mendelian-inheritance-error checking is capable of detection of additional errors, such as mutations and null alleles. Neither strategy is able to detect all errors. Inheritance checking of the commercial marker data identified that the results contained 0.13% human errors and 0.12% other errors (0.25% total error), whereas concordance checking found 0.16% human errors. Similarly, Mendelian-inheritance-error checking of the custom-set data identified 1.37% errors, compared with 2.38% human errors identified by concordance checking. A greater variety of error types were detected by Mendelian-inheritance-error checking than by duplication of samples or by independent reanalysis of gels. These data suggest that Mendelian-inheritance-error checking is a worthwhile strategy for both types of genotyping data, whereas fine-mapping studies benefit more from concordance checking than do studies using commercial marker data. Maximization of error identification increases the likelihood of linkage when complex diseases are analyzed.     

SNPs in stress-responsive rice genes: validation, genotyping, functional relevance and population structure

ABSTRACT: BACKGROUND: Single nucleotide polymorphism (SNP) validation and large-scale genotyping are required to maximize the use of DNA sequence variation and determine the functional relevance of candidate genes for complex stress tolerance traits through genetic association in rice. We used the bead array platform-based Illumina GoldenGate assay to validate and genotype SNPs in a select set of stress-responsive genes to understand their functional relevance and study the population structure in rice. RESULTS: Of the 384 putative SNPs assayed, we successfully validated and genotyped 362 (94.3%). Of these 325 (84.6%) showed polymorphism among the 91 rice genotypes examined. Physical distribution, degree of allele sharing, admixtures and introgression, and amino acid replacement of SNPs in 263 abiotic and 62 biotic stress-responsive genes provided clues for identification and targeted mapping of trait-associated genomic regions. We assessed the functional and adaptive significance of validated SNPs in a set of contrasting drought tolerant upland and sensitive lowland rice genotypes by correlating their allelic variation with amino acid sequence alterations in catalytic domains and three-dimensional secondary protein structure encoded by stress-responsive genes. We found a strong genetic association among SNPs in the nine stress-responsive genes with upland and lowland ecological adaptation. Higher nucleotide diversity was observed in indica accessions compared with other rice sub-populations based on different population genetic parameters. The inferred ancestry of 16% among rice genotypes was derived from admixed populations with the maximum between upland aus and wild Oryza species. CONCLUSIONS: SNPs validated in biotic and abiotic stress-responsive rice genes can be used in association analyses to identify candidate genes and develop functional markers for stress tolerance in rice.     

Construction and validation of a prototype microarray for efficient and high-throughput genotyping of angiosperms

Until recently, the identification of plants relied on conventional techniques, such as morphological, anatomical and chemical profiling, that are often inefficient or unfeasible in certain situations. Extensive literature exists describing the use of polymerase chain reaction (PCR) DNA-based identification techniques, which offer a reliable platform, but their broad application is often limited by a low throughput. However, hybridization-based microarray technology represents a rapid and high-throughput tool for genotype identification at a molecular level. Using an innovative technique, a 'Subtracted Diversity Array' (SDA) of 376 features was constructed from a pooled genomic DNA library of 49 angiosperm species, from which pooled non-angiosperm genomic DNA was subtracted. Although not the first use of a subtraction technique for genotyping, the SDA method was superior in accuracy, sensitivity and efficiency, and showed high-throughput capacity and broad application. The SDA technique was validated for potential genotyping use, and the results indicated a successful subtraction of non-angiosperm DNA. Statistical analysis of the polymorphic features from the pilot study enabled the establishment of accurate phylogenetic relationships, confirming the potential use of the SDA technique for genotyping. Further, the technique substantially enriched the presence of polymorphic sequences; 68% were polymorphic when using the array to differentiate six angiosperm clades (Asterids, Rosids, Caryophyllids, Ranunculids, Monocots and Eumagnoliids). The 'proof of concept' experiments demonstrate the potential of establishing a highly informative, reliable and high-throughput microarray-based technique for novel application to sequence independent genotyping of major angiosperm clades.     

A geminivirus-induced gene silencing system for gene function validation in cassava

We have constructed an African cassava mosaic virus (ACMV) based gene-silencing vector as a reverse genetics tool for gene function analysis in cassava. The vector carrying a fragment from the Nicotiana tabacumsulfur gene (su), encoding one unit of the chloroplast enzyme magnesium chelatase, was used to induce the silencing of the cassava orthologous gene resulting in yellow–white spots characteristic of the inhibition of su expression. This result suggests that well developed sequence databases from model plants like Arabidopsis thaliana, Nicotiana benthamiana, N. tabacum, Lycopersicon esculentum and others could be used as a major source of information and sequences for functional genomics in cassava. Furthermore, a fragment of the cassava CYP79D2endogenous gene, sharing 89% homology with CYP79D1endogenous gene was inserted into the ACMV vector. The resultant vector was inducing the down regulation of the expression of these two genes which catalyze the first-dedicated step in the synthesis of linamarin, the major cyanogenic glycoside in cassava. At 21 days post-inoculation (dpi), a 76% reduction of linamarin content was observed in silenced leaves. Using transgenic plants expressing antisense RNA of CYP79D1and CYP79D2, Siritunga and Sayre (2003) obtained several lines with a reduction level varying from 60% to 94%. This result provides the first example of direct comparison of the efficiency of a virus-induced gene silencing (VIGS) system and the transgenic approach for suppression of a biosynthetic pathway. The ACMV VIGS system will certainly be a complement and in some cases an alternative to the transgenic approach, for gene discovery and gene function analysis in cassava.     

鸭脂联素基因单核苷酸多态性检测及群体遗传分析

以昆山麻鸭、樱桃谷鸭、高邮鸭、荆江麻鸭、金定鸭, 山麻鸭、龙白鸭和白羽番鸭8个鸭品种为实验材料, 根据鸭脂联素基因开放阅读框序列设计5对引物, 用PCR-SSCP方法进行单核苷酸多态性分析,并对不同品种群体进行群体遗传学分析。结果发现引物4扩增片段上共存有7个单碱基突变, 第430、457、523处的G-A 、A-G、T-C单碱基突变导致第144、153、175个氨基酸分别由丙氨酸（A）变为苏氨酸（T）、异亮氨酸（I）变为缬氨酸(V)、酪氨酸(Y)变为组氨酸(H); 而C507T, T540C, C576T和C597T 4个单碱基突变为沉默突变。鸭群中表现出AA、AB、AC、BB、BC、CC、DD、DE 8种基因型。8种基因型在8个鸭品种间分布存在极显著的差异(P<0.01)。除金定鸭外, 其他品种均处于Hardy-Weinberg平衡状态。群体遗传分析表明金定鸭的纯合度最高, 高邮鸭最低, 其他各群体的纯合度较相近; 金定鸭为低度多态, 高邮鸭为高度多态, 其他品种为中度多态。表明鸭脂联素基因不同品种中具有丰富的单核苷酸多态性, 可以进一步作为候选基因来分析其与脂肪性状的相关性。     

A high-throughput SNuPE assay for genotyping SNPs in the flanking regions of Zea mays sequence tagged simple sequence repeats

Information on genetic relationships among individuals is of importanceto maize breeders for line and hybrid development. Estimates on the geneticsimilarity of breeding materials is best obtained using DNA markers. SingleNucleotide Polymorphisms and small insertions/deletions are both emerging as anew generation of markers, due to their abundance and amenability to fullyautomated genotyping. Application of SNPs, for example in genetic mappingprojects or breeding programs, involves the analysis of a large number ofsamples, and therefore requires rapid, inexpensive, and highly automatedmethodsto genotype the sequence variants. Towards this, we have applied a highthroughput Single Nucleotide Primer Extension assay to assess 23 polymorphicbase variations at five microsatellite loci in 22 inbred maize lines, as wellasin a mapping population of two of the inbred lines. Using a MegaBACE automatedgenotyper, we are able to assay more than 1248 (96 × 13) samples in asingle 75 minute run. The SNuPE method successfully genotyped the base pairvariations of interest in all maize lines. It was also found that primerscontaining polymorphisms themselves could be used to genotype the samples. Thistechnique allows the rapid production of valuable high throughput informationongenetic relationships among maize varieties. We further demonstrate the utilityof this method, using it to successfully map one of the microsatellite loci.     

Identification of a large dataset of SNPs in Larimichthys polyactis using high-throughput 2b-RAD sequencing

The small yellow croaker (Larimichthys polyactis) is one of the most commercially exploited marine fishes along the coast of the Yellow-Bohai Sea and the East China Sea. In this study, we used next-generation high-throughput 2b-RAD-seq technology to identify novel SNPs in L. polyactis. We scored a total of 1 374 008 putative SNPs genome-wide. Further filtration yielded a final dataset of 6457 high-quality SNPs. These SNP markers presented sufficient power in detecting genetic distinction between two wild samples from the Yellow-Bohai Sea and one captive sample from the East China Sea, and inbreeding reflected by strong heterozygote deficiency within the samples; some of the genetic polymorphisms are diagnostic among the samples and related to biological functions. The panel of SNPs could be used as powerful tools in further population genetic and stock assessment research in L. polyactis as it has relatively scarce genomic resources. The findings from this study will advance our understanding of population and functional genomics for facilitating fishery resource management and developing desirable characteristics for the benefit of culture and farming of L. polyactis.     

New insights into the Lake Chad Basin population structure revealed by high-throughput genotyping of mitochondrial DNA coding SNPs

Background

Located in the Sudan belt, the Chad Basin forms a remarkable ecosystem, where several unique agricultural and pastoral techniques have been developed. Both from an archaeological and a genetic point of view, this region has been interpreted to be the center of a bidirectional corridor connecting West and East Africa, as well as a meeting point for populations coming from North Africa through the Saharan desert.

Methodology/Principal Findings

Samples from twelve ethnic groups from the Chad Basin (n = 542) have been high-throughput genotyped for 230 coding region mitochondrial DNA (mtDNA) Single Nucleotide Polymorphisms (mtSNPs) using Matrix-Assisted Laser Desorption/Ionization Time-Of-Flight (MALDI-TOF) mass spectrometry. This set of mtSNPs allowed for much better phylogenetic resolution than previous studies of this geographic region, enabling new insights into its population history. Notable haplogroup (hg) heterogeneity has been observed in the Chad Basin mirroring the different demographic histories of these ethnic groups. As estimated using a Bayesian framework, nomadic populations showed negative growth which was not always correlated to their estimated effective population sizes. Nomads also showed lower diversity values than sedentary groups.

Conclusions/Significance

Compared to sedentary population, nomads showed signals of stronger genetic drift occurring in their ancestral populations. These populations, however, retained more haplotype diversity in their hypervariable segments I (HVS-I), but not their mtSNPs, suggesting a more ancestral ethnogenesis. Whereas the nomadic population showed a higher Mediterranean influence signaled mainly by sub-lineages of M1, R0, U6, and U5, the other populations showed a more consistent sub-Saharan pattern. Although lifestyle may have an influence on diversity patterns and hg composition, analysis of molecular variance has not identified these differences. The present study indicates that analysis of mtSNPs at high resolution could be a fast and extensive approach for screening variation in population studies where labor-intensive techniques such as entire genome sequencing remain unfeasible.     

Identification of high utility SNPs for population assignment and traceability purposes in the pig using high-throughput sequencing

The objectives of this study were to develop breed-specific single nucleotide polymorphisms (SNPs) in five pig breeds sequenced with Illumina's Genome Analyzer and to investigate their usefulness for breed assignment purposes. DNA pools were prepared for Duroc, Landrace, Large White, Pietrain and Wild Boar. The total number of animals used for sequencing was 153. SNP discovery was performed by aligning the filtered reads against Build 7 of the pig genome. A total of 313,964 high confidence SNPs were identified and analysed for the presence of breed-specific SNPs (defined in this context as SNPs for which one of the alleles was detected in only one breed). There were 29,146 putative breed-specific SNPs identified, of which 4441 were included in the PorcineSNP60 beadchip. Upon re-examining the genotypes obtained using the beadchip, 193 SNPs were confirmed as being breed specific. These 193 SNPs were subsequently used to assign an additional 490 individuals from the same breeds, using the sequenced individuals as reference populations. In total, four breed assignment tests were performed. Results showed that for all methods tested 99% of the animals were correctly assigned, with an average probability of assignment of at least 99.2%, indicating the high utility of breed-specific markers for breed assignment and traceability. This study provides a blueprint for the way next-generation sequencing technologies can be used for the identification of breed-specific SNPs, as well as evidence that these SNPs may be a powerful tool for breed assignment and traceability of animal products to their breeds of origin.     

Identification and validation of QTL for Sclerotinia midstalk rot resistance in sunflower by selective genotyping

Midstalk rot, caused by Sclerotinia sclerotiorum (Lib.) de Bary, is an important cause of yield loss in sunflower (Helianthus annuus L.). Objectives of this study were to: (1) estimate the number, genomic positions and genetic effects of quantitative trait loci (QTL) for resistance to midstalk rot in line TUB-5-3234, derived from an interspecific cross; (2) determine congruency of QTL between this line and other sources of resistance; and (3) make inferences about the efficiency of selective genotyping (SG) in detecting QTL conferring midstalk rot resistance in sunflower. Phenotypic data for three resistance (stem lesion, leaf lesion and speed of fungal growth) and two morphological (leaf length and leaf length with petiole) traits were obtained from 434 F₃ families from cross CM625 (susceptible) × TUB-5-3234 (resistant) under artificial infection in field experiments across two environments. The SG was applied by choosing the 60 most resistant and the 60 most susceptible F₃ families for stem lesion. For genotyping of the respective F₂ plants, 78 simple sequence repeat markers were used. Genotypic variances were highly significant for all traits. Heritabilities and genotypic correlations between resistance traits were moderate to high. Three to four putative QTL were detected for each resistance trait explaining between 40.8% and 72.7% of the genotypic variance ( ). Two QTL for stem lesion showed large genetic effects and corroborated earlier findings from the cross NDBLOS_sel (resistant) × CM625 (susceptible). Our results suggest that SG can be efficiently used for QTL detection and the analysis of congruency for resistance genes across populations.     

Single-track sequencing for genotyping of multiple SNPs in the N-acetyltransferase 1 (NAT1) gene

Background  

Fast, cheap and reliable methods are needed to identify large populations, which may be at risk in relation to environmental exposure. Polymorphisms in NAT1 (N-acetyl transferase) may be suitable markers to identify individuals at risk.     

Identification and expression analysis of MinD gene involved in plastid division in cassava

Cassava is a tropical crop known for its starchy root and excellent properties. Considering that starch biosynthesis in the amyloplast is affected by its division, it appears conceivable that the regulation of plastid division plays an important role in starch accumulation. As a member of the Min system genes, MinD participated in the spatial regulation of the position of the plastid division site.In our studies, sequence analysis and phylogenetic analysis showed that MeMinD has been highly conserved during the evolutionary process. Subcellular localisation indicated that MeMinD carries a chloroplast transit peptide and was localised in the chloroplast. Overexpression of MeMinD resulted in division site misplacement and filamentous formation in E. coli, indicating that MeMinD protein was functional across species. MeMinD exhibited different spatial and temporal expression patterns which was highly expressed in the source compared to that in the sink organ.     

Identification and in silico analysis of functional SNPs of the BRCA1 gene

Single-nucleotide polymorphisms (SNPs) play a major role in the understanding of the genetic basis of many complex human diseases. Also, the genetics of human phenotype variation could be understood by knowing the functions of these SNPs. It is still a major challenge to identify the functional SNPs in a disease-related gene. In this work, we have analyzed the genetic variation that can alter the expression and the function of the BRCA1 gene using computational methods. Of the total 477 SNPs, 65 were found to be nonsynonymous (ns) SNPs. Among the 14 SNPs in the untranslated region, 4 were found in the 5' and 10 were found in the 3' untranslated region (UTR). It was found that 16.9% of the nsSNPs were damaging, by both the SIFT and the PolyPhen servers. The UTR Resource tool suggested that 2 of 4 SNPs in the 5' UTR and 3 of 10 SNPs in the 3' UTR might change the protein expression levels. We identified major mutations from proline to serine at positions 1776 and 1812 of the native protein of the BRCA1 gene. From a comparison of the stabilizing residues of the native and mutant proteins, we propose that an nsSNP (rs1800751) could be an important candidate for the breast cancer caused by the BRCA1 gene.     

Primer design for PCR and sequencing in high-throughput analysis of SNPs

To achieve high-throughput analysis of allele frequencies in human SNPs, we have developed automated methodsfor designing PCR and DNA sequencing primers. We found we could run the PCR assays at quite stringent, uniform conditions. The design process used freely available databases, including dbSNP, SNPper, and TSC, and publicly available software including RepeatMasker and Primer3. We describe parameters for the software and other considerations that increase experimental success. As anticipated. some assays filed at the design stage due primarily to the genomic locations of repetitive sequences, extreme GC content regions, or lack of sufficient sequence. However, over 23,000 assays, including 96% of those recently analyzed, have been experimentally successfuL Similar design methods could be usedfor PCR assays in any organism with substantial available sequence.     

Identification of three SNPs in the porcine myostatin gene (MSTN)

Thirteen pairs of primers were designed for the entire porcine MSTN gene to enable PCR amplification for the detection of single nucleotide polymorphisms (SNPs) by a PCR-SSCP approach. Altogether 96.5% (1089/1128) of the encoding regions and 971 bp of the non-coding regions were screened. A total of three polymorphisms were identified with PCR-SSCP. They were located in the promoter, intron one and exon three regions of the gene. These polymorphisms were then confirmed to be point mutations (T --> A transversion, G --> A transition and C --> T transition respectively) by sequencing. Allele frequencies were determined for all three SNPs in several different pig breed populations. The polymorphisms were found to be rare in Western breeds, but much more common in Chinese breeds. Whether they have any relationship with the marked difference in lean meat mass between Western and Chinese breeds requires further study.     

Selecting additional tag SNPs for tolerating missing data in genotyping

Background  

Recent studies have shown that the patterns of linkage disequilibrium observed in human populations have a block-like structure, and a small subset of SNPs (called tag SNPs) is sufficient to distinguish each pair of haplotype patterns in the block. In reality, some tag SNPs may be missing, and we may fail to distinguish two distinct haplotypes due to the ambiguity caused by missing data.     

A high-throughput method for genotyping S-RNase alleles in apple

We present a new efficient screening tool for detection of S-alleles in apple. The protocol using general and multiplexed primers for PCR reaction and fragment detection on an automatized capillary DNA sequencer exposed a higher number of alleles than any previous studies. Analysis of alleles is made on basis of three individual fragment sizes making the allele interpretation highly accurate. The method was employed to genotype 432 Malus accessions and exposed 25 different S-alleles in a selection of Malus domestica cultivars of mainly Danish origin (402 accessions) as well as a selection of other Malus species (30 accessions). The allele S3 (28 %) was the most common among the Danish cultivars followed by S1 and S7 (both 27 %). The alleles S36 and S40 not previously reported from M. domestica were found in 6 and 17 cultivars, respectively. Complete allelic composition was found in 91 % of the 369 diploid accessions and in 86 % of the 63 triploids concerned. We further identified a relatively high frequency of S33 and S34, which has not been considered by most previous studies. The protocol presented here is easy to adopt and saves both time and work effort compared to previous methods. The robustness is illustrated by the great accuracy and a high number of S-alleles presented.