Standardisation of nomenclature for dog mtDNA D-loop: a prerequisite for launching a Canis familiaris database

Characterization of Nuclear Mitochondrial Insertions in Canine Genome Assemblies

BACKGROUND

The presence of mitochondrial sequences in the nuclear genome (Numts) confounds analyses of mitochondrial sequence variation, and is a potential source of false positives in disease studies. To improve the analysis of mitochondrial variation in canines, we completed a systematic assessment of Numt content across genome assemblies, canine populations and the carnivore lineage.

RESULTS

Centering our analysis on the UU_Cfam_GSD_1.0/canFam4/Mischka assembly, a commonly used reference in dog genetic variation studies, we found a total of 321 Numts located throughout the nuclear genome and encompassing the entire sequence of the mitochondria. A comparison with 14 canine genome assemblies identified 63 Numts with presence-absence dimorphism among dogs, wolves, and a coyote. Furthermore, a subset of Numts were maintained across carnivore evolutionary time (arctic fox, polar bear, cat), with eight sequences likely more than 10 million years old, and shared with the domestic cat. On a population level, using structural variant data from the Dog10K Consortium for 1879 dogs and wolves, we identified 11 Numts that are absent in at least one sample, as well as 53 Numts that are absent from the Mischka assembly.

CONCLUSIONS

We highlight scenarios where the presence of Numts is a potentially confounding factor and provide an annotation of these sequences in canine genome assemblies. This resource will aid the identification and interpretation of polymorphisms in both somatic and germline mitochondrial studies in canines.

Genomic Diversity and Runs of Homozygosity in Bernese Mountain Dogs

Bernese mountain dogs are a large dog breed formed in the early 1900s in Switzerland. While originally farm dogs that were used for pulling carts, guarding, and driving cattle, today they are considered multi-purpose companion and family dogs. The breed is predisposed to several complex diseases, such as histiocytic sarcoma, degenerative myelopathy, or hip dysplasia. Using whole-genome sequencing (WGS) data, we assessed the genomic architecture of 33 unrelated dogs from four countries: France, Sweden, Switzerland, and the United States. Analysis of runs of homozygosity (ROH) identified 12,643 ROH with an average length of 2.29 Mb and an average inbreeding coefficient of 0.395. Multidimensional scaling analysis of the genetic relatedness revealed limited clustering of European versus USA dogs, suggesting exchanges of breeding stock between continents. Furthermore, only two mtDNA haplotypes were detected in the 33 studied dogs, both of which are widespread throughout multiple dog breeds. WGS-based ROH analyses revealed several fixed or nearly fixed regions harboring discreet morphological trait-associated as well as disease-associated genetic variants. Several genes involved in the regulation of immune cells were found in the ROH shared by all dogs, which is notable in the context of the breed’s strong predisposition to hematopoietic cancers. High levels of inbreeding and relatedness, strongly exaggerated in the last 30 years, have likely led to the high prevalence of specific genetic disorders in this breed.

Evaluation of the tRNA-Leu (UUR) gene haplotype profile observed in canine mammary gland tumours based on comparative analysis with the MT-TL1 human gene

The aetiology and pathogenesis of many canine tumours are likely to be similar to cancers found in humans. This study aimed to evaluate a plausible link between changes in the tRNA-Leu (UUR) gene and the carcinogenesis process in dogs with mammary gland tumours. The whole mitochondrial DNA (mtDNA) isolated from blood and tumour tissues of 13 dogs with malignant mammary gland tumours was sequenced. The present work is the first report showing that some polymorphisms might occur at the corresponding positions in the human and canine mtDNA genome, which in turn may provoke similar deleterious effects. The homology between the human MT-TL1 and canine tRNA-Leu (UUR) genes was 84%. After resequencing of the whole mitochondrial DNA genome with the use of the NGS technology, two polymorphisms in two haplotypes were identified: m.2683G>A (observed in 18 out of 27 samples) and m.2678_2679insG (27 out of 27 samples). The m.2683G>A polymorphism corresponded to a deleterious change at m.3243A>G, which is linked with MELAS (Mitochondrial Encephalomyopathy, Lactic Acidosis, Stroke-like episodes) syndrome and with different types of cancers in humans as well. The comparative analysis of MT-TL1 and tRNA-Leu (UUR) led us to hypothesise that the m.2678_2679insG and m.2683G>A polymorphisms might influence the dog’s condition and might be linked with tumourigenesis, as observed in humans.

Genomic diversity and population structure of the Leonberger dog breed

BACKGROUND

Leonberger is a giant dog breed formed in the 1850s in Germany. Its post-World War II popularity has resulted in a current global population of ~ 30,000 dogs. The breed has predispositions to neurodegenerative disorders and cancer, which is likely due in large part to limited genetic diversity. However, to date there is no scientific literature on the overall demography and genomic architecture of this breed.

RESULTS

We assessed extensive pedigree records, SNP array genotype data, and whole-genome sequences (WGS) on 142,072, 1203 and 39 Leonberger dogs, respectively. Pedigree analyses identified 22 founder animals and revealed an apparent popular sire effect. The average pedigree-based inbreeding coefficient of 0.29 and average kinship of 0.31 show a dramatic loss of genetic diversity. The observed average life span decreased over time from 9.4 years in 1989 to 7.7 years in 2004. A global health survey confirmed a high prevalence of cancer and neurological disorders. Analysis of SNP-based runs of homozygosity (ROH) identified 125,653 ROH with an average length of 5.88 Mb, and confirmed an average inbreeding coefficient of 0.28. Genome-wide filtering of the WGS data revealed 28 non-protein-changing variants that were present in all Leonberger individuals and a list of 22 potentially pathogenic variants for neurological disorders of which 50% occurred only in Leonbergers and 50% occurred rarely in other breeds. Furthermore, one of the two mtDNA haplogroups detected was present in one dog only.

