The sequence and de novo assembly of the giant panda genome

Generation and characterization of giant panda induced pluripotent stem cells

The giant panda (Ailuropoda melanoleuca) stands as a flagship and umbrella species, symbolizing global biodiversity. While traditional assisted reproductive technology faces constraints in safeguarding the genetic diversity of giant pandas, induced pluripotent stem cells (iPSCs) known for their capacity to differentiate into diverse cells types, including germ cells, present a transformative potential for conservation of endangered animals. In this study, primary fibroblast cells were isolated from the giant panda, and giant panda iPSCs (GPiPSCs) were generated using a non-integrating episomal vector reprogramming method. Characterization of these GPiPSCs revealed their state of primed pluripotency and demonstrated their potential for differentiation. Furthermore, we innovatively formulated a species-specific chemically defined FACL medium and unraveled the intricate signaling pathway networks responsible for maintaining the pluripotency and fostering cell proliferation of GPiPSCs. This study provides key insights into rare species iPSCs, offering materials for panda characteristics research and laying the groundwork for in vitro giant panda gamete generation, potentially aiding endangered species conservation.

Late Pleistocene polar bear genomes reveal the timing of allele fixation in key genes associated with Arctic adaptation

The polar bear (Ursus maritimus) occupies a relatively narrow ecological niche, with many traits adapted for cold temperatures, movement across snow, ice and open water, and for consuming highly lipid-dense prey species. The divergence of polar bears from brown bears (Ursus arctos) and their adaptation to their Arctic lifestyle is a well-known example of rapid evolution. Previous research investigating whole genomes uncovered twelve key genes that are highly differentiated between polar and brown bears, show signatures of selection in the polar bear lineage, and are associated with polar bear adaptation to the Arctic environment. Further research suggested fixed derived alleles in these genes arose from selection on both standing variation and de novo mutations in the evolution of polar bears. Here, we reevaluate these findings based on a larger and geographically more representative dataset of 119 polar bears and 135 brown bears, and assess the timing of derived allele fixation in polar bears by incorporating the genomes of two Late Pleistocene individuals (aged 130-100,000 years old and 100-70,000 years old). In contrast with previous results, we found no evidence of derived alleles fixed in present-day polar bears within the key genes arising from de novo mutation. Most derived alleles fixed in present-day polar bears were also fixed in the Late Pleistocene polar bears, suggesting selection occurred prior to 70,000 years ago. However, some derived alleles fixed in present-day polar bears were not fixed in the two Late Pleistocene polar bears, including at sites within APOB, LYST, and TTN. These three genes are associated with cardiovascular function, metabolism, and pigmentation, suggesting selection may have acted on different loci at different times.

The Chinese longsnout catfish genome provides novel insights into the feeding preference and corresponding metabolic strategy of carnivores

Fish show variation in feeding habits to adapt to complex environments. However, the genetic basis of feeding preference and the corresponding metabolic strategies that differentiate feeding habits remain elusive. Here, by comparing the whole genome of a typical carnivorous fish (Leiocassis longirostris Günther) with that of herbivorous fish, we identify 250 genes through both positive selection and rapid evolution, including taste receptor taste receptor type 1 member 3 (tas1r3) and trypsin We demonstrate that tas1r3 is required for carnivore preference in tas1r3-deficient zebrafish and in a diet-shifted grass carp model. We confirm that trypsin correlates with the metabolic strategies of fish with distinct feeding habits. Furthermore, marked alterations in trypsin activity and metabolic profiles are accompanied by a transition of feeding preference in tas1r3-deficient zebrafish and diet-shifted grass carp. Our results reveal a conserved adaptation between feeding preference and corresponding metabolic strategies in fish, and provide novel insights into the adaptation of feeding habits over the evolution course.

Research progress on formation mechanism of pearl

Pearls are deeply cherished for their rich color and gorgeous luster, and their quality directly affects their value. Currently, the evaluation of pearl quality is mainly based on four aspects: color, shape, size and smoothness. The quality of pearls is influenced by a variety of factors, categorized into internal factors, such as the structural composition of the nacreous layer and genetic factors of the mussels, and external factors, including the aquaculture environment. Existing research results indicates that genetic factors are the dominant factor controlling the pearl quality. However, the macromolecules such as metal ions, organic pigments and various physical and chemical factors in the aquaculture water environment will also significantly impact pearl quality. Among these, matrix proteins are organic macromolecules found in the nacreous layer that play an important role in pearl quality. They participate in the deposition of calcium carbonate and the construction of the organic framework, affecting the pearls' size and shape. The color of pearls is influenced by the deposition of metal ions, the transport of organic pigments and the regulation of microstructure.

Comparative transcriptome and methylome of polar bears, giant and red pandas reveal diet-driven adaptive evolution

Epigenetic regulation plays an important role in the evolution of species adaptations, yet little information is available on the epigenetic mechanisms underlying the adaptive evolution of bamboo-eating in both giant pandas (Ailuropoda melanoleuca) and red pandas (Ailurus fulgens). To investigate the potential contribution of epigenetic to the adaptive evolution of bamboo-eating in giant and red pandas, we performed hepatic comparative transcriptome and methylome analyses between bamboo-eating pandas and carnivorous polar bears (Ursus maritimus). We found that genes involved in carbohydrate, lipid, amino acid, and protein metabolism showed significant differences in methylation and expression levels between the two panda species and polar bears. Clustering analysis of gene expression revealed that giant pandas did not form a sister group with the more closely related polar bears, suggesting that the expression pattern of genes in livers of giant pandas and red pandas have evolved convergently driven by their similar diets. Compared to polar bears, some key genes involved in carbohydrate metabolism and biological oxidation and cholesterol synthesis showed hypomethylation and higher expression in giant and red pandas, while genes involved in fat digestion and absorption, fatty acid metabolism, lysine degradation, resistance to lipid peroxidation and detoxification showed hypermethylation and low expression. Our study elucidates the special nutrient utilization mechanism of giant pandas and red pandas and provides some insights into the molecular mechanism of their adaptive evolution of bamboo feeding. This has important implications for the breeding and conservation of giant pandas and red pandas.

Genome-Wide Identification of Specific Genetic Loci Common to Sheep and Goat

Sheep and goat may become carriers of some zoonotic diseases. They are important livestock and experimental model animals for human beings. The fast and accurate identification of genetic materials originating from sheep and goat can prevent and inhibit the spread of some zoonotic diseases, monitor market product quality, and maintain the stability of animal husbandry and food industries. This study proposed a methodology for identifying sheep and goat common specific sites from a genome-wide perspective. A total of 150 specific sites were selected from three data sources, including the coding sequences of single copy genes from nine species (sheep, goat, cow, pig, dog, horse, human, mouse, and chicken), the dbSNPs for these species, and human 100-way alignment data. These 150 sites exhibited low intraspecific heterogeneity in the resequencing data of 1450 samples from five species (sheep, goat, cow, pig, and chicken) and high interspecific divergence in the human 100-way alignment data after quality control. The results were proven to be reliable at the data level. Using the process proposed in this study, specific sites of other species can be screened, and genome-level species identification can be performed using the screened sites.

A new chromosome-scale genome of wild Brassica oleracea provides insights into the domestication of Brassica crops

The cultivated diploid Brassica oleracea is an important vegetable crop, but the genetic basis of its domestication remains largely unclear in the absence of high-quality reference genomes of wild B. oleracea. Here, we report the first chromosome-level assembly of the wild Brassica oleracea L. W03 genome (total genome size, 630.7 Mb; scaffold N50, 64.6 Mb). Using the newly assembled W03 genome, we constructed a gene-based B. oleracea pangenome and identified 29 744 core genes, 23 306 dispensable genes, and 1896 private genes. We re-sequenced 53 accessions, representing six potential wild B. oleracea progenitor species. The results of the population genomic analysis showed that the wild B. oleracea populations had the highest level of diversity and represents the most closely related population to modern-day horticultural B. oleracea. In addition, the WUSCHEL gene was found to play a decisive role in domestication and to be involved in cauliflower and broccoli curd formation. We also illustrate the loss of disease-resistance genes during selection for domestication. Our results provide new insights into the domestication of B. oleracea and will facilitate the future genetic improvement of Brassica crops.

A chromosome-scale genome of Rhus chinensis Mill. provides new insights into plant-insect interaction and gallotannins biosynthesis

