Whole-genome sequencing of Berkshire (European native pig) provides insights into its origin and domestication

The role of G protein-coupled receptors and their ligands in animal domestication

The domestication of plants and animals has resulted in one of the most significant cultural and socio-economical transitions in human history. Domestication of animals, including human-supervised reproduction, largely uncoupled particular animal species from their natural, evolutionary history driven by environmental and ecological factors. The primary motivations for domesticating animals were, and still are, producing food and materials (e.g. meat, eggs, honey or milk products, wool, leather products, jewelry and medication products) to support plowing in agriculture or in transportation (e.g. horse, cattle, camel and llama) and to facilitate human activities (for hunting, rescuing, therapeutic aid, guarding behavior and protecting or just as a companion). In recent years, decoded genetic information from more than 40 domesticated animal species have become available; these studies have identified genes and mutations associated with specific physiological and behavioral traits contributing to the complex genetic background of animal domestication. These breeding-altered genomes provide insights into the regulation of different physiological areas, including information on links between e.g. endocrinology and behavior, with important pathophysiological implications (e.g. for obesity and cancer), extending the interest in domestication well beyond the field. Several genes that have undergone selection during domestication and breeding encode specific G protein-coupled receptors, a class of membrane-spanning receptors involved in the regulation of a number of overarching functions such as reproduction, development, body homeostasis, metabolism, stress responses, cognition, learning and memory. Here we summarize the available literature on variations in G protein-coupled receptors and their ligands and how these have contributed to animal domestication.

A chromosome-level genome of Chenghua pig provides new insights into the domestication and local adaptation of pigs

As a country with abundant genetic resources of pigs, the domestication history of pigs in China and the adaptive evolution of Chinese pig breeds at different latitudes have rarely been elucidated at the genome-wide level. To fill this gap, we first assembled a high-quality chromosome-level genome of the Chenghua pig and used it as a benchmark to analyse the genomes of 272 samples from three genera of three continents. The divergence of the three species belonging to three genera, Phacochoerus africanus, Potamochoerus porcus, and Sus scrofa, was assessed. The introgression of pig breeds redefined that the migration routes were basically from southern China to central and southwestern China, then spread to eastern China, arrived in northern China, and finally reached Europe. The domestication of pigs in China occurred ∼12,000 years ago, earlier than the available Chinese archaeological domestication evidence. In addition, FBN1 and NR6A1 were identified in our study as candidate genes related to extreme skin thickness differences in Eurasian pig breeds and adaptive evolution at different latitudes in Chinese pig breeds, respectively. Our study provides a new resource for the pig genomic pool and refines our understanding of pig genetic diversity, domestication, migration, and adaptive evolution at different latitudes.

Haplotype-resolved 3D chromatin architecture of the hybrid pig

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Integrative Analysis of Transcriptomic and Lipidomic Profiles Reveals a Differential Subcutaneous Adipose Tissue Mechanism among Ningxiang Pig and Berkshires, and Their Offspring

Adipose tissue composition contributes greatly to the quality and nutritional value of meat. Transcriptomic and lipidomic techniques were used to investigate the molecular mechanisms of the differences in fat deposition in Ningxiang pigs, Berkshires and F1 offspring. Transcriptomic analysis identified 680, 592, and 380 DEGs in comparisons of Ningxiang pigs vs. Berkshires, Berkshires vs. F1 offspring, and Ningxiang pigs vs. F1 offspring. The lipidomic analysis screened 423, 252, and 50 SCLs in comparisons of Ningxiang pigs vs. Berkshires, Berkshires vs. F1 offspring, and Ningxiang pigs vs. F1 offspring. Lycine, serine, and the threonine metabolism pathway, fatty acid biosynthesis and metabolism-related pathways were significantly enriched in comparisons of Berkshires vs. Ningxiang pigs and Berkshires vs. F1 offspring. The DEGs (PHGDH, LOC110256000) and the SCLs (phosphatidylserines) may have a great impact on the glycine, serine, and the threonine metabolism pathway. Moreover, the DEGs (FASN, ACACA, CBR4, SCD, ELOV6, HACD2, CYP3A46, CYP2B22, GPX1, and GPX3) and the SCLs (palmitoleic acid, linoleic acid, arachidonic acid, and icosadienoic acid) play important roles in the fatty acid biosynthesis and metabolism of fatty acids. Thus, the difference in fat deposition among Ningxiang pig, Berkshires, and F1 offspring may be caused by differences in the expression patterns of key genes in multiple enriched KEGG pathways. This research revealed multiple lipids that are potentially available biological indicators and screened key genes that are potential targets for molecular design breeding. The research also explored the molecular mechanisms of the difference in fat deposition among Ningxiang pig, Berkshires, and F1 pigs, and provided an insight into selection for backfat thickness and the fat composition of adipose tissue for future breeding strategies.

Backcrossing to different parents produced two distinct hybrid species

Repeated homoploid hybrid speciation (HHS) events with the same parental species have rarely been reported. In this study, we used population transcriptome data to test paraphyly and HHS events in the conifer Picea brachytyla. Our analyses revealed non-sister relationships for two lineages of P. brachytyla, with the southern lineage being placed within the re-circumscribed P. likiangensis species complex (PLSC) and P. brachytyla sensu stricto (s.s.) consisted solely of the northern lineage, forming a distinct clade that is paratactic to both the PLSC and P. wilsonii. Our phylogenetic and coalescent analyses suggested that P. brachytyla s.s. arose from HHS between the ancestor of the PLSC before its diversification and P. wilsonii through an intermediate hybrid lineage at an early stage and backcrossing to the ancestral PLSC. Additionally, P. purpurea shares the same parents and an extinct lineage with P. brachytyla s.s. but backcrossing to the other parent, P. wilsonii at a later stage. We reveal the first case that backcrossing to different parents of the same extinct hybrid lineage produced two different hybrid species. Our results highlight the existence of more reticulate evolution during species diversification in the spruce genus and more complex homoploid hybrid events than previously identified.

Environment and genotype predict the genomic nature of domestication of salmonids as revealed by gene expression

Billions of salmonids are produced annually by artificial reproduction for harvest and conservation. Morphologically, behaviourally and physiologically these fish differ from wild-born fish, including in ways consistent with domestication. Unlike most studied domesticates, which diverged from wild ancestors millennia ago, salmonids offer a tractable model for early-stage domestication. Here, we review a fundamental mechanism for domestication-driven differences in early-stage domestication, differentially expressed genes (DEGs), in salmonids. We found 34 publications examining DEGs under domestication driven by environment and genotype, covering six species, over a range of life-history stages and tissues. Three trends emerged. First, domesticated genotypes have increased expression of growth hormone and related metabolic genes, with differences magnified under artificial environments with increased food. Regulatory consequences of these DEGs potentially drive overall DEG patterns. Second, immune genes are often DEGs under domestication and not simply owing to release from growth-immune trade-offs under increased food. Third, domesticated genotypes exhibit reduced gene expression plasticity, with plasticity further reduced in low-complexity environments typical of production systems. Recommendations for experimental design improvements, coupled with tissue-specific expression and emerging analytical approaches for DEGs present tractable avenues to understand the evolution of domestication in salmonids and other species.

Population Structure and Selection Signatures Underlying Domestication Inferred from Genome-Wide Copy Number Variations in Chinese Indigenous Pigs

