Comparative genomics of koala,cattle and sheep strains of Chlamydia pecorum

Background

Chlamydia pecorum is an important pathogen of domesticated livestock including sheep, cattle and pigs. This pathogen is also a key factor in the decline of the koala in Australia. We sequenced the genomes of three koala C. pecorum strains, isolated from the urogenital tracts and conjunctiva of diseased koalas. The genome of the C. pecorum VR629 (IPA) strain, isolated from a sheep with polyarthritis, was also sequenced.

Results

Comparisons of the draft C. pecorum genomes against the complete genomes of livestock C. pecorum isolates revealed that these strains have a conserved gene content and order, sharing a nucleotide sequence similarity > 98%. Single nucleotide polymorphisms (SNPs) appear to be key factors in understanding the adaptive process. Two regions of the chromosome were found to be accumulating a large number of SNPs within the koala strains. These regions include the Chlamydia plasticity zone, which contains two cytotoxin genes (toxA and toxB), and a 77 kbp region that codes for putative type III effector proteins. In one koala strain (MC/MarsBar), the toxB gene was truncated by a premature stop codon but is full-length in IPTaLE and DBDeUG. Another five pseudogenes were also identified, two unique to the urogenital strains C. pecorum MC/MarsBar and C. pecorum DBDeUG, respectively, while three were unique to the koala C. pecorum conjunctival isolate IPTaLE. An examination of the distribution of these pseudogenes in C. pecorum strains from a variety of koala populations, alongside a number of sheep and cattle C. pecorum positive samples from Australian livestock, confirmed the presence of four predicted pseudogenes in koala C. pecorum clinical samples. Consistent with our genomics analyses, none of these pseudogenes were observed in the livestock C. pecorum samples examined. Interestingly, three SNPs resulting in pseudogenes identified in the IPTaLE isolate were not found in any other C. pecorum strain analysed, raising questions over the origin of these point mutations.

Conclusions

The genomic data revealed that variation between C. pecorum strains were mainly due to the accumulation of SNPs, some of which cause gene inactivation. The identification of these genetic differences will provide the basis for further studies to understand the biology and evolution of this important animal pathogen.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-667) contains supplementary material, which is available to authorized users.     

Quantitative PCR method for detection of mycoplasma spp. DNA in nasal lavage samples from the desert tortoise (Gopherus agassizii)

Mycoplasma agassizii and M. testudineum have been associated with upper respiratory tract disease (URTD) in the threatened desert tortoise (Gopherus agassizii). Because microbiological culture methods have proven difficult to employ in wild desert tortoises, our goal was to develop a sensitive and specific qPCR method for detecting and quantifying mycoplasma DNA in nasal lavage fluid collected in the field. Primers for 16S ribosomal RNA gene sequences specific for M. agassizii and M. testudineum were designed, together with primers that recognize conserved sequences of both microorganisms. Standard curves generated with DNA extracted from known numbers of mycoplasma cells revealed a lower detection limit of approximately 5 fg. The qPCR method did not recognize normal flora DNA, and nasal lavage fluid contained no interfering substances. Nasal lavage samples collected from 20 captive desert tortoises housed at the Desert Tortoise Conservation Center (Clark County, Nevada, USA) revealed the presence of M. agassizii DNA in 100% of the tortoises. Concentrations ranged from a low of 6 pg ml^− 1 to a high of 72,962 pg ml^− 1. Only one of the tortoises was positive for M. testudineum. Interestingly, not all of the qPCR positive tortoises showed evidence of seroconversion, suggesting that they were colonized but not infected. This new quantitative method will provide a critical tool for managing threatened populations of the desert tortoise.     

Analysis of partial 16S rRNA nucleotide sequences of Chlamydia pecorum and C. psittaci

Partial 16S nucleotide sequences of Chlamydia psittaci isolates S26/3 (abortion), P94/1 (pigeon) and Chlamydia pecorum isolates W73 (enteric) and E58 (encephalomyelitis) were determined. Analysis of these data indicates very high levels of interspecies sequence conservation, with C. psittaci being more closely related to C. pecorum than to C. pneumoniae or C. trachomatis. Restriction enzyme analysis of nucleotide sequences indicated that BslI can be used to clearly distinguish C. psittaci and C. pecorum isolates. Psittacine and non-psittacine (pigeon) avian isolates of C. psittaci were distinguished using MaeI.     