Rhus chinensis Mill., an economically valuable Anacardiaceae species, is parasitized by the galling aphid Schlechtendalia chinensis, resulting in the formation of the Chinese gallnut (CG). Here, we report a chromosomal-level genome assembly of R. chinensis, with a total size of 389.40 Mb and scaffold N50 of 23.02 Mb. Comparative genomic and transcriptome analysis revealed that the enhanced structure of CG and nutritional metabolism contribute to improving the adaptability of R. chinensis to S. chinensis by supporting CG and galling aphid growth. CG was observed to be abundant in hydrolysable tannins (HT), particularly gallotannin and its isomers. Tandem repeat clusters of dehydroquinate dehydratase/shikimate dehydrogenase (DQD/SDH) and serine carboxypeptidase-like (SCPL) and their homologs involved in HT production were determined as specific to HT-rich species. The functional differentiation of DQD/SDH tandem duplicate genes and the significant contraction in the phenylalanine ammonia-lyase (PAL) gene family contributed to the accumulation of gallic acid and HT while minimizing the production of shikimic acid, flavonoids, and condensed tannins in CG. Furthermore, we identified one UDP glucosyltransferase (UGT84A), three carboxylesterase (CXE), and six SCPL genes from conserved tandem repeat clusters that are involved in gallotannin biosynthesis and hydrolysis in CG. We then constructed a regulatory network of these genes based on co-expression and transcription factor motif analysis. Our findings provide a genomic resource for the exploration of the underlying mechanisms of plant-galling insect interaction and highlight the importance of the functional divergence of tandem duplicate genes in the accumulation of secondary metabolites.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

Genomics comparisons of three chromosome-level mudskipper genome assemblies reveal molecular clues for water-to-land evolution and adaptation

INTRODUCTION

Mudskippers are a large group of amphibious fishes that have developed many morphological and physiological capacities to live on land. Genomics comparisons of chromosome-level genome assemblies of three representative mudskippers, Boleophthalmus pectinirostris (BP), Periophthalmus magnuspinnatus (PM) and P. modestus (PMO), may be able to provide novel insights into the water-to-land evolution and adaptation.

METHODS

Two chromosome-level genome assemblies for BP and PM were respectively sequenced by an integration of PacBio, Nanopore and Hi-C sequencing. A series of standard assembly and annotation pipelines were subsequently performed for both mudskippers. We also re-annotated the PMO genome, downloaded from NCBI, to obtain a redundancy-reduced annotation. Three-way comparative analyses of the three mudskipper genomes in a large scale were carried out to discover detailed genomic differences, such as different gene sizes, and potential chromosomal fission and fusion events. Comparisons of several representative gene families among the three amphibious mudskippers and some other teleosts were also performed to find some molecular clues for terrestrial adaptation.

RESULTS

We obtained two high-quality haplotype genome assemblies with 23 and 25 chromosomes for BP and PM respectively. We also found two specific chromosome fission events in PM. Ancestor chromosome analysis has discovered a common fusion event in mudskipper ancestor. This fusion was then retained in all the three mudskipper species. A loss of some SCPP (secretory calcium-binding phosphoprotein) genes were identified in the three mudskipper genomes, which could lead to reduction of scales for a part-time terrestrial residence. The loss of aanat1a gene, encoding an important enzyme (arylalkylamine N-acetyltransferase 1a, AANAT1a) for dopamine metabolism and melatonin biosynthesis, was confirmed in PM but not in PMO (as previously reported existence in BP), suggesting a better air vision of PM than both PMO and BP. Such a tiny variation within the genus Periophthalmus exemplifies to prove a step-by-step evolution for the mudskippers' water-to-land adaptation.

CONCLUSION

These high-quality mudskipper genome assemblies will become valuable genetic resources for in-depth discovery of genomic evolution for the terrestrial adaptation of amphibious fishes.

Natural mating ability is associated with gut microbiota composition and function in captive male giant pandas

The issue of poor sexual performance of some male giant pandas seriously impairs the growth and the genetic diversity of the captive population, yet there is still no clear understanding of the cause of the loss of this ability and its underlying mechanism. In this study, we analyzed the gut microbiota and its function in 72 fecal samples obtained from 20 captive male giant pandas, with an equal allocation between individuals capable and incapable of natural mating. Additionally, we investigated fecal hormone levels and behavioral differences between the two groups. A correlation analysis was then conducted among these factors to explore the influencing factors of their natural mating ability. The results showed significant differences in the composition of gut microbiota between the two groups of male pandas. The capable group had significantly higher abundance of Clostridium sensu stricto 1 (p adjusted = .0021, GLMM), which was positively correlated with fatty acid degradation and two-component system functions (Spearman, p adjusted < .05). Additionally, the capable group showed higher gene abundance in gut microbiota function including purine and pyrimidine metabolism and galactose metabolism, as well as pathways related to biological processes such as ribosome and homologous recombination (DEseq2, p adjusted < .05). We found no significant differences in fecal cortisol and testosterone levels between the two groups, and no difference was found in their behavior either. Our study provides a theoretical and practical basis for further studying the behavioral degradation mechanisms of giant pandas and other endangered mammal species.

The Potential of Co-Evolution and Interactions of Gut Bacteria-Phages in Bamboo-Eating Pandas: Insights from Dietary Preference-Based Metagenomic Analysis

Bacteria and phages are two of the most abundant biological entities in the gut microbiome, and diet and host phylogeny are two of the most critical factors influencing the gut microbiome. A stable gut bacterial community plays a pivotal role in the host's physiological development and immune health. A phage is a virus that directly infects bacteria, and phages' close associations and interactions with bacteria are essential for maintaining the stability of the gut bacterial community and the entire microbial ecosystem. Here, we utilized 99 published metagenomic datasets from 38 mammalian species to investigate the relationship (diversity and composition) and potential interactions between gut bacterial and phage communities and the impact of diet and phylogeny on these communities. Our results highlight the co-evolutionary potential of bacterial-phage interactions within the mammalian gut. We observed a higher alpha diversity in gut bacteria than in phages and identified positive correlations between bacterial and phage compositions. Furthermore, our study revealed the significant influence of diet and phylogeny on mammalian gut bacterial and phage communities. We discovered that the impact of dietary factors on these communities was more pronounced than that of phylogenetic factors at the order level. In contrast, phylogenetic characteristics had a more substantial influence at the family level. The similar omnivorous dietary preference and closer phylogenetic relationship (family Ursidae) may contribute to the similarity of gut bacterial and phage communities between captive giant panda populations (GPCD and GPYA) and omnivorous animals (OC; including Sun bear, brown bear, and Asian black bear). This study employed co-occurrence microbial network analysis to reveal the potential interaction patterns between bacteria and phages. Compared to other mammalian groups (carnivores, herbivores, and omnivores), the gut bacterial and phage communities of bamboo-eating species (giant pandas and red pandas) exhibited a higher level of interaction. Additionally, keystone species and modular analysis showed the potential role of phages in driving and maintaining the interaction patterns between bacteria and phages in captive giant pandas. In sum, gaining a comprehensive understanding of the interaction between the gut microbiota and phages in mammals is of great significance, which is of great value in promoting healthy and sustainable mammals and may provide valuable insights into the conservation of wildlife populations, especially endangered animal species.

Evolution of Endogenous Retroviruses in the Subfamily of Caprinae

The interest in endogenous retroviruses (ERVs) has been fueled by their impact on the evolution of the host genome. In this study, we used multiple pipelines to conduct a de novo exploration and annotation of ERVs in 13 species of the Caprinae subfamily. Through analyses of sequence identity, structural organization, and phylogeny, we defined 28 ERV groups within Caprinae, including 19 gamma retrovirus groups and 9 beta retrovirus groups. Notably, we identified four recent and potentially active groups prevalent in the Caprinae genomes. Additionally, our investigation revealed that most long noncoding genes (lncRNA) and protein-coding genes (PC) contain ERV-derived sequences. Specifically, we observed that ERV-derived sequences were present in approximately 75% of protein-coding genes and 81% of lncRNA genes in sheep. Similarly, in goats, ERV-derived sequences were found in approximately 74% of protein-coding genes and 75% of lncRNA genes. Our findings lead to the conclusion that the majority of ERVs in the Caprinae genomes can be categorized as fossils, representing remnants of past retroviral infections that have become permanently integrated into the genomes. Nevertheless, the identification of the Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25 groups indicates the presence of relatively recent and potentially active ERVs in these genomes. These particular groups may contribute to the ongoing evolution of the Caprinae genome. The identification of putatively active ERVs in the Caprinae genomes raises the possibility of harnessing them for future genetic marker development.