Single nucleotide polymorphism was widely used to perform genetic and evolution research in pigs. However, little is known about the effect of copy number variation (CNV) on characteristics in pigs. This study performed a genome-wide comparison of CNVs between Wannan black pigs (WBP) and Asian wild boars (AWB), using whole genome resequencing data. By using Manta, we detected in total 28,720 CNVs that covered approximately 1.98% of the pig genome length. We identified 288 selected CNVs (top 1%) by performing Fst statistics. Functional enrichment analyses for genes located in selected CNVs were found to be muscle related (NDN, TMOD4, SFRP1, and SMYD3), reproduction related (GJA1, CYP26B1, WNT5A, SRD5A2, PTPN11, SPEF2, and CCNB1), residual feed intake (RFI) related (MAP3K5), and ear size related (WIF1). This study provides essential information on selected CNVs in Wannan black pigs for further research on the genetic basis of the complex phenotypic and provides essential information for direction in the protection and utilization of Wannan black pig.

Lost pigs of Angola: Whole genome sequencing reveals unique regions of selection with emphasis on metabolism and feed efficiency

Angola, in the western coast of Africa, has been through dramatic social events that have led to the near-disappearance of native swine populations, and the recent introduction of European exotic breeds has also contributed to the erosion of this native swine repertoire. In an effort to investigate the genetic basis of native pigs in Angola (ANG) we have generated whole genomes from animals of a remote local pig population in Huambo province, which we have compared with 78 genomes of European and Asian pig breeds as well as European and Asian wild boars that are currently in public domain. Analyses of population structure showed that ANG pigs grouped within the European cluster and were clearly separated from Asian pig breeds. Pairwise F_ST ranged from 0.14 to 0.26, ANG pigs display lower levels of genetic differentiation towards European breeds. Finally, we have identified candidate regions for selection using a complementary approach based on various methods. All results suggest that selection towards feed efficiency and metabolism has occurred. Moreover, all analysis identified CDKAL1 gene, which is related with insulin and cholesterol metabolism, as a candidate gene overlapping signatures of selection unique to ANG pigs. This study presents the first assessment of the genetic relationship between ANG pigs and other world breeds and uncovers selection signatures that may indicate adaptation features unique to this important genetic resource.

Whole-genome sequencing reveals the artificial selection and local environmental adaptability of pigeons (Columba livia)

To meet human needs, domestic pigeons (Columba livia) with various phenotypes have been bred to provide genetic material for our research on artificial selection and local environmental adaptation. Seven pigeon breeds were resequenced and can be divided into commercial varieties (Euro-pigeon, Shiqi, Shen King, Taishen, and Silver King), ornamental varieties (High Fliers), and local varieties (Tarim pigeon). Phylogenetic analysis based on population resequencing showed that one group contained local breeds and ornamental pigeons from China, whereas all commercial varieties were clustered together. It is revealed that the traditional Chinese ornamental pigeon is a branch of Tarim pigeon. Runs of homozygosity (ROH) and linkage disequilibrium (LD) analyses revealed significant differences in the genetic diversity of the three types of pigeons. Genome sweep analysis revealed that the selected genes of commercial breeds were related to body size, reproduction, and plumage color. The genomic imprinting genes left by the ornamental pigeon breeds were mostly related to special human facial features and muscular dystrophy. The Tarim pigeon has evolved genes related to chemical ion transport, photoreceptors, oxidative stress, organ development, and olfaction in order to adapt to local environmental stress. This research provides a molecular basis for pigeon genetic resource evaluation and genetic improvement and suggests that the understanding of adaptive evolution should integrate the effects of various natural environmental characteristics.

Genomic signatures of a major adaptive event in the pathogenic fungus Melampsora larici-populina

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host–pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype–phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

Deletion in KRT75L4 linked to frizzle feather in Xiushui Yellow Chickens

Bird feathers are the product of interactions between natural and artificial selection. Feather-related traits are important for chicken selection and breeding. Frizzle feather is characterized by the abnormally development of feathers in chickens. In the current study, frizzle feather characteristics were observed in a local breed called Xiushui Yellow Chicken in Jiangxi, China. To determine the molecular mechanisms that underlie frizzle feather in Xiushui Yellow Chicken, four populations of three breeds (Xiushui Yellow Chicken with frizzle feathers, Xiushui Yellow Chicken with normal feathers, Guangfeng White-Ear Yellow Chicken, and Ningdu Yellow Chicken) were selected for whole-genome resequencing. Using a comparative genome strategy and genome-wide association study, a missense mutation (g.5281494A>G) and a 15-bp deletion (g.5285437-5285451delGATGCCGGCAGGACG) in KRT75L4 were identified as candidate mutations associated with frizzle feather in Xiushui Yellow Chicken. Based on genotyping performed in a large Xiushui Yellow Chicken population, the g.5285437-5285451delGATGCCGGCAGGACG mutation in KRT75L4 was confirmed as the putative causative mutation of frizzle feather. These results deepen the understanding of the molecular mechanisms responsible for frizzle feather, as well as facilitating the molecular detection and selection of the feather phenotype in Xiushui Yellow Chickens.

Development and validation of simultaneous identification of 26 mammalian and poultry species by a multiplex assay

A multiplex PCR assay was developed to simultaneously identify 22 mammalian species (alpaca, Asiatic black bear, Bactrian camel, brown rat, cat, cattle, common raccoon, dog, European rabbit, goat, horse, house mouse, human, Japanese badger, Japanese wild boar, masked palm civet, pig, raccoon dog, red fox, sheep, Siberian weasel, and sika deer) and four poultry species (chicken, domestic turkey, Japanese quail, and mallard), even from a biological sample containing a DNA mixture of multiple species. The assay was designed to identify species through multiplex PCR and capillary electrophoresis, with a combination of amplification of a partial region of the mitochondrial D-loop by universal primer sets and a partial region of the cytochrome b (cyt b) gene by species-specific primer sets. The assay was highly sensitive, with a detection limit of 100 copies of mitochondrial DNA. The assay's ability to identify species from complex DNA mixtures was demonstrated using an experimental sample consisting of 10 species. Efficacy, accuracy, and reliability of the assay were validated for use in forensic analysis with the guidelines of Scientific Working Group on DNA Analysis Methods (SWGDAM). The multiplex PCR assay developed in this study enables cost-effective, highly sensitive, and simultaneous species identification without massively parallel sequencing (MPS) platforms. Thus, the technique described is straightforward and suitable for routine forensic investigations.

Genomic and transcriptomic analyses of Citrus sinensis varieties provide insights into Valencia orange fruit mastication trait formation

Valencia orange (Citrus sinensis Osbeck) (VO) is a type of late-ripening sweet orange whose ripening occurs 4 to 5 months later than that of the mid-ripening common sweet orange (CO). Notably, the mastication trait of VO fruit is inferior to that of CO fruit. To date, how inferior pulp mastication trait forms in VO has not been determined. In this study, 13 VO varieties and 12 CO varieties were subjected to whole-genome resequencing. A total of 2.98 million SNPs were identified from 25 varieties, and a SNP molecular marker was developed to distinguish VO and CO. Moreover, 144 and 141 genes identified by selective sweep analysis were selected during VO and CO evolution, respectively. Based on gene functional enrichment analysis, most of the selected VO genes were related to the stress response and lignin biosynthesis. Simultaneously, we comparatively analyzed the transcriptome profiles of peel and pulp tissues among three VO varieties and three CO varieties, and the results demonstrated differences in lignin biosynthesis between VO and CO fruits. Furthermore, coexpression network analysis was performed to identify hub genes of lignin-related and variety-specific networks, which included CsERF74, CsNAC25, CsHSFB3, CsSPL4/13, etc. Overall, this study provides important insights into the mastication trait formation of Valencia orange fruit.

Genomic Variation Landscape of the Model Salt Cress Eutrema salsugineum

