An approach to high-throughput genotyping

Limitations of methods for measuring the concentration of human genomic DNA and oligonucleotide samples

We compared different methods (absorbance, fluorescent dye-binding, and digital PCR) for measuring the concentrations of human genomic DNA from cultured cells and absorbance measurements of a synthetic DNA oligonucleotide. NIST Standard Reference Material (SRM) 2082, a pathlength absorbance standard, was used to benchmark the absorbance measurements done with microvolume spectrophotometers and a microvolume plate reader. Control absorbance values were measured on a high accuracy spectrophotometer and a NIST calibrated pathlength cuvette. Measurements of the human genomic DNA sample were done with several types of fluorescent dye binding assays using different DNA calibrators. The fluorescent dye binding methods gave different results for genomic DNA depending upon the type of DNA calibrator and the fluorescent dye that was used. The human genomic DNA sample was also characterized by using six different droplet digital PCR assays (amplicons located on different chromosomes) to measure the average copy number. Conversion of the digital PCR data to copy numbers was sensitive to the droplet size used for calculations and conversion to mass concentration was dependent upon the molecular weight of the human genome used for the calculations. The results from the different methods were compared and the caveats for each measurement method were discussed.

A genome-wide scan points to a susceptibility locus for bipolar disorder on chromosome 12

Our previous results pointed to a putative gene for susceptibility to bipolar affective disorder located on the chromosomal region 12q23-q24 that segregated in the Saguenay-Lac-St-Jean population of Quebec. We report here results from a second genome-wide scan based on the analysis of 380 polymorphic microsatellite markers. For the purpose of this analysis, an additional 18 families were recruited from the Saguenay-Lac-St-Jean region and pooled to our previous sample to improve its statistical power, giving a total of 394 sampled individuals. This work confirms the presence of a susceptibility locus for affective disorder on chromosome 12q24 with parametric LOD score value of 3.35 at D12S378 when pedigrees were broken into nuclear families and analysed under a recessive segregation model. This result was supported by neighbouring markers and by a LOD score value of 5.05 at D12S378 under model-free analysis. Other regions of lower interest were indicated on chromosomes 2, 5, 7, 9, 10, 17 and 20.

The Family Investigation of Nephropathy and Diabetes (FIND): design and methods

BACKGROUND

The Family Investigation of Nephropathy and Diabetes (FIND) is a multicenter study designed to identify genetic determinants of diabetic nephropathy. It is conducted in eight U.S. clinical centers and a coordinating center, and with four ethnic groups (European Americans, African Americans, Mexican Americans, and American Indians). Two strategies are used to localize susceptibility genes: a family-based linkage study and a case-control study using mapping by admixture linkage disequilibrium (MALD).

METHODS

In the family-based study, probands with diabetic nephropathy are recruited with their parents and selected siblings. Linkage analyses will be conducted to identify chromosomal regions containing genes that influence the development of diabetic nephropathy or related quantitative traits such as serum creatinine concentration, urinary albumin excretion, and plasma glucose concentrations. Regions showing evidence of linkage will be examined further with both genetic linkage and association studies to identify genes that influence diabetic nephropathy or related traits. Two types of MALD studies are being done. One is a case-control study of unrelated individuals of Mexican American heritage in which both cases and controls have diabetes, but only the case has nephropathy. The other is a case-control study of African American patients with nephropathy (cases) and their spouses (controls) unaffected by diabetes and nephropathy; offspring are genotyped when available to provide haplotype data.

CONCLUSIONS

Identification of genes that influence susceptibility to diabetic nephropathy will lead to a better understanding of how nephropathy develops. This should eventually lead to improved treatment and prevention.

A robust, low- to medium-throughput prnp genotyping system in sheep

BACKGROUND

In many countries breeding programs for resistance to scrapie in sheep are established. Therefore, the demand on genotyping capacities of the polymorphisms of the prion protein gene (prnp) relevant to presently known disease associations and EU regulations is steadily increasing. Most published typing methods are not well suited for routine typing of large sample numbers in smaller service laboratories for different reasons: they require partly manual data processing, sophisticated and sensitive protocols, high efforts regarding time and manpower, multiple step reactions or substantial hardware investments. To overcome these drawbacks, we developed a prnp typing method that is based on a 'multiplex amplification refractory mutation system' (ARMS) reaction.

METHODS

In this study we combined the amplification refractory mutation system (ARMS) with standard fluorescent based fragment length analyses method to develop a prnp genotyping method (PRNP ARMS).

RESULTS

