Evolutionary analysis of a large mtDNA translocation (numt) into the nuclear genome of the Panthera genus species

Mitogenomic analysis of a late Pleistocene jaguar from North America

The jaguar (Panthera onca) is the largest living cat species native to the Americas and one of few large American carnivorans to have survived into the Holocene. However, the extent to which jaguar diversity declined during the end-Pleistocene extinction event remains unclear. For example, Pleistocene jaguar fossils from North America are notably larger than the average extant jaguar, leading to hypotheses that jaguars from this continent represent a now-extinct subspecies (Panthera onca augusta) or species (Panthera augusta). Here, we used a hybridization capture approach to recover an ancient mitochondrial genome from a large, late Pleistocene jaguar from Kingston Saltpeter Cave, Georgia, United States, which we sequenced to 26-fold coverage. We then estimated the evolutionary relationship between the ancient jaguar mitogenome and those from other extinct and living large felids, including multiple jaguars sampled across the species' current range. The ancient mitogenome falls within the diversity of living jaguars. All sampled jaguar mitogenomes share a common mitochondrial ancestor ~400 thousand years ago, indicating that the lineage represented by the ancient specimen dispersed into North America from the south at least once during the late Pleistocene. While genomic data from additional and older specimens will continue to improve understanding of Pleistocene jaguar diversity in the Americas, our results suggest that this specimen falls within the variation of extant jaguars despite the relatively larger size and geographic location and does not represent a distinct taxon.

From Caves to the Savannah, the Mitogenome History of Modern Lions (Panthera leo) and Their Ancestors

Lions (Panthera leo) play a crucial ecological role in shaping and maintaining fragile ecosystems within Africa. Conservation efforts should focus on genetic variability within wild populations when considering reintroduction attempts. We studied two groups of lions from two conservation sites located in Zambia and Zimbabwe to determine their genetic make-up, information that is usually unknown to the sites. In this study, we analysed 17 specimens for cytb and seven microsatellite markers to ascertain family relationships and genetic diversity previously obtained by observational studies. We then produced a standardised haplogroup phylogeny using all available entire mitogenomes, as well as calculating a revised molecular clock. The modern lion lineage diverged ~151 kya and was divided into two subspecies, both containing three distinct haplogroups. We confirm that Panthera leo persica is not a subspecies, but rather a haplogroup of the northern P.l. leo that exited Africa at least ~31 kya. The progenitor to all lions existed ~1.2 Mya, possibly in SE Africa, and later exited Africa and split into the two cave lion lineages ~175 kya. Species demography is correlated to major climactic events. We now have a detailed phylogeny of lion evolution and an idea of their conservation status given the threat of climate change.

Defining cat mitogenome variation and accounting for numts via multiplex amplification and Nanopore sequencing

Hair shed by domestic cats is a potentially useful source of forensic evidence. Analysable hair DNA is predominantly mitochondrial, but the recent domestication history of cats means that mtDNA diversity is low. A 402-bp control region segment is usually sequenced, defining only a small number of distinct haplotypes in populations. Previously, we used a long-amplicon approach to sequence whole mitogenomes in a sample of blood DNAs from 119 UK cats, greatly increasing observed diversity and reducing random match probabilities. To exploit this variation for forensic analysis, we here describe a multiplex system that amplifies the cat mitogenome in 60 overlapping amplicons of mean length 360 bp, followed by Nanopore sequencing. Variants detected in multiplex sequence data from unrooted hair completely mirror those from long-amplicon data from blood from the same individuals. However, applying the multiplex to matched blood DNA reveals additional sequence variants which derive from the major feline nuclear mitochondrial insertion sequence (numt), which covers 7.9 kb of the 17-kb mitogenome and exists in multiple tandem copies. We use long-amplicon Nanopore sequencing to investigate numt variation in a set of cats, together with an analysis of published genome sequences, and show that numt arrays are variable in both structure and sequence, thus providing a potential source of uncertainty when nuclear DNA predominates in a sample. Forensic application of the multiplex was demonstrated by matching hairs from a cat with skeletal remains from its putative mother, both of which shared a globally common haplotype at the control region. The random match probability in this case with the CR 402-bp segment was 0.21 and this decreased to 0.03 when considering the whole mitogenome. The developed multiplex and sequencing approach, when applied to cat hair where nuclear DNA is scarce, can provide a reliable and highly discriminating source of forensic genetic evidence from a single hair. The confounding effect of numt co-amplification in degraded samples where mixed sequences are observed can be mitigated by variant phasing, and by comparison with numt sequence diversity data, such as those presented here.

Ancient DNA reveals genetic admixture in China during tiger evolution

The tiger (Panthera tigris) is a charismatic megafauna species that originated and diversified in Asia and probably experienced population contraction and expansion during the Pleistocene, resulting in low genetic diversity of modern tigers. However, little is known about patterns of genomic diversity in ancient populations. Here we generated whole-genome sequences from ancient or historical (100-10,000 yr old) specimens collected across mainland Asia, including a 10,600-yr-old Russian Far East specimen (RUSA21, 8× coverage) plus six ancient mitogenomes, 14 South China tigers (0.1-12×) and three Caspian tigers (4-8×). Admixture analysis showed that RUSA21 clustered within modern Northeast Asian phylogroups and partially derived from an extinct Late Pleistocene lineage. While some of the 8,000-10,000-yr-old Russian Far East mitogenomes are basal to all tigers, one 2,000-yr-old specimen resembles present Amur tigers. Phylogenomic analyses suggested that the Caspian tiger probably dispersed from an ancestral Northeast Asian population and experienced gene flow from southern Bengal tigers. Lastly, genome-wide monophyly supported the South China tiger as a distinct subspecies, albeit with mitochondrial paraphyly, hence resolving its longstanding taxonomic controversy. The distribution of mitochondrial haplogroups corroborated by biogeographical modelling suggested that Southwest China was a Late Pleistocene refugium for a relic basal lineage. As suitable habitat returned, admixture between divergent lineages of South China tigers took place in Eastern China, promoting the evolution of other northern subspecies. Altogether, our analysis of ancient genomes sheds light on the evolutionary history of tigers and supports the existence of nine modern subspecies.