Western blot can distinguish natural and acquired antibodies to Mycoplasma agassizii in the desert tortoise (Gopherus agassizii)

Mycoplasma agassizi has been identified as a cause of upper respiratory tract disease (URTD) in the threatened Mojave population of the desert tortoise (Gopherus agassizii), and anti-M. agassizii antibodies have been found by ELISA in as many as 15% of these animals across their geographic range. Here we report that a cohort of 16 egg-reared desert tortoises never exposed to M. agassizii had ELISA antibody titers to this organism that overlapped with titers obtained from some M. agassizii-infected tortoises. These natural antibodies were predominantly of the IgM class. Western blots of plasma from these non-infected tortoises produced a characteristic banding pattern against M. agassizii antigens. A group of 38 wild-caught desert tortoises was tested by ELISA, and although some of these tortoises had antibody titers significantly higher than the non-infected tortoises, there was considerable overlap at the lower titer levels. However, Western blot analysis revealed distinct banding patterns that could readily distinguish between the non-infected tortoises and tortoises with acquired antibodies, regardless of ELISA antibody titers. We conclude that desert tortoises have natural antibodies to M. agassizii that can compromise the determination of infection status by ELISA. However, the Western blot technique can distinguish between natural and acquired antibody patterns and can be used to confirm the diagnosis of M. agassizii infections in the desert tortoise.     

Chlamydia buteonis,a new Chlamydia species isolated from a red-shouldered hawk

Chlamydiaceae are obligate intracellular bacterial pathogens for humans and animals. A recent study highlighted that a Chlamydiaceae intermediary between C. psittaci and C. abortus can infect hawks. Here, an isolate was obtained upon passage of cloacal and conjunctival sac material collected from a female hatch-year red-shouldered hawk (Buteo lineatus) in cultured cells. The diseased bird, one of 12 birds housed in a rehabilitation center, developed conjunctivitis and later died. Swabs from both sites tested positive for Chlamydia using the QuickVue Chlamydia test. The isolate, named RSHA, tested negative in qPCR assays specific for C. psittaci and C. abortus, respectively. Analysis of the 16S rRNA, 23S rRNA and whole genome sequences as well as MLST, ANIb and TETRA values reveal that C. psittaci and C. abortus are the closest relatives of RSHA. However, the overall results strongly suggest a phylogenetic intermediate position between these two species. Therefore, we propose the introduction of a new species designated Chlamydia buteonis with RSHA^T as the type strain.     

Genome sequencing and comparative analysis of three Chlamydia pecorum strains associated with different pathogenic outcomes

Background

Chlamydia pecorum is the causative agent of a number of acute diseases, but most often causes persistent, subclinical infection in ruminants, swine and birds. In this study, the genome sequences of three C. pecorum strains isolated from the faeces of a sheep with inapparent enteric infection (strain W73), from the synovial fluid of a sheep with polyarthritis (strain P787) and from a cervical swab taken from a cow with metritis (strain PV3056/3) were determined using Illumina/Solexa and Roche 454 genome sequencing.

Results

Gene order and synteny was almost identical between C. pecorum strains and C. psittaci. Differences between C. pecorum and other chlamydiae occurred at a number of loci, including the plasticity zone, which contained a MAC/perforin domain protein, two copies of a >3400 amino acid putative cytotoxin gene and four (PV3056/3) or five (P787 and W73) genes encoding phospholipase D. Chlamydia pecorum contains an almost intact tryptophan biosynthesis operon encoding trpABCDFR and has the ability to sequester kynurenine from its host, however it lacks the genes folA, folKP and folB required for folate metabolism found in other chlamydiae. A total of 15 polymorphic membrane proteins were identified, belonging to six pmp families. Strains possess an intact type III secretion system composed of 18 structural genes and accessory proteins, however a number of putative inc effector proteins widely distributed in chlamydiae are absent from C. pecorum. Two genes encoding the hypothetical protein ORF663 and IncA contain variable numbers of repeat sequences that could be associated with persistence of infection.

Conclusions

Genome sequencing of three C. pecorum strains, originating from animals with different disease manifestations, has identified differences in ORF663 and pseudogene content between strains and has identified genes and metabolic traits that may influence intracellular survival, pathogenicity and evasion of the host immune system.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-23) contains supplementary material, which is available to authorized users.     

