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 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.

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.

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.

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.

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.

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.

Genetic convergence in the adaptation of dogs and humans to the high-altitude environment of the tibetan plateau.

The high-altitude hypoxic environment represents one of the most extreme challenges for mammals. Previous studies of humans on the Tibetan plateau and in the Andes Mountains have identified statistical signatures of selection in different sets of loci. Here, we first measured the hemoglobin levels in village dogs from Tibet and those from Chinese lowlands. We found that the hemoglobin levels are very similar between the two groups, suggesting that Tibetan dogs might share similar adaptive strategies as the Tibetan people. Through a whole-genome sequencing approach, we have identified EPAS1 and HBB as candidate genes for the hypoxic adaptation on the Tibetan plateau. The population genetic analysis shows a significant convergence between humans and dogs in Tibet. The similarities in the sets of loci that exhibit putative signatures of selection and the hemoglobin levels between humans and dogs of the same environment, but not between human populations in different regions, suggests an extraordinary landscape of convergent evolution between human beings and their best friend on the Tibetan plateau.

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.

Retrieving Y chromosomal haplogroup trees using GWAS data.

Phylogenetically informative Y chromosomal single-nucleotide polymorphisms (Y-SNPs) integrated in DNA chips have not been sufficiently explored in most genome-wide association studies (GWAS). Herein, we introduce a pipeline to retrieve Y-SNP data. We introduce the software YTool (http://mitotool.org/ytool/) to handle conversion, filtering, and annotation of the data. Genome-wide SNP data from populations in Myanmar are used to construct a haplogroup tree for 117 Y chromosomes based on 369 high-confidence Y-SNPs. Parallel genotyping and published resequencing data of Y chromosomes confirm the validity of our pipeline. We apply this strategy to the CEU HapMap data set and construct a haplogroup tree with 107 Y-SNPs from 39 individuals. The retrieved Y-SNPs can discern the parental genetic structure of populations. Given the massive quantity of data from GWAS, this method facilitates future investigations of Y chromosome diversity.

Origin of Chinese goldfish and sequential loss of genetic diversity accompanies new breeds.

BACKGROUND: Goldfish, Carassius auratus, have experienced strong anthropogenic selection during their evolutionary history, generating a tremendous extent of morphological variation relative to that in native Carassius. To locate the geographic origin of goldfish, we analyzed nucleotide sequences from part of the control region (CR) and the entire cytochrome b (Cytb) mitochondrial DNA genes for 234 goldfish and a large series of native specimens. Four important morphological characteristics used in goldfish taxonomy-body shape, dorsal fin, eye shape, and tailfin-were selected for hypothesis-testing to identify those that better correspond to evolutionary history. PRINCIPAL FINDING: Haplotypes of goldfish rooted in two sublineages (C5 and C6), which contained the haplotypes of native C. a. auratus from southern China. Values of F(ST) and N(m) revealed a close relationship between goldfish and native C. a. auratus from the lower Yangtze River. An extraordinary, stepwise loss of genetic diversity was detected from native fish to goldfish and from Grass-goldfish relative to other breeds. Significantly negative results for the tests of Tajima's D and Fu and Li's D* and F* were identified in goldfish, including the Grass breed. The results identified eye-shape as being the least informative character for grouping goldfish with respect to their evolutionary history. Fisher's exact test identified matrilineal constraints on domestication. CONCLUSIONS: Chinese goldfish have a matrilineal origin from native southern Chinese C. a. auratus, especially the lineages from the lower Yangtze River. Anthropogenic selection of the native Carassius eliminated aesthetically unappealing goldfish and this action appeared to be responsible for the stepwise decrease in genetic diversity of domesticated goldfish, a process similar to that reported for the domestication of pigs, rice, and maize. The three-breed taxonomy--Grass-goldfish, Egg-goldfish, and Wen-goldfish--better reflected the history of domestication.

Ecological genetics of Chinese rhesus macaque in response to mountain building: all things are not equal.

BACKGROUND: Pliocene uplifting of the Qinghai-Tibetan Plateau (QTP) and Quaternary glaciation may have impacted the Asian biota more than any other events. Little is documented with respect to how the geological and climatological events influenced speciation as well as spatial and genetic structuring, especially in vertebrate endotherms. Macaca mulatta is the most widely distributed non-human primate. It may be the most suitable model to test hypotheses regarding the genetic consequences of orogenesis on an endotherm. METHODOLOGY AND PRINCIPAL FINDINGS: Using a large dataset of maternally inherited mitochondrial DNA gene sequences and nuclear microsatellite DNA data, we discovered two maternal super-haplogroups exist, one in western China and the other in eastern China. M. mulatta formed around 2.31 Ma (1.51-3.15, 95%), and divergence of the two major matrilines was estimated at 1.15 Ma (0.78-1.55, 95%). The western super-haplogroup exhibits significant geographic structure. In contrast, the eastern super-haplogroup has far greater haplotypic variability with little structure based on analyses of six variable microsatellite loci using Structure and Geneland. Analysis using Migrate detected greater gene flow from WEST to EAST than vice versa. We did not detect signals of bottlenecking in most populations. CONCLUSIONS: Analyses of the nuclear and mitochondrial datasets obtained large differences in genetic patterns for M. mulatta. The difference likely reflects inheritance mechanisms of the maternally inherited mtDNA genome versus nuclear biparentally inherited STRs and male-mediated gene flow. Dramatic environmental changes may be responsible for shaping the matrilineal history of macaques. The timing of events, the formation of M. mulatta, and the divergence of the super-haplogroups, corresponds to both the uplifting of the QTP and Quaternary climatic oscillations. Orogenesis likely drove divergence of western populations in China, and Pleistocene glaciations are likely responsible for genetic structuring in the eastern super-haplogroup via geographic isolation and secondary contact.

Mitochondrial DNA haplogroup Y is associated to Leigh syndrome in Chinese population.

Although Leigh syndrome (LS) is a well characterized clinical mitochondrial disorder; the exact mutation is not found in all cases and it is not clear whether matrilineal background has contributed to this disease. To address this issue, we extensively studied and compared the haplogroup composition of a sample of 171 Chinese LS patients with that of 1597 controls. Our results show that haplogroup Y may increase the risk of LS in Chinese by 2.867 fold (95% CI=1.135-7.240, P=0.020). Haplogroup B5 has also this trend (1.737 fold, 95% CI=0.961-3.139), but with a borderline P-value (P=0.065). Both haplogroups belong to macro-haplogroup N and share a common reverse mutation on nucleotide position 10398 (A10398G). In fact, the combined haplogroup N with 10398G is also associated with an increased risk for LS (OR=1.882, 95% CI=1.134-3.124, P=0.013).

Patrilineal perspective on the Austronesian diffusion in Mainland Southeast Asia.

The Cham people are the major Austronesian speakers of Mainland Southeast Asia (MSEA) and the reconstruction of the Cham population history can provide insights into their diffusion. In this study, we analyzed non-recombining region of the Y chromosome markers of 177 unrelated males from four populations in MSEA, including 59 Cham, 76 Kinh, 25 Lao, and 17 Thai individuals. Incorporating published data from mitochondrial DNA (mtDNA), our results indicated that, in general, the Chams are an indigenous Southeast Asian population. The origin of the Cham people involves the genetic admixture of the Austronesian immigrants from Island Southeast Asia (ISEA) with the local populations in MSEA. Discordance between the overall patterns of Y chromosome and mtDNA in the Chams is evidenced by the presence of some Y chromosome lineages that prevail in South Asians. Our results suggest that male-mediated dispersals via the spread of religions and business trade might play an important role in shaping the patrilineal gene pool of the Cham people.

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.