Development of short-target primers for species identification in biological studies of Carnivora

Noninvasive genetic sampling greatly facilitates studies on the genetics, ecology, and conservation of threatened species. Species identification is often a prerequisite for noninvasive sampling-based biological studies. Due to the low quantity and quality of genomic DNA from noninvasive samples, high-performance short-target PCR primers are necessary for DNA barcoding applications. The order Carnivora is characterized by an elusive habit and threatened status. In this study, we developed three pairs of short-target primers for identifying Carnivora species. The COI279 primer pair was suitable for samples with better DNA quality. The COI157a and COI157b primer pairs performed well for noninvasive samples and reduced the interference of nuclear mitochondrial pseudogenes (numts). COI157a could effectively identify samples from Felidae, Canidae, Viverridae, and Hyaenidae, while COI157b could be applied to samples from Ursidae, Ailuridae, Mustelidae, Procyonidae, and Herpestidae. These short-target primers will facilitate noninvasive biological studies and efforts to conserve Carnivora species.

Lions and brown bears colonized North America in multiple synchronous waves of dispersal across the Bering Land Bridge

The Bering Land Bridge connecting North America and Eurasia was periodically exposed and inundated by oscillating sea levels during the Pleistocene glacial cycles. This land connection allowed the intermittent dispersal of animals, including humans, between Western Beringia (far northeast Asia) and Eastern Beringia (northwest North America), changing the faunal community composition of both continents. The Pleistocene glacial cycles also had profound impacts on temperature, precipitation and vegetation, impacting faunal community structure and demography. While these palaeoenvironmental impacts have been studied in many large herbivores from Beringia (e.g., bison, mammoths, horses), the Pleistocene population dynamics of the diverse guild of carnivorans present in the region are less well understood, due to their lower abundances. In this study, we analyse mitochondrial genome data from ancient brown bears (Ursus arctos; n = 103) and lions (Panthera spp.; n = 39), two megafaunal carnivorans that dispersed into North America during the Pleistocene. Our results reveal striking synchronicity in the population dynamics of Beringian lions and brown bears, with multiple waves of dispersal across the Bering Land Bridge coinciding with glacial periods of low sea levels, as well as synchronous local extinctions in Eastern Beringia during Marine Isotope Stage 3. The evolutionary histories of these two taxa underline the crucial biogeographical role of the Bering Land Bridge in the distribution, turnover and maintenance of megafaunal populations in North America.

Comparison of detection methods and genome quality when quantifying nuclear mitochondrial insertions in vertebrate genomes

The integration of mitochondrial genome fragments into the nuclear genome is well documented, and the transfer of these mitochondrial nuclear pseudogenes (numts) is thought to be an ongoing evolutionary process. With the increasing number of eukaryotic genomes available, genome-wide distributions of numts are often surveyed. However, inconsistencies in genome quality can reduce the accuracy of numt estimates, and methods used for identification can be complicated by the diverse sizes and ages of numts. Numts have been previously characterized in rodent genomes and it was postulated that they might be more prevalent in a group of voles with rapidly evolving karyotypes. Here, we examine 37 rodent genomes, and an additional 26 vertebrate genomes, while also considering numt detection methods. We identify numts using DNA:DNA and protein:translated-DNA similarity searches and compare numt distributions among rodent and vertebrate taxa to assess whether some groups are more susceptible to transfer. A combination of protein sequence comparisons (protein:translated-DNA) and BLASTN genomic DNA searches detect 50% more numts than genomic DNA:DNA searches alone. In addition, higher-quality RefSeq genomes produce lower estimates of numts than GenBank genomes, suggesting that lower quality genome assemblies can overestimate numts abundance. Phylogenetic analysis shows that mitochondrial transfers are not associated with karyotypic diversity among rodents. Surprisingly, we did not find a strong correlation between numt counts and genome size. Estimates using DNA: DNA analyses can underestimate the amount of mitochondrial DNA that is transferred to the nucleus.

Rapid Macrosatellite Evolution Promotes X-Linked Hybrid Male Sterility in a Feline Interspecies Cross

The sterility or inviability of hybrid offspring produced from an interspecific mating result from incompatibilities between parental genotypes that are thought to result from divergence of loci involved in epistatic interactions. However, attributes contributing to the rapid evolution of these regions also complicates their assembly, thus discovery of candidate hybrid sterility loci is difficult and has been restricted to a small number of model systems. Here we reported rapid interspecific divergence at the DXZ4 macrosatellite locus in an interspecific cross between two closely related mammalian species: the domestic cat (Felis silvestris catus) and the Jungle cat (Felis chaus). DXZ4 is an interesting candidate due to its structural complexity, copy number variability, and described role in the critical yet complex biological process of X-chromosome inactivation. However, the full structure of DXZ4 was absent or incomplete in nearly every available mammalian genome assembly given its repetitive complexity. We compared highly continuous genomes for three cat species, each containing a complete DXZ4 locus, and discovered that the felid DXZ4 locus differs substantially from the human ortholog, and that it varies in copy number between cat species. Additionally, we reported expression, methylation, and structural conformation profiles of DXZ4 and the X chromosome during stages of spermatogenesis that have been previously associated with hybrid male sterility. Collectively, these findings suggest a new role for DXZ4 in male meiosis and a mechanism for feline interspecific incompatibility through rapid satellite divergence.

