Genome of the Giant Panda Roundworm Illuminates Its Host Shift and Parasitic Adaptation

Baylisascaris schroederi, a roundworm (ascaridoid) parasite specific to the bamboo-feeding giant panda (Ailuropoda melanoleuca), represents a leading cause of mortality in wild giant panda populations. Here, we present a 293-megabase chromosome-level genome assembly of B. schroederi to infer its biology, including host adaptations. Comparative genomics revealed an evolutionary trajectory accompanied by host-shift events in ascaridoid parasite lineages after host separations, suggesting their potential for transmission and rapid adaptation to new hosts. Genomic and anatomical lines of evidence, including expansion and positive selection of genes related to the cuticle and basal metabolisms, indicate that B. schroederi undergoes specific adaptations to survive in the sharp-edged bamboo-enriched gut of giant pandas by structurally increasing its cuticle thickness and efficiently utilizing host nutrients through gut parasitism. Additionally, we characterized the secretome of B. schroederi and predicted potential drug and vaccine targets for new control strategies. Overall, this genome resource provides new insights into the host adaptation of B. schroederi to the giant panda as well as the host-shift events in ascaridoid parasite lineages. Our findings on the unique biology of B. schroederi will also aid in the development of prevention and treatment measures to protect giant panda populations from roundworm parasitism.     

Major histocompatibility complex alleles associated with parasite susceptibility in wild giant pandas

Major histocompatibility complex (MHC) polymorphism is thought to be driven by antagonistic coevolution between pathogens and hosts, mediated through either overdominance or frequency-dependent selection. However, investigations under natural conditions are still rare for endangered mammals which often exhibit depleted variation, and the mechanism of selection underlying the maintenance of characteristics remains a considerable debate. In this study, 87 wild giant pandas were used to investigate MHC variation associated with parasite load. With the knowledge of the MHC profile provided by the genomic data of the giant panda, seven DRB1, seven DQA1 and eight DQA2 alleles were identified at each single locus. Positive selection evidenced by a significantly higher number of non-synonymous substitutions per non-synonymous codon site relative to synonymous substitutions per synonymous codon site could only be detected at the DRB1 locus, which leads to the speculation that DRB1 may have a more important role in dealing with parasite infection for pandas. Coprological analyses revealed that 55.17% of individuals exhibited infection with 1–2 helminthes and 95.3% of infected pandas carried Baylisascaris shroederi. Using a generalized linear model, we found that Aime-DRB1*10 was significantly associated with parasite infection, but no resistant alleles could be detected. MHC heterozygosity of the pandas was found to be uncorrelated with the infection status or the infection intensity. These results suggested that the possible selection mechanisms in extant wild pandas may be frequency dependent rather than being determined by overdominance selection. Our findings could guide the candidate selection for the ongoing reintroduction or translocation of pandas.     

Diet drives convergent evolution of gut microbiomes in bamboo-eating species

Gut microbiota plays a critical role in host physiology and health. The coevolution between the host and its gut microbes facilitates animal adaptation to its specific ecological niche. Multiple factors such as host diet and phylogeny modulate the structure and function of gut microbiota. However, the relative contribution of each factor in shaping the structure of gut microbiota remains unclear. The giant(Ailuropoda melanoleuca) and red(Ailurus styani) pandas belong to different families of order Carnivora. They have evolved as obligate bamboo-feeders and can be used as a model system for studying the gut microbiome convergent evolution. Here, we compare the structure and function of gut microbiota of the two pandas with their carnivorous relatives using 16S rRNA and metagenome sequencing. We found that both panda species share more similarities in their gut microbiota structure with each other than each species shares with its carnivorous relatives. This indicates that the specialized herbivorous diet rather than host phylogeny is the dominant driver of gut microbiome convergence within Arctoidea.Metagenomic analysis revealed that the symbiotic gut microbiota of both pandas possesses a high level of starch and sucrose metabolism and vitamin B12 biosynthesis. These findings suggest a diet-driven convergence of gut microbiomes and provide new insight into host-microbiota coevolution of these endangered species.     