By optimised primer design it was possible to type the 4 relevant single nucleotide polymorphisms (SNPs) in the prnp simultaneously in one multiplex reaction. Automated fragment length analysis enabled automated allele designation. Suitability of the PRNP ARMS for routine application was proven by typing samples with known genotypes and larger sample numbers from half-sib families.

CONCLUSION

The ARMS PRNP typing method established in this study is universally suited for a broad range of typing projects with different requirements. It provides an efficient and inexpensive diagnostic mutation analysis that will improve the quality of prnp genotyping compared with other low-cost methods. It can be implemented by most molecular genetic laboratories using standard equipment.

A novel method for automatic genotyping of microsatellite markers based on parametric pattern recognition

Genetic mapping of loci affecting complex phenotypes in human and other organisms is presently being conducted on a very large scale, using either microsatellite or single nucleotide polymorphism (SNP) markers and by partly automated methods. A critical step in this process is the conversion of the instrument output into genotypes, both a time-consuming and error prone procedure. Errors made during this calling of genotypes will dramatically reduce the ability to map the location of loci underlying a phenotype. Accurate methods for automatic genotype calling are therefore important. Here, we describe novel algorithms for automatic calling of microsatellite genotypes using parametric pattern recognition. The analysis of microsatellite data is complicated both by the occurrence of stutter bands, which arise from Taq polymerase misreading the number of repeats, and additional bands derived form the non-template dependent addition of a nucleotide to the 3' end of the PCR products. These problems, together with the fact that the lengths of two alleles in a heterozygous individual may differ by only two nucleotides, complicate the development of an automated process. The novel algorithms markedly reduce the need for manual editing and the frequency of miscalls, and compares very favourably with commercially available software for automatic microsatellite genotyping.

Undetected genotyping errors cause apparent overtransmission of common alleles in the transmission/disequilibrium test

The transmission/disequilibrium test (TDT), a family-based test of linkage and association, is a popular and intuitive statistical test for studies of complex inheritance, as it is nonparametric and robust to population stratification. We carried out a literature search and located 79 significant TDT-derived associations between a microsatellite marker allele and a disease. Among these, there were 31 (39%) in which the most common allele was found to exhibit distorted transmission to affected offspring, implying that the allele may be associated with either susceptibility to or protection from a disease. In 27 of these 31 studies (87%), the most common allele appeared to be overtransmitted to affected offspring (a risk factor), and, in the remaining 4 studies, the most common allele appeared to be undertransmitted (a protective factor). In a second literature search, we identified 92 case-control studies in which a microsatellite marker allele was found to have significantly different frequencies in case and control groups. Of these, there were 37 instances (40%) in which the most common allele was involved. In 12 of these 37 studies (32%), the most common allele was enriched in cases relative to controls (a risk factor), and, in the remaining 25 studies, the most common allele was enriched in controls (a protective factor). Thus, the most common allele appears to be a risk factor when identified through the TDT, and it appears to be protective when identified through case-control analysis. To understand this phenomenon, we incorporated an error model into the calculation of the TDT statistic. We show that undetected genotyping error can cause apparent transmission distortion at markers with alleles of unequal frequency. We demonstrate that this distortion is in the direction of overtransmission for common alleles. Therefore, we conclude that undetected genotyping errors may be contributing to an inflated false-positive rate among reported TDT-derived associations and that genotyping fidelity must be increased.

High throughput genotyping technologies for pharmacogenomics

Genetic differences between individuals play a role in determining susceptibility to diseases as well as in drug response. The challenge today is first to discover the range of genetic variability in the human population and then to define the particular gene variants, or alleles, that contribute to clinically important outcomes. Consequently, high throughput, automated methods are being developed that allow rapid scoring of microsatellite alleles and single nucleotide polymorphisms (SNPs). Many detection technologies are being used to accomplish this goal, including electrophoresis, standard fluorescence, fluorescence polarization, fluorescence resonance energy transfer, and mass spectrometry. SNP alleles may be distinguished by any one of several methods, including single nucleotide primer extension, allele-specific hybridization, allele-specific primer extension, oligonucleotide ligation assay, and invasive signal amplification. Newer methods require less sample manipulation, increase sensitivity, allow more flexibility, and decrease reagent costs. Recent developments show promise for continuing these trends by combining amplification and detection steps and providing flexible, miniaturized platforms for genotyping.

Toward high-throughput genotyping: dynamic and automatic software for manipulating large-scale genotype data using fluorescently labeled dinucleotide markers