Conservation Genomics and Metagenomics of Giant and Red Pandas in the Wild

Giant pandas and red pandas are endangered species with similar specialized bamboo diet and partial sympatric distribution in China. Over the last two decades, the rapid development of genomics and metagenomics research on these species has enriched our knowledge of their biology, ecology, physiology, genetics, and evolution, which is crucial and useful for their conservation. We describe the evolutionary history, endangerment processes, genetic diversity, and population structure of wild giant pandas and two species of red pandas (Chinese and Himalayan red pandas). In addition, we explore how genomics and metagenomics studies have provided insight into the convergent adaptation of pandas to the specialized bamboo diet. Finally, we discuss how these findings are applied to effective conservation management of giant and red pandas in the wild and in captivity to promote the long-term persistence of these species.

Mozambican Coffea accessions from Ibo and Quirimba Islands: identification and geographical distribution

Mozambique does not have a tradition of farming Coffea arabica or Coffea canephora, the two species that dominate the worldwide coffee market. However, native coffee plants have been growing spontaneously and in some cases cultivated in the Ibo and Quirimba islands in the north of the country and Inhambane province in the south. Historically there has been confusion over the precise taxonomic classification of these indigenous coffee plants, with different botanists identifying the species as C. racemosa, C. zanguebariae or various synonyms of both. The present research aims to clarify the subject and provide new information on these little-described coffee species which may prove valuable as new breeding material for future cultivars, something that is sorely needed to face the present and future challenges of coffee production. Leaf samples were collected from 40 accessions from Ibo Island, Quirimba Island and Inhambane province. The samples were sequenced by whole-genome technology and WGS reads were filtered to identify relevant SNP variants. Diversity among the samples was assessed by PCA, and a phylogenetic tree including several Coffea species was built using additional data available in public databases. Experimental data confirm the presence of C. zanguebariae as the only coffee species present in both Ibo and Quirimba Islands, while it appears that C. racemosa is exclusive to the southern Inhambane province. The present research provides the most detailed analysis so far on the genetic identity of the traditional Mozambican coffee crops. This is the prerequisite for undertaking further scientific studies on these almost unknown coffee species and for starting agronomic development programs for the economic revival of Ibo and Quirimba islands based on coffee cultivation. Furthermore, these species could provide much-needed genetic material for the breeding of new hybrids with the two main commercial coffee species.

Multiple independent losses of the biosynthetic pathway for two tropane alkaloids in the Solanaceae family

Hyoscyamine and scopolamine (HS), two valuable tropane alkaloids of significant medicinal importance, are found in multiple distantly related lineages within the Solanaceae family. Here we sequence the genomes of three representative species that produce HS from these lineages, and one species that does not produce HS. Our analysis reveals a shared biosynthetic pathway responsible for HS production in the three HS-producing species. We observe a high level of gene collinearity related to HS synthesis across the family in both types of species. By introducing gain-of-function and loss-of-function mutations at key sites, we confirm the reduced/lost or re-activated functions of critical genes involved in HS synthesis in both types of species, respectively. These findings indicate independent and repeated losses of the HS biosynthesis pathway since its origin in the ancestral lineage. Our results hold promise for potential future applications in the artificial engineering of HS biosynthesis in Solanaceae crops.

Impact of Holocene environmental change on the evolutionary ecology of an Arctic top predator

The Arctic is among the most climatically sensitive environments on Earth, and the disappearance of multiyear sea ice in the Arctic Ocean is predicted within decades. As apex predators, polar bears are sentinel species for addressing the impact of environmental variability on Arctic marine ecosystems. By integrating genomics, isotopic analysis, morphometrics, and ecological modeling, we investigate how Holocene environmental changes affected polar bears around Greenland. We uncover reductions in effective population size coinciding with increases in annual mean sea surface temperature, reduction in sea ice cover, declines in suitable habitat, and shifts in suitable habitat northward. Furthermore, we show that west and east Greenlandic polar bears are morphologically, and ecologically distinct, putatively driven by regional biotic and genetic differences. Together, we provide insights into the vulnerability of polar bears to environmental change and how the Arctic marine ecosystem plays a vital role in shaping the evolutionary and ecological trajectories of its inhabitants.

Genome wide copy number variations using Porcine 60K SNP Beadchip in Landlly pigs

In the present study, Porcine 60K SNP genotype data from 69 Landlly pigs were used to explore Copy Number Variations (CNVs) across the autosomes. A total of 386 CNVs were identified using Hidden Markov Model (HMM) in PennCNV software, which were subsequently aggregated to 115 CNV regions (CNVRs). Among the total detected CNVRs, 58 gain, 49 were loss type while remaining 8 events were both gain and loss types. Identified CNVRs covered 12.5 Mb (0.55%) of Sus scrofa reference 11.1 genome. Comparison of our results with previous investigations on pigs revealed that approximately 75% CNVRs were novel, which may be due to differences in genetic background, environment and implementation of artificial selection in Landlly pigs. Functional annotation and pathway analysis showed the significant enrichment of 267 well-annotated Sus scrofa genes in CNVRs. These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity. Additionally, KIT and FUT1 were two major genes detected on CNVR in our population. This investigation provided a comprehensive overview of CNV distribution in the Indian porcine genome for the first time, which may be useful for further investigating the association of important quantitative traits in Landlly pigs.Highlights115 CNVRs were identified in 69 Landlly pig population.Approximately 75% detected CNVRs were novel for Landlly population.Significant enrichment of 267 well-annotated Sus scrofa genes observed in these CNVRs.These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity.Comprehensive CNV map in the Indian porcine genome developed for the first time.

Integrated transcriptomic and proteomic analyses reveal the mechanism of easy acceptance of artificial pelleted diets during food habit domestication in Largemouth bass (Micropterus salmoides)

Acceptance of artificial pelleted diets contributes to increasing the cultured areas and output of carnivorous fish. However, the mechanism of acceptance of artificial pelleted diets remains largely unknown. In this study, the easy acceptance of artificial pelleted diets (EAD) group and the not easy acceptance of artificial pelleted diets (NAD) group of Largemouth bass (Micropterus salmoides) were divided based on the ratios of stomach weight/body weight (SB) after 0.5 h feeding, which was bigger than 18% in the EAD group and ranged from 8 to 12% in the NAD group. Through transcriptome and proteome sequencing, a total of 2463 differentially expressed genes (DEGs) and 230 differentially expressed proteins (DEPs) were identified, respectively. Integrated analyses of transcriptome and proteome data revealed that 152 DEPs were matched with the corresponding DEGs (named co-DEGs-DEPs), and 54 co-DEGs-DEPs were enriched in 16 KEGG pathways, including the metabolic pathways, steroid biosynthesis, fatty acid biosynthesis, etc. Furthermore, 3 terpenoid backbone biosynthesis-related genes (Hmgcr, Hmgcs, and Fdps) in metabolic pathways, 10 steroid biosynthesis-related genes (Fdft1, Sqle, Lss, Cyp51a1, Tm7sf2, Nsdhl, Hsd17b7, Dhcr24, Sc5d, and Dhcr7), and 3 fatty acid biosynthesis-related genes (Acaca, Fasn, and Ascl) were all up-regulated in the EAD group, suggesting that the lipid metabolism pathway and steroid biosynthesis pathway play important roles in early food habit domestication in Largemouth bass. In addition, the detection results of randomly selected 15 DEGs and 15 DEPs indicated that both transcriptome and proteome results in the study were reliable. Our study provides useful information for further research on the mechanisms of food habit domestication in fish.

Unravelling the genomic features, phylogeny and genetic basis of tooth ontogenesis in Characiformes through analysis of four genomes

Characiformes is a diverse and evolutionarily significant order of freshwater fish encompassing over 2,300 species. Despite its diversity, our understanding of Characiformes' evolutionary relationships and adaptive mechanisms is limited due to insufficient genome sequences. In this study, we sequenced and assembled the genomes of four Characiformes species, three of which were chromosome-level assemblies. Our analyses revealed dynamic changes in gene family evolution, repeat sequences and variations in chromosomal collinearity within these genomes. With the assembled genomes, we were not only able to elucidate the evolutionary relationship of the four main orders in Otophysi but also indicated Characiformes as the paraphyletic group. Comparative genomic analysis with other available fish genomes shed light on the evolution of genes related to tooth development in Characiformes. Notably, variations in the copy number of secretory calcium-binding phosphoproteins (SCPP) genes were observed among different orders of Otophysi, indicating their potential contribution to the diversity of tooth types. Our study offers invaluable genome sequences and novel insights into Characiformes' evolution, paving the way for further genomic and evolutionary research in fish.

