Isolation, genome analysis and comparison of a novel parainfluenza virus 5 from a Siberian tiger (Panthera tigris).

Paramyxoviruses are important pathogens affecting various animals, including mammals and humans. Parainfluenza virus 5 (PIV5)-a member of the family Paramyxoviridae-is a major threat to the health of mammals and humans. However, studies on terrestrial wild animals infected with PIV5 are scanty. In this study, we utilized reverse transcription PCR to detect PIV5 infection in the visceral organ tissues of a Siberian tiger (Panthera tigris ssp. altaica) with vomiting, diarrhea, and dyspnea before its death. A novel PIV5 (named SR strain) with a slowly progressive cytopathic effect was isolated in Vero cells and validated using a transmission electron microscope. Full-length sequencing and analysis revealed that the whole genome of the PIV5 SR strain contained 15,246 nucleotides (nt) and seven non-overlapping genes (3'-N-V/P-M-F-SH-HN-L-5') encoding eight proteins. Phylogenetic analysis of three PIV5 strains identified in the same zoo confirmed that PIV5 strains SR and ZJQ-221 shared the closest genetic relationship as they were clustered in the same branch, while the recently found Siberian tiger strain SZ2 kept a certain distance and formed a relatively unique branch. Furthermore, mutations of nt and amino acids (aa) between strains ZJQ-221, SR, and SZ2 were identified. In summary, we report the identification and genomic characterization of a novel PIV5 strain SR isolated in a Siberian tiger, which may help future research on interspecific transmission mechanisms.

Molecular characterization of a novel Spiruromorpha species in wild Chinese pangolin by mitogenome sequence analysis.

Pangolins are susceptible to a variety of gastrointestinal nematodes due to their burrowing lifestyle and feeding habits, and few parasitic nematodes have been reported. Here, a Chinese pangolin with old wounds on its leg and tail was rescued from the Heyuan City, Guangdong Province. The cox1 and SSU rRNA of the worms from the intestine of the Chinese pangolin had the highest sequence identity of 89.58% and 97.95% to the species in the infraorder Spiruromorpha. The complete mitogenome of the worm was further assembled by next-generation sequencing, with a size of 13,708 bp and a GC content of 25.6%. The worm mitogenome had the highest sequence identity of 78.56% to that of Spirocerca lupi, sharing the same gene arrangement with S. lupi and some species in other families under Spiruromorpha. However, the mitogenome between the worm and S. lupi showed differences in codon usage of PCGs, sequences of NCR, and tRNA secondary structures. Phylogenetic analysis showed that the worm mitogenome was clustered with S. lupi in the family Thelaziidae to form a separate branch. However, it is still difficult to identify the worm in the family Thelaziidae because the species in the family Thelaziidae are confused, specifically S. lupi and Thelazia callipaeda in the family Thelaziidae were separated and grouped with species from other families. Thus, the parasitic nematode from the Chinese pangolin may be a novel species in Spiruromorpha and closely related to S. lupi. This study enriches the data on gastrointestinal nematodes in the Chinese pangolin.

Genomic evolution of island birds from the view of the Swinhoe's pheasant (Lophura swinhoii).

Island endemic birds account for the majority of extinct vertebrates in the past few centuries. To date, the evolutionary characteristics of island endemic bird's is poorly known. In this research, we de novo assembled a high-quality chromosome-level reference genome for the Swinhoe's pheasant, which is a typical endemic island bird. Results of collinearity tests suggest rapid ancient chromosome rearrangement that may have contributed to the initial species radiation within Phasianidae, and a role for the insertions of CR1 transposable elements in rearranging chromosomes in Phasianidae. During the evolution of the Swinhoe's pheasant, natural selection positively selected genes involved in fecundity and body size functions, at both the species and population levels, which reflect genetic variation associated with island adaptation. We further tested for variation in population genomic traits between the Swinhoe's pheasant and its phylogenetically closely related mainland relative the silver pheasant, and found higher levels of genetic drift and inbreeding in the Swinhoe's pheasant genome. Divergent demographic histories of insular and mainland bird species during the last glacial period may reflect the differing impact of insular and continental climates on the evolution of species. Our research interprets the natural history and population genetic characteristics of the insular endemic bird the Swinhoe's pheasant, at a genome-wide scale, provides a broader perspective on insular speciation, and adaptive evolution and contributes to the genetic conservation of island endemic birds.

