The Porcine and Chicken Innate DNA Sensing cGAS-STING-IRF Signaling Axes Exhibit Differential Species Specificity.

The innate immune DNA sensing cyclic GMP-AMP synthase (cGAS)-stimulator of IFN genes (STING) signaling pathway plays a key role in host antiviral function. Although the cGAS-STING pathway has been extensively studied, the cGAS-STING signaling in livestock and poultry is not well understood, and whether the species specificity exists is still unknown. In this study, we found that porcine and chicken STING, but not cGAS, exhibit species differences in regulation of IFN; that is, porcine (p)STING mediates good induction of IFN in mammalian cells and low IFN induction in chicken DF-1 cells; on the contrary, chicken (ch)STING mediates IFN induction only in chicken cells but not in mammalian cells. Furthermore, it was found that the motifs pLxIS of pSTING and pLxVS of chSTING are responsible for the species disparity, with the IFN activity of pSTING and chSTING exchanged by swapping the two pLxI/VS motifs. The pLxI/VS motifs mediated the interactions of various STING with downstream IFN regulatory factors (IRFs), reflecting the species-specific pIRF3 and chIRF7. Next, the STING, IRFs, and STING-IRFs were reconstituted in porcine and chicken macrophages that were genetically knocked out for STING and/or IRFs by the CRISPR-Cas9 approach. The results showed that pSTING plus pIRF3 or chIRF7 are able to induce IFN; however, chSTING plus chIRF7 but not pIRF3 are able to induce IFN, suggesting that pIRF3 is specific and stringent, which underlies the inability of chSTING to induce IFN in mammalian cells. In summary, our findings reveal the differential species specificity in the cGAS-STING pathway and the underlying mechanisms, thus providing valuable insights on the cGAS-STING-IRF signaling axis for comparative immunology.

Identification of a Monoclonal Antibody against Porcine Deltacoronavirus Membrane Protein.

Porcine deltacoronavirus (PDCoV) is an emerging virus that poses a significant threat to the global swine industry. Its membrane (M) protein is crucial for virion assembly and virus-host interactions. We selected the hydrophilic region of M protein for prokaryotic expression, purification, and recombinant protein production. Utilizing hybridoma technology, we prepared the monoclonal antibody (mAb) 24-A6 against M protein. The mAb 24-A6 was shown to be suitable for use in immunofluorescence assays, western blotting, and immunoprecipitation, with specificity for PDCoV and no cross-reactivity with other five porcine viruses. The M protein was observed to be expressed as early as 3 h after PDCoV infection, increasing its expression over the duration of infection. Notably, the antigenic epitope of the M protein identified as 103SPESRL108 recognized by mAb 24-A6 was found within a conserved structural domain (SWWSFNPETNNL) of the coronavirus M protein, indicating a crucial overlap between a functionally important viral assembly region and a region recognized by the immune system. Our findings provide valuable insights into mAb 24-A6 targeting the antigenic epitope of M protein and may contribute to the development of diagnostic tools for PDCoV infection and fundamental research into the function of PDCoV M protein.

ncRNAs regulate bovine adipose tissue deposition.

Lipid metabolism, which encompasses synthesis and degradation of lipids, is critical for a wide range of cellular functions, including structural and morphological properties of organelles, energy storage, signalling, and the stability and function of membrane proteins. Adipose tissue is a dynamic tissue type that performs a lot of significant physiological functions, including secretion, and is involved in maintaining homeostasis and in regulatory roles of other tissues based on paracrine or endocrine. More recently, several classes of non-coding RNAs (ncRNAs), such as long non-coding RNA (lncRNA), microRNA (miRNA) and circular RNA (circRNA), have been discovered in adipocytes, and they act as critical regulators of gene expression in adipogenesis and regulate adipogenesis through multiple pathways. In the present paper, we discussed several classes of non-coding RNAs and summarized the latest research on the regulatory role of ncRNAs in bovine adipogenesis. We gave examples for known modes of action to look forward to providing reference information future scientific research in cattle breeding.

