Unlocking the story in the swab: A new genotyping assay for the amphibian chytrid fungus Batrachochytrium dendrobatidis

One of the most devastating emerging pathogens of wildlife is the chytrid fungus, Batrachochytrium dendrobatidis (Bd), which affects hundreds of amphibian species around the world. Genomic data from pure Bd cultures have advanced our understanding of Bd phylogenetics, genomic architecture and mechanisms of virulence. However, pure cultures are laborious to obtain and whole‐genome sequencing is comparatively expensive, so relatively few isolates have been genetically characterized. Thus, we still know little about the genetic diversity of Bd in natural systems. The most common noninvasive method of sampling Bd from natural populations is to swab amphibian skin. Hundreds of thousands of swabs have been collected from amphibians around the world, but Bd DNA collected via swabs is often low in quality and/or quantity. In this study, we developed a custom Bd genotyping assay using the Fluidigm Access Array platform to amplify 192 carefully selected regions of the Bd genome. We obtained robust sequence data for pure Bd cultures and field‐collected skin swabs. This new assay has the power to accurately discriminate among the major Bd clades, recovering the basic tree topology previously revealed using whole‐genome data. Additionally, we established a critical value for initial Bd load for swab samples (150 Bd genomic equivalents) above which our assay performs well. By leveraging advances in microfluidic multiplex PCR technology and the globally distributed resource of amphibian swab samples, noninvasive skin swabs can now be used to address critical spatial and temporal questions about Bd and its effects on declining amphibian populations.     

Genomic transition to pathogenicity in chytrid fungi

Understanding the molecular mechanisms of pathogen emergence is central to mitigating the impacts of novel infectious disease agents. The chytrid fungus Batrachochytrium dendrobatidis (Bd) is an emerging pathogen of amphibians that has been implicated in amphibian declines worldwide. Bd is the only member of its clade known to attack vertebrates. However, little is known about the molecular determinants of - or evolutionary transition to - pathogenicity in Bd. Here we sequence the genome of Bd's closest known relative - a non-pathogenic chytrid Homolaphlyctis polyrhiza (Hp). We first describe the genome of Hp, which is comparable to other chytrid genomes in size and number of predicted proteins. We then compare the genomes of Hp, Bd, and 19 additional fungal genomes to identify unique or recent evolutionary elements in the Bd genome. We identified 1,974 Bd-specific genes, a gene set that is enriched for protease, lipase, and microbial effector Gene Ontology terms. We describe significant lineage-specific expansions in three Bd protease families (metallo-, serine-type, and aspartyl proteases). We show that these protease gene family expansions occurred after the divergence of Bd and Hp from their common ancestor and thus are localized to the Bd branch. Finally, we demonstrate that the timing of the protease gene family expansions predates the emergence of Bd as a globally important amphibian pathogen.     

Short-term exposure to warm microhabitats could explain amphibian persistence with Batrachochytrium dendrobatidis

Environmental conditions can alter the outcomes of symbiotic interactions. Many amphibian species have declined due to chytridiomycosis, caused by the pathogenic fungus Batrachochytrium dendrobatidis (Bd), but many others persist despite high Bd infection prevalence. This indicates that Bd's virulence is lower, or it may even be a commensal, in some hosts. In the Australian Wet Tropics, chytridiomycosis extirpated Litoria nannotis from high-elevation rain forests in the early 1990 s. Although the species is recolonizing many sites, no population has fully recovered. Litoria lorica disappeared from all known sites in the early 1990 s and was thought globally extinct, but a new population was discovered in 2008, in an upland dry forest habitat it shares with L. nannotis. All frogs of both species observed during three population censuses were apparently healthy, but most carried Bd. Frogs perch on sun-warmed rocks in dry forest streams, possibly keeping Bd infections below the lethal threshold attained in cooler rain forests. We tested whether short-term elevated temperatures can hamper Bd growth in vitro over one generation (four days). Simulating the temperatures available to frogs on strongly and moderately warmed rocks in dry forests, by incubating cultures at 33°C for one hour daily, reduced Bd growth below that of Bd held at 15°C constantly (representing rain forest habitats). Even small decreases in the exponential growth rate of Bd on hosts may contribute to the survival of frogs in dry forests.     

Susceptibility of amphibians to chytridiomycosis is associated with MHC class II conformation