CONCLUSIONS

The increasing size of the Leonberger population has been accompanied by a considerable loss of genetic diversity after the bottleneck that occurred in the 1940s due to the intensive use of popular sires resulting in high levels of inbreeding. This might explain the high prevalence of certain disorders; however, genomic data provide no evidence for fixed coding variants that explain these predispositions. The list of candidate causative variants for polyneuropathy needs to be further evaluated. Preserving the current genetic diversity is possible by increasing the number of individuals for breeding while restricting the number of litters per sire/dam. In addition, outcrossing would help optimize long-term genetic diversity and contribute to the sustainability and health of the population.

Who killed my dog? Use of forensic genetics to investigate an enigmatic case

Genetic testing of animal biological material has become a valuable tool in forensic investigations, and it is successfully used to identify unknown crime perpetrators, to unmask food frauds, or to clarify cases of animal attacks on humans or other animals. When DNA profiling is not possible due to inadequate amounts of nuclear DNA, mitochondrial DNA (mtDNA) testing is the only viable alternative, as in the case of shed hair samples. In this case, a dog was allegedly killed by wild animals while being hosted in a boarding house. Extraneous hair fragments recovered from the dog’s mouth and paws were subjected to genetic analysis: the cytochrome b gene located on mtDNA was amplified and sequenced in order to determine the species responsible for the killing. The mtDNA analysis provided evidence that the dog was killed by other dogs, thus unmasking a false wild animal attack and putting the case in an entirely different perspective.

Canis mtDNA HV1 database: a web-based tool for collecting and surveying Canis mtDNA HV1 haplotype in public database

BACKGROUND

Canine and wolf mitochondrial DNA haplotypes, which can be used for forensic or phylogenetic analyses, have been defined in various schemes depending on the region analyzed. In recent studies, the 582 bp fragment of the HV1 region is most commonly used. 317 different canine HV1 haplotypes have been reported in the rapidly growing public database GenBank. These reported haplotypes contain several inconsistencies in their haplotype information. To overcome this issue, we have developed a Canis mtDNA HV1 database. This database collects data on the HV1 582 bp region in dog mitochondrial DNA from the GenBank to screen and correct the inconsistencies. It also supports users in detection of new novel mutation profiles and assignment of new haplotypes.

DESCRIPTION

The Canis mtDNA HV1 database (CHD) contains 5567 nucleotide entries originating from 15 subspecies in the species Canis lupus. Of these entries, 3646 were haplotypes and grouped into 804 distinct sequences. 319 sequences were recognized as previously assigned haplotypes, while the remaining 485 sequences had new mutation profiles and were marked as new haplotype candidates awaiting further analysis for haplotype assignment. Of the 3646 nucleotide entries, only 414 were annotated with correct haplotype information, while 3232 had insufficient or lacked haplotype information and were corrected or modified before storing in the CHD. The CHD can be accessed at http://chd.vnbiology.com . It provides sequences, haplotype information, and a web-based tool for mtDNA HV1 haplotyping. The CHD is updated monthly and supplies all data for download.

CONCLUSIONS

The Canis mtDNA HV1 database contains information about canine mitochondrial DNA HV1 sequences with reconciled annotation. It serves as a tool for detection of inconsistencies in GenBank and helps identifying new HV1 haplotypes. Thus, it supports the scientific community in naming new HV1 haplotypes and to reconcile existing annotation of HV1 582 bp sequences.

Phylogeography of the Koala, (Phascolarctos cinereus), and Harmonising Data to Inform Conservation

The Australian continent exhibits complex biogeographic patterns but studies of the impacts of Pleistocene climatic oscillation on the mesic environments of the Southern Hemisphere are limited. The koala (Phascolarctos cinereus), one of Australia’s most iconic species, was historically widely distributed throughout much of eastern Australia but currently represents a complex conservation challenge. To better understand the challenges to koala genetic health, we assessed the phylogeographic history of the koala. Variation in the maternally inherited mitochondrial DNA (mtDNA) Control Region (CR) was examined in 662 koalas sampled throughout their distribution. In addition, koala CR haplotypes accessioned to Genbank were evaluated and consolidated. A total of 53 unique CR haplotypes have been isolated from koalas to date (including 15 haplotypes novel to this study). The relationships among koala CR haplotypes were indicative of a single Evolutionary Significant Unit and do not support the recognition of subspecies, but were separated into four weakly differentiated lineages which correspond to three geographic clusters: a central lineage, a southern lineage and two northern lineages co-occurring north of Brisbane. The three geographic clusters were separated by known Pleistocene biogeographic barriers: the Brisbane River Valley and Clarence River Valley, although there was evidence of mixing amongst clusters. While there is evidence for historical connectivity, current koala populations exhibit greater structure, suggesting habitat fragmentation may have restricted female-mediated gene flow. Since mtDNA data informs conservation planning, we provide a summary of existing CR haplotypes, standardise nomenclature and make recommendations for future studies to harmonise existing datasets. This holistic approach is critical to ensuring management is effective and small scale local population studies can be integrated into a wider species context.

Defect of the mitochondrial DNA hypervariable region as a risk factor for canine mammary tumour