Evolution of Transmissible Gastroenteritis Virus (TGEV): A Codon Usage Perspective.

Transmissible gastroenteritis virus (TGEV) is a coronavirus associated with diarrhea and high mortality in piglets. To gain insight into the evolution and adaptation of TGEV, a comprehensive analysis of phylogeny and codon usage bias was performed. The phylogenetic analyses of maximum likelihood and Bayesian inference displayed two distinct genotypes: genotypes I and II, and genotype I was classified into subtypes Ia and Ib. The compositional properties revealed that the coding sequence contained a higher number of A/U nucleotides than G/C nucleotides, and that the synonymous codon third position was A/U-enriched. The principal component analysis based on the values of relative synonymous codon usage (RSCU) showed the genotype-specific codon usage patterns. The effective number of codons (ENC) indicated moderate codon usage bias in the TGEV genome. Dinucleotide analysis showed that CpA and UpG were over-represented and CpG was under-represented in the coding sequence of the TGEV genome. The analyses of Parity Rule 2 plot, ENC-plot, and neutrality plot displayed that natural selection was the dominant evolutionary driving force in shaping codon usage preference in genotypes Ia and II. In addition, natural selection played a major role, while mutation pressure had a minor role in driving the codon usage bias in genotype Ib. The codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) analyses suggested that genotype I might be more adaptive to pigs than genotype II. Current findings contribute to understanding the evolution and adaptation of TGEV.

A rare novel mutation in TECTA causes autosomal dominant nonsyndromic hearing loss in a Mongolian family.

BACKGROUND: The genetic basis of autosomal dominant nonsyndromic hearing loss is complex. Genetic factors are responsible for approximately 50% of cases with congenital hearing loss. However, no previous studies have documented the clinical phenotype and genetic basis of autosomal dominant nonsyndromic hearing loss in Mongolians. METHODS: In this study, we performed exon capture sequencing of a Mongolian family with hereditary hearing loss and identified a novel mutation in TECTA gene, which encodes α -tectorin, a major component of the inner ear extracellular matrix that contacts the specialized sensory hair cells. RESULTS: The novel G → T missense mutation at nucleotide 6016 results in a substitution of amino acid aspartate at 2006 with tyrosine (Asp2006Tyr) in a highly conserved zona pellucida (ZP) domain of α-tectorin. The mutation is not found in control subjects from the same family with normal hearing and a genotype-phenotype correlation is observed. CONCLUSION: A novel missense mutation c.6016 G > T (p.Asp2006Tyr) of TECTA gene is a characteristic TECTA-related mutation which causes autosomal dominant nonsyndromic hearing loss. Our result indicated that mutation in TECTA gene is responsible for the hearing loss in this Mongolian family.

The Evolution and Global Spatiotemporal Dynamics of Senecavirus A.