A multilocus DNA mini-barcode assay to identify twenty vertebrate wildlife species.

The world faces significant challenges in preserving the diversity of vertebrate species due to wildlife crimes. DNA barcoding, an effective molecular marker for insufficient nuclear DNA, is an authentic and quick identification technique to trace the origin of seized samples in forensic investigations. Here, we present a multiplex assay capable of identifying twenty vertebrate wildlife species utilizing twenty species-specific primers that target short fragments of the mitochondrial Cyt b, COI, 16S rRNA, and 12S rRNA genes. The assay achieved strong species specificity and sensitivity with a detection limit as low as 5 pg of DNA input. Additionally, it effectively discriminated a minor contributor (≥1%) from binary mixtures and successfully identified of noninvasive samples, inhibited DNA samples, artificially degraded DNA samples, and case samples, demonstrating a sensitive, robust, practical and easily interpretable tool in screening, and investigating forensic wildlife crimes.

Comparative genomics reveals the hybrid origin of a macaque group.

Although species can arise through hybridization, compelling evidence for hybrid speciation has been reported only rarely in animals. Here, we present phylogenomic analyses on genomes from 12 macaque species and show that the fascicularis group originated from an ancient hybridization between the sinica and silenus groups ~3.45 to 3.56 million years ago. The X chromosomes and low-recombination regions exhibited equal contributions from each parental lineage, suggesting that they were less affected by subsequent backcrossing and hence could have played an important role in maintaining hybrid integrity. We identified many reproduction-associated genes that could have contributed to the development of the mixed sexual phenotypes characteristic of the fascicularis group. The phylogeny within the silenus group was also resolved, and functional experimentation confirmed that all extant Western silenus species are susceptible to HIV-1 infection. Our study provides novel insights into macaque evolution and reveals a hybrid speciation event that has occurred only very rarely in primates.

Pathogenicity, tissue tropism and potential vertical transmission of SARSr-CoV-2 in Malayan pangolins.

Malayan pangolin SARS-CoV-2-related coronavirus (SARSr-CoV-2) is closely related to SARS-CoV-2. However, little is known about its pathogenicity in pangolins. Using CT scans we show that SARSr-CoV-2 positive Malayan pangolins are characterized by bilateral ground-glass opacities in lungs in a similar manner to COVID-19 patients. Histological examination and blood gas tests are indicative of dyspnea. SARSr-CoV-2 infected multiple organs in pangolins, with the lungs the major target, and histological expression data revealed that ACE2 and TMPRSS2 were co-expressed with viral RNA. Transcriptome analysis indicated that virus-positive pangolins were likely to have inadequate interferon responses, with relative greater cytokine and chemokine activity in the lung and spleen. Notably, both viral RNA and viral proteins were detected in three pangolin fetuses, providing initial evidence for vertical virus transmission. In sum, our study outlines the biological framework of SARSr-CoV-2 in pangolins, revealing striking similarities to COVID-19 in humans.

Virus diversity, wildlife-domestic animal circulation and potential zoonotic viruses of small mammals, pangolins and zoo animals.

Wildlife is reservoir of emerging viruses. Here we identified 27 families of mammalian viruses from 1981 wild animals and 194 zoo animals collected from south China between 2015 and 2022, isolated and characterized the pathogenicity of eight viruses. Bats harbor high diversity of coronaviruses, picornaviruses and astroviruses, and a potentially novel genus of Bornaviridae. In addition to the reported SARSr-CoV-2 and HKU4-CoV-like viruses, picornavirus and respiroviruses also likely circulate between bats and pangolins. Pikas harbor a new clade of Embecovirus and a new genus of arenaviruses. Further, the potential cross-species transmission of RNA viruses (paramyxovirus and astrovirus) and DNA viruses (pseudorabies virus, porcine circovirus 2, porcine circovirus 3 and parvovirus) between wildlife and domestic animals was identified, complicating wildlife protection and the prevention and control of these diseases in domestic animals. This study provides a nuanced view of the frequency of host-jumping events, as well as assessments of zoonotic risk.

Population genomic analysis provides evidence of the past success and future potential of South China tiger captive conservation.