To efficiently manipulate large amounts of genotype data generated with fluorescently labeled dinucleotide markers, we developed a Microsoft database management system, named. offers several advantages. First, it accommodates the dynamic nature of the accumulations of genotype data during the genotyping process; some data need to be confirmed or replaced by repeat lab procedures. By using, the raw genotype data can be imported easily and continuously and incorporated into the database during the genotyping process that may continue over an extended period of time in large projects. Second, almost all of the procedures are automatic, including autocomparison of the raw data read by different technicians from the same gel, autoadjustment among the allele fragment-size data from cross-runs or cross-platforms, autobinning of alleles, and autocompilation of genotype data for suitable programs to perform inheritance check in pedigrees. Third, provides functions to track electrophoresis gel files to locate gel or sample sources for any resultant genotype data, which is extremely helpful for double-checking consistency of raw and final data and for directing repeat experiments. In addition, the user-friendly graphic interface of renders processing of large amounts of data much less labor-intensive. Furthermore, has built-in mechanisms to detect some genotyping errors and to assess the quality of genotype data that then are summarized in the statistic reports automatically generated by. The can easily handle >500,000 genotype data entries, a number more than sufficient for typical whole-genome linkage studies. The modules and programs we developed for the can be extended to other database platforms, such as Microsoft SQL server, if the capability to handle still greater quantities of genotype data simultaneously is desired.

A schizophrenia-susceptibility locus at 6q25, in one of the world's largest reported pedigrees

We have completed a genome scan of a 12-generation, 3,400-member pedigree with schizophrenia. Samples from 210 individuals were collected from the pedigree. We performed an "affecteds-only" genome-scan analysis using 43 members of the pedigree. The affected individuals included 29 patients with schizophrenia, 10 with schizoaffective disorders, and 4 with psychosis not otherwise specified. Two sets of white-European allele frequencies were used-one from a Swedish control population (46 unrelated individuals) and one from the pedigree (210 individuals). All analyses pointed to the same region: D6S264, located at 6q25.2, showed a maximum LOD score of 3.45 when allele frequencies in the Swedish control population were used, compared with a maximum LOD score of 2.59 when the pedigree's allele frequencies were used. We analyzed additional markers in the 6q25 region and found a maximum LOD score of 6.6 with marker D6S253, as well as a 6-cM haplotype (markers D6S253-D6S264) that segregated, after 12 generations, with the majority of the affected individuals. Multipoint analysis was performed with the markers in the 6q25 region, and a maximum LOD score of 7.7 was obtained. To evaluate the significance of the genome scan, we simulated the complete analysis under the assumption of no linkage. The results showed that a LOD score >2.2 should be considered as suggestive of linkage, whereas a LOD score >3.7 should be considered as significant. These results suggest that a common ancestral region was inherited by the affected individuals in this large pedigree.

High prevalence of the very rare Wilson disease gene mutation Leu708Pro in the Island of Gran Canaria (Canary Islands, Spain): a genetic and clinical study

The molecular basis of Wilson disease (WD), an autosomal recessive disorder, is the presence of mutations in the ATP7B gene, a copper transporting ATPase. Hospital records indicated a higher prevalence of WD (1 in 2,600) in some counties in the northeastern region of the island of Gran Canaria (Canary Islands, Spain) that was around 10-fold higher than that described for European populations (1 in 30,000). The ATP7B gene was analyzed for mutations in 24 affected subjects, revealing a high prevalence of the rare Leu708Pro mutation present in 12 homozygous and 7 heterozygous individuals. In these patients, who constitute one of the largest described cohorts of WD homozygotes, we found a variable clinical presentation of the disease, although the biochemical picture was homogenous and characteristic, thereby confirming that the Leu708Pro change is indeed a mutation associated with WD. Haplotype analysis of subjects homozygous for the Leu708Pro mutation showed a conserved shared region smaller than 1 centimorgan (cM), and the region of linkage disequilibrium between the Leu708Pro mutation and neighboring microsatellite markers extended approximately 4.6 cM. When comparing the amount of linkage disequilibrium versus genetic distance from the disease mutation, it was estimated that a common ancestral Leu708Pro chromosome may have been introduced in Gran Canaria over 56 generations ago, dating it back to pre-Hispanic times. The prevalence, and the tight geographical distribution of the Leu708Pro chromosome suggests that the Canary Islands can be considered a genetic isolate for linkage disequilibrium studies.

Lack of evidence for linkage to chromosomes 13 and 8 for schizophrenia and schizoaffective disorder