De novo assembly and annotation of the singing mouse genome

BACKGROUND

Developing genomic resources for a diverse range of species is an important step towards understanding the mechanisms underlying complex traits. Specifically, organisms that exhibit unique and accessible phenotypes-of-interest allow researchers to address questions that may be ill-suited to traditional model organisms. We sequenced the genome and transcriptome of Alston's singing mouse (Scotinomys teguina), an emerging model for social cognition and vocal communication. In addition to producing advertisement songs used for mate attraction and male-male competition, these rodents are diurnal, live at high-altitudes, and are obligate insectivores, providing opportunities to explore diverse physiological, ecological, and evolutionary questions.

RESULTS

Using PromethION, Illumina, and PacBio sequencing, we produced an annotated genome and transcriptome, which were validated using gene expression and functional enrichment analyses. To assess the usefulness of our assemblies, we performed single nuclei sequencing on cells of the orofacial motor cortex, a brain region implicated in song coordination, identifying 12 cell types.

CONCLUSIONS

These resources will provide the opportunity to identify the molecular basis of complex traits in singing mice as well as to contribute data that can be used for large-scale comparative analyses.

The unique gut microbiome of giant pandas involved in protein metabolism contributes to the host's dietary adaption to bamboo

BACKGROUND

The gut microbiota of the giant panda (Ailuropoda melanoleuca), a global symbol of conservation, are believed to be involved in the host's dietary switch to a fibrous bamboo diet. However, their exact roles are still largely unknown.

RESULTS

In this study, we first comprehensively analyzed a large number of gut metagenomes giant pandas (n = 322), including 98 pandas sequenced in this study with deep sequencing (Illumina) and third-generation sequencing (nanopore). We reconstructed 408 metagenome-assembled genomes (MAGs), and 148 of which (36.27%) were near complete. The most abundant MAG was classified as Streptococcus alactolyticus. A pairwise comparison of the metagenomes and meta-transcriptomes in 14 feces revealed genes involved in carbohydrate metabolism were lower, but those involved in protein metabolism were greater in abundance and expression in giant pandas compared to those in herbivores and omnivores. Of note, S. alactolyticus was positively correlated to the KEGG modules of essential amino-acid biosynthesis. After being isolated from pandas and gavaged to mice, S. alactolyticus significantly increased the relative abundance of essential amino acids in mice jejunum.

CONCLUSIONS

The study highlights the unique protein metabolic profiles in the giant panda's gut microbiome. The findings suggest that S. alactolyticus is an important player in the gut microbiota that contributes to the giant panda's dietary adaptation by more involvement in protein rather than carbohydrate metabolism. Video Abstract.

Comparative study of the digestion and metabolism related genes' expression changes during the postnatal food change in different dietary mammals

The changes in the expression of genes related to digestion and metabolism may be various in different dietary mammals from juvenile to adult, especially, the giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens), which were once carnivores but have shifted to being specialized bamboo eaters, are unique features of their changes are more unclear. To elucidate the changing patterns of gene expression related to digestion and metabolism from juvenile to adult in different dietary mammals, we performed transcriptome analysis of the liver or pancreas in giant and red pandas, herbivorous rabbits (Oryctolagus cuniculus) and macaques (Macaca mulatta), carnivorous ferrets (Mustela putorius furo), and omnivorous mice (Mus musculus) from juvenile to adult. During the transition from juvenile to adulthood, giant and red pandas, as well as rabbits and macaques, show significant upregulation of key genes for carbohydrate metabolism, such as starch hydrolysis and sucrose metabolism, and unsaturated fatty acid metabolism, such as linoleic acid, while there is no significant difference in the expression of key genes for fatty acid β-oxidation. A large number of amino acid metabolism related genes were upregulated in adult rabbits and macaques compared to juveniles. While adult giant and red pandas mainly showed upregulation of key genes for arginine synthesis and downregulation of key genes for arginine and lysine degradation. In adult stages, mouse had significantly higher expression patterns in key genes for starch hydrolysis and sucrose metabolism, as well as lipid and protein metabolism. In contrast to general expectations, genes related to lipid, amino acid and protein metabolism were significantly higher expressed in adult group of ferrets, which may be related to their high metabolic levels. Our study elucidates the pattern of changes in the expression of genes related to digestion and metabolism from juvenile to adult in different dietary mammals, with giant and red pandas showing adaptations associated with specific nutritional limitations of bamboo.

Age-related alterations in metabolome and microbiome provide insights in dietary transition in giant pandas

We conducted UPLC-MS-based metabolomics, 16S rRNA, and metagenome sequencing on the fecal samples of 44 captive giant pandas (Ailuropoda melanoleuca) from four age groups (i.e., Cub, Young, Adult, and Old) to comprehensively understand age-related changes in the metabolism and gut microbiota of giant pandas. We characterized the metabolite profiles of giant pandas based on 1,376 identified metabolites, with 152 significantly differential metabolites (SDMs) found across the age groups. We found that the metabolites and the composition/function of the gut microbiota changed in response to the transition from a milk-dominant diet in panda cubs to a bamboo-specific diet in young and adult pandas. Lipid metabolites such as choline and hippuric acid were enriched in the Cub group, and many plant secondary metabolites were significantly higher in the Young and Adult groups, while oxidative stress and inflammatory related metabolites were only found in the Old group. However, there was a decrease in the α-diversity of gut microbiota in adult and old pandas, who exclusively consume bamboo. The abundance of bacteria related to the digestion of cellulose-rich food, such as Firmicutes, Streptococcus, and Clostridium, significantly increased from the Cub to the Adult group, while the abundance of beneficial bacteria such as Faecalibacterium, Sarcina, and Blautia significantly decreased. Notably, several potential pathogenic bacteria had relatively high abundances, especially in the Young group. Metagenomic analysis identified 277 CAZyme genes including cellulose degrading genes, and seven of the CAZymes had abundances that significantly differed between age groups. We also identified 237 antibiotic resistance genes (ARGs) whose number and diversity increased with age. We also found a significant positive correlation between the abundance of bile acids and gut bacteria, especially Lactobacillus and Bifidobacterium. Our results from metabolome, 16S rRNA, and metagenome data highlight the important role of the gut microbiota-bile acid axis in the regulation of age-related metabolism and provide new insights into the lipid metabolism of giant pandas. IMPORTANCE The giant panda is a member of the order Carnivora but is entirely herbivorous. The giant panda's specialized diet and related metabolic mechanisms have not been fully understood. It is therefore crucial to investigate the dynamic changes in metabolites as giant pandas grow and physiologically adapt to their herbivorous diet. This study conducted UPLC-MS-based metabolomics 16S rRNA, and metagenome sequencing on the fecal samples of captive giant pandas from four age groups. We found that metabolites and the composition/function of gut microbiota changed in response to the transition from a milk-dominant diet in cubs to a bamboo-specific diet in young and adult pandas. The metabolome, 16S rRNA, and metagenome results highlight that the gut microbiota-bile acid axis has an important role in the regulation of age-related metabolism, and our study provides new insights into the lipid metabolism of giant pandas.

Tuning the optoelectronic properties of acridine-triphenylamine (ACR-TPA) based novel hole transporting material for high efficiency perovskite and organic solar cell

In this research, five distinct small donor molecules (designated as ACR-TPA-X1, ACR-TPA-X2, ACR-TPA-X3, ACR-TPA-X4, ACR-TPA-X5) are constructed by replacing the methoxy groups on both sides of the model molecule (ACR-TPA-R) with thiophene bridged acceptor moieties. We have used the B3LYP/6-31G (d,p) model for our computational studies. Our model molecule's morphological alteration has resulted in a lowered E_g of 1.77-2.51 eV as compared to model (ACR-TPA-R=3.84 eV). ACR-TPA-X2 investigated the λ_max at 776 nm. ACR-TPA-X4 was found to be most miscible with dichloromethane (DCM). The greatest V_OC(1.21 eV) was observed in ACR-TPA-X1. Among all of the variants, ACR-TPA-X1 had the highest PCE (23.42%). It was found that ACR-TPA-X4 had the highest electron mobility (0.00370 eV) and ACR-TPA-X5 had the highest hole mobility (0.00324 eV) of all the materials examined. The findings prove the worth of the methods used, paving the way for the development of effective small donors for OSCs and HTMs for PSCs.

