Relaxation of selective constraint on dog mitochondrial DNA following domestication

Analysis of Canis mitochondrial DNA demonstrates high concordance between the control region and ATPase genes

Background

Phylogenetic studies of wild Canis species have relied heavily on the mitochondrial DNA control region (mtDNA CR) to infer species relationships and evolutionary lineages. Previous analyses of the CR provided evidence for a North American evolved eastern wolf (C. lycaon), that is more closely related to red wolves (C. rufus) and coyotes (C. latrans) than grey wolves (C. lupus). Eastern wolf origins, however, continue to be questioned. Therefore, we analyzed mtDNA from 89 wolves and coyotes across North America and Eurasia at 347 base pairs (bp) of the CR and 1067 bp that included the ATPase6 and ATPase8 genes. Phylogenies and divergence estimates were used to clarify the evolutionary history of eastern wolves, and regional comparisons of nonsynonomous to synonomous substitutions (dN/dS) at the ATPase6 and ATPase8 genes were used to elucidate the potential role of selection in shaping mtDNA geographic distribution.

Results

We found high concordance across analyses between the mtDNA regions studied. Both had a high percentage of variable sites (CR = 14.6%; ATP = 9.7%) and both phylogenies clustered eastern wolf haplotypes monophyletically within a North American evolved lineage apart from coyotes. Divergence estimates suggest the putative red wolf sequence is more closely related to coyotes (D_xyCR = 0.01982 ± 0.00494 SD; D_xyATP = 0.00332 ± 0.00097 SD) than the eastern wolf sequences (D_xyCR = 0.03047 ± 0.00664 SD; D_xyATP = 0.00931 ± 0.00205 SD). Neutrality tests on both genes were indicative of the population expansion of coyotes across eastern North America, and dN/dS ratios suggest a possible role for purifying selection in the evolution of North American lineages. dN/dS ratios were higher in European evolved lineages from northern climates compared to North American evolved lineages from temperate regions, but these differences were not statistically significant.

Conclusions

These results demonstrate high concordance between coding and non-coding regions of mtDNA, and provide further evidence that the eastern wolf possessed distinct mtDNA lineages prior to recent coyote introgression. Purifying selection may have influenced North American evolved Canis lineages, but detection of adaptive selection in response to climate is limited by the power of current statistical tests. Increased sampling and development of alternative analytical tools will be necessary to disentangle demographic history from processes of natural selection.

Impacts of recent cultivation on genetic diversity pattern of a medicinal plant, Scutellaria baicalensis (Lamiaceae)

Background

Cultivation of medicinal plants is not only a means for meeting current and future demands for large volume production of plant-based drug and herbal remedies, but also a means of relieving harvest pressure on wild populations. Scutellaria baicalensis Georgi (Huang-qin or Chinese skullcap) is a very important medicinal plant in China. Over the past several decades, wild resource of this species has suffered rapid declines and large-scale cultivation was initiated to meet the increasing demand for its root. However, the genetic impacts of recent cultivation on S. baicalensis have never been evaluated. In this study, the genetic diversity and genetic structure of 28 wild and 22 cultivated populations were estimated using three polymorphic chloroplast fragments. The objectives of this study are to provide baseline data for preserving genetic resource of S. baicalensis and to evaluate the genetic impacts of recent cultivation on medicinal plants, which may be instructive to future cultivation projects of traditional Chinese medicinal plants.

Results

