Search tips
Search criteria 


Logo of tswjThe Scientific World Journal
ScientificWorldJournal. 2012; 2012: 365104.
Published online 2012 June 18. doi:  10.1100/2012/365104
PMCID: PMC3385616

Mixed Sequence Reader: A Program for Analyzing DNA Sequences with Heterozygous Base Calling


The direct sequencing of PCR products generates heterozygous base-calling fluorescence chromatograms that are useful for identifying single-nucleotide polymorphisms (SNPs), insertion-deletions (indels), short tandem repeats (STRs), and paralogous genes. Indels and STRs can be easily detected using the currently available Indelligent or ShiftDetector programs, which do not search reference sequences. However, the detection of other genomic variants remains a challenge due to the lack of appropriate tools for heterozygous base-calling fluorescence chromatogram data analysis. In this study, we developed a free web-based program, Mixed Sequence Reader (MSR), which can directly analyze heterozygous base-calling fluorescence chromatogram data in .abi file format using comparisons with reference sequences. The heterozygous sequences are identified as two distinct sequences and aligned with reference sequences. Our results showed that MSR may be used to (i) physically locate indel and STR sequences and determine STR copy number by searching NCBI reference sequences; (ii) predict combinations of microsatellite patterns using the Federal Bureau of Investigation Combined DNA Index System (CODIS); (iii) determine human papilloma virus (HPV) genotypes by searching current viral databases in cases of double infections; (iv) estimate the copy number of paralogous genes, such as β-defensin 4 (DEFB4) and its paralog HSPDP3.

1. Introduction

The detection of genomic variations is important in studying the relationships between causative genes and diseases and the relationships between predisposing genes and complex trait diseases, such as type 2 diabetes, coronary heart disease, and cancers [14]. Structural genomic variations also provide important information about both genetic diversity and human evolution [5]. Human genomic variations include single-nucleotide polymorphisms (SNPs), variable number of tandem repeats (VNTRs), short tandem repeats (STRs, or microsatellites), and copy number variations (CNVs) [6]. Among these genomic variants, there are currently 51,810,853 reference SNPs for the human genome, which include 6,516,668 indel sequences and 5,214 microsatellite markers, according to dbSNP Build 135. For instance, a recent genome sequence of a human individual revealed 292,102 heterozygous indel events and 559,473 homozygous indels [7].

Genomic variants are frequently identified with heterozygous base-calling fluorescence chromatogram data generated from the direct sequencing of genomic PCR products using the dye-terminator method with Applied Biosystems (ABIs) autosequencers, such as models 3700 or 3730. Several groups have developed programs to analyze heterozygous base-calling fluorescence chromatogram data. For instance, ShiftDetector is a program for detecting shift mutations and calculating a probability score for sequences reconstructed from  .abi sequence files [8]. The Indelligent program uses dynamic programming optimization to decode the heterozygous indels with International Union of Pure and Applied Chemistry DNA code system (IUPAC code) [9, 10]. The CHILD (CHromatogram In/Del Location and Detection) program was specifically designed to detect indels in DNA mixtures where one variant is rare, and it can also estimate the ratio of two variants [11]. Generally speaking, all of the currently available programs can be applied for the analysis of indel genomic variants even without reference sequences. However, the alignment of heterozygous base-calling fluorescent data with a reference database can be used to detect the physical position of indel within the genome. Even so, some heterozygous indels may not be easily visualized (Figures 1(b-4) and 2(a)).

Figure 1
Workflow and user interfaces (input and result output) of the Mixed Sequences Reader (MSR). A flowchart describing the MSR workflow is shown in the left-side panel (a); the input interface of MSR is shown in the right-side panel (b); the result output ...
Figure 2
Experimental confirmation of the Indel identified by Mixed Sequence Reader. A 9-bp insertion at chromosome 7:55249011 was detected by Mixed Sequence Reader from a directly sequenced PCR product. (a) PCR direct sequencing chromatography trace, (b) MSR ...

