Subtype distribution and expression of the koala retrovirus in the Japanese zoo koala population

We investigated the proviral copies and RNA expression in koala retrovirus (KoRV)-infected koalas. To ascertain any variation in viral load by institution, age, sex, or body condition score, we quantified KoRV proviral DNA and RNA loads in captive koalas (n = 37) reared in Japanese zoos. All koalas were positive for KoRV genes (pol, LTRs, and env of KoRV-A) in genomic DNA (gDNA), and 91.89% were positive for the pol gene in RNA. In contrast, the distribution rates of KoRV-B, KoRV-C, KoRV-D, and KoRV-F env genes in gDNA were 94.59%, 27.03%, 67.57%, and 54.05%, respectively. A wide inter-individual variation and/or a significant inter-institutional difference in proviral DNA (p < 0.0001) and RNA (p < 0.001) amounts (copies/10³ koala β-actin copies) were observed in Awaji Farm England Hill Zoo koalas, which were obtained from southern koala populations, suggesting exogenous incorporation of KoRV in these koalas. Significant (p < 0.05) age differences were noted in KoRV RNA load (p < 0.05) and median total RNA load (p < 0.001), with loads higher in younger koalas (joeys and juveniles). Thus, the current study provides the distribution of KoRV subtypes in Japanese zoo koala populations and identifies several additional risk factors (sex, age, and body condition) associated with KoRV expression.