Thirty-two haplotypes of S. baicalensis (HapA-Y and Hap1-7) were identified when three chloroplast spacers were combined. These haplotypes constituted a shallow gene tree without obvious clusters for cultivated populations, suggesting multiple origins of cultivated S. baicalensis. Cultivated populations (h_T = 0.832) maintained comparable genetic variation with wild populations (h_T = 0.888), indicating a slight genetic bottleneck due to multiple origins of cultivation. However, a substantial amount of rare alleles (10 out of 25 haplotypes within wild populations) lost during the course of S. baicalensis cultivation. The genetic differentiation for cultivated group (G_ST = 0.220) was significantly lower than that of wild group (G_ST = 0.701). Isolation by distance analysis showed that the effect of geographical isolation on genetic structure was significant in wild populations (r = 0.4346, P < 0.0010), but not in cultivated populations (r = 0.0599, P = 0.2710). These genetic distribution patterns suggest that a transient cultivation history and the extensive seed change among different geographical areas during the course of S. baicalensis cultivation.

Conclusions

Although cultivated S. baicalensis maintains comparable genetic diversity relative to wild populations, recent cultivation has still imposed profound impacts on genetic diversity patterns of the cultivated S. baicalensis populations, i.e., the loss of rare alleles and homogenization of cultivated populations. This study suggests that conservation-by-cultivation is an effective means for protecting genetic resources of S. baicalensis, however, the wild resources still need to be protected in situ and the evolutionary consequences of extensive seed exchange mediated by human being should be monitored carefully.

A complete mitochondrial genome of wheat (Triticum aestivum cv. Chinese Yumai), and fast evolving mitochondrial genes in higher plants

Plant mitochondrial genomes, encoding necessary proteins involved in the system of energy production, play an important role in the development and reproduction of the plant. They occupy a specific evolutionary pattern relative to their nuclear counterparts. Here, we determined the winter wheat (Triticum aestivum cv. Chinese Yumai) mitochondrial genome in a length of 452 and 526 bp by shotgun sequencing its BAC library. It contains 202 genes, including 35 known protein-coding genes, three rRNA and 17 tRNA genes, as well as 149 open reading frames (ORFs; greater than 300 bp in length). The sequence is almost identical to the previously reported sequence of the spring wheat (T. aestivum cv. Chinese Spring); we only identified seven SNPs (three transitions and four transversions) and 10 indels (insertions and deletions) between the two independently acquired sequences, and all variations were found in non-coding regions. This result confirmed the accuracy of the previously reported mitochondrial sequence of the Chinese Spring wheat. The nucleotide frequency and codon usage of wheat are common among the lineage of higher plant with a high AT-content of 58%. Molecular evolutionary analysis demonstrated that plant mitochondrial genomes evolved at different rates, which may correlate with substantial variations in metabolic rate and generation time among plant lineages. In addition, through the estimation of the ratio of non-synonymous to synonymous substitution rates between orthologous mitochondrion-encoded genes of higher plants, we found an accelerated evolutionary rate that seems to be the result of relaxed selection.

Feline non-repetitive mitochondrial DNA control region database for forensic evidence

The domestic cat is the one of the most popular pets throughout the world. A by-product of owning, interacting with, or being in a household with a cat is the transfer of shed fur to clothing or personal objects. As trace evidence, transferred cat fur is a relatively untapped resource for forensic scientists. Both phenotypic and genotypic characteristics can be obtained from cat fur, but databases for neither aspect exist. Because cats incessantly groom, cat fur may have nucleated cells, not only in the hair bulb, but also as epithelial cells on the hair shaft deposited during the grooming process, thereby generally providing material for DNA profiling. To effectively exploit cat hair as a resource, representative databases must be established. The current study evaluates 402 bp of the mtDNA control region (CR) from 1394 cats, including cats from 25 distinct worldwide populations and 26 breeds. Eighty-three percent of the cats are represented by 12 major mitotypes. An additional 8.0% are clearly derived from the major mitotypes. Unique sequences are found in 7.5% of the cats. The overall genetic diversity for this data set is 0.8813±0.0046 with a random match probability of 11.8%. This region of the cat mtDNA has discriminatory power suitable for forensic application worldwide.

Complex population structure in African village dogs and its implications for inferring dog domestication history