Eutrema salsugineum has long been used as the model for examining salt and other abiotic stress in plants. In addition to the forward genetics approaches widely used in the lab, natural variations undoubtedly will provide a rich genetic resource for studying molecular mechanisms underlying the stress tolerance and local adaptation of this species. We used 90 resequencing whole genomes of natural populations of this species across its Asian and North American distributions to detect the selection signals for genes involved in salt and other stresses at the species-range level and local distribution. We detected selection signals for genes involved in salt and other abiotic tolerance at the species level. In addition, several cold-induced and defense genes showed selection signals due to local adaptation in North America-NE Russia or northern China, respectively. These variations and findings provide valuable resources for further deciphering genetic mechanisms underlying the stress tolerance and local adaptations of this model species.

Maternal lineage of Nicobari pig (Sus scrofa nicobaricus) correlated with migration of Nicobarese, a native tribal population of Andaman and Nicobar Islands, India

Nicobari pig is reared by Nicobarese, a native tribal population of Andaman and Nicobar Islands. Nicobari pig has maintained its genetic identity due to geographical isolation. This communication is the first report on maternal inheritance of Nicobari pigs. DNA polymorphism data showed seven haplotypes. D-loop sequence information and mitogenome analysis were able to earmark Nicobari pigs to Asian clade. The domestication process of pigs and its expansion pattern help to understand human migration pattern. Based on this hypothesis, this communication elucidates the probable origin of Nicobarese. Earlier studies indicated that Nicobarese had genetic affinities to races distributed in China, Malaysia and Thailand. Our data on maternal inheritance of Nicobari pig correlates with the data on migration of Nicobarese. Moreover, we could establish a novel connection of Nicobarese with people of Northeastern parts of India, Philippines and Vietnam through phylogenetic signal and geographical provenance of Nicobari pig. We further concluded that migration of Nicobarese happened during Western route of migration (WRM) ∼4000 years before present. Therefore, we propose one wave hypothesis of peopling of Nicobar based on our study and existence of Ausrtroasiatic language, Mon-Khmer in these islands.

Whole-genome resequencing provides insights into the population structure and domestication signatures of ducks in eastern China

Background

Duck is an ancient domesticated animal with high economic value, used for its meat, eggs, and feathers. However, the origin of indigenous Chinese ducks remains elusive. To address this question, we performed whole-genome resequencing to first explore the genetic relationship among variants of these domestic ducks with their potential wild ancestors in eastern China, as well as understand how the their genomes were shaped by different natural and artificial selective pressures.

Results

Here, we report the resequencing of 60 ducks from Chinese spot-billed ducks (Anas zonorhyncha), mallards (Anas platyrhnchos), Fenghua ducks, Shaoxing ducks, Shanma ducks and Cherry Valley Pekin ducks of eastern China (ten from each population) at an average effective sequencing depth of ~ 6× per individual. The results of population and demographic analysis revealed a deep phylogenetic split between wild (Chinese spot-billed ducks and mallards) and domestic ducks. By applying selective sweep analysis, we identified that several candidate genes, important pathways and GO categories associated with artificial selection were functionally related to cellular adhesion, type 2 diabetes, lipid metabolism, the cell cycle, liver cell proliferation, and muscle functioning in domestic ducks.

Conclusion

Genetic structure analysis showed a close genetic relationship of Chinese spot-billed ducks and mallards, which supported that Chinese spot-billed ducks contributed to the breeding of domestic ducks. During the long history of artificial selection, domestic ducks have developed a complex biological adaptation to captivity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07710-2.

The History of Lentil (Lens culinaris subsp. culinaris) Domestication and Spread as Revealed by Genotyping-by-Sequencing of Wild and Landrace Accessions

Protein-rich legumes accompanied carbohydrate-rich cereals since the beginning of agriculture and yet their domestication history is not as well understood. Lentil (Lens culinaris Medik. subsp. culinaris) was first cultivated in Southwest Asia (SWA) 8000-10,000 years ago but archeological evidence is unclear as to how many times it may have been independently domesticated, in which SWA region(s) this may have happened, and whether wild species within the Lens genus have contributed to the cultivated gene pool. In this study, we combined genotyping-by-sequencing (GBS) of 190 accessions from wild (67) and domesticated (123) lentils from the Old World with archeological information to explore the evolutionary history, domestication, and diffusion of lentils to different environments. GBS led to the discovery of 87,647 single-nucleotide polymorphisms (SNPs), which allowed us to infer the phylogeny of genus Lens. We confirmed previous studies proposing four groups within it. The only gene flow detected was between cultivated varieties and their progenitor (L. culinaris subsp. orientalis) albeit at very low levels. Nevertheless, a few putative hybrids or naturalized cultivars were identified. Within cultivated lentil, we found three geographic groups. Phylogenetics, population structure, and archeological data coincide in a scenario of protracted domestication of lentils, with two domesticated gene pools emerging in SWA. Admixed varieties are found throughout their range, suggesting a relaxed selection process. A small number of alleles involved in domestication and adaptation to climatic variables were identified. Both novel mutation and selection on standing variation are presumed to have played a role in adaptation of lentils to different environments. The results presented have implications for understanding the process of plant domestication (past), the distribution of genetic diversity in germplasm collections (present), and targeting genes in breeding programs (future).

Selective sweep analysis reveals extensive parallel selection traits between large white and Duroc pigs

In the process of pig genetic improvement, different commercial breeds have been bred for the same purpose, improving meat production. Most of the economic traits, such as growth and fertility, have been selected similarly despite the discrepant selection pressure, which is known as parallel selection. Here, 28 whole-genome sequencing data of Danish large white pigs with an approximately 25-fold depth each were generated, resulting in about 12 million high-quality SNPs for each individual. Combined with the sequencing data of 27 Duroc and 23 European wild boars, we investigated the parallel selection of Danish large white and Duroc pigs using two complementary methods, Fst and iHS. In total, 67 candidate regions were identified as the signatures of parallel selection. The genes in candidate regions of parallel selection were mainly associated with sensory perception, growth rate, and body size. Further functional annotation suggested that the striking consistency of the terms may be caused by the polygenetic basis of quantitative traits, and revealing the complex genetic basis of parallel selection. Besides, some unique terms were enriched in population-specific selection regions, such as the limb development-related terms enriched in Duroc-specific selection regions, suggesting unique selections of breed-specific selected traits. These results will help us better understand the parallel selection process of different breeds. Moreover, we identified several potential causal SNPs that may contribute to the pig genetic breeding process.

Genomic and Phenotypic Divergence in Wild Barley Driven by Microgeographic Adaptation

Microgeographic adaptation is a fundamental driving force of evolution, but the underlying causes remain undetermined. Here, the phenotypic, genomic and transcriptomic variations of two wild barley populations collected from sharply divergent and adjacent micro-geographic sites to identify candidate genes associated with edaphic local adaptation are investigated. Common garden and reciprocal transplant studies show that large phenotypic differentiation and local adaptation to soils occur between these populations. Genetic, phylogenetic and admixture analyses based on population resequencing show that significant genetic divergences occur between basalt and chalk populations. These divergences are consistent with the phenotypic variations observed in the field. Genome sweep analyses reveal 162.7 Mb of selected regions driven by edaphic local adaptation, in which 445 genes identified, including genes associated with root architecture, metal transport/detoxification, and ABA signaling. When the phenotypic, genomic and transcriptomic data are combined, HvMOR, encoding an LBD transcription factor, is determined to be the vital candidate for regulating the root architecture to adapt to edaphic conditions at the microgeographic scale. This study provides new insights into the genetic basis of edaphic adaptation and demonstrates that edaphic factors may contribute to the evolution and speciation of barley.

Examining the effect of feralization on craniomandibular morphology in pigs, Sus scrofa (Artiodactyla: Suidae)