Early Pleistocene origin and extensive intra-species diversity of the extinct cave lion

The cave lion is an extinct felid that was widespread across the Holarctic throughout the Late Pleistocene. Its closest extant relative is the lion (Panthera leo), but the timing of the divergence between these two taxa, as well as their taxonomic ranking are contentious. In this study we analyse 31 mitochondrial genome sequences from cave lion individuals that, through a combination of ¹⁴C and genetic tip dating, are estimated to be from dates extending well into the mid-Pleistocene. We identified two deeply diverged and well-supported reciprocally monophyletic mitogenome clades in the cave lion, and an additional third distinct lineage represented by a single individual. One of these clades was restricted to Beringia while the other was prevalent across western Eurasia. These observed clade distributions are in line with previous observations that Beringian and European cave lions were morphologically distinct. The divergence dates for these lineages are estimated to be far older than those between extant lions subspecies. By combining our radiocarbon tip-dates with a split time prior that takes into account the most up-to-date fossil stem calibrations, we estimated the mitochondrial DNA divergence between cave lions and lions to be 1.85 Million ya (95% 0.52- 2.91 Mya). Taken together, these results support previous hypotheses that cave lions existed as at least two subspecies during the Pleistocene, and that lions and cave lions were distinct species.

Spatial and temporal genetic diversity of the Texas kangaroo rat, Dipodomys elator (Rodentia: Heteromyidae)

The Texas kangaroo rat (Dipodomys elator) is listed as a threatened species in Texas because of its scarcity and small geographic range. We assessed patterns of genetic diversity in D. elator that could affect extinction risk or influence management decisions. Specific objectives included: 1) document levels of genetic diversity, 2) document the degree and patterns of genetic divergence among localities, and 3) compare levels of genetic diversity between different time periods at the same locality. Portions of the mitochondrial genome (mtDNA; control region, cytochrome c oxidase subunit I, and cytochrome b) were sequenced and nuclear microsatellites were examined. Low mtDNA diversity was observed, which could be explained by an historical, species-wide genetic bottleneck. In contrast, microsatellites exhibited ample variation, and analyses were conducted using data from 11 loci and four populations (designated Quanah, Iowa Park, Vernon, and Harrold). Allelic diversity and heterozygosity were similar between populations and temporal samples. Estimates of effective population size (Ne) ranged from 5 to 856, depending on method and population, with Iowa Park showing consistently lower values than Quanah. All methods addressing population structure indicated that the Iowa Park population was divergent from the others, with Vernon and Harrold showing a somewhat intermediate relationship but with a closer affiliation with Quanah than Iowa Park, despite their closer proximity to Iowa Park. This pattern did not conform to isolation by distance, thus genetic drift appears to have played a greater role than gene flow in establishing genetic structure. There was much less difference between temporal samples compared to geographic samples, indicating that genetic drift has had only minimal impacts in shifting allelic frequencies over the time periods examined (17–36 years).

Genetic analysis of African lions (Panthera leo) in Zambia support movement across anthropogenic and geographical barriers

The Luangwa Valley in eastern Zambia is a transverse offshoot of the Great Rift Valley system. This region appears to have an isolating effect as evidenced by suspected endemic subspecies, such as the Cookson's wildebeest and Thornicroft's giraffe. Recent mitochondrial DNA studies demonstrated that African lions in Zambia consist of two highly diverse eastern and western sub-populations. Herein, we report nuclear and mitochondrial DNA results from 409 lions that support this population substructure across Zambia but proposes only partial isolation of the Luangwa Valley with more movement between the populations than previously thought. Population assignment analysis identifies two populations with little evidence of admixture assigning lions to either the eastern or western sub-populations. A high occurrence of private alleles and clear evidence for a Wahlund effect further justify the presence of a highly structured population. But, while mitochondrial DNA analysis still shows little to no matrilineal gene flow (FST = 0.53) between sub-populations, microsatellite analysis suggests there is gene flow (FST = 0.04) with low but significant isolation-by-distance and an average of 6 migrants per generation. Evidence of isolation-by-distance is also found in factorial correspondence analysis with the Lower Zambezi National Park and eastern corridor clusters overlapping isolated clusters of the Luangwa Valley and western sub-population. From this evidence, the Luangwa Valley appears separated from the western sub-population with some dispersal through the southern regions of the eastern sub-population. Both the eastern and western sub-populations have high heterozygosity (0.68 and 0.69, respectively) and genetic diversity (0.47 and 0.50, respectively) values, indicative of genetically healthy populations.

Possibility of numt co-amplification from gigantic genome of Orthoptera: testing efficiency of standard PCR protocol in producing orthologous COI sequences

