Large-Scale Chromosomal Changes Lead to Genome-Level Expression Alterations, Environmental Adaptation, and Speciation in the Gayal (Bos frontalis).

Determining the functional consequences of karyotypic changes is invariably challenging because evolution tends to obscure many of its own footprints, such as accumulated mutations, recombination events, and demographic perturbations. Here, we describe the assembly of a chromosome-level reference genome of the gayal (Bos frontalis) thereby revealing the structure, at base-pair-level resolution, of a telo/acrocentric-to-telo/acrocentric Robertsonian translocation (2;28) (T/A-to-T/A rob[2;28]). The absence of any reduction in the recombination rate or genetic introgression within the fusion region of gayal served to challenge the long-standing view of a role for fusion-induced meiotic dysfunction in speciation. The disproportionate increase noted in the distant interactions across pro-chr2 and pro-chr28, and the change in open-chromatin accessibility following rob(2;28), may, however, have led to the various gene expression irregularities observed in the gayal. Indeed, we found that many muscle-related genes, located synthetically on pro-chr2 and pro-chr28, exhibited significant changes in expression. This, combined with genome-scale structural variants and expression alterations in genes involved in myofibril composition, may have driven the rapid sarcomere adaptation of gayal to its rugged mountain habitat. Our findings not only suggest that large-scale chromosomal changes can lead to alterations in genome-level expression, thereby promoting both adaptation and speciation, but also illuminate novel avenues for studying the relationship between karyotype evolution and speciation.

Genomic Analyses Reveal Potential Independent Adaptation to High Altitude in Tibetan Chickens.

Much like other indigenous domesticated animals, Tibetan chickens living at high altitudes (2,200-4,100 m) show specific physiological adaptations to the extreme environmental conditions of the Tibetan Plateau, but the genetic bases of these adaptations are not well characterized. Here, we assembled a de novo genome of a Tibetan chicken and resequenced whole genomes of 32 additional chickens, including Tibetan chickens, village chickens, game fowl, and Red Junglefowl, and found that the Tibetan chickens could broadly be placed into two groups. Further analyses revealed that several candidate genes in the calcium-signaling pathway are possibly involved in adaptation to the hypoxia experienced by these chickens, as these genes appear to have experienced directional selection in the two Tibetan chicken populations, suggesting a potential genetic mechanism underlying high altitude adaptation in Tibetan chickens. The candidate selected genes identified in this study, and their variants, may be useful targets for clarifying our understanding of the domestication of chickens in Tibet, and might be useful in current breeding efforts to develop improved breeds for the highlands.

Population variation revealed high-altitude adaptation of Tibetan mastiffs.

With the assistance of their human companions, dogs have dispersed into new environments during the expansion of human civilization. Tibetan Mastiff (TM), a native of the Tibetan Plateau, was derived from the domesticated Chinese native dog and, like Tibetans, has adapted to the extreme environment of high altitude. Here, we genotyped genome-wide single-nucleotide polymorphisms (SNPs) from 32 TMs and compared them with SNPs from 20 Chinese native dogs and 14 gray wolves (Canis lupus). We identified 16 genes with signals of positive selection in the TM, with 12 of these candidate genes associated with functions that have roles in adaptation to high-altitude adaptation, such as EPAS1, SIRT7, PLXNA4, and MAFG that have roles in responses to hypoxia. This study provides important information on the genetic diversity of the TM and potential mechanisms for adaptation to hypoxia.

Mitochondrial DNA variation of Nigerian domestic helmeted guinea fowl.