Feralization is the process by which domestic animals return to the wild and produce self-sustaining populations. It is often considered as a model in understanding the permanence of morphological changes associated with domestication; however, it is still unclear how much the release of anthropogenic selective pressures affects domestic traits. Here, we assessed the influence of feralization on the domestic morphological traits acquired through selective breeding using craniomandibular differences in shape and size between populations of feral pigs, wild boar and domestic pigs, using landmark-based geometric morphometrics. Our results suggest that numerous cranial and mandibular traits associated with domestication still exist in feral specimens, corroborating that domestication-induced changes in the shape of morphological elements are broadly maintained in feral populations. This is not the case for size variations, however, as the cranium is significantly smaller in feral pigs than in domesticated breeds, which could be due to the selective pressures associated with founding events. Our exploratory study, therefore, underlines the complexity of feral population history, the intricate influence of variations in genetic diversity, and novel selection pressures in the morphology of these groups. Future studies will need to expand the sample to take into account the diversity of morphotypes.

Genome-wide investigation of a Korean synthetic breed, Woori-Heukdon using the Illumina PorcineSNP60K BeadChip

BACKGROUND

Woori-Heukdon (KWH) is a Korean synthetic pig breed generated using Chookjin-Duroc (KCD), Chookjin-Chamdon (KCC), and their crossbreds. Currently, there is a severe lack of studies investigating the Korean breed populations including wild boars (KWB) throughout the genome.

OBJECTIVE

This study was performed to investigate the genetic characteristics of Korean pig populations at the genome-wide level.

METHODS

Using the SNP dataset derived from genotyped and downloaded datasets using the Illumina PorcineSNP60K BeadChip, we compared the genomes of 532 individuals derived from 23 pig breeds to assess the genetic diversity, inbreeding coefficient, genetic differentiation, and population structure.

RESULTS

KWB showed the lowest average expected heterozygosity (H_E = 0.1904), while KWH showed the highest genetic diversity (H_E = 0.02859) among Korean populations. We verified that the genetic composition of KWH, showing USD of 74.8% and KCC of 25.2% in ADMIXTURE analysis. In population structure analyses, KCC was consistently shown to be separated from other pig populations. In addition, we observed gene flow from Western pigs to a part of Chinese populations.

CONCLUSION

This study showed that Korean native pigs, KCC have genetic differences in comparison with Chinese and Western pigs; despite some historical records and recent genetic studies, we could not find any clear evidence that KCC was significantly influenced by Chinese or Western breeds in this study. We also verified the theoretical genomic composition of KWH at the molecular level in structure analyses. To our knowledge, this is the first genomic study to investigate the genomic characteristics of KWH and KCC.

Phylogeography of feral Monteiro pig in the Brazilian Pantanal Ecosystem

The Monteiro is a feral pig found in the Brazilian Pantanal ecosystem. The goal of this research is to generate data and knolewdge related to animal populations wich can be used for management and development of an in vitro conservation program for animal resourses at Pantanal ecosystem. The present study evaluated animals sampled from 10 distinct locations within the region, using 19 microsatellite markers (N = 189) and the control region of mitochondrial DNA (mtDNA) (N = 392). Low genetic differences were found between populations with the microsatellite data. The F_ST range was between 0.009 and 0.063 (p-value < 0.05). The Mantel test corroborated with previous results, as low correlations between genetic and geographic distances were observed (r² = 0.2309, p = 0.06). Bayesian analysis for genetic structure identification placed the Monteiro pigs into three main clusters (MOB, Pop 1 and all others Pantanal populations). Most of the Monteiro pigs share a single European haplotype as seen by mtDNA analyses. This haplotype is not exclusive, as it is shared with other swine populations (commercial and other locally adapted breeds). Monteiro populations from different geographic locations within Pantanal are not isolated and can be considered as a large unique population. Since animals roam freely to seek food and water, or even due to seasonal flooding of their habitat, the Monteiro populations presented absence of major genetic structure and evidence of high gene flow. These results can be used to create a management plan and in situ and ex situ conservation program for conservation and use of the Monteiro breed in the Pantanal ecosystem.

Whole genome sequence analysis reveals genetic structure and X-chromosome haplotype structure in indigenous Chinese pigs

Chinese indigenous pigs exhibit considerable phenotypic diversity, but their population structure and the genetic basis of agriculturally important traits need further exploration. Here, we sequenced the whole genomes of 24 individual pigs representing 22 breeds distributed throughout China. For comparison with European and commercial breeds (one pig per breed), we included seven published pig genomes with our new genomes for analyses. Our results showed that breeds grouped together based on morphological classifications are not necessarily more genetically similar to each other than to breeds from other groups. We found that genetic material from European pigs likely introgressed into five Chinese breeds. We have identified two new subpopulations of domestic pigs that encompass morphology-based criteria in China. The Southern Chinese subpopulation comprises the classical South Chinese Type and part of the Central China Type. In contrast, the Northern Chinese subpopulation comprises the North China Type, the Lower Yangtze River Basin Type, the Southwest Type, the Plateau Type, and the remainder of the Central China Type. Eight haplotypes and two recombination sites were identified within a conserved 40.09 Mb linkage-disequilibrium (LD) block on the X chromosome. Potential candidate genes (LEPR, FANCC, COL1A1, and PCCA) influencing body size were identified. Our findings provide insights into the phylogeny of Chinese indigenous pig breeds and benefit gene mining efforts to improve major economic traits.

Transcriptomic comparison of liver tissue between Anqing six-end-white pigs and Yorkshire pigs based on RNA sequencing

Chinese indigenous pig and Western commercial pig breeds show different patterns of lipid metabolism, fat deposition, and fatty acid composition; for these reasons, they have become vitally important models of energy metabolism and obesity in humans. To compare the mechanisms underlying lipid metabolism between Yorkshire pigs (lean type) and Anqing six-end-white pigs (obese type), the liver transcriptomes of six castrated boars with a body weight of approximately 100 kg (three Yorkshire and three Anqing) were analyzed by RNA-seq. The total number of reads produced for each liver sample ranged from 47.05 to 62.6 million. Among 362 differentially expressed genes, 142 were up-regulated and 220 were down-regulated in Anqing six-end-white pigs. Based on these data, 79 GO terms were significantly enriched. The top 10 (the 10 with lowest corrected P-value) significantly enriched GO terms were identified, including lipid metabolic process and carboxylic acid metabolic process. Pathway analysis revealed three significantly enriched KEGG pathways including PPAR signaling pathway, steroid hormone biosynthesis, and retinol metabolism. Based on protein–protein interaction networks, multiple genes responsible for lipid metabolism were identified, such as PCK1, PPARA, and CYP7A1, and these were considered promising candidate genes that could affect porcine liver lipid metabolism and fat deposition. Our results provide abundant transcriptomic information that will be useful for animal breeding and biomedical research.

Genome-wide analysis reveals molecular convergence underlying domestication in 7 bird and mammals

Background

In response to ecological niche of domestication, domesticated mammals and birds developed adaptively phenotypic homoplasy in behavior modifications like fearlessness, altered sociability, exploration and cognition, which partly or indirectly result in consequences for economic productivity. Such independent adaptations provide an excellent model to investigate molecular mechanisms and patterns of evolutionary convergence driven by artificial selection.

Results