Long-term cold, freezing and drought: overlapping and specific regulatory mechanisms and signal transduction in tea plant (Camellia sinensis (L.) Kuntze)

Introduction

Low temperatures and drought are two main environmental constraints reducing the yield and geographical distribution of horticultural crops worldwide. Understanding the genetic crosstalk between stress responses has potential importance for crop improvement.

Methods

In this study, Illumina RNA-seq and Pac-Bio genome resequencing were used to annotate genes and analyze transcriptome dynamics in tea plants under long-term cold, freezing, and drought.

Results

The highest number of differentially expressed genes (DEGs) was identified under long-term cold (7,896) and freezing (7,915), with 3,532 and 3,780 upregulated genes, respectively. The lowest number of DEGs was observed under 3-day drought (47) and 9-day drought (220), with five and 112 genes upregulated, respectively. The recovery after the cold had 6.5 times greater DEG numbers as compared to the drought recovery. Only 17.9% of cold-induced genes were upregulated by drought. In total, 1,492 transcription factor genes related to 57 families were identified. However, only 20 transcription factor genes were commonly upregulated by cold, freezing, and drought. Among the 232 common upregulated DEGs, most were related to signal transduction, cell wall remodeling, and lipid metabolism. Co-expression analysis and network reconstruction showed 19 genes with the highest co-expression connectivity: seven genes are related to cell wall remodeling (GATL7, UXS4, PRP-F1, 4CL, UEL-1, UDP-Arap, and TBL32), four genes are related to calcium-signaling (PXL1, Strap, CRT, and CIPK6), three genes are related to photo-perception (GIL1, CHUP1, and DnaJ11), two genes are related to hormone signaling (TTL3 and GID1C-like), two genes are involved in ROS signaling (ERO1 and CXE11), and one gene is related to the phenylpropanoid pathway (GALT6).

Discussion

Based on our results, several important overlapping mechanisms of long-term stress responses include cell wall remodeling through lignin biosynthesis, o-acetylation of polysaccharides, pectin biosynthesis and branching, and xyloglucan and arabinogalactan biosynthesis. This study provides new insight into long-term stress responses in woody crops, and a set of new target candidate genes were identified for molecular breeding aimed at tolerance to abiotic stresses.

Population genomic analysis provides evidence of the past success and future potential of South China tiger captive conservation

BACKGROUND

Among six extant tiger subspecies, the South China tiger (Panthera tigris amoyensis) once was widely distributed but is now the rarest one and extinct in the wild. All living South China tigers are descendants of only two male and four female wild-caught tigers and they survive solely in zoos after 60 years of effective conservation efforts. Inbreeding depression and hybridization with other tiger subspecies were believed to have occurred within the small, captive South China tiger population. It is therefore urgently needed to examine the genomic landscape of existing genetic variation among the South China tigers.

RESULTS

In this study, we assembled a high-quality chromosome-level genome using long-read sequences and re-sequenced 29 high-depth genomes of the South China tigers. By combining and comparing our data with the other 40 genomes of six tiger subspecies, we identified two significantly differentiated genomic lineages among the South China tigers, which harbored some rare genetic variants introgressed from other tiger subspecies and thus maintained a moderate genetic diversity. We noticed that the South China tiger had higher F_ROH values for longer runs of homozygosity (ROH > 1 Mb), an indication of recent inbreeding/founder events. We also observed that the South China tiger had the least frequent homozygous genotypes of both high- and moderate-impact deleterious mutations, and lower mutation loads than both Amur and Sumatran tigers. Altogether, our analyses indicated an effective genetic purging of deleterious mutations in homozygous states from the South China tiger, following its population contraction with a controlled increase in inbreeding based on its pedigree records.

CONCLUSIONS

The identification of two unique founder/genomic lineages coupled with active genetic purging of deleterious mutations in homozygous states and the genomic resources generated in our study pave the way for a genomics-informed conservation, following the real-time monitoring and rational exchange of reproductive South China tigers among zoos.

The chromosome-level genome of Cherax quadricarinatus

Red claw crayfish (Cherax quadricarinatus) is an aquatic crustacean with considerable potential for the commercial culture and an ideal model for studying the mechanism of sex determination. To provide better genomic resources, we assembled a chromosome-level genome with a size of 5.26 Gb and contig N50 of 144.33 kb. Nearly 90% of sequences were anchored to 100 chromosomes, which represents the high-quality crustacean genome with the largest number of chromosomes ever reported. The genome contained 78.69% repeat sequences and 20,460 protein-coding genes, of which 82.40% were functionally annotated. This chromosome-scale genome would be a valuable reference for assemblies of other complex genomes and studies of evolution in crustaceans.

Tailoring the solar cell efficiency of Y-series based non-fullerene acceptors through end cap modification

Y-series-based non-fullerene acceptors (NFAs) have achieved significant deliberation by chemists and physicists because the promising optical and photochemical properties associated with high-performance OSCs can be further tuned through end-capped modification. In this work, such modifications of Y-series benzothiadiazole-based NFAs were accomplished theoretically to propose new acceptors for photovoltaic cells (PVCs). The recently synthesized Y-series non-fullerene acceptor m-BTP-PhC6 was taken as a reference acceptor. We designed five new acceptors (BTP1-BTP5) through the structural modification at both ends of acceptor groups and evaluated their performance by applying DFT and TD-DFT. The newly engineered molecules exhibited a narrower bandgap (E_g) than the reference (R) resulting in better intramolecular charge transfer (ICT). Further, the designed acceptors expressed the maximum absorption in the region of 600-800 nm revealing a redshift in their absorption spectrum. Low excitation energy and low exciton binding energy were noted for designed acceptors confirming them as better candidates for high PCE of solar cells. Low reorganizational energy for the mobility of holes and electrons was also observed for the designed molecules, indicating improved charge transfer properties. The newly tailored acceptor BTP4 was found to be the promising candidate among all acceptors because of lower bandgap, lower exciton binding energy, reorganizational energy, and redshift of the absorption spectrum. The complex analysis of BTP4 with donor polymer PTB7-Th and PM6 was executed at the same DFT level. Furthermore, FMOs studies showed relatively rich electron density in the acceptor groups of LUMO as compared to the reference molecule. The overall theoretical results of this study showed that the designed acceptors played a productive and effective role in uplifting the efficiency of fullerene-free energy devices.

Designing of anthracene-arylamine hole transporting materials for organic and perovskite solar cells

This study focuses on the creation of 5 small donor molecules (A102W1-A102W5) by substituting the one-sided methoxy group of model (A102R) with different thiophene bridged acceptor moieties. B3LYP/6-31**G (d,p) model has been employed for computational analysis. The best miscibility was found for A102W3 in dichloromethane (DCM) solvent, where its λ_max was also found to be at 753 nm, its E_g was found to be 1.55 eV as well as dipole moment in DCM was 21.47 D. The percentage of PCE among all the variants was greatest for A102W2 (25.31%). The electron reorganization energy shown by A102W4 was 0.00470 eV, whereas the hole reorganization energy investigated in A102W2 was 0.00586 eV representing their maximum electron and hole mobility respectively amongst all. Results validate the value of specified techniques, opening a new door to create efficient small donors for OSCs and HTMs for PSCs.

Whole-genome assembly and annotation for the little yellow croaker (Larimichthys polyactis) provide insights into the evolution of hermaphroditism and gonochorism

The evolutionary direction of gonochorism and hermaphroditism is an intriguing mystery to be solved. The special transient hermaphroditic stage makes the little yellow croaker (Larimichthys polyactis) an appealing model for studying hermaphrodite formation. However, the origin and evolutionary relationship between of L. polyactis and Larimichthys crocea, the most famous commercial fish species in East Asia, remain unclear. Here, we report the sequence of the L. polyactis genome, which we found is ~706 Mb long (contig N50 = 1.21 Mb and scaffold N50 = 4.52 Mb) and contains 25,233 protein-coding genes. Phylogenomic analysis suggested that L. polyactis diverged from the common ancestor, L. crocea, approximately 25.4 million years ago. Our high-quality genome assembly enabled comparative genomic analysis, which revealed several within-chromosome rearrangements and translocations, without major chromosome fission or fusion events between the two species. The dmrt1 gene was identified as the male-specific gene in L. polyactis. Transcriptome analysis showed that the expression of dmrt1 and its upstream regulatory gene (rnf183) were both sexually dimorphic. Rnf183, unlike its two paralogues rnf223 and rnf225, is only present in Larimichthys and Lates but not in other teleost species, suggesting that it originated from lineage-specific duplication or was lost in other teleosts. Phylogenetic analysis shows that the hermaphrodite stage in male L. polyactis may be explained by the sequence evolution of dmrt1. Decoding the L. polyactis genome not only provides insight into the genetic underpinnings of hermaphrodite evolution, but also provides valuable information for enhancing fish aquaculture.

