Phylogeographic study of the Bufo gargarizans species complex, with emphasis on Northeast Asia.

We conduct a phylogeographic and population genetic study of the Asiatic toad (Bufo gargarizans) to understand its evolutionary history, and the influence of geology and climate. A total of 292 individuals from 94 locations were genotyped for two mitochondrial loci (cytb, ND2) and five nuclear introns (Sox9-2, Rho-3, CCNB2-3, UCH-2, and DBI-2), and we performed a suite of phylogenetic, population genetic, and divergence dating analyses. The phylogenetic trees constructed using mitochondrial loci inferred B. gargarizans being divided into two major groups: China mainland and Northeast Asia (Northeast China, Russia, and Korean Peninsula). As with previous studies of this species, we recover population genetic structure not tied to geographic region. Additionally, we discover a new genetic clade restricted to Northeast Asia that points towards the Korean Peninsula being a glacial refugium during the Pleistocene. The weak phylogeographic pattern of B. gargarizans is likely the result of multiple biological, anthropogenic, and historical factors - robust dispersal abilities as a consequence of physiological adaptations, human translocation, geologic activity, and glacial cycles of the Pleistocene. We highlight the complex geologic and climatic history of Northeast Asia and encourage further research to understand its impact on the biodiversity in the region.

Asia-wide phylogeography of wild boar (Sus scrofa) based on mitochondrial DNA and Y-chromosome: Revising the migration routes of wild boar in Asia.

Genetics of pigs has been well studied in Europe and Asia, but most of previous studies of molecular phylogeny of Sus scrofa have been based on sequences of both wild and domestic forms. In this study we analysed genetic traits of Sus scrofa from 13 regions in Asia (including previously undisclosed Eastern Caucasus and Trans-Baikal regions) using purely wild boar samples. Mitochondrial control region and Y-chromosome genes (AMELY & USP9Y) were employed to resolve phylogeographic relationships. We discussed spatio-temporal dynamics of wild boar distribution and compared molecular data to morphological and cytogenetic data on wild boar variability and taxonomy. A total of 51 haplotypes were detected in mtDNA control region and five haplotypes were found in combined sequences of Y-chromosome genes. The phylogeography of Asia-wide wild boars supported a hypothesis of migration from South-East Asia to South Asia, followed by migration to East and West Asia. We present a hypothesis about independent dispersal of wild boars into West Asia from South and North-East Asia. Mitochondrial DNA phylogeny generally fits the morphologically based intraspecies taxonomy. Distribution of chromosomal variants of wild boar presently does not show clear correlation with mtDNA clades.

Tandem duplication of mitochondrial DNA in the black-faced spoonbill, Platalea minor.

Mitochondrial (mt) heteroplasmy in the control region (CR) of the black-faced spoonbill was investigated using LA-PCR. To avoid amplification of transpositioned nuclear genome fragment from mtDNA (numt), PCR product of the almost-complete mitochondrial genome was amplified using primers designed to anneal on the COIII gene. Then nested LA-PCR product was amplified between the cyt b and 12S rRNA genes using the almost-complete mitochondrial genome PCR product as a template. Nucleotide sequencing revealed tandem duplication composed of two units. The first contains cyt b-1, tRNA(Thr)-1, tRNA(Pro)-1, ND6-1, tRNA(Glu)-1 and CR1, and the second consists of cyt b-2, tRNA(Thr)-2, tRNA(Pro)-2, ND6-2, tRNA(Glu)-2 and CR2, followed by tRNA(Phe) and 12S rRNA. The duplicated cyt b-2 sequence coincided with 499 bp at the 3' end of cyt b-1. With the exception of the CR, the other genes in the duplicated sequence were identical to the original corresponding gene. Even though both CR1 and CR2 contain functional blocks, such as a poly-C site, a goose hairpin and a TAS structure in Domain I, the 3' end of CR1 was followed by a 112 bp sequence (non-coding region) that was not found in CR2 or in sequence homology analysis of similar genes. Meanwhile, CR2 ended in a complicated repeat sequence. The 5' franking region in the Domain I (Region A) and the 3' franking region in the Domain I (Region B) of the two CRs evolve in quite different manners: Region A was highly variable between CR1 and CR2 in the same individuals, while Region B was almost identical between them, which indicates concerted evolution.