Janet S. Ziegle

Application of automated DNA sizing technology for genotyping microsatellite loci

Highly polymorphic microsatellite loci offer great promise for gene mapping studies, but fulfillment of this potential will require substantial improvements in methods for accurate and efficient genotyping. Here, we report a genotyping method based on fluorescently labeled PCR primers and size characterization of PCR products using an automated DNA fragment analyzer. We capitalize on the availability of three distinct fluorescent dyes to label uniquely loci that overlap in size, and this innovation increases by threefold the number of loci that can be analyzed simultaneously. We label size standards with a fourth dye and combine these with the microsatellite PCR products in each gel lane. Computer programs provide very rapid and accurate sizing of microsatellite alleles and efficient data management. In addition, fluorescence signals are linear over a much greater range of intensity than conventional autoradiography. This facilitates multiplexing of loci (since signal intensities often vary greatly) and helps distinguish major peaks from artifacts, thereby improving genotyping accuracy.

A comprehensive linkage map of the dog genome.

We have leveraged the reference sequence of a boxer to construct the first complete linkage map for the domestic dog. The new map improves access to the dogs unique biology, from human disease counterparts to fascinating evolutionary adaptations. The map was constructed with ∼3000 microsatellite markers developed from the reference sequence. Familial resources afforded 450 mostly phase-known meioses for map assembly. The genotype data supported a framework map with ∼1500 loci. An additional ∼1500 markers served as map validators, contributing modestly to estimates of recombination rate but supporting the framework content. Data from ∼22,000 SNPs informing on a subset of meioses supported map integrity. The sex-averaged map extended 21 M and revealed marked region- and sex-specific differences in recombination rate. The map will enable empiric coverage estimates and multipoint linkage analysis. Knowledge of the variation in recombination rate will also inform on genomewide patterns of linkage disequilibrium (LD), and thus benefit association, selective sweep, and phylogenetic mapping approaches. The computational and wet-bench strategies can be applied to the reference genome of any nonmodel organism to assemble a de novo linkage map.

The prevalence of folate-remedial MTHFR enzyme variants in humans.

Studies of rare, inborn metabolic diseases establish that the phenotypes of some mutations in vitamin-dependent enzymes can be suppressed by supplementation of the cognate vitamin, which restores function of the defective enzyme. To determine whether polymorphisms exist that more subtly affect enzymes yet are augmentable in the same way, we sequenced the coding region of a prototypical vitamin-dependent enzyme, methylenetetrahydrofolate reductase (MTHFR), from 564 individuals of diverse ethnicities. All nonsynonymous changes were evaluated in functional in vivo assays in Saccharomyces cerevisiae to identify enzymatic defects and folate remediability of impaired alleles. We identified 14 nonsynonymous changes: 11 alleles with minor allele frequencies <1% and 3 common alleles (A222V, E429A, and R594Q). Four of 11 low-frequency alleles affected enzyme function, as did A222V. Of the five impaired alleles, four could be restored to normal functionality by elevating intracellular folate levels. All five impaired alleles mapped to the N-terminal catalytic domain of the enzyme, whereas changes in the C-terminal regulatory domain had little effect on activity. Impaired activity correlated with the phosphorylation state of MTHFR, with more severe mutations resulting in lower abundance of the phosphorylated protein. Significantly, diploid yeast heterozygous for mutant alleles were impaired for growth, particularly with lower folate supplementation. These results suggested that multiple less-frequent alleles, in aggregate, might significantly contribute to metabolic dysfunction. Furthermore, vitamin remediation of mutant enzymes may be a common phenomenon in certain domains of proteins.

Y-chromosome analysis reveals genetic divergence and new founding native lineages in Athapaskan- and Eskimoan-speaking populations

For decades, the peopling of the Americas has been explored through the analysis of uniparentally inherited genetic systems in Native American populations and the comparison of these genetic data with current linguistic groupings. In northern North America, two language families predominate: Eskimo-Aleut and Na-Dene. Although the genetic evidence from nuclear and mtDNA loci suggest that speakers of these language families share a distinct biological origin, this model has not been examined using data from paternally inherited Y chromosomes. To test this hypothesis and elucidate the migration histories of Eskimoan- and Athapaskan-speaking populations, we analyzed Y-chromosomal data from Inuvialuit, Gwich’in, and Tłįchǫ populations living in the Northwest Territories of Canada. Over 100 biallelic markers and 19 chromosome short tandem repeats (STRs) were genotyped to produce a high-resolution dataset of Y chromosomes from these groups. Among these markers is an SNP discovered in the Inuvialuit that differentiates them from other Aboriginal and Native American populations. The data suggest that Canadian Eskimoan- and Athapaskan-speaking populations are genetically distinct from one another and that the formation of these groups was the result of two population expansions that occurred after the initial movement of people into the Americas. In addition, the population history of Athapaskan speakers is complex, with the Tłįchǫ being distinct from other Athapaskan groups. The high-resolution biallelic data also make clear that Y-chromosomal diversity among the first Native Americans was greater than previously recognized.