First performing population genomic and brain transcriptional comparisons in 68 wild and domesticated chickens, we revealed evolutionary trajectories, genetic architectures and physiologic bases of adaptively behavioral alterations. To extensively decipher molecular convergence on behavioral changes thanks to domestication, we investigated selection signatures in hundreds of genomes and brain transcriptomes across chicken and 6 other domesticated mammals. Although no shared substitution was detected, a common enrichment of the adaptive mutations in regulatory sequences was observed, presenting significance to drive adaptations. Strong convergent pattern emerged at levels of gene, gene family, pathway and network. Genes implicated in neurotransmission, semaphorin, tectonic protein and modules regulating neuroplasticity were central focus of selection, supporting molecular repeatability of homoplastic behavior reshapes. Genes at nodal positions in trans-regulatory networks were preferably targeted. Consistent down-regulation of majority brain genes may be correlated with reduced brain size during domestication. Up-regulation of splicesome genes in chicken rather mammals highlights splicing as an efficient way to evolve since avian-specific genomic contraction of introns and intergenics. Genetic burden of domestication elicits a general hallmark. The commonly selected genes were relatively evolutionary conserved and associated with analogous neuropsychiatric disorders in human, revealing trade-off between adaption to life with human at the cost of neural changes affecting fitness in wild.

Conclusions

After a comprehensive investigation on genomic diversity and evolutionary trajectories in chickens, we revealed basis, pattern and evolutionary significance of molecular convergence in domesticated bird and mammals, highlighted the genetic basis of a compromise on utmost adaptation to the lives with human at the cost of high risk of neurophysiological changes affecting animals’ fitness in wild.

Whole Genome Resequencing Reveals Selection Signatures Associated With Important Traits in Ethiopian Indigenous Goat Populations

Ethiopia is considered as the main gateway for the introduction of livestock species, including goat, to the African continent. Ethiopian goats are characterized by their unique adaptive ability, and different physical characteristics in terms of morphology, body size, coat colors, and other important traits. The comparative population genomic analysis provides useful genomic information associated with important traits. Whole-genome resequencing of 44 Ethiopian indigenous goats produced 16 million single-nucleotide polymorphisms (SNPs) as well as 123,577 insertions and deletions. Specifically, 11,137,576, 10,760,581, 10,833,847, 12,229,657 and 10,749,996 putative SNPs were detected in Abergelle, Afar, Begait, Central Highland and Meafure goat populations, respectively. In this study, we used population differentiation (F_ST) and pooled heterozygosity (H_P) Cbased approaches. From the F_ST analysis, we identified 480 outlier windows. The H_P approach detected 108 and 205 outlier windows for Abergelle, and Begait, respectively. About 11 and 5 genes under selective signals were common for both approaches that were associated with important traits. After genome annotation, we found 41 Gene ontology (GO) terms (12 in biological processes, 8 in cellular components and 11 in the molecular function) and 10 Kyoto Encyclopedia of Genes and Genomes pathways. Several of the candidate genes are involved in the reproduction, body weight, fatty acids, and disease related traits. Our investigation contributes to deliver valuable genetic information and paves the way to design conservation strategy, breed management, genetic improvement, and utilization programs. The genomic resources generated in the study will offer an opportunity for further investigations.

Optimal sequencing depth design for whole genome re-sequencing in pigs

BACKGROUND

As whole-genome sequencing is becoming a routine technique, it is important to identify a cost-effective depth of sequencing for such studies. However, the relationship between sequencing depth and biological results from the aspects of whole-genome coverage, variant discovery power and the quality of variants is unclear, especially in pigs. We sequenced the genomes of three Yorkshire boars at an approximately 20X depth on the Illumina HiSeq X Ten platform and downloaded whole-genome sequencing data for three Duroc and three Landrace pigs with an approximately 20X depth for each individual. Then, we downsampled the deep genome data by extracting twelve different proportions of 0.05, 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9 paired reads from the original bam files to mimic the sequence data of the same individuals at sequencing depths of 1.09X, 2.18X, 3.26X, 4.35X, 6.53X, 8.70X, 10.88X, 13.05X, 15.22X, 17.40X, 19.57X and 21.75X to evaluate the influence of genome coverage, the variant discovery rate and genotyping accuracy as a function of sequencing depth. In addition, SNP chip data for Yorkshire pigs were used as a validation for the comparison of single-sample calling and multisample calling algorithms.

RESULTS

Our results indicated that 10X is an ideal practical depth for achieving plateau coverage and discovering accurate variants, which achieved greater than 99% genome coverage. The number of false-positive variants was increased dramatically at a depth of less than 4X, which covered 95% of the whole genome. In addition, the comparison of multi- and single-sample calling showed that multisample calling was more sensitive than single-sample calling, especially at lower depths. The number of variants discovered under multisample calling was 13-fold and 2-fold higher than that under single-sample calling at 1X and 22X, respectively. A large difference was observed when the depth was less than 4.38X. However, more false-positive variants were detected under multisample calling.

CONCLUSIONS

Our research will inform important study design decisions regarding whole-genome sequencing depth. Our results will be helpful for choosing the appropriate depth to achieve the same power for studies performed under limited budgets.

Discovery of Genomic Characteristics and Selection Signatures in Korean Indigenous Goats Through Comparison of 10 Goat Breeds

Indigenous breeds develop their own genomic characteristics by adapting to local environments or cultures over long periods of time. Most of them are not particularly productive in commercial terms, but they have abilities to survive in harsh environments or tolerate to specific diseases. Their adaptive characteristics play an important role as genetic materials for improving commercial breeds. As a step toward this goal, we analyzed the genome of Korean indigenous goats within 10 goat breeds. We collected 136 goat individuals by sequencing 46 new goats and employing 90 publicly available goats. Our whole-genome data was comprised of three indigenous breeds (Korean indigenous goat, Iranian indigenous goat, and Moroccan indigenous goat; n = 29, 18, 20), six commercial breeds (Saanen, Boer, Anglo-Nubian, British Alpine, Alpine, and Korean crossbred; n = 16, 11, 5, 5, 2, 13), and their ancestral species (Capra aegagrus; n = 17). We identified that the Iranian indigenous goat and the Moroccan indigenous goat have relatively similar genomic characteristics within a large category of genomic diversity but found that the Korean indigenous goat has unique genomic characteristics distinguished from the other nine breeds. Through population analysis, we confirmed that these characteristics have resulted from a near-isolated environment with strong genetic drift. The Korean indigenous goat experienced a severe genetic bottleneck upon entering the Korean Peninsula about 2,000 years ago, and has subsequently rarely experienced genetic interactions with other goat breeds. From selection analysis and gene-set enrichment analysis, we revealed selection signals for Salmonella infection and cardiomyopathy in the genome of the Korean indigenous goat. These adaptive characteristics were further identified with genomic-based evidence. We uncovered genomic regions of selective sweeps in the LBP and BPI genes (Salmonella infection) and the TTN and ITGB6 genes (cardiomyopathy), among several candidate genes. Our research presents unique genomic characteristics and distinctive selection signals of the Korean indigenous goat based on the extensive comparison. Although the adaptive traits require further validation through biological experiments, our findings are expected to provide a direction for future biodiversity conservation strategies and to contribute another option to genomic-based breeding programmes for improving the viability of Capra hircus.

Screening for signatures of selection of Tianzhu white yak using genome-wide re-sequencing

The Tianzhu white yak, a domestic yak indigenous to the Qilian Mountains, migrated inland from the Qinghai-Tibet Plateau. Specific ecological and long-term artificial selection influenced the evolution of its pure white coat and physiological characteristics. Therefore, it is not only a natural population that represents a genomic selective region of environmental adaptability but is also an animal model for studying the pigmentation of the yak coat. A total of 24 261 829 variants, including 22 445 252 SNPs, were obtained from 29 yaks by genome-wide re-sequencing. According to the results of a selective sweep analysis of Tianzhu white yak in comparison to Tibetan yaks, nine candidate genes under selection in Tianzhu white yak were identified by combining π, Tajima's D, πA/πB and F_ST statistics, with threshold standards of 5%. These genes include PDCD1, NUP210, ABCG8, NEU4, LOC102287650, D2HGDH, COL4A1, RTP5 and HDAC11. Five of the nine genes were classified into 12 molecular signaling pathways, and most of these signaling pathways are involved in environmental information processing, organismal systems and metabolism. A majority of these genes has not been implicated in previous studies of yak coat color and high-altitude animals. Our findings are helpful not only for explaining the molecular mechanism of yak coat pigmentation but also for exploring the genetic changes in Tianzhu white yak due to environmental adaptation.