We analyzed genetic diversity of 215 mitochondrial DNA (mtDNA) D-loop sequences from seven populations of domesticated helmeted guinea fowl (Numida meleagris) in Nigeria and compared that with results of samples collected in Kenya (n = 4) and China (n = 22). In total, 241 sequences were assigned to 22 distinct haplotypes. Haplotype diversity in Nigeria was 0.693 ± 0.022. The network grouped most matrilines into two main haplogroups: A and B. There was an absence of a geographic signal, and two haplotypes dominated across all locations with the exception of the Kebbi population in the northwest of the country; AMOVA also confirmed this observation (FST  = 0.035). The low genetic diversity may be a result of recent domestication, whereas the lack of maternal genetic structure likely suggests the extensive genetic intermixing within the country. Additionally, the differentiation of the Kebbi population may be due to a certain demographic history and/or artificial selection that shaped its haplotype profile. The current data do not permit us to make further conclusions; therefore, more research evidence from genetics and archaeology is still required.

Large-scale mtDNA screening reveals a surprising matrilineal complexity in east Asia and its implications to the peopling of the region.

In order to achieve a thorough coverage of the basal lineages in the Chinese matrilineal pool, we have sequenced the mitochondrial DNA (mtDNA) control region and partial coding region segments of 6,093 mtDNAs sampled from 84 populations across China. By comparing with the available complete mtDNA sequences, 194 of those mtDNAs could not be firmly assigned into the available haplogroups. Completely sequencing 51 representatives selected from these unclassified mtDNAs identified a number of novel lineages, including five novel basal haplogroups that directly emanate from the Eurasian founder nodes (M and N). No matrilineal contribution from the archaic hominid was observed. Subsequent analyses suggested that these newly identified basal lineages likely represent the genetic relics of modern humans initially peopling East Asia instead of being the results of gene flow from the neighboring regions. The observation that most of the newly recognized mtDNA lineages have already differentiated and show the highest genetic diversity in southern China provided additional evidence in support of the Southern Route peopling hypothesis of East Asians. Specifically, the enrichment of most of the basal lineages in southern China and their rather ancient ages in Late Pleistocene further suggested that this region was likely the genetic reservoir of modern humans after they entered East Asia.

Lactase persistence may have an independent origin in Tibetan populations from Tibet, China.

Milk consumption is prevalent in daily diets of Tibetans. To digest the milk sugar lactose, lactase persistence (LP) should be required. However, little is known about the genetic basis of LP in Tibetans. We screened 495 Tibetan individuals for five previously reported single-nucleotide polymorphisms (SNPs): -13907C/G (rs41525747), -13910C/T (rs4988235), -13915T/G (rs41380347), -14010G/C and -22018G/A (rs182549), which are associated with the LP in populations from a vast region surrounding Tibet. The five SNPs were nearly absent in Tibetan populations, suggesting LP likely to have an independent origin in Tibetans rather than to be introduced via gene flow from neighboring populations. We identified three novel SNPs (-13838G/A, -13906T/A and -13908C/T) in Tibetans. In particular, -13838G/A might be functional as it is located in the binding motif for HNF4α that acts as a transcription factor for intestinal gene expression. To investigate the potential association of this variant with LP, further detailed studies are required in the future.

Revisiting the role of the Himalayas in peopling Nepal: insights from mitochondrial genomes.

Himalayas was believed to be a formidably geographical barrier between South and East Asia. The observed high frequency of the East Eurasian paternal lineages in Nepal led some researchers to suggest that these lineages were introduced into Nepal from Tibet directly; however, it is also possible that the East Eurasian genetic components might trace their origins to northeast India where abundant East Eurasian maternal lineages have been detected. To trace the origin of the Nepalese maternal genetic components, especially those of East Eurasian ancestry, and then to better understand the role of the Himalayas in peopling Nepal, we have studied the matenal genetic composition extensively, especially the East Eurasian lineages, in Nepalese and its surrounding populations. Our results revealed the closer affinity between the Nepalese and the Tibetans, specifically, the Nepalese lineages of the East Eurasian ancestry generally are phylogenetically closer with the ones from Tibet, albeit a few mitochondrial DNA haplotypes, likely resulted from recent gene flow, were shared between the Nepalese and northeast Indians. It seems that Tibet was most likely to be the homeland for most of the East Eurasian in the Nepalese. Taking into account the previous observation on Y chromosome, now it is convincing that bearer of the East Eurasian genetic components had entered Nepal across the Himalayas around 6 kilo years ago (kya), a scenario in good agreement with the previous results from linguistics and archeology.