Mitochondrial DNA has been the preferential genome biodiversity studies. However, several factors contribute to its inadequacy. Numts constitute one of the main complications that prevent obtaining orthologous mitochondrial sequences. Orthoptera have been a model group in numt studies because of their huge genome size. In this study we aimed to; (i) test efficiency of standard PCR protocol in producing orthologous sequences of cytochrome C oxidase, (ii) study presence/absence of numts in several unstudied Orthoptera species, (iii) test if there is a threshold between the length of mtDNA targeted for amplification and possibility of encountering numts, and (iv) estimate reliability of the sequences in databases in light of these findings. For these aims we studied 38 species of Orthoptera representing different sublineages and genome sizes. DNA extracted from each sample was used to amplify five different fragments of COI region by standard PCR protocol. Sequenced PCR amplicons were checked for numt possibility by several different numt criteria. No sequences without numt signs were obtained for the first fragment. The number of samples with numt signs for the other four fragments differed between the suborders Ensifera and Caelifera. The percentage of samples with numt signs was higher in Caelifera than Ensifera for all fragments. The numt percentage considerably decreased for the longest two fragments. Numts are more prevalent in families with larger genome size. We arrived at the following conclusions: (i) numts are common in all members of Orthoptera, but, their prevalence differs among intra-lineages, especially more prevalent in Caelifera, (ii) there seems a correlation between numt rate and genome size, (iii) there is no threshold to avoid numt co-amplification, but, a 1,000 bp length may be a threshold for Ensifera, (iv) Folmer region of COI doesn't seem an appropriate marker for animal barcoding. Additionally, a phylogenetic tree produced from the numt sequences of fragment four detected in genus Anterastes suggested a paleonumt gained in generic ancestor a 3.5-4 times slower divergence rate for numt sequences.

First Report of a Mitochondrial Pseudogene in Agnathan Vertebrates (Cyclostomata: Petromyzontidae)

We report herein the characterization of a nuclear paralog of a fragment of the mitochondrial genome (a numt) in two closely related species of lampreys (Ichthyomyzon spp.). Although numts have been characterized in several vertebrate taxa, numts have yet to be reported for fishes in general. Given the phylogenetic position of lampreys relative to other vertebrates, the presence of numts within the lamprey genome is either evidence of an ancestral trait lost in other fishes but uniquely retained in agnathans and amniotes, or (more intriguingly) a product of the genome rearrangements these animals undergo during development.

Occurrence of mitochondrial CO1 pseudogenes in Neocalanus plumchrus (Crustacea: Copepoda): Hybridization indicated by recombined nuclear mitochondrial pseudogenes

A portion of the mitochondrial cytochrome c oxidase I gene was sequenced using both genomic DNA and complement DNA from three planktonic copepod Neocalanus species (N. cristatus, N. plumchrus, and N. flemingeri). Small but critical sequence differences in CO1 were observed between gDNA and cDNA from N. plumchrus. Furthermore, careful observation revealed the presence of recombination between sequences in gDNA from N. plumchrus. Moreover, a chimera of the N. cristatus and N. plumchrus sequences was obtained from N. plumchrus gDNA. The observed phenomena can be best explained by the preferential amplification of the nuclear mitochondrial pseudogenes from gDNA of N. plumchrus. Two conclusions can be drawn from the observations. First, nuclear mitochondrial pseudogenes are pervasive in N. plumchrus. Second, a mating between a female N. cristatus and a male N. plumchrus produced viable offspring, which further backcrossed to a N. plumchrus individual. These observations not only demonstrate intriguing mating behavior in these species, but also emphasize the importance of careful interpretation of species marker sequences amplified from gDNA.

Mitochondrial Haplotype Diversity in Zambian Lions: Bridging a Gap in the Biogeography of an Iconic Species

Analysis of DNA sequence diversity at the 12S to 16S mitochondrial genes of 165 African lions (Panthera leo) from five main areas in Zambia has uncovered haplotypes which link Southern Africa with East Africa. Phylogenetic analysis suggests Zambia may serve as a bridge connecting the lion populations in southern Africa to eastern Africa, supporting earlier hypotheses that eastern-southern Africa may represent the evolutionary cradle for the species. Overall gene diversity throughout the Zambian lion population was 0.7319 +/- 0.0174 with eight haplotypes found; three haplotypes previously described and the remaining five novel. The addition of these five novel haplotypes, so far only found within Zambia, nearly doubles the number of haplotypes previously reported for any given geographic location of wild lions. However, based on an AMOVA analysis of these haplotypes, there is little to no matrilineal gene flow (Fst = 0.47) when the eastern and western regions of Zambia are considered as two regional sub-populations. Crossover haplotypes (H9, H11, and Z1) appear in both populations as rare in one but common in the other. This pattern is a possible result of the lion mating system in which predominately males disperse, as all individuals with crossover haplotypes were male. The determination and characterization of lion sub-populations, such as done in this study for Zambia, represent a higher-resolution of knowledge regarding both the genetic health and connectivity of lion populations, which can serve to inform conservation and management of this iconic species.

Looking back on a decade of barcoding crustaceans

Species identification represents a pivotal component for large-scale biodiversity studies and conservation planning but represents a challenge for many taxa when using morphological traits only. Consequently, alternative identification methods based on molecular markers have been proposed. In this context, DNA barcoding has become a popular and accepted method for the identification of unknown animals across all life stages by comparison to a reference library. In this review we examine the progress of barcoding studies for the Crustacea using the Web of Science data base from 2003 to 2014. All references were classified in terms of taxonomy covered, subject area (identification/library, genetic variability, species descriptions, phylogenetics, methods, pseudogenes/numts), habitat, geographical area, authors, journals, citations, and the use of the Barcode of Life Data Systems (BOLD). Our analysis revealed a total number of 164 barcoding studies for crustaceans with a preference for malacostracan crustaceans, in particular Decapoda, and for building reference libraries in order to identify organisms. So far, BOLD did not establish itself as a popular informatics platform among carcinologists although it offers many advantages for standardized data storage, analyses and publication.

