Mid-Pleistocene Transitions Forced Himalayan ibex to Evolve Independently after Split into an Allopatric Refugium.

Pleistocene glaciations had profound impact on the spatial distribution and genetic makeup of species in temperate ecosystems. While the glacial period trapped several species into glacial refugia and caused abrupt decline in large populations, the interglacial period facilitated population growth and range expansion leading to allopatric speciation. Here, we analyzed 40 genomes of four species of ibex and found that Himalayan ibex in the Pamir Mountains evolved independently after splitting from its main range about 0.1 mya following the Pleistocene species pump concept. Demographic trajectories showed Himalayan ibex experienced two historic bottlenecks, one each c. 0.8-0.5 mya and c. 50-30 kya, with an intermediate large population expansion c. 0.2-0.16 mya coinciding with Mid-Pleistocene Transitions. We substantiate with multi-dimensional evidence that Himalayan ibex is an evolutionary distinct phylogenetic species of Siberian ibex which need to be prioritized as Capra himalayensis for taxonomic revision and conservation planning at a regional and global scale.

​Comparative genomics of Balto, a famous historic dog, captures lost diversity of 1920s sled dogs.

We reconstruct the phenotype of Balto, the heroic sled dog renowned for transporting diphtheria antitoxin to Nome, Alaska, in 1925, using evolutionary constraint estimates from the Zoonomia alignment of 240 mammals and 682 genomes from dogs and wolves of the 21st century. Balto shares just part of his diverse ancestry with the eponymous Siberian husky breed. Balto's genotype predicts a combination of coat features atypical for modern sled dog breeds, and a slightly smaller stature. He had enhanced starch digestion compared with Greenland sled dogs and a compendium of derived homozygous coding variants at constrained positions in genes connected to bone and skin development. We propose that Balto's population of origin, which was less inbred and genetically healthier than that of modern breeds, was adapted to the extreme environment of 1920s Alaska.

Pangenome obtained by long-read sequencing of 11 genomes reveal hidden functional structural variants in pigs.

Long-read sequencing (LRS) facilitates both the genome assembly and the discovery of structural variants (SVs). Here, we built a graph-based pig pangenome by incorporating 11 LRS genomes with an average of 94.01% BUSCO completeness score, revealing 206-Mb novel sequences. We discovered 183,352 nonredundant SVs (63% novel), representing 12.12% of the reference genome. By genotyping SVs in an additional 196 short-read sequencing samples, we identified thousands of population stratified SVs. Particularly, we detected 7,568 Tibetan specific SVs, some of which demonstrate significant population differentiation between Tibetan and low-altitude pigs, which might be associated with the high-altitude hypoxia adaptation in Tibetan pigs. Further integrating functional genomic data, the most promising candidate genes within the SVs that might contribute to the high-altitude hypoxia adaptation were discovered. Overall, our study generates a benchmark pangenome resource for illustrating the important roles of SVs in adaptive evolution, domestication, and genetic improvement of agronomic traits in pigs.

Markhor-derived Introgression of a Genomic Region Encompassing PAPSS2 Confers High-altitude Adaptability in Tibetan Goats.

Understanding the genetic mechanism of how animals adapt to extreme conditions is fundamental to determine the relationship between molecular evolution and changing environments. Goat is one of the first domesticated species and has evolved rapidly to adapt to diverse environments, including harsh high-altitude conditions with low temperature and poor oxygen supply but strong ultraviolet radiation. Here, we analyzed 331 genomes of domestic goats and wild caprid species living at varying altitudes (high > 3000 m above sea level and low < 1200 m), along with a reference-guided chromosome-scale assembly (contig-N50: 90.4 Mb) of a female Tibetan goat genome based on PacBio HiFi long reads, to dissect the genetic determinants underlying their adaptation to harsh conditions on the Qinghai-Tibetan Plateau (QTP). Population genomic analyses combined with genome-wide association studies (GWAS) revealed a genomic region harboring the 3'-phosphoadenosine 5'-phosphosulfate synthase 2 (PAPSS2) gene showing strong association with high-altitude adaptability (PGWAS = 3.62 × 10-25) in Tibetan goats. Transcriptomic data from 13 tissues revealed that PAPSS2 was implicated in hypoxia-related pathways in Tibetan goats. We further verified potential functional role of PAPSS2 in response to hypoxia in PAPSS2-deficient cells. Introgression analyses suggested that the PAPSS2 haplotype conferring the high-altitude adaptability in Tibetan goats originated from a recent hybridization between goats and a wild caprid species, the markhor (Capra falconeri). In conclusion, our results uncover a hitherto unknown contribution of PAPSS2 to high-altitude adaptability in Tibetan goats on QTP, following interspecific introgression and natural selection.

Unveiling the functional and evolutionary landscape of RNA editing in chicken using genomics and transcriptomics.