The pathogenic chytrid fungus Batrachochytrium dendrobatidis (Bd) can cause precipitous population declines in its amphibian hosts. Responses of individuals to infection vary greatly with the capacity of their immune system to respond to the pathogen. We used a combination of comparative and experimental approaches to identify major histocompatibility complex class II (MHC-II) alleles encoding molecules that foster the survival of Bd-infected amphibians. We found that Bd-resistant amphibians across four continents share common amino acids in three binding pockets of the MHC-II antigen-binding groove. Moreover, strong signals of selection acting on these specific sites were evident among all species co-existing with the pathogen. In the laboratory, we experimentally inoculated Australian tree frogs with Bd to test how each binding pocket conformation influences disease resistance. Only the conformation of MHC-II pocket 9 of surviving subjects matched those of Bd-resistant species. This MHC-II conformation thus may determine amphibian resistance to Bd, although other MHC-II binding pockets also may contribute to resistance. Rescuing amphibian biodiversity will depend on our understanding of amphibian immune defence mechanisms against Bd. The identification of adaptive genetic markers for Bd resistance represents an important step forward towards that goal.     

Mimulus is an emerging model system for the integration of ecological and genomic studies

The plant genus Mimulus is rapidly emerging as a model system for studies of evolutionary and ecological functional genomics. Mimulus contains a wide array of phenotypic, ecological and genomic diversity. Numerous studies have proven the experimental tractability of Mimulus in laboratory and field studies. Genomic resources currently under development are making Mimulus an excellent system for determining the genetic and genomic basis of adaptation and speciation. Here, we introduce some of the phenotypic and genetic diversity in the genus Mimulus and highlight how direct genetic studies with Mimulus can address a wide spectrum of ecological and evolutionary questions. In addition, we present the genomic resources currently available for Mimulus and discuss future directions for research. The integration of ecology and genetics with bioinformatics and genome technology offers great promise for exploring the mechanistic basis of adaptive evolution and the genetics of speciation.     

From genotype to phenotype in bovine functional genomics

In past 5 years, the promise that came with genome sequencing has revolutionized the functional genomics research field at unprecedented manner. It would soon know what all known genes do, particularly genes involved in genetic improvement of animal health and increase food animal production. With the availability of full bovine genomic sequence, yet we still have a lot of daunting tasks on 'genotype-to-phenotype problem' particularly about the phenotypic variations and trying to predict what genes are likely to be involved, and improved integrated interactive database. This article outlined and discussed about the current status of bovine functional genomics, recent development in bovine genome databases particularly in annotation of bovine genome, bovine quantitative trait loci database and its potential impact to unveil the from genotype-to-phenotype problem.     

Genetic evidence for a high diversity and wide distribution of endemic strains of the pathogenic chytrid fungus Batrachochytrium dendrobatidis in wild Asian amphibians

Population declines and extinctions of amphibians have been attributed to the chytrid fungus Batrachochytrium dendrobatidis (Bd), especially one globally emerging recombinant lineage (‘Bd‐GPL’). We used PCR assays that target the ribosomal internal transcribed spacer region (ITS) of Bd to determine the prevalence and genetic diversity of Bd in South Korea, where Bd is widely distributed but is not known to cause morbidity or mortality in wild populations. We isolated Korean Bd strains from native amphibians with low infection loads and compared them to known worldwide Bd strains using 19 polymorphic SNP and microsatellite loci. Bd prevalence ranged between 12.5 and 48.0%, in 11 of 17 native Korean species, and 24.7% in the introduced bullfrog Lithobates catesbeianus. Based on ITS sequence variation, 47 of the 50 identified Korean haplotypes formed a group closely associated with a native Brazilian Bd lineage, separated from the Bd‐GPL lineage. However, multilocus genotyping of three Korean Bd isolates revealed strong divergence from both Bd‐GPL and the native Brazilian Bd lineages. Thus, the ITS region resolves genotypes that diverge from Bd‐GPL but otherwise generates ambiguous phylogenies. Our results point to the presence of highly diversified endemic strains of Bd across Asian amphibian species. The rarity of Bd‐GPL‐associated haplotypes suggests that either this lineage was introduced into Korea only recently or Bd‐GPL has been outcompeted by native Bd strains. Our results highlight the need to consider possible complex interactions among native Bd lineages, Bd‐GPL and their associated amphibian hosts when assessing the spread and impact of Bd‐GPL on worldwide amphibian populations.     

Fine mapping of the Bsr1 barley stripe mosaic virus resistance gene in the model grass Brachypodium distachyon