Short tandem repeats (STRs) are the most important markers in forensic genetic analysis, and several commercial kits for STR analysis are available [12, 13]. Two-nucleotide repeats are the most prevalent STRs, while repeats with more nucleotides (n = 3~6) are less common in the human genome. Among the more than 5000 human STR markers, the Federal Bureau of Investigation (FBI) has included thirteen loci in its Combined DNA Index System (CODIS) database, which contains information from more than 6,384,379 individuals [1416]. Current STR genotyping uses multiplex PCR with fluorescent STR primers to amplify genomic regions containing VNTRs or STRs, and the PCR products are then separated with capillary electrophoresis, which distinguishes fragments by length but does not display the actual sequence [17]. Nevertheless, sequencing methods remain very useful for analysis of STRs because they reveal the actual sequences of STRs [18], although the number may be too ambiguous to interpret (Figure 3(c)). To resolve this ambiguity, DNA cloning is often required to identify the different sequences of the two alleles. Currently, no program is able to analyze microsatellite repeat units (or CODIS) directly using the chromatogram trace data, even though the chromatogram trace profiles of heterozygous microsatellites are similar to those of heterozygous indels with 2 or more nucleotide deletions (Supplementary Figure  1 in Supplementary Materials available online at doi:10.1100/2012/365104).