Mitochondrial genome evidence reveals successful Late Paleolithic settlement on the Tibetan Plateau.

Due to its numerous environmental extremes, the Tibetan Plateau--the world's highest plateau--is one of the most challenging areas of modern human settlement. Archaeological evidence dates the earliest settlement on the plateau to the Late Paleolithic, while previous genetic studies have traced the colonization event(s) to no earlier than the Neolithic. To explore whether the genetic continuity on the plateau has an exclusively Neolithic time depth, we studied mitochondrial DNA (mtDNA) genome variation within 6 regional Tibetan populations sampled from Tibet and neighboring areas. Our results confirm that the vast majority of Tibetan matrilineal components can trace their ancestry to Epipaleolithic and Neolithic immigrants from northern China during the mid-Holocene. Significantly, we also identified an infrequent novel haplogroup, M16, that branched off directly from the Eurasian M founder type. Its nearly exclusive distribution in Tibetan populations and ancient age (>21 kya) suggest that M16 may represent the genetic relics of the Late Paleolithic inhabitants on the plateau. This partial genetic continuity between the Paleolithic inhabitants and the contemporary Tibetan populations bridges the results and inferences from archaeology, history, and genetics.

Reticulate evolution: frequent introgressive hybridization among Chinese hares (genus lepus) revealed by analyses of multiple mitochondrial and nuclear DNA loci.

BACKGROUND: Interspecific hybridization may lead to the introgression of genes and genomes across species barriers and contribute to a reticulate evolutionary pattern and thus taxonomic uncertainties. Since several previous studies have demonstrated that introgressive hybridization has occurred among some species within Lepus, therefore it is possible that introgressive hybridization events also occur among Chinese Lepus species and contribute to the current taxonomic confusion. RESULTS: Data from four mtDNA genes, from 116 individuals, and one nuclear gene, from 119 individuals, provides the first evidence of frequent introgression events via historical and recent interspecific hybridizations among six Chinese Lepus species. Remarkably, the mtDNA of L. mandshuricus was completely replaced by mtDNA from L. timidus and L. sinensis. Analysis of the nuclear DNA sequence revealed a high proportion of heterozygous genotypes containing alleles from two divergent clades and that several haplotypes were shared among species, suggesting repeated and recent introgression. Furthermore, results from the present analyses suggest that Chinese hares belong to eight species. CONCLUSION: This study provides a framework for understanding the patterns of speciation and the taxonomy of this clade. The existence of morphological intermediates and atypical mitochondrial gene genealogies resulting from frequent hybridization events likely contribute to the current taxonomic confusion of Chinese hares. The present study also demonstrated that nuclear gene sequence could offer a powerful complementary data set with mtDNA in tracing a complete evolutionary history of recently diverged species.

Tracing the origins of Hakka and Chaoshanese by mitochondrial DNA analysis.

Hakka and Chaoshanese are two unique Han populations residing in southern China but with northern Han (NH) cultural traditions and linguistic influences. Although most of historical records indicate that both populations migrated from northern China in the last two thousand years, no consensus on their origins has been reached so far. To shed more light on the origins of Hakka and Chaoshanese, mitochondrial DNAs (mtDNAs) of 170 Hakka from Meizhou and 102 Chaoshanese from Chaoshan area, Guangdong Province, were analyzed. Our results show that some southern Chinese predominant haplogroups, e.g. B, F, and M7, have relatively high frequencies in both populations. Although median network analyses show that Hakka/Chaoshanese share some haplotypes with NH, interpopulation comparison reveals that both populations show closer affinity with southern Han (SH) populations than with NH. In consideration of previous results from nuclear gene (including Y chromosome) research, it is likely that matrilineal landscapes of both Hakka and Chaoshanese have largely been shaped by the local people during their migration southward and/or later colonization in southern China, and factors such as cultural assimilation, patrilocality, and even sex-bias in the immigrants might have played important roles during the process.