The aim of this study was to identify mutations in the hypervariable region of mitochondrial DNA in canine mammary tumours and to determine their association with the process of neoplastic transformation. A total of 93 biological samples, including blood as well as normal and neoplastic tissue samples from 31 dogs with diagnosed malignant canine mammary tumours were analysed. DNA extraction, amplification and sequencing of the D-loop as well as bioinformatic and statistical analyses were performed. In the mitochondrial D-loop sequence, 26 polymorphic loci and 5 mutations were identified. For the first time, D-loop length heteroplasmy was detected in dogs with mammary tumours. The malignancy grade exerted no effect on the presence of nucleotide changes. A statistically significant association between the presence of mutations and polymorphisms and the size of dogs was demonstrated. The 100% frequency of length heteroplasmy may imply that this is a hotspot mutation of canine mammary tumour.

Identification of additional mitochondrial DNA mutations in canine mast cell tumours

BACKGROUND

Research has revealed the presence of somatic mutations in mitochondrial DNA (mtDNA) of certain types of tumours. As this has not been studied for canine mast cell tumours, the aim of this study was to identify mutations in the hypervariable region of mtDNA in mast cell tumours in dogs and determine their association with the process of neoplastic transformation.

RESULTS

Samples from 17 dogs with histopathologically confirmed mast cell tumours were analysed. The samples consisted of tumour tissues (n = 17), normal tissues (n = 17), and blood (n = 17). Amplicons of the displacement loop (D-loop) were sequenced and the obtained nucleotide sequences were subjected to bioinformatics analyses. Somatic mutations were detected in seven positions of the D-loop nucleotide sequences in 47 % of the dogs, while polymorphisms were identified in 94 % of the dogs. Most of these changes were homoplasmic, while heteroplasmy was detected in two individuals. Six new haplotypes were established as being characteristic for canine mast cell tumours. There was no association between the presence of the mutations and sex, haplotype, or malignancy grade assessed in 3 and 2-grade scales.

CONCLUSIONS

Differences in the frequency of somatic mutations imply their direct association with the neoplastic transformation. However, their functional consequences and clinical significance are not clear. The mutations may be used for diagnosis and prognosis of canine mast cell tumours in the future.

Complete mitochondrial genome database and standardized classification system for Canis lupus familiaris

To contribute to the complete mitogenome database of the species Canis lupus familiaris and shed more light on its origin, we have sequenced mitochondrial genomes of 120 modern dogs from worldwide populations. Together with all the previously published mitogenome sequences of acceptable quality, we have reconstructed a global phylogenetic tree of 555 C. l. familiaris mitogenomes and standardized haplogroup nomenclature. The phylogenetic tree presented here and available online at http://clf.mtdna.tree.cm.umk.pl/ could be further used by forensic and evolutionary geneticists as well cynologists, for data quality control and unambiguous haplogroup classification. Our in-depth phylogeographic analysis of all C. l. familiaris mitogenomes confirmed that domestic dogs may have originated in East Asia during the Mesolithic and Upper Paleolithic time periods and started to expand to other parts of the world during Neolithic times.

Mitochondrial DNA haplogroup phylogeny of the dog: Proposal for a cladistic nomenclature

Canis lupus familiaris mitochondrial DNA analysis has increased in recent years, not only for the purpose of deciphering dog domestication but also for forensic genetic studies or breed characterization. The resultant accumulation of data has increased the need for a normalized and phylogenetic-based nomenclature like those provided for human maternal lineages. Although a standardized classification has been proposed, haplotype names within clades have been assigned gradually without considering the evolutionary history of dog mtDNA. Moreover, this classification is based only on the D-loop region, proven to be insufficient for phylogenetic purposes due to its high number of recurrent mutations and the lack of relevant information present in the coding region. In this study, we design 1) a refined mtDNA cladistic nomenclature from a phylogenetic tree based on complete sequences, classifying dog maternal lineages into haplogroups defined by specific diagnostic mutations, and 2) a coding region SNP analysis that allows a more accurate classification into haplogroups when combined with D-loop sequencing, thus improving the phylogenetic information obtained in dog mitochondrial DNA studies.

Local mitochondrial DNA haplotype databases needed for domestic dog populations that have experienced founder effect

Biological material from pets is often collected as evidence from crime scenes. Due to sample type and quality, mitochondrial DNA (mtDNA) is frequently evaluated to identify the potential contributor. MtDNA has a lower discriminatory power than nuclear DNA with multiple individuals in a population potentially carrying the same mtDNA sequence, or haplotype. The frequency distribution of mtDNA haplotypes in a population must be known in order to determine the evidentiary value of a match between crime scene evidence and the potential contributor of the biological material. This is especially important in geographic areas that include remote and/or isolated populations where founder effect may have lead to a decrease in genetic diversity and a non-random distribution of haplotypes relative to the population at large. Here we compared the haplotype diversity in dogs from the noncontiguous states of Alaska and Hawaii relative to the contiguous United States (US). We report a greater proportion of dogs carrying an A haplotype in Alaska relative to any other US population. Significant variation in the distribution of haplotype frequencies was discovered when comparing the haplotype diversity of dogs in Hawaii to that of the continental US. Each of these regions exhibits reduced genetic diversity relative to the contiguous US, likely due to founder effect. We recommend that specific databases be created to accurately represent the mitochondrial haplotype diversity in these remote areas. Furthermore, our work demonstrates the importance of local surveys for populations that may have experienced found effect.

Length heteroplasmy of the polyC-polyT-polyC stretch in the dog mtDNA control region