The evolutionary and functional features of RNA editing are well studied in mammals, cephalopods, and insects, but not in birds. Here, we integrated transcriptomic and whole-genomic analyses to exhaustively characterize the expansive repertoire of adenosine-to-inosine (A-to-I) RNA editing sites (RESs) in the chicken. In addition, we investigated the evolutionary status of the chicken editome as a potential mechanism of domestication. We detected the lowest editing level in the liver of chickens, compared to muscles in humans, and found higher editing activity and specificity in the brain than in non-neural tissues, consistent with the brain's functional complexity. To a certain extent, specific editing activity may account for the specific functions of tissues. Our results also revealed that sequences critical to RES secondary structures remained conserved within avian evolution. Furthermore, the RNA editome was shaped by purifying selection during chicken domestication and most RESs may have served as a selection pool for a few functional RESs involved in chicken domestication, including evolution of nervous and immune systems. Regulation of RNA editing in chickens by adenosine deaminase acting on RNA (ADAR) enzymes may be affected by non-ADAR factors whose expression levels changed widely after ADAR knockdown. Collectively, we provide comprehensive lists of candidate RESs and non-ADAR-editing regulators in the chicken, thus contributing to our current understanding of the functions and evolution of RNA editing in animals.

The Genetic Echo of the Tarim Mummies in Modern Central Asians.

The diversity of Central Asians has been shaped by multiple migrations and cultural diffusion. Although ancient DNA studies have revealed the demographic changes of the Central Asian since the Bronze Age, the contribution of the ancient populations to the modern Central Asian remains opaque. Herein, we performed high-coverage sequencing of 131 whole genomes of Indo-European-speaking Tajik and Turkic-speaking Kyrgyz populations to explore their genomic diversity and admixture history. By integrating the ancient DNA data, we revealed more details of the origins and admixture history of Central Asians. We found that the major ancestry of present-day Tajik populations can be traced back to the admixture of the Bronze Age Bactria-Margiana Archaeological Complex and Andronovo-related populations. Highland Tajik populations further received additional gene flow from the Tarim mummies, an isolated ancient North Eurasian-related population. The West Eurasian ancestry of Kyrgyz is mainly derived from Historical Era populations in Xinjiang of China. Furthermore, the recent admixture signals detected in both Tajik and Kyrgyz are ascribed to the expansions of Eastern Steppe nomadic pastoralists during the Historical Era.

A polar bear paleogenome reveals extensive ancient gene flow from polar bears into brown bears.

Polar bears (Ursus maritimus) and brown bears (Ursus arctos) are sister species possessing distinct physiological and behavioural adaptations that evolved over the last 500,000 years. However, comparative and population genomics analyses have revealed that several extant and extinct brown bear populations have relatively recent polar bear ancestry, probably as the result of geographically localized instances of gene flow from polar bears into brown bears. Here, we generate and analyse an approximate 20X paleogenome from an approximately 100,000-year-old polar bear that reveals a massive prehistoric admixture event, which is evident in the genomes of all living brown bears. This ancient admixture event was not visible from genomic data derived from living polar bears. Like more recent events, this massive admixture event mainly involved unidirectional gene flow from polar bears into brown bears and occurred as climate changes caused overlap in the ranges of the two species. These findings highlight the complex reticulate paths that evolution can take within a regime of radically shifting climate.

Genome-wide scan for selection signatures and genes related to heat tolerance in domestic chickens in the tropical and temperate regions in Asia.

Heat stress is one of the major environmental stressors challenging the global poultry industry. Identifying the genes responsible for heat tolerance is fundamentally important for direct breeding programs. To uncover the genetic basis underlying the ambient temperature adaptation of chickens, we analyzed a total of 59 whole genomes from indigenous chickens that inhabit South Asian tropical regions and temperate regions from Northern China. We applied FST and π-ratio to scan selective sweeps and identified 34 genes with a signature of positive selection in chickens from tropical regions. Several of these genes are functionally implicated in metabolism (FABP2, RAMP3, SUGCT, and TSHR) and vascular smooth muscle contractility (CAMK2), and they may be associated with adaptation to tropical regions. In particular, we found a missense mutation in thyroid-stimulating hormone receptor (41020238:G>A) that shows significant differences in allele frequency between the chicken populations of the two regions. To evaluate whether the missense mutation in TSHR could enhance the heat tolerance of chickens, we constructed segregated chicken populations and conducted heat stress experiments using homozygous mutations (AA) and wild-type (GG) chickens. We found that GG chickens exhibited significantly higher concentrations of alanine aminotransferase, lactate dehydrogenase, and creatine kinase than AA chickens under heat stress (35 ± 1°C) conditions (P < 0.05). These results suggest that TSHR (41020238:G>A) can facilitate heat tolerance and adaptation to higher ambient temperature conditions in tropical climates. Overall, our results provide potential candidate genes for molecular breeding of heat-tolerant chickens.