The ND18 strain of Barley stripe mosaic virus (BSMV) infects several lines of Brachypodium distachyon, a recently developed model system for genomics research in cereals. Among the inbred lines tested, Bd3-1 is highly resistant at 20 to 25 °C, whereas Bd21 is susceptible and infection results in an intense mosaic phenotype accompanied by high levels of replicating virus. We generated an F(6:7) recombinant inbred line (RIL) population from a cross between Bd3-1 and Bd21 and used the RILs, and an F(2) population of a second Bd21 × Bd3-1 cross to evaluate the inheritance of resistance. The results indicate that resistance segregates as expected for a single dominant gene, which we have designated Barley stripe mosaic virus resistance 1 (Bsr1). We constructed a genetic linkage map of the RIL population using SNP markers to map this gene to within 705 Kb of the distal end of the top of chromosome 3. Additional CAPS and Indel markers were used to fine map Bsr1 to a 23 Kb interval containing five putative genes. Our study demonstrates the power of using RILs to rapidly map the genetic determinants of BSMV resistance in Brachypodium. Moreover, the RILs and their associated genetic map, when combined with the complete genomic sequence of Brachypodium, provide new resources for genetic analyses of many other traits.     

There Is No Evidence for a Temporal Link between Pathogen Arrival and Frog Extinctions in North-Eastern Australia

Pathogen spread can cause population declines and even species extinctions. Nonetheless, in the absence of tailored monitoring schemes, documenting pathogen spread can be difficult. In the case of worldwide amphibian declines the best present understanding is that the chytrid fungus Batrachochytrium dendrobatidis (Bd) has recently spread, causing amphibian declines and extinction in the process. However, good evidence demonstrating pathogen arrival followed by amphibian decline is rare, and analysis of putative evidence is often inadequate. Here we attempt to examine the relationship between Bd arrival and amphibian decline across north-eastern Australia, using sites where a wave-like pattern of amphibian decline was first noticed and at which intensive research has since been conducted. We develop an analytical framework that allows rigorous estimation of pathogen arrival date, which can then be used to test for a correlation between the time of pathogen arrival and amphibian decline across sites. Our results show that, with the current dataset, the earliest possible arrival date of Bd in north-eastern Australia is completely unresolved; Bd could have arrived immediately before sampling commenced or may have arrived thousands of years earlier, the present data simply cannot say. The currently available data are thus insufficient to assess the link between timing of pathogen arrival and population decline in this part of the world. This data insufficiency is surprising given that there have been decades of research on chytridiomycosis in Australia and that there is a general belief that the link between Bd arrival and population decline is well resolved in this region. The lack of data on Bd arrival currently acts as a major impediment to determining the role of environmental factors in driving the global amphibian declines, and should be a major focus of future research.     

Investigating Differences across Host Species and Scales to Explain the Distribution of the Amphibian Pathogen Batrachochytrium dendrobatidis

Many pathogens infect more than one host species, and clarifying how these different hosts contribute to pathogen dynamics can facilitate the management of pathogens and can lend insight into the functioning of pathogens in ecosystems. In this study, we investigated a suite of native and non-native amphibian hosts of the pathogen Batrachochytrium dendrobatidis (Bd) across multiple scales to identify potential mechanisms that may drive infection patterns in the Colorado study system. Specifically, we aimed to determine if: 1) amphibian populations vary in Bd infection across the landscape, 2) amphibian community composition predicts infection (e.g., does the presence or abundance of any particular species influence infection in others?), 3) amphibian species vary in their ability to produce infectious zoospores in a laboratory infection, 4) heterogeneity in host ability observed in the laboratory scales to predict patterns of Bd prevalence in the landscape. We found that non-native North American bullfrogs (Lithobates catesbeianus) are widespread and have the highest prevalence of Bd infection relative to the other native species in the landscape. Additionally, infection in some native species appears to be related to the density of sympatric L. catesbeianus populations. At the smaller host scale, we found that L. catesbeianus produces more of the infective zoospore stage relative to some native species, but that this zoospore output does not scale to predict infection in sympatric wild populations of native species. Rather, landscape level infection relates most strongly to density of hosts at a wetland as well as abiotic factors. While non-native L. catesbeianus have high levels of Bd infection in the Colorado Front Range system, we also identified Bd infection in a number of native amphibian populations allopatric with L. catesbeianus, suggesting that multiple host species are important contributors to the dynamics of the Bd pathogen in this landscape.     

Amphibian chytrid fungus and ranaviruses in the Northwest Territories, Canada