Odintifier--A computational method for identifying insertions of organellar origin from modern and ancient high-throughput sequencing data based on haplotype phasing

Background

Cellular organelles with genomes of their own (e.g. plastids and mitochondria) can pass genetic sequences to other organellar genomes within the cell in many species across the eukaryote phylogeny. The extent of the occurrence of these organellar-derived inserted sequences (odins) is still unknown, but if not accounted for in genomic and phylogenetic studies, they can be a source of error. However, if correctly identified, these inserted sequences can be used for evolutionary and comparative genomic studies. Although such insertions can be detected using various laboratory and bioinformatic strategies, there is currently no straightforward way to apply them as a standard organellar genome assembly on next-generation sequencing data. Furthermore, most current methods for identification of such insertions are unsuitable for use on non-model organisms or ancient DNA datasets.

Results

We present a bioinformatic method that uses phasing algorithms to reconstruct both source and inserted organelle sequences. The method was tested in different shotgun and organellar-enriched DNA high-throughput sequencing (HTS) datasets from ancient and modern samples. Specifically, we used datasets from lions (Panthera leo ssp. and Panthera leo leo) to characterize insertions from mitochondrial origin, and from common grapevine (Vitis vinifera) and bugle (Ajuga reptans) to characterize insertions derived from plastid genomes. Comparison of the results against other available organelle genome assembly methods demonstrated that our new method provides an improvement in the sequence assembly.

Conclusion

Using datasets from a wide range of species and different levels of complexity we showed that our novel bioinformatic method based on phasing algorithms can be used to achieve the next two goals: i) reference-guided assembly of chloroplast/mitochondrial genomes from HTS data and ii) identification and simultaneous assembly of odins. This method represents the first application of haplotype phasing for automatic detection of odins and reference-based organellar genome assembly.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0682-1) contains supplementary material, which is available to authorized users.

Phylogenetic relationship of Eurasian lynx (Lynx lynx) revealed by complete mitochondrial genome

The Eurasian lynx (Lynx lynx) is an Endangered species in northeast China. We first obtained muscle sample, extracted the sample DNA and sequenced the whole mtDNA genome of lynx from northeast China. We reconstructed the phylogenetic tree of Eurasian lynx and 10 other most closely related Felidae species. This lynx's complete mitogenome is 17 054bp in length, includes 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and one control region. The phylogenetic tree confirmed previous research results.

Genetic ancestry of the extinct Javan and Bali tigers

The Bali (Panthera tigris balica) and Javan (P. t. sondaica) tigers are recognized as distinct tiger subspecies that went extinct in the 1940s and 1980s, respectively. Yet their genetic ancestry and taxonomic status remain controversial. Following ancient DNA procedures, we generated concatenated 1750bp mtDNA sequences from 23 museum samples including 11 voucher specimens from Java and Bali and compared these to diagnostic mtDNA sequences from 122 specimens of living tiger subspecies and the extinct Caspian tiger. The results revealed a close genetic affinity of the 3 groups from the Sunda Islands (Bali, Javan, and Sumatran tigers P. t. sumatrae). Bali and Javan mtDNA haplotypes differ from Sumatran haplotypes by 1-2 nucleotides, and the 3 island populations define a monophyletic assemblage distinctive and equidistant from other mainland subspecies. Despite this close phylogenetic relationship, no mtDNA haplotype was shared between Sumatran and Javan/Bali tigers, indicating little or no matrilineal gene flow among the islands after they were colonized. The close phylogenetic relationship among Sunda tiger subspecies suggests either recent colonization across the islands, or else a once continuous tiger population that had subsequently isolated into different island subspecies. This supports the hypothesis that the Sumatran tiger is the closest living relative to the extinct Javan and Bali tigers.

Revealing the maternal demographic history of Panthera leo using ancient DNA and a spatially explicit genealogical analysis

Background

Understanding the demographic history of a population is critical to conservation and to our broader understanding of evolutionary processes. For many tropical large mammals, however, this aim is confounded by the absence of fossil material and by the misleading signal obtained from genetic data of recently fragmented and isolated populations. This is particularly true for the lion which as a consequence of millennia of human persecution, has large gaps in its natural distribution and several recently extinct populations.

Results

We sequenced mitochondrial DNA from museum-preserved individuals, including the extinct Barbary lion (Panthera leo leo) and Iranian lion (P. l. persica), as well as lions from West and Central Africa. We added these to a broader sample of lion sequences, resulting in a data set spanning the historical range of lions. Our Bayesian phylogeographical analyses provide evidence for highly supported, reciprocally monophyletic lion clades. Using a molecular clock, we estimated that recent lion lineages began to diverge in the Late Pleistocene. Expanding equatorial rainforest probably separated lions in South and East Africa from other populations. West African lions then expanded into Central Africa during periods of rainforest contraction. Lastly, we found evidence of two separate incursions into Asia from North Africa, first into India and later into the Middle East.

Conclusions

We have identified deep, well-supported splits within the mitochondrial phylogeny of African lions, arguing for recognition of some regional populations as worthy of independent conservation. More morphological and nuclear DNA data are now needed to test these subdivisions.