BACKGROUND: Among six extant tiger subspecies, the South China tiger (Panthera tigris amoyensis) once was widely distributed but is now the rarest one and extinct in the wild. All living South China tigers are descendants of only two male and four female wild-caught tigers and they survive solely in zoos after 60 years of effective conservation efforts. Inbreeding depression and hybridization with other tiger subspecies were believed to have occurred within the small, captive South China tiger population. It is therefore urgently needed to examine the genomic landscape of existing genetic variation among the South China tigers. RESULTS: In this study, we assembled a high-quality chromosome-level genome using long-read sequences and re-sequenced 29 high-depth genomes of the South China tigers. By combining and comparing our data with the other 40 genomes of six tiger subspecies, we identified two significantly differentiated genomic lineages among the South China tigers, which harbored some rare genetic variants introgressed from other tiger subspecies and thus maintained a moderate genetic diversity. We noticed that the South China tiger had higher FROH values for longer runs of homozygosity (ROH > 1 Mb), an indication of recent inbreeding/founder events. We also observed that the South China tiger had the least frequent homozygous genotypes of both high- and moderate-impact deleterious mutations, and lower mutation loads than both Amur and Sumatran tigers. Altogether, our analyses indicated an effective genetic purging of deleterious mutations in homozygous states from the South China tiger, following its population contraction with a controlled increase in inbreeding based on its pedigree records. CONCLUSIONS: The identification of two unique founder/genomic lineages coupled with active genetic purging of deleterious mutations in homozygous states and the genomic resources generated in our study pave the way for a genomics-informed conservation, following the real-time monitoring and rational exchange of reproductive South China tigers among zoos.

Biological characteristics of two mesenchymal stem cell cultures isolated from the umbilical cord and adipose tissue of a neonatal common hippo (Hippopotamus amphibius).

Mesenchymal stem cells (MSCs) are multipotent adult stem cells and can be isolated from many tissues of the body. Due to their potentials to treat various diseases and be applied in animal breeding, MSCs have been isolated and identified regarding their biological properties. Common hippos (Hippopotamus amphibius) are a vulnerable species and yet the cryopreservation of their genetic materials is scare. In this study, we successfully established two MSC cultures (UC-MSCs and AT-MSCs) from the umbilical cord and adipose tissue of a neonatal common hippo and comparatively described their features. Both UC-MSCs and AT-MSCs showed fibroblastoid morphology and could be continuously passaged for over 17 passages without dramatic signs of senescence. The cell cultures had normal chromosome composition, say, 17 pairs of autosomes and 1 pair of X chromosomes. UC-MSCs and AT-MSCs displayed similar gene expression profiles. They were positive for CD34, CD45, CD73, CD90 and CD105 and negative for HLA-DR. They demonstrated stemness maintenance by expression of classical stem cell markers. UC-MSCs and AT-MSCs manifested different differentiation potentials into other cell lineages. In summary, these two cell cultures demonstrated the essential properties of mesenchymal stem cells and might play a role in the future research.

Heteroplasmic mitochondrial genomes of a Raillietina tapeworm in wild Pangolin.

BACKGROUND: Raillietina species belong to the family Davaineidae, which parasitizes in a wide variety of mammals and birds, causing stunted growth, lethargy, emaciation, and digestive tract obstruction. However, only a limited number of Raillietina species have been identified in wild animals. METHODS: We analyzed and annotated the complete mitochondrial (mt) genome of a worm from the intestine of a wild pangolin using Illumina sequencing of whole genomic DNA. RESULTS: These findings showed the presence of two mtDNA sequences in Raillietina sp., designated as mt1 and mt2, with the lengths of 14,331 bp and 14,341 bp, respectively. Both the mts genomes of Raillietina sp. comprised 36 genes, containing 12 protein-coding genes (PCGs), 2 ribosomal RNAs, and 22 transfer RNAs. Gene arrangements of both mt genomes of Raillietina sp. were similar to those of most flatworms, except for taeniids, which shift positions between tRNAL1 and tRNAS2 genes. Twenty of 22 tRNA secondary structures of Raillietina sp. had a typical cloverleaf structure similar to Raillietina tetragona. Sequence differences between the mt1 and mt2 genomes were 4.4%, and this difference arises from the mtDNA heteroplasmic mutations. Moreover, heteroplasmic mtDNA mutations were detected in PCGs, tRNAs, rRNAs, NCRs, and intergenes, but the highest proportion of heteroplasmy of 79.0% was detected in PCGs, indicating the occurrence of mtDNA heteroplasmy in Raillietina sp. To our knowledge, this is the first report of mtDNA heteroplasmy in tapeworm parasites. Phylogenetic analyses of 18S rRNA, ITS2, and 12 PCG sequences demonstrated that the worm was clustered with other Raillietina species in the Davaneidae family. CONCLUSIONS: We found a novel Raillietina species in wild pangolin with the existence of mitochondrial DNA heteroplasmy. Thus, these findings provide insights into the heterogeneity of the mt genome in parasitic cestodes, and mt genome data contributes to the understanding of pangolin-parasitic cestodes in terms of their molecular biology, epidemiology, diagnosis, and taxonomy.