Profile of microRNA in Giant Panda Blood: A Resource for Immune-Related and Novel microRNAs

The giant panda (Ailuropoda melanoleuca) is one of the world’s most beloved endangered mammals. Although the draft genome of this species had been assembled, little was known about the composition of its microRNAs (miRNAs) or their functional profiles. Recent studies demonstrated that changes in the expression of miRNAs are associated with immunity. In this study, miRNAs were extracted from the blood of four healthy giant pandas and sequenced by Illumina next generation sequencing technology. As determined by miRNA screening, a total of 276 conserved miRNAs and 51 novel putative miRNAs candidates were detected. After differential expression analysis, we noticed that the expressions of 7 miRNAs were significantly up-regulated in young giant pandas compared with that of adults. Moreover, 2 miRNAs were up-regulated in female giant pandas and 1 in the male individuals. Target gene prediction suggested that the miRNAs of giant panda might be relevant to the expressions of 4,602 downstream genes. Subseuqently, the predicted target genes were conducted to KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis and we found that these genes were mainly involved in host immunity, including the Ras signaling pathway, the PI3K-Akt signaling pathway, and the MAPK signaling pathway. In conclusion, our results provide the first miRNA profiles of giant panda blood, and the predicted functional analyses may open an avenue for further study of giant panda immunity.     

Coevolution between parasite virulence and host life-history traits

Epidemiological models generally explore the evolution of parasite life-history traits, namely, virulence and transmission, against a background of constant host life-history traits. However, life-history models have predicted the evolution of host traits in response to parasitism. The coevolution of host and parasite life-history traits remains largely unexplored. We present an epidemiological model, based on resource allocation theory, that provides an analysis of the coevolution between host reproductive effort and parasite virulence. This model allows for hosts with either a fixed (i.e., genetic) or conditional (i.e., a phenotypically plastic) response to parasitism. It also considers superinfections. We show that parasitism always favors increased allocation to host reproduction, but because of epidemiological feedbacks, the evolutionarily stable host reproductive effort does not always increase with parasite virulence. Superinfection drives the evolution of parasite virulence and acts on the evolution of the host through parasite evolution, generally leading to higher host reproductive effort. Coevolution, as opposed to cases where only one of the antagonists evolves, may generate correlations between host and parasite life-history traits across environmental gradients affecting the fecundity or the survival of the host. Our results provide a theoretical framework against which experimental coevolution outcomes or field observations can be contrasted.     

The soybean cyst nematode, Heterodera glycines: a genetic model system for the study of plant-parasitic nematodes

Despite advances in understanding plant responses to nematode infection, little information exists regarding parasitic mechanisms. Recently, it has become possible to perform genetic analysis of soybean cyst nematode. Integration of classic and reverse genetics and genomic approaches for the parasite, with host genetics and genomics will expand our knowledge of nematode parasitism.     

Relatedness and spatial proximity as determinants of host–parasite interactions in the brood parasitic Barrow's goldeneye (Bucephala islandica)

Recent studies, which have found evidence for kin-biased egg donation, have sparked interest in re-assessing the parasitic nature of conspecific brood parasitism (CBP). Since host–parasite kinship is essential for mutual benefits to arise from CBP, we explored the role of relatedness in determining the behaviour of conspecific nest parasites and their hosts in nesting female Barrow's goldeneyes ( Bucephala islandica ), a duck in which CBP is common. The results revealed that the amount of parasitism increased with host–parasite relatedness, the effect of which was independent of geographical proximity of host and parasite nests. Proximity per se was also positively associated with the amount of parasitism. Furthermore, while hosts appeared to reduce their clutch size as a response to the presence of parasitic eggs, the magnitude of host clutch reduction also tended to increase with increasing relatedness to the parasite. Hence, our results indicate that both relatedness and spatial proximity are important determinants of CBP, and that host clutch reduction may be an adaptation to nest parasitism, modulated by host–parasite relatedness. Taken together, the results provide a demonstration that relatedness influences host and parasite behaviour in Barrow's goldeneyes, resulting in kin-biased egg donation.     

Bacterial pore-forming proteins as anthelmintics

