DNA Commission of the International Society for Forensic Genetics: Recommendations on the validation of software programs performing biostatistical calculations for forensic genetics applications.

The use of biostatistical software programs to assist in data interpretation and calculate likelihood ratios is essential to forensic geneticists and part of the daily case work flow for both kinship and DNA identification laboratories. Previous recommendations issued by the DNA Commission of the International Society for Forensic Genetics (ISFG) covered the application of bio-statistical evaluations for STR typing results in identification and kinship cases, and this is now being expanded to provide best practices regarding validation and verification of the software required for these calculations. With larger multiplexes, more complex mixtures, and increasing requests for extended family testing, laboratories are relying more than ever on specific software solutions and sufficient validation, training and extensive documentation are of upmost importance. Here, we present recommendations for the minimum requirements to validate bio-statistical software to be used in forensic genetics. We distinguish between developmental validation and the responsibilities of the software developer or provider, and the internal validation studies to be performed by the end user. Recommendations for the software provider address, for example, the documentation of the underlying models used by the software, validation data expectations, version control, implementation and training support, as well as continuity and user notifications. For the internal validations the recommendations include: creating a validation plan, requirements for the range of samples to be tested, Standard Operating Procedure development, and internal laboratory training and education. To ensure that all laboratories have access to a wide range of samples for validation and training purposes the ISFG DNA commission encourages collaborative studies and public repositories of STR typing results.

Alu insertion polymorphisms and the genetic structure of human populations from the Caucasus.

An analysis of 8 Alu insertion loci (ACE, TPA25, PV92, APO, FXIIIB, D1, A25, B65) has been carried out in six populations from the Caucasus, including Indo-European-speaking Armenians; Altaic-speaking Azerbaijanians; North Caucasian-speaking Cherkessians, Darginians, and Ingushians; and South Caucasian (Kartvelian)-speaking Georgians. The Caucasus populations exhibit low levels of within-population variation and high levels of between-population differentiation, with the average Fst value for the Caucasus of 0.113, which is almost as large as the Fst value of 0.157 for worldwide populations. Maximum likelihood tree and principal coordinate analyses both group the Caucasus populations with European populations. Neither geographic nor linguistic relationships appear to explain the genetic relationships of Caucasus populations. Instead, it appears as if they have been small and relatively isolated, and hence genetic drift has been the dominant influence on the genetic structure of Caucasus populations.

dbSNP: the NCBI database of genetic variation.

In response to a need for a general catalog of genome variation to address the large-scale sampling designs required by association studies, gene mapping and evolutionary biology, the National Center for Biotechnology Information (NCBI) has established the dbSNP database [S.T.Sherry, M.Ward and K. Sirotkin (1999) Genome Res., 9, 677-679]. Submissions to dbSNP will be integrated with other sources of information at NCBI such as GenBank, PubMed, LocusLink and the Human Genome Project data. The complete contents of dbSNP are available to the public at website: http://www.ncbi.nlm.nih.gov/SNP. The complete contents of dbSNP can also be downloaded in multiple formats via anonymous FTP at ftp://ncbi.nlm.nih.gov/snp/.

Reading between the LINEs: human genomic variation induced by LINE-1 retrotransposition.