High genetic diversity of East Asian village dogs has recently been used to argue for an East Asian origin of the domestic dog. However, global village dog genetic diversity and the extent to which semiferal village dogs represent distinct, indigenous populations instead of admixtures of various dog breeds has not been quantified. Understanding these issues is critical to properly reconstructing the timing, number, and locations of dog domestication. To address these questions, we sampled 318 village dogs from 7 regions in Egypt, Uganda, and Namibia, measuring genetic diversity >680 bp of the mitochondrial D-loop, 300 SNPs, and 89 microsatellite markers. We also analyzed breed dogs, including putatively African breeds (Afghan hounds, Basenjis, Pharaoh hounds, Rhodesian ridgebacks, and Salukis), Puerto Rican street dogs, and mixed breed dogs from the United States. Village dogs from most African regions appear genetically distinct from non-native breed and mixed-breed dogs, although some individuals cluster genetically with Puerto Rican dogs or United States breed mixes instead of with neighboring village dogs. Thus, African village dogs are a mosaic of indigenous dogs descended from early migrants to Africa, and non-native, breed-admixed individuals. Among putatively African breeds, Pharaoh hounds, and Rhodesian ridgebacks clustered with non-native rather than indigenous African dogs, suggesting they have predominantly non-African origins. Surprisingly, we find similar mtDNA haplotype diversity in African and East Asian village dogs, potentially calling into question the hypothesis of an East Asian origin for dog domestication.

Contrast features of CpG islands in the promoter and other regions in the dog genome

The recent release of the domestic dog genome provides us with an ideal opportunity to investigate dog-specific genomic features. In this study, we performed a systematic analysis of CpG islands (CGIs), which are often considered gene markers, in the dog genome. Relative to the human and mouse genomes, the dog genome has a remarkably large number of CGIs and high CGI density, which is contributed by its noncoding sequences. Surprisingly, the dog genome has fewer CGIs associated with the promoter regions of genes than the human or the mouse. Further examination of functional features of dog-human-mouse homologous genes suggests that the dog might have undergone a faster erosion rate of promoter-associated CGIs than the human or mouse. Some genetic or genomic factors such as local recombination rate and karyotype may be related to the unique dog CGI features.

Mitochondrial genome DNA analysis of the domestic dog: identifying informative SNPs outside of the control region

While the mitochondrial control region has proven successful for human forensic evaluations by indicating ethnic origin, domestic dogs (Canis lupus familiaris) of seemingly unrelated breeds often form large groups based on identical control region sequences. In an attempt to break up these large haplotype groups, we have analyzed the remaining c. 15,484 base pairs of the canine mitochondrial genome for 79 dogs and used phylogenetic and population genetic methods to search for additional variability in the form of single nucleotide polymorphisms (SNPs). We have identified 356 SNPs and 65 haplotypes in the remainder of the mitochondrial genome excluding the control region. The exclusion capacity was found to be 0.018. The mitochondrial control region was also evaluated for the same 79 dogs. The signals from the different fragments do not conflict, but instead support one another and provide a larger fragment of DNA that can be analyzed as forensic evidence.

Relaxation of selective constraints on avian mitochondrial DNA following the degeneration of flight ability

The evolution of flight is the most important feature of birds, and this ability has helped them become one of the most successful groups of vertebrates. However, some species have independently lost their ability to fly. The degeneration of flight ability is a long process, and some species remain transitional locomotive models. Most of the energy required for locomotion is supplied by mitochondria via oxidative phosphorylation. Thus, rapidly flying birds should require a more energy efficient metabolism than weakly flying or flightless species. Therefore, we speculated that evolutionary constraints acted on the mtDNA of birds with different locomotive abilities. To test this hypothesis, we compared 76 complete avian mitochondrial genomes. Weakly flying and flightless birds, compared to rapidly flying birds, accumulated more nonsynonymous nucleotide substitutions relative to synonymous substitutions. Even after controlling for mutation rate, this trend remained significant. This finding was further tested for its generality by examining 214 complete mammalian mitochondrial genomes. The same as birds, a negative correlation was also found for the K(a)/K(s) ratio and locomotive speed. Our results demonstrated that, in addition to the previously described role for effective population size, functional constraints due to locomotion play an important role in the evolution of mtDNA.