Twelve natural estrogens in urines of six threatened or endangered mammalian species in Zoo Park: implications and their potential risk.

This work was the first to report twelve natural estrogens (NEs) in the urines of six threatened or endangered mammalians in a Zoo Park of Guangzhou (i.e., panda, gorilla, elephant, African lion, jaguar, and leopard). Ten out of twelve NEs were detected at least in one urine sample of the six mammalians studied, including the four major NEs (i.e., estrone (E1), 17ß-estradiol (E2), 17α-estradiol (αE2), estriol (E3)), and six other NEs (i.e., 4-hydroxyestrone (4OHE1), 2-hydroxyestradiol (2OHE2), 4-hydroxyestradiol (4OHE2), 16α-hydroxyestrone (16α-OHE1), 16ketoestradiol (16ketoE2), and 17epiestriol (17epiE3)). The six studied mammalians, ranked in the order of high to low urinary concentration of total NEs, were jaguar, African lion, gorilla, elephant, panda, and leopard, with respective urinary concentrations of 110.4, 86.4, 71.4, 66.0, 55.9, and 52.8 ng/mL. According to the average urinary concentration of NE in the six mammalians ranked from high to low, the top five NEs detected were 16α-OHE1, 4OHE1, E1, E3, and 17epiE3, respectively. These clearly indicated the occurrence of NEs other than the four major types in urines of animals in a Zoo Park. Moreover, the daily excretion rates of the five detected NEs by one elephant ranged from 1162-2254 µg/d with a total daily excretion rate of 8260 µg/d, suggesting that the total urinary excretion of NEs by one adult elephant was equivalent to that by 170 premenopausal women or 506 adult men. Consequently, it appears from this study that NEs in the urines of zoo animals should be considered an emerging source of NEs.

Molecular detection of a novel Ancylostoma sp. by whole mtDNA sequence from pangolin Manis javanica.

BACKGROUND: Ancylostoma species are hematophagous parasites that cause chronic hemorrhage in various animals and humans. Pangolins, also known as scaly anteaters, are mammals that live in soil environments where they are readily exposed to soil-borne parasitic nematodes. However, only a limited number of helminth species have been identified in this animal host so far. METHODS: Ancylostoma sp. was isolated from a wild pangolin, and the complete mitochondrial (mt) genome of Ancylostoma sp. was obtained by Illumina sequencing of total genomic DNA. RESULTS: The circular complete mt genome that was assembled had a total length of 13,757 bp and comprised 12 protein-coding genes (PCGs), 22 transfer ribosomal RNAs, two ribosomal RNAs (rRNAs), two non-coding regions and one AT-rich region, but lacked the gene coding for ATPase subunit 8 (atp8). The overall AT content of the mt genome of Ancylostoma sp. was 76%, which is similar to that of other nematodes. The PCGs used two start codons (ATT and TTG) and three stop codons (TAA, TAG, and T). The nucleotide identity of the 12 PCGs ranged from 83.1% to 89.7% and had the highest sequence identity with Ancylostoma caninum among species in the Ancylostomatidae family. Also, the pangolin-derived Ancylostoma sp. lacked repeat sequences in the non-coding regions and in the unique sequence of the short non-coding regions, which differentiated it from other Ancylostoma species. In addition, phylogenetic analyses of 18S rRNA and mtDNA sequences revealed that the Ancylostoma sp. was positioned in a separate branch in the subfamily Ancylostomatinae along with other Ancylostoma species. CONCLUSIONS: The Ancylostoma sp. isolated from a pangolin in this study was identified as a possible new Ancylostoma species. The identification of this Ancylostoma sp. from pangolin enriches our knowledge of the species in the Ancylostomatidae family and provides information that will lead to a better understanding of the taxonomy, diagnostics, and biology of hookworms.

Comparative study of gut microbiota from captive and confiscated-rescued wild pangolins.