Population Differentiation of Southern Indian Male Lineages Correlates with Agricultural Expansions Predating the Caste System

Previous studies that pooled Indian populations from a wide variety of geographical locations, have obtained contradictory conclusions about the processes of the establishment of the Varna caste system and its genetic impact on the origins and demographic histories of Indian populations. To further investigate these questions we took advantage that both Y chromosome and caste designation are paternally inherited, and genotyped 1,680 Y chromosomes representing 12 tribal and 19 non-tribal (caste) endogamous populations from the predominantly Dravidian-speaking Tamil Nadu state in the southernmost part of India. Tribes and castes were both characterized by an overwhelming proportion of putatively Indian autochthonous Y-chromosomal haplogroups (H-M69, F-M89, R1a1-M17, L1-M27, R2-M124, and C5-M356; 81% combined) with a shared genetic heritage dating back to the late Pleistocene (10–30 Kya), suggesting that more recent Holocene migrations from western Eurasia contributed <20% of the male lineages. We found strong evidence for genetic structure, associated primarily with the current mode of subsistence. Coalescence analysis suggested that the social stratification was established 4–6 Kya and there was little admixture during the last 3 Kya, implying a minimal genetic impact of the Varna (caste) system from the historically-documented Brahmin migrations into the area. In contrast, the overall Y-chromosomal patterns, the time depth of population diversifications and the period of differentiation were best explained by the emergence of agricultural technology in South Asia. These results highlight the utility of detailed local genetic studies within India, without prior assumptions about the importance of Varna rank status for population grouping, to obtain new insights into the relative influences of past demographic events for the population structure of the whole of modern India.

Multiplex single-nucleotide polymorphism typing of the human Y chromosome using TaqMan probes

BackgroundThe analysis of human Y-chromosome variation in the context of population genetics and forensics requires the genotyping of dozens to hundreds of selected single-nucleotide polymorphisms (SNPs). In the present study, we developed a 121-plex (121 SNPs in a single array) TaqMan array capable of distinguishing most haplogroups and subhaplogroups on the Y-chromosome human phylogeny in Europe.ResultsWe present data from 264 samples from several European areas and ethnic groups. The array developed in this study shows >99% accuracy of assignation to the Y human phylogeny (with an average call rate of genotypes >96%).ConclusionsWe have created and evaluated a robust and accurate Y-chromosome multiplex which minimises the possible errors due to mixup when typing the same sample in several independent reactions.

Genome-wide signatures of male-mediated migration shaping the Indian gene pool

Multiple questions relating to contributions of cultural and demographical factors in the process of human geographical dispersal remain largely unanswered. India, a land of early human settlement and the resulting diversity is a good place to look for some of the answers. In this study, we explored the genetic structure of India using a diverse panel of 78 males genotyped using the GenoChip. Their genome-wide single-nucleotide polymorphism (SNP) diversity was examined in the context of various covariates that influence Indian gene pool. Admixture analysis of genome-wide SNP data showed high proportion of the Southwest Asian component in all of the Indian samples. Hierarchical clustering based on admixture proportions revealed seven distinct clusters correlating to geographical and linguistic affiliations. Convex hull overlay of Y-chromosomal haplogroups on the genome-wide SNP principal component analysis brought out distinct non-overlapping polygons of F*-M89, H*-M69, L1-M27, O2a-M95 and O3a3c1-M117, suggesting a male-mediated migration and expansion of the Indian gene pool. Lack of similar correlation with mitochondrial DNA clades indicated a shared genetic ancestry of females. We suggest that ancient male-mediated migratory events and settlement in various regional niches led to the present day scenario and peopling of India.

Aboriginal Australian mitochondrial genome variation – an increased understanding of population antiquity and diversity

