Abstract
Reptarenaviruses cause boid inclusion body disease (BIBD), a fatal disease particularly impacting captive boa constrictor collections. The development of cytoplasmic inclusion bodies (IBs) comprising reptarenavirus nucleoprotein (NP) in many cell types of affected snakes is characteristic of BIBD. However, snakes can harbor reptarenaviruses without showing IBs, hence representing carriers and a potential source of transmission. The RNA genome of reptarenaviruses comprises a small (S) and a large (L) segment, and the snakes with BIBD commonly carry a swarm of reptarenavirus segments. To design sensitive and reliable tools for the diagnosis of reptarenavirus infection in snake colonies, we used metatranscriptomics to determine the reptarenavirus segments present in a large boa constrictor breeding colony. The analysis identified one reptarenavirus S segment and three L segments in the colony. The sequence data served to design real-time reverse transcription-PCR (RT-PCR) targeting the found S segment. This allowed us to identify all infected animals and to quantify the S segment RNA levels, which we found to correlate with the presence of IBs. We further found a positive correlation between the number of L segments and the S segment RNA level, which could suggest that L segment excess also contributes to the IB formation. Information on cohousing of the snakes showed a clear association of reptarenavirus infection with cohousing in general and cohousing with infected animals. Information on breeding and offspring confirmed that vertical transmission occurred. Furthermore, our data suggest that some animals might be able to clear the infection or at least exhibit transient or intermittent viremia. IMPORTANCE Boid inclusion body disease (BIBD) is caused by reptarenavirus infection, and while reptarenavirus nucleoprotein is the main component of the inclusion bodies (IBs) characteristic of BIBD, not all reptarenavirus-infected snakes demonstrate IBs in their cells. Identification of infected individuals is critical for controlling the spread of the disease; however, the genetic divergence of reptarenaviruses complicates reverse transcription-PCR (RT-PCR)-based diagnostics. Here, we tested a next-generation-sequencing-based approach to establish a tailored "colony-specific" set of diagnostic tools for the detection of reptarenavirus small (S) and large (L) genome segments. With this approach, we could demonstrate that an S-segment-specific RT-PCR is highly effective in identifying the infected individuals. We further found the S segment RNA level to positively correlate with the presence of IBs and the number of L segments, which could direct future studies to identify the BIBD pathogenetic mechanisms.