Determination of genetic associations between indels in 11 candidate genes and milk composition traits in Chinese Holstein population

Background

We have previously identified 11 promising candidate genes for milk composition traits by resequencing the whole genomes of 8 Holstein bulls with extremely high and low estimated breeding values for milk protein and fat percentages (high and low groups), including FCGR2B, CENPE, RETSAT, ACSBG2, NFKB2, TBC1D1, NLK, MAP3K1, SLC30A2, ANGPT1 and UGDH those contained 25 indels between high and low groups. In this study, the purpose was to further examine whether these candidates have significant genetic effects on milk protein and fat traits.

Results

With PCR product sequencing, 13 indels identified by whole genome resequencing were successfully genotyped. With association analysis in 769 Chinese Holstein cows, we found that the indel in FCGR2B was significantly associated with milk yield, protein yield and protein percentage (P = 0.0041 to 0.0297); five indels in CENPE and one indel in MAP3K1 were markedly relevant to milk yield, fat yield and protein yield (P < 0.0001 to 0.0073); polymorphism in RETSAT was evidently associated with milk yield, fat yield, protein yield and protein percentage (P = 0.0001 to 0.0237); variant in ACSBG2 affected fat yield and protein percentage (P = 0.0088 and 0.0052); one indel in TBC1D1 was with respect to fat percentage and protein percentage (P = 0.0224 and 0.0209). Significant associations were shown between indels in NLK and protein yield and protein percentage (P = 0.0012 to 0.0257); variant in UGDH was related to the milk yield (P = 0.0312). The two exonic indels in FCGR2B and CENPE were predicted to change the mRNA and protein secondary structures, and resulted in the corresponding protein dysfunction.

Conclusion

Our findings presented here provide the first evidence for the associations of eight functional genes with milk yield and composition traits in dairy cattle.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0751-y) contains supplementary material, which is available to authorized users.

Genome-Wide Signatures of Selection Detection in Three South China Indigenous Pigs

South China indigenous pigs are famous for their superior meat quality and crude feed tolerance. Saba and Baoshan pigs without saddleback were located in the high-altitude area of Yunnan Province, while Tunchang and Ding’an pigs with saddleback were located in the low-altitude area of Hainan Province. Although these pigs are different in appearance, the underlying genetic differences have not been investigated. In this study, based on the single-nucleotide polymorphism (SNP) genotypes of 124 samples, both the cross-population extended haplotype homozygosity (XP-EHH) and the fixation index (F_ST) statistic were used to identify potential signatures of selection in these pig breeds. We found nine potential signatures of selection detected simultaneously by two methods, annotated 22 genes in Hainan pigs, when Baoshan pigs were used as the reference group. In addition, eleven potential signatures of selection detected simultaneously by two methods, annotated 24 genes in Hainan pigs compared with Saba pigs. These candidate genes were most enriched in GO: 0048015~phosphatidylinositol-mediated signaling and ssc00604: Glycosphingolipid biosynthesis—ganglio series. These selection signatures were likely to overlap with quantitative trait loci associated with meat quality traits. Furthermore, one potential selection signature, which was associated with different coat color, was detected in Hainan pigs. These results contribute to a better understanding of the underlying genetic architecture of South China indigenous pigs.

Whole-Genome Sequencing of Three Native Cattle Breeds Originating From the Northernmost Cattle Farming Regions

Northern Fennoscandia and the Sakha Republic in the Russian Federation represent the northernmost regions on Earth where cattle farming has been traditionally practiced. In this study, we performed whole-genome sequencing to genetically characterize three rare native breeds Eastern Finncattle, Western Finncattle and Yakutian cattle adapted to these northern Eurasian regions. We examined the demographic history, genetic diversity and unfolded loci under natural or artificial selection. On average, we achieved 13.01-fold genome coverage after mapping the sequencing reads on the bovine reference genome (UMD 3.1) and detected a total of 17.45 million single nucleotide polymorphisms (SNPs) and 1.95 million insertions-deletions (indels). We observed that the ancestral species (Bos primigenius) of Eurasian taurine cattle experienced two notable prehistorical declines in effective population size associated with dramatic climate changes. The modern Yakutian cattle exhibited a higher level of within-population variation in terms of number of SNPs and nucleotide diversity than the contemporary European taurine breeds. This result is in contrast to the results of marker-based cattle breed diversity studies, indicating assortment bias in previous analyses. Our results suggest that the effective population size of the ancestral Asiatic taurine cattle may have been higher than that of the European cattle. Alternatively, our findings could indicate the hybrid origins of the Yakutian cattle ancestries and possibly the lack of intensive artificial selection. We identified a number of genomic regions under selection that may have contributed to the adaptation to the northern and subarctic environments, including genes involved in disease resistance, sensory perception, cold adaptation and growth. By characterizing the native breeds, we were able to obtain new information on cattle genomes and on the value of the adapted breeds for the conservation of cattle genetic resources.

In silico approaches to discover the functional impact of non-synonymous single nucleotide polymorphisms in selective sweep regions of the Landrace genome

OBJECTIVE

The aim of this study was to discover the functional impact of non-synonymous single nucleotide polymorphisms (nsSNPs) that were found in selective sweep regions of the Landrace genome.

METHODS

Whole-genome re-sequencing data were obtained from 40 pigs, including 14 Landrace, 16 Yorkshire, and 10 wild boars, which were generated with the Illumina HiSeq 2000 platform. The nsSNPs in the selective sweep regions of the Landrace genome were identified, and the impacts of these variations on protein function were predicted to reveal their potential association with traits of the Landrace breed, such as reproductive capacity.

RESULTS

Total of 53,998 nsSNPs in the mapped regions of pigs were identified, and among them, 345 nsSNPs were found in the selective sweep regions of the Landrace genome which were reported previously. The genes featuring these nsSNPs fell into various functional categories, such as reproductive capacity or growth and development during the perinatal period. The impacts of amino acid sequence changes by nsSNPs on protein function were predicted using two in silico SNP prediction algorithms, i.e., sorting intolerant from tolerant and polymorphism phenotyping v2, to reveal their potential roles in biological processes that might be associated with the reproductive capacity of the Landrace breed.

CONCLUSION

The findings elucidated the domestication history of the Landrace breed and illustrated how Landrace domestication led to patterns of genetic variation related to superior reproductive capacity. Our novel findings will help understand the process of Landrace domestication at the genome level and provide SNPs that are informative for breeding.

Population genomic analysis reveals that homoploid hybrid speciation can be a lengthy process

An increasing number of species are thought to have originated by homoploid hybrid speciation (HHS), but in only a handful of cases are details of the process known. A previous study indicated that Picea purpurea, a conifer in the Qinghai-Tibet Plateau (QTP), originated through HHS from P. likiangensis and P. wilsonii. To investigate this origin in more detail, we analysed transcriptome data for 114 individuals collected from 34 populations of the three Picea species from their core distributions in the QTP. Phylogenetic, principal component and admixture analyses of nuclear SNPs showed the species to be delimited genetically and that P. purpurea was admixed with approximately 60% of its ancestry derived from P. wilsonii and 40% from P. likiangensis. Coalescent simulations revealed the best-fitting model of origin involved formation of an intermediate hybrid lineage between P. likiangensis and P. wilsonii approximately 6 million years ago (mya), which backcrossed to P. wilsonii to form P. purpurea approximately one mya. The intermediate hybrid lineage no longer exists and is referred to as a "ghost" lineage. Our study emphasizes the power of population genomic analysis combined with coalescent analysis for reconstructing the stages involved in the origin of a homoploid hybrid species over an extended period. In contrast to other studies, we show that these stages can in some instances span a relatively long period of evolutionary time.