Hybridisation, paternal leakage and mitochondrial DNA linearization in three anomalous fish (Scombridae)

Using mitochondrial DNA for species identification and population studies assumes that the genome is maternally inherited, circular, located in the cytoplasm and lacks recombination. This study explores the mitochondrial genomes of three anomalous mackerel. Complete mitochondrial genome sequencing plus nuclear microsatellite genotyping of these fish identified them as Scomberomorus munroi (spotted mackerel). Unlike normal S. munroi, these three fish also contained different linear, mitochondrial genomes of Scomberomorus semifasciatus (grey mackerel). The results are best explained by hybridisation, paternal leakage and mitochondrial DNA linearization. This unusual observation may provide an explanation for mtDNA outliers in animal population studies.

Nuclear mitochondrial pseudogenes in Austinograea alayseae hydrothermal vent crabs (Crustacea: Bythograeidae): effects on DNA barcoding

Members of the brachyuran crab family, Bythograeidae, are among the most abundant and common crabs in vent fields. However, their identification based on morphological characteristics often leads to incorrect species recognition due to a lack of taxonomic factors and the existence of sibling (or cryptic) species. For these reasons, we used DNA barcoding for vent crabs using mitochondrial cytochrome c oxidase subunit 1 (CO1). However, several nuclear mitochondrial pseudogenes (Numts) were amplified from Austinograea alayseae Guinot, 1990, using universal primers (Folmer primers). The Numts were characterized in six haplotypes, with 13.58-14.11% sequence divergence from A. alayseae, a higher nonsynonymous substitution ratio than true CO1, and the formation of an independent clade in bythograeids. In a neighbour-joining tree, the origin of the Numts would be expected to incorporate into the nucleus at an ancestral node of Austinograea, and they mutated more slowly in the nucleus than CO1 in the mitochondria. This evolutionary process may have resulted in the higher binding affinity of Numts for the Folmer primers than CO1. In the present study, we performed long PCR for the amplification of CO1 in A. alayseae. We also present evidence that Numts can introduce serious ambiguity into DNA barcoding, including overestimating the number of species in bythograeids. These results may help in conducting taxonomic studies using mitochondrial genes from organisms living in hydrothermal vent fields.

[Summary of phylogeny in family Felidae of Carnivora]

Felidae (cats) is one of the strict carnivorous groups in the order Carnivora, many of which are most familiar and spectacular to us. They are the top predators in the world. Thirty-six of 37 living cat species are considered as either "endangered" or "threatened". The relationships among species of the family Felidae, which evolved recently and rapidly, are difficult to resolve, and have been the subject of debate. Construction of a reliable Felidae phylogeny will be of evolutionarily significance and conservation value. In this paper, we summarized phylogeny of Felidae, including cytological, morphological and molecular evidence, and pointed out the existing phylogenetic problems. This review is expected to guide future researches of Felidae phylogeny, and to lay a theoretic foundation for the protection of this animal group.

The complete mitochondrial genome analysis of the tiger (Panthera tigris)

The complete mitochondrial genomes of five tiger samples from three subspecies (P. t. sumatrae, P. t. altica, and P. t. tigris) were successfully obtained by using 26 specifically designed Panthera-specific primer sets. The genome organization and gene arrangement of the five tiger samples were similar to each other; however polymorphic tandem repeat sequences were observed in the control region (CR). This led to a difference in the genome lengths obtained from these five samples with an average size of 16,994 bp for the five tiger mitochondrial genomes. The nucleotide base composition was on average as follows: A, 31.8%; T, 27.0%; C, 26.6%; G, 14.6% and exhibited compositional asymmetry. Most of tiger mitochondrial genome characteristics are similar to those of other common vertebrate species; however, some distinctive features were observed in the CR. First, the repetitive sequence 2 (RS 2) contained two repeat units of 80 bp and the first 15 bp of what would be the third repeat motif. The repetitive sequence 3 (RS 3) contained 47-50 repeat motifs of a shorter 8 bp (ACGTAYAC)(n). Second, length heteroplasmy polycystosine (poly-C) stretches was observed at the end of the HV I locus in all tiger samples.

Complete sequences of mitochondria genomes of Aedes aegypti and Culex quinquefasciatus and comparative analysis of mitochondrial DNA fragments inserted in the nuclear genomes

We present complete sequences of the mitochondrial genomes for two important mosquitoes, Aedes aegypti and Culex quinquefasciatus, that are major vectors of dengue virus and lymphatic filariasis, respectively. The A. aegypti mitochondrial genome is 16,655 bp in length and that of C. quinquefasciatus is 15,587 bp, yet both contain 13 protein coding genes, 22 transfer RNA (tRNA) genes, one 12S ribosomal RNA (rRNA) gene, one 16S rRNA gene and a control region (CR) in the same order. The difference in the genome size is due to the difference in the length of the control region. We also analyzed insertions of nuclear copies of mtDNA-like sequences (NUMTs) in a comparative manner between the two mosquitoes. The NUMT sequences occupy ~0.008% of the A. aegypti genome and ~0.001% of the C. quinquefasciatus genome. Several NUMTs were found localized in the introns of predicted protein coding genes in both genomes (32 genes in A. aegypti but only four in C. quinquefasciatus). None of these NUMT-containing genes had an ortholog between the two species or had paralogous copies within a genome that was also NUMT-containing. It was further observed that the NUMT-containing genes were relatively longer but had lower GC content compared to the NUMT-less paralogous copies. Moreover, stretches of homologies are present among the genic and non-genic NUMTs that may play important roles in genomic rearrangement of NUMTs in these genomes. Our study provides new insights on understanding the roles of nuclear mtDNA sequences in genome complexities of these mosquitoes.

