Using next-generation sequencing technology alone, we have successfully generated and assembled a draft sequence of the giant panda genome. The assembled contigs (2.25 gigabases (Gb)) cover approximately 94% of the whole genome, and the remaining gaps (0.05 Gb) seem to contain carnivore-specific repeats and tandem repeats. Comparisons with the dog and human showed that the panda genome has a lower divergence rate. The assessment of panda genes potentially underlying some of its unique traits indicated that its bamboo diet might be more dependent on its gut microbiome than its own genetic composition. We also identified more than 2.7 million heterozygous single nucleotide polymorphisms in the diploid genome. Our data and analyses provide a foundation for promoting mammalian genetic research, and demonstrate the feasibility for using next-generation sequencing technologies for accurate, cost-effective and rapid de novo assembly of large eukaryotic genomes.

The giant panda genome codes for all necessary enzymes associated with a carnivorous digestive system but lacks genes for enzymes needed to digest cellulose, the principal component of their bamboo diet. It has been posited that this iconic species must therefore possess microbial symbionts capable of metabolizing cellulose, but these symbionts have remained undetected. Here we examined 5,522 prokaryotic ribosomal RNA gene sequences in wild and captive giant panda fecal samples. We found lower species richness of the panda microbiome than of mammalian microbiomes for herbivores and nonherbivorous carnivores. We detected 13 operational taxonomic units closely related to Clostridium groups I and XIVa, both of which contain taxa known to digest cellulose. Seven of these 13 operational taxonomic units were unique to pandas compared with other mammals. Metagenomic analysis using ∼37-Mbp contig sequences from gut microbes recovered putative genes coding two cellulose-digesting enzymes and one hemicellulose-digesting enzyme, cellulase, β-glucosidase, and xylan 1,4-β-xylosidase, in Clostridium group I. Comparing glycoside hydrolase profiles of pandas with those of herbivores and omnivores, we found a moderate abundance of oligosaccharide-degrading enzymes for pandas (36%), close to that for humans (37%), and the lowest abundance of cellulases and endohemicellulases (2%), which may reflect low digestibility of cellulose and hemicellulose in the panda's unique bamboo diet. The presence of putative cellulose-digesting microbes, in combination with adaptations related to feeding, physiology, and morphology, show that giant pandas have evolved a number of traits to overcome the anatomical and physiological challenge of digesting a diet high in fibrous matter.

T cells become dysfunctional when they encounter self antigens or are exposed to chronic infection or to the tumour microenvironment1. The function of T cells is tightly regulated by a combinational co-stimulatory signal, and dominance of negative co-stimulation results in T cell dysfunction2. However, the molecular mechanisms that underlie this dysfunction remain unclear. Here, using an in vitro T cell tolerance induction system in mice, we characterize genome-wide epigenetic and gene expression features in tolerant T cells, and show that they are distinct from effector and regulatory T cells. Notably, the transcription factor NR4A1 is stably expressed at high levels in tolerant T cells. Overexpression of NR4A1 inhibits effector T cell differentiation, whereas deletion of NR4A1 overcomes T cell tolerance and exaggerates effector function, as well as enhancing immunity against tumour and chronic virus. Mechanistically, NR4A1 is preferentially recruited to binding sites of the transcription factor AP-1, where it represses effector-gene expression by inhibiting AP-1 function. NR4A1 binding also promotes acetylation of histone 3 at lysine 27 (H3K27ac), leading to activation of tolerance-related genes. This study thus identifies NR4A1 as a key general regulator in the induction of T cell dysfunction, and a potential target for tumour immunotherapy.

The panda lineage dates back to the late Miocene and ultimately leads to only one extant species, the giant panda (Ailuropoda melanoleuca). Although global climate change and anthropogenic disturbances are recognized to shape animal population demography their contribution to panda population dynamics remains largely unknown. We sequenced the whole genomes of 34 pandas at an average 4.7-fold coverage and used this data set together with the previously deep-sequenced panda genome to reconstruct a continuous demographic history of pandas from their origin to the present. We identify two population expansions, two bottlenecks and two divergences. Evidence indicated that, whereas global changes in climate were the primary drivers of population fluctuation for millions of years, human activities likely underlie recent population divergence and serious decline. We identified three distinct panda populations that show genetic adaptation to their environments. However, in all three populations, anthropogenic activities have negatively affected pandas for 3,000 years.