The insertion of mobile elements into the genome represents a new class of genetic markers for the study of human evolution. Long interspersed elements (LINEs) have amplified to a copy number of about 100,000 over the last 100 million years of mammalian evolution and comprise approximately 15% of the human genome. The majority of LINE-1 (L1) elements within the human genome are 5' truncated copies of a few active L1 elements that are capable of retrotransposition. Some of the young L1 elements have inserted into the human genome so recently that populations are polymorphic for the presence of an L1 element at a particular chromosomal location. L1 insertion polymorphisms offer several advantages over other types of polymorphisms for human evolution studies. First, they are typed by rapid, simple, polymerase chain reaction (PCR)-based assays. Second, they are stable polymorphisms that rarely undergo deletion. Third, the presence of an L1 element represents identity by descent, because the probability is negligible that two different young L1 repeats would integrate independently between the exact same two nucleotides. Fourth, the ancestral state of L1 insertion polymorphisms is known to be the absence of the L1 element, which can be used to root plots/trees of population relationships. Here we report the development of a PCR-based display for the direct identification of dimorphic L1 elements from the human genome. We have also developed PCR-based assays for the characterization of six polymorphic L1 elements within the human genome. PCR analysis of human/rodent hybrid cell line DNA samples showed that the polymorphic L1 elements were located on several different chromosomes. Phylogenetic analysis of nonhuman primate DNA samples showed that all of the recently integrated "young" L1 elements were restricted to the human genome and absent from the genomes of nonhuman primates. Analysis of a diverse array of human populations showed that the allele frequencies and level of heterozygosity for each of the L1 elements was variable. Polymorphic L1 elements represent a new source of identical-by-descent variation for the study of human evolution. [The sequence data described in this paper have been submitted to the GenBank data library under accession nos. AF242435-AF242451.]

Long DOP-PCR of rare archival anthropological samples.

The application of molecular DNA technologies to anthropological questions has meant that rare or archival samples of human remains, including blood, hair, and bone, can now be used as a source of material for genetic analysis. Often, these samples are irreplaceable, and/or yield very small quantities of DNA, so methods for preamplifying as much of the whole genome as possible would greatly enhance their usefulness. DOP-PCR (degenerate oligonucleotide-primed polymerase chain reaction) is an amplification method that uses a degenerate primer and very low initial annealing temperatures to amplify the whole genome. We adapted a published DOP-PCR protocol to long PCR enzyme and amplification conditions. The effectiveness of these modifications was tested by PCR amplification of DOP-PCR products at a mixture of genomic targets including 66 different microsatellites, 11 Alu insertion polymorphisms, and variable-length segments of the human lipoprotein lipase gene (LPL). The selected microsatellite markers were chosen to represent every chromosome, with expected product sizes ranging from 150 base pairs to 8,000 base pairs in length, while the 22 Alu insertion polymorphisms were selected to reveal biases in the recovery of alleles of different sizes. To determine nucleotide sequence variation, 2 kilobases (kb) of the LPL gene in 30 Mongolian individuals were sequenced. All gene-specific targets from DOP-PCR product template were amplified. No unexpected polymorphisms in the sequence results attributable to the DOP-PCR step were found, and 93% to 95% of Alu genotypes that have been amplified from total genomic DNA were replicated. The incorrect typings were all due to the preferential amplification of the shorter of two possible alleles in individuals heterozygous for an Alu insertion and were all correctly typed on subsequent reamplification of the gene-specific PCR products. This method of whole-genome amplification promises to be an efficient way to maximize the genetic use of rare anthropological samples.

dbSNP: a database of single nucleotide polymorphisms.

In response to a need for a general catalog of genome variation to address the large-scale sampling designs required by association studies, gene mapping and evolutionary biology, the National Cancer for Biotechnology Information (NCBI) has established the dbSNP database. Submissions to dbSNP will be integrated with other sources of information at NCBI such as GenBank, PubMed, LocusLink and the Human Genome Project data. The complete contents of dbSNP are available to the public at website: http://www.ncbi.nlm.nih.gov/SNP. Submitted SNPs can also be downloaded via anonymous FTP at ftp://ncbi.nlm.nih.gov/snp/

Use of molecular variation in the NCBI dbSNP database.

While high quality information regarding variation in genes is currently available in locus-specific or specialized mutation databases, the need remains for a general catalog of genome variation to address the large-scale sampling designs required by association studies, gene mapping, and evolutionary biology. In response to this need, the National Center for Biotechnology Information (NCBI) has established the dbSNP database http://ncbi. nlm.nih.gov/SNP/ to serve as a generalized, central variation database. Submissions to dbSNP will be integrated with other sources of information at NCBI such as GenBank, PubMed, LocusLink, and the Human Genome Project data, and the complete contents of dbSNP are available to the public via anonymous FTP. Hum Mutat 15:68-75, 2000. Published 2000 Wiley-Liss, Inc.