The seminiferous epithelium cycle and its duration in different breeds of dog (Canis familiaris)

Testis structure and function in dogs are relatively poorly investigated. The aim of the present study was to carry out a comparative investigation of the stages of the seminiferous epithelium cycle and its duration in different breeds of dog. Fifty-six sexually mature dogs (mongrel, n = 12; pinscher, n = 12; beagle, n = 5; American pit bull, n = 9; poodle, n = 12; and Labrador retriever, n = 6) were analysed. Intratesticular injections of tritiated thymidine were given to determine the duration of spermatogenesis. Orchiectomy was performed at different time periods following injection (1 h, 2 and 4 weeks). Testis fragments were embedded in plastic and routinely prepared for histological and autoradiographic evaluations. Eight stages were characterized based on the acrosome system. Significant (P < 0.05) differences were found for the frequencies of the different stages characterized (except Stages V, VI and VIII), particularly for the mongrel. Stage IV (when spermiation occurs) was the most frequent in all six breeds (~25%), whereas Stages II and VIII were the least frequent (< 8%). Each spermatogenic cycle and the total duration of spermatogenesis lasted 13.73 +/- 0.03 and 61.9 +/- 0.14 days, respectively, for the mongrel, poodle, pinscher, beagle, and Labrador retriever. These values were approximately 10% lower (P < 0.03) for the American pit bull (12.55 +/- 0.26 and 56.5 +/- 1.17 days, respectively). To our knowledge, this is the first comprehensive study to perform a careful investigation of stage frequencies and seminiferous epithelium cycle duration in this very important domestic species.

Origin and phylogenetic analysis of Tibetan Mastiff based on the mitochondrial DNA sequence

At present, the Tibetan Mastiff is the oldest and most ferocious dog in the world. However, the origin of the Tibetan Mastiff and its phylogenetic relationship with other large breed dogs such as Saint Bernard are unclear. In this study, the primers were designed according to the mitochondrial genome sequence of the domestic dog, and the 2,525 bp mitochondrial sequence, containing the whole sequence of Cytochrome b, tRNA-Thr, tRNA-Pro, and control region of the Tibetan Mastiff, was obtained. Using grey wolves and coyotes as outgroups, the Tibetan Mastiff and 12 breeds of domestic dogs were analyzed in phylogenesis. Tibetan Mastiff, domestic dog breeds, and grey wolves were clustered into a group and coyotes were clustered in a group separately. This indicated that the Tibetan Mastiff and the other domestic dogs originated from the grey wolf, and the Tibetan Mastiff belonged to Carnivora, Canidae, Canis, Canis lupus, Canis lupus familiaris on the animal taxonomy. In domestic dogs, the middle and small breed dogs were clustered at first; German Sheepdog, Swedish Elkhound, and Black Russian Terrier were clustered into one group, and the Tibetan Mastiff, Old English Sheepdog, Leonberger, and Saint Bernard were clustered in another group. This confirmed the viewpoint that many of the famous large breed dogs worldwide such as Saint Bernard possibly had the blood lineage of the Tibetan Mastiff, based on the molecular data. According to the substitution rate, we concluded that the approximate divergence time between Tibetan Mastiff and grey wolf was 58,000 years before the present (YBP), and the approximate divergence time between other domestic dogs and grey wolf was 42,000 YBP, demonstrating that the time of origin of the Tibetan Mastiff was earlier than that of the other domestic dogs.

Monophyletic origin of domestic bactrian camel (Camelus bactrianus) and its evolutionary relationship with the extant wild camel (Camelus bactrianus ferus)