Previously, the mitochondrial control region of 214 Belgian dogs was sequenced. Analysis of this data indicated length heteroplasmy of the polyT stretch in the polyC-polyT-polyC stretch from positions 16661 to 16674. Nine polyC-polyT-polyC haplotype combinations were observed, consisting of seven major haplotypes (highest signal intensity) combined with minor haplotypes (lower signal intensity) one T shorter than the major haplotype in all but three dogs. The longer the polyT stretch, the smaller was the difference in signal intensity between the major and minor haplotype peaks. Additional sequencing, cloning, and PCR trap experiments were performed to further study the intra-individual variation of this mitochondrial DNA (mtDNA) region. Cloning experiments demonstrated that the proportion of clones displaying the minor haplotypes also increased with the length of the polyT stretch. Clone amplification showed that in vitro polymerase errors might contribute to the length heteroplasmy of polyT stretches with at least 10 Ts. Although major and minor polyC-polyT-polyC haplotypes did not differ intra-individually within and between tissues in this study, interpretation of polyT stretch variation should be handled with care in forensic casework.

Coding region SNP analysis to enhance dog mtDNA discrimination power in forensic casework

The high population frequencies of three control region haplotypes contribute to the low discrimination power of the dog mtDNA control region. It also diminishes the evidential power of a match with one of these haplotypes in forensic casework. A mitochondrial genome study of 214 Belgian dogs suggested 26 polymorphic coding region sites that successfully resolved dogs with the three most frequent control region haplotypes. In this study, three SNP assays were developed to determine the identity of the 26 informative sites. The control region of 132 newly sampled dogs was sequenced and added to the study of 214 dogs. The assays were applied to 58 dogs of the haplotypes of interest, which confirmed their suitability for enhancing dog mtDNA discrimination power. In the Belgian population study of 346 dogs, the set of 26 sites divided the dogs into 25 clusters of mtGenome sequences with substantially lower population frequency estimates than their control region sequences. In case of a match with one of the three control region haplotypes, using these three SNP assays in conjunction with control region sequencing would augment the exclusion probability of dog mtDNA analysis from 92.9% to 97.0%.

Reconstructing the colonization history of lost wolf lineages by the analysis of the mitochondrial genome

The grey wolves (Canis lupus) originally inhabited major parts of the Northern hemisphere, but many local populations became extinct. Two lineages of wolves in Japan, namely, Japanese or Honshu (C. l. hodophilax) and Ezo or Hokkaido (C. l. hattai) wolves, rapidly went extinct between 100 and 120years ago. Here we analyse the complete mitochondrial genome sequences from ancient specimens and reconstruct the colonization history of the two extinct subspecies. We show a unique status of Japanese wolves in wolf phylogeny, suggesting their long time separation from other grey wolf populations. Japanese wolves appeared to have colonized the Japanese archipelago in the Late Pleistocene (ca. 25,000-125,000years ago). By contrast, Ezo wolves, which are clearly separated from Japanese wolves in phylogeny, are likely to have arrived at Japan relatively recently (<14,000years ago). Interestingly, their colonization history to Japan tallies well with the dynamics of wolf populations in Europe and America during the last several millennia. Our analyses suggest that at least several thousands of wolves once inhabited in the Japanese archipelago. Our analyses also show that an enigmatic clade of domestic dogs is likely to have originated from rare admixture events between male dogs and female Japanese wolves.

Forensic informativity of ~3000 bp of coding sequence of domestic dog mtDNA

The discriminatory power of the noncoding control region (CR) of domestic dog mitochondrial DNA alone is relatively low. The extent to which the discriminatory power could be increased by analyzing additional highly variable coding regions of the mitochondrial genome (mtGenome) was therefore investigated. Genetic variability across the mtGenome was evaluated by phylogenetic analysis, and the three most variable ~1 kb coding regions identified. We then sampled 100 Swedish dogs to represent breeds in accordance with their frequency in the Swedish population. A previously published dataset of 59 dog mtGenomes collected in the United States was also analyzed. Inclusion of the three coding regions increased the exclusion capacity considerably for the Swedish sample, from 0.920 for the CR alone to 0.964 for all four regions. The number of mtDNA types among all 159 dogs increased from 41 to 72, the four most frequent CR haplotypes being resolved into 22 different haplotypes.

Dog mitochondrial genome sequencing to enhance dog mtDNA discrimination power in forensic casework

A Belgian dog population sample and several population studies worldwide have confirmed that only a limited number of mtDNA control region haplotypes is observed in the majority of dogs. The high population frequency of these haplotypes negatively impacts both the exclusion probability of dog mtDNA analysis and the evidential value of a match with one of these haplotypes in casework. Variation within the mtDNA coding region was explored to improve the discrimination power of dog mtDNA analysis. In the current study, the entire mitochondrial genome of 161 dogs was sequenced applying a quality assured strategy and resulted in a total of 119 different mitochondrial genome sequences. Our research was focused on those dogs with the six most common control region haplotypes from a previous Belgian population study. We identified 33 informative SNPs that successfully divide the six most common control region haplotypes into 32 clusters of mitochondrial genome sequences. Determining the identity of these 33 polymorphic sites in addition to control region sequencing in case of a match with one of these 6 control region haplotypes could augment the exclusion probability of forensic dog mtDNA analysis from 92.5% to 97.5%.

Population structure and genetic differentiation of livestock guard dog breeds from the Western Balkans

