Search

Thoracolumbar vertebra (TLV) and rib primordium (RP) development is a common evolutionary feature across vertebrates, although whole-organism analysis of TLV and RP gene expression dynamics has been lacking. Here we investigated the single-cell transcriptomic landscape of thoracic vertebra (TV), lumbar vertebra (LV), and RP cells from a pig embryo at 27 days post-fertilization (dpf) and identified six cell types with distinct gene-expression signatures. In-depth dissection of the gene-expression dynamics and RNA velocity revealed a coupled process of osteogenesis and angiogenesis during TLV and rib development. Further analysis of cell-type-specific and strand-specific expression uncovered the extremely high levels of Hoxa10 3'-untranslated region (UTR) sequence specific to osteoblast of LV cells, which may function as anti-Hoxa10-antisense by counteracting the Hoxa10-antisense effect to determine TLV transition. Thus, this work provides a valuable resource for understanding embryonic osteogenesis and angiogenesis underlying vertebrate TLV and RP development at the cell-type-specific resolution, which serves as a comprehensive view on the transcriptional profile of animal embryo development.

Obsessive-compulsive disorder (OCD) represents a heterogeneous collection of diseases with diverse levels of phenotypic, genetic, and etiologic variability, making it difficult to identify the underlying genetic and biological mechanisms in humans. Domestic dogs exhibit several OCD-like behaviors. Using continuous circling as a representative phenotype for OCD, we screened two independent dog breeds, the Belgian Malinois and Kunming Dog and subsequently sequenced ten circling dogs and ten unaffected dogs for each breed. Using population differentiation analyses, we identified 11 candidate genes in the extreme tail of the differentiated regions between cases and controls. These genes overlap significantly with genes identified in a genome wide association study (GWAS) of human OCD, indicating strong convergence between humans and dogs. Through gene expressional analysis and functional exploration, we found that two candidate OCD risk genes, PPP2R2B and ADAMTSL3, affected the density and morphology of dendritic spines. Therefore, changes in dendritic spine may underlie some common biological and physiological pathways shared between humans and dogs. Our study revealed an unprecedented level of convergence in OCD shared between humans and dogs, and highlighted the importance of using domestic dogs as a model species for many human diseases including OCD.

Understanding the zoonotic origin and evolution history of SARS-CoV-2 will provide critical insights for alerting and preventing future outbreaks. A significant gap remains for the possible role of pangolins as a reservoir of SARS-CoV-2 related coronaviruses (SC2r-CoVs). Here, we screened SC2r-CoVs in 172 samples from 163 pangolin individuals of four species, and detected positive signals in muscles of four Manis javanica and, for the first time, one M. pentadactyla. Phylogeographic analysis of pangolin mitochondrial DNA traced their origins from Southeast Asia. Using in-solution hybridization capture sequencing, we assembled a partial pangolin SC2r-CoV (pangolin-CoV) genome sequence of 22 895 bp (MP20) from the M. pentadactyla sample. Phylogenetic analyses revealed MP20 was very closely related to pangolin-CoVs that were identified in M. javanica seized by Guangxi Customs. A genetic contribution of bat coronavirus to pangolin-CoVs via recombination was indicated. Our analysis revealed that the genetic diversity of pangolin-CoVs is substantially higher than previously anticipated. Given the potential infectivity of pangolin-CoVs, the high genetic diversity of pangolin-CoVs alerts the ecological risk of zoonotic evolution and transmission of pathogenic SC2r-CoVs.

Over the last several hundred years, donkeys have adapted to high-altitude conditions on the Tibetan Plateau. Interestingly, the kiang, a closely related equid species, also inhabits this region. Previous reports have demonstrated the importance of specific genes and adaptive introgression in divergent lineages for adaptation to hypoxic conditions on the Tibetan Plateau. Here, we assessed whether donkeys and kiangs adapted to the Tibetan Plateau via the same or different biological pathways and whether adaptive introgression has occurred. We assembled a de novo genome from a kiang individual and analyzed the genomes of five kiangs and 93 donkeys (including 24 from the Tibetan Plateau). Our analyses suggested the existence of a strong hard selective sweep at the EPAS1 locus in kiangs. In Tibetan donkeys, however, another gene, i.e., EGLN1, was likely involved in their adaptation to high altitude. In addition, admixture analysis found no evidence for interspecific gene flow between kiangs and Tibetan donkeys. Our findings indicate that despite the short evolutionary time scale since the arrival of donkeys on the Tibetan Plateau, as well as the existence of a closely related species already adapted to hypoxia, Tibetan donkeys did not acquire adaptation via admixture but instead evolved adaptations via a different biological pathway.

Genomic imprinting often results in parent-of-origin specific differential expression of maternally and paternally inherited alleles and plays an essential role in mammalian development and growth. Mammalian genomic imprinting has primarily been studied in mice and humans, with only limited information available for pigs. To systematically characterize this phenomenon and evaluate imprinting status between different species, we investigated imprinted genes on a genome-wide scale in pig brain tissues. Specifically, we performed bioinformatics analysis of high-throughput sequencing results from parental genomes and offspring transcriptomes of hybrid crosses between Duroc and Diannan small-ear pigs. We identified 11 paternally and five maternally expressed imprinted genes in pigs with highly stringent selection criteria. Additionally, we found that the KCNQ1 and IGF2R genes, which are related to development, displayed a different imprinting status in pigs compared with that in mice and humans. This comprehensive research should help improve our knowledge on genomic imprinting in pigs and highlight the potential use of imprinted genes in the pig breeding field.

