A First Y-Chromosomal Haplotype Network to Investigate Male-Driven Population Dynamics in Domestic and Wild Bactrian Camels

Camelus knoblochi genome reveals the complex evolutionary history of Old World camels

Extant Old World camels (genus Camelus) contributed to the economic and cultural exchanges between the East and West for thousands of years.1,2 Although many remains have been unearthed,3,4,5 we know neither whether the prevalent hybridization observed between extant Camelus species2,6,7 also occurred between extinct lineages and the ancestors of extant Camelus species nor why some populations became extinct while others survived. To investigate these questions, we generated paleogenomic and stable isotope data from an extinct two-humped camel species, Camelus knoblochi. We find that in the mitochondrial phylogeny, all C. knoblochi form a paraphyletic group that nests within the diversity of modern, wild two-humped camels (Camelus ferus). In contrast, they are clearly distinguished from both wild and domesticated (Camelus bactrianus) two-humped camels on the nuclear level. Moreover, the divergence pattern of the three camel species approximates a trifurcation, because the most common topology is only slightly more frequent than the two other possible topologies. This mito-nuclear phylogenetic discordance likely arose due to interspecific gene flow between all three species, suggesting that interspecific hybridization is not exclusive to modern camels but a recurrent phenomenon throughout the evolutionary history of the genus Camelus. These results suggest that the genomic complexity of Old World camels' evolutionary history is underestimated when considering data from only modern species. Finally, we find that C. knoblochi populations began declining prior to the last glacial maximum and, by integrating palaeoecological evidence and stable isotope data, suggest that this was likely due to failure to adapt to a changing environment.

Inferences about the population history of Rangifer tarandus from Y chromosome and mtDNA phylogenies

Reindeer, called caribou in North America, has a circumpolar distribution and all extant populations belong to the same species (Rangifer tarandus). It has survived the Holocene thanks to its immense adaptability and successful coexistence with humans in different forms of hunting and herding cultures. Here, we examine the paternal and maternal history of Rangifer based on robust Y-chromosomal and mitochondrial DNA (mtDNA) trees representing Eurasian tundra reindeer, Finnish forest reindeer, Svalbard reindeer, Alaska tundra caribou, and woodland caribou. We first assembled Y-chromosomal contigs, representing 1.3 Mb of single-copy Y regions. Based on 545 Y-chromosomal and 458 mtDNA SNPs defined in 55 males, maximum parsimony trees were created. We observed two well separated clades in both phylogenies: the "EuroBeringian clade" formed by animals from Arctic Islands, Eurasia, and a few from North America and the "North American clade" formed only by caribou from North America. The time calibrated Y tree revealed an expansion and dispersal of lineages across continents after the Last Glacial Maximum. We show for the first time unique paternal lineages in Svalbard reindeer and Finnish forest reindeer and reveal a circumscribed Y haplogroup in Fennoscandian tundra reindeer. The Y chromosome in domesticated reindeer is markedly diverse indicating that several male lineages have undergone domestication and less intensive selection on males. This study places R. tarandus onto the list of species with resolved Y and mtDNA phylogenies and builds the basis for studies of the distribution and origin of paternal and maternal lineages in the future.

Refining the evolutionary tree of the horse Y chromosome

The Y chromosome carries information about the demography of paternal lineages, and thus, can prove invaluable for retracing both the evolutionary trajectory of wild animals and the breeding history of domesticates. In horses, the Y chromosome shows a limited, but highly informative, sequence diversity, supporting the increasing breeding influence of Oriental lineages during the last 1500 years. Here, we augment the primary horse Y-phylogeny, which is currently mainly based on modern horse breeds of economic interest, with haplotypes (HT) segregating in remote horse populations around the world. We analyze target enriched sequencing data of 5 Mb of the Y chromosome from 76 domestic males, together with 89 whole genome sequenced domestic males and five Przewalski's horses from previous studies. The resulting phylogeny comprises 153 HTs defined by 2966 variants and offers unprecedented resolution into the history of horse paternal lineages. It reveals the presence of a remarkable number of previously unknown haplogroups in Mongolian horses and insular populations. Phylogenetic placement of HTs retrieved from 163 archaeological specimens further indicates that most of the present-day Y-chromosomal variation evolved after the domestication process that started around 4200 years ago in the Western Eurasian steppes. Our comprehensive phylogeny significantly reduces ascertainment bias and constitutes a robust evolutionary framework for analyzing horse population dynamics and diversity.

McAN: a novel computational algorithm and platform for constructing and visualizing haplotype networks