Livestock guard dog (LGD) breeds from the Western Balkans are a good example of how complex genetic diversity pattern observed in dog breeds has been shaped by transition in dog breeding practices. Despite their common geographical origin and relatively recent formal recognition as separate breeds, the Karst Shepherd, Sarplaninac and Tornjak show distinct population dynamics, assessed by pedigree, microsatellite and mtDNA data. We genotyped 493 dogs belonging to five dog breeds using a set of 18 microsatellite markers and sequenced mtDNA from 94 dogs from these breeds. Different demographic histories of the Karst Shepherd and Tornjak breeds are reflected in the pedigree data with the former breed having more unbalanced contributions of major ancestors and a realized effective population size of less than 20 animals. The highest allelic richness was found in Sarplaninac (5.94), followed by Tornjak (5.72), whereas Karst Shepherd dogs exhibited the lowest allelic richness (3.33). Similarly, the highest mtDNA haplotype diversity was found in Sarplaninac, followed by Tornjak and Karst Shepherd, where only one haplotype was found. Based on FST differentiation values and high percentages of animals correctly assigned, all breeds can be considered genetically distinct. However, using microsatellite data, common ancestry between the Karst Shepherd and Sarplaninac could not be reconstructed, despite pedigree and mtDNA evidence of their historical admixture. Using neighbour-joining, STRUCTURE or DAPC methods, Sarplaninac and Caucasian Shepherd breeds could not be separated and additionally showed close proximity in the NeighborNet tree. STRUCTURE analysis of the Tornjak breed demonstrated substructuring, which needs further investigation. Altogether, results of this study show that the official separation of these dog breeds strongly affected the resolution of genetic differentiation and thus suggest that the relationships between breeds are not only determined by breed relatedness, but in small populations even more importantly by stochastic effects.

Reviewing population studies for forensic purposes: Dog mitochondrial DNA

The identification of dog hair through mtDNA analysis has become increasingly important in the last 15 years, as it can provide associative evidence connecting victims and suspects. The evidential value of an mtDNA match between dog hair and its potential donor is determined by the random match probability of the haplotype. This probability is based on the haplotype’s population frequency estimate. Consequently, implementing a population study representative of the population relevant to the forensic case is vital to the correct evaluation of the evidence. This paper reviews numerous published dog mtDNA studies and shows that many of these studies vary widely in sampling strategies and data quality. Therefore, several features influencing the representativeness of a population sample are discussed. Moreover, recommendations are provided on how to set up a dog mtDNA population study and how to decide whether or not to include published data. This review emphasizes the need for improved dog mtDNA population data for forensic purposes, including targeting the entire mitochondrial genome. In particular, the creation of a publicly available database of qualitative dog mtDNA population studies would improve the genetic analysis of dog traces in forensic casework.

Validation of two canine STR multiplex-assays following the ISFG recommendations for non-human DNA analysis

To gain general acceptance forensic DNA testing in animals needs to improve standardization of analysis methods and data interpretation. Recently, the International Society of Forensic Genetics (ISFG) took particular care of this topic by publishing recommendations for forensic non-human DNA analysis following the successful example of human DNA analysis in order to provide a basis for harmonization of the still existing inter-laboratory variability. By following these recommendations we demonstrate the performance of two short tandem repeat (STR) multiplexes for forensic identity testing of canine biological material. Thirteen STRs and two sex-specific markers were selected and validated according to the ISFG guidelines. Population genetic parameters were calculated based on 295 dog samples collected in Austria (124) and Germany (171). A repeat-based nomenclature of the mainly tetrameric STRs and corresponding allelic ladders are presented. All 146 different alleles included in the ladders were sequenced for correct allele calling. Additionally, a canine cell line (DH82-D3167) was evaluated as standard reference material.

Evaluating purifying selection in the mitochondrial DNA of various mammalian species

Mitochondrial DNA (mtDNA), the circular DNA molecule inside the mitochondria of all eukaryotic cells, has been shown to be under the effect of purifying selection in several species. Traditional testing of purifying selection has been based simply on ratios of nonsynonymous to synonymous mutations, without considering the relative age of each mutation, which can be determined by phylogenetic analysis of this non-recombining molecule. The incorporation of a mutation time-ordering from phylogeny and of predicted pathogenicity scores for nonsynonymous mutations allow a quantitative evaluation of the effects of purifying selection in human mtDNA. Here, by using this additional information, we show that purifying selection undoubtedly acts upon the mtDNA of other mammalian species/genera, namely Bos sp., Canis lupus, Mus musculus, Orcinus orca, Pan sp. and Sus scrofa. The effects of purifying selection were comparable in all species, leading to a significant major proportion of nonsynonymous variants with higher pathogenicity scores in the younger branches of the tree. We also derive recalibrated mutation rates for age estimates of ancestors of these various species and proposed a correction curve in order to take into account the effects of selection. Understanding this selection is fundamental to evolutionary studies and to the identification of deleterious mutations.

Genetic characterization, at the mitochondrial and nuclear DNA levels, of five Canary Island dog breeds

Many studies presenting genetic analysis of dog breeds have been conducted without the inclusion of island dog breeds, although isolation can be one of the main factors in their origin. Here we report the genetic analysis at the nuclear and mitochondrial DNA levels of five Canary Island dog breeds (Canarian Warren Hound, Canary Island Mastiff, Garafiano Shepherd, La Palma Rat-Hunter and El Hierro Wolfhound) to fill this gap and, at the same time, genetically characterize these breeds. We identified 168 alleles in autosomal microsatellites and 16 mitochondrial haplotypes. Observed and expected heterozygosities ranged from 0.556 to 0.783 and from 0.737 to 0.943 respectively. Furthermore, three haplotypes were newly described and exclusive to a particular breed (A17+ in the Canary Island Mastiff; A33+ in the Canarian Warren Hound; Bi in the La Palma Rat-Hunter). The outcome of our analyses also revealed different breed histories consistent with historical documents and hypothetical origin designations. Although mtDNA haplotypes showed poor breed discriminating power, autosomal markers allowed a clear clustering of each single population. We expect that our results, together with further analyses, will help to make the population histories of island dog breeds clearer.