Recurrent outbreaks of senecavirus A (SVA)-associated vesicular disease have led to a large number of infected pigs being culled and has caused considerable economic losses to the swine industry. Although SVA was discovered 2 decades ago, knowledge about the evolutionary and transmission histories of SVA remains unclear. Herein, we performed an integrated analysis of the recombination, phylogeny, selection, and spatiotemporal dynamics of SVA. Phylogenetic analysis demonstrated that SVA diverged into two main branches, clade I (pre-2007 strains) and clade II (post-2007 strains). Importantly, analysis of selective strength showed that clade II was evolving under relaxed selection compared with clade I. Positive selection analysis identified 27 positive selective sites, most of which are located on the outer surface of capsid protomer or on the important functional domains of nonstructure proteins. Bayesian phylodynamics suggested that the estimated time to the most recent common ancestor of SVA was around 1986, and the estimated substitution rate of SVA was 3.3522 × 10-3 nucleotide substitutions/site/year. Demographic history analysis revealed that the effective population size of SVA has experienced a gradually increasing trend with slight fluctuation until 2017 followed by a sharp decline. Notably, Bayesian phylogeographic analysis inferred that Brazil might be the source of SVA's global transmission since 2015. In summary, these data illustrated that the ongoing evolution of SVA drove the lineage-specific innovation and potentially phenotypically important variation. Our study sheds new light on the fundamental understanding of SVA evolution and spread history. IMPORTANCE Recurrent outbreaks and global epidemics of senecavirus A-associated vesicular disease have caused heavy economic losses and have threatened the development of the pig industry. However, the question of where the virus comes from has been one of the biggest puzzles due to the stealthy nature of the virus. Consequently, tracing the source, evolution, and transmission pattern of SVA is a very challenging task. Based on the most comprehensive analysis, we revealed the origin time, rapid evolution, epidemic dynamics, and selection of SVA. We observed two main genetic branches, clade I (pre-2007 strains) and clade II (post-2007 strains), and described the epidemiological patterns of SVA in different countries. We also first identified Brazil as the source of SVA's global transmission since 2015. Findings in this study provide important implications for the control and prevention of the virus.

Emerging of H5N6 Subtype Influenza Virus with 129-Glycosylation Site on Hemagglutinin in Poultry in China Acquires Immune Pressure Adaption.

For an investigation into the effects of glycosylation site modification on hemagglutinin (HA) on the biological characteristics of the H5N6 subtype avian influenza virus (AIV), the HA sequences of H5N6 AIVs from Global Initiative on Sharing All Influenza Data (GISAID) and the isolates in China were analyzed for genetic evolution and glycosylation site patterns. Eight recombinant H5N6 AIVs with different glycosylation site patterns were constructed, and their biological characteristics were determined. The results showed that H5N6 AIVs containing a 129-glycosylation site on HA are becoming prevalent strains in China. Acquisition of the 129-glycosylation site on the HA of H5N6 AIVs increased thermostability, decreased pH stability, and attenuated pathogenicity and contact transmission in chickens. Most importantly, H5N6 AIVs escaped the neutralization activity of the Re-8-like serum antibody. Our findings reveal that H5N6 AIVs containing the 129-glycosylation site affect antigenicity and have become prevalent strains in China. IMPORTANCE H5N6 avian influenza viruses (AIVs) were first reported in 2013 and have spread throughout many countries. In China, compulsory vaccine inoculation has been adopted to control H5 subtype avian influenza. However, the effect of vaccination on the antigenic drift of H5N6 AIVs remains unknown. Here, we found that H5N6 AIVs with the 129-glycosylation site on hemagglutinin were the dominant strains in poultry in China. The neutralization assay of the serum antibody against the H5 subtype vaccine Re-8 showed a significantly lower neutralization activity against H5N6 AIVs with the 129-glycosylation site compared to that against H5N6 AIVs without the 129-glycosylation site, indicating that the 129-glycosylation site may be a crucial molecular marker for immune evasion.

Genome-Wide Reassortment Analysis of Influenza A H7N9 Viruses Circulating in China during 2013-2019.