Crystal (Cry) proteins are made by the Gram-positive bacterium Bacillus thuringiensis (Bt). Cry proteins are pore-forming proteins and are the most widely used biological insecticides in the world. Our laboratory found some Cry proteins are highly effective against a broad range of nematodes (roundworms). Here, we discuss our results of Cry protein activity against intestinal roundworms. Both Cry5B and Cry21A have therapeutic activities against infections of the roundworm Heligmosomoides polygyrus bakeri in mice. Cry5B also shows highly therapeutic activity against Ancylostoma ceylanicum infection in hamsters. A. ceylanicum is a minor hookworm parasite of humans, and it is closely related to the more prevalent Ancylostoma duodenale. In addition, Cry proteins show excellent combinatorial therapeutic properties with nicotinic acetylcholine receptor (nAChR) agonists, one of the two classes of compounds approved by the World Health Organization for the treatment for intestinal roundworms in humans. Given their non-toxicity to humans and their broad spectrum of nematicidal action, Cry proteins show great potential as next-generation anthelmintics.     

Subtilase activity in intrusive cells mediates haustorium maturation in parasitic plants

Parasitic plants that infect crops are devastating to agriculture throughout the world. These parasites develop a unique inducible organ called the haustorium that connects the vascular systems of the parasite and host to establish a flow of water and nutrients. Upon contact with the host, the haustorial epidermal cells at the interface with the host differentiate into specific cells called intrusive cells that grow endophytically toward the host vasculature. Following this, some of the intrusive cells re-differentiate to form a xylem bridge (XB) that connects the vasculatures of the parasite and host. Despite the prominent role of intrusive cells in host infection, the molecular mechanisms mediating parasitism in the intrusive cells remain poorly understood. In this study, we investigated differential gene expression in the intrusive cells of the facultative parasite Phtheirospermum japonicum in the family Orobanchaceae by RNA-sequencing of laser-microdissected haustoria. We then used promoter analyses to identify genes that are specifically induced in intrusive cells, and promoter fusions with genes encoding fluorescent proteins to develop intrusive cell-specific markers. Four of the identified intrusive cell-specific genes encode subtilisin-like serine proteases (SBTs), whose biological functions in parasitic plants are unknown. Expression of SBT inhibitors in intrusive cells inhibited both intrusive cell and XB development and reduced auxin response levels adjacent to the area of XB development. Therefore, we propose that subtilase activity plays an important role in haustorium development in P. japonicum.

Subtilases specifically expressed in intrusive cells regulate auxin-mediated host–parasite connections in the parasitic plant Phtheirospermum japonicum.     

Geographic pattern of antibiotic resistance genes in the metagenomes of the giant panda

The rise in infections by antibiotic-resistant bacteria poses a serious public health problem worldwide. The gut microbiome of animals is a reservoir for antibiotic resistance genes (ARGs). However, the correlation between the gut microbiome of wild animals and ARGs remains controversial. Here, based on the metagenomes of giant pandas (including three wild populations from the Qinling, Qionglai and Xiaoxiangling Mountains, and two major captive populations from Yaan and Chengdu), we investigated the potential correlation between the constitution of the gut microbiome and the composition of ARGs across the different geographic locations and living environments. We found that the types of ARGs were correlated with gut microbiome composition. The NMDS cluster analysis using Jaccard distance of the ARGs composition of the gut microbiome of wild giant pandas displayed a difference based on geographic location. Captivity also had an effect on the differences in ARGs composition. Furthermore, we found that the Qinling population exhibited profound dissimilarities of both gut microbiome composition and ARGs (the highest proportion of Clostridium and vancomycin resistance genes) when compared to the other wild and captive populations studies, which was supported by previous giant panda whole-genome sequencing analysis. In this study, we provide an example of a potential consensus pattern regarding host population genetics, symbiotic gut microbiome and ARGs. We revealed that habitat isolation impacts the ARG structure in the gut microbiome of mammals. Therefore, the difference in ARG composition between giant panda populations will provide some basic information for their conservation and management, especially for captive populations.     

Strong specificity in the interaction between parasitoids and symbiont‐protected hosts