Pangolins are among the most critically endangered animals due to widespread poaching and worldwide trafficking. Captive breeding is considered to be one way to protect them and increase the sizes of their populations. However, comparative studies of captive and wild pangolins in the context of gut microbiota are rare. Here, the gut microbiome of captive and confiscated-rescued wild pangolins is compared, and the effects of different periods of captivity and captivity with and without antibiotic treatment are considered. We show that different diets and periods of captivity, as well as the application of antibiotic therapy, can alter gut community composition and abundance in pangolins. Compared to wild pangolins, captive pangolins have an increased capacity for chitin and cellulose/hemicellulose degradation, fatty acid metabolism, and short-chain fatty acid synthesis, but a reduced ability to metabolize exogenous substances. In addition to increasing the ability of the gut microbiota to metabolize nutrients in captivity, captive breeding imposes some risks for survival by resulting in a greater abundance of antibiotic resistance genes and virulence factors in captive pangolins than in wild pangolins. Our study is important for the development of guidelines for pangolin conservation, including health assessment, disease prevention, and rehabilitation of wild pangolin populations.

Covariation of the Fecal Microbiome with Diet in Nonpasserine Birds.

Opportunistic feeding and multiple other environment factors can modulate the gut microbiome, and bias conclusions, when wild animals are used for studying the influence of phylogeny and diet on their gut microbiomes. Here, we controlled for these other confounding factors in our investigation of the magnitude of the effect of diet on the gut microbiome assemblies of nonpasserine birds. We collected fecal samples, at one point in time, from 35 species of birds in a single zoo as well as 6 species of domestic poultry from farms in Guangzhou city to minimize the influences from interfering factors. Specifically, we describe 16S rRNA amplicon data from 129 fecal samples obtained from 41 species of birds, with additional shotgun metagenomic sequencing data generated from 16 of these individuals. Our data show that diets containing native starch increase the abundance of Lactobacillus in the gut microbiome, while those containing plant-derived fiber mainly enrich the level of Clostridium Greater numbers of Fusobacteria and Proteobacteria are detected in carnivorous birds, while in birds fed a commercial corn-soybean basal diet, a stronger inner-connected microbial community containing Clostridia and Bacteroidia was enriched. Furthermore, the metagenome functions of the microbes (such as lipid metabolism and amino acid synthesis) were adapted to the different food types to achieve a beneficial state for the host. In conclusion, the covariation of diet and gut microbiome identified in our study demonstrates a modulation of the gut microbiome by dietary diversity and helps us better understand how birds live based on diet-microbiome-host interactions.IMPORTANCE Our study identified food source, rather than host phylogeny, as the main factor modulating the gut microbiome diversity of nonpasserine birds, after minimizing the effects of other complex interfering factors such as weather, season, and geography. Adaptive evolution of microbes to food types formed a dietary-microbiome-host interaction reciprocal state. The covariation of diet and gut microbiome, including the response of microbiota assembly to diet in structure and function, is important for health and nutrition in animals. Our findings help resolve the major modulators of gut microbiome diversity in nonpasserine birds, which had not previously been well studied. The diet-microbe interactions and cooccurrence patterns identified in our study may be of special interest for future health assessment and conservation in birds.

Identification of Amblyomma javanense and detection of tick-borne Ehrlichia spp. in confiscated Malayan Pangolins.

Due to habitat destruction and illegal hunting and trade, the number of pangolins has been sharply reduced. To protect pangolins from extinction, relevant departments are combined and active action have been taken. A total of 21 confiscated Malayan pangolins were rescued in 2019, but died continuously for unknown reasons. This study aimed to investigate the reasons for the death of these pangolin and rescue them. 19 of the 21 confiscated pangolins had ticks on their body integument. A total of 303 ticks were collected and identified as Amblyomma javanense (A. javanense) according to their morphology and the sequences of 16S rRNA and internal transcribed spacer 2 (ITS2). There were multi-organ damages in the dead pangolins, especially congestion and hemorrhage in lung, heart and kidney and inflammation of which were observed using HE staining. Pathogens' nucleic acid detection showed ticks were only positive for Ehrlichia spp, with 56.7% positive rate of collected ticks (127/224), which was further confirmed in tissues from dead pangolins. Our findings confirm that ehrlichiosis caused by Ehrlichia spp. from A. javanense might accelerate the confiscated pangolin's death. More attention should be payed to tick-elimination work and the diagnoses and treatment of tick-borne diseases in the follow-up rescue operation.

Author Correction: 863 genomes reveal the origin and domestication of chicken.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