Identification of single nucleotide polymorphisms within the mtDNA genome of the domestic dog to discriminate individuals with common HVI haplotypes

We sequenced the entire ∼16 kb canine mitochondrial genome (mtGenome) of 100 unrelated domestic dogs (Canis lupus familiaris) and compared these to 246 published sequences to assess hypervariable region I (HVI) haplotype frequencies. We then used all available sequences to identify informative single nucleotide polymorphisms (SNPs) outside of the control region for use in further resolving mtDNA haplotypes corresponding to common HVI haplotypes. Haplotype frequencies in our data set were highly correlated with previous ones (e.g., F(ST)=0.02, r=0.90), suggesting the total data set reasonably reflected the broader dog population. A total of 128 HVI haplotypes was represented. The 10 most common HVI haplotypes (n=184 dogs) represented 53.3% of the sample. We identified a total 71 SNPs in the mtGenomes (external to the control region) that resolved the 10 most common HVI haplotypes into 63 mtGenome subhaplotypes. The random match probability of the dataset based solely on the HVI sequence was 4%, whereas the random match probability of the mtGenome subhaplotypes was <1%. Thus, the panel of 71 SNPs identified in this study represents a useful forensic tool to further resolve the identity of individual dogs from mitochondrial DNA (mtDNA).

Sequence analysis of the canine mitochondrial DNA control region from shed hair samples in criminal investigations

In recent years, evidence from domestic dogs has increasingly been analyzed by forensic DNA testing. Especially, canine hairs have proved most suitable and practical due to the high rate of hair transfer occurring between dogs and humans. Starting with the description of a contamination-free sample handling procedure, we give a detailed workflow for sequencing hypervariable segments (HVS) of the mtDNA control region from canine evidence. After the hair material is lysed and the DNA extracted by Phenol/Chloroform, the amplification and sequencing strategy comprises the HVS I and II of the canine control region and is optimized for DNA of medium-to-low quality and quantity. The sequencing procedure is based on the Sanger Big-dye deoxy-terminator method and the separation of the sequencing reaction products is performed on a conventional multicolor fluorescence detection capillary electrophoresis platform. Finally, software-aided base calling and sequence interpretation are addressed exemplarily.

Belgian canine population and purebred study for forensics by improved mitochondrial DNA sequencing

In canine population studies for forensics, the mitochondrial DNA is profiled by sequencing the two hyper variable regions, HV1 and HV2 of the control region. In a first effort to create a Belgian population database some samples showed partially poor sequence quality. We demonstrated that a nuclear pseudogene was co-amplified with the mtDNA control region. Using a new combination of primers this adverse result was no longer observed and sequencing quality was improved. All former samples with poor sequence data were reanalyzed. Furthermore, the forensic canine population study was extended to 208 breed and mixed dogs. In total, 58 haplotypes were identified, resulting in an exclusion capacity of 0.92. The profile distribution of the Belgian population sample was not significantly different from those observed in population studies of three other countries. In addition to the total population study 107 Belgian registered pedigree dogs of six breeds were profiled. Per breed, the obtained haplotypes were supplemented with those from population and purebred studies. The combined data revealed that some haplotypes were more or less prominent present in particular dog breeds. The statistically significant differences in haplotype distribution between breeds and population sample can have consequences on mtDNA databasing and matching probabilities in forensics.

ISFG: recommendations regarding the use of non-human (animal) DNA in forensic genetic investigations

The use of non-human DNA typing in forensic science investigations, and specifically that from animal DNA, is ever increasing. The term animal DNA in this document refers to animal species encountered in a forensic science examination but does not include human DNA. Non-human DNA may either be: the trade and possession of a species, or products derived from a species, which is contrary to legislation; as evidence where the crime is against a person or property; instances of animal cruelty; or where the animal is the offender. The first instance is addressed by determining the species present, and the other scenarios can often be addressed by assigning a DNA sample to a particular individual organism. Currently there is little standardization of methodologies used in the forensic analysis of animal DNA or in reporting styles. The recommendations in this document relate specifically to animal DNA that is integral to a forensic science investigation and are not relevant to the breeding of animals for commercial purposes. This DNA commission was formed out of discussions at the International Society for Forensic Genetics 23rd Congress in Buenos Aires to outline recommendations on the use of non-human DNA in a forensic science investigation. Due to the scope of non-human DNA typing that is possible, the remit of this commission is confined to animal DNA typing only.

Segregation of point mutation heteroplasmy in the control region of dog mtDNA studied systematically in deep generation pedigrees

Heteroplasmy, the presence of two or more variants in an organism, may render mitochondrial DNA (mtDNA)-based individual identification challenging in forensic analysis. However, the variation of heteroplasmic proportions and the segregation of heteroplasmic variants through generations and within families have not been systematically described at a large scale in animals such as the domestic dog. Therefore, we performed the largest study to date in domestic dogs and screened a 582-bp-long fragment of the mtDNA control region in 180 individuals in 58 pedigrees for signs of heteroplasmy. We identified three pedigrees (5.17%) with heteroplasmic point mutations. To follow the segregation of the point mutations, we then analyzed 131 samples from these three independent pedigrees and found significant differences in heteroplasmy between generations and among siblings. Frequently (10% of cases), the proportion of one base changed from 0–10% to 80–90% (as judged from Sanger electropherograms) between generations and varied to a similar extent among siblings. We included also a literature review of heteroplasmic and potential mutational hot spot positions in the studied region which showed that all heteroplasmic positions appear to be mutational hot spots. Thus, although heteroplasmy may be used to increase the significance of a match in forensic case work, it may also cause erroneous exclusion of related individuals because of sharp switches from one state to the other within a single generation or among siblings especially in the presented mutational hot spots.