Coevolution between hosts and parasites may promote the maintenance of genetic variation in both antagonists by negative frequency‐dependence if the host–parasite interaction is genotype‐specific. Here we tested for specificity in the interaction between parasitoids (Lysiphlebus fabarum) and aphid hosts (Aphis fabae) that are protected by a heritable defensive endosymbiont, the γ‐proteobacterium Hamiltonella defensa. Previous studies reported a lack of genotype specificity between unprotected aphids and parasitoids, but suggested that symbiont‐conferred resistance might exhibit a higher degree of specificity. Indeed, in addition to ample variation in host resistance as well as parasitoid infectivity, we found a strong aphid clone‐by‐parasitoid line interaction on the rates of successful parasitism. This genotype specificity appears to be mediated by H. defensa, highlighting the important role that endosymbionts can play in host–parasite coevolution.     

Cuckoos, cowbirds and hosts: adaptations, trade-offs and constraints

The interactions between brood parasitic birds and their host species provide one of the best model systems for coevolution. Despite being intensively studied, the parasite-host system provides ample opportunities to test new predictions from both coevolutionary theory as well as life-history theory in general. I identify four main areas that might be especially fruitful: cuckoo female gentes as alternative reproductive strategies, non-random and nonlinear risks of brood parasitism for host individuals, host parental quality and targeted brood parasitism, and differences and similarities between predation risk and parasitism risk. Rather than being a rare and intriguing system to study coevolutionary processes, I believe that avian brood parasites and their hosts are much more important as extreme cases in the evolution of life-history strategies. They provide unique examples of trade-offs and situations where constraints are either completely removed or particularly severe.     

Egg color variation,but not egg rejection behavior,changes in a cuckoo host breeding in the absence of brood parasitism

Interactions between parasitic cuckoos and their songbird hosts form a classical reciprocal “arms race,” and are an excellent model for understanding the process of coevolution. Changes in host egg coloration via the evolution of interclutch variation in egg color or intraclutch consistency in egg color are hypothesized counter adaptations that facilitate egg recognition and thus limit brood parasitism. Whether these antiparasitism strategies are maintained when the selective pressure of parasitism is relaxed remains debated. However, introduced species provide unique opportunities for testing the direction and extent of natural selection on phenotypic trait maintenance and variation. Here, we investigated egg rejection behavior and egg color polymorphism in the red‐billed leiothrix (Leiothrix lutea), a common cuckoo (Cuculus canorus) host, in a population introduced to Hawaii 100 years ago (breeding without cuckoos) and a native population in China (breeding with cuckoos). We found that egg rejection ability was equally strong in both the native and the introduced populations, but levels of interclutch variation and intraclutch consistency in egg color in the native population were higher than in the introduced population. This suggests that egg rejection behavior in hosts can be maintained in the absence of brood parasitism and that egg appearance is maintained by natural selection as a counter adaptation to brood parasitism. This study provides rare evidence that host antiparasitism strategies can change under parasite‐relaxed conditions and reduced selection pressure.     

Similar yet different: co-analysis of the genetic diversity and structure of an invasive nematode parasite and its invasive mammalian host

Animal parasitic nematodes can cause serious diseases and their emergence in new areas can be an issue of major concern for biodiversity conservation and human health. Their ability to adapt to new environments and hosts is likely to be affected by their degree of genetic diversity, with gene flow between distinct populations counteracting genetic drift and increasing effective population size. The raccoon roundworm (Baylisascaris procyonis), a gastrointestinal parasite of the raccoon (Procyon lotor), has increased its global geographic range after being translocated with its host. The raccoon has been introduced multiple times to Germany, but not all its populations are infected with the parasite. While fewer introduced individuals may have led to reduced diversity in the parasite, admixture between different founder populations may have counteracted genetic drift and bottlenecks. Here, we analyse the population genetic structure of the roundworm and its raccoon host at the intersection of distinct raccoon populations infected with B. procyonis. We found evidence for two parasite clusters resulting from independent introductions. Both clusters exhibited an extremely low genetic diversity, suggesting small founding populations subjected to inbreeding and genetic drift with no, or very limited, genetic influx from population admixture. Comparison of the population genetic structures of both host and parasite suggested that the parasite spread to an uninfected raccoon founder population. On the other hand, an almost perfect match between cluster boundaries also suggested that the population genetic structure of B. procyonis has remained stable since its introduction, mirroring that of its raccoon host.     