Haplotype networks are graphs used to represent evolutionary relationships between a set of taxa and are characterized by intuitiveness in analyzing genealogical relationships of closely related genomes. We here propose a novel algorithm termed McAN that considers mutation spectrum history (mutations in ancestry haplotype should be contained in descendant haplotype), node size (corresponding to sample count for a given node) and sampling time when constructing haplotype network. We show that McAN is two orders of magnitude faster than state-of-the-art algorithms without losing accuracy, making it suitable for analysis of a large number of sequences. Based on our algorithm, we developed an online web server and offline tool for haplotype network construction, community lineage determination, and interactive network visualization. We demonstrate that McAN is highly suitable for analyzing and visualizing massive genomic data and is helpful to enhance the understanding of genome evolution. Availability: Source code is written in C/C++ and available at https://github.com/Theory-Lun/McAN and https://ngdc.cncb.ac.cn/biocode/tools/BT007301 under the MIT license. Web server is available at https://ngdc.cncb.ac.cn/bit/hapnet/. SARS-CoV-2 dataset are available at https://ngdc.cncb.ac.cn/ncov/. Contact: songshh@big.ac.cn (Song S), zhaowm@big.ac.cn (Zhao W), baoym@big.ac.cn (Bao Y), zhangzhang@big.ac.cn (Zhang Z), ybxue@big.ac.cn (Xue Y).

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY, DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes Y1A, Y1B, Y2A and Y2B with a marked geographic partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and 7 wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during the migrations into northern Europe, eastern and southern Asia and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19^th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographic range.

What's in a name? Common name misuse potentially confounds the conservation of the wild camel Camelus ferus

Common names allow species diversity to be acknowledged by experts and non-specialists alike; they are descriptors with both scientific and cultural implications. However, a lack of clarity when using a common name could risk altering perceptions of threatened species. This is the case for the Critically Endangered wild camel Camelus ferus, which, despite extensive evidence of its species status, is frequently referred to in English as wild Bactrian camel. However, the wild camel (Mongolian: хавтгай, khavtgai; Chinese: 野骆驼, ye luo tuo) is not a wild version of the domestic Bactrian camel Camelus bactrianus but a separate species near extinction, with an estimated population of c. 950. Failure to clearly separate Bactrian and wild camels in name risks masking the plight of the few remaining wild camels with the visible abundance of the domesticated species. Here we advocate the use of the accurate English common name wild camel for C. ferus ideally alongside its Indigenous names to correctly represent its cultural and conservation importance.

Haplogroup Prediction Using Y-Chromosomal Short Tandem Repeats in the General Population of Bosnia and Herzegovina

Human Y-chromosomal haplogroups are an important tool used in population genetics and forensic genetics. A conventional method used for Y haplogroup assignment is based on a set of Y-single nucleotide polymorphism (SNP) markers deployed, which exploits the low mutation rate nature of these markers. Y chromosome haplogroups can be successfully predicted from Y-short tandem repeat (STR) markers using different software packages, and this method gained much attention recently due to its labor-, time-, and cost-effectiveness. The present study was based on the analysis of a total of 480 adult male buccal swab samples collected from different regions of Bosnia and Herzegovina. Y haplogroup prediction was performed using Whit Athey’s Haplogroup Predictor, based on haplotype data on 23 Y-STR markers contained within the PowerPlex® Y23 kit. The results revealed the existence of 14 different haplogroups, with I2a, R1a, and E1b1b being the most prevalent with frequencies of 43.13, 14.79, and 14.58%, respectively. Compared to the previously published studies on Bosnian-Herzegovinian population based on Y-SNP and Y-STR data, this study represents an upgrade of molecular genetic data with a significantly larger number of samples, thus offering more accurate results and higher probability of detecting rare haplogroups.

An 8.22 Mb Assembly and Annotation of the Alpaca (Vicugna pacos) Y Chromosome

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca (Vicugna pacos) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

Paternal Origins and Migratory Episodes of Domestic Sheep

The domestication and subsequent global dispersal of livestock are crucial events in human history, but the migratory episodes during the history of livestock remain poorly documented [1-3]. Here, we first developed a set of 493 novel ovine SNPs of the male-specific region of Y chromosome (MSY) by genome mapping. We then conducted a comprehensive genomic analysis of Y chromosome, mitochondrial DNA, and whole-genome sequence variations in a large number of 595 rams representing 118 domestic populations across the world. We detected four different paternal lineages of domestic sheep and resolved, at the global level, their paternal origins and differentiation. In Northern European breeds, several of which have retained primitive traits (e.g., a small body size and short or thin tails), and fat-tailed sheep, we found an overrepresentation of MSY lineages y-HC and y-HB, respectively. Using an approximate Bayesian computation approach, we reconstruct the demographic expansions associated with the segregation of primitive and fat-tailed phenotypes. These results together with archaeological evidence and historical data suggested the first expansion of early domestic hair sheep and the later expansion of fat-tailed sheep occurred ∼11,800-9,000 years BP and ∼5,300-1,700 years BP, respectively. These findings provide important insights into the history of migration and pastoralism of sheep across the Old World, which was associated with different breeding goals during the Neolithic agricultural revolution.