Evidence of evolutionary history and selective sweeps in the genome of Meishan pig reveals its genetic and phenotypic characterization

Background

Meishan is a pig breed indigenous to China and famous for its high fecundity. The traits of Meishan are strongly associated with its distinct evolutionary history and domestication. However, the genomic evidence linking the domestication of Meishan pigs with its unique features is still poorly understood. The goal of this study is to investigate the genomic signatures and evolutionary evidence related to the phenotypic traits of Meishan via large-scale sequencing.

Results

We found that the unique domestication of Meishan pigs occurred in the Taihu Basin area between the Majiabang and Liangzhu Cultures, during which 300 protein-coding genes have underwent positive selection. Notably, enrichment of the FoxO signaling pathway with significant enrichment signal and the harbored gene IGF1R were likely associated with the high fertility of Meishan pigs. Moreover, NFKB1 exhibited strong selective sweep signals and positively participated in hyaluronan biosynthesis as the key gene of NF-kB signaling, which may have resulted in the wrinkled skin and face of Meishan pigs. Particularly, three population-specific synonymous single-nucleotide variants occurred in PYROXD1, MC1R, and FAM83G genes; the T305C substitution in the MCIR gene explained the black coat of the Meishan pigs well. In addition, the shared haplotypes between Meishan and Duroc breeds confirmed the previous Asian-derived introgression and demonstrated the specific contribution of Meishan pigs.

Conclusions

These findings will help us explain the unique genetic and phenotypic characteristics of Meishan pigs and offer a plausible method for their utilization of Meishan pigs as valuable genetic resources in pig breeding and as an animal model for human wrinkled skin disease research.

Integrated analysis of methylome, transcriptome and miRNAome of three pig breeds

AIM

Integrated analysis of methylome and transcriptome may help understand the molecular basis of the different breeds with different traits of commercial interest.

MATERIALS & METHODS

We obtained the first genome-wide methylome with single-base resolution, miRNAome and transcriptome from three swine breeds.

RESULTS

We displayed the landscape of the three omics in the whole-genome level. Integrated outcomes of methylome with genetic selection, miRNAome and transcriptome are also provided. Finally, we identified 11 candidate differentially methylated genes associated with phenotype variance in pigs.

CONCLUSION

DNA methylation not only suppresses transcriptome but also miRNAome. The different -omics data have complicated interaction in directly or indirectly and exhibited close relations with the distinct phenotypic traits of growth, disease resistance and energy metabolism.

Whole Genome Sequencing of Hulunbuir Short-Tailed Sheep for Identifying Candidate Genes Related to the Short-Tail Phenotype

The Hulunbuir short-tailed sheep (Ovis aries) is a breed native to China, in which the short-tail phenotype is the result of artificial and natural selection favoring a specific set of genetic mutations. Here, we analyzed the genetic differences between short-tail and normal-tail phenotypes at the genomic level. Selection signals were identified in genome-wide sequences. From 16 sheep, we identified 72,101,346 single nucleotide polymorphisms. Selection signals were detected based on the fixation index and heterozygosity. Seven genomic regions under putative selection were identified, and these regions contained nine genes. Among these genes, T was the strongest candidate as T is related to vertebral development. In T, a nonsynonymous mutation at c.G334T resulted in p.G112W substitution. We inferred that the c.G334T mutation in T leads to functional changes in Brachyury-encoded by this gene-resulting in the short-tail phenotype. Our findings provide a valuable insight into the development of the short-tail phenotype in sheep and other short-tailed animals.

Identification of selection signals by large-scale whole-genome resequencing of cashmere goats

Inner Mongolia and Liaoning cashmere goats are two outstanding Chinese multipurpose breeds that adapt well to the semi-arid temperate grassland. These two breeds are characterized by their soft cashmere fibers, thus making them great models to identify genomic regions that are associated with cashmere fiber traits. Whole-genome sequencing of 70 cashmere goats produced more than 5.52 million single-nucleotide polymorphisms and 710,600 short insertions and deletions. Further analysis of these genetic variants showed some population-specific molecular markers for the two cashmere goat breeds that are otherwise phenotypically similar. By analyzing F_ST and θ_π outlier values, we identified 135 genomic regions that were associated with cashmere fiber traits within the cashmere goat populations. These selected genomic regions contained genes, which are potential involved in the production of cashmere fiber, such as FGF5, SGK3, IGFBP7, OXTR, and ROCK1. Gene ontology enrichment analysis of identified short insertions and deletions also showed enrichment in keratinocyte differentiation and epidermal cell differentiation. These findings demonstrate that this genomic resource will facilitate the breeding of cashmere goat and other Capra species in future.

Genetic and phenotypic characterization of the novel mouse substrain C57BL/6N Korl with increased body weight

In inbred mouse lines, there is generally little genetic difference between individuals. This small genetic variability facilitates carrying out research on minute changes of various traits and the gene pool. Also, characterizing the diversity and detecting selective genetic and phenotypic signatures are crucial to understanding the genomic basis of a population and to identify specific patterns of evolutionary change. In this study, we investigated the underlying genetic profiles of a newly developed mouse strain, C57BL/6NKorl (Korl), established through sibling mating over 30 generations. To analyse the distinctive genomic features of Korl mice, we used whole-genome sequencing from six samples, which were compared to those of other C57BL/6N-based mouse strains. Korl strain-specific polymorphisms were identified and signatures of a selective sweep were detected. In particular, the candidate genes related to the increased body weight of the Korl strain were identified. Establishment of the genetic profile of Korl mice can provide insight into the inbreeding-induced changes to the gene pool, and help to establish this strain as a useful model for practical and targeted research purposes.

Self-domestication in Homo sapiens: Insights from comparative genomics

This study identifies and analyzes statistically significant overlaps between selective sweep screens in anatomically modern humans and several domesticated species. The results obtained suggest that (paleo-)genomic data can be exploited to complement the fossil record and support the idea of self-domestication in Homo sapiens, a process that likely intensified as our species populated its niche. Our analysis lends support to attempts to capture the "domestication syndrome" in terms of alterations to certain signaling pathways and cell lineages, such as the neural crest.

Effect of single nucleotide polymorphism on the total number of piglets born per parity of three different pig breeds

Objective

To determine the effects of genomic breeding values (GBV) and single nucleotide polymorphisms (SNP) on the total number of piglets born (TNB) in 3 pig breeds (Berkshire, Landrace, and Yorkshire).

Methods

After collecting genomic information (Porcine SNP BeadChip) and phenotypic TNB records for each breed, the effects of GBV and SNP were estimated by using single step best linear unbiased prediction (ssBLUP) method.

Results