Genomic consequences of population decline in critically endangered pangolins and their demographic histories.

Pangolins are among the most critically endangered animals due to heavy poaching and worldwide trafficking. However, their demographic histories and the genomic consequences of their recent population declines remain unknown. We generated high-quality de novo reference genomes for critically endangered Malayan (Manis javanica, MJ) and Chinese (M. pentadactyla, MP) pangolins and re-sequencing population genomic data from 74 MJs and 23 MPs. We recovered the population identities of illegally traded pangolins and previously unrecognized genetic populations that should be protected as evolutionarily distinct conservation units. Demographic reconstruction suggested environmental changes have resulted in a population size fluctuation of pangolins. Additionally, recent population size declines due to human activities have resulted in an increase in inbreeding and genetic load. Deleterious mutations were enriched in genes related to cancer/diseases and cholesterol homeostasis, which may have increased their susceptibility to diseases and decreased their survival potential to adapt to environmental changes and high-cholesterol diets. This comprehensive study provides not only high-quality pangolin reference genomes, but also valuable information concerning the driving factors of long-term population size fluctuations and the genomic impact of recent population size declines due to human activities, which is essential for pangolin conservation management and global action planning.

Potential dual expansion of domesticated donkeys revealed by worldwide analysis on mitochondrial sequences.

Molecular studies on donkey mitochondrial sequences have clearly defined two distinct maternal lineages involved in domestication. However, domestication histories of these two lineages remain enigmatic. We therefore compared several population characteristics between these two lineages based on global sampling, which included 171 sequences obtained in this study (including Middle Asian, East Asian, and African samples) plus 536 published sequences (including European, Asian, and African samples). The two lineages were clearly separated from each other based on whole mitochondrial genomes and partial non-coding displacement loop (D-loop) sequences, respectively. The Clade I lineage experienced an increase in population size more than 8 000 years ago and shows a complex haplotype network. In contrast, the population size of the Clade II lineage has remained relatively constant, with a simpler haplotype network. Although the distribution of the two lineages was almost equal across the Eurasian mainland, they still presented discernible but complex geographic bias in most parts of Africa, which are known as their domestication sites. Donkeys from sub-Saharan Africa tended to descend from the Clade I lineage, whereas the Clade II lineage was dominant along the East and North coasts of Africa. Furthermore, the migration routes inferred from diversity decay suggested different expansion across China between the two lineages. Altogether, these differences indicated non-simultaneous domestication of the two lineages, which was possibly influenced by the response of pastoralists to the desertification of the Sahara and by the social expansion and trade of ancient humans in Northeast Africa, respectively.

Mitochondrial genomes uncover the maternal history of the Pamir populations.

The Pamirs, among the world's highest mountains in Central Asia, are one of homelands with the most extreme high altitude for several ethnic groups. The settlement history of modern humans on the Pamirs remains still opaque. Herein, we have sequenced the mitochondrial DNA (mtDNA) genomes of 382 individuals belonging to eight populations from the Pamirs and the surrounding lowlands in Central Asia. We construct the Central Asian (including both highlanders and lowlanders) mtDNA haplogroup tree at the highest resolution. All the matrilineal components are assigned into the defined mtDNA haplogroups in East and West Eurasians. No basal lineages that directly emanate from the Eurasian founder macrohaplogroups M, N, and R are found. Our data support the origin of Central Asian being the result of East-West Eurasian admixture. The coalescence ages for more than 93% mtDNA lineages in Central Asians are dated after the last glacial maximum (LGM). The post-LGM and/or later dispersals/admixtures play dominant roles in shaping the maternal gene pool of Central Asians. More importantly, our analyses reveal the mtDNA heterogeneity in the Pamir highlanders, not only between the Turkic Kyrgyz and the Indo-European Tajik groups, but also among three highland Tajiks. No evidence supports positive selection or relaxation of selective constraints in the mtDNAs of highlanders as compared to that of lowlanders. Our results suggest a complex history for the peopling of Pamirs by multiple waves of migrations from various genetic resources during different time scales.