Revealing evolution of tropane alkaloid biosynthesis by analyzing two genomes in the Solanaceae family

Tropane alkaloids (TAs) are widely distributed in the Solanaceae, while some important medicinal tropane alkaloids (mTAs), such as hyoscyamine and scopolamine, are restricted to certain species/tribes in this family. Little is known about the genomic basis and evolution of TAs biosynthesis and specialization in the Solanaceae. Here, we present chromosome-level genomes of two representative mTAs-producing species: Atropa belladonna and Datura stramonium. Our results reveal that the two species employ a conserved biosynthetic pathway to produce mTAs despite being distantly related within the nightshade family. A conserved gene cluster combined with gene duplication underlies the wide distribution of TAs in this family. We also provide evidence that branching genes leading to mTAs likely have evolved in early ancestral Solanaceae species but have been lost in most of the lineages, with A. belladonna and D. stramonium being exceptions. Furthermore, we identify a cytochrome P450 that modifies hyoscyamine into norhyoscyamine. Our results provide a genomic basis for evolutionary insights into the biosynthesis of TAs in the Solanaceae and will be useful for biotechnological production of mTAs via synthetic biology approaches.

Application of a core genome sequence typing (cgMLST) pipeline for surveillance of Clostridioides difficile in China

Introduction

Clostridioides difficile (C. difficile) is a nosocomial bacterial pathogen that causes antibiotic-associated diarrhea mediated by cellular exotoxins secreted into the intestine during bacterial growth. Multilocus sequence typing (MLST) and PCR ribotyping are the main molecular typing for C. difficile. Whole genome sequencing (WGS) core genome multilocus sequence typing (cgMLST) was developed for genetic evolution and outbreak investigation of C. difficile with higher precision and accuracy.

Methods

A total of 699 whole (complete and draft) genome sequences of distinct C. difficile strains were used in this study to identify core gene set (2469 core genes) and the cgMLST scheme for the phylogeny analysis of C. difficile. This cgMLST pipeline was then carried the Chinese Pathogen Identification Net (China PIN) for surveillance of C. difficile in China. Within the China PIN, 195 WGS of C. difficile and an outbreak of CDI with 12 WGS of C. difficile were used to evaluate the cgMLST pipeline.

Results

The result displayed that mostly tested C. difficile isolates could be successfully divided into 5 classic clades and the outbreak event was also successfully identified.

Discussion

The results are meaningful and offer a practicable pipeline for a national-wide surveillance of C. difficile in China.

Comparative Analysis of Olfactory Receptor Repertoires Sheds Light on the Diet Adaptation of the Bamboo-Eating Giant Panda Based on the Chromosome-Level Genome

The giant panda (Ailuropoda melanoleuca) is the epitome of a flagship species for wildlife conservation and also an ideal model of adaptive evolution. As an obligate bamboo feeder, the giant panda relies on the olfaction for food recognition. The number of olfactory receptor (OR) genes and the rate of pseudogenes are the main factors affecting the olfactory ability of animals. In this study, we used the chromosome-level genome of the giant panda to identify OR genes and compared the genome sequences of OR genes with five other Ursidae species (spectacled bear (Tremarctos ornatus), American black bear (Ursus americanus), brown bear (Ursus arctos), polar bear (Ursus maritimus) and Asian black bear (Ursus thibetanus)). The giant panda had 639 OR genes, including 408 functional genes, 94 partial OR genes and 137 pseudogenes. Among them, 222 OR genes were detected and distributed on 18 chromosomes, and chromosome 8 had the most OR genes. A total of 448, 617, 582, 521 and 792 OR genes were identified in the spectacled bear, American black bear, brown bear, polar bear and Asian black bear, respectively. Clustering analysis based on the OR protein sequences of the six species showed that the OR genes distributed in 69 families and 438 subfamilies based on sequence similarity, and the six mammals shared 72 OR gene subfamilies, while the giant panda had 31 unique OR gene subfamilies (containing 35 genes). Among the 35 genes, there are 10 genes clustered into 8 clusters with 10 known human OR genes (OR8J3, OR51I1, OR10AC1, OR1S2, OR1S1, OR51S1, OR4M1, OR4M2, OR51T1 and OR5W2). However, the kind of odor molecules can be recognized by the 10 known human OR genes separately, which needs further research. The phylogenetic tree showed that 345 (about 84.56%) functional OR genes were clustered as Class-II, while only 63 (about 15.44%) functional OR genes were clustered as Class-I, which required further and more in-depth research. The potential odor specificity of some giant panda OR genes was identified through the similarity to human protein sequences. Sequences similar to OR2B1, OR10G3, OR11H6 and OR11H7P were giant panda-specific lacking, which may be related to the transformation and specialization from carnivore to herbivore of the giant panda. Since our reference to flavoring agents comes from human research, the possible flavoring agents from giant panda-specific OR genes need further investigation. Moreover, the conserved motifs of OR genes were highly conserved in Ursidae species. This systematic study of OR genes in the giant panda will provide a solid foundation for further research on the olfactory function and variation of the giant panda.

Chromosome-level genome assembly for takin (Budorcas taxicolor) provides insights into its taxonomic status and genetic diversity

The takin (Budorcas taxicolor) is one of the largest bovid herbivores in the subfamily Caprinae. The takin is at high risk of extinction, but its taxonomic status and genetic diversity remain unclear. In this study, we constructed the first reference genome of Bu. taxicolor using PacBio long High-Fidelity reads and Hi-C technology. The assembled genome is ~2.95 Gb with a contig N50 of 68.05 Mb, which were anchored onto 25+XY chromosomes. We found that the takin was more closely related to muskox than to other Caprinae species. Compared to the common ancestral karyotype of bovidae (2n = 60), we found the takin (2n = 52) experienced four chromosome fusions and one large translocation. Furthermore, we resequenced nine golden takins from the main distribution area, the Qinling Mountains, and identified 3.3 million single nucleotide polymorphisms. The genetic diversity of takin was very low (θπ = 0.00028 and heterozygosity =0.00038), among the lowest detected in domestic and wild mammals. Takin genomes showed a high inbreeding coefficient (FROH =0.217), suggesting severe inbreeding depression. The demographic history showed that the effective population size of takins declined significantly from ~100,000 years ago. Our results provide valuable information for protection of takins and insights into their evolution.

Bamboo Plant Part Preference Affects the Nutrients Digestibility and Intestinal Microbiota of Geriatric Giant Pandas

Bamboo part preference plays a critical role in influencing the nutrient utilization and gastrointestinal microbiota composition of captive giant pandas. However, the effects of bamboo part consumption on the nutrient digestibility and gut microbiome of geriatric giant pandas remain unknown. A total of 11 adult and 11 aged captive giant pandas were provided with bamboo shoots or bamboo leaves in the respective single-bamboo-part consumption period, and the nutrient digestibility and fecal microbiota of both adult and aged giant pandas in each period were evaluated. Bamboo shoot ingestion increased the crude protein digestibility and decreased the crude fiber digestibility of both age groups. The fecal microbiome of the bamboo shoot-fed giant pandas exhibited greater alpha diversity indices and significantly different beta diversity index than the bamboo leaf-fed counterparts regardless of age. Bamboo shoot feeding significantly changed the relative abundance of predominant taxa at both phylum and genus levels in adult and geriatric giant pandas. Bamboo shoot-enriched genera were positively correlated with crude protein digestibility and negatively correlated with crude fiber digestibility. Taken together, these results suggest that bamboo part consumption dominates over age in affecting the nutrient digestibility and gut microbiota composition of giant pandas.

The Gut-Prostate Axis: A New Perspective of Prostate Cancer Biology through the Gut Microbiome

