Whole-genome resequencing reveals loci under selection during chicken domestication.

Domestic animals are excellent models for genetic studies of phenotypic evolution. They have evolved genetic adaptations to a new environment, the farm, and have been subjected to strong human-driven selection leading to remarkable phenotypic changes in morphology, physiology and behaviour. Identifying the genetic changes underlying these developments provides new insight into general mechanisms by which genetic variation shapes phenotypic diversity. Here we describe the use of massively parallel sequencing to identify selective sweeps of favourable alleles and candidate mutations that have had a prominent role in the domestication of chickens (Gallus gallus domesticus) and their subsequent specialization into broiler (meat-producing) and layer (egg-producing) chickens. We have generated 44.5-fold coverage of the chicken genome using pools of genomic DNA representing eight different populations of domestic chickens as well as red jungle fowl (Gallus gallus), the major wild ancestor. We report more than 7,000,000 single nucleotide polymorphisms, almost 1,300 deletions and a number of putative selective sweeps. One of the most striking selective sweeps found in all domestic chickens occurred at the locus for thyroid stimulating hormone receptor (TSHR), which has a pivotal role in metabolic regulation and photoperiod control of reproduction in vertebrates. Several of the selective sweeps detected in broilers overlapped genes associated with growth, appetite and metabolic regulation. We found little evidence that selection for loss-of-function mutations had a prominent role in chicken domestication, but we detected two deletions in coding sequences that we suggest are functionally important. This study has direct application to animal breeding and enhances the importance of the domestic chicken as a model organism for biomedical research.

Ultra-deep sequencing of mouse mitochondrial DNA: mutational patterns and their origins.

Somatic mutations of mtDNA are implicated in the aging process, but there is no universally accepted method for their accurate quantification. We have used ultra-deep sequencing to study genome-wide mtDNA mutation load in the liver of normally- and prematurely-aging mice. Mice that are homozygous for an allele expressing a proof-reading-deficient mtDNA polymerase (mtDNA mutator mice) have 10-times-higher point mutation loads than their wildtype siblings. In addition, the mtDNA mutator mice have increased levels of a truncated linear mtDNA molecule, resulting in decreased sequence coverage in the deleted region. In contrast, circular mtDNA molecules with large deletions occur at extremely low frequencies in mtDNA mutator mice and can therefore not drive the premature aging phenotype. Sequence analysis shows that the main proportion of the mutation load in heterozygous mtDNA mutator mice and their wildtype siblings is inherited from their heterozygous mothers consistent with germline transmission. We found no increase in levels of point mutations or deletions in wildtype C57Bl/6N mice with increasing age, thus questioning the causative role of these changes in aging. In addition, there was no increased frequency of transversion mutations with time in any of the studied genotypes, arguing against oxidative damage as a major cause of mtDNA mutations. Our results from studies of mice thus indicate that most somatic mtDNA mutations occur as replication errors during development and do not result from damage accumulation in adult life.

Genetic origin of the Swedish Sami inferred from HLA class I and class II allele frequencies.

Sami of northern Scandinavia are genetic outliers among European populations and their origin has been difficult to determine. In order to study the genetic origin of the Swedish Sami, we have performed high-resolution typing of the class I HLA-A and -B loci and the class II DRB1, DQB1 and DQA1 loci in the northern and southern Swedish Sami. Several of the common class I alleles in Sami (B*0702, B*1501, B*4002 and A*0301) are found at high frequency in other European populations. However, a number of class I and class II alleles (B*4001, A*2402, DRB1*0901 and DRB1*1101) in the Swedish Sami are characteristic of Asian populations. Admixture analyses indicate that 87% of the Sami gene pool is of European origin and that the Asian contribution is 13%. Our HLA analyses indicate a higher proportion of Asian ancestry in the Sami than shown by previous genetic studies.

Polymorphisms in 9q32 and TSCOT are linked to cervical cancer in affected sib-pairs with high mean age at diagnosis.