The Role of Endangered Species Reintroduction in Ecosystem Restoration: Tortoise–Cactus Interactions on Española Island, Galápagos

We evaluated the role that endangered species reintroduction efforts can play in the larger context of ecosystem restoration. To do so, we examined interactions between endangered giant tortoises (Geochelone nigra hoodensis), currently being reintroduced to Isla Española, Galápagos, and an arboreal cactus (Opuntia megasperma var. megasperma), which is itself endangered and a keystone resource for many animals on the island. We collected information on spatial patterns of occurrence of cacti, tortoises, and woody vegetation and compared recruitment of juvenile cacti in areas occupied versus unoccupied by tortoises. Reintroduced tortoises appeared to suppress cactus recruitment near the few remaining adult cacti at the study site, but facilitate it at longer distances, with tortoise–cactus interactions mediated by the presence of woody vegetation, which likely alters tortoise movements and thereby patterns of cactus seed dispersal. The net effect of tortoises on cacti appeared to be positive insofar as tortoise presence was associated with greater recruitment of juveniles into cactus populations. Our study provides support for reintroducing endangered reptiles and other animals to aid ecosystem restoration in areas where they might once have played an important role in grazing upon and dispersing plants.     

Biogeography of Parasitic Nematode Communities in the Galápagos Giant Tortoise: Implications for Conservation Management

The Galápagos giant tortoise is an icon of the unique, endemic biodiversity of Galápagos, but little is known of its parasitic fauna. We assessed the diversity of parasitic nematode communities and their spatial distributions within four wild tortoise populations comprising three species across three Galápagos islands, and consider their implication for Galápagos tortoise conservation programmes. Coprological examinations revealed nematode eggs to be common, with more than 80% of tortoises infected within each wild population. Faecal samples from tortoises within captive breeding centres on Santa Cruz, Isabela and San Cristobal islands also were examined. Five different nematode egg types were identified: oxyuroid, ascarid, trichurid and two types of strongyle. Sequencing of the 18S small-subunit ribosomal RNA gene from adult nematodes passed with faeces identified novel sequences indicative of rhabditid and ascaridid species. In the wild, the composition of nematode communities varied according to tortoise species, which co-varied with island, but nematode diversity and abundance were reduced or altered in captive-reared animals. Evolutionary and ecological factors are likely responsible for the variation in nematode distributions in the wild. This possible species/island-parasite co-evolution has not been considered previously for Galápagos tortoises. We recommend that conservation efforts, such as the current Galápagos tortoise captive breeding/rearing and release programme, be managed with respect to parasite biogeography and host-parasite co-evolutionary processes in addition to the biogeography of the host.     

Are captive tortoises a reservoir for conservation? An assessment of genealogical affiliation of captive Gopherus agassizii to local, wild populations

The conservation of tortoises poses a unique situation because several threatened species are commonly kept as pets within their native ranges. Thus, there is potential for captive populations to be a reservoir for repatriation efforts. We assess the utility of captive populations of the threatened Agassiz’s desert tortoise (Gopherus agassizii) for recovery efforts based on genetic affinity to local areas. We collected samples from 130 captive desert tortoises from three desert communities: two in California (Ridgecrest and Joshua Tree) and the Desert Tortoise Conservation Center (Las Vegas) in Nevada. We tested all samples for 25 short tandem repeats and sequenced 1,109 bp of the mitochondrial genome. We compared captive genotypes to a database of 1,258 Gopherus samples, including 657 wild caught G. agassizii spanning the full range of the species. We conducted population assignment tests to determine the genetic origins of the captive individuals. For our total sample set, only 44 % of captive individuals were assigned to local populations based on genetic units derived from the reference database. One individual from Joshua Tree, California, was identified as being a Morafka’s desert tortoise, G. morafkai, a cryptic species which is not native to the Mojave Desert. Our data suggest that captive desert tortoises kept within the native range of G. agassizii cannot be presumed to have a genealogical affiliation to wild tortoises in their geographic proximity. Precautions should be taken before considering the release of captive tortoises into the wild as a management tool for recovery.     

Asymptomatic Endemic Chlamydia pecorum Infections Reduce Growth Rates in Calves by up to 48 Percent

