Genomic Evidence for the Complex Evolutionary History of Macaques (Genus Macaca).

The genus Macaca is widely distributed, occupies a variety of habitats, shows diverse phenotypic characteristics, and is one of the best-studied genera of nonhuman primates. Here, we reported five re-sequencing Macaca genomes, including one M. cyclopis, one M. fuscata, one M. thibetana, one M. silenus, and one M. sylvanus. Together with published genomes of other macaque species, we combined 20 genome sequences of 10 macaque species to investigate the gene introgression and genetic differences among the species. The network analysis of the SNV-fragment trees indicates a reticular phylogeny of macaque species. Combining the results from various analytical methods, we identified extensive ancient introgression events among macaque species. The multiple introgression signals between different species groups were also observed, such as between fascicularis group species and silenus group species. However, gene flow signals between fascicularis and sinica group were not as strong as those between fascicularis group and silenus group. On the other hand, the unidirect gene flow in M. arctoides probably occurred between the progenitor of M. arctoides and the common ancestor of fascicularis group. Our study also shows that the genetic backgrounds and genetic diversity of different macaques vary dramatically among species, even among populations of the same species. In conclusion, using whole genome sequences and multiple methods, we have studied the evolutionary history of the genus Macaca and provided evidence for extensive introgression among the species.

Single-cell analysis of the cellular landscape of vulvar melanoma provides new insight for immunotherapy administration.

BACKGROUND: Vulvar and vaginal melanoma (VuM & VaM) is a rare gynecologic malignancy with high mortality but low effectiveness to checkpoint immunotherapy compared to cutaneous melanoma. This article aims to elucidate the role of the disordered immune microenvironment in cancer progression in VuM. METHODS: At first, this article applied single-cell RNA sequencing (scRNA-seq) to the VuM obtained from a 68-year-old female patient, and constructed a single-cell atlas of VuM consist of 12,243 single cells. Then this article explores the genomic complexity and core signal channel in VuM microenvironment. RESULTS: This article provides new insights about the pathogenesis of VuM based on single-cell resolution data. It was found that the activation of CD8+ T cell contributed to induce tumor angiogenesis and immune escape, and the activation of the antigen-presenting molecular function participated in melanoma metastasis. CONCLUSION: This article provided new insights into underlining VuM molecular regulation and potential signaling involved in immunotherapy, which would benefit the clinical practice and administration.

Metagenomic analysis of oral and intestinal microbiome of patients during the initial stage of orthodontic treatment.

INTRODUCTION: This prospective study analyzed changes in the oral and intestinal microbiomes in patients before and after fixed orthodontic treatment, elucidating the impacts of fixed orthodontic treatment on patient health and metabolism. METHODS: Metagenomic analysis was conducted on stool, dental plaque, and saliva samples from 10 fixed orthodontic patients. All the samples were sequenced with Illumina NovaSeq 6000 with a paired-end sequencing length of 150 bp. Identification of taxa in metagenomes and functional annotation of genes of the microbiota were performed using the data after quality control. Clinical periodontal parameters, including the gingiva index, plaque index, and pocket probing depth, were examined at each time point in triplicates. Patients also received a table to record their oral hygiene habits of brushing, flossing, and dessert consumption frequency over 1 month. RESULTS: The brushing and flossing times per day of patients were significantly increased after treatment compared with baseline. The number of times a patient ate dessert daily was also fewer after treatment than at baseline. In addition, the plaque index decreased significantly, whereas the pH value of saliva, gingiva index, and pocket probing depth did not change. No significant differences were observed between the participants before and after orthodontic treatment regarding alpha-diversity analysis of the gut, dental plaque, or saliva microbiota. However, on closer analysis, periodontal disease-associated bacteria levels in the oral cavity remain elevated. Alterations in gut microbiota were also observed after orthodontic treatment. CONCLUSIONS: The richness and diversity of the microbiome did not change significantly during the initial stage of fixed orthodontic treatment. However, the levels of periodontal disease-associated bacteria increased.