Genomic basis for skin phenotype and cold adaptation in the extinct Steller's sea cow.

Steller's sea cow, an extinct sirenian and one of the largest Quaternary mammals, was described by Georg Steller in 1741 and eradicated by humans within 27 years. Here, we complement Steller's descriptions with paleogenomic data from 12 individuals. We identified convergent evolution between Steller's sea cow and cetaceans but not extant sirenians, suggesting a role of several genes in adaptation to cold aquatic (or marine) environments. Among these are inactivations of lipoxygenase genes, which in humans and mouse models cause ichthyosis, a skin disease characterized by a thick, hyperkeratotic epidermis that recapitulates Steller's sea cows' reportedly bark-like skin. We also found that Steller's sea cows' abundance was continuously declining for tens of thousands of years before their description, implying that environmental changes also contributed to their extinction.

Genome Sequencing of a Gray Wolf from Peninsular India Provides New Insights into the Evolution and Hybridization of Gray Wolves.

The gray wolf (Canis lupus) is among the few large carnivores that survived the Late Pleistocene megafaunal extinctions. Thanks to their complex history of admixture and extensive geographic range, the number of gray wolf subspecies and their phylogenetic relationships remain poorly understood. Here, we perform whole-genome sequencing of a gray wolf collected from peninsular India that was phenotypically distinct from gray wolves outside India. Genomic analyses reveal that the Indian gray wolf is an evolutionarily distinct lineage that diverged from other extant gray wolf lineages ∼110 thousand years ago. Demographic analyses suggest that the Indian wolf population declined continuously decline since separating from other gray wolves and, today, has exceptionally low genetic diversity. We also find evidence for pervasive and mosaic gene flow between the Indian wolf and African canids including African wolf, Ethiopian wolf, and African wild dog despite their current geographical separation. Our results support the hypothesis that the Indian subcontinent was a Pleistocene refugium and center of diversification and further highlight the complex history of gene flow that characterized the evolution of gray wolves.

NOD1 Is Associated With the Susceptibility of Pekin Duck Flock to Duck Hepatitis A Virus Genotype 3.

Duck viral hepatitis (DVH) is an acute, highly lethal infectious disease of ducklings that causes huge losses in the duck industry. Duck hepatitis A virus genotype 3 (DHAV-3) has been one of the most prevalent DVH pathogen in the Asian duck industry in recent years. Here, we investigated the genetic basis of the resistance and susceptibility of ducks to DVH by comparing the genomes and transcriptomes of a resistant Pekin duck flock (Z8) and a susceptible Pekin duck flock (SZ7). Our comparative genomic and transcriptomic analyses suggested that NOD1 showed a strong signal of association with DVH susceptibility in ducks. Then, we found that NOD1 showed a significant expression difference between the livers of susceptible and resistant individuals after infection with DHAV-3, with higher expression in the SZ7 flock. Furthermore, suppression and overexpression experiments showed that the number of DHAV-3 genomic copies in primary duck hepatocytes was influenced by the expression level of NOD1. In addition, in situ RNAscope analysis showed that the localization of NOD1 and DHAV-3 in liver cells was consistent. Altogether, our data suggested that NOD1 was likely associated with DHAV-3 susceptibility in ducks, which provides a target for future investigations of the pathogenesis of DVH.

Large-scale genomic analysis reveals the genetic cost of chicken domestication.

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.

Genomes reveal selective sweeps in kiang and donkey for high-altitude adaptation.

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.

Intraoperative Oxygen Concentration and Postoperative Delirium After Laparoscopic Gastric and Colorectal Malignancies Surgery: A Randomized, Double-Blind, Controlled Trial.

PURPOSE: Postoperative delirium (POD) is common in elderly patients undergoing laparoscopic surgery for gastric and colorectal malignancies. POD may be affected by different fraction of inspired oxygen (FiO2). The purpose of this study was to compare the effects of different FiO2 on POD. PATIENTS AND METHODS: A randomized, double-blind controlled trial was performed in Qingdao Municipal Hospital Affiliated to Qingdao University. A total of 662 patients aged 65 to 85 years old underwent isolated laparoscopic radical gastrectomy, radical resection of colon cancer, or radical resection of rectal cancer only. A random number table method was used to divide the patients into two groups: 40% FiO2 (group A) and 80% FiO2 (group B). The primary endpoint was the incidence of POD, which was assessed by the Confusion Assessment Method (CAM) twice daily during the first 7 postoperative days, and POD severity was measured by the Memorial Delirium Assessment Scale (MDAS). The secondary endpoints were the intraoperative regional cerebral oxygen saturation (rSO2), Bispectral (BIS) index, invasive arterial blood pressure (IABP), oxygen saturation (SpO2), end-tidal carbon dioxide partial pressure (PETCO2), the number of atelectasis cases and visual analogue scale (VAS) scores on days 1-7 after surgery. RESULTS: The incidence of POD was 19.37% (122/630), including 20.38% (64/314) in group A and 18.35% (58/316) in group B. No statistical significance was found in the incidence of POD between the two groups (P > 0.05); compared with group B, SpO2, rSO2 and PaO2 decreased at T2 to T4 time point (P < 0.01), and the incidence of postoperative atelectasis decreased (P < 0.05) in group A. CONCLUSION: The incidence of POD was not significantly affected by different FiO2 and the incidence of postoperative atelectasis was decreased at low FiO2.