Intracellular Chlamydia (C.) bacteria cause in cattle some acute but rare diseases such as abortion, sporadic bovine encephalomyelitis, kerato-conjunctivitis, pneumonia, enteritis and polyarthritis. More frequent, essentially ubiquitous worldwide, are low-level, asymptomatic chlamydial infections in cattle. We investigated the impact of these naturally acquired infections in a cohort of 51 female Holstein and Jersey calves from birth to 15 weeks of age. In biweekly sampling, we measured blood/plasma markers of health and infection and analyzed their association with clinical appearance and growth in dependence of chlamydial infection intensity as determined by mucosal chlamydial burden or contemporaneous anti-chlamydial plasma IgM. Chlamydia 23S rRNA gene PCR and ompA genotyping identified only C. pecorum (strains 1710S, Maeda, and novel strain Smith3v8) in conjunctival and vaginal swabs. All calves acquired the infection but remained clinically asymptomatic. High chlamydial infection associated with reduction of body weight gains by up to 48% and increased conjunctival reddening (P<10⁻⁴). Simultaneously decreased plasma albumin and increased globulin (P<10⁻⁴) suggested liver injury by inflammatory mediators as mechanisms for the growth inhibition. This was confirmed by the reduction of plasma insulin like growth factor-1 at high chlamydial infection intensity (P<10⁻⁴). High anti-C. pecorum IgM associated eight weeks later with 66% increased growth (P = 0.027), indicating a potential for immune protection from C. pecorum-mediated growth depression. The worldwide prevalence of chlamydiae in livestock and their high susceptibility to common feed-additive antibiotics suggests the possibility that suppression of chlamydial infections may be a major contributor to the growth promoting effect of feed-additive antibiotics.     

The Dominance of Introduced Plant Species in the Diets of Migratory Galapagos Tortoises Increases with Elevation on a Human‐Occupied Island

The distribution of resources and food selection are fundamental to the ecology, life history, physiology, population dynamics, and conservation of animals. Introduced plants are changing foraging dynamics of herbivores in many ecosystems often with unknown consequences. Galapagos tortoises, like many herbivores, undertake migrations along elevation gradients driven by variability in vegetation productivity which take them into upland areas dominated by introduced plants. We sought to characterize diet composition of two species of Galapagos tortoises, focussing on how the role of introduced forage species changes over space and the implications for tortoise conservation. We quantified the distribution of tortoises with elevation using GPS telemetry. Along the elevation gradient, we quantified the abundance of introduced and native plant species, estimated diet composition by recording foods consumed by tortoises, and assessed tortoise physical condition from body weights and blood parameter values. Tortoises ranged between 0 and 429 m in elevation over which they consumed at least 64 plant species from 26 families, 44 percent of which were introduced species. Cover of introduced species and the proportion of introduced species in tortoise diets increased with elevation. Introduced species were positively selected for by tortoises at all elevations. Tortoise physical condition was either consistent or increased with elevation at the least biologically productive season on Galapagos. Santa Cruz tortoises are generalist herbivores that have adapted their feeding behavior to consume many introduced plant species that has likely made a positive contribution to tortoise nutrition. Some transformed habitats that contain an abundance of introduced forage species are compatible with tortoise conservation.     

Nutrition and husbandry conditions of Palearctic tortoises (Testudo spp.) in captivity

Mediterranean and Russian tortoises (Testudo spp.) are popular companion animals (pets), despite ongoing controversy concerning privately keeping reptiles. The arguments used during these controversial discussions have often been based on outdated facts. Therefore, a survey was developed to evaluate the current population structure, husbandry conditions, diet regime, and health status of Testudo species in captivity. More than 75% of the 1075 respondents housed their tortoises in an outdoor enclosure containing a greenhouse or cold frame, which is considered the most species-appropriate way of husbandry. Of the respondents, 67.7% fed their tortoises with the optimum diet of more than 80% grasses and weeds during the summer vegetation period. Only 8.2% of respondents owned a tortoise with a diagnosed disease. According to the results, the likelihood of tortoises developing pyramidal growth syndrome, which can be used as an indicator of the quality of tortoise husbandry, was high in tortoises kept in a terrarium and/or fed a diet of less than 80% grasses and weeds in summer. This likelihood varied among species, with a higher incidence in Hermann’s tortoises (Testudo hermanni).     

Multi Locus Sequence Typing of Chlamydia Reveals an Association between Chlamydia psittaci Genotypes and Host Species