Metagenomic analysis reveals distinct changes in the gut microbiome of obese Chinese children.

BACKGROUND: The prevalence of obese children in China is increasing, which poses a great challenge to public health. Gut microbes play an important role in human gut health, and changes in gut status are closely related to obesity. However, how gut microbes contribute to obesity in children remains unclear. In our study, we performed shotgun metagenomic sequencing of feces from 23 obese children, 8 overweight children and 22 control children in Chengdu, Sichuan, China. RESULTS: We observed a distinct difference in the gut microbiome of obese children and that of controls. Compared with the controls, bacterial pathogen Campylobacter rectus was significantly more abundant in obese children. In addition, functional annotation of microbial genes revealed that there might be gut inflammation in obese children. The guts of overweight children might belong to the transition state between obese and control children due to a gradient in relative abundance of differentially abundant species. Finally, we compared the gut metagenomes of obese Chinese children and obese Mexican children and found that Trichuris trichiura was significantly more abundant in the guts of obese Mexican children. CONCLUSIONS: Our results contribute to understanding the changes in the species and function of intestinal microbes in obese Chinese children.

Gene expressions between obligate bamboo-eating pandas and non-herbivorous mammals reveal converged specialized bamboo diet adaptation.

BACKGROUND: It is inevitable to change the function or expression of genes during the environmental adaption of species. Both the giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens) belong to Carnivora and have developed similar adaptations to the same dietary switch to bamboos at the morphological and genomic levels. However, the genetic adaptation at the gene expression level is unclear. Therefore, we aimed to examine the gene expression patterns of giant and red panda convergent specialized bamboo-diets. We examined differences in liver and pancreas transcriptomes between the two panda species and other non-herbivorous species. RESULTS: The clustering and PCA plots suggested that the specialized bamboo diet may drive similar expression shifts in these two species of pandas. Therefore, we focused on shared liver and pancreas DEGs (differentially expressed genes) in the giant and red panda relative to other non-herbivorous species. Genetic convergence occurred at multiple levels spanning carbohydrate metabolism, lipid metabolism, and lysine degradation. The shared adaptive convergence DEGs in both organs probably be an evolutionary response to the high carbohydrate, low lipid and lysine bamboo diet. Convergent expression of those nutrient metabolism-related genes in both pandas was an intricate process and subjected to multi-level regulation, including DNA methylation and transcription factor. A large number of lysine degradation and lipid metabolism related genes were hypermethylated in promoter regions in the red panda. Most genes related to carbohydrate metabolism had reduced DNA methylation with increased mRNA expression in giant pandas. Unlike the red panda, the core gene of the lysine degradation pathway (AASS) doesn't exhibit hypermethylation modification in the giant panda, and dual-luciferase reporter assay showed that transcription factor, NR3C1, functions as a transcriptional activator in AASS transcription through the binding to AASS promoter region. CONCLUSIONS: Our results revealed the adaptive expressions and regulations of the metabolism-related genes responding to the unique nutrients in bamboo food and provided data accumulation and research hints for the future revelation of complex mechanism of two pandas underlying convergent adaptation to a specialized bamboo diet.

Multi-omics analysis reveals changes in tryptophan and cholesterol metabolism before and after sexual maturation in captive macaques.