Phenotypic convergence between distantly related taxa often mirrors adaptation to similar selective pressures and may be driven by genetic convergence. The giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens) belong to different families in the order Carnivora, but both have evolved a specialized bamboo diet and adaptive pseudothumb, representing a classic model of convergent evolution. However, the genetic bases of these morphological and physiological convergences remain unknown. Through de novo sequencing the red panda genome and improving the giant panda genome assembly with added data, we identified genomic signatures of convergent evolution. Limb development genes DYNC2H1 and PCNT have undergone adaptive convergence and may be important candidate genes for pseudothumb development. As evolutionary responses to a bamboo diet, adaptive convergence has occurred in genes involved in the digestion and utilization of bamboo nutrients such as essential amino acids, fatty acids, and vitamins. Similarly, the umami taste receptor gene TAS1R1 has been pseudogenized in both pandas. These findings offer insights into genetic convergence mechanisms underlying phenotypic convergence and adaptation to a specialized bamboo diet.de novo genome | phenotype convergence | amino acid convergence | positive selection | pseudogenization

Systemic lupus erythematosus (SLE), a worldwide autoimmune disease with high heritability, shows differences in prevalence, severity and age of onset among different ancestral groups. Previous genetic studies have focused more on European populations, which appear to be the least affected. Consequently, the genetic variations that underlie the commonalities, differences and treatment options in SLE among ancestral groups have not been well elucidated. To address this, we undertake a genome-wide association study, increasing the sample size of Chinese populations to the level of existing European studies. Thirty-eight novel SLE-associated loci and incomplete sharing of genetic architecture are identified. In addition to the human leukocyte antigen (HLA) region, nine disease loci show clear ancestral differences and implicate antibody production as a potential mechanism for differences in disease manifestation. Polygenic risk scores perform significantly better when trained on ancestry-matched data sets. These analyses help to reveal the genetic basis for disparities in SLE among ancestral groups.

The carnivoran giant panda has a specialized bamboo diet, to which its alimentary tract is poorly adapted. Measurements of daily energy expenditure across five captive and three wild pandas averaged 5.2 megajoules (MJ)/day, only 37.7% of the predicted value (13.8 MJ/day). For the wild pandas, the mean was 6.2 MJ/day, or 45% of the mammalian expectation. Pandas achieve this exceptionally low expenditure in part by reduced sizes of several vital organs and low physical activity. In addition, circulating levels of thyroid hormones thyroxine (T4) and triiodothyronine (T3) averaged 46.9 and 64%, respectively, of the levels expected for a eutherian mammal of comparable size. A giant panda-unique mutation in the DUOX2 gene, critical for thyroid hormone synthesis, might explain these low thyroid hormone levels. A combination of morphological, behavioral, physiological, and genetic adaptations, leading to low energy expenditure, likely enables giant pandas to survive on a bamboo diet.

Most patients with cancer do not develop durable antitumor responses after programmed cell death protein 1 (PD-1) or programmed cell death ligand 1(PD-L1) checkpoint inhibition monotherapy because of an ephemeral reversal of T cell dysfunction and failure to promote long-lasting immunological T cell memory. Activating costimulatory pathways to induce stronger T cell activation may improve the efficacy of checkpoint inhibition and lead to durable antitumor responses. We performed single-cell RNA sequencing of more than 2000 tumor-infiltrating CD8 + T cells in mice receiving both PD-1 and GITR (glucocorticoid-induced tumor necrosis factor receptor-related protein) antibodies and found that this combination synergistically enhanced the effector function of expanded CD8 + T cells by restoring the balance of key homeostatic regulators CD226 and T cell immunoreceptor with Ig and ITIM domains (TIGIT), leading to a robust survival benefit. Combination therapy decreased CD8 + T cell dysfunction and induced a highly proliferative precursor effector memory T cell phenotype in a CD226-dependent manner.PD-1 inhibition rescued CD226 activity by preventing PD-1-Src homology region 2 (SHP2) dephosphophorylation of the CD226 intracellular domain, whereas GITR agonism decreased TIGIT expression. Unmasking the molecular pathways driving durable antitumor responses will be essential to the development of rational approaches to optimizing cancer immunotherapy.