863 genomes reveal the origin and domestication of chicken.

Despite the substantial role that chickens have played in human societies across the world, both the geographic and temporal origins of their domestication remain controversial. To address this issue, we analyzed 863 genomes from a worldwide sampling of chickens and representatives of all four species of wild jungle fowl and each of the five subspecies of red jungle fowl (RJF). Our study suggests that domestic chickens were initially derived from the RJF subspecies Gallus gallus spadiceus whose present-day distribution is predominantly in southwestern China, northern Thailand and Myanmar. Following their domestication, chickens were translocated across Southeast and South Asia where they interbred locally with both RJF subspecies and other jungle fowl species. In addition, our results show that the White Leghorn chicken breed possesses a mosaic of divergent ancestries inherited from other subspecies of RJF. Despite the strong episodic gene flow from geographically divergent lineages of jungle fowls, our analyses show that domestic chickens undergo genetic adaptations that underlie their unique behavioral, morphological and reproductive traits. Our study provides novel insights into the evolutionary history of domestic chickens and a valuable resource to facilitate ongoing genetic and functional investigations of the world's most numerous domestic animal.

Isolation of SARS-CoV-2-related coronavirus from Malayan pangolins.

The current outbreak of coronavirus disease-2019 (COVID-19) poses unprecedented challenges to global health1. The new coronavirus responsible for this outbreak-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-shares high sequence identity to SARS-CoV and a bat coronavirus, RaTG132. Although bats may be the reservoir host for a variety of coronaviruses3,4, it remains unknown whether SARS-CoV-2 has additional host species. Here we show that a coronavirus, which we name pangolin-CoV, isolated from a Malayan pangolin has 100%, 98.6%, 97.8% and 90.7% amino acid identity with SARS-CoV-2 in the E, M, N and S proteins, respectively. In particular, the receptor-binding domain of the S protein of pangolin-CoV is almost identical to that of SARS-CoV-2, with one difference in a noncritical amino acid. Our comparative genomic analysis suggests that SARS-CoV-2 may have originated in the recombination of a virus similar to pangolin-CoV with one similar to RaTG13. Pangolin-CoV was detected in 17 out of the 25 Malayan pangolins that we analysed. Infected pangolins showed clinical signs and histological changes, and circulating antibodies against pangolin-CoV reacted with the S protein of SARS-CoV-2. The isolation of a coronavirus from pangolins that is closely related to SARS-CoV-2 suggests that these animals have the potential to act as an intermediate host of SARS-CoV-2. This newly identified coronavirus from pangolins-the most-trafficked mammal in the illegal wildlife trade-could represent a future threat to public health if wildlife trade is not effectively controlled.

Some maternal lineages of domestic horses may have origins in East Asia revealed with further evidence of mitochondrial genomes and HVR-1 sequences.

OBJECTIVES: There are large populations of indigenous horse (Equus caballus) in China and some other parts of East Asia. However, their matrilineal genetic diversity and origin remained poorly understood. Using a combination of mitochondrial DNA (mtDNA) and hypervariable region (HVR-1) sequences, we aim to investigate the origin of matrilineal inheritance in these domestic horses. METHODS: To investigate patterns of matrilineal inheritance in domestic horses, we conducted a phylogenetic study using 31 de novo mtDNA genomes together with 317 others from the GenBank. In terms of the updated phylogeny, a total of 5,180 horse mitochondrial HVR-1 sequences were analyzed. RESULTS: Eightteen haplogroups (Aw-Rw) were uncovered from the analysis of the whole mitochondrial genomes. Most of which have a divergence time before the earliest domestication of wild horses (about 5,800 years ago) and during the Upper Paleolithic (35-10 KYA). The distribution of some haplogroups shows geographic patterns. The Lw haplogroup contained a significantly higher proportion of European horses than the horses from other regions, while haplogroups Jw, Rw, and some maternal lineages of Cw, have a higher frequency in the horses from East Asia. The 5,180 sequences of horse mitochondrial HVR-1 form nine major haplogroups (A-I). We revealed a corresponding relationship between the haplotypes of HVR-1 and those of whole mitochondrial DNA sequences. The data of the HVR-1 sequences also suggests that Jw, Rw, and some haplotypes of Cw may have originated in East Asia while Lw probably formed in Europe. CONCLUSIONS: Our study supports the hypothesis of the multiple origins of the maternal lineage of domestic horses and some maternal lineages of domestic horses may have originated from East Asia.