Aboriginal Australians represent one of the oldest continuous cultures outside Africa, with evidence indicating that their ancestors arrived in the ancient landmass of Sahul (present-day New Guinea and Australia) ~55 thousand years ago. Genetic studies, though limited, have demonstrated both the uniqueness and antiquity of Aboriginal Australian genomes. We have further resolved known Aboriginal Australian mitochondrial haplogroups and discovered novel indigenous lineages by sequencing the mitogenomes of 127 contemporary Aboriginal Australians. In particular, the more common haplogroups observed in our dataset included M42a, M42c, S, P5 and P12, followed by rarer haplogroups M15, M16, N13, O, P3, P6 and P8. We propose some major phylogenetic rearrangements, such as in haplogroup P where we delinked P4a and P4b and redefined them as P4 (New Guinean) and P11 (Australian), respectively. Haplogroup P2b was identified as a novel clade potentially restricted to Torres Strait Islanders. Nearly all Aboriginal Australian mitochondrial haplogroups detected appear to be ancient, with no evidence of later introgression during the Holocene. Our findings greatly increase knowledge about the geographic distribution and phylogenetic structure of mitochondrial lineages that have survived in contemporary descendants of Australia’s first settlers.

Mitochondrial DNA diversity of present-day Aboriginal Australians and implications for human evolution in Oceania

Aboriginal Australians are one of the more poorly studied populations from the standpoint of human evolution and genetic diversity. Thus, to investigate their genetic diversity, the possible date of their ancestors’ arrival and their relationships with neighboring populations, we analyzed mitochondrial DNA (mtDNA) diversity in a large sample of Aboriginal Australians. Selected mtDNA single-nucleotide polymorphisms and the hypervariable segment haplotypes were analyzed in 594 Aboriginal Australians drawn from locations across the continent, chiefly from regions not previously sampled. Most (~78%) samples could be assigned to mtDNA haplogroups indigenous to Australia. The indigenous haplogroups were all ancient (with estimated ages >40 000 years) and geographically widespread across the continent. The most common haplogroup was P (44%) followed by S (23%) and M42a (9%). There was some geographic structure at the haplotype level. The estimated ages of the indigenous haplogroups range from 39 000 to 55 000 years, dates that fit well with the estimated date of colonization of Australia based on archeological evidence (~47 000 years ago). The distribution of mtDNA haplogroups in Australia and New Guinea supports the hypothesis that the ancestors of Aboriginal Australians entered Sahul through at least two entry points. The mtDNA data give no support to the hypothesis of secondary gene flow into Australia during the Holocene, but instead suggest long-term isolation of the continent.

A robust targeted sequencing approach for low input and variable quality DNA from clinical samples

Next-generation deep sequencing of gene panels is being adopted as a diagnostic test to identify actionable mutations in cancer patient samples. However, clinical samples, such as formalin-fixed, paraffin-embedded specimens, frequently provide low quantities of degraded, poor quality DNA. To overcome these issues, many sequencing assays rely on extensive PCR amplification leading to an accumulation of bias and artifacts. Thus, there is a need for a targeted sequencing assay that performs well with DNA of low quality and quantity without relying on extensive PCR amplification. We evaluate the performance of a targeted sequencing assay based on Oligonucleotide Selective Sequencing, which permits the enrichment of genes and regions of interest and the identification of sequence variants from low amounts of damaged DNA. This assay utilizes a repair process adapted to clinical FFPE samples, followed by adaptor ligation to single stranded DNA and a primer-based capture technique. Our approach generates sequence libraries of high fidelity with reduced reliance on extensive PCR amplification—this facilitates the accurate assessment of copy number alterations in addition to delivering accurate single nucleotide variant and insertion/deletion detection. We apply this method to capture and sequence the exons of a panel of 130 cancer-related genes, from which we obtain high read coverage uniformity across the targeted regions at starting input DNA amounts as low as 10 ng per sample. We demonstrate the performance using a series of reference DNA samples, and by identifying sequence variants in DNA from matched clinical samples originating from different tissue types.Cancer diagnostics: Targeted DNA sequencing for low-quality tumor samplesA new DNA sequencing technology enables comprehensive genetic analyses of poor-quality tumor samples. Hanlee Ji from Stanford University in California, USA, together with colleagues from a company he cofounded called TOMA Biosciences, tested the performance of a targeted sequencing assay known as oligonucleotide-selective sequencing (OS-Seq). They used the “in-solution” version of OS-Seq, which involves a pre-processing step to remove any damaged DNA and then sequences target regions of the genome to look for duplications, insertions or deletions of DNA segments. Using archival specimens (which often contain low quantities of degraded DNA) from patients with lung and colorectal cancer, the researchers showed they could detect sequence variants in a panel of 130 cancer-related genes. The findings suggest the OS-Seq assay could help inform treatment decisions for cancer patients, even with clinical specimens of low quality.