Multicenter analysis and a rapid screening model to predict early novel coronavirus pneumonia using a random forest algorithm.

ABSTRACT: Early determination of coronavirus disease 2019 (COVID-19) pneumonia from numerous suspected cases is critical for the early isolation and treatment of patients.The purpose of the study was to develop and validate a rapid screening model to predict early COVID-19 pneumonia from suspected cases using a random forest algorithm in China.A total of 914 initially suspected COVID-19 pneumonia in multiple centers were prospectively included. The computer-assisted embedding method was used to screen the variables. The random forest algorithm was adopted to build a rapid screening model based on the training set. The screening model was evaluated by the confusion matrix and receiver operating characteristic (ROC) analysis in the validation.The rapid screening model was set up based on 4 epidemiological features, 3 clinical manifestations, decreased white blood cell count and lymphocytes, and imaging changes on chest X-ray or computed tomography. The area under the ROC curve was 0.956, and the model had a sensitivity of 83.82% and a specificity of 89.57%. The confusion matrix revealed that the prospective screening model had an accuracy of 87.0% for predicting early COVID-19 pneumonia.Here, we developed and validated a rapid screening model that could predict early COVID-19 pneumonia with high sensitivity and specificity. The use of this model to screen for COVID-19 pneumonia have epidemiological and clinical significance.

Genomic Analyses of Unveil Helmeted Guinea Fowl (Numida meleagris) Domestication in West Africa.

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.

Integrating Genomic and Transcriptomic Data to Reveal Genetic Mechanisms Underlying Piao Chicken Rumpless Trait.

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 remain unclear. In this study, we comprehensively employed comparative transcriptomic and genomic analyses to unravel potential genetic underpinnings of rumplessness in 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, including Irx4, Il18, Hspb2, and Cryab. Retrieval of quantitative trait loci (QTL) and gene functions implies 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 changes in gene activity in mesenchymal stem cells of the tail bud. The ectopic activity could eventually result in tail truncation by impeding differentiation and proliferation of the stem cells. Our study provides fundamental insights into early initiation and genetic basis of the rumpless phenotype in Piao chicken.

A rapid screening model for early predicting novel coronavirus pneumonia in Zhejiang Province of China: a multicenter study.

Novel coronavirus pneumonia (NCP) has been widely spread in China and several other countries. Early finding of this pneumonia from huge numbers of suspects gives clinicians a big challenge. The aim of the study was to develop a rapid screening model for early predicting NCP in a Zhejiang population, as well as its utility in other areas. A total of 880 participants who were initially suspected of NCP from January 17 to February 19 were included. Potential predictors were selected via stepwise logistic regression analysis. The model was established based on epidemiological features, clinical manifestations, white blood cell count, and pulmonary imaging changes, with the area under receiver operating characteristic (AUROC) curve of 0.920. At a cut-off value of 1.0, the model could determine NCP with a sensitivity of 85% and a specificity of 82.3%. We further developed a simplified model by combining the geographical regions and rounding the coefficients, with the AUROC of 0.909, as well as a model without epidemiological factors with the AUROC of 0.859. The study demonstrated that the screening model was a helpful and cost-effective tool for early predicting NCP and had great clinical significance given the high activity of NCP.

Revisiting the evolutionary history of domestic and wild ducks based on genomic analyses.

Although domestic ducks have been important poultry species throughout human history, their origin remains enigmatic, with mallards and/or Chinese spot-billed ducks being proposed as the direct wild ancestor(s) of domestic ducks. Here, we analyzed 118 whole genomes from mallard, Chinese spot-billed, and domestic ducks to reconstruct their evolutionary history. We found pervasive introgression patterns among these duck populations. Furthermore, we showed that domestic ducks separated from mallard and Chinese spot-billed ducks nearly 38 thousand years ago (kya) and 54 kya, respectively, which is considerably outside the time period of presumed duck domestication. Thus, our results suggest that domestic ducks may have originated from another wild duck population that is currently undefined or unsampled, rather than from present-day mallard and/or Chinese spot-billed ducks, as previously thought. Overall, this study provides new insight into the complex evolution of ducks.