基于线粒体12S对中国两大山系大熊猫西氏贝蛔虫种群遗传多样性的分析

大熊猫是中国特有的珍稀濒危物种,而西氏贝蛔虫是大熊猫体内最为常见的肠道寄生虫,对野生和圈养大熊猫危害极大。考虑到大熊猫因分布区域的差异而形成了不同的亚种以及寄生虫与宿主间广泛的协同进化关系,西氏贝蛔虫是否也存在与大熊猫相适应的亚种分化一直是野生动物学家极其关注和热议的话题。为此,本文选择中国两大山系（岷山和邛崃）大熊猫种群体内共计34株西氏贝蛔虫虫体样本进行种群遗传多态性研究。利用PCR技术扩增出了岷山（14株）和邛崃（20株）西氏贝蛔虫的线粒体12S基因全序列并对其做了遗传多样性分析。结果表明：（1）34个样本包含9个单倍型,呈现出一个高单倍性多样性和低核苷酸多样性的特点;（2）负的Tajima's D和Fu's Fs中性检验值及“多峰型”的种群歧点分布图暗示种群不久前曾经历过突增长的现象;（3）低的种群间的分化系数和高的基因流表明两个地理种群间未形成显著的遗传分化;（4）系统发育树和单倍型网络图表明两山系种群分布无区域特异性。因此,岷山和邛崃山系的大熊猫体内的西氏贝蛔虫种群遗传变异性较低,分化不明显。该发现不仅暗示了西氏贝蛔虫与其宿主（大熊猫）的进化不同步,而且还为不同区域大熊猫西氏贝蛔虫病的监控提供了理论依据。     

Genomic evidence for population‐specific responses to co‐evolving parasites in a New Zealand freshwater snail

Reciprocal co‐evolving interactions between hosts and parasites are a primary source of strong selection that can promote rapid and often population‐ or genotype‐specific evolutionary change. These host–parasite interactions are also a major source of disease. Despite their importance, very little is known about the genomic basis of co‐evolving host–parasite interactions in natural populations, especially in animals. Here, we use gene expression and sequence evolution approaches to take critical steps towards characterizing the genomic basis of interactions between the freshwater snail Potamopyrgus antipodarum and its co‐evolving sterilizing trematode parasite, Microphallus sp., a textbook example of natural coevolution. We found that Microphallus‐infected P. antipodarum exhibit systematic downregulation of genes relative to uninfected P. antipodarum. The specific genes involved in parasite response differ markedly across lakes, consistent with a scenario where population‐level co‐evolution is leading to population‐specific host–parasite interactions and evolutionary trajectories. We also used an F_ST‐based approach to identify a set of loci that represent promising candidates for targets of parasite‐mediated selection across lakes as well as within each lake population. These results constitute the first genomic evidence for population‐specific responses to co‐evolving infection in the P. antipodarum‐Microphallus interaction and provide new insights into the genomic basis of co‐evolutionary interactions in nature.     

The evolutionary significance of parasitism: do parasite‐driven genetic dynamics occur ex silico?

It has long been recognized that reciprocal antagonism might lock host and parasite populations into a process of constant change, adapting and reacting in open‐ended coevolution. A significant body of theory supports this intuition: dynamic genetic polymorphisms are a common outcome of computer simulations of host–parasite coevolution. These in silico experiments have also shown that dynamical interactions could be responsible for high levels of genetic diversity in host populations, and even be the principle determinant of rates of genetic recombination and sexuality. The evolutionary significance of parasitism depends on the strength and prevalence of parasite‐mediated selection in nature. Here I appraise whether parasitism is a pervasive agent of evolutionary change by detailing empirical evidence for selection. Although there is considerable evidence of genetic variation for resistance, and hence the potential for selection, direct observation of parasite‐driven genetic change is lacking.