Reassortment with the H9N2 virus gave rise to the zoonotic H7N9 avian influenza virus (AIV), which caused more than five outbreak waves in humans, with high mortality. The frequent exchange of genomic segments between H7N9 and H9N2 has been well-documented. However, the reassortment patterns have not been described and are not yet fully understood. Here, we used phylogenetic analyses to investigate the patterns of intersubtype and intrasubtype/intralineage reassortment across the eight viral segments. The H7N9 virus and its progeny frequently exchanged internal genes with the H9N2 virus but rarely with the other AIV subtypes. Before beginning the intrasubtype/intralineage reassortment analyses, five Yangtze River Delta (YRD A-E) and two Pearl River Delta (PRD A-B) clusters were divided according to the HA gene phylogeny. The seven reset segment genes were also nomenclatured consistently. As revealed by the tanglegram results, high intralineage reassortment rates were determined in waves 2-3 and 5. Additionally, the clusters of PB2 c05 and M c02 were the most dominant in wave 5, which could have contributed to the onset of the largest H7N9 outbreak in 2016-2017. Meanwhile, a portion of the YRD-C cluster (HP H7N9) inherited their PB2, PA, and M segments from the co-circulating YRD-E (LP H7N9) cluster during wave 5. Untanglegram results revealed that the reassortment rate between HA and NA was lower than HA with any of the other six segments. A multidimensional scaling plot revealed a robust genetic linkage between the PB2 and PA genes, indicating that they may share a co-evolutionary history. Furthermore, we observed relatively more robust positive selection pressure on HA, NA, M2, and NS1 proteins. Our findings demonstrate that frequent reassortment, particular reassorted patterns, and adaptive mutations shaped the H7N9 viral genetic diversity and evolution. Increased surveillance is required immediately to better understand the current state of the HP H7N9 AIV.

Disruption of interferon-ß production by the Npro of atypical porcine pestivirus.

Atypical porcine pestivirus (APPV) is an emerging porcine virus that threatens global swine production. Pestiviruses can prevent interferon (IFN) production to avoid the host innate immune response, and the Npro viral protein can play a critical role. Knowledge of the host immune response to APPV infection is limited. Here, we showed that the IFN-ß production was suppressed by APPV-Npro and the IFN regulatory factor 3 (IRF3) promoter activity stimulated by adaptor molecules of the IFN-ß signaling pathway was also inhibited in the APPV-Npro-expressed cells. The APPV-Npro was able to interact with IRF3 and interfere the phosphorylation of IRF3, indicated that the IFN-ß antagonism of APPV-Npro mainly depended on blocking IRF3 activity. To identify the functional region of APPV-Npro, a panel of truncated APPV-Npro was constructed, and its influence on the IRF3 activation was investigated. The results showed that the N-terminal 31-51 amino acids of APPV-Npro were mainly associated with inhibition of the IFN-ß response. Taken together, this is the first study focusing on elucidating the function of APPV protein by revealing a novel mechanism of Npro in disruption of host IFN-ß production, which will enlighten future study in addressing APPV pathogenesis and immune evasion.

Genetic and antigenic diversity of H7N9 highly pathogenic avian influenza virus in China.

Avian influenza virus (AIV) H7N9 that emerged in 2013 in eastern China is a novel zoonotic agent mainly circulating in poultry without clinical signs but causing severe disease with high fatality in humans in more than 5 waves. Since the emergence of highly pathogenic (HP) H7N9 variants in 2016, it has induced heavy losses in the poultry industry leading to the implementation of an intensive nationwide vaccination program at the end of wave 5 (September 2017). To characterize the ongoing evolution of H7N9 AIV, we conducted analyses of H7N9 glycoprotein genes obtained from 2013 to 2019. Bayesian analyses revealed a decreasing population size of HP H7N9 variants post wave 5. Phylogenetic topologies revealed that two novel small subclades were formed and carried several fixed amino acid mutations that were along HA and NA phylogenetic trees since wave 5. Some of the mutations were located at antigenic sites or receptor binding sites. The antigenic analysis may reveal a significant antigenic drift evaluated by hemagglutinin inhibition (HI) assay and the antigenicity of H7N9 AIV might evolute in large leaps in wave 7. Molecular simulations found that the mutations (V135T, S145P, and L226Q) around the HA receptor pocket increased the affinity to α2,3-linked sialic acid (SIA) while decreased to α2,6-linked SIA. Altered affinity may suggest that HP H7N9 variations aggravate the pathogenicity to poultry but lessen the threat to public health. Selection analyses showed that the HP H7N9 AIV experienced an increasing selection pressure since wave 5, and the national implementation of vaccination might intensify the role of natural selection during the evolution waves 6 and 7. In summary, our data provide important insights about the genetic and antigenic diversity of circulating HP H7N9 viruses from 2017 to 2019. Enhanced surveillance is urgently warranted to understand the current situation of HP H7N9 AIV.