Rhesus macaques (Macaca mulatta, RMs) are widely used in sexual maturation studies due to their high genetic and physiological similarity to humans. However, judging sexual maturity in captive RMs based on blood physiological indicators, female menstruation, and male ejaculation behavior can be inaccurate. Here, we explored changes in RMs before and after sexual maturation based on multi-omics analysis and identified markers for determining sexual maturity. We found that differentially expressed microbiota, metabolites, and genes before and after sexual maturation showed many potential correlations. Specifically, genes involved in spermatogenesis (TSSK2, HSP90AA1, SOX5, SPAG16, and SPATC1) were up-regulated in male macaques, and significant changes in gene (CD36), metabolites (cholesterol, 7-ketolithocholic acid, and 12-ketolithocholic acid), and microbiota (Lactobacillus) related to cholesterol metabolism were also found, suggesting the sexually mature males have stronger sperm fertility and cholesterol metabolism compared to sexually immature males. In female macaques, most differences before and after sexual maturity were related to tryptophan metabolism, including changes in IDO1, IDO2, IFNGR2, IL1Β, IL10, L-tryptophan, kynurenic acid (KA), indole-3-acetic acid (IAA), indoleacetaldehyde, and Bifidobacteria, indicating that sexually mature females exhibit stronger neuromodulation and intestinal immunity than sexually immature females. Cholesterol metabolism-related changes (CD36, 7-ketolithocholic acid, 12-ketolithocholic acid) were also observed in female and male macaques. Exploring differences before and after sexual maturation through multi-omics, we identified potential biomarkers of sexual maturity in RMs, including Lactobacillus (for males) and Bifidobacterium (for females) valuable for RM breeding and sexual maturation research.

Pyfastx: a robust Python package for fast random access to sequences from plain and gzipped FASTA/Q files.