Figure 3
Experimental confirmation of one CODIS-STR locus identified by the Mixed Sequencer Reader. (a) The FBI CODIS Core STR Loci map (adapted from The D13S317 locus (chromosome 13:55144219-55144352 region) was ...

Some genotypes of HPV are oncogenic. Therefore, routine screening and genotyping of HPV in women are crucial for prevention of cervical cancer. Current HPV detection methods are either using PCR techniques or HPV genotyping array system [19]. About 90% of cervical cancer tissues are infected by a single HPV genotype, while about 10% of specimens by two or more types of HPV. In the single HPV-infected samples, the viral genotype can be easily identified by genotyping array or PCR. For those with mixed infection by double types of HPV, mixed chromatogram traces are observed, in a similar way to those of SNP and indel sequences (Figure 4(b)).

Figure 4
Application of MSR in identifying double-infection of HPV in a cervical cancer sample. (a) A sample infected with two strains of HPV was selected according to HPV array results. The HPV sequences in the specimen were amplified with a pair of PCR primers ...

CNVs, another type of genomic variation, are segments of DNA with variable copy numbers; the length of a single CNV may range from one kilobase to several megabases [20]. One well-known CNV occurs in the β-defensin 4 locus (DEFB4), which is known to have copy numbers that range from two to seven. The CNV of the defensin genes is associated with increased susceptibility to infectious diseases, autoimmune diseases, inflammatory disorders, Crohn's disease, and certain cancers [21, 22]. Array-based comparative genomic hybridization (aCGH) is used to detect genomic CNVs. Other methods for CNV validation include multiplex amplifiable probe hybridization combined with restriction enzyme digest variant ratios (MAPH/REDVRs) [23, 24], multiplex ligation-dependent probe amplification (MLPA) [25, 26], paralog ratio test (PRT) [24, 27, 28], and real-time polymerase chain reaction [29]. The PRT method uses a pair of primers to amplify two paralogous sequences, which are then separated with capillary electrophoresis. By detecting the ratio between the chromatographic intensities of the two PCR products, PRT can estimate the copy number of the CNV [24, 27, 28]. However, the chromatogram traces of the PCR products derived from paralogous genes are often too heterozygous to be analyzed with any currently available programs (Figure 5(b)).

Figure 5
Experimental confirmation of HSPDP3 and DEFB4 copy numbers identified by Mixed Sequence Reader. (a) The primers used to amplify the DEF4B and HSPDP3 genes in the human genome. The LRi ratio between DEFB4 and HSPDP3 was calculated for each chromatography ...

To address the above issues, we have developed a program, Mixed Sequencer Reader (MSR), which can be used to identify indels, microsatellite copy numbers, and CODIS combinations. On the basis of the Indelligent method, the heterozygous sequences are identified as two distinct sequences, which are further aligned with reference genomic sequences to provide information about the physical position of indels, copy number or types of STR, and paralogous genes. We also applied this program to identify double HPV infections in cervical cancer tissues and estimate the copy number of paralogous human genes (e.g., the β-defensin 4 gene DEFB4 and its paralog HSPDP3). The software is freely available at

2. Materials and Methods

2.1. Mixed Sequence Reader (MSR)

The MSR program was developed to detect heterogeneity in chromatographic traces, determine the physical positions of the detected variants in the human genome, and identify the type of genomic variation present. The algorithm used in MSR was modified from that of Indelligent, but MSR is designed to use reference database alignment. The analytic steps used in MSR are described below and shown in Figure 1(a).

2.1.1. Importing the DNA Sequence Chromatographs and Selecting the Reference Database

The imported files are chromatography traces in the  .abi format. The base peak positions, quality values, and four channel signals (A, C, G, T) recorded in the  .abi file are extracted and analyzed to identify the major and the minor signals at each base location (Figures 1(b-1) and 1(b-4)). Users can select the desired reference database (Figure 1(b-2)).

2.1.2. Defining the Log Ratio of Intensity (LRi)

The LRi value for each base position is defined as the log ratio of the chromatographic intensities of two combined sequences. The formula for LRi is:

LRi=log2(major fluorescence intensityminor fluorescence intensity).

If the sequence position contained only one major band (no heterogeneity in the chromatographic trace), the value of LRi should be infinite (log 2(1/0)). If the DNA sequences contain two heterogeneous chromatographic traces with equal intensities, the LRi should be 0 (log 2(1/1)). After the  .abi file is imported, users can define an LRi cut-off value (Figure 1(b-3)), or the MSR program can automatically set the cut-off value using smooth LRi (Figure 1(c-1), green line). For example, in the results shown in Figure 1(c-1), the shift of signal intensity at 212 bp was detected by the MSR program and considered an “indel”-type heterozygosity. The LRi values of “indel” sequences are higher than the cutoff value for sequences without heterogeneous fluorescence chromatography traces, so the LRi line drops when heterogeneous fluorescence traces are identified. When there is no obvious shift of the LRi line, the sequences might be a “mixed” type, and the LRi cutoff value is automatically set at 2.0 by the MSR program. The sequences of either “indel” or “mixed” heterozygosities are then displayed in IUPAC code (Figure 1(c-3)).

2.1.3. Decomposing the IUPAC Code Using the Indelligent Algorithm

In the Indelligent algorithm, dynamic programming is used to convert the IUPAC code into two nucleotides (i.e., M is converted into A/C, W into A/T, Y into C/T, K into G/T, and S into G/C) [9, 10]. For ambiguous bases that cannot be decomposed with Indelligent, major and minor sequences are assigned according to the intensities of the corresponding fluorescence signals (Figure 1(c-4)).

2.1.4. Finding the DNA Sequence of the Most Possible Variances (MPVs) by BLAST

The major and minor DNA sequences derived from Section 2.1.3 were further analyzed by BLASTN against the databases that were built into the MSR program (GRCh37 primary reference assembly or CODIS, HPV, DEFB4/HSPDP3 reference sequences) or user-defined reference sequences (Figure 1(b-2)) to detect the most possible variances (Figure 1(c-5)). The major and minor sequences were BLASTed against reference sequences to obtain the major and minor MPVs. Once the physical positions of the major and minor MPVs were identified, the sequences were categorized as indels (Figure 1(c-6)) or considered as a “mixed” type. For “mixed” type heterogeneities, MSR continues with the following procedures.

2.1.5. Deriving the Optimal Combinations

Each top M (major) MPV and top N (minor) MPV are combined pairwise into IUPAC code, resulting in (M × N) combinations. Each IUPAC code combination is aligned to the original signal IUPAC code (Figure 1(c-3)), and the (M × N) IUPAC code combinations with the highest scores are identified as the optimal MPV combinations by the MSR program (Figure 1(c-7)). Using the selected reference database, such as the default MSR human genome sequence database or a user-imported database, the MSR program identifies the genotypes of the mixed sequences.

2.1.6. Calculating the Ratio between the Major and Minor Sequences

For all heterogeneous bases, the LRi values were calculated as described in Section 2.1.2. MSR calculated the medium LRi value of heterogeneous nucleotides as the “Sequence Mix Ratio” (example shown in Figure 4(d)). The Sequence Mix Ratio is proportional to the signal ratio between the major and minor sequences.

2.2. DNA Cloning

To validate DNA sequences, PCR products amplified from the experimental samples were selected for cloning experiments to confirm the variations. The PCR products were cloned into the pCRII-TOPO cloning system (Invitrogen) according to the manufacturer's recommendations. The cloned DNA sequences were then analyzed with an ABI 3770 autosequencer.

2.3. Determination of the β-Defensin 4 (DEFB4) Copy Number

All DNA samples used in this study were unlinked from clinical information, and the DNA collection was approved by the Institutional Review Board (IRB) of Chang Gung Memorial Hospital (CGMH) (#99-0229B, IRB#100-2900A3). Genomic DNA samples from 100 normal individuals were tested for CNVs of the paralogous DEFB4/HSPDP3 genes. DEFB4 and its paralog, the HSPDP3 pseudogene, were amplified from 50 genomic DNA specimens. PCR using a previously reported pair of primers [27] amplifies two products of similar size, one from the DEFB4 gene on chromosome 8 and the other from the HSPDP3 pseudogene on chromosome 5 (the copy number of HSPDP3 is always 2). These PCR products were directly sequenced (see Figure 5(b) for an example). For each peak in the chromatogram trace, we calculated the LRi value, which indicated the DEFB4/HSPDP3 ratio. For each sample, the median values of all heterogeneous sequences were used to determine the copy number of the DEFB4 gene. The k-means method was used to partition the different copy numbers.

3. Results

3.1. The Mixed Sequences Reader Interface

The data import interface of the Mixed Sequence Reader is divided into three parts (Figure 1(b)): (1) data import, (2) reference database selection, and (3) MSR parameter settings. At the data import step, users can import an  .abi file into MSR (Figure 1(b-1)). Users can also test the performance of MSR using 260 sample  .abi files (Figure 1(b-1)), some of which were experimentally validated in this study. After the  .abi files are imported, users can preview the chromatography data by moving the cursor in the right side panel (Figure 1(b-4)). Users can then either select the reference database to use (GRCh37, CODIS, HPV, DEFB4/HSPDP3) or import their own reference sequences in FASTA format (Figure 1(b-2)). Then, users can define the LRi cutoff value, sequence type and specify the ignored head and tail sequence lengths (Figure 1(b-3)). Then, users click the “Run” button to execute the MSR program. The results of a sample analysis are shown in Figure 1(c).

3.2. Indel and Short Tandem Repeat Sequences

We first analyzed indel sequences with MSR (Figure 2(b)). From the NCBI dbSNP database, we selected 6 indel sites of high heterozygosity, PCR amplified the corresponding genomic DNA, and directly sequenced the PCR products to obtain  .abi files. Fourteen  .abi files were successfully analyzed by MSR (Supplementary Table  1). To validate the MSR results, the predicted indel sequences were confirmed by cloning the PCR products, followed by sequencing of the plasmid DNAs. For example, one heterozygous chromatography trace was predicted to contain a nine nucleotide insertion at human chr7:55249011 (Figure 2(b)). Although similar results were also obtained when the sequences were analyzed with Indelligent or other programs, the MSR was able to identify the physical position of the indel by BLASTing the sequences against the GRCh37 database (Figure 2(b)). The experimental results confirmed the results predicted by MSR (Figure 2(a)).

We also selected several simple short tandem repeats (STRs) from the NCBI dbSNP database. Three STR-containing regions of the genome were amplified by PCR, sequenced, and analyzed with MSR (13  .abi files of STR in Supplementary Table  2). Among the 13  .abi files analyzed, the PCR products corresponding to 2  .abi files were further cloned and validated by single-plasmid DNA sequencing. The number of each repeats were also defined by MSR (Supplementary Figure  1).

3.3. The Repeat Structure of a Short Tandem Repeat in the CODIS Database

Because CODIS is the largest STR database currently available, the repeat patterns of the 13 loci in CODIS were documented (Figure 3(a)) [14, 15]. The D13S317 locus ([TATC]7–15) was amplified by PCR using specific primers (Supplementary Table  3), and the resulting PCR products were directly sequenced (Figures 3(b) and 3(c)). According to the MSR results, the D13S317 locus genotype in the tested individual should be [TATC]8/[TATC]9 (Figure 3(d)). The PCR products were further cloned and sequenced, and the results confirmed the MSR prediction (Figure 3(e)). Similar results were obtained when the same  .abi file was analyzed using the Indelligent program.

Some STRs in the CODIS database contain complex repeat structures, such as the D21S11 locus with the 65-repeat structure [TCTA]4–11[TCTG]3–14[TCTA]0–3[TA]0-1 [TCTA]3TCA[TCTA]2 TCCATA [TCTA]6–15. The alignment of the two sequences in the sample 211-CODIS-D21S11-3  .ab1 file detected two 4-bp gap structures (Supplementary Figure  2). Because of the complexity of the CODIS STR repeat structure, its pattern was not easily solved using the Indelligent program without a reference database. However, the optimal STR repeat structure was successfully identified by the MSR program using the reference databases. A total of 29 sequences for 10 CODIS STR sites were analyzed with MSR, and the results are shown in Supplementary Table  4.

3.4. Infection with Two Genotypes of HPV in Cervical Cancer Specimens

Another type of “mixed” sequence is represented by viral coinfection in the same specimen. For example, double HPV genotypes within one sample could be analyzed by EasyChip using probes designed to amplify the variable regions that are unique to each HPV genotype [19]. To validate the microarray data, dually infected specimens were amplified at the L1 region of HPV and the resulting PCR products were sequenced directly (Figure 4(c)). The alignment of sequences representing the different genotypes of HPV revealed both indels and SNPs (Figures 4(b) and 4(c)). One of the samples was identified by MSR as coinfected with HPV-33 and HPV-81, compatible with the EasyChip results (Figure 4(e)). In addition, the ratio between the major (HPV-33) and minor (HPV-81) sequences calculated by MSR could be used to estimate the relative ratio of the two genotypes of HPV. The ratio between HPV-33 and HPV-81 was 3.4 : 1. Genotype-specific PCR confirmed the results of MSR prediction for all 7 specimens that were defined by EasyChip as infected with multiple HPV strains (Supplementary Table  5).

3.5. Copy Number Variations of Paralogous Genes

We also applied MSR to detect the copy number of the DEFB4 gene, which is well known for its multiple copy number variations. The PCR primers used amplified DEFB4 and its paralog, the pseudogene HSPDP3 (Figure 5(a)) [27]; thus, the chromatography traces generated by PCR direct sequencing were very heterozygous (Figure 5(b)) because the chromatographs contained sequences from both chromosome 5 (HSPDP3 as reference, copy number n = 2) and chromosome 8 (DEFB4, often variable copy numbers). The LRi value of each chromatography peak was used to estimate the ratio between DEFB4/HSPDP3 and calculate the DEFB4 copy number (Figure 5(a), right-side panel). The med(LRi) was defined as the median LRi value in all heterogeneous sequences. The scatter plot of the med(LRi) values of the DNA sequences amplified by the forward and reverse primers for 98 individual DNA samples is shown in Figure 5(d). The five groups shown in different colors were partitioned with the K-means clustering algorithm. To confirm these results, 5 specimens were randomly selected and the PCR products were cloned into the pCRII-TOPO cloning system for single-plasmid sequencing (Figure 5(d), circles; Supplementary Table  6). For each specimen, at least 20 clones were sequenced to calculate the HSPDP3 and DEFB4 ratio. All of the cloning-sequencing results were identical to the results predicted by MSR.

4. Discussion

Structural variations in the human genome are clinically important [6]. For instance, CNVs in DEFB4 are associated with susceptibility to infectious disease, autoimmune, inflammatory disorders, and even cancers [5, 20]. The copy numbers of the genes CCL3 L1, CCL4 L1, and TBC1D3 vary between 0 and 10 copies, and such variations have been found to be associated with susceptibility to HIV-1 [30, 31]. In addition, microsatellite markers are used as indicators for global genome stability and are especially useful in genomic research of cancer [32]. In short, structural variations shape the genome and determine disease susceptibility at the individual level. Therefore, analytical methods that can detect structural genomic variations are acutely needed to study the relationship between genomic variations and disease.

The MSR program introduced in this study can be used to directly analyze heterozygous base-calling chromatographs to detect multiple structural variations, including SNPs, microsatellites, and CNVs, in the human genome. The fluorescence intensity of chromatographs has already been used to detect SNPs [33], but we have extended this analysis of heterozygous base-calling chromatography to explore more structural variations in the human genome. The accuracy of the MSR predictions was validated by other methods. Our analyses show that MSR can also be used to identify double infections by different genotypes of human papilloma virus (HPV) in cervical cancer tissues [34, 35]. It is worth testing whether MSR can be used to determine the presence of multiple viral infections in other cancer tissues [34, 35]. The ability to identify dual viral infections is limited because the reference sequences for these viruses often mutate rapidly, especially for RNA viruses [36]. Therefore, to improve this program for the aforementioned applications, we would need more clinical specimens to challenge this software.

The MSR program does have some limitations, however. First, the majority of structural variations in the human genome, which can be readily identified with MSR, are small indels. On the other hand, deletions of large fragments are likely better analyzed by array-based comparative genomic hybridization methods (aCGH) [37]. Second, some of the heterozygosity identified by base-calling chromatography may be caused by the formation of DNA secondary structures (likely in the GC-rich or AT-rich regions) that result in band compression; these sequences are difficult to analyze with MSR. Third, the STRs in the human genome are not yet fully characterized. Therefore, the MSR may not be able to predict all STRs. In this study, we used some types of STRs as examples of the potential of MSR, but more samples are needed to analyze all 13 core loci in the CODIS database and other STRs. Fourth, MSR and all other web sources only provide tools to analyze or predict structural variations in the human genome. We should rely on experimental data to confirm these results. Fifth, the MSR program is designed to read only fluorescence chromatography tracings derived from the ABI 3730 autosequencer, but not for those generated by the recently developed ultrahigh-throughput sequencers, such as Roche 454, ABI SoLid, and Illumina/Solexa.

A comparison of the features of the MSR program and other currently available DNA sequence analysis software is summarized in Table 1. All of the available programs can detect indels. Most of the programs can read files in  .abi format, with the exception of Indelligent, which only processes IUPAC code data. Only MSR can map sequence data to a reference database and report the most possible mixed sequence. MSR and CHILD can estimate the ratio of sequences in a mixture. Only PolyScan can process next-generation sequencing data [38].

Table 1
Feature comparison between MSR and other currently available DNA sequencing software.

In conclusion, we have developed a user-friendly web-based program, Mixed Sequences Reader (MSR), to analyze heterozygous fluorescent chromatographs derived from an autosequencer. Using this program, several types of human genomic variations, including SNPs, indels, and CNVs of microsatellites or genes can be detected from a single DNA sequence read. Furthermore, MSR is useful for detecting viral infection with double genotypes in clinical specimens.

Supplementary Material

Supplementary Table 1: The filename, LRi cutoff value, and physical position in samples used for InDel analysis in this study.

Supplementary Table 2: The filename, LRi cutoff value and repeat structure in samples used for microsatellite analysis in this study.

Supplementary Table 3: The primer paris designed for amplification of 10 CODIS locus in this study.

Supplementary Table 4: The filename, LRi cutoff value and repeat structure in samples analyzing CODIS pattern in this study.

Supplementary Table 5: The filename, LRi cutoff value, and HPV genotype in dual infection cervical cancers samples in this study.

Supplementary Table 6: The copy number ratio of DEFB4/HSPDP3 predicted by MSR program and T/A cloning experiment.

Supplementary Figure 1: (A) Chromatography trace of one Indel [CA] STR repeat structure and those of plasmid cloning, (B) MSR analyzed results are shown as the IUPAC code and decomposed sequence.

Supplementary Figure 2: (A) Chromatography trace of a combined two CODIS D21S11 Locus subtype, (B) there were two 4-bp gaps in the alignment of two subtypes in CODIS D21S11 Locus.

Authors' Contribution

Chun-Tien Chang and Chi-Neu Tsai contributed equally to this paper.

Conflict of Interests

The authors declare that they have no conflict of interests.


This work was supported by the National Research Program for Genomic Medicine (NSC-97-3112-B-001-020 to Y.-S. Lee and C.-H. Chen); The National Science Council, Taiwan (NSC-98-2320-B-182-034-MY3 to C.-N. Tsai and NSC-100-2221-E-126-011-MY3 to C. Y. Tang); Chang-Gung Memorial Hospital (CMRPD190571 to C.-N. Tsai); and the Department of Health, Taiwan (DOH99-TD-I-111-TM013, DOH99-TD-C-111-006 to T.-H. Wang).


1. Janssens AC, van Duijn CM. Genome-based prediction of common diseases: advances and prospects. Human Molecular Genetics. 2008;17(2):R166–R173. [PubMed]
2. Manolio TA. Genomewide association studies and assessment of the risk of disease. The New England Journal of Medicine. 2010;363(2):166–176. [PubMed]
3. Menashe I, Maeder D, Garcia-Closas M, et al. Pathway analysis of breast cancer genome-wide association study highlights three pathways and one canonical signaling cascade. Cancer Research. 2010;70(11):4453–4459. [PMC free article] [PubMed]
4. Wacholder S, et al. Performance of common genetic variants in breast-cancer risk models. The New England Journal of Medicine. 2010;362(11):986–993. [PMC free article] [PubMed]
5. Conrad DF, Pinto D, Redon R, et al. Origins and functional impact of copy number variation in the human genome. Nature. 2010;464(7289):704–712. [PMC free article] [PubMed]
6. Nakamura Y. DNA variations in human and medical genetics: 25 years of my experience. Journal of Human Genetics. 2009;54(1):1–8. [PubMed]
7. Levy S, et al. The diploid genome sequence of an individual human. PLoS Biology. 2007;5(10, article e254) [PMC free article] [PubMed]
8. Seroussi E, Ron M, Kedra D. ShiftDetector: detection of shift mutations. Bioinformatics. 2002;18(8):1137–1138. [PubMed]
9. Dmitriev DA, Rakitov RA. Decoding of superimposed traces produced by direct sequencing of heterozygous indels. PLoS Computational Biology. 2008;4(7) Article ID e1000113. [PMC free article] [PubMed]
10. Dmitriev DA, Rakitov RA. Indelligent v.1.2, 2008,
11. Zhidkov I, Cohen R, Geifman N, Mishmar D, Rubin E. CHILD: a new tool for detecting low-abundance insertions and deletions in standard sequence traces. Nucleic Acids Research. 2011;39(7, article e47) [PMC free article] [PubMed]
12. Budowle B, Masibay A, Anderson SJ, et al. STR primer concordance study. Forensic Science International. 2001;124(1):47–54. [PubMed]
13. Cotton EA, Allsop RF, Guest JL, et al. Validation of the AMPFlSTR SGM Plus system for use in forensic casework. Forensic Science International. 2000;112(2-3):151–161. [PubMed]
14. Butler JM. Genetics and genomics of core short tandem repeat loci used in human identity testing. Journal of Forensic Sciences. 2006;51(2):253–265. [PubMed]
15. Butler JM. Short tandem repeat typing technologies used in human identity testing. BioTechniques. 2007;43(4):2–5. [PubMed]
16. The Federal Bureau of Investigation. Combined DNA Index System (CODIS)
17. Butler JM, Buel E, Crivellente F, McCord BR. Forensic DNA typing by capillary electrophoresis using the ABI Prism 310 and 3100 genetic analyzers for STR analysis. Electrophoresis. 2004;25(10-11):1397–1412. [PubMed]
18. Divne AM, Edlund H, Allen M. Forensic analysis of autosomal STR markers using Pyrosequencing. Forensic Science International. 2010;4(2):122–129. [PubMed]
19. Lin CY, Chao A, Yang YC, et al. Human papillomavirus typing with a polymerase chain reaction-based genotyping array compared with type-specific PCR. Journal of Clinical Virology. 2008;42(4):361–367. [PubMed]
20. Cook EH, Jr., Scherer SW. Copy-number variations associated with neuropsychiatric conditions. Nature. 2008;455(7215):919–923. [PubMed]
21. Gersemann M, Wehkamp J, Fellermann K, Stange EF. Crohn’s disease-defect in innate defence. World Journal of Gastroenterology. 2008;14(36):5499–5503. [PMC free article] [PubMed]
22. Groth M, Wiegand C, Szafranski K, et al. Both copy number and sequence variations affect expression of human DEFB4. Genes and Immunity. 2010;11(6):458–466. [PubMed]
23. Aldred PMR, Hollox EJ, Armour JAL. Copy number polymorphism and expression level variation of the human α-defensin genes DEFA1 and DEFA3. Human Molecular Genetics. 2005;14(14):2045–2052. [PubMed]
24. Hollox EJ, Huffmeier U, Zeeuwen PLJM, et al. Psoriasis is associated with increased β-defensin genomic copy number. Nature Genetics. 2008;40(1):23–25. [PMC free article] [PubMed]
25. Groth M, Szafranski K, Taudien S, et al. High-resolution mapping of the 8p23.1 beta-defensin cluster reveals strictly concordant copy number variation of all genes. Human Mutation. 2008;29(10):1247–1254. [PubMed]
26. Schouten JP, McElgunn CJ, Waaijer R, Zwijnenburg D, Diepvens F, Pals G. Relative quantification of 40 nucleic acid sequences by multiplex ligation-dependent probe amplification. Nucleic acids research. 2002;30(12):p. e57. [PMC free article] [PubMed]
27. Armour JAL, Palla R, Zeeuwen PLJM, Heijer MD, Schalkwijk J, Hollox EJ. Accurate, high-throughput typing of copy number variation using paralogue ratios from dispersed repeats. Nucleic Acids Research. 2007;35(3, article e19) [PMC free article] [PubMed]
28. Deutsch S, Choudhury U, Merla G, Howald C, Sylvan A, Antonarakis SE. Detection of aneuploidies by paralogous sequence quantification. Journal of Medical Genetics. 2004;41(12):908–915. [PMC free article] [PubMed]
29. Linzmeier RM, Ganz T. Copy number polymorphisms are not a common feature of innate immune genes. Genomics. 2006;88(1):122–126. [PubMed]
30. Townson JR, L.F. Barcellos, and RJNibbs. Gene copy number regulates the production of the human chemokine CCL3-L1. European Journal of Immunology. 2002;32(10):3016–3026. [PubMed]
31. Gonzalez E, Kulkarni H, Bolivar H, et al. The influence of CCL3L1 gene-containing segmental duplications on HIV-1/AIDS susceptibility. Science. 2005;307(5714):1434–1440. [PubMed]
32. Perry GH. The evolutionary significance of copy number variation in the human genome. Cytogenetic and Genome Research. 2009;123(1–4):283–287. [PMC free article] [PubMed]
33. Ngamphiw C, Kulawonganunchai S, Assawamakin A, Jenwitheesuk E, Tongsima S. VarDetect: a nucleotide sequence variation exploratory tool. BMC Bioinformatics. 2008;9(12, article S9) [PMC free article] [PubMed]
34. Lai CH, Chang CJ, Huang HJ, et al. Role of human papillomavirus genotype in prognosis of early-stage cervical cancer undergoing primary surgery. Journal of Clinical Oncology. 2007;25(24):3628–3634. [PubMed]
35. Yang SW, Lee YS, Chen TA, Wu CJ, Tsai CN. Human papillomavirus in oral leukoplakia is no prognostic indicator of malignant transformation. Cancer Epidemiology. 2009;33(2):118–122. [PubMed]
36. Duffy S, Shackelton LA, Holmes EC. Rates of evolutionary change in viruses: patterns and determinants. Nature Reviews Genetics. 2008;9(4):267–276. [PubMed]
37. Lee YS, Chao A, Chao AS, et al. CGcgh: a tool for molecular karyotyping using DNA microarray-based comparative genomic hybridization (array-CGH) Journal of Biomedical Science. 2008;15(6):687–696. [PubMed]
38. Chen K, McLellan MD, Ding L, et al. PolyScan: an automatic indel and SNP detection approach to the analysis of human resequencing data. Genome Research. 2007;17(5):659–666. [PubMed]

Articles from The Scientific World Journal are provided here courtesy of Hindawi Limited