Assembly of a high-resolution map of the Acadian Usher syndrome region and localization of the nuclear EF-hand acidic gene.

Usher syndrome type 1C (USH1C) occurs in a small population of Acadian descendants from southwestern Louisiana. Linkage and linkage disequilibrium analyses localize USH1C to chromosome 11p between markers D11S1397 and D11S1888, an interval of less than 680 kb. Here, we refine the USH1C linkage to a region less than 400 kb, between genetic markers D11S1397 and D11S1890. Using 17 genetic markers from this interval, we have isolated a contiguous set of 60 bacterial artificial chromosomes (BACs) that span the USH1C critical region. Exon trapping of BAC clones from this region resulted in the recovery of an exon of the nuclear EF-hand acidic (NEFA) gene. However, DNA sequence analysis of the NEFA cDNA from lymphocytes of affected individuals provided no evidence of mutation, making structural mutations in the NEFA protein unlikely as the cellular cause of Acadian Usher syndrome.

Genetic traces of ancient demography.

Patterns of gene differences among humans contain information about the demographic history of our species. Haploid loci like mitochondrial DNA and the nonrecombining part of the Y chromosome show a pattern indicating expansion from a population of only several thousand during the late middle or early upper Pleistocene. Nuclear short tandem repeat loci also show evidence of this expansion. Both mitochondrial DNA and the Y chromosome coalesce within the last several hundred thousand years, and they cannot provide information about the population before their coalescence. Several nuclear loci are informative about our ancestral population size during nearly the whole Pleistocene. They indicate a small effective size, on the order of 10,000 breeding individuals, throughout this time period. This genetic evidence denies any version of the multiregional model of modern human origins. It implies instead that our ancestors were effectively a separate species for most of the Pleistocene.

Alu evolution in human populations: using the coalescent to estimate effective population size.

There are estimated to be approximately 1000 members of the Ya5 Alu subfamily of retroposons in humans. This subfamily has a distribution restricted to humans, with a few copies in gorillas and chimpanzees. Fifty-seven Ya5 elements were previously cloned from a HeLa-derived randomly sheared total genomic library, sequenced, and screened for polymorphism in a panel of 120 unrelated humans. Forty-four of the 57 cloned Alu repeats were monomorphic in the sample and 13 Alu repeats were dimorphic for insertion presence/absence. The observed distribution of sample frequencies of the 13 dimorphic elements is consistent with the theoretical expectation for elements ascertained in a single diploid cell line. Coalescence theory is used to compute expected total pedigree branch lengths for monomorphic and dimorphic elements, leading to an estimate of human effective population size of approximately 18,000 during the last one to two million years.

Identification and characterization of two polymorphic Ya5 Alu repeats.

Two new polymorphic Alu elements (HS2.25 and HS4.14) belonging to the young (Ya5/8) subfamily of human-specific Alu repeats have been identified. DNA sequence analysis of both Alu repeats revealed that each Alu repeat had a long 3'-oligo-dA-rich tail (41 and 52 nucleotides in length) and a low level of random mutations. HS2.25 and HS4.14 were flanked by short precise direct repeats of 8 and 14 nucleotides in length, respectively. HS2.25 was located on human chromosome 13, and HS4.14 on chromosome 1. Both Alu elements were absent from the orthologous positions within the genomes of non-human primates, and were highly polymorphic in a survey of twelve geographically diverse human groups.

Evolutionary history of the COII/tRNALys intergenic 9 base pair deletion in human mitochondrial DNAs from the Pacific.