Genomic signatures of domestication in Old World camels

Domestication begins with the selection of animals showing less fear of humans. In most domesticates, selection signals for tameness have been superimposed by intensive breeding for economical or other desirable traits. Old World camels, conversely, have maintained high genetic variation and lack secondary bottlenecks associated with breed development. By re-sequencing multiple genomes from dromedaries, Bactrian camels, and their endangered wild relatives, here we show that positive selection for candidate genes underlying traits collectively referred to as 'domestication syndrome' is consistent with neural crest deficiencies and altered thyroid hormone-based signaling. Comparing our results with other domestic species, we postulate that the core set of domestication genes is considerably smaller than the pan-domestication set - and overlapping genes are likely a result of chance and redundancy. These results, along with the extensive genomic resources provided, are an important contribution to understanding the evolutionary history of camels and the genomic features of their domestication.

Whole-genome sequencing of 128 camels across Asia reveals origin and migration of domestic Bactrian camels

The domestic Bactrian camels were treated as one of the principal means of locomotion between the eastern and western cultures in history. However, whether they originated from East Asia or Central Asia remains elusive. To address this question, we perform whole-genome sequencing of 128 camels across Asia. The extant wild and domestic Bactrian camels show remarkable genetic divergence, as they were split from dromedaries. The wild Bactrian camels also contribute little to the ancestry of domestic ones, although they share close habitat in East Asia. Interestingly, among the domestic Bactrian camels, those from Iran exhibit the largest genetic distance and the earliest split from all others in the phylogeny, despite evident admixture between domestic Bactrian camels and dromedaries living in Central Asia. Taken together, our study support the Central Asian origin of domestic Bactrian camels, which were then immigrated eastward to Mongolia where native wild Bactrian camels inhabit.

Ming, Yuan et al. performed whole-genome sequencing on 128 wild and domesticated Bactrian camels across Asia. They show that wild and domestic Bactrian camels are genetically diverged from dromedaries, and that wild camels contributed little to domestic camel ancestry despite sharing a habitat in East Asia.

Population genetic analysis of the domestic Bactrian camel in China by RAD-seq

Restriction site-associated DNA sequencing (RAD-seq) is one of the most effective high-throughput sequencing technologies for SNP development and utilization and has been applied to studying the origin and evolution of various species. The domestic Bactrian camels play an important role in economic trade and cultural construction. They are precious species resources and indispensable animals in China's agricultural production. Recently, the rapid development of modern transportation and agriculture, and the deterioration of the environment have led to a sharp decline in the number of camels. Although there have been some reports on the evolution history of the domestic Bactrian camel in China, the origin, evolutionary relationship, and genetic diversity of the camels are unclear due to the limitations of sample size and sequencing technology. Therefore, 47 samples of seven domestic Bactrian camel species from four regions (Inner Mongolia, Gansu, Qinghai, and Xinjiang) were prepared for RAD-seq analysis to study the evolutionary relationship and genetic diversity. In addition, seven domestic Bactrian camel species are located in different ecological zones, forming different characteristics and having potential development value. A total of 6,487,849 SNPs were genotyped. On the one hand, the filtered SNP information was used to conduct polymorphism mapping construction, LD attenuation analysis, and nucleotide diversity analysis. The results showed that the number of SNPs in Dongjiang camel was the highest, the LD coefficient decayed the fastest, and the nucleotide diversity was the highest. It indicates that Dongjiang camel has the highest genetic diversity. On the other hand, the filtered SNPs information was used to construct the phylogenetic tree, and F ST analysis, inbreeding coefficient analysis, principal component analysis, and population structure analysis were carried out. The results showed that Nanjiang camel and Beijiang camels grouped together, and the other five Bactrian camel populations gathered into another branch. It may be because the mountains in the northern part of Xinjiang and the desert in the middle isolate the two groups from the other five groups.

Old World camels in a modern world - a balancing act between conservation and genetic improvement

Old World camels have served humans in cross‐continental caravans, transporting people and goods, connecting different cultures and providing milk, meat, wool and draught since their domestication around 3000–6000 years ago. In a world of modern transport and fast connectivity, these beasts of burden seem to be out‐dated. However, a growing demand for sustainable milk and meat production, especially in countries affected by climate change and increasing desertification, brings dromedaries (Camelus dromedarius) and Bactrian camels (Camelus bactrianus) back onstage and into the focus of animal breeders and scientists. In this review on the molecular genetics of these economically important species we give an overview about the evolutionary history, domestication and dispersal of Old World camels, whereas highlighting the need for conservation of wild two‐humped camels (Camelus ferus) as an evolutionarily unique and highly endangered species. We provide cutting‐edge information on the current molecular resources and on‐going sequencing projects. We cannot emphasise enough the importance of balancing the need for improving camel production traits with maintaining the genetic diversity in two domestic species with specific physiological adaptation to a desert environment.