DNA barcoding meets molecular scatology: short mtDNA sequences for standardized species assignment of carnivore noninvasive samples

Although species assignment of scats is important to study carnivore biology, there is still no standardized assay for the identification of carnivores worldwide, which would allow large-scale routine assessments and reliable cross-comparison of results. Here, we evaluate the potential of two short mtDNA fragments [ATP6 (126 bp) and cytochrome oxidase I gene (COI) (187 bp)] to serve as standard markers for the Carnivora. Samples of 66 species were sequenced for one or both of these segments. Alignments were complemented with archival sequences and analysed with three approaches (tree-based, distance-based and character-based). Intraspecific genetic distances were generally lower than between-species distances, resulting in diagnosable clusters for 86% (ATP6) and 85% (COI) of the species. Notable exceptions were recently diverged species, most of which could still be identified using diagnostic characters and uniqueness of haplotypes or by reducing the geographic scope of the comparison. In silico analyses were also performed for a 110-bp cytochrome b (cytb) segment, whose identification success was lower (70%), possibly due to the smaller number of informative sites and/or the influence of misidentified sequences obtained from GenBank. Finally, we performed case studies with faecal samples, which supported the suitability of our two focal markers for poor-quality DNA and allowed an assessment of prey DNA co-amplification. No evidence of prey DNA contamination was found for ATP6, while some cases were observed for COI and subsequently eliminated by the design of more specific primers. Overall, our results indicate that these segments hold good potential as standard markers for accurate species-level identification in the Carnivora.

Applying molecular genetic tools to tiger conservation

The utility of molecular genetic approaches in conservation of endangered taxa is now commonly recognized. Over the past decade, conservation genetic analyses based on mitochondrial DNA sequencing and microsatellite genotyping have provided powerful tools to resolve taxonomy uncertainty of tiger subspecies, to define conservation units, to reconstruct phylogeography and demographic history, to examine the genetic ancestry of extinct subspecies, to assess population genetic status non-invasively, and to verify genetic background of captive tigers worldwide. The genetic status of tiger subspecies and populations and implications for developing strategies for the survival of this charismatic species both in situ and ex situ are discussed.

Confirmation of germ-line transmission in the red fluorescence protein (RFP) transgenic cloned male cat

The production of transgenic animals is highly desirable for biotechnology and basic research. We investigated the reproductive ability of a red fluorescence protein (RFP) transgenic cloned male cat (RFP TG cat) by natural mating with a domestic female cat. The RFP expression levels were examined in early embryogenesis, tissues from 45-day-old two fetuses, and four RFP TG cat offspring. The RFP gene was detected in tissue samples from one dead kitten, including several organs and the skin. Also, under a fluorescent light source, we were able to directly detect the RFP expression of in in vitro-produced blastocysts derived with sperm from the RFP TG cat. These results indicate that the RFP TG cat exhibits normal reproductive fertility, stable germ-line transmission of the RFP transgene, and characteristic RFP expression in its offspring. We isolated feline neural progenitor cells from a 45-day-old fetus derived from the natural mating of the RFP TG cat with a domestic female cat. Isolated brain and retinal progenitor cells were successfully passaged at least four times post isolation (day 23), and showed a high RFP expression level. This method of producing genetically modified cloned cats will be important for generating biomedical models of human diseases.

Chimeric mitochondrial minichromosomes of the human body louse, Pediculus humanus: evidence for homologous and non-homologous recombination

The mitochondrial (mt) genome of the human body louse, Pediculus humanus, consists of 18 minichromosomes. Each minichromosome is 3 to 4 kb long and has 1 to 3 genes. There is unequivocal evidence for recombination between different mt minichromosomes in P. humanus. It is not known, however, how these minichromosomes recombine. Here, we report the discovery of eight chimeric mt minichromosomes in P. humanus. We classify these chimeric mt minichromosomes into two groups: Group I and Group II. Group I chimeric minichromosomes contain parts of two different protein-coding genes that are from different minichromosomes. The two parts of protein-coding genes in each Group I chimeric minichromosome are joined at a microhomologous nucleotide sequence; microhomologous nucleotide sequences are hallmarks of non-homologous recombination. Group II chimeric minichromosomes contain all of the genes and the non-coding regions of two different minichromosomes. The conserved sequence blocks in the non-coding regions of Group II chimeric minichromosomes resemble the "recombination repeats" in the non-coding regions of the mt genomes of higher plants. These repeats are essential to homologous recombination in higher plants. Our analyses of the nucleotide sequences of chimeric mt minichromosomes indicate both homologous and non-homologous recombination between minichromosomes in the mitochondria of the human body louse.

Nuclear copies of mitochondrial genes: another problem for ancient DNA

The application of ancient DNA techniques is subject to many problems caused primarily by low quality and by low quantity of DNA. For these reasons most studies employing ancient DNA rely on the characterization of mitochondrial DNA, which is present in many more copies per cell than nuclear DNA and hence more copies are likely to survive. We used universal and taxon specific mitochondrial primers to amplify DNA from museum specimens, and found many instances where the amplification of nuclear copies of the mitochondrial gene (numts) instead of the targeted mitochondrial fragment had occurred. Furthermore, the likelihood of amplifying numts increased dramatically when universal primers were utilized. Here we suggest that ancient DNA practitioners must consider the possibility that numts can be amplified at higher rates than previously thought. This is another complication for ancient DNA studies, but it also suggests that more extensive inclusion of nuclear markers in ancient DNA studies should be feasible.