Length changes in human mitochondrial DNA (mtDNA) are potentially useful markers for inferring the evolutionary history of populations. One such length change is a nine base pair (9-bp) deletion that is located in the intergenic region between the COII gene and the Lysine tRNA gene (COII/tRNALys intergenic region). This deletion has been used as a genetic marker to trace descent from peoples of East Asian origin. A geographic cline of the deletion frequency across modern Pacific Islander populations suggests that the deletion may be useful for tracing prehistoric Polynesian origins and affinities. Mitochondrial DNA sequence variation within two variable segments of the control region (CR) permits a number of inferences regarding the evolutionary history of the 9-bp deletion that cannot be determined from frequency data alone. We obtained CR sequences from 74 mtDNAs with the 9-bp deletion from Indonesia, coastal Papua New Guinea (PNG), and American Samoa. Phylogenetic and pairwise distribution analysis of these CR sequences pooled with previously published CR sequences reveals that the deletion arose independently in Africa and Asia and suggests possible multiple origins of the deletion in Asia. A clinal increase of the frequency of the 9-bp deletion across the three Pacific populations is associated with a decrease in CR sequence diversity, consistent with founder events. Furthermore, analysis of pairwise difference distributions indicates an expansion time of proto-Polynesians that began 5,500 yr ago from Southeast Asia. These results are consistent with the express train model of Polynesian origins.

Mismatch distributions of mtDNA reveal recent human population expansions.

Although many genetic studies of human evolution have tried to make distinctions between the replacement and the multiregional evolution hypotheses, current methods and data have not resolved the issue. However, new advances in nucleotide divergence theory can complement these investigations with a description of human demographic behavior during the late Middle and Upper Paleolithic (approximately the last 250,000 years). Restriction fragment length polymorphism (RFLP) and DNA sequence analyses of human mitochondrial DNAs (mtDNAs) from 25 ethnic and racial groups indicate that significant expansions occurred during the late Middle and Upper Paleolithic in 23 of the 25 populations examined. Estimates for the individual group expansion times are consistently less than 100,000 years ago with a mean expansion time of approximately 40,000 years ago. The dramatic expansions suggested by these data occurred well after modern human anatomy appeared, approximately 100,000 years ago, but are concordant with archeological evidence for the expansion of modern human technology, approximately 50,000 years ago.

New approaches to dating suggest a recent age for the human mtDNA ancestor.

The most critical and controversial feature of the African origin hypothesis of human mitochondrial DNA (mtDNA) evolution is the relatively recent age of about 200 ka inferred for the human mtDNA ancestor. If this age is wrong, and the actual age instead approaches 1 million years ago, then the controversy abates. Reliable estimates of the age of the human mtDNA ancestor and the associated standard error are therefore crucial. However, more recent estimates of the age of the human ancestor rely on comparisons between human and chimpanzee mtDNAs that may not be reliable and for which standard errors are difficult to calculate. We present here two approaches for deriving an intraspecific calibration of the rate of human mtDNA sequence evolution that allow standard errors to be readily calculated. The estimates resulting from these two approaches for the age of the human mtDNA ancestor (and approximate 95% confidence intervals) are 133 (63-356) and 137 (63-416) ka ago. These results provide the strongest evidence yet for a relatively recent origin of the human mtDNA ancestor.

The structure of human mitochondrial DNA variation.

Restriction analysis of mitochondrial DNA (mtDNA) of 3065 humans from 62 geographic samples identified 149 haplotypes and 81 polymorphic sites. These data were used to test several aspects of the evolutionary past of the human species. A dendrogram depicting the genetic relatedness of all haplotypes shows that the native African populations have the greatest diversity and, consistent with evidence from a variety of sources, suggests an African origin for our species. The data also indicate that two individuals drawn at random from the entire sample will differ at approximately 0.4% of their mtDNA nucleotide sites, which is somewhat higher than previous estimates. Human mtDNA also exhibits more interpopulation heterogeneity (GST = 0.351 +/- 0.025) than does nuclear DNA (GST = 0.12). Moreover, the virtual absence of intermediate levels of linkage disequilibrium between pairs of sites is consistent with the absence of genetic recombination and places constraints on the rate of mutation. Tests of the selective neutrality of mtDNA variation, including the Ewens-Watterson and Tajima tests, indicate a departure in the direction consistent with purifying selection, but this departure is more likely due to the rapid growth of the human population and the geographic heterogeneity of the variation. The lack of a good fit to neutrality poses problems for the estimation of times of coalescence from human mtDNA data.