Spatiotemporal Associations and Molecular Evolution of Highly Pathogenic Avian Influenza A H7N9 Virus in China from 2017 to 2021.

Highly pathogenic (HP) H7N9 avian influenza virus (AIV) emerged in China in 2016. HP H7N9 AIV caused at least 33 human infections and has been circulating in poultry farms continuously since wave 5. The genetic divergence, geographic patterns, and hemagglutinin adaptive and parallel molecular evolution of HP H7N9 AIV in China since 2017 are still unclear. Here, 10 new strains of HP H7N9 AIVs from October 2019 to April 2021 were sequenced. We found that HP H7N9 was primarily circulating in Northern China, particularly in the provinces surrounding the Bohai Sea (Liaoning, Hebei, and Shandong) since wave 6. Of note, HP H7N9 AIV phylogenies exhibit a geographical structure compatible with high levels of local transmission after unidirectional rapid geographical expansion towards the north of China in 2017. In addition, we showed that two major subclades were continually expanding with the viral population size undergoing a sharp increase after 2018 with an obvious seasonal tendency. Notably, the hemagglutinin gene showed signs of parallel evolution and positive selection. Our research sheds light on the current epidemiology, evolution, and diversity of HP H7N9 AIV that can help prevent and control the spreading of HP H7N9 AIV.

Integrative Analysis Revealing Human Heart-Specific Genes and Consolidating Heart-Related Phenotypes.

Elucidating expression patterns of heart-specific genes is crucial for understanding developmental, physiological, and pathological processes of the heart. The aim of the present study is to identify functionally and pathologically important heart-specific genes by performing the Ingenuity Pathway Analysis (IPA). Through a median-based analysis of tissue-specific gene expression based on the Genotype-Tissue Expression (GTEx) data, we identified 56 genes with heart-specific or elevated expressions in the heart (heart-specific/enhanced), among which three common heart-specific/enhanced genes and four atrial appendage-specific/enhanced genes were unreported regarding the heart. Differential expression analysis further revealed 225 differentially expressed genes (DEGs) between atrial appendage and left ventricle. Our integrative analyses of those heart-specific/enhanced genes and DEGs with IPA revealed enriched heart-related traits and diseases, consolidating evidence of relationships between these genes and heart function. Our reports on comprehensive identification of heart-specific/enhanced genes and DEGs and their relation to pathways associated with heart-related traits and diseases provided molecular insights into essential regulators of cardiac physiology and pathophysiology and potential new therapeutic targets for heart diseases.

Phylogenetic and codon usage analysis of atypical porcine pestivirus (APPV).

Atypical porcine pestivirus (APPV) has been identified as the main causative agent for congenital tremor (CT) type A-II in piglets, which is threatening the health of the global swine herd. However, the evolution of APPV remains largely unknown. In this study, phylogenetic analysis showed that APPV could be divided into three phylogroups (I, II, and III). Phylogroups I and II included viral strains from China, while phylogroup III contained strains from Europe, North America, and Asia. Phylogroups I and II are tentatively thought to be of Chinese origin. Next, compositional property analysis revealed that a high frequency of nucleotide A and A-end codons was used in the APPV genome. Intriguingly, the analysis of preferred codons revealed that the AGA[Arg] and AGG[Arg] were overrepresented. Dinucleotide CC was found to be overrepresented, and dinucleotide CG was underrepresented. Furthermore, it was found that the weak codon usage bias of APPV was mainly dominated by selection pressures versus mutational forces. The codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) analyses showed that the codon usage patterns of phylogroup II and III were more similar to the one of a pig than phylogroup I, suggesting that phylogroup II and III may be more adaptive to pigs. Overall, this study provides insights into APPV evolution through phylogeny and codon usage pattern analysis.