Chlamydia comprises a group of obligate intracellular bacterial parasites responsible for a variety of diseases in humans and animals, including several zoonoses. Chlamydia trachomatis causes diseases such as trachoma, urogenital infection and lymphogranuloma venereum with severe morbidity. Chlamydia pneumoniae is a common cause of community-acquired respiratory tract infections. Chlamydia psittaci, causing zoonotic pneumonia in humans, is usually hosted by birds, while Chlamydia abortus, causing abortion and fetal death in mammals, including humans, is mainly hosted by goats and sheep. We used multi-locus sequence typing to asses the population structure of Chlamydia. In total, 132 Chlamydia isolates were analyzed, including 60 C. trachomatis, 18 C. pneumoniae, 16 C. abortus, 34 C. psittaci and one of each of C. pecorum, C. caviae, C. muridarum and C. felis. Cluster analyses utilizing the Neighbour-Joining algorithm with the maximum composite likelihood model of concatenated sequences of 7 housekeeping fragments showed that C. psittaci 84/2334 isolated from a parrot grouped together with the C. abortus isolates from goats and sheep. Cluster analyses of the individual alleles showed that in all instances C. psittaci 84/2334 formed one group with C. abortus. Moving 84/2334 from the C. psittaci group to the C. abortus group resulted in a significant increase in the number of fixed differences and elimination of the number of shared mutations between C. psittaci and C. abortus. C. psittaci M56 from a muskrat branched separately from the main group of C. psittaci isolates. C. psittaci genotypes appeared to be associated with host species. The phylogentic tree of C. psittaci did not follow that of its host bird species, suggesting host species jumps. In conclusion, we report for the first time an association between C. psittaci genotypes with host species.     

Low‐coverage reduced representation sequencing reveals subtle within‐island genetic structure in Aldabra giant tortoises

Aldabrachelys gigantea (Aldabra giant tortoise) is one of only two giant tortoise species left in the world and survives as a single wild population of over 100,000 individuals on Aldabra Atoll, Seychelles. Despite this large current population size, the species faces an uncertain future because of its extremely restricted distribution range and high vulnerability to the projected consequences of climate change. Captive‐bred A. gigantea are increasingly used in rewilding programs across the region, where they are introduced to replace extinct giant tortoises in an attempt to functionally resurrect degraded island ecosystems. However, there has been little consideration of the current levels of genetic variation and differentiation within and among the islands on Aldabra. As previous microsatellite studies were inconclusive, we combined low‐coverage and double‐digest restriction‐associated DNA (ddRAD) sequencing to analyze samples from 33 tortoises (11 from each main island). Using 5426 variant sites within the tortoise genome, we detected patterns of within‐island population structure, but no differentiation between the islands. These unexpected results highlight the importance of using genome‐wide genetic markers to capture higher‐resolution genetic structure to inform future management plans, even in a seemingly panmictic population. We show that low‐coverage ddRAD sequencing provides an affordable alternative approach to conservation genomic projects of non‐model species with large genomes.     

Tracing Genetic Lineages of Captive Desert Tortoises in Arizona

Abstract: We genotyped 180 captive desert tortoises (Gopherus agassizii) from Kingman (n = 45), Phoenix (n = 113), and Tucson (n = 22), Arizona, USA, to determine if the genetic lineage of captives is associated with that of wild tortoises in the local area (Sonoran Desert). We tested all samples for 16 short tandem repeats and sequenced 1,109 base pairs of mitochondrial DNA (mtDNA). To determine genetic origin, we performed assignment tests against a reference database of 997 desert tortoise samples collected throughout the Mojave and Sonoran Deserts. We found that >40% of our Arizona captive samples were genetically of Mojave Desert or hybrid origin, with the percentage of individuals exhibiting the Mojave genotype increasing as the sample locations approached the California, USA, border. In Phoenix, 11.5% were Sonoran–Mojave crosses, and 8.8% were hybrids between desert tortoise and Texas tortoise (G. berlandieri). Our findings present many potential implications for wild tortoises in the Sonoran Desert of Arizona. Escaped or released captive tortoises with Mojave or hybrid genotypes have the potential to affect the genetic composition of Sonoran wild populations. Genotyping captive desert tortoises could be used to inform the adoption process, and thereby provide additional protection to native desert-tortoise populations in Arizona.     

Novel molecular markers of <Emphasis Type="Italic">Chlamydia pecorum</Emphasis> genetic diversity in the koala (<Emphasis Type="Italic">Phascolarctos cinereus</Emphasis>)