A previous report [Blouin et al., 1998: Nat Genet 20:70-73] suggesting linkage to chromosomes 13q32 and 8p21 in families with schizophrenia led us to investigate these regions in a large set of 301 multiplex families with schizophrenia. Multipoint analyses failed to reveal evidence for linkage to any portion of chromosome 13, while only a weakly positive score was present on 8p using the identical marker reported in the earlier report. Failure to confirm the Blouin et al claims in a substantially larger cohort adds emphasis to the inconsistency of the findings concerning linkage in schizophrenia. Am. J. Med. Genet. (Neuropsychiatr. Genet.) 96:235-239, 2000.

Using quality measures to facilitate allele calling in high-throughput genotyping

Currently, the main limitation in high-throughput microsatellite genotyping is the required manual editing of allele calls. Even though programs for automated allele calling have been available for several years, they have limited capability because accurate data could only be assured by manual inspection of the electropherograms for confirmation. Here we describe the development of a parametric approach to allele call quality control that eliminates much of the time required for manual editing of the data. This approach was implemented in an editing tool, Decode-GT, that works downstream of the allele calling program, TrueAllele (TA). Decode-GT reads the output data from TA, displays the underlying electropherograms for the genotypes, and sorts the allele calls into three categories: good, bad, and ambiguous. It discards the bad calls, accepts the good calls, and suggests that the user inspect the ambiguous calls, thereby reducing dependence on manual editing. For the categorization we use the following parameters: (1) the quality value for each allele call from TrueAllele; (2) the peak height of the alleles; and (3) the size of the peak shift needed to move peaks into the nearest bin. Here we report how we optimized the parameters such that the size of the ambiguous category was minimized, and both the number of miscalled genotypes in the good category and the useable genotypes in the bad category were negligible. This approach reduces the manual editing time and results in <1% miscalls.

Loss of Heterozygosity Analysis Using Whole Genome Amplification, Cell Sorting, and Fluorescence-Based PCR

Loss of heterozygosity (LOH) is a common genetic lesion found in many human neoplasms. Extending investigation of LOH to large-scale clinical and public health science studies has proven difficult because of the small size and cellular and genetic heterogeneity of human neoplasms, in addition to the challenges associated with increasing throughput. Our approach to LOH analysis was developed using clinical biopsy samples from patients with Barrett’s esophagus (BE) and uses flow cytometric cell sorting to increase sample purity, whole genome amplification to increase sample amount, and automated fluorescent genotyping to increase sample throughput. This approach allows LOH assessment at 20 loci in DNA extracted from 1000 flow-purified cells while maintaining accurate and reproducible allele ratios compared with the standard method of using genomic DNA. This method of analysis should allow accurate, reproducible determination of allele ratios in a variety of human tumors and premalignant conditions.

Characterization of SYBR Gold nucleic acid gel stain: a dye optimized for use with 300-nm ultraviolet transilluminators

The highest sensitivity nucleic acid gel stains developed to date are optimally excited using short-wavelength ultraviolet or visible light. This is a disadvantage for laboratories equipped only with 306- or 312-nm UV transilluminators. We have developed a new unsymmetrical cyanine dye that overcomes this problem. This new dye, SYBR Gold nucleic acid gel stain, has two fluorescence excitation maxima when bound to DNA, one centered at approximately 300 nm and one at approximately 495 nm. We found that when used with 300-nm transillumination and Polaroid black-and-white photography, SYBR Gold stain is more sensitive than ethidium bromide, SYBR Green I stain, and SYBR Green II stain for detecting double-stranded DNA, single-stranded DNA, and RNA. SYBR Gold stain's superior sensitivity is due to the high fluorescence quantum yield of the dye-nucleic acid complexes ( approximately 0.7), the dye's large fluorescence enhancement upon binding to nucleic acids ( approximately 1000-fold), and its capacity to more fully penetrate gels than do the SYBR Green gel stains. We found that SYBR Gold stain is as sensitive as silver staining for detecting DNA-with a single-step staining procedure. Finally, we found that staining nucleic acids with SYBR Gold stain does not interfere with subsequent molecular biology protocols.

Quantitative-trait locus for specific language and reading deficits on chromosome 6p

Reading disability (RD), or dyslexia, is a complex cognitive disorder manifested by difficulties in learning to read, in otherwise normal individuals. Individuals with RD manifest deficits in several reading and language skills. Previous research has suggested the existence of a quantitative-trait locus (QTL) for RD on the short arm of chromosome 6. In the present study, RD subjects' performance in several measures of word recognition and component skills of orthographic coding, phonological decoding, and phoneme awareness were individually subjected to QTL analysis, with a new sample of 126 sib pairs, by means of a multipoint mapping method and eight informative DNA markers on chromosome 6 (D6S461, D6S276, D6S105, D6S306, D6S258, D6S439, D6S291, and D6S1019). The results indicate significant linkage across a distance of at least 5 cM for deficits in orthographic (LOD = 3.10) and phonological (LOD = 2.42) skills, confirming previous findings.