The heritability estimates for TNB in Berkshire, Landrace, and Yorkshire breeds were 0.078, 0.107, and 0.121, respectively. The breeding value estimates for TNB in Berkshire, Landrace, and Yorkshire breeds were in the range of −1.34 to 1.47 heads, −1.79 to 1.87 heads, and −2.60 to 2.94 heads, respectively. Of sows having records for TNB, the reliability of breeding value for individuals with SNP information was higher than that for individuals without SNP information. Distributions of the SNP effects on TNB did not follow gamma distribution. Most SNP effects were near zero. Only a few SNPs had large effects. The numbers of SNPs with absolute value of more than 4 standard deviations in Berkshire, Landrace, and Yorkshire breeds were 11, 8, and 19, respectively. There was no SNP with absolute value of more than 5 standard deviations in Berkshire or Landrace. However, in Yorkshire, four SNPs (ASGA 0089457, ASGA0103374, ALGA0111816, and ALGA0098882) had absolute values of more than 5 standard deviations.

Conclusion

There was no common SNP with large effect among breeds. This might be due to the large genetic composition differences and the small size of reference population. For the precise evaluation of genetic performance of individuals using a genomic selection method, it may be necessary to establish the appropriate size of reference population.

Genome-wide sequence variations between wild and cultivated tomato species revisited by whole genome sequence mapping

Background

Cultivated tomato (Solanum lycopersicum L.) is the second most important vegetable crop after potato and a member of thirteen interfertile species of Solanum genus. Domestication and continuous selection for desirable traits made cultivated tomato species susceptible to many stresses as compared to the wild species. In this study, we analyzed and compared the genomes of wild and cultivated tomato accessions to identify the genomic regions that encountered changes during domestication.

Results

Analysis was based on SNP and InDel mining of twentynine accessions of twelve wild tomato species and forty accessions of cultivated tomato. Percentage of common SNPs among the accessions within a species corresponded with the reproductive behavior of the species. SNP profiles of the wild tomato species within a phylogenetic subsection varied with their geographical distribution. Interestingly, the ratio of genic SNP to total SNPs increased with phylogenetic distance of the wild tomato species from the domesticated species, suggesting that variations in gene-coding region play a major role in speciation. We retrieved 2439 physical positions in 1594 genes including 32 resistance related genes where all the wild accessions possessed a common wild variant allele different from all the cultivated accessions studied. Tajima’s D analysis predicted a very strong purifying selection associated with domestication in nearly 1% of its genome, half of which is contributed by chromosome 11. This genomic region with a low Tajima’s D value hosts a variety of genes associated with important agronomic trait such as, fruit size, tiller number and wax deposition.

Conclusion

Our analysis revealed a broad-spectrum genetic base in wild tomato species and erosion of that in cultivated tomato due to recurrent selection for agronomically important traits. Identification of the common wild variant alleles and the genomic regions undergoing purifying selection during cultivation would facilitate future breeding program by introgression from wild species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3822-3) contains supplementary material, which is available to authorized users.

The effects of Pleistocene climate change on biotic differentiation in a montane songbird clade from Wallacea

The role of Pleistocene Ice Age in tropical diversification is poorly understood, especially in archipelagos, in which glaciation-induced sea level fluctuations may lead to complicated changes in land distribution. To assess how Pleistocene land bridges may have facilitated gene flow in tropical archipelagos, we investigated patterns of diversification in the rarely-collected rusty-bellied fantail Rhipidura teysmanni (Passeriformes: Rhipiduridae) complex from Wallacea using a combination of bioacoustic traits and whole-genome sequencing methods (dd-RADSeq). We report a biogeographic leapfrog pattern in the vocalizations of these birds, and uncover deep genomic divergence among island populations despite the presence of intermittent land connections between some. We demonstrate how rare instances of genetic introgression have affected the evolution of this species complex, and document the presence of double introgressive mitochondrial sweeps, highlighting the dangers of using only mitochondrial DNA in evolutionary research. By applying different tree inference approaches, we demonstrate how concatenation methods can give inaccurate results when investigating divergence in closely-related taxa. Our study highlights high levels of cryptic avian diversity in poorly-explored Wallacea, elucidates complex patterns of Pleistocene climate-mediated diversification in an elusive montane songbird, and suggests that Pleistocene land bridges may have accounted for limited connectivity among montane Wallacean biota.

Genomewide association studies for hematological traits and T lymphocyte subpopulations in a Duroc × Erhualian F resource population

It has been shown that hematological traits can act as important indicators of immune function in both humans and livestock. T lymphocytes are key components of the adaptive immune system, playing a critical role in immune response. To identify genomic regions affecting hematological traits and T lymphocyte subpopulations, we performed both a SNP-based genomewide association study (GWAS) and a haplotype analysis for 20 hematological traits and 8 T cell subpopulations at 3 different time points (d 20, 33, and 35) in a Duroc × Erhualian F intercross population. Bonferroni correction was used to calculate the threshold -values for suggestive and 5% genomewide significance levels. In total, for SNP-based GWAS, we detected 96 significant SNP, including 15 genomewide-significant SNP, associated with 23 hematological traits and 234 significant SNP, including 27 genomewide-significant SNP, associated with 8 T cell subpopulations. Meanwhile, we identified 563 significant SNP, including 7 genomewide-significant SNP, associated with 5 hematological traits and 2,407 significant SNP, including 1,261 genomewide-significant SNP, associated with 8 T cell subpopulations by haplotype analysis. Among the significant regions detected, we propose both the () gene and the () gene on SSC3 as plausible candidate genes associated with CD/CD T lymphocytes at d 20. The findings provide insights into the basis of molecular mechanisms that are involved with immune response in the domestic pig and would aid further identification of causative mutations.

Whole-genome scanning for the litter size trait associated genes and SNPs under selection in dairy goat (Capra hircus)

Dairy goats are one of the most utilized domesticated animals in China. Here, we selected extreme populations based on differential fecundity in two Laoshan dairy goat populations. Utilizing deep sequencing we have generated 68.7 and 57.8 giga base of sequencing data, and identified 12,458,711 and 12,423,128 SNPs in the low fecundity and high fecundity groups, respectively. Following selective sweep analyses, a number of loci and candidate genes in the two populations were scanned independently. The reproduction related genes CCNB2, AR, ADCY1, DNMT3B, SMAD2, AMHR2, ERBB2, FGFR1, MAP3K12 and THEM4 were specifically selected in the high fecundity group whereas KDM6A, TENM1, SWI5 and CYM were specifically selected in the low fecundity group. A sub-set of genes including SYCP2, SOX5 and POU3F4 were localized both in the high and low fecundity selection windows, suggesting that these particular genes experienced strong selection with lower genetic diversity. From the genome data, the rare nonsense mutations may not contribute to fecundity, whereas nonsynonymous SNPs likely play a predominant role. The nonsynonymous exonic SNPs in SETDB2 and CDH26 which were co-localized in the selected region may take part in fecundity traits. These observations bring us a new insights into the genetic variation influencing fecundity traits within dairy goats.

Genomic analysis reveals selection in Chinese native black pig

Identification of genomic signatures that help reveal mechanisms underlying desirable traits in domesticated pigs is of significant biological, agricultural and medical importance. To identify the genomic footprints left by selection during domestication of the Enshi black pig, a typical native and meat-lard breed in China, we generated about 72-fold coverage of the pig genome using pools of genomic DNA representing three different populations of Enshi black pigs from three different locations. Combining this data with the available whole genomes of 13 Chinese wild boars, we identified 417 protein-coding genes embedded in the selected regions of Enshi black pigs. These genes are mainly involved in developmental and metabolic processes, response to stimulus, and other biological processes. Signatures of selection were detected in genes involved in body size and immunity (RPS10 and VASN), lipid metabolism (GSK3), male fertility (INSL6) and developmental processes (TBX19). These findings provide a window into the potential genetic mechanism underlying development of desirable phenotypes in Enshi black pigs during domestication and subsequent artificial selection. Thus, our results illustrate how domestication has shaped patterns of genetic variation in Enshi black pigs and provide valuable genetic resources that enable effective use of pigs in agricultural production.