Background  

Chlamydia pecorum is an obligate intracellular bacterium and the causative agent of reproductive and ocular disease in several animal hosts including koalas, sheep, cattle and goats. C. pecorum strains detected in koalas are genetically diverse, raising interesting questions about the origin and transmission of this species within koala hosts. While the ompA gene remains the most widely-used target in C. pecorum typing studies, it is generally recognised that surface protein encoding genes are not suited for phylogenetic analysis and it is becoming increasingly apparent that the ompA gene locus is not congruent with the phylogeny of the C. pecorum genome. Using the recently sequenced C. pecorum genome sequence (E58), we analysed 10 genes, including ompA, to evaluate the use of ompA as a molecular marker in the study of koala C. pecorum genetic diversity.     

The origin of Indian Star tortoises (Geochelone elegans) based on nuclear and mitochondrial DNA analysis: A story of rescue and repatriation

The Indian Star tortoise (Geochelone elegans) belongs to the family Testunidae and is distributed in southwest India and Sri Lanka. In addition to facing loss of its natural habitat, the species is also illegally traded as food and as an exotic pet internationally. Here we report DNA-based analyses for identification and repatriation of these tortoises into their natural habitat. We have attempted to establish the geographical origin of these tortoises rescued from smugglers, by comparing the microsatellite and mitochondrial markers of rescued animals with animals of known provenance. Star tortoises exhibited strong genetic structure in India. The populations from western India were genetically distinct at microsatellite and mitochondrial loci from southern populations. The rescued individuals had similar multilocus genotypes and mitochondrial DNA haplotypes as the reference individuals from south India. However, the precise geographic origin of many of the rescued samples remains unresolved, because we could not assign them to southern populations and the Neighbor-Joining cluster analysis indicated that some of rescued tortoises formed distinct clusters. These data strongly suggest that the rescued group of tortoises is composed of a mix of individuals from differentiated source populations that are probably located in southern India and possibly Sri Lanka. Our study provides valuable information based on molecular markers for the assessment of genetic diversity in Indian Star tortoises.     

Body size development of captive and free‐ranging Leopard tortoises (Geochelone pardalis)

The growth and weight development of Leopard tortoise hatchings (Geochelone pardalis) kept at the Al Wabra Wildlife Preservation (AWWP), Qatar, was observed for more than four years, and compared to data in literature for free‐ranging animals on body weight or carapace measurements. The results document a distinctively faster growth in the captive animals. Indications for the same phenomenon in other tortoise species (Galapagos giant tortoises, G. nigra; Spur‐thighed tortoises, Testudo graeca; Desert tortoises, Gopherus agassizi) were found in the literature. The cause of the high growth rate most likely is the constant provision with highly digestible food of low fiber content. Increased growth rates are suspected to have negative consequences such as obesity, high mortality, gastrointestinal illnesses, renal diseases, “pyramiding,” fibrous osteodystrophy or metabolic bone disease. The apparently widespread occurrence of high growth rates in intensively managed tortoises underlines how easily ectothermic animals can be oversupplemented with nutrients. Zoo Biol 29:517–525, 2010. © 2009 Wiley‐Liss, Inc.     

Potential Facilitation Between a Commensal and a Pathogenic Microbe in a Wildlife Disease

We assessed the potential for microbial interactions influencing a well-documented host–pathogen system. Mycoplasma agassizii is the known etiological agent of upper respiratory tract disease in Mojave desert tortoises (Gopherus agassizii), but disease in wild animals is extremely heterogeneous. For example, a much larger proportion of animals harbor M. agassizii than those that develop disease. With the availability of a new quantitative PCR assay for a microbe that had previously been implicated in disease, Pasteurella testudinis, we tested 389 previously collected samples of nasal microbes from tortoise populations across the Mojave desert. We showed that P. testudinis is a common commensal microbe. However, we did find that its presence was associated with higher levels of M. agassizii among the tortoises positive for this pathogen. The best predictor of P. testudinis prevalence in tortoise populations was average size of tortoises, suggesting that older populations have higher levels of P. testudinis. The prevalence of co-infection in populations was associated with the prevalence of URTD, providing additional evidence for an indirect interaction between the two microbes and inflammatory disease. We showed that URTD, like many chronic, polymicrobial diseases involving mucosal surfaces, shows patterns of a polymicrobial etiology.