A cryptic mitochondrial DNA link between North European and West African dogs.

Domestic dogs have an ancient origin and a long history in Africa. Nevertheless, the timing and sources of their introduction into Africa remain enigmatic. Herein, we analyse variation in mitochondrial DNA (mtDNA) D-loop sequences from 345 Nigerian and 37 Kenyan village dogs plus 1530 published sequences of dogs from other parts of Africa, Europe and West Asia. All Kenyan dogs can be assigned to one of three haplogroups (matrilines; clades): A, B, and C, while Nigerian dogs can be assigned to one of four haplogroups A, B, C, and D. None of the African dogs exhibits a matrilineal contribution from the African wolf (Canis lupus lupaster). The genetic signal of a recent demographic expansion is detected in Nigerian dogs from West Africa. The analyses of mitochondrial genomes reveal a maternal genetic link between modern West African and North European dogs indicated by sub-haplogroup D1 (but not the entire haplogroup D) coalescing around 12,000 years ago. Incorporating molecular anthropological evidence, we propose that sub-haplogroup D1 in West African dogs could be traced back to the late-glacial dispersals, potentially associated with human hunter-gatherer migration from southwestern Europe.

Origin, genetic diversity, and population structure of Chinese domestic sheep.

To characterize the origin, genetic diversity, and phylogeographic structure of Chinese domestic sheep, we here analyzed a 531-bp fragment of mtDNA control region of 449 Chinese autochthonous sheep from 19 breeds/populations from 13 geographic regions, together with previously reported 44 sequences from Chinese indigenous sheep. Phylogenetic analysis showed that all three previously defined lineages A, B, and C were found in all sampled Chinese sheep populations, except for the absence of lineage C in four populations. Network profiles revealed that the lineages B and C displayed a star-like phylogeny with the founder haplotype in the centre, and that two star-like subclades with two founder haplotypes were identified in lineage A. The pattern of genetic variation in lineage A, together with the divergence time between the two central founder haplotypes suggested that two independent domestication events have occurred in sheep lineage A. Considerable mitochondrial diversity was observed in Chinese sheep. Weak structuring was observed either among Chinese indigenous sheep populations or between Asian and European sheep and this can be attributable to long-term strong gene flow induced by historical human movements. The high levels of intra-population diversity in Chinese sheep and the weak phylogeographic structuring indicated three geographically independent domestication events have occurred and the domestication place was not only confined to the Near East, but also occurred in other regions.

Identifying molecular signatures of hypoxia adaptation from sex chromosomes: A case for Tibetan Mastiff based on analyses of X chromosome.

Genome-wide studies on high-altitude adaptation have received increased attention as a classical case of organismal evolution under extreme environment. However, the current genetic understanding of high-altitude adaptation emanated mainly from autosomal analyses. Only a few earlier genomic studies paid attention to the allosome. In this study, we performed an intensive scan of the X chromosome of public genomic data generated from Tibetan Mastiff (TM) and five other dog populations for indications of high-altitude adaptation. We identified five genes showing signatures of selection on the X chromosome. Notable among these genes was angiomotin (AMOT), which is related to the process of angiogenesis. We sampled additional 11 dog populations (175 individuals in total) at continuous altitudes in China from 300 to 4,000 meters to validate and test the association between the haplotype frequency of AMOT gene and altitude adaptation. The results suggest that AMOT gene may be a notable candidate gene for the adaptation of TM to high-altitude hypoxic conditions. Our study shows that X chromosome deserves consideration in future studies of adaptive evolution.

Comparative population genomics reveals genetic basis underlying body size of domestic chickens.