Comprehensive Analysis of Codon Usage on Porcine Astrovirus.

Porcine astrovirus (PAstV), associated with mild diarrhea and neurological disease, is transmitted in pig farms worldwide. The purpose of this study is to elucidate the main factors affecting codon usage to PAstVs. Phylogenetic analysis showed that the subtype PAstV-5 sat at the bottom of phylogenetic tree, followed by PAstV-3, PAstV-1, PAstV-2, and PAstV-4, indicating that the five existing subtypes (PAstV1-PAstV5) may be formed by multiple differentiations of PAstV ancestors. A codon usage bias was found in the PAstVs-2,3,4,5 from the analyses of effective number of codons (ENC) and relative synonymous codon usage (RSCU). Nucleotides A/U are more frequently used than nucleotides C/G in the genome CDSs of the PAstVs-3,4,5. Codon usage patterns of PAstV-5 are dominated by mutation pressure and natural selection, while natural selection is the main evolutionary force that affects the codon usage pattern of PAstVs-2,3,4. The analyses of codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) showed the codon usage similarities between the PAstV and animals might contribute to the broad host range and the cross-species transmission of astrovirus. Our results provide insight into understanding the PAstV evolution and codon usage patterns.

Identification of a Novel Imprinted Transcript in the Porcine GNAS Complex Locus Using Methylome and Transcriptome of Parthenogenetic Fetuses.

Genomic imprinting in domestic animals contributes to the variance of performance traits. However, research remains to be done on large-scale detection of epigenetic landscape of porcine imprinted loci including the GNAS complex locus. The purpose of this study was to generate porcine parthenogenetic fetuses and comprehensively identify imprinting patterns of the GNAS locus in transcript levels. To this end, both normally fertilized and bimaternal (uniparental) parthenogenetic porcine fetuses were generated, and whole genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) were performed to construct methylome and transcriptome, respectively. Differentially methylated regions (DMRs) between the fetuses were identified through methylome analysis, and parental-origin-specific expression patterns of transcripts were examined with transcriptome. As a result, three major DMRs were identified: paternally methylated Nesp DMR, maternally methylated Nespas-Gnasxl DMR, and maternally methylated Exon1B-Exon1A DMR. Parental-origin-specific expressions of those five DMR-affected transcripts were found, including a novel imprinted transcript, Exon1B, in pigs. In conclusion, using parthenotes, parental-origin-specific imprinting patterns in the porcine GNAS locus was comprehensively identified, and our approach paves the way for the discovery of novel imprinted genes and loci in a genomic context across species.

Integrative Analysis Revealing Human Adipose-Specific Genes and Consolidating Obesity Loci.

Identification of adipose-specific genes has contributed to an understanding of mechanisms underlying adipocyte development and obesity. Herein, our analyses of the recent Genotype-Tissue Expression (GTEx) database revealed 38 adipose-specific/enhanced protein coding genes, among which 3 genes were novel adipose-specific, and 414 highly differentially expressed genes (DEGs) between subcutaneous and omental adipose depots. By integrative analyses of genome-wide association studies (GWASs), 14 adipose-specific/enhanced genes and 60 DEGs were found to be associated with obesity-related traits and diseases, consolidating evidence for contribution of these genes to the regional fat distribution and obesity phenotypes. In addition, expression of HOXC cluster was up-regulated in subcutaneous adipose tissue, and the majority of the HOXB cluster was expressed highly in omental adipose tissue, indicating differential expression patterns of HOX clusters in adipose depots. Our findings on the distinct gene expression profiles in adipose tissue and their relation to obesity provide an important foundation for future functional biological studies and therapeutic targets in obesity and associated diseases.

Leptin attenuates the growth of rabbit mesenchymal stem cells via the extracellular signal-regulated kinase signaling pathway.