FASTA and FASTQ are the most widely used biological data formats that have become the de facto standard to exchange sequence data between bioinformatics tools. With the avalanche of next-generation sequencing data, the amount of sequence data being deposited and accessed in FASTA/Q formats is increasing dramatically. However, the existing tools have very low efficiency at random retrieval of subsequences due to the requirement of loading the entire index into memory. In addition, most existing tools have no capability to build index for large FASTA/Q files because of the limited memory. Furthermore, the tools do not provide support to randomly accessing sequences from FASTA/Q files compressed by gzip, which is extensively adopted by most public databases to compress data for saving storage. In this study, we developed pyfastx as a versatile Python package with commonly used command-line tools to overcome the above limitations. Compared to other tools, pyfastx yielded the highest performance in terms of building index and random access to sequences, particularly when dealing with large FASTA/Q files with hundreds of millions of sequences. A key advantage of pyfastx over other tools is that it offers an efficient way to randomly extract subsequences directly from gzip compressed FASTA/Q files without needing to uncompress beforehand. Pyfastx can easily be installed from PyPI (https://pypi.org/project/pyfastx) and the source code is freely available at https://github.com/lmdu/pyfastx.

Transcriptomics and metagenomics of common cutworm (Spodoptera litura) and fall armyworm (Spodoptera frugiperda) demonstrate differences in detoxification and development.

BACKGROUND: Spodoptera litura is an important polyphagous pest that causes significant damage to the agricultural sector. We performed RNA-seq of 15 S. litura individuals from larval (fifth and sixth instar larvae), chrysalis, and adult developmental stages. We also compared the S. litura transcriptome data with Spodoptera frugiperda across the same developmental stages, which was sequenced in our previous study. RESULTS: A total of 101,885 differentially expressed transcripts (DETs) were identified in S. litura. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that S. litura may undergo active xenobiotic and detoxifying metabolism during its larval and adult stages, which may explain difficulties with current population control measures. We also found that DETs of single-copy orthologous genes between S. litura and S. frugiperda were involved in basic metabolism and development. However, energy and metabolic processes genes had a higher expression in S. litura, whereas nervous and olfactory function genes had a higher expression in S. frugiperda. Metagenomics analysis in larval S. litura and S. frugiperda revealed that microbiota participate in the detoxification and metabolism processes, but the relative abundance of detoxification-related microbiota was more abundant in S. frugiperda. Transcriptome results also confirmed the detoxification-related pathway of S. frugiperda was more abundant than in S. litura. CONCLUSIONS: Significant changes at transcriptional level were identified during the different development stages of S. litura. Importantly, we also identified detoxification associated genes and gut microbiota between S. litura and S. frugiperda at different developmental stages, which will be valuable in revealing possible mechanisms of detoxification and development in these two lepidopterans.

Sex-specific gene expression in the blood of four primates.

Blood is an important non-reproductive tissue, but little is known about the sex-specific gene expressions in the blood. Therefore, we investigated sex-specific gene expression differences in the blood tissues of four primates, rhesus macaques (Macaca mulatta), Tibetan macaques (M. thibetana), yellow baboons (Papio cynocephalus), and humans. We identified seven sex-specific differentially expressed genes (SDEGs) in each non-human primate and 31 SDEGs in humans. The four primates had only one common SDEG, MAP7D2. In humans, immune-related SDEGs were identified as up-regulated, but also down-regulated in females. We also found that most of the X-Y gene pairs had similar expression levels between species, except pair EIF1AY/EIF1AX. The expression level of X-Y gene pairs of rhesus and Tibetan macaques showed no significant differential expression levels, while humans had six significant XY-biased and three XX-biased X-Y gene pairs. Our observed sex differences in blood should increase understanding of sex differences in primate blood tissue.

Age-related gene expression and DNA methylation changes in rhesus macaque.

Aging is a very complicated biological process that can change gene expressions. The Chinese rhesus macaque (Macaca mulatta lasiota; CR) is closely related to humans. We explored gene expression with increasing age and DNA methylation changes in young and old CRs. Results showed blood transcriptome and DNA methylome significantly changed from young to old CRs. The age-associated differentially expressed genes (DEGs) and differentially methylated regions (DMRs) were associated with age-related biological features, such as immunity, blood coagulation, and biosynthetic process. The measurements of coagulation indicators confirmed old CRs had shorter coagulation time than young CRs, and the activities of coagulation factor II (FII) and factor VIII (FVIII) were enhanced in old CRs. Humans and CRs exhibited the same enhanced blood coagulation with age phenotype. Our study found aging is a critical factor affecting gene expression in CRs, and also provided new insights into the blood coagulation changes in non-human primates.

Characterization of the primate TRIM gene family reveals the recent evolution in primates.

The tripartite motif (TRIM) gene family encodes diverse distinct proteins that play important roles in many biological processes. However, the molecular evolution and phylogenetic relationships of TRIM genes in primates are still elusive. We performed a genomic approach to identify and characterize TRIM genes in human and other six primate genomes. In total, 537 putative functional TRIM genes were identified and TRIM members varied among primates. A neighbor joining (NJ) tree based on the protein sequences of 82 human TRIM genes indicates seven TRIM groups, which is consistent with the results based on the architectural motifs. Many TRIM gene duplication events were identified, indicating a recent expansion of TRIM family in primate lineages. Interestingly, the chimpanzee genome shows the greatest TRIM gene expansion among the primates; however, its congeneric species, bonobo, has the least number of TRIM genes and no duplication event. Moreover, we identified a ~ 200 kb deletion on chromosome 11 of bonobos that results in a loss of cluster3 TRIM genes. The loss of TRIM genes might have occurred within the last 2 mys. Analysis of positive selection recovered 9 previously reported and 21 newly identified positively selected TRIM genes. In particular, most positive selected sites are located in the B30.2 domains. Our results have provided new insight into the evolution of primate TRIM genes and will broaden our understanding on the functions of the TRIM family.

Genomic evidence for the Old divergence of Southern European wolf populations.

The grey wolf (Canis lupus) is one of the most widely distributed mammals in which a variety of distinct populations have been described. However, given their currently fragmented distribution and recent history of human-induced population decline, little is known about the events that led to their differentiation. Based on the analysis of whole canid genomes, we examined the divergence times between Southern European wolf populations and their ancient demographic history. We found that all present-day Eurasian wolves share a common ancestor ca 36 000 years ago, supporting the hypothesis that all extant wolves derive from a single population that subsequently expanded after the Last Glacial Maximum. We also estimated that the currently isolated European populations of the Iberian Peninsula, Italy and the Dinarics-Balkans diverged very closely in time, ca 10 500 years ago, and maintained negligible gene flow ever since. This indicates that the current genetic and morphological distinctiveness of Iberian and Italian wolves can be attributed to their isolation dating back to the end of the Pleistocene, predating the recent human-induced extinction of wolves in Central Europe by several millennia.

Comprehensive analysis of lncRNA and mRNA expression changes in Tibetan chicken lung tissue between three developmental stages.

The Tibetan chicken is a native Chinese breed that lives at high elevations and has adapted to the extreme environmental conditions of the Tibetan Plateau. However, its hypoxic adaptation at the gene expression level is unclear. Here, we sequenced nine lung transcriptomes of the Tibetan chicken at three developmental stages (5 and 42 weeks and 4.5 years). A total of 1.02 billion clean reads were obtained. We identified 16 012 mRNAs and 6898 lncRNAs. The expression of mRNA showed that nine samples were significantly divided into three clusters, with higher correlation and closer relationship between the 5 and 42 week groups. We identified 399 differentially expressed genes (DEGs) between the 5 and 42 week groups, 3532 DEGs between the 5 week and 4.5 year groups, and 3909 DEGs between the 42 week and 4.5 year groups. The up-regulated DEGs in the 5 week group, compared with 42 week and 4.5 year groups, were enriched in GO terms associated with growth and development, whereas the up-regulated DEGs in the 4.5 year group were mainly enriched in many metabolic-related categories. Moreover, the enrichment results with up-regulated DEGs in the 5 and/or 42 week groups, compared with the 4.5 year group, were associated with hypoxic adaptation, such as oxygen transport, oxygen binding and oxygen carrier activity, and calcium signaling pathway. In addition, we identified 978 high-correlation lncRNA and protein-coding gene pairs, and 524 significant neighboring protein-coding genes were also DEGs. Our results provide new insights into gene expression of lung tissue in Tibetan chickens during the aging process.

The first draft genome of Lophophorus: A step forward for Phasianidae genomic diversity and conservation.

The monal genus (Lophophorus) is a branch of Phasianidae and its species inhabit the high-altitude mountains of the Qinghai-Tibet Plateau. The Chinese monal, L. lhuysii, is a threatened endemic bird of China that possesses high-altitude adaptability, diversity of plumage color and potentially low reproductive life history. This is the first study to describe the monal genome using next generation sequencing technology. The Chinese monal genome size is 1.01 Gb, with 16,940 protein-coding genes. Gene annotation yielded 100.93 Mb (9.97%) repeat elements, 785 ncRNA, 5,465,549 bp (0.54%) SSR and 15,550 (92%) genes in public databases. Compared to other birds and mammals, the genome evolution analysis showed numerous expanded gene families and positive selected genes involved in high-altitude adaptation, especially related to the adaptation of low temperature and hypoxia. Consequently, this gene data can be used to investigate the molecular evolution of high-altitude adaptation in future bird research. Our first published genome of the genus Lophophorus will be integral for the study of monal population genetic diversity and conservation, genomic evolution and Galliformes species differentiation in the Qinghai-Tibetan Plateau.

Worldwide patterns of genomic variation and admixture in gray wolves.

The gray wolf (Canis lupus) is a widely distributed top predator and ancestor of the domestic dog. To address questions about wolf relationships to each other and dogs, we assembled and analyzed a data set of 34 canine genomes. The divergence between New and Old World wolves is the earliest branching event and is followed by the divergence of Old World wolves and dogs, confirming that the dog was domesticated in the Old World. However, no single wolf population is more closely related to dogs, supporting the hypothesis that dogs were derived from an extinct wolf population. All extant wolves have a surprisingly recent common ancestry and experienced a dramatic population decline beginning at least ∼30 thousand years ago (kya). We suggest this crisis was related to the colonization of Eurasia by modern human hunter-gatherers, who competed with wolves for limited prey but also domesticated them, leading to a compensatory population expansion of dogs. We found extensive admixture between dogs and wolves, with up to 25% of Eurasian wolf genomes showing signs of dog ancestry. Dogs have influenced the recent history of wolves through admixture and vice versa, potentially enhancing adaptation. Simple scenarios of dog domestication are confounded by admixture, and studies that do not take admixture into account with specific demographic models are problematic.

Genome-wide analysis reveals the genomic features of the turkey vulture (Cathartes aura) as a scavenger.

Cathartidae is a small family of large-bodied carrion-feeding birds, of which the turkey vulture (Cathartes aura, Cathartidae) is the most widespread distributed. To investigate the chemoreception system, detoxification system, and immune system in the turkey vulture, we compared its genome to 14 other avian genomes. Comparative genomics demonstrated the expansion in the chemoreception system, especially the olfactory receptors, while the genes in the detoxification system of the turkey vulture did not show apparent expansion. We identified five positively selected genes associated with the immune system in the turkey vulture, which was likely to strengthen the immune defense against pathogenic invasion. Functional enrichment analysis indicated that many positively selected genes were involved in the regulation of immune system processes, implying important reorganization of the immune system in the turkey vulture. The turkey vulture-specific missense mutations were found in one positively selected gene (BCL6), and all the missense mutations were classified as deleterious by PolyPhen-2, possibly contributing to immune adaptation to the carrion feeding. Furthermore, we identified four turkey vulture-specific missense mutations in three ß-defensin genes of the turkey vulture, which was an indispensable part in the innate immunity (a natural barrier against invasive microbes including bacteria, fungi, and viruses). Our genomic analyses in the turkey vulture provided insights into the genetic signatures of the adaptation to the carrion feeding.

Demographically-Based Evaluation of Genomic Regions under Selection in Domestic Dogs.

Controlling for background demographic effects is important for accurately identifying loci that have recently undergone positive selection. To date, the effects of demography have not yet been explicitly considered when identifying loci under selection during dog domestication. To investigate positive selection on the dog lineage early in the domestication, we examined patterns of polymorphism in six canid genomes that were previously used to infer a demographic model of dog domestication. Using an inferred demographic model, we computed false discovery rates (FDR) and identified 349 outlier regions consistent with positive selection at a low FDR. The signals in the top 100 regions were frequently centered on candidate genes related to brain function and behavior, including LHFPL3, CADM2, GRIK3, SH3GL2, MBP, PDE7B, NTAN1, and GLRA1. These regions contained significant enrichments in behavioral ontology categories. The 3rd top hit, CCRN4L, plays a major role in lipid metabolism, that is supported by additional metabolism related candidates revealed in our scan, including SCP2D1 and PDXC1. Comparing our method to an empirical outlier approach that does not directly account for demography, we found only modest overlaps between the two methods, with 60% of empirical outliers having no overlap with our demography-based outlier detection approach. Demography-aware approaches have lower-rates of false discovery. Our top candidates for selection, in addition to expanding the set of neurobehavioral candidate genes, include genes related to lipid metabolism, suggesting a dietary target of selection that was important during the period when proto-dogs hunted and fed alongside hunter-gatherers.

Ancient hybridization and admixture in macaques (genus Macaca) inferred from whole genome sequences.

The evolutionary history of the stump-tailed macaque (Macaca arctoides) and its genetic relationship to other macaques is a subject of continuing controversy. Here, we have reported the first genome sequences of two stump-tailed macaques and one Assamese macaque (M. assamensis). Additionally, we have investigated the genetic diversity between macaque species and analyzed ancient hybridization events. Genome-wide analyses demonstrated that the stump-tailed macaque is more closely related to sinica species than to fascicularis/mulatta species. This topology contradicts the mitochondrial sequence-based phylogeny that places the stump-tailed macaque into the fascicularis/mulatta group. However, our results further show that stump-tailed macaques have genetic backgrounds distinct from sinica species, and present evidence of gene flows with rhesus macaques. We suggest that an ancient introgression occurred after stump-tailed macaques diverged from sinica species. The distinct gene flow between proto-arctoides and proto-mulatta resulted in the transfer of rhesus macaque-type mitochondria into proto-arctoides. The rhesus macaque-type mitochondria remained in populations because of genetic drift during the bottleneck. The PSMC results and morphological and geographic evidence are consistent with the mitochondria capture pattern in the stump-tailed macaque. The molecular clock estimates suggest that the mitochondrial transference into stump-tailed macaques occurred 0.4-1.4 million years ago. Furthermore, we detected extensive admixtures between different macaque species, indicating that gene flow has played an important role in the evolutionary history of the genus Macaca.

Comparative genome-wide survey of single nucleotide variation uncovers the genetic diversity and potential biomedical applications among six Macaca species.

Macaca is of great importance in evolutionary and biomedical research. Aiming at elucidating genetic diversity patterns and potential biomedical applications of macaques, we characterized single nucleotide variations (SNVs) of six Macaca species based on the reference genome of Macaca mulatta. Using eight whole-genome sequences, representing the most comprehensive genomic SNV study in Macaca to date, we focused on discovery and comparison of nonsynonymous SNVs (nsSNVs) with bioinformatic tools. We observed that SNV distribution patterns were generally congruent among the eight individuals. Outlier tests of nsSNV distribution patterns detected 319 bins with significantly distinct genetic divergence among macaques, including differences in genes associated with taste transduction, homologous recombination, and fat and protein digestion. Genes with specific nsSNVs in various macaques were differentially enriched for metabolism pathways, such as glycolysis, protein digestion and absorption. On average, 24.95% and 11.67% specific nsSNVs were putatively deleterious according to PolyPhen2 and SIFT4G, respectively, among which the shared deleterious SNVs were located in 564⁻1981 genes. These genes displayed enrichment signals in the 'obesity-related traits' disease category for all surveyed macaques, confirming that they were suitable models for obesity related studies. Additional enriched disease categories were observed in some macaques, exhibiting promising potential for biomedical application. Positively selected genes identified by PAML in most tested Macaca species played roles in immune and nervous system, growth and development, and fat metabolism. We propose that metabolism and body size play important roles in the evolutionary adaptation of macaques.

Similar genomic proportions of copy number variation within gray wolves and modern dog breeds inferred from whole genome sequencing.

BACKGROUND: Whole genome re-sequencing data from dogs and wolves are now commonly used to study how natural and artificial selection have shaped the patterns of genetic diversity. Single nucleotide polymorphisms, microsatellites and variants in mitochondrial DNA have been interrogated for links to specific phenotypes or signals of domestication. However, copy number variation (CNV), despite its increasingly recognized importance as a contributor to phenotypic diversity, has not been extensively explored in canids. RESULTS: Here, we develop a new accurate probabilistic framework to create fine-scale genomic maps of segmental duplications (SDs), compare patterns of CNV across groups and investigate their role in the evolution of the domestic dog by using information from 34 canine genomes. Our analyses show that duplicated regions are enriched in genes and hence likely possess functional importance. We identify 86 loci with large CNV differences between dogs and wolves, enriched in genes responsible for sensory perception, immune response, metabolic processes, etc. In striking contrast to the observed loss of nucleotide diversity in domestic dogs following the population bottlenecks that occurred during domestication and breed creation, we find a similar proportion of CNV loci in dogs and wolves, suggesting that other dynamics are acting to particularly select for CNVs with potentially functional impacts. CONCLUSIONS: This work is the first comparison of genome wide CNV patterns in domestic and wild canids using whole-genome sequencing data and our findings contribute to study the impact of novel kinds of genetic changes on the evolution of the domestic dog.