Transportin mediates nuclear entry of DNA in vertebrate systems

Delivery of DNA to the cell nucleus is an essential step in many types of viral infection, transfection, gene transfer by the plant pathogen Agrobacterium tumefaciens and in strategies for gene therapy. Thus, the mechanism by which DNA crosses the nuclear pore complex (NPC) is of great interest. Using nuclei reconstituted in vitro in Xenopus egg extracts, we previously studied DNA passage through the nuclear pores using a single-molecule approach based on optical tweezers. Fluorescently labeled DNA molecules were also seen to accumulate within nuclei. Here we find that this import of DNA relies on a soluble protein receptor of the importin family. To identify this receptor, we used different pathway-specific cargoes in competition studies as well as pathway-specific dominant negative inhibitors derived from the nucleoporin Nup153. We found that inhibition of the receptor transportin suppresses DNA import. In contrast, inhibition of importin beta has little effect on the nuclear accumulation of DNA. The dependence on transportin was fully confirmed in assays using permeabilized HeLa cells and a mammalian cell extract. We conclude that the nuclear import of DNA observed in these different vertebrate systems is largely mediated by the receptor transportin. We further report that histones, a known cargo of transportin, can act as an adaptor for the binding of transportin to DNA.

In situ population structure and ex situ representation of the endangered Amur tiger

The Amur tiger (Panthera tigris altaica) is a critically endangered felid that suffered a severe demographic contraction in the 1940s. In this study, we sampled 95 individuals collected throughout their native range to investigate questions relative to population genetic structure and demographic history. Additionally, we sampled targeted individuals from the North American ex situ population to assess the genetic representation found in captivity. Population genetic and Bayesian structure analyses clearly identified two populations separated by a development corridor in Russia. Despite their well-documented 20th century decline, we failed to find evidence of a recent population bottleneck, although genetic signatures of a historical contraction were detected. This disparity in signal may be due to several reasons, including historical paucity in population genetic variation associated with postglacial colonization and potential gene flow from a now extirpated Chinese population. Despite conflicting signatures of a bottleneck, our estimates of effective population size (N(e) = 27-35) and N(e)/N ratio (0.07-0.054) were substantially lower than the only other values reported for a wild tiger population. Lastly, the extent and distribution of genetic variation in captive and wild populations were similar, yet gene variants persisted ex situ that were lost in situ. Overall, our results indicate the need to secure ecological connectivity between the two Russian populations to minimize loss of genetic diversity and overall susceptibility to stochastic events, and support a previous study suggesting that the captive population may be a reservoir of gene variants lost in situ.

Mitochondrial phylogeography illuminates the origin of the extinct caspian tiger and its relationship to the amur tiger

The Caspian tiger (Panthera tigris virgata) flourished in Central Asian riverine forest systems in a range disjunct from that of other tigers, but was driven to extinction in 1970 prior to a modern molecular evaluation. For over a century naturalists puzzled over the taxonomic validity, placement, and biogeographic origin of this enigmatic animal. Using ancient-DNA (aDNA) methodology, we generated composite mtDNA haplotypes from twenty wild Caspian tigers from throughout their historic range sampled from museum collections. We found that Caspian tigers carry a major mtDNA haplotype differing by only a single nucleotide from the monomorphic haplotype found across all contemporary Amur tigers (P. t. altaica). Phylogeographic analysis with extant tiger subspecies suggests that less than 10,000 years ago the Caspian/Amur tiger ancestor colonized Central Asia via the Gansu Corridor (Silk Road) from eastern China then subsequently traversed Siberia eastward to establish the Amur tiger in the Russian Far East. The conservation implications of these findings are far reaching, as the observed genetic depletion characteristic of modern Amur tigers likely reflects these founder migrations and therefore predates human influence. Also, due to their evolutionary propinquity, living Amur tigers offer an appropriate genetic source should reintroductions to the former range of the Caspian tiger be implemented.

Development and testing of an optimized method for DNA-based identification of jaguar (Panthera onca) and puma (Puma concolor) faecal samples for use in ecological and genetic studies

The elusive nature and endangered status of most carnivore species imply that efficient approaches for their non-invasive sampling are required to allow for genetic and ecological studies. Faecal samples are a major potential source of information, and reliable approaches are needed to foster their application in this field, particularly in areas where few studies have been conducted. A major obstacle to the reliable use of faecal samples is their uncertain species-level identification in the field, an issue that can be addressed with DNA-based assays. In this study we describe a sequence-based approach that efficiently distinguishes jaguar versus puma scats, and that presents several desirable properties: (1) considerably high amplification and sequencing rates; (2) multiple diagnostic sites reliably differentiating the two focal species; (3) high information content that allows for future application in other carnivores; (4) no evidence of amplification of prey DNA; and (5) no evidence of amplification of a nuclear mitochondrial DNA insertion known to occur in the jaguar. We demonstrate the reliability and usefulness of this approach by evaluating 55 field-collected samples from four locations in the highly fragmented Atlantic Forest biome of Brazil and Argentina, and document the presence of one or both of these endangered felids in each of these areas.