Pathogens can cause serious declines in host species, and knowing where pathogens associated with host declines occur facilitates understanding host-pathogen ecology. Suspected drivers of global amphibian declines include infectious diseases, with 2 pathogens in particular, Batrachochytrium dendrobatidis (Bd) and ranaviruses, causing concern. We explored the host range and geographic distribution of Bd and ranaviruses in the Taiga Plains ecoregion of the Northwest Territories, Canada, in 2007 and 2008. Both pathogens were detected, greatly extending their known geographic distributions. Ranaviruses were widespread geographically, but found only in wood frogs. In contrast, Bd was found at a single site, but was detected in all 3 species of amphibians in the survey area (wood frogs, boreal chorus frogs, western toads). The presence of Bd in the Northwest Territories is not congruent with predicted distributions based on niche models, even though findings from other studies at northern latitudes are consistent with those same models. Unexpectedly, we also found evidence that swabs routinely used to collect samples for Bd screening detected fewer infections than toe clips. Our use and handling of the swabs was consistent with other studies, and the cause of the apparent lack of integrity of swabs is unknown. The ranaviruses detected in our study were confirmed to be Frog Virus 3 by sequence analysis of a diagnostic 500 bp region of the major capsid protein gene. It is unknown whether Bd or ranaviruses are recent arrivals to the Canadian north. However, the genetic analyses required to answer that question can inform larger debates about the origin of Bd in North America as well as the potential effects of climate change and industrial development on the distributions of these important amphibian pathogens.     

Amphibian skin fungal communities vary across host species and do not correlate with infection by a pathogenic fungus

Amphibian population declines caused by the fungus Batrachochytrium dendrobatidis (Bd) have prompted studies on the bacterial community that resides on amphibian skin. However, studies addressing the fungal portion of these symbiont communities have lagged behind. Using ITS1 amplicon sequencing, we examined the fungal portion of the skin microbiome of temperate and tropical amphibian species currently coexisting with Bd in nature. We assessed cooccurrence patterns between bacterial and fungal OTUs using a subset of samples for which bacterial 16S rRNA gene amplicon data were also available. We determined that fungal communities were dominated by members of the phyla Ascomycota and Basidiomycota, and also by Chytridiomycota in the most aquatic amphibian species. Alpha diversity of the fungal communities differed across host species, and fungal community structure differed across species and regions. However, we did not find a correlation between fungal diversity/community structure and Bd infection, though we did identify significant correlations between Bd and specific OTUs. Moreover, positive bacterial–fungal cooccurrences suggest that positive interactions between these organisms occur in the skin microbiome. Understanding the ecology of amphibian skin fungi, and their interactions with bacteria will complement our knowledge of the factors influencing community assembly and the overall function of these symbiont communities.     

Batrachochytrium dendrobatidis shows high genetic diversity and ecological niche specificity among haplotypes in the Maya Mountains of Belize

The amphibian pathogen Batrachochytrium dendrobatidis (Bd) has been implicated in amphibian declines around the globe. Although it has been found in most countries in Central America, its presence has never been assessed in Belize. We set out to determine the range, prevalence, and diversity of Bd using quantitative PCR (qPCR) and sequencing of a portion of the 5.8 s and ITS1-2 regions. Swabs were collected from 524 amphibians of at least 26 species in the protected areas of the Maya Mountains of Belize. We sequenced a subset of 72 samples that had tested positive for Bd by qPCR at least once; 30 samples were verified as Bd. Eight unique Bd haplotypes were identified in the Maya Mountains, five of which were previously undescribed. We identified unique ecological niches for the two most broadly distributed haplotypes. Combined with data showing differing virulence shown in different strains in other studies, the 5.8 s - ITS1-2 region diversity found in this study suggests that there may be substantial differences among populations or haplotypes. Future work should focus on whether specific haplotypes for other genomic regions and possibly pathogenicity can be associated with haplotypes at this locus, as well as the integration of molecular tools with other ecological tools to elucidate the ecology and pathogenicity of Bd.     

Ecological genomics in full colour