Body size is the most important economic trait for animal production and breeding. Several hundreds of loci have been reported to be associated with growth trait and body weight in chickens. The loci are mapped to large genomic regions due to the low density and limited number of genetic markers in previous studies. Herein, we employed comparative population genomics to identify genetic basis underlying the small body size of Yuanbao chicken (a famous ornamental chicken) based on 89 whole genomes. The most significant signal was mapped to the BMP10 gene, whose expression was upregulated in the Yuanbao chicken. Overexpression of BMP10 induced a significant decrease in body length by inhibiting angiogenic vessel development in zebrafish. In addition, three other loci on chromosomes 1, 2, and 24 were also identified to be potentially involved in the development of body size. Our results provide a paradigm shift in identification of novel loci controlling body size variation, availing a fast and efficient strategy. These loci, particularly BMP10, add insights into ongoing research of the evolution of body size under artificial selection and have important implications for future chicken breeding.

Genomic analysis of snub-nosed monkeys (Rhinopithecus) identifies genes and processes related to high-altitude adaptation.

The snub-nosed monkey genus Rhinopithecus includes five closely related species distributed across altitudinal gradients from 800 to 4,500 m. Rhinopithecus bieti, Rhinopithecus roxellana, and Rhinopithecus strykeri inhabit high-altitude habitats, whereas Rhinopithecus brelichi and Rhinopithecus avunculus inhabit lowland regions. We report the de novo whole-genome sequence of R. bieti and genomic sequences for the four other species. Eight shared substitutions were found in six genes related to lung function, DNA repair, and angiogenesis in the high-altitude snub-nosed monkeys. Functional assays showed that the high-altitude variant of CDT1 (Ala537Val) renders cells more resistant to UV irradiation, and the high-altitude variants of RNASE4 (Asn89Lys and Thr128Ile) confer enhanced ability to induce endothelial tube formation in vitro. Genomic scans in the R. bieti and R. roxellana populations identified signatures of selection between and within populations at genes involved in functions relevant to high-altitude adaptation. These results provide valuable insights into the adaptation to high altitude in the snub-nosed monkeys.

Ancient inland human dispersals from Myanmar into interior East Asia since the Late Pleistocene.

Given the existence of plenty of river valleys connecting Southeast and East Asia, it is possible that some inland route(s) might have been adopted by the initial settlers to migrate into the interior of East Asia. Here we analyzed mitochondrial DNA (mtDNA) HVS variants of 845 newly collected individuals from 14 Myanmar populations and 5,907 published individuals from 115 populations from Myanmar and its surroundings. Enrichment of basal lineages with the highest genetic diversity in Myanmar suggests that Myanmar was likely one of the differentiation centers of the early modern humans. Intriguingly, some haplogroups were shared merely between Myanmar and southwestern China, hinting certain genetic connection between both regions. Further analyses revealed that such connection was in fact attributed to both recent gene flow and certain ancient dispersals from Myanmar to southwestern China during 25-10 kya, suggesting that, besides the coastal route, the early modern humans also adopted an inland dispersal route to populate the interior of East Asia.

DNA methylation signatures of long intergenic noncoding RNAs in porcine adipose and muscle tissues.

Long intergenic noncoding RNAs (lincRNAs) are one of the major unexplored components of genomes. Here we re-analyzed a published methylated DNA immunoprecipitation sequencing (MeDIP-seq) dataset to characterize the DNA methylation pattern of pig lincRNA genes in adipose and muscle tissues. Our study showed that the methylation level of lincRNA genes was higher than that of mRNA genes, with similar trends observed in comparisons of the promoter, exon or intron regions. Different methylation pattern were observed across the transcription start sites (TSS) of lincRNA and protein-coding genes. Furthermore, an overlap was observed between many lincRNA genes and differentially methylated regions (DMRs) identified among different breeds of pigs, which show different fat contents, sexes and anatomic locations of tissues. We identify a lincRNA gene, linc-sscg3623, that displayed differential methylation levels in backfat between Min and Large White pigs at 60 and 120 days of age. We found that a demethylation process occurred between days 150 and 180 in the Min and Large White pigs, which was followed by remethylation between days 180 and 210. These results contribute to our understanding of the domestication of domestic animals and identify lincRNA genes involved in adipogenesis and muscle development.