Solar energy, which is essential for the origin and evolution of all life forms on Earth, can be objectively recorded through attributes such as climatic ambient temperature (CAT), ultraviolet radiation (UVR), and sunlight duration (SD). These attributes have specific geographical variations and may cause different adaptation traits. However, the adaptation profile of each attribute and the selective role of solar energy as a whole during human evolution remain elusive. Here, we performed a genome-wide adaptation study with respect to CAT, UVR, and SD using the Human Genome Diversity Project-Centre Etude Polymorphism Humain (HGDP-CEPH) panel data. We singled out CAT as the most important driving force with the highest number of adaptive loci (6 SNPs at the genome-wide 1 × 10−7 level; 401 at the suggestive 1 × 10−5 level). Five of the six genome-wide significant adaptation SNPs were successfully replicated in an independent Chinese population (N = 1395). The corresponding 316 CAT adaptation genes were mostly involved in development and immunity. In addition, 265 (84%) genes were related to at least one genome-wide association study (GWAS)-mapped human trait, being significantly enriched in anthropometric loci such as those associated with body mass index (χ2; P

This study aimed to investigate to the effects of dietary CP levels and protease supplementation on growth performance, carcass traits, meat quality, nutrients utilization, and standardized ileal digestibility of amino acid in Pekin ducks fed a complex diet. A total of 960 14-day-old male ducks were weighed and randomly allotted to a 2 × 5 factorial arrangement of 10 treatments with 6 replicate pens per treatment and 16 ducks per pen fed to 49 D of age. Experimental factors included five dietary CP levels ranging from 13.5 to 17.5% and with or without protease (200 mg/kg) supplementation. Between day 28 to 34, the digestible and metabolizable trials were performed. Significant CP × protease interactions (P

Classic concepts of genetic (gene) diversity (heterozygosity) such as Nei & Li’s nucleotide diversity were defined within a population context. Although variations are often measured in population context, the basic carriers of variation are individuals. Hence, measuring variations such as SNP of an individual against a reference genome, which has been ignored previously, is certainly in its own right. Indeed, similar practice has been a tradition in community ecology, where the basic unit of diversity measure is individual community sample. We propose to use Renyi’s-entropy-based Hill numbers to define individual-level genetic diversity and similarity and demonstrate the definitions with the SNP (single nucleotide polymorphism) datasets from the 1000-Genomes Project. Hill numbers, derived from Renyi’s entropy (of which Shannon’s entropy is a special case), have found widely applications including measuring the quantum information entanglement and ecological diversity. The demonstrated individual-level SNP diversity not only complements the existing population-level genetic diversity concepts, but also offers building blocks for comparative genetic analysis at higher levels. The concept of individual covers, but is not limited to, individual chromosome, region of chromosome, gene cluster(s), or whole genome. Similarly, the SNP can be replaced by other structural variants or mutation types such as indels.

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

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

Cis-regulatory elements play an important role in the development of traits and disease in organisms (Ma et al., 2020; Woolfe et al., 2005) and their annotation could facilitate genetic studies. The Encyclopedia of DNA Elements (ENCODE) (Davis et al., 2018) and Functional Annotation of Animal Genomes (FAANG) (FAANG Consortium et al., 2015) offer pioneering data on regulatory elements in several species. Currently, however, regulatory element annotation data remain limited for most organisms. In this study, we developed a tool (OrthReg) for annotating conserved orthologous cis-regulatory elements in targeted genomes using an annotated reference genome. Cross-species validation of this annotation tool using human and mouse ENCODE data confirmed the robustness of this strategy. To explore the efficiency of the tool, we annotated the pig genome and identified more than 28 million regulatory annotation records using the reference human ENCODE data. With this regulatory annotation, some putative regulatory non-coding variants were identified within domestication sweeps in European and East Asian pigs. Thus, this tool can utilize data produced by ENCODE, FAANG, and similar projects, and can be easily extended to customized experimental data. The extensive application of this tool will help to identify informative single nucleotide polymorphisms (SNPs) in post-genome-wide association studies and resequencing analysis of organisms with limited regulatory annotation data.

Background Despite decades of research, the horse domestication scenario in East Asia remains poorly understood. Results The study identified 16 haplogroups with fine-scale phylogenetic resolution using mitochondrial genomes of 317 horse samples. The time to the most recent common ancestor of the 16 haplogroups ranges from [0.8–3.1] thousand years ago (KYA) to [7.9–27.1] KYA. With combined analyses of the mitochondrial control region for 35 extant Przewalski’s horses, 3544 modern and 203 ancient horses across the world, researchers provide evidence for that East Asian prevalent haplogroups Q and R were indigenously domesticated or they were involved in numerous distinct genetic components from wild horses in the southern part of East Asia. These events of haplotypes Q and R occurred during 4.7 to 16.3 KYA and 2.1 to 11.5 KYA, respectively. The diffusion of preponderant European haplogroups L from west to East Asia is consistent with the external gene input. Furthermore, genetic differences were detected between northern East Asia and southern East Asia cohorts by Principal Component Analysis, Analysis of Molecular Variance test, the χ2 test and phylogeographic analyses. Conclusions All results suggest a complex picture of horse domestication, as well as geographic pattern in East Asia. Both local origin and external input occurred in East Asia horse populations. And besides, there are at least two different domestication or hybridization centers in East Asia.

Domestication of the helmeted guinea fowl (HGF; Numida meleagris) in Africa remains elusive. Here we report a high-quality de novo genome assembly for domestic HGF generated by long- and short-reads sequencing together with optical and chromatin interaction mapping. Using this assembly as the reference, we performed population genomic analyses for newly sequenced whole-genomes for 129 birds from Africa, Asia, and Europe, including domestic animals (n = 89), wild progenitors (n = 34), and their closely related wild species (n = 6). Our results reveal domestication of HGF in West Africa around 1,300-5,500 years ago. Scanning for selective signals characterized the functional genes in behavior and locomotion changes involved in domestication of HGF. The pleiotropy and linkage in genes affecting plumage color and fertility were revealed in the recent breeding of Italian domestic HGF. In addition to presenting a missing piece to the jigsaw puzzle of domestication in poultry, our study provides valuable genetic resources for researchers and breeders to improve production in this species.

Coexistence and cooperation between dogs and humans over thousands of years have supported convergent evolutionary processes in the two species. Previous studies found that Eurasian dogs evolved into a distinct geographic cluster. In this study, we used the genomes of 242 European dogs, 38 Southeast Asian indigenous (SEAI) dogs, and 41 gray wolves to identify adaptation of European dogs . We report 86 unique positively selected genes in European dogs, among which is LCT (lactase). LCT encodes lactase, which is fundamental for the digestion of lactose. We found that an A-to-G mutation (chr19:38,609,592) is almost fixed in Middle Eastern and European dogs. The results of two-dimensional site frequency spectrum (2D SFS) support that the mutation is under soft sweep . We inferred that the onset of positive selection of the mutation is shorter than 6,535 years and behind the well-developed dairy economy in central Europe. It increases the expression of LCT by reducing its binding with ZEB1, which would enhance dog’s ability to digest milk-based diets. Our study uncovers the genetic basis of convergent evolution between humans and dogs with respect to diet, emphasizing the import of the dog as a biomedical model for studying mechanisms of the digestive system.

Although DNA array-based approaches for genome wide association studies (GWAS) permit the collection of thousands of low-cost genotypes, it is often at the expense of resolution and completeness, as SNP chip technologies are ultimately limited by SNPs chosen during array development. An alternative low-cost approach is low-pass whole genome sequencing (WGS) followed by imputation. Rather than relying on high levels of genotype confidence at a set of select loci, low-pass WGS and imputation relies on the combined information from millions of randomly sampled low confidence genotypes. To investigate low-pass WGS and imputation in the dog, we assessed accuracy and performance by downsampling 97 high-coverage (>15x) WGS datasets from 51 different breeds to approximately 1x coverage, simulating low-pass WGS. Using a reference panel of 676 dogs from 91 breeds, genotypes were imputed from the downsampled data and compared to a truth set of genotypes generated from high coverage WGS. Using our truth set, we optimized a variant quality filtering strategy that retained approximately 80% of 14M imputed sites and lowered the imputation error rate from 3.0% to 1.5%. Seven million sites remained with a MAF > 5% and an average imputation quality score of 0.95. Finally, we simulated the impact of imputation errors on outcomes for case-control GWAS, where small effect sizes were most impacted and medium to large effect sizes were minorly impacted. These analyses provide best practice guidelines for study design and data post-processing of low-pass WGS imputed genotypes in dogs.

Study question To investigate the different metabolomic profiling in serum between pregnant and non-pregnant women during early implantation period. Summary answer Metabolomics of progesterone-related hormones enhances from ET day3 for pregnancy women compared with non-pregnancy women. What is known already Metabolomics is based on high-throughput analytical methods to identify and quantify metabolites. Compared to other omics study, metabolomics is the closest one to the phenotype, allowing the observation of dynamic changes in phenotype at specific timepoints. So far there is no published work about the metabolomics profile in human early implantation period. Study design, size, duration: Study design: comparative study. Size: 14 pregnancy women and 14 non-pregnancy women. duration: time-course. Participants/materials, setting, methods Participants: pregnancy women and unpregnancy women after embryo transfer (ET). Setting: university-based study. Methods: Peripheral blood were collected at ET day0, 3, 6 and 9. metabolomic profiling in serum by platforms of capillary electrophoresis-mass spectrometry (CE-MS) and liquid chromatography–mass spectrometry (LC-MS). Main results and the role of chance There were no statistical difference of the age, BMI, basal FSH level, endometrium thickness on the day of embryo transfer, distribution of primary and secondary fertility, embryo transfer cycle as well as the infertile types between the two groups. After deleting those with over 50% missing data, we finally have 310 metabolites into statistical analysis. Among the 310 metabolite, lipid metabolites account the largest percentage, nearly half of all metabolites. The second biggest class of metabolites in our data was organic acids. Combined results in repeated measurement ANOVA (RM-ANOVA) and ANOVA-simultaneous component analysis (ASCA) as well as multivariate empirical Bayes time-series analysis (MEBA), we finally found that progesterone-related hormones were the most important metabolites for the whole time-series data. Those significant metabolites showed a significant down regulation from ET day0 to ET day3 and up regulation from ET day3 to ET day9. Limitations, reasons for caution we have limited sample size for this study and further validation is necessary for confirmation. Wider implications of the findings: The phenomenon of upregulation of progesterone-related hormones from day3 in pregnancy group might be related to the embryo-originated hcg. Because the embryo has entered into endometrium at day3 and produced cytokines, hcg and other interaction with endometrium. Trial registration number NA

The genus Dichaetophora is currently comprised of 67 formally described species assigned into five species groups, i.e., agbo, tenuicauda, acutissima, sinensis and trilobita. In the present study, we challenged to delimit species from a huge amount of samples of Dichaetophora and allied taxa (the genus Mulgravea and the subgenus Dudaica of Drosophila) collected from a wide range of the Oriental and east Palearctic regions. We first sorted all specimens into morpho-species, which were tentatively classified into the species groups of Dichaetophora, Mulgravea and Dudaica based on the morphological diagnoses for these supraspecific taxa. Then, we selected representative specimen(s) for each morpho-species and subjected them to barcoding of COI (the cytochrome c oxidase subunit I gene) sequences. The applied ABGD (automatic barcode gap discovery) algorithm estimated a total of 222 MOTUs (molecular operational taxonomic units). Although most of the morpho-species were recognized as individual MOTUs, some others were divided into multiple MOTUs or combined into a single MOTU, suggesting the existence of cryptic sibling species, subspecies or intraspecific genetic varieties. Out of the 222 MOTUs, 88 representing the supraspecific taxa of Dichaetophora, Mulgravea and Dudaica were selected, along with 33 species from major genera and subgenera of Drosophila in the tribe Drosophilini, as in-group (four species from the tribe Colocasiomyini as out-group) for phylogenetic reconstruction. We analyzed a dataset of concatenated nucleotide sequences of 12 nuclear gene markers by the maximum likelihood and Bayesian inference methods. As a result, the three focal taxa (i.e., Dichaetophora, Mulgravea and Dudaica) formed a clade, which we called the “pan-Dichaetophora”. Within this large clade, the agbo, tenuicauda, sinensis and trilobita groups of Dichaetophora, Mulgravea and Dudaica were recovered as monophyletic groups, but Dichaetophora and its acutissima group were regarded as paraphyletic. In addition, two clusters were recognized among ungrouped MOTUs of Dichaetophora. Thus, the present study has uncovered some problems in the taxonomy of the pan-Dichaetophora. Solving such problems will be challenged in Part II of this serial work.

Thoracolumbar vertebra (TLV) and rib primordium (RP) development is a common evolutionary feature across vertebrates although whole-organism analysis of TLV and RP gene expression dynamics has been lacking. Here we investigated the single-cell transcriptomic landscape of thoracic vertebra (TV), lumbar vertebra (LV), and RP cells from a pig embryo at 27 days post-fertilization (dpf) and identified six cell types with distinct gene-expression signatures. In-depth dissection of the gene-expression dynamics and RNA velocity revealed a coupled process of osteogenesis and angiogenesis during TLV and rib development. Further analysis of cell-type-specific and strand-specific expression uncovered the extremely high levels of HOXA10 3’-UTR sequence specific to osteoblast of LV cells, which may function as anti-HOXA10-antisense by counteracting the HOXA10-antisense effect to determine TLV transition. Thus, this work provides a valuable resource for understanding embryonic osteogenesis and angiogenesis underlying vertebrate TLV and RP development at the cell-type-specific resolution, which serves as a comprehensive view on the transcriptional profile of animal embryo development.

The spread of COVID-19, the disease caused by SARS-CoV-2, has progressed into a global pandemic. To date, thousands of genetic variants have been identified across strains of SARS-CoV-2 isolated from worldwide patients. However, there is still little direct evidence linking viral variants and clinical features. Based on two tightly linked SNPs, we previously divided SARS-CoV-2 into two major lineages: the ancestral ‘S lineage’ (U8,782 and C28,144) and the derived ‘L lineage’ (C8,782 and U28,144). Here, we identified SARS-CoV-2 lineages from 271 COVID-19 patients during the early outbreak of this pandemic in Wuhan, including 73 S- and 198 L-lineage patients. The S-lineage patients exhibited significantly worse clinical severity than the L-lineage patients, even after excluding other risk factors. This study suggests that SARS-CoV-2 lineage may provide useful clinical information for the management of COVID-19 and support the argument that viral clades should be analyzed as a function of clinical severity.

Domestication of the helmeted guinea fowl (HGF; Numida meleagris) in Africa remains elusive. Here we report a high-quality de novo genome assembly for domestic HGF generated by long- and short-reads sequencing together with optical and chromatin interaction mapping. Using this assembly as the reference, we performed population genomic analyses for newly sequenced whole-genomes for 129 birds from Africa, Asia, and Europe, including domestic animals (n = 89), wild progenitors (n = 34), and their closely related wild species (n = 6). Our results reveal domestication of HGF in West Africa around 1,300–5,500 years ago. Scanning for selective signals characterized the functional genes in behavior and locomotion changes involved in domestication of HGF. The pleiotropy and linkage in genes affecting plumage color and fertility were revealed in the recent breeding of Italian domestic HGF. In addition to presenting a missing piece to the jigsaw puzzle of domestication in poultry, our study provides valuable genetic resources for researchers and breeders to improve production in this species.

The phylogeography of the European wild boar was mainly determined by postglacial recolonization patterns from Mediterranean refugia after the last ice age. Here we present the first analysis of SNP polymorphism within the complete mtDNA genome of West Russian (n = 8), European (n = 64), and North African (n = 5) wild boar. Our analyses provided evidence of unique lineages in the East‐Caucasian (Dagestan) region and in Central Italy. A phylogenetic analysis revealed that these lineages are basal to the other European mtDNA sequences. We also show close connection between the Western Siberian and Eastern European populations. Also, the North African samples were clustered with the Iberian population. Phylogenetic trees and migration modeling revealed a high proximity of Dagestan sequences to those of Central Italy and suggested possible gene flow between Western Asia and Southern Europe which was not directly related to Northern and Central European lineages. Our results support the presence of old maternal lineages in two Southern glacial refugia (i.e., Caucasus and the Italian peninsula), as a legacy of an ancient wave of colonization of Southern Europe from an Eastern origin.

Tibetan pig is native to the Qinghai-Tibet Plateau and has adapted to the high-altitude environmental condition such as hypoxia. However, its origin and genetic mechanisms underlying high-altitude adaptation still remain controversial and enigmatic. Herein, we analyze 229 genomes of wild and domestic pigs from Eurasia, including 63 Tibetan pigs, and detect 49.6 million high-quality variants. Phylogenomic and structure analyses show that Tibetan pigs have a close relationship with low-land domestic pigs in China, implying a common domestication origin. Positively selected genes in Tibetan pigs involved in high-altitude physiology, such as hypoxia, cardiovascular systems, UV damage, DNA repair. Three of loci with strong signals of selection are associated with EPAS1, CYP4F2, and THSD7A genes, related to hypoxia and circulation. We validated four non-coding mutations nearby EPAS1 and CYP4F2 showing reduced transcriptional activity in Tibetan pigs. A high-frequency missense mutation is found in THSD7A (Lys561Arg) in Tibetan pigs. The selective sweeps in Tibetan pigs was found in association with selection against non-coding variants, indicating an important role of regulatory mutations in Tibetan pig evolution. This study is important in understanding the evolution of Tibetan pigs and advancing our knowledge on animal adaptation to high-altitude environments.

Indigenous chickens at the Swahili coast and other traditional migratory corridors in Kenya represent important populations that are inconclusively characterized. Using a comprehensive dataset of Kenyan indigenous chickens and additional mined data of chickens from 8 African and 5 Asian countries, we performed univariate and multivariate assessments to uncover the underlying phenotypic and morphometric variability. Kenyan indigenous chickens expressed differentiation of several qualitative and quantitative traits, both among different counties in the Swahili coast, and among coastal, western, and northern migratory corridors. There was a substantial population stratification of these chickens, particularly distinctive clustering of chickens from Marsabit, Lamu, and Kilifi counties. The pooled dataset further clarified a closer phenotypic and morphometric proximity of chickens within different geographical regions. We additionally revealed a preponderance of bantam and rumpless traits to hot and humid locales, and feathered shanks to cooler regions. Currently, most chicken breeding programs in developing countries rely on phenotypic and morphometric properties. Hence, the high chicken diversity and population stratification observed in our study, possibly shaped by natural and artificial selective pressures, reveal opportunities for complementary phenotypic and genotypic assessments to identify resources for effective breed improvement and conservation strategies of indigenous chickens in the tropics.

Viviparous (live-bearing) vertebrates have evolved repeatedly within otherwise oviparous (egg-laying) clades. Over two-thirds of these changes in vertebrate reproductive parity mode happened in squamate reptiles, where the transition has happened between 98 and 129 times. The transition from oviparity to viviparity requires numerous physiological, morphological, and immunological changes to the female reproductive tract, including eggshell reduction, delayed oviposition, placental development for supply of water and nutrition to the embryo by the mother, enhanced gas exchange, and suppression of maternal immune rejection of the embryo. We performed genomic and transcriptomic analyses of a closely related oviparous–viviparous pair of lizards ( Phrynocephalus przewalskii and Phrynocephalus vlangalii ) to examine these transitions. Expression patterns of maternal oviduct through reproductive development of the egg and embryo differ markedly between the two species. We found changes in expression patterns of appropriate genes that account for each of the major aspects of the oviparity to viviparity transition. In addition, we compared the gene sequences in transcriptomes of four oviparous–viviparous pairs of lizards in different genera ( Phrynocephalus , Eremias , Scincella , and Sphenomorphus ) to look for possible gene convergence at the sequence level. We discovered low levels of convergence in both amino acid replacement and evolutionary rate shift. This suggests that most of the changes that produce the oviparity–viviparity transition are changes in gene expression, so occasional reversals to oviparity from viviparity may not be as difficult to achieve as has been previously suggested.

Background The ability to digest dietary lactose is associated with lactase persistence (LP) in the intestinal lumen in human. The genetic basis of LP has been investigated in many populations in the world. Iran has a long history of pastoralism and the daily consumption of dairy products; thus, we aim to assess how LP has evolved in the Iranian population. We recruited 400 adult individuals from seven Iranian ethnic groups, from whom we investigated their lactose tolerance and screened the genetic variants in their lactase gene locus. Results The LP frequency distribution ranged from 0 to 29.9% in the seven Iranian ethnic groups with an average value of 9.8%. The variants, − 13910*T and − 22018*A, were significantly associated with LP phenotype in Iranians. We found no evidence of hard selective sweep for − 13910*T and − 22018*A in Persians, the largest ethnic group of Iran. The extremely low frequency of − 13915*G in the Iranian population challenged the view that LP distribution in Iran resulted from the demic diffusion, especially mediated by the spread of Islam, from the Arabian Peninsula. Conclusions Our results indicate the distribution of LP in seven ethnic groups across the Iranian plateau. Soft selective sweep rather than hard selective sweep played a substantial role in the evolution of LP in Iranian populations.

Speciation is a process whereby the evolution of reproductive barriers leads to isolated species. Although many studies have addressed large-effect genetic footprints in the advanced stages of speciation, the genetics of reproductive isolation in nascent stage of speciation remains unclear. Here, we show that pig domestication offers an interesting model for studying the early stages of speciation in great details. Pig breeds have not evolved the large X-effect of hybrid incompatibility commonly observed between “good species.” Instead, deleterious epistatic interactions among multiple autosomal loci are common. These weak Dobzhansky–Muller incompatibilities confer partial hybrid inviability with sex biases in crosses between European and East Asian domestic pigs. The genomic incompatibility is enriched in pathways for angiogenesis, androgen receptor signaling and immunity, with an observation of many highly differentiated cis-regulatory variants. Our study suggests that partial hybrid inviability caused by pervasive but weak interactions among autosomal loci may be a hallmark of nascent speciation in mammals.

It has been reported that cyclin-dependent kinase like 3 (CDKL3) plays a crucial role in cell proliferation and migration in several cancers. However, the function of CDKL3 in triple-negative breast cancer (TNBC) is still unclear. In the present study, immunohistochemistry (IHC) was conducted to detect the CDKL3 expression. CCK-8, flow cytometry, Transwell assays, and mice xenograft models, were performed to explore the roles of CDKL3 on the proliferation and migration of TNBC in vitro and in vivo. Besides, protein chip analysis was used to screen the potential pathways, which was further confirmed by promoter activity assay, western blotting, and CCK-8 assay. Our findings reveal a high expression of CDKL3 in TNBC tissues, which is closely related to a poor prognosis of patients with TNBC. In TNBC cells, CDKL3 knockdown inhibits cell proliferation and migration, whereas CDKL3 overexpression has exactly the opposite effect. Consistently, CDKL3 knockdown induces cell apoptosis in vitro but suppresses tumor growth in vivo. Furthermore, CDKL3 knockdown increases p53 expression and reduces cell viability, and these effects are significantly weakened by the p53 inhibitor, PFT-α. In conclusion, the current study highlights that CDKL3 promotes TNBC progressions via regulating the p53 signaling pathway, suggesting that CDKL3 is a novel therapeutic target for TNBC treatment.

The African cattle provide unique genetic resources shaped up by both diverse tropical environmental conditions and human activities, the assessment of their genetic diversity will shade light on the mechanism of their remarkable adaptive capacities. We therefore analyzed the genetic diversity of cattle samples from Nigeria using both maternal and paternal DNA markers. Nigerian cattle can be assigned to 80 haplotypes based on the mitochondrial DNA (mtDNA) D-loop sequences and haplotype diversity was 0.985 + 0.005. The network showed two major matrilineal clustering: the dominant cluster constituting the Nigerian cattle together with other African cattle while the other clustered Eurasian cattle. Paternal analysis indicates only zebu haplogroup in Nigerian cattle with high genetic diversity 1.000 ± 0.016 compared to other cattle. There was no signal of maternal genetic structure in Nigerian cattle population, which may suggest an extensive genetic intermixing within the country. The absence of Bos indicus maternal signal in Nigerian cattle is attributable to vulnerability bottleneck of mtDNA lineages and concordance with the view of male zebu genetic introgression in African cattle. Our study shades light on the current genetic diversity in Nigerian cattle and population history in West Africa.

The domestic dog has evolved to be an important biomedical model for studies regarding the genetic basis of disease, morphology and behavior. Genetic studies in the dog have relied on a draft reference genome of a purebred female boxer dog named “Tasha” initially published in 2005. Derived from a Sanger whole genome shotgun sequencing approach coupled with limited clone-based sequencing, the initial assembly and subsequent updates have served as the predominant resource for canine genetics for 15 years. While the initial assembly produced a good-quality draft, as with all assemblies produced at the time, it contained gaps, assembly errors and missing sequences, particularly in GC-rich regions, which are found at many promoters and in the first exons of protein-coding genes. Here, we present Dog10K_Boxer_Tasha_1.0, an improved chromosome-level highly contiguous genome assembly of Tasha created with long-read technologies that increases sequence contiguity &gt, 100-fold, closes &gt, 23,000 gaps of the CanFam3.1 reference assembly and improves gene annotation by identifying &gt, 1200 new protein-coding transcripts. The assembly and annotation are available at NCBI under the accession GCF_000002285.5.

Background Species domestication is generally characterized by the exploitation of high-impact mutations through processes that involve complex shifting demographics of domesticated species. These include not only inbreeding and artificial selection that may lead to the emergence of evolutionary bottlenecks, but also post-divergence gene flow and introgression. Although domestication potentially affects the occurrence of both desired and undesired mutations, the way wild relatives of domesticated species evolve and how expensive the genetic cost underlying domestication is remain poorly understood. Here, we investigated the demographic history and genetic load of chicken domestication. Results We analyzed a dataset comprising over 800 whole genomes from both indigenous chickens and wild jungle fowls. We show that despite having a higher genetic diversity than their wild counterparts (average π, 0.00326 vs. 0.00316), the red jungle fowls, the present-day domestic chickens experienced a dramatic population size decline during their early domestication. Our analyses suggest that the concomitant bottleneck induced 2.95% more deleterious mutations across chicken genomes compared with red jungle fowls, supporting the “cost of domestication” hypothesis. Particularly, we find that 62.4% of deleterious SNPs in domestic chickens are maintained in heterozygous states and masked as recessive alleles, challenging the power of modern breeding programs to effectively eliminate these genetic loads. Finally, we suggest that positive selection decreases the incidence but increases the frequency of deleterious SNPs in domestic chicken genomes. Conclusion This study reveals a new landscape of demographic history and genomic changes associated with chicken domestication and provides insight into the evolutionary genomic profiles of domesticated animals managed under modern human selection.

The COVID-19 pandemic poses a great threat to human society. SARS-CoV-2 is mainly transmitted through social contact; however, it is highly debated whether cold-chain related transmission has occurred and can be identified in the epidemic areas of COVID-19. Here, we provide a new method and distinguish two transmission routes by detecting a lineage-specific reduction of SARS-CoV-2 mutation rate. After analyzing 1,610,125 SARS-CoV-2 genomic sequences, we find that two outbreaks in Xinfadi-Beijing and Auckland are cold-chain related and respectively caused by two mutation-dormant variants. A Dalian outbreak in July 2020 and a Yingkou outbreak ten months later are epidemiologically connected and derived from a cold-chain related variant. Mutation-dormant variants are detected during the spread of spike D614G variant and the Delta Variant of Concern. Cold-chain contaminations repeatedly caused by epidemiologically connected patients are also found and have resulted in infections. Moreover, the COVID-19 outbreak in Wuhan is likely to be cold-chain related. A systematic identification reveals that the frequency of cold-chain related transmission is in the order of magnitude of 0.1-10%. Our results indicate that that cold-chain related transmission is rare but happens globally.

The rise and expansion of Tibetan Empire in the 7th to 9th centuries AD affected the course of history across East Eurasia, but the genetic impact of Tibetans on surrounding populations remains undefined. We sequenced 60 genomes for four populations from Pakistan and Tajikistan to explore their demographic history. We showed that the genomes of Balti people from Baltistan comprised 22.6–26% Tibetan ancestry. We inferred a single admixture event and dated it to about 39–21 generations ago, a period that postdated the conquest of Baltistan by the ancient Tibetan Empire. The analyses of mitochondrial DNA, Y, and X chromosome data indicated that both ancient Tibetan males and females were involved in the male-biased dispersal. Given the fact that the Balti people adopted Tibetan language and culture in history, our study suggested the impact of Tibetan Empire on Baltistan involved dominant cultural and minor demic diffusion.

Background The reconstruction of phylogenetic relationships for allopolyploids using genomic data is challenging because hybridization and polyploidy can blur history. In cyprinids, one allopolyploidization event involves goldfish and common carp, yet little is known about the origins of other cyprinid polyploid lineages, including their maternal and paternal ancestral lineages. Results Herein we employ 10 HOX genes, genomic and transcriptomic data of representative species from seven subfamilies in Cyprinidae to investigate the origins of polyploid lineages. Analyses of the nuclear and mitochondria genomes reveal that the Schizothoracinae and Cyprininae share the same maternal common ancestor, and identify the candidate genes for identifying maternal gene-copies from Cyprininae. Gene-trees show a close relationship between Sinocyclocheilus grahami and Cyprininae, and indicate that they share the same whole genome duplication (WGD) event about 12.10 ~ 14.03 Ma. Another allopolyploidization event involves Torinae; one duplication clustering with Schizothoracinae and Schizopygopsinae identifies them as paternal siblings. Further, Labeoninae has a history of recurrent hybridization, which is supported by both gene-trees from genomes and HOX genes. Conclusions Collectively, the diverse WGD history makes Cyprinidae a candidate system for investigating the origin of hybridization and polyploidization in vertebrates. Further investigations should enlarge the representatives, combine morphological traits, and reconstruct gene-trees based on whole-genome markers.

The Himalaya are among the youngest and highest mountains in the world, but the exact timing of their uplift and origins of their biodiversity are still in debate. The Himalayan region is a relatively small area but with exceptional diversity and endemism. One common hypothesis to explain the rich montane diversity is uplift-driven diversification—that orogeny creates conditions favoring rapid in situ speciation of resident lineages. We test this hypothesis in the Himalayan region using amphibians and reptiles, two environmentally sensitive vertebrate groups. In addition, analysis of diversification of the herpetofauna provides an independent source of information to test competing geological hypotheses of Himalayan orogenesis. We conclude that the origins of the Himalayan herpetofauna date to the early Paleocene, but that diversification of most groups was concentrated in the Miocene. There was an increase in both rates and modes of diversification during the early to middle Miocene, together with regional interchange (dispersal) between the Himalaya and adjacent regions. Our analyses support a recently proposed stepwise geological model of Himalayan uplift beginning in the Paleocene, with a subsequent rapid increase of uplifting during the Miocene, finally giving rise to the intensification of the modern South Asian Monsoon.

Objective The data provided herein represent the whole-genome resequencing data related to three wolves and three Iranian local dogs. The understanding of genome evolution during animal domestication is an interesting subject in genome biology. Dog is an excellent model for understanding of domestication due to its considerable variety of behavioral and physical traits. The Zagros area of current day Iran has been identified as one of the initial centers of animal domestication. The availability of the complete genome sequences of Iranian local canids can be a valuable resource for researchers to address questions and testing hypotheses on the dog domestication process. Data description We collected blood samples from six Iranian local canids including two hunting dogs (Saluki breed), a mastiff dog (Qahderijani ecotype) and three wolves. We extracted genomic DNA from blood samples. Sequence data were produced using the Illumina HiSeq 2500 system. All sequence data are available in the National Genomics Data Center (NGDC), Genome Sequence Archive (GSA) database under the accession of CRA001324 and the National Center for Biotechnology Information (NCBI) under the accession of PRJNA639312. The short-read sequences with the mean depth of 16X were aligned to the dog reference genome (CanFam3.1) and achieved 99% coverage of the reference assembly. The obtained information from this experiment will be useful in evolutionary biology.

Understanding the mutational and evolutionary dynamics of SARS-CoV-2 is essential for treating COVID-19 and the development of a vaccine. Here, we analyzed publicly available 15,818 assembled SARS-CoV-2 genome sequences, along with 2,350 raw sequence datasets sampled worldwide. We investigated the distribution of inter-host single nucleotide polymorphisms (inter-host SNPs) and intra-host single nucleotide variations (iSNVs). Mutations have been observed at 35.6% (10,649/29,903) of the bases in the genome. The substitution rate in some protein coding regions is higher than the average in SARS-CoV-2 viruses, and the high substitution rate in some regions might be driven to escape immune recognition by diversifying selection. Both recurrent mutations and human-to-human transmission are mechanisms that generate fitness advantageous mutations. Furthermore, the frequency of three mutations (S protein, F400L; ORF3a protein, T164I; and ORF1a protein, Q6383H) has gradual increased over time on lineages, which provides new clues for the early detection of fitness advantageous mutations. Our study provides theoretical support for vaccine development and the optimization of treatment for COVID-19. We call researchers to submit raw sequence data to public databases.

Despite the substantial role that chickens have played in human societies across the world, both the geographic and temporal origins of their domestication remain controversial. To address this issue, we analyzed 863 genomes from a worldwide sampling of chickens and representatives of all four species of wild jungle fowl and each of the five subspecies of red jungle fowl (RJF). Our study suggests that domestic chickens were initially derived from the RJF subspecies Gallus gallus spadiceus whose present-day distribution is predominantly in southwestern China, northern Thailand and Myanmar. Following their domestication, chickens were translocated across Southeast and South Asia where they interbred locally with both RJF subspecies and other jungle fowl species. In addition, our results show that the White Leghorn chicken breed possesses a mosaic of divergent ancestries inherited from other subspecies of RJF. Despite the strong episodic gene flow from geographically divergent lineages of jungle fowls, our analyses show that domestic chickens undergo genetic adaptations that underlie their unique behavioral, morphological and reproductive traits. Our study provides novel insights into the evolutionary history of domestic chickens and a valuable resource to facilitate ongoing genetic and functional investigations of the world’s most numerous domestic animal.

Background The pathophysiological roles of cystic fibrosis transmembrane-conductance regulator (CFTR) Cl− channels in the regulation of blood pressure (BP) remain controversial. Here we studied the function of CFTR Cl− channels in regulation of BP and in the high-fructose-salt-diet (HFSD) induced hypertension in mice. Methods The systolic, diastolic and mean BP (SBP, DBP and MBP, respectively) were continuously monitored from unrestricted conscious wild-type (cftr+/+) FVB and CFTR-knockout (cftr−/−) mice (8-week old, male). HFSD (64.7% fructose, 2% NaCl water) or control normal starch diet (CNSD, 58.9% corn starch, 0 NaCl water) was given for 8 weeks and vascular Doppler were performed. Real-time PCR and Western blot were used to examine mRNA and protein expression, respectively. Results The aortic stiffness, daytime and nighttime SBP, DBP, and MBP of the cftr−/− mice were significantly higher than those in the age- and gender-matched cftr+/+ mice, which is consistent with the findings of increased vascular resistance in cystic fibrosis patients. The aortic stiffness, daytime and nighttime SBP, DBP, and MBP of cftr+/+ mice fed with HFSD were all significantly higher than those fed with CNSD. Importantly, HFSD caused a significant decrease in mRNA and protein expression of WINK1, WINK4 and CFTR in aorta and mesenteric arteries, but not in the kidney, corroborating that HSFD-induced downregulation of WINKs and loss of CFTR function specifically in the arteries may mediate the increased BP. Conclusions CFTR regulates peripheral arterial resistance and BP in vivo. HFSD-induced CFTR downregulation specifically in the arteries may be a novel mechanism for hypertension.

As the only woody resurrection plant, Myrothamnus flabellifolia has a strong tolerance to drought and can survive long-term in a desiccated environment. However, the molecular mechanisms related to the stress tolerance of M. flabellifolia are largely unknown, and few tolerance-related genes previously identified had been functionally characterized. WRKYs are a group of unique and complex plant transcription factors, and have reported functions in diverse biological processes, especially in the regulation of abiotic stress tolerances, in various species. However, little is known about their roles in response to abiotic stresses in M. flabellifolia. In this study, we characterized a dehydration-inducible WRKY transcription factor gene, MfWRKY17, from M. flabellifolia. MfWRKY17 shows high degree of homology with genes from Vitis vinifera and Vitis pseudoreticulata, belonging to group II of the WRKY family. Unlike known WRKY17s in other organisms acting as negative regulators in biotic or abiotic stress responses, overexpression of MfWRKY17 in Arabidopsis significantly increased drought and salt tolerance. Further investigations indicated that MfWRKY17 participated in increasing water retention, maintaining chlorophyll content, and regulating ABA biosynthesis and stress-related gene expression. These results suggest that MfWRKY17 possibly acts as a positive regulator of stress tolerance in the resurrection plant M. flabellifolia.

Genetic introgression not only provides material for adaptive evolution but also confounds our understanding of evolutionary history. This is particularly true for canids, a species complex in which genome sequencing and analysis has revealed a complex history of admixture and introgression. Here, we sequence 19 new whole genomes from high-altitude Tibetan and Himalayan wolves and dogs and combine these into a larger data set of 166 whole canid genomes. Using these data, we explore the evolutionary history and adaptation of these and other canid lineages. We find that Tibetan and Himalayan wolves are closely related to each other, and that ∼39% of their nuclear genome is derived from an as-yet-unrecognized wolf-like lineage that is deeply diverged from living Holarctic wolves and dogs. The EPAS1 haplotype, which is present at high frequencies in Tibetan dog breeds and wolves and confers an adaptive advantage to animals living at high altitudes, was probably derived from this ancient lineage. Our study underscores the complexity of canid evolution and demonstrates how admixture and introgression can shape the evolutionary trajectories of species.

Allotetraploid cyprinids have a unique strategy for balancing subgenomic stabilization and diversification in rediploidization., A persistent enigma is the rarity of polyploidy in animals, compared to its prevalence in plants. Although animal polyploids are thought to experience deleterious genomic chaos during initial polyploidization and subsequent rediploidization processes, this hypothesis has not been tested. We provide an improved reference-quality de novo genome for allotetraploid goldfish whose origin dates to ~15 million years ago. Comprehensive analyses identify changes in subgenomic evolution from asymmetrical oscillation in goldfish and common carp to diverse stabilization and balanced gene expression during continuous rediploidization. The homoeologs are coexpressed in most pathways, and their expression dominance shifts temporally during embryogenesis. Homoeolog expression correlates negatively with alternation of DNA methylation. The results show that allotetraploid cyprinids have a unique strategy for balancing subgenomic stabilization and diversification. Rediploidization process in these fishes provides intriguing insights into genome evolution and function in allopolyploid vertebrates.

Background Yellow-feathered chickens (YFCs) have a long history in China. They are well-known for the nutritional and commercial importance attributable to their yellow color phenotype. Currently, there is a huge paucity in knowledge of the genetic determinants responsible for phenotypic and biochemical properties of these iconic chickens. This study aimed to uncover the genetic structure and the molecular underpinnings of the YFCs trademark coloration. Results The whole-genomes of 100 YFCs from 10 major traditional breeds and 10 Huaibei partridge chickens from China were re-sequenced. Comparative population genomics based on autosomal single nucleotide polymorphisms (SNPs) revealed three geographically based clusters among the YFCs. Compared to other Chinese indigenous chicken genomes incorporated from previous studies, a closer genetic proximity within YFC breeds than between YFC breeds and other chicken populations is evident. Through genome-wide scans for selective sweeps, we identified RALY heterogeneous nuclear ribonucleoprotein (RALY), leucine rich repeat containing G protein-coupled receptor 4 (LGR4), solute carrier family 23 member 2 (SLC23A2), and solute carrier family 2 member 14 (SLC2A14), besides the classical beta-carotene dioxygenase 2 (BCDO2), as major candidates pigment determining genes in the YFCs. Conclusion We provide the first comprehensive genomic data of the YFCs. Our analyses show phylogeographical patterns among the YFCs and potential candidate genes giving rise to the yellow color trait of the YFCs. This study lays the foundation for further research on the genome-phenotype cross-talks that define important poultry traits and for formulating genetic breeding and conservation strategies for the YFCs.

Background Body size is considered as one of the most fundamental properties of an organism. Due to intensive breeding and artificial selection throughout the domestication history, horses exhibit striking variations for heights at withers and body sizes. Debao pony (DBP), a famous Chinese horse, is known for its small body size and lives in Guangxi mountains of southern China. In this study, we employed comparative population genomics to study the genetic basis underlying the small body size of DBP breed based on the whole genome sequencing data. To detect genomic signatures of positive selection, we applied three methods based on population comparison, fixation index (FST), cross population composite likelihood ratio (XP-CLR) and nucleotide diversity (θπ), and further analyzed the results to find genomic regions under selection for body size-related traits. Results A number of protein-coding genes in windows with the top 1% values of FST (367 genes), XP-CLR (681 genes), and log2 (θπ ratio) (332 genes) were identified. The most significant signal of positive selection was mapped to the NELL1 gene, probably underlies the body size and development traits, and may also have been selected for short stature in the DBP population. In addition, some other loci on different chromosomes were identified to be potentially involved in the development of body size. Conclusions Results of our study identified some positively selected genes across the horse genome, which are possibly involved in body size traits. These novel candidate genes may be useful targets for clarifying our understanding of the molecular basis of body size and as such they should be of great interest for future research into the genetic architecture of relevant traits in horse breeding program.

BackgroundPolydactyly is one of the most common congenital limb dysplasia in many animal species. Although preaxial polydactyly (PPD) has been comprehensively studied in humans as a common abnormality, the genetic variations in other animal species have not been fully understood. Herein, we focused on the pig, as an even-toed ungulate mammal model with its unique advantages in medical and genetic researches, two PPD families consisting of four affected and 20 normal individuals were sequenced.ResultsOur results showed that the PPD in the sampled pigs were not related to previously reported variants. A strong association was identified atABCC4and it encodes a transmembrane protein involved in ciliogenesis. We found that the affected and normal individuals were highly differentiated atABCC4, and all the PPD individuals shared long haplotype stretches as compared with the unaffected individuals. A highly differentiated missense mutation (I85T) inABCC4was observed at a residue from a transmembrane domain highly conserved among a variety of organisms.ConclusionsThis study reportsABCC4as a new candidate gene and identifies a missense mutation for PPD in pigs. Our results illustrate a putative role of ciliogenesis process in PPD, coinciding with an earlier observation of ciliogenesis abnormality resulting in pseudo-thumb development in pandas. These results expand our knowledge on the genetic variations underlying PPD in animals.

Piao chicken, a rare Chinese native poultry breed, lacks primary tail structures, such as pygostyle, caudal vertebra, uropygial gland and tail feathers. So far, the molecular mechanisms underlying tail absence in this breed have remained unclear. We employed comprehensive comparative transcriptomic and genomic analyses to unravel potential genetic underpinnings of rumplessness in the Piao chicken. Our results reveal many biological factors involved in tail development and several genomic regions under strong positive selection in this breed. These regions contain candidate genes associated with rumplessness, includingIRX4, IL-18, HSPB2, andCRYAB. Retrieval of quantitative trait loci (QTL) and gene functions implied that rumplessness might be consciously or unconsciously selected along with the high-yield traits in Piao chicken. We hypothesize that strong selection pressures on regulatory elements might lead to gene activity changes in mesenchymal stem cells of the tail bud and eventually result in tail truncation by impeding differentiation and proliferation of the stem cells. Our study provides fundamental insights into early initiation and genetic bases of the rumpless phenotype in Piao chicken.

BackgroundHybridisation and introgression play key roles in the evolutionary history of animal species. They are commonly observed within several orders in wild birds. The domestic chickenGallus gallus domesticusis the most common livestock species. More than 65 billion chickens are raised annually to produce meat and 80 million metric tons of egg for global human consumption by the commercial sector. Unravelling the origin of its genetic diversity has major application for sustainable breeding improvement programmes.ResultsIn this study, we report genome-wide analyses for signatures of introgression between indigenous domestic village chicken and the four wildGallusspecies. We first assess the genome-wide phylogeny and divergence time across the genusGallus. Genome-wide sequence divergence analysis supports a sister relationship between the Grey junglefowlG. sonneratiiand Ceylon junglefowlG. lafayettii. Both species form a clade that is sister to the Red junglefowlG. gallus, with the Green junglefowlG. variusthe most ancient lineage within the genus. We reveal extensive bidirectional introgression between the Grey junglefowl and the domestic chicken and to a much lesser extent with the Ceylon junglefowl. We identify a single case of Green junglefowl introgression. These introgressed regions include genes with biological functions related to development and immune system.ConclusionsOur study shows that while the Red junglefowl is the main ancestral species, introgressive hybridisation episodes have impacted the genome and contributed to the diversity of the domestic chicken, although likely at different levels across its geographic range.