The evolutionary relationship between the domestic bactrian camel and the extant wild two-humped camel and the factual origin of the domestic bactrian camel remain elusive. We determined the sequence of mitochondrial cytb gene from 21 camel samples, including 18 domestic camels (three Camelus bactrianus xinjiang, three Camelus bactrianus sunite, three Camelus bactrianus alashan, three Camelus bactrianus red, three Camelus bactrianus brown and three Camelus bactrianus normal) and three wild camels (Camelus bactrianus ferus). Our phylogenetic analyses revealed that the extant wild two-humped camel may not share a common ancestor with the domestic bactrian camel and they are not the same subspecies at least in their maternal origins. Molecular clock analysis based on complete mitochondrial genome sequences indicated that the sub-speciation of the two lineages had begun in the early Pleistocene, about 0.7 million years ago. According to the archaeological dating of the earliest known two-humped camel domestication (5000-6000 years ago), we could conclude that the extant wild camel is a separate lineage but not the direct progenitor of the domestic bactrian camel. Further phylogenetic analysis suggested that the bactrian camel appeared monophyletic in evolutionary origin and that the domestic bactrian camel could originate from a single wild population. The data presented here show how conservation strategies should be implemented to protect the critically endangered wild camel, as it is the last extant form of the wild tribe Camelina.

Molecular evolution of the antiretroviral TRIM5 gene

In 2004, the first report of TRIM5alpha as a cellular antiretroviral factor triggered intense interest among virologists, particularly because some primate orthologs of TRIM5alpha have activity against HIV. Since that time, a complex and eventful evolutionary history of the TRIM5 locus has emerged. A review of the TRIM5 literature constitutes a veritable compendium of evolutionary phenomena, including elevated rates of nonsynonymous substitution, divergence in subdomains due to short insertions and deletions, expansions and contractions in gene copy number, pseudogenization, balanced polymorphism, trans-species polymorphism, convergent evolution, and the acquisition of new domains by exon capture. Unlike most genes, whose history is dominated by long periods of purifying selection interspersed with rare instances of genetic innovation, analysis of restriction factor loci is likely to be complicated by the unpredictable and more-or-less constant influence of positive selection. In this regard, the molecular evolution and population genetics of restriction factor loci most closely resemble patterns that have been documented for immunity genes, such as class I and II MHC genes, whose products interact directly with microbial targets. While the antiretroviral activity encoded by TRIM5 provides plausible mechanistic hypotheses for these unusual evolutionary observations, evolutionary analyses have reciprocated by providing significant insights into the structure and function of the TRIM5alpha protein. Many of the lessons learned from TRIM5 should be applicable to the study of other restriction factor loci, and molecular evolutionary analysis could facilitate the discovery of new antiviral factors, particularly among the many TRIM genes whose functions remain as yet unidentified.

Quantification of adaptive evolution of genes expressed in avian brain and the population size effect on the efficacy of selection

Whether protein evolution is mainly due to fixation of beneficial alleles by positive selection or to random genetic drift has remained a contentious issue over the years. Here, we use two genomewide polymorphism data sets collected from chicken populations, together with divergence data from >5,000 chicken-zebra finch gene orthologs expressed in brain, to assess the amount of adaptive evolution in protein-coding genes of birds. First, we show that estimates of the fixation index (FI, the ratio of fixed nonsynonymous-to-synonymous changes over the ratio of the corresponding polymorphisms) are highly dependent on the character of the underlying data sets. Second, by using polymorphism data from high-frequency alleles, to avoid the confounding effect of slightly deleterious mutations segregating at low frequency, we estimate that about 20% of amino acid changes have been brought to fixation through positive selection during avian evolution. This estimate is intermediate to that obtained in humans (lower) and flies as well as bacteria (higher), and is consistent with population genetics theory that stipulates a positive relationship between the efficiency of selection and the effective population size. Further, by comparing the FIs for common and all alleles, we estimate that approximately 20% of nonsynonymous variation segregating in chicken populations represent slightly deleterious mutations, which is less than in Drosophila. Overall, these results highlight the link between the effective population size and positive as well as negative selection.