Identification of forensically informative SNPs in the domestic dog mitochondrial control region

Dog hair is often found at crime scenes either due to the dog's involvement in the crime or secondary transfer. As little nuclear DNA is present in shed hair, a 1000 base pair fragment of the mitochondrial control region (mtCR) from 552 dogs was assessed for forensically useful sequence variation. Through pairwise alignment to a standard reference sequence, existing haplotypes were further described and 36 new haplotypes and 24 new single nucleotide polymorphisms were identified. The probability of exclusion was found to be 0.957. Breeds were found to have similar sequences, although not identical. No genetic basis was found for grouping dogs by either purebred or mixed or geographic location within the continental United States. Our research demonstrates that the domestic dog mtCR has not been thoroughly surveyed for sequence variation and that a single database comprised of purebred and mixed breed dogs is sufficient for the continental United States.

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.

Forensic analysis of dog (Canis lupus familiaris) mitochondrial DNA sequences: an inter-laboratory study of the GEP-ISFG working group

A voluntary collaborative exercise aiming at the mitochondrial analysis of canine biological samples was carried out in 2006-2008 by the Non-Human Forensic Genetics Commission of the Spanish and Portuguese Working Group (GEP) of the International Society for Forensic Genetics (ISFG). The participating laboratories were asked to sequence two dog samples (one bloodstain and one hair sample) for the mitochondrial D-loop region comprised between positions 15,372 and 16,083 using suggested primers and PCR conditions, and to compare their results against a reference sequence. Twenty-one participating laboratories reported a total of 67.5% concordant results, 15% non-concordant results, and 17.5% no results. The hair sample analysis presented more difficulty to the participants than the bloodstain analysis, with a high percentage (29%) failing to obtain a result. The high level of participation showed the interest of the community in the analysis of dog forensic samples but the results reveal that crucial methodological issues need to be addressed and further training is required in order to respond proficiently to the demands of forensic casework.

Forensic species identification based on size variation of mitochondrial DNA hypervariable regions

In this study, two new systems for species identification were developed based on size variation of mitochondrial DNA hypervariable regions among animals: one was a conventional method using non-fluorescent primer sets and agarose gel electrophoresis and the other was an automatic method using fluorescent primer sets and capillary electrophoresis. DNA samples from 18 mammal, four birds, and 19 fish species were amplified using three primer sets specific for mammals, birds, and fishes, respectively. The differences in the sizes of the polymerase chain reaction (PCR) products, ranging from about 350 to 900 bp, permitted us to identify species. These systems were successfully applied to various specimens from several criminal cases. In unknown samples, which were different in size from reference DNA markers, sequencing of the PCR products and subsequent BLAST analysis helped to identify species. Furthermore, the sequence data provided us with information on individuals. Because these species identification methods are very simple, easy, rapid, and exact, they are useful in the field of forensic science.

An assessment of the clonality of the components of canine mixed mammary tumours by mitochondrial DNA analysis

The aim of this study was to investigate if mutations in the mitochondrial DNA (mtDNA) D-loop fragment control region of canine mammary mixed tumours could be used as clonal markers that identified the cell population of origin. Ten benign mixed mammary tumours and nine carcinomas arising from benign mixed tumours were microdissected and DNA from epithelial and mesenchymal tumour cells and from normal mammary tissue was examined for sequence variations in a fragment of the hypervariable control region. Identical sequence variants in both the epithelial and mesenchymal components (as well as in the corresponding normal tissue) were found in 80% of the benign mixed tumours and in 89% of the carcinomas arising from benign mixed tumours suggesting a shared clonal origin. The distinctive sequence alterations identified in the epithelial and mesenchymal components of 15.8% of all 19 tumours examined, suggests the possibility that a minority of mammary tumours are polyclonal in origin or that early clonal divergence occurs. Increased mutation within the mtDNA D-loop fragment of mixed tumour components was not observed.

Molecular characterization of the canine mitochondrial DNA control region for forensic applications

The canine mitochondrial DNA (mtDNA) control region of 133 dogs living in the area around Innsbruck, Austria was sequenced. A total of 40 polymorphic sites were observed in the first hypervariable segment and 15 in the second, which resulted in the differentiation of 40 distinct haplotypes. We observed five nucleotide positions that were highly polymorphic within different haplogroups, and they represent good candidates for mtDNA screening. We found five point heteroplasmic positions; all located in HVS-I and a polythymine region in HVS-II, the latter often being associated with length heteroplasmy. In contrast to human mtDNA, the canine control region contains a hypervariable 10 nucleotide repeat region, which is located between the two hypervariable regions. In our population sample, we observed eight different repeat types, which we characterized by direct sequencing and fragment length analysis. The discrimination power of the canine mtDNA control region was 0.93, not taking the polymorphic repeat region into consideration.

Mitochondrial DNA analysis of the domestic dog: control region variation within and among breeds