Obesity and a high-fat diet are risk factors associated with prostate cancer, and lifestyle, especially diet, impacts the gut microbiome. The gut microbiome plays important roles in the development of several diseases, such as Alzheimer's disease, rheumatoid arthritis, and colon cancer. The analysis of feces from patients with prostate cancer by 16S rRNA sequencing has uncovered various associations between altered gut microbiomes and prostate cancer. Gut dysbiosis caused by the leakage of gut bacterial metabolites, such as short-chain fatty acids and lipopolysaccharide results in prostate cancer growth. Gut microbiota also play a role in the metabolism of androgen which could affect castration-resistant prostate cancer. Moreover, men with high-risk prostate cancer share a specific gut microbiome and treatments such as androgen-deprivation therapy alter the gut microbiome in a manner that favors prostate cancer growth. Thus, implementing interventions aiming to modify lifestyle or altering the gut microbiome with prebiotics or probiotics may curtail the development of prostate cancer. From this perspective, the "Gut-Prostate Axis" plays a fundamental bidirectional role in prostate cancer biology and should be considered when screening and treating prostate cancer patients.

Next-generation development and application of codon model in evolution

To date, numerous nucleotide, amino acid, and codon substitution models have been developed to estimate the evolutionary history of any sequence/organism in a more comprehensive way. Out of these three, the codon substitution model is the most powerful. These models have been utilized extensively to detect selective pressure on a protein, codon usage bias, ancestral reconstruction and phylogenetic reconstruction. However, due to more computational demanding, in comparison to nucleotide and amino acid substitution models, only a few studies have employed the codon substitution model to understand the heterogeneity of the evolutionary process in a genome-scale analysis. Hence, there is always a question of how to develop more robust but less computationally demanding codon substitution models to get more accurate results. In this review article, the authors attempted to understand the basis of the development of different types of codon-substitution models and how this information can be utilized to develop more robust but less computationally demanding codon substitution models. The codon substitution model enables to detect selection regime under which any gene or gene region is evolving, codon usage bias in any organism or tissue-specific region and phylogenetic relationship between different lineages more accurately than nucleotide and amino acid substitution models. Thus, in the near future, these codon models can be utilized in the field of conservation, breeding and medicine.

Comparative transcriptome analyses of immune responses to LPS in peripheral blood mononuclear cells from the giant panda, human, mouse, and monkey

Gram-negative bacteria are major pathogens that can cause illnesses in giant pandas. Lipopolysaccharides (LPS), components of Gram-negative bacteria, can activate immune responses in mammals (i.e., humans and mice) through recognition by toll-like receptors (TLRs). However, the giant pandas' immune response to LPS stimulation and the differences between the giant panda and other mammals are not fully known. In this study, we administrated peripheral blood mononuclear cells (PBMCs) from giant pandas, humans, C57BL/6 mice, and rhesus monkeys by LPS treatment at 6 h followed by RNA sequencing (RNA-seq), respectively, with control of non-stimulation. KEGG analyses of differentially expressed genes (DEGs) pathways indicated that LPS could activate the classic signaling pathway of NF-κB in PBMCs from those four tested species. Thus, similar to the other three species, NF-κB is an LPS-responsive regulator of innate immune responses in giant pandas. Furthermore, the expression patterns of adapter genes, inflammatory cytokine genes, chemokines, interferon genes, cytokine genes related to cell growth and development, costimulatory molecules, Th1/Th2 cytokine genes, Th17 cytokine genes, Th9, and Th22 cytokine genes were compared among giant pandas and three other species. Our data indicated that in addition to the similar expression patterns of certain genes among giant pandas and other species, the unique expression pattern response to LPS in giant pandas was also discovered. Furthermore, Th9, Th17, and Th22 cells might be involved in the response to LPS in giant pandas at this tested time point. This study reveals that LPS-induced immune responses have different sensitivities and response timelines in giant pandas compared with other mammals. This study facilitates further understanding of the role of the TLR signaling pathway and the immune system in giant pandas, which might be helpful for disease prevention and protection.

Phased Genome Assemblies

The ultimate goal of de novo assembly of reads sequenced from a diploid individual is the separate reconstruction of the sequences corresponding to the two copies of each chromosome. Unfortunately, the allele linkage information needed to perform phased genome assemblies has been difficult to generate. Hence, most current genome assemblies are a haploid mixture of the two underlying chromosome copies present in the sequenced individual. Sequencing technologies providing long (20 kb) and accurate reads are the basis to generate phased genome assemblies. This chapter provides a brief overview of the main milestones in traditional genome assembly, focusing on the bioinformatic techniques developed to generate haplotype information from different specialized protocols. Using these techniques as a knowledge background, the chapter reviews the current algorithms to generate phased assemblies from long reads with low error rates. Current techniques perform haplotype-aware error correction steps to increase the quality of the raw reads. In addition, variations on the traditional overlap-layout-consensus (OLC) graph have been developed in an effort to eliminate edges between reads sequenced from different chromosome copies. This allows for large presence-absence variants between the chromosome copies to be taken into account. The development of these algorithms, along with the improved sequencing technologies has been crucial to finish chromosome-level assemblies of complex genomes.

Dynamic Monitoring of Changes in Fecal Flora of Giant Pandas in Mice: Co-Occurrence Network Reconstruction

Giant pandas are uniquely vulnerable mammals in western China. It is important to develop an animal model to explore the intestinal flora of giant pandas to understand the relationship between digestive diseases and flora. Existing animal models of intestinal flora focus on human flora-associated animals, such as mice, and there is a very limited amount of knowledge regarding giant panda flora-associated animals. To fill this gap, fecal microorganisms from giant pandas were transplanted into pseudosterile and germfree mice using single and multiple gavages. Fecal samples were collected from mice at four time points after transplantation for microbial community analysis. We determined that compared to pseudosterile mice, the characteristics of intestinal flora in pandas were better reproduced in germfree mice. There was no significant difference in microbial diversity between germfree mice and giant panda gut microbes from day 3 to day 21. Germfree mice at the phylum level possessed large amounts of Firmicutes and Proteobacteria, and at the genus level, Escherichia-Shigella, Clostridium sensu stricto 1, and Streptococcus dominated the intestinal flora structure. The microbial community co-occurrence network based on indicator species indicated that germfree mice transplanted with fecal bacteria tended to form a microbial community co-occurrence network similar to that of giant pandas, while pseudosterile mice tended to restore the microbial community co-occurrence network originally present in these mice. Our data are helpful for the study of giant panda flora-associated animals and provide new insights for the in vitro study of giant panda intestinal flora. IMPORTANCE The giant panda is a unique vulnerable mammal in western China, and its main cause of death is digestive system diseases regardless of whether these animals are in the wild or in captivity. The relationship between the intestinal flora and the host exerts a significant impact on the nutrition and health of the giant pandas. However, the protected status of the giant panda has made in vivo, repeatable, and large-sample sampling studies of their intestinal flora difficult. This greatly hinders the research depth of the giant panda intestinal flora from the source. The development and utilization of specific animal models to simulate the structure and characteristics of the intestinal flora provide another means to deal with these research limitations. However, current research examining giant panda flora-associated animals is limited. This study is the first to reveal dynamic changes in the fecal flora of giant pandas in mice after transplantation.

Using Blood Transcriptome Analysis to Determine the Changes in Immunity and Metabolism of Giant Pandas with Age

A low reproductive rate coupled with human activities has endangered the giant panda, a species endemic to southwest China. Although giant pandas feed almost exclusively on bamboo, they retain carnivorous traits and suffer from carnivorous diseases. Additionally, their immune system is susceptible to aging, resulting in a reduced ability to respond to diseases. This study aimed to determine the genes and pathways expressed differentially with age in blood tissues. The differentially expressed genes in different age groups of giant pandas were identified by RNA-seq. The elderly giant pandas had many differentially expressed genes compared with the young group (3 years old), including 548 upregulated genes and 401 downregulated genes. Further, functional enrichment revealed that innate immune upregulation and adaptive immune downregulation were observed in the elderly giant pandas compared with the young giant pandas. Meanwhile, the immune genes in the elderly giant pandas changed considerably, including genes involved in innate immunity and adaptive immunity such as PLSCR1, CLEC7A, CCL5, CCR9, and EPAS1. Time series analysis found that giant pandas store glycogen by prioritizing fat metabolism at age 11, verifying changes in the immune system. The results reported in this study will provide a foundation for further research on disease prevention and the energy metabolism of giant pandas.

Role of feeding specialization in taste receptor loss: insights from sweet and umami receptor evolution in Carnivora