Colour patterns in animals have long offered an opportunity to observe adaptive traits in natural populations. Colour plays myriad roles in interactions within and among species, from reproductive signalling to predator avoidance, leading to multiple targets of natural and sexual selection and opportunities for diversification. Understanding the genetic and developmental underpinnings of variation in colour promises a fuller understanding of these evolutionary processes, but the path to unravelling these connections can be arduous. The advent of genomic techniques suitable for nonmodel organisms is now beginning to light the way. Two new studies in this issue of Molecular Ecology use genomic sequencing of laboratory crosses to map colour traits in cichlid fishes, a remarkably diverse group in which coloration has played a major role in diversification. They illustrate how genomic approaches, specifically RAD sequencing, can rapidly identify both simple and more complex genetic variation underlying ecologically important traits. In the first, Henning et al. ( 2014 ) detect a single locus that appears to control in a Mendelian fashion the presence of horizontal stripes, a trait that has evolved in numerous cichlid lineages. In the second, Albertson et al. ( 2014 ) identify several genes and epistatic interactions affecting multiple colour traits, as well as a novel metric describing integration across colour traits. Albertson et al. ( 2014 ) go further, by quantifying differential expression of parental alleles at a candidate locus and by relating differentiation among natural populations at mapped loci to trait divergence. Herein lies the promise of ecological genomics – efficiently integrating genetic mapping of phenotypes with population genomic data to both identify functional genes and unravel their evolutionary history. These studies offer guidance on how genomic techniques can be tailored to a research question or study system, and they also add to the growing body of empirical examples addressing basic questions about how ecologically important traits evolve in natural populations.     

Modelling the amphibian chytrid fungus spread by connectivity analysis: towards a national monitoring network in Italy

The emerging amphibian disease, Batrachochytrium dendrobatidis (Bd), is driving population declines worldwide and even species extinctions in Australia, South and Central America. In order to mitigate effects of Bd on amphibian populations, high-exposed areas should be identified at the local scale and effective conservation measures should be planned at the national level. This assessment is actually lacking in the Mediterranean basin, and in particular in Italy, one of the most relevant amphibian diversity hotspots in the entire region. In this study, we reviewed the available information on Bd in Italy, and conducted a 5-year molecular screening on 1274 individual skin swabs belonging to 18 species. Overall, we found presence of Bd in 13 species and in a total of 56 known occurrence locations for peninsular Italy and Sardinia. We used these occurrence locations and climate data to model habitat suitability of Bd for current and future climatic scenarios. We then employed electric circuit theory to model landscape permeability to the diffusion of Bd, using a resistance map. With this procedure, we were able to model, for the first time, the diffusion pathways of Bd at the landscape scale, characterising the main future pathways towards areas with a high probability of Bd occurrence. Thus, we identified six national protected areas that will become pivotal for a nationally-based strategic plan in order to monitor, mitigate and possibly contrast Bd diffusion in Italy.

     

Patterns of relatedness and genetic diversity inferred from whole genome sequencing of archival blood fluke miracidia (Schistosoma japonicum)

Genomic approaches hold great promise for resolving unanswered questions about transmission patterns and responses to control efforts for schistosomiasis and other neglected tropical diseases. However, the cost of generating genomic data and the challenges associated with obtaining sufficient DNA from individual schistosome larvae (miracidia) from mammalian hosts have limited the application of genomic data for studying schistosomes and other complex macroparasites. Here, we demonstrate the feasibility of utilizing whole genome amplification and sequencing (WGS) to analyze individual archival miracidia. As an example, we sequenced whole genomes of 22 miracidia from 11 human hosts representing two villages in rural Sichuan, China, and used these data to evaluate patterns of relatedness and genetic diversity. We also down-sampled our dataset to test how lower coverage sequencing could increase the cost effectiveness of WGS while maintaining power to accurately infer relatedness. Collectively, our results illustrate that population-level WGS datasets are attainable for individual miracidia and represent a powerful tool for ultimately providing insight into overall genetic diversity, parasite relatedness, and transmission patterns for better design and evaluation of disease control efforts.     

A guide to the genomics of ecological speciation in natural animal populations

Interest in ecological speciation is growing, as evidence accumulates showing that natural selection can lead to rapid divergence between subpopulations. However, whether and how ecological divergence can lead to the buildup of reproductive isolation remains under debate. What is the relative importance of natural selection vs. neutral processes? How does adaptation generate reproductive isolation? Can ecological speciation occur despite homogenizing gene flow? These questions can be addressed using genomic approaches, and with the rapid development of genomic technology, will become more answerable in studies of wild populations than ever before. In this article, we identify open questions in ecological speciation theory and suggest useful genomic methods for addressing these questions in natural animal populations. We aim to provide a practical guide for ecologists interested in incorporating genomic methods into their research programs. An increased integration between ecological research and genomics has the potential to shed novel light on the origin of species.     

Effects of nutrient supplementation on host‐pathogen dynamics of the amphibian chytrid fungus: a community approach

相似文献