The mitochondrial DNA (mtDNA) control regions of 125 domestic dogs (Canis familiaris) encompassing 43 breeds, as well as one coyote and two wolves were sequenced and subsequently examined for sequence variation in an effort to construct a reference dog mtDNA data set for forensic analysis. Forty informative variable sites were identified that described 45 haplotypes, 29 of which were observed only once. Substantial variation was found both within and among breeds in the mtDNA derived from tissue, indicating that analysis of the mtDNA derived from dog hairs could be a valuable, discriminating piece of evidence in forensic investigations. The dog data set single nucleotide polymorphisms (SNPs) ranged from having one to six changes on a phylogenetic tree. On average, there were 1.9 character changes for each variable position on the tree. The most variable sites (with four or more changes each, listed from the most changes to the fewest) observed were 15,639 (L=6), 16,672 (L=5), 15,955 (L=4), 15,627 (L=3), 16,431 (L=3), and 16,439 (L=3). These sites were consistent with other reports on variable positions in the dog mtDNA genome. A total of 26 SNPs were chosen to best identify all major clusters in the domestic dog data set. The descriptive analyses revealed that this data set is similar to other published canine data sets and further demonstrates that this domestic dog data set is a useful resource for forensic applications. This reference data set has been compiled and validated against the published dog genetic literature with an aim to aid forensic investigations that seek to incorporate mtDNA sequences and SNPs from trace evidence such as dog hair.

Forensic informativity of domestic dog mtDNA control region sequences

We have analysed the genetic information to be obtained from analysis of mitochondrial DNA (mtDNA) in domestic dogs studying the exclusion capacity in different populations and the correlation between mtDNA types and breeds or types of dogs. The exclusion capacities for a 573 bp sequence of the mitochondrial control region was between 0.86 and 0.95 for dogs in Sweden, the UK, Germany, Japan and China. The direct correlation between mtDNA type and breed, type of dog, and geographical origin of breed was generally low, but in some cases certain mtDNA types were overrepresented in one breed, and for wider groupings such as morphologically similar breeds, some mtDNA types were in many cases found in a distinct group of breeds, often originating from the same geographic region. This type of information may be used as an indication of the breed and, with some degree of probability, to include or exclude certain breeds from being the source of evidence materials.

The mtDNA catalogue of all Portuguese autochthonous goat (Capra hircus) breeds: high diversity of female lineages at the western fringe of European distribution

Goat is believed to be the first true livestock domesticated and, apart from its historical importance, keeps playing an essential economic role in very diverse human societies. We have analysed the female gene pool of all Portuguese autochthonous breeds (Bravia, Serrana, Charnequeira, Serpentina and Algarvia) through the mtDNA HVI sequencing of 288 unrelated animals sampled throughout the country. All breeds proved to be extremely diverse (average haplotype diversity of 0.977), in contrast with the Portuguese peripheral geographic situation in the distribution range of the species. Moreover, observed genetic distances between breeds do not correlate with microgeography inside Portugal. These observations are consistent with recurrent refreshment of the breeding stock through the introduction of exotic animals. Fitting the new data into the still loosely defined female genetic pool landscape of goats, all Portuguese animals, one sample excepted (belonging to Bravia), are classified into haplogroup A and haplotype sharing is geographically very sparse, including a Far East match. Our results confirm that goats stand out among most of domesticates as exceptionally diverse and showing an unparalleled degree of mobility of animals (at least females) used for reproductive purposes.

MtDNA diversity among four Portuguese autochthonous dog breeds: a fine-scale characterisation

Background

The picture of dog mtDNA diversity, as obtained from geographically wide samplings but from a small number of individuals per region or breed, has revealed weak geographic correlation and high degree of haplotype sharing between very distant breeds. We aimed at a more detailed picture through extensive sampling (n = 143) of four Portuguese autochthonous breeds – Castro Laboreiro Dog, Serra da Estrela Mountain Dog, Portuguese Sheepdog and Azores Cattle Dog-and comparatively reanalysing published worldwide data.

Results

Fifteen haplotypes belonging to four major haplogroups were found in these breeds, of which five are newly reported. The Castro Laboreiro Dog presented a 95% frequency of a new A haplotype, while all other breeds contained a diverse pool of existing lineages. The Serra da Estrela Mountain Dog, the most heterogeneous of the four Portuguese breeds, shared haplotypes with the other mainland breeds, while Azores Cattle Dog shared no haplotypes with the other Portuguese breeds.

A review of mtDNA haplotypes in dogs across the world revealed that: (a) breeds tend to display haplotypes belonging to different haplogroups; (b) haplogroup A is present in all breeds, and even uncommon haplogroups are highly dispersed among breeds and continental areas; (c) haplotype sharing between breeds of the same region is lower than between breeds of different regions and (d) genetic distances between breeds do not correlate with geography.

Conclusion

MtDNA haplotype sharing occurred between Serra da Estrela Mountain dogs (with putative origin in the centre of Portugal) and two breeds in the north and south of the country-with the Castro Laboreiro Dog (which behaves, at the mtDNA level, as a sub-sample of the Serra da Estrela Mountain Dog) and the southern Portuguese Sheepdog. In contrast, the Azores Cattle Dog did not share any haplotypes with the other Portuguese breeds, but with dogs sampled in Northern Europe. This suggested that the Azores Cattle Dog descended maternally from Northern European dogs rather than Portuguese mainland dogs. A review of published mtDNA haplotypes identified thirteen non-Portuguese breeds with sufficient data for comparison. Comparisons between these thirteen breeds, and the four Portuguese breeds, demonstrated widespread haplotype sharing, with the greatest diversity among Asian dogs, in accordance with the central role of Asia in canine domestication.