Controversy and misunderstanding surround the role of feeding specialization in taste receptor loss in vertebrates. We refined and tested the hypothesis that this loss is caused by feeding specializations. Specifically, feeding specializations were proposed to trigger time-dependent process of taste receptor loss through deprivation of benefit of using the receptor's gustatory function. We propose that this process may be accelerated by abiotic environmental conditions or decelerated/stopped because of extragustatory functions of the receptor's protein(s). As test case we used evolution of the sweet (TAS1R2+TAS1R3) and umami (TAS1R1+TAS1R3) receptors in Carnivora (dogs, cats, and kin). We predicted these receptors' absence/presence using data on presence/absence of inactivating mutations in these receptors' genes and data from behavioral sweet/umami preference tests. We identified 20 evolutionary events of sweet (11) or umami (9) receptor loss. These events affected species with feeding specializations predicted to favor sweet/umami receptor loss (27 and 22 species, respectively). All species with feeding habits predicted to favor sweet/umami receptor retention (11 and 24, respectively) were found to retain that receptor. Six species retained the sweet (5) or umami (1) receptor despite feeding specialization predicted to favor loss of that receptor, which can be explained by the time dependence of sweet/umami receptor loss process and the possible decelerating effect of TAS1R extragustatory functions so that the sweet/umami receptor process is ongoing in these species. Our findings support the idea that feeding specialization leads to taste receptor loss and is the main if not only triggering factor for evolutionary loss of taste receptors.

Whole genome sequencing of the fast-swimming Southern bluefin tuna (Thunnus maccoyii)

The economically important Southern bluefin tuna (Thunnus maccoyii) is a world-famous fast-swimming fish, but its genomic information is limited. Here, we performed whole genome sequencing and assembled a draft genome for Southern bluefin tuna, aiming to generate useful genetic data for comparative functional prediction. The final genome assembly is 806.54 Mb, with scaffold and contig N50 values of 3.31 Mb and 67.38 kb, respectively. Genome completeness was evaluated to be 95.8%. The assembled genome contained 23,403 protein-coding genes and 236.1 Mb of repeat sequences (accounting for 29.27% of the entire assembly). Comparative genomics analyses of this fast-swimming tuna revealed that it had more than twice as many hemoglobin genes (18) as other relatively slow-moving fishes (such as seahorse, sunfish, and tongue sole). These hemoglobin genes are mainly localized in two big clusters (termed as "MNˮ and "LAˮ respectively), which is consistent with other reported fishes. However, Thr39 of beta-hemoglobin in the MN cluster, conserved in other fishes, was mutated as cysteine in tunas including the Southern bluefin tuna. Since hemoglobins are reported to transport oxygen efficiently for aerobic respiration, our genomic data suggest that both high copy numbers of hemoglobin genes and an adjusted function of the beta-hemoglobin may support the fast-swimming activity of tunas. In summary, we produced a primary genome assembly and predicted hemoglobin-related roles for the fast-swimming Southern bluefin tuna.

An improved, chromosome-level genome of the giant panda (Ailuropoda melanoleuca)

BACKGROUND

The iconic giant panda (Ailuropoda melanoleuca), as both a flagship and umbrella species endemic to China, is a world famous symbol for wildlife conservation. The giant panda has several specific biological traits and holds a relatively small place in evolution. A high-quality genome of the giant panda is key to understanding the biology of this vulnerable species.

FINDINGS

We generated a 2.48-Gb chromosome-level genome (GPv1) of the giant panda named "Jing Jing" with a contig N50 of 28.56 Mb and scaffold N50 of 134.17 Mb, respectively. The total length of chromosomes (n = 21) was 2.39-Gb, accounting for 96.4% of the whole genome. Compared with the previously published four genomes of the giant panda, our genome is characterized by the highest completeness and the correct sequence orientation. A gap-free and 850 kb length of immunoglobulin heavy-chain gene cluster was manually annotated in close proximity to the telomere of chromosome 14. Additionally, we developed an algorithm to predict the centromere position of each chromosome. We also constructed a complete chromatin structure for "Jing Jing", which includes inter-chromosome interaction pattern, A/B compartment, topologically associated domain (TAD), TAD-clique and promoter-enhancer interaction (PEI).

CONCLUSIONS

We presented an improved chromosome-level genome and complete chromatin structure for the giant panda. This is a valuable resource for the future genetic and genomic studies on giant panda.

Chromosome-level genome assembly for the Aldabra giant tortoise enables insights into the genetic health of a threatened population

BACKGROUND

The Aldabra giant tortoise (Aldabrachelys gigantea) is one of only two giant tortoise species left in the world. The species is endemic to Aldabra Atoll in Seychelles and is listed as Vulnerable on the International Union for Conservation of Nature Red List (v2.3) due to its limited distribution and threats posed by climate change. Genomic resources for A. gigantea are lacking, hampering conservation efforts for both wild and ex situpopulations. A high-quality genome would also open avenues to investigate the genetic basis of the species' exceptionally long life span.

FINDINGS

We produced the first chromosome-level de novo genome assembly of A. gigantea using PacBio High-Fidelity sequencing and high-throughput chromosome conformation capture. We produced a 2.37-Gbp assembly with a scaffold N50 of 148.6 Mbp and a resolution into 26 chromosomes. RNA sequencing-assisted gene model prediction identified 23,953 protein-coding genes and 1.1 Gbp of repetitive sequences. Synteny analyses among turtle genomes revealed high levels of chromosomal collinearity even among distantly related taxa. To assess the utility of the high-quality assembly for species conservation, we performed a low-coverage resequencing of 30 individuals from wild populations and two zoo individuals. Our genome-wide population structure analyses detected genetic population structure in the wild and identified the most likely origin of the zoo-housed individuals. We further identified putatively deleterious mutations to be monitored.

CONCLUSIONS

We establish a high-quality chromosome-level reference genome for A. gigantea and one of the most complete turtle genomes available. We show that low-coverage whole-genome resequencing, for which alignment to the reference genome is a necessity, is a powerful tool to assess the population structure of the wild population and reveal the geographic origins of ex situ individuals relevant for genetic diversity management and rewilding efforts.

HSDatabase-a database of highly similar duplicate genes from plants, animals, and algae

Gene duplication is an important evolutionary mechanism capable of providing new genetic material, which in some instances can help organisms adapt to various environmental conditions. Recent studies, for example, have indicated that highly similar duplicate genes (HSDs) are aiding adaptation to extreme conditions via gene dosage. However, for most eukaryotic genomes HSDs remain uncharacterized, partly because they can be hard to identify and categorize efficiently and effectively. Here, we collected and curated HSDs in nuclear genomes from various model animals, land plants and algae and indexed them in an online, open-access sequence repository called HSDatabase. Currently, this database contains 117 864 curated HSDs from 40 distinct genomes; it includes statistics on the total number of HSDs per genome as well as individual HSD copy numbers/lengths and provides sequence alignments of the duplicate gene copies. HSDatabase also allows users to download sequences of gene copies, access genome browsers, and link out to other databases, such as Pfam and Kyoto Encyclopedia of Genes and Genomes. What is more, a built-in Basic Local Alignment Search Tool option is available to conveniently explore potential homologous sequences of interest within and across species. HSDatabase has a user-friendly interface and provides easy access to the source data. It can be used on its own for comparative analyses of gene duplicates or in conjunction with HSDFinder, a newly developed bioinformatics tool for identifying, annotating, categorizing and visualizing HSDs. Database URL: http://hsdfinder.com/database/.

Roles of Species-Specific Legumains in Pathogenicity of the Pinewood Nematode Bursaphelenchus xylophilus

Peptidases are very important to parasites, which have central roles in parasite biology and pathogenesis. In this study, by comparative genome analysis, genome-wide peptidase diversities among plant-parasitic nematodes are estimated. We find that genes encoding cysteine peptidases in family C13 (legumain) are significantly abundant in pine wood nematodes Bursaphelenchus genomes, compared to those in other plant-parasitic nematodes. By phylogenetic analysis, a clade of B. xylophilus-specific legumain is identified. RT-qPCR detection shows that these genes are highly expressed at early stage during the nematode infection process. Utilizing transgene technology, cDNAs of three species-specific legumain were introduced into the Arabidopsis γvpe mutant. Functional complementation assay shows that these B. xylophilus legumains can fully complement the activity of Arabidopsis γVPE to mediate plant cell death triggered by the fungal toxin FB1. Secretory activities of these legumains are experimentally validated. By comparative transcriptome analysis, genes involved in plant cell death mediated by legumains are identified, which enrich in GO terms related to ubiquitin protein transferase activity in category molecular function, and response to stimuli in category biological process. Our results suggest that B. xylophilu-specific legumains have potential as effectors to be involved in nematode-plant interaction and can be related to host cell death.