When stimulated, mesenchymal stem cells (MSCs) may differentiate into chondroblasts, adipocytes or osteoblasts. Leptin is an adipocyte-derived hormone, which regulates food intake and glucose homeostasis. The aim of the present study was to identify the potential role of mitogen-activated protein kinase in the leptin-induced growth of rabbit bone MSCs (rBMSCs). Various concentrations of leptin were used to culture rBMSCs and the viability of cells was observed as well as alterations in the phosphorylation state of extracellular signal-regulated kinase 1/2 (ERK1/2), c-Jun N-terminal kinase and p38. It was revealed that the growth of leptin-treated rBMSCs was primarily inhibited by phosphorylated ERK1/2, which was mediated by the leptin receptor. In conclusion, the results of the present study demonstrated that leptin inhibits the growth of rBMSCs principally via the ERK1/2 signaling pathway.

Whole-genome sequencing of 175 Mongolians uncovers population-specific genetic architecture and gene flow throughout North and East Asia.

The genetic variation in Northern Asian populations is currently undersampled. To address this, we generated a new genetic variation reference panel by whole-genome sequencing of 175 ethnic Mongolians, representing six tribes. The cataloged variation in the panel shows strong population stratification among these tribes, which correlates with the diverse demographic histories in the region. Incorporating our results with the 1000 Genomes Project panel identifies derived alleles shared between Finns and Mongolians/Siberians, suggesting that substantial gene flow between northern Eurasian populations has occurred in the past. Furthermore, we highlight that North, East, and Southeast Asian populations are more aligned with each other than these groups are with South Asian and Oceanian populations.

The genome of a Mongolian individual reveals the genetic imprints of Mongolians on modern human populations.

Mongolians have played a significant role in modern human evolution, especially after the rise of Genghis Khan (1162[?]-1227). Although the social cultural impacts of Genghis Khan and the Mongolian population have been well documented, explorations of their genome structure and genetic imprints on other human populations have been lacking. We here present the genome of a Mongolian male individual. The genome was de novo assembled using a total of 130.8-fold genomic data produced from massively parallel whole-genome sequencing. We identified high-confidence variation sets, including 3.7 million single nucleotide polymorphisms (SNPs) and 756,234 short insertions and deletions. Functional SNP analysis predicted that the individual has a pathogenic risk for carnitine deficiency. We located the patrilineal inheritance of the Mongolian genome to the lineage D3a through Y haplogroup analysis and inferred that the individual has a common patrilineal ancestor with Tibeto-Burman populations and is likely to be the progeny of the earliest settlers in East Asia. We finally investigated the genetic imprints of Mongolians on other human populations using different approaches. We found varying degrees of gene flows between Mongolians and populations living in Europe, South/Central Asia, and the Indian subcontinent. The analyses demonstrate that the genetic impacts of Mongolians likely resulted from the expansion of the Mongolian Empire in the 13th century. The genome will be of great help in further explorations of modern human evolution and genetic causes of diseases/traits specific to Mongolians.

Camelid genomes reveal evolution and adaptation to desert environments.

Bactrian camel (Camelus bactrianus), dromedary (Camelus dromedarius) and alpaca (Vicugna pacos) are economically important livestock. Although the Bactrian camel and dromedary are large, typically arid-desert-adapted mammals, alpacas are adapted to plateaus. Here we present high-quality genome sequences of these three species. Our analysis reveals the demographic history of these species since the Tortonian Stage of the Miocene and uncovers a striking correlation between large fluctuations in population size and geological time boundaries. Comparative genomic analysis reveals complex features related to desert adaptations, including fat and water metabolism, stress responses to heat, aridity, intense ultraviolet radiation and choking dust. Transcriptomic analysis of Bactrian camels further reveals unique osmoregulation, osmoprotection and compensatory mechanisms for water reservation underpinned by high blood glucose levels. We hypothesize that these physiological mechanisms represent kidney evolutionary adaptations to the desert environment. This study advances our understanding of camelid evolution and the adaptation of camels to arid-desert environments.