A genome-wide search for schizophrenia susceptibility genes

We completed a systematic genome-wide search for evidence of loci linked to schizophrenia using a collection of 70 pedigrees containing multiple affected individuals according to three phenotype classifications: schizophrenia only (48 pedigrees; 70 sib-pairs); schizophrenia plus schizoaffective disorder (70 pedigrees; 101 sib-pairs); and a broad category consisting of schizophrenia, schizoaffective disorder, paranoid or schizotypal personality disorder, psychosis not otherwise specified (NOS), delusional disorder, and brief reactive psychosis (70 pedigrees; 111 sib-pairs). All 70 families contained at least one individual affected with chronic schizophrenia according to DSM-III-R criteria. Three hundred and thirty-eight markers spanning the genome were typed in all pedigrees for an average resolution of 10.5 cM (range, 0-31 cM) and an average heterozygosity of 74.3% per marker. The data were analyzed using multipoint nonparametric allele-sharing and traditional two-point lod score analyses using dominant and recessive, affecteds-only models. Twelve chromosomes (1, 2, 4, 5, 8, 10, 11, 12, 13, 14, 16, and 22) had at least one region with a nominal P value <0.05, and two of these chromosomes had a nominal P value <0.01 (chromosomes 13 and 16), using allele-sharing tests in GENEHUNTER. Five chromosomes (1, 2, 4, 11, and 13) had at least one marker with a lod score >2.0, allowing for heterogeneity. These regions will be saturated with additional markers and investigated in a new, larger set of families to test for replication.

Linkage of bipolar affective disorder to chromosome 18 markers in a new pedigree series

Several groups have reported evidence suggesting linkage of bipolar affective disorder (BPAD) to chromosome 18. We have reported data from 28 pedigrees that showed linkage to marker loci on 18p and to loci 40 cM distant on 18q. Most of the linkage evidence derived from families with affected phenotypes in only the paternal lineage and from marker alleles transmitted on the paternal chromosome. We now report results from a series of 30 new pedigrees (259 individuals) genotyped for 13 polymorphic markers spanning chromosome 18. Subjects were interviewed by a psychiatrist and were diagnosed by highly reliable methods. Genotypes were generated with automated technology and were scored blind to phenotype. Affected sib pairs showed excess allele sharing at the 18q markers D18S541 and D18S38. A parent-of-origin effect was observed, but it was not consistently paternal. No robust evidence of linkage was detected for markers elsewhere on chromosome 18. Multipoint nonparametric linkage analysis in the new sample combined with the original sample of families supports linkage on chromosome 18q, but the susceptibility gene is not well localized.

Emerging technologies from the Human Genome Project for understanding susceptibility and risk

The new technologies from the Human Genome Program provide exceptional opportunities for surveying and measuring human exposure, as well as determining susceptibility on an individual-by-individual basis. These new technologies will soon enable rapid screening of populations at risk, as well as the broader public, for a variety of genes known to be associated with increased risk. These include specific oncogenes, tumor suppressor genes and DNA repair enzymes. Use of these technologies also presents a number of ethical issues, both in screening and in use of the information about individuals. Overall, the use of rapid genotyping technologies will introduce a specificity and possible group identifiers that will present new challenges to the determination of risk within the EPA mandate.

A simple method for automated allele binning in microsatellite markers

High-throughput fluorescent genotyping requires a considerable amount of automation for accurate and efficient processing of genetic markers. Automated DNA sequencers and corresponding software products are commercially available that contribute substantially to increased throughput rates for large-scale genotyping projects. However, some conceptually simple tasks still require time-consuming manual intervention that imposes bottlenecks on throughput capacity. One of these tasks is the conversion of imprecise DNA fragment sizes determined by commercial software programs to the underlying discrete alleles that the sizes represent. Here we describe a simple method for assigning allele sizes into their appropriate allele "bins" using least-squares minimization procedures. The method requires no special treatment of family data on plates, internal/external size standards, or electropherogram data manipulation. Tests of the method using the ABI 373A automated DNA sequencer and accompanying Genescan/Genotyper software resulted in accurate automatic classification of all alleles in >80% of 208 markers analyzed, with the remaining 20% being appropriately identified as requiring additional attention to laboratory conditions. Specific characteristics of different markers, including differences in PCR product size and inexact repeat lengths (e.g., 1. 9 bp for a dinucleotide repeat), are accommodated by the method and their properties discussed.