A complete mitochondrial genome sequence of the wild two-humped camel (Camelus bactrianus ferus): an evolutionary history of camelidae

Background

The family Camelidae that evolved in North America during the Eocene survived with two distinct tribes, Camelini and Lamini. To investigate the evolutionary relationship between them and to further understand the evolutionary history of this family, we determined the complete mitochondrial genome sequence of the wild two-humped camel (Camelus bactrianus ferus), the only wild survivor of the Old World camel.

Results

The mitochondrial genome sequence (16,680 bp) from C. bactrianus ferus contains 13 protein-coding, two rRNA, and 22 tRNA genes as well as a typical control region; this basic structure is shared by all metazoan mitochondrial genomes. Its protein-coding region exhibits codon usage common to all mammals and possesses the three cryptic stop codons shared by all vertebrates. C. bactrianus ferus together with the rest of mammalian species do not share a triplet nucleotide insertion (GCC) that encodes a proline residue found only in the nd1 gene of the New World camelid Lama pacos. This lineage-specific insertion in the L. pacos mtDNA occurred after the split between the Old and New World camelids suggests that it may have functional implication since a proline insertion in a protein backbone usually alters protein conformation significantly, and nd1 gene has not been seen as polymorphic as the rest of ND family genes among camelids. Our phylogenetic study based on complete mitochondrial genomes excluding the control region suggested that the divergence of the two tribes may occur in the early Miocene; it is much earlier than what was deduced from the fossil record (11 million years). An evolutionary history reconstructed for the family Camelidae based on cytb sequences suggested that the split of bactrian camel and dromedary may have occurred in North America before the tribe Camelini migrated from North America to Asia.

Conclusion

Molecular clock analysis of complete mitochondrial genomes from C. bactrianus ferus and L. pacos suggested that the two tribes diverged from their common ancestor about 25 million years ago, much earlier than what was predicted based on fossil records.

Origin of mitochondrial DNA diversity of domestic yaks

Background

The domestication of plants and animals was extremely important anthropologically. Previous studies have revealed a general tendency for populations of livestock species to include deeply divergent maternal lineages, indicating that they were domesticated in multiple, independent events from genetically discrete wild populations. However, in water buffalo, there are suggestions that a similar deep maternal bifurcation may have originated from a single population. These hypotheses have rarely been rigorously tested because of a lack of sufficient wild samples. To investigate the origin of the domestic yak (Poephagus grunnies), we analyzed 637 bp of maternal inherited mtDNA from 13 wild yaks (including eight wild yaks from a small population in west Qinghai) and 250 domesticated yaks from major herding regions.

Results

The domestic yak populations had two deeply divergent phylogenetic groups with a divergence time of > 100,000 yrs BP. We here show that haplotypes clustering with two deeply divergent maternal lineages in domesticated yaks occur in a single, small, wild population. This finding suggests that all domestic yaks are derived from a single wild gene pool. However, there is no clear correlation of the mtDNA phylogenetic clades and the 10 morphological types of sampled yaks indicating that the latter diversified recently. Relatively high diversity was found in Qinghai and Tibet around the current wild distribution, in accordance with previous suggestions that the earliest domestications occurred in this region. Conventional molecular clock estimation led to an unrealistic early dating of the start of the domestication. However, Bayesian estimation of the coalescence time allowing a relaxation of the mutation rate are better in agreement with a domestication during the Holocene as supported by archeological records.

Conclusion

The information gathered here and the previous studies of other animals show that the demographic histories of domestication of livestock species were highly diverse despite the common general feature of deeply divergent maternal lineages. The results further suggest that domestication of local wild prey ungulate animals was a common occurrence during the development of human civilization following the postglacial colonization in different locations of the world, including the high, arid Qinghai-Tibetan Plateau.