Cervical cancer is a multifactorial disease influenced by both environmental and genetic factors. We have previously found linkage to 9q32 in a genomewide scan of affected sib-pairs (ASPs) with cervical cancer and to the thymic stromal co-transporter (TSCOT), a candidate gene in this region. Here we examined the contribution of 9q32 and TSCOT to cervical cancer susceptibility using at larger material of 641 ASPs, 278 of which were included in the earlier genome-scan. Since heritable forms of cancer frequently show stronger genetic effects in families with early onset of disease, we stratified the ASPs into two groups based on mean age at diagnosis (MAAD) within sib-pairs. Surprisingly, ASPs with high MAAD (30.5-47.5 years) showed increased sharing at all microsatellite markers at 9q31.1-33.1 and linkage signals of up to MLS = 2.74 for TSCOT SNPs, while ASPs with low MAAD (19-30 years) showed no deviation from random genetic sharing (MLS = 0.00). The difference in allelic sharing between the two MAAD strata was significant (P < 0.005) and is not likely to be explained by the HLA haplotype, a previously known genetic susceptibility factor for cervical cancer. Our results indicate locus heterogeneity in the susceptibility to cervical cancer between the two strata, with polymorphisms in the 9q32 region mainly showing an effect in women with high MAAD.

mtDB: Human Mitochondrial Genome Database, a resource for population genetics and medical sciences.

The mitochondrial genome, contained in the subcellular mitochondrial network, encodes a small number of peptides pivotal for cellular energy production. Mitochondrial genes are highly polymorphic and cataloguing existing variation is of interest for medical scientists involved in the identification of mutations causing mitochondrial dysfunction, as well as for population genetics studies. Human Mitochondrial Genome Database (mtDB) (http://www.genpat.uu.se/mtDB) has provided a comprehensive database of complete human mitochondrial genomes since early 2000. At this time, owing to an increase in the number of published complete human mitochondrial genome sequences, it became necessary to provide a web-based database of human whole genome and complete coding region sequences. As of August 2005 this database contains 2104 sequences (1544 complete genome and 560 coding region) available to download or search for specific polymorphisms. Of special interest to medical researchers and population geneticists evaluating specific positions is a complete list of (currently 3311) mitochondrial polymorphisms among these sequences. Recent expansions in the capabilities of mtDB include a haplotype search function and the ability to identify and download sequences carrying particular variants.

A recent genetic link between Sami and the Volga-Ural region of Russia.

The genetic origin of the Sami is enigmatic and contributions from Continental Europe, Eastern Europe and Asia have been proposed. To address the evolutionary history of northern and southern Swedish Sami, we have studied their mtDNA haplogroup frequencies and complete mtDNA genome sequences. While the majority of mtDNA diversity in the northern Swedish, Norwegian and Finnish Sami is accounted for by haplogroups V and U5b1b1, the southern Swedish Sami have other haplogroups and a frequency distribution similar to that of the Continental European population. Stratification of the southern Sami on the basis of occupation indicates that this is the result of recent admixture with the Swedish population. The divergence time for the Sami haplogroup V sequences is 7600 YBP (years before present), and for U5b1b1, 5500 YBP amongst Sami and 6600 YBP amongst Sami and Finns. This suggests an arrival in the region soon after the retreat of the glacial ice, either by way of Continental Europe and/or the Volga-Ural region. Haplogroup Z is found at low frequency in the Sami and Northern Asian populations but is virtually absent in Europe. Several conserved substitutions group the Sami Z lineages strongly with those from Finland and the Volga-Ural region of Russia, but distinguish them from Northeast Asian representatives. This suggests that some Sami lineages shared a common ancestor with lineages from the Volga-Ural region as recently as 2700 years ago, indicative of a more recent contribution of people from the Volga-Ural region to the Sami population.

Lifestyle, genetics, and disease in Sami.

AIM: To present a summary of the lifestyle, genetic origin, diet, and disease in the population of Sami, indigenous people of northern Fennoscandia. METHOD: A survey of the available scientific literature and preliminary results from our own study of the Swedish Sami population. RESULTS: The Sami probably have a heterogeneous genetic origin, with a major contribution of continental or Eastern European tribes and a smaller contribution from Asia. The traditional Sami diet, high in animal products, persists in Sami groups still involved with reindeer herding, but others have adopted a diet typical of Western cultures. Early reports indicated a lower prevalence of heart disease and most cancers, except stomach cancer. Recent studies have not found a lower risk of heart disease, but have consistently shown an overall reduced cancer risk. Sami have been reported to share some specific health-related genetic polymorphisms with other European populations, but none that would explain the observed differences in disease risk. CONCLUSION: The genetic structure of the Sami population makes it suitable for studies of the genetic and environmental factors influencing the development of common diseases. The difference in incidence of heart disease between studies may reflect the ongoing transition from a traditional to a more Westernized lifestyle. The ability to compare population segments with different lifestyles, combined with the genetic structure of the population, creates unusual possibilities for studies of the genetic and environmental factors involved in the development of common disease.

SNP frequency estimation using massively parallel sequencing of pooled DNA.

Resequencing of genomic regions that have been implicated by linkage and/or association studies to harbor genetic susceptibility loci represents a necessary step to identify causal variants. Massively parallel sequencing (MPS) offers the possibility of SNP discovery and frequency determination among pooled DNA samples. The strategies of pooling DNA samples and pooling PCR amplicons generated from individual DNA samples were evaluated, and both were found to return accurate estimates of SNP frequencies across varying levels of sequence coverage.

Evaluation of mitochondrial DNA coding region assays for increased discrimination in forensic analysis.

There is an increasing trend to use mitochondrial DNA (mtDNA) analysis in criminal investigations where only limited amounts of DNA are available. However, analysis of the mtDNA control region has the drawback of low discrimination power, due to the lack of recombination that results from uniparental (maternal) inheritance. As a strategy to increase discrimination, a number of typing assays detecting variation in the mitochondrial coding region have been developed. In this study, several of these assays are evaluated for their discriminatory capacity using data obtained from 495 complete Caucasian mtDNA sequences. In order to add a local geographic perspective to this evaluation, we have also sequenced and analysed the entire mtDNA from 20 individuals of Swedish origin. We find that the coding region assays are very useful for resolving sequences with identical HVI/HVII regions. The best-performing coding region assay was able to discriminate 46% of the resolvable sequences, compared to 20-30% for the other coding region assays we evaluated.

Rate variation between mitochondrial domains and adaptive evolution in humans.

It has been proposed that as a result of human adaptation to different climates, mitochondrial genes have been affected by natural selection. To further study the selective pressure on human mitochondrial DNA, we have analysed polymorphism at the gene, domain and nucleotide site level in four geographic regions. The ratio of non-synonymous relative to synonymous substitutions is elevated for ATP6 in North Asia, ND3 in Africa and cytb in Europe relative to other mitochondrial genes. In addition, non-synonymous substitutions appear nearly five times more frequently outside core functional domains as compared with within functional domains. Therefore, a large part of the rate variation between mitochondrial genes is explained by differences in the length of the functional domain relative to the length of the gene. The accumulation of non-synonymous substitutions in the ATP6 gene among sequences from North Asia has been a gradual process over many thousands of years, consistent with relaxed purifying selection rather than recent strong selective sweeps. It is not necessary to invoke adaptive selection to explain the pattern of polymorphism in mitochondrial genes, when the most parsimonious explanation is relaxed purifying selection and the action of genetic drift.

Mitochondrial genomics identifies major haplogroups in Aboriginal Australians.

We classified diversity in eight new complete mitochondrial genome sequences and 41 partial sequences from living Aboriginal Australians into five haplogroups. Haplogroup AuB belongs to global lineage M, and AuA, AuC, AuD, and AuE to N. Within N, we recognize subdivisions, assigning AuA to haplogroup S, AuD to haplogroup O, AuC to P4, and AuE to P8. On available evidence, (S)AuA and (M)AuB are widespread in Australia. (P4)AuC is found in the Riverine region of western New South Wales, and was identified by others in northern Australia. (O)AuD and (P8)AuE were clearly identified only from central Australia. Our eight Australian full mt genome sequences, combined with 20 others (Ingman and Gyllensten 2003 Genome Res. 13:1600-1606) and compared with full mt genome sequences from regions to the north that include Papua New Guinea, Malaya, and Andaman and Nicobar Islands, show that ancestral connections between regions are deep and limited to clustering at the level of the N and M macrohaplogroups. The Australian-specific distribution of the five haplogroups identified indicates genetic isolation over a long period. Ancestral connections within Australia are deeper than those reflected by known linguistic or culturally based affinities. Applying a coalescence analysis to a gene tree for the coding regions of the eight genomic sequences, we made estimates of time depth that support a continuity of presence for the descendants of a founding population already established by 40,000 years ago.

Mitochondrial genome variation and evolutionary history of Australian and New Guinean aborigines.

To study the evolutionary history of the Australian and New Guinean indigenous peoples, we analyzed 101 complete mitochondrial genomes including populations from Australia and New Guinea as well as from Africa, India, Europe, Asia, Melanesia, and Polynesia. The genetic diversity of the Australian mitochondrial sequences is remarkably high and is similar to that found across Asia. This is in contrast to the pattern seen in previously described Y-chromosome data where an Australia-specific haplotype was found at high frequency. The mitochondrial genome data indicate that Australia was colonized between 40 and 70 thousand years ago, either by a single migration from a heterogeneous source population or by multiple movements of smaller groups occurring over a period of time. Some Australian and New Guinea sequences form clades, suggesting the possibility of a joint colonization and/or admixture between the two regions.