The evolution of brain neuron numbers in amniotes.

SignificanceThe evolution of brain processing capacity has traditionally been inferred from data on brain size. However, similarly sized brains of distantly related species can differ in the number and distribution of neurons, their basic computational units. Therefore, a finer-grained approach is needed to reveal the evolutionary paths to increased cognitive capacity. Using a new, comprehensive dataset, we analyzed brain cellular composition across amniotes. Compared to reptiles, mammals and birds have dramatically increased neuron numbers in the telencephalon and cerebellum, which are brain parts associated with higher cognition. Astoundingly, a phylogenetic analysis suggests that as few as four major changes in neuron-brain scaling in over 300 million years of evolution pave the way to intelligence in endothermic land vertebrates.

Identification of potentially zoonotic parasites in captive orangutans and semi-captive mandrills: Phylogeny and morphological comparison.

Cysts and trophozoites of vestibuliferid ciliates and larvae of Strongyloides were found in fecal samples from captive orangutans Pongo pygmaeus and P. abelii from Czech and Slovak zoological gardens. As comparative material, ciliates from semi-captive mandrills Mandrillus sphinx from Gabon were included in the study. Phylogenetic analysis of the detected vestibuliferid ciliates using ITS1-5.8s-rRNA-ITS2 and partial 18S ribosomal deoxyribonucleic acid (rDNA) revealed that the ciliates from orangutans are conspecific with Balantioides coli lineage A, while the ciliates from mandrills clustered with Buxtonella-like ciliates from other primates. Morphological examination of the cysts and trophozoites using light microscopy did not reveal differences robust enough to identify the genera of the ciliates. Phylogenetic analysis of detected L1 larvae of Strongyloides using partial cox1 revealed Strongyloides stercoralis clustering within the cox1 lineage A infecting dogs, humans, and other primates. The sequences of 18S rDNA support these results. As both B. coli and S. stercoralis are zoonotic parasites and the conditions in captive and semi-captive settings may facilitate transmission to humans, prophylactic measures should reflect the findings.

GLOBAL RETROSPECTIVE REVIEW OF SEVERE ACUTE RESPIRATORY SYNDROME SARS COV-2 INFECTIONS IN NONDOMESTIC FELIDS: MARCH 2020-FEBRUARY 2021.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections in nondomestic felids have been documented in North America, South America, Africa, Europe, and Asia. Between March 2020 and February 2021, at nine institutions across three continents, infection was confirmed in 16 tigers (Panthera tigris), 14 lions (Panthera leo), three snow leopards (Panthera uncia), one cougar (Puma concolor), and one Amur leopard cat (Prionailurus bengalensis euptilurus) ranging from 2 to 21 yr old (average, 10 yr). Infection was suspected in an additional 12 tigers, 4 lions, and 9 cougars. Clinical signs (in order of most to least common) included coughing, ocular and/or nasal discharge, wheezing, sneezing, decreased appetite, lethargy, diarrhea, and vomiting. Most felids recovered uneventfully, but one geriatric tiger with comorbidities developed severe dyspnea and neurologic signs necessitating euthanasia. Clinical signs lasted 1-19 d (average, 8 d); one tiger was asymptomatic. Infection was confirmed by various methods, including antigen tests and/or polymerase chain reaction (PCR) of nasal or oral swabs, tracheal wash, and feces, or virus isolation from feces or tracheal wash. Infection status and resolution were determined by testing nasal swabs from awake animals, fecal PCR, and observation of clinical signs. Shedding of fecal viral RNA was significantly longer than duration of clinical signs. Postinfection seropositivity was confirmed by four institutions including 11 felids (5 lions, 6 tigers). In most instances, asymptomatic or presymptomatic keepers were the presumed or confirmed source of infection, although in some instances the infection source remains uncertain. Almost all infections occurred despite using cloth facemasks and disposable gloves when in proximity to the felids and during food preparation. Although transmission may have occurred during momentary lapses in personal protective equipment compliance, it seems probable that cloth masks are insufficient at preventing transmission of SARS-CoV-2 from humans to nondomestic felids. Surgical or higher grade masks may be warranted when working with nondomestic felids.

Reverse-zoonotic transmission of SARS-CoV-2 lineage alpha (B.1.1.7) to great apes and exotic felids in a zoo in the Czech Republic.

We report an outbreak of SARS-CoV-2 lineage alpha in gorillas and felid species in a zoo in Prague, Czech Republic. The course of illness and clinical signs are described, as are the results of characterization of these particular SARS-CoV-2 variants by next-generation sequencing and phylogenetic analysis. The putative transmission routes are also discussed.

Bioacoustics as a tool to monitor the estrus cycle in a female slow loris (Nycticebus sp.).

Evidence-based husbandry leading to increased reproductive success and strengthening of ex situ populations of slow lorises (Nycticebus spp.) living in zoos is highly important. Better fulfillment of their social needs is one of the main priorities in achieving these objectives. We performed 21-month long acoustic monitoring of a zoo-kept female slow loris (Nycticebus sp.) housed either singly or with one of two males to examine whether her estrus cycle potentially could be detected based on her vocal activity. We found a regular cycle of remarkably increased whistle production that lasted approximately 31.2 days in a nocturnal exhibit and approximately 39 days in an off-show room. The regular cycle of increased vocal activity corresponded to a previously described estrus cycle of slow lorises and was observable in the presence of both males as well as when the female was housed singly. Additionally, vaginal smears collected from the female close to the peak period of her vocal activity showed signs of proestrus and estrus. The acoustic properties of the whistles, specifically that they did not overlap with or were loud enough to exceed background noise commonly occurring in zoos, made them perfect candidates for analyses involving automatic processing of a large number of recordings. We conclude that bioacoustics represents a promising, completely noninvasive and relatively easily applicable tool that allows detection and anticipation of the estrus cycle in some females, thus improving the social management of slow lorises living in zoos.

Genetic diversity, evolution and selection in the major histocompatibility complex DRB and DQB loci in the family Equidae.

BACKGROUND: The mammalian Major Histocompatibility Complex (MHC) is a genetic region containing highly polymorphic genes with immunological functions. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. The MHC class II sub-region contains genes expressed in antigen presenting cells. The antigen binding site is encoded by the second exon of genes encoding antigen presenting molecules. The exon 2 sequences of these MHC genes have evolved under the selective pressure of pathogens. Interspecific differences can be observed in the class II sub-region. The family Equidae includes a variety of domesticated, and free-ranging species inhabiting a range of habitats exposed to different pathogens and represents a model for studying this important part of the immunogenome. While equine MHC class II DRA and DQA loci have received attention, the genetic diversity and effects of selection on DRB and DQB loci have been largely overlooked. This study aimed to provide the first in-depth analysis of the MHC class II DRB and DQB loci in the Equidae family. RESULTS: Three DRB and two DQB genes were identified in the genomes of all equids. The genes DRB2, DRB3 and DQB3 showed high sequence conservation, while polymorphisms were more frequent at DRB1 and DQB1 across all species analyzed. DQB2 was not found in the genome of the Asiatic asses Equus hemionus kulan and E. h. onager. The bioinformatic analysis of non-zero-coverage-bases of DRB and DQB genes in 14 equine individual genomes revealed differences among individual genes. Evidence for recombination was found for DRB1, DRB2, DQB1 and DQB2 genes. Trans-species allele sharing was identified in all genes except DRB1. Site-specific selection analysis predicted genes evolving under positive selection both at DRB and DQB loci. No selected amino acid sites were identified in DQB3. CONCLUSIONS: The organization of the MHC class II sub-region of equids is similar across all species of the family. Genomic sequences, along with phylogenetic trees suggesting effects of selection as well as trans-species polymorphism support the contention that pathogen-driven positive selection has shaped the MHC class II DRB/DQB sub-regions in the Equidae.

Comparative Study of the Bush Dog (Speothos venaticus) Karyotype and Analysis of Satellite DNA Sequences and Their Chromosome Distribution in Six Species of Canidae.

The bush dog (Speothos venaticus, 2n = 74) is a near threatened species taxonomically classified among South American canids. We revised the bush dog karyotype and performed a comparative sequence analysis of satellite and satellite-like DNAs in 6 canids: the bush dog, domestic dog (Canis familiaris, 2n = 78), grey wolf (C. lupus, 2n = 78), Chinese raccoon dog (Nyctereutes procyonoides procyonoides, 2n = 54+B), red fox (Vulpes vulpes, 2n = 34+B), and arctic fox (V. lagopus, 2n = 48-50) to specify the species position among Canidae. Using FISH with painting and BAC probes, we found that the distribution of canid evolutionarily conserved chromosome segments in the bush dog karyotype is similar to that of the domestic dog and grey wolf. The bush dog karyotype differs by 2 acrocentric chromosome pairs formed by tandem fusions of the canine (29;34) and (26;35) orthologues. An interstitial signal of the telomeric probe was observed in the (26;35) fusion site in the bush dog indicating a recent evolutionary origin of this rearrangement. Sequences and hybridisation patterns of satellite DNAs were compared, and a phylogenetic tree of the 6 canid species was constructed which confirmed the bush dog position close to the wolf-like canids, and apart from the raccoon dog and foxes.

Seroprevalence of Neospora caninum and Toxoplasma gondii in exotic ruminants and camelids in the Czech Republic.

Neospora caninum and Toxoplasma gondii are the protozoan parasites with definitive hosts from order Carnivora. Due to vertical transmission, both parasites can cause abortions and neonatal mortality that lead to significant productive and economic losses in the domestic ruminants. The aim of this study was to describe N. caninum and T. gondii seroprevalence in the group of frequently farmed captive exotic ruminants (n = 184) including Bovidae (barbary sheep, bezoar goat, common eland, American bison, water buffalo, and yak) and Camelidae (bactrian camel, guanaco, llama, and alpaca). Antibodies were tested by indirect fluorescent antibody test (IFAT) and enzyme-linked immunosorbent assay (ELISA). Higher prevalence of T. gondii antibodies (31% in IFAT and 52% in ELISA) was detected compared to N. caninum (24% in IFAT and 17% in cELISA). Mixed infection was found in 18 (10%) and 22 (12%) animals by IFAT and ELISA, respectively. Higher seroprevalence of both N. caninum and T. gondii was found in Camelidae compared to Bovidae. To author knowledge, this is the first detection of T. gondii and N. caninum antibodies in common elands and bezoar goats.

First Description of the Karyotype and Sex Chromosomes in the Komodo Dragon (Varanus komodoensis).

The Komodo dragon (Varanus komodoensis) is the largest lizard in the world. Surprisingly, it has not yet been cytogenetically examined. Here, we present the very first description of its karyotype and sex chromosomes. The karyotype consists of 2n = 40 chromosomes, 16 macrochromosomes and 24 microchromosomes. Although the chromosome number is constant for all species of monitor lizards (family Varanidae) with the currently reported karyotype, variability in the morphology of the macrochromosomes has been previously documented within the group. We uncovered highly differentiated ZZ/ZW sex microchromosomes with a heterochromatic W chromosome in the Komodo dragon. Sex chromosomes have so far only been described in a few species of varanids including V. varius, the sister species to Komodo dragon, whose W chromosome is notably larger than that of the Komodo dragon. Accumulations of several microsatellite sequences in the W chromosome have recently been detected in 3 species of monitor lizards; however, these accumulations are absent from the W chromosome of the Komodo dragon. In conclusion, although varanids are rather conservative in karyotypes, their W chromosomes exhibit substantial variability at the sequence level, adding further evidence that degenerated sex chromosomes may represent the most dynamic genome part.

Defining management units for European captive aardvarks.

The Aardvark (Orycteropus afer) is a very unique, but relatively widespread African mammal. Although some morphological variation has been observed between forest and savannah populations and among different African regions, they are all considered as a single species. However, no modern taxonomic revision is available. All captive aardvarks in Europe are believed to stem from wild born animals from Namibia, but recently several new wild-caught aardvarks from Tanzania have been integrated into the captive population. This raises the question, whether these specimens should be interbred with the existing captive population or whether there is a risk of outbreeding depression. We studied the genetic structure of the captive populations by sequencing two mitochondrial genes (cytochrome b and 16S rRNA) to assess the degree of genetic differentiation between the two source regions. Our data suggest that the aardvarks kept in European zoos belong to the same phylogenetic (mitochondrial) lineage as the differentiation in the two studied mitochondrial markers was extremely low. A more comprehensive analysis of a larger sample with well documented origin (covering the complete geographic range) and with more sensitive genetic markers is needed to infer any final conclusions concerning the aardvark's taxonomy and identification of suitable aardvark management units.

Supernumerary Marker Chromosome Identified in Asian Elephant (Elephas maximus).

We identified a small, supernumerary marker chromosome (sSMC) in two phenotypically normal Asian elephants (Elephas maximus): a female (2n = 57,XX,+mar) and her male offspring (2n = 57,XY,+mar). sSMCs are defined as structurally abnormal chromosomes that cannot be identified by conventional banding analysis since they are usually small and often lack distinct banding patterns. Although current molecular techniques can reveal their origin, the mechanism of their formation is not yet fully understood. We determined the origin of the marker using a suite of conventional and molecular cytogenetic approaches that included (a) G- and C-banding, (b) AgNOR staining, (c) preparation of a DNA clone using laser microdissection of the marker chromosome, (d) FISH with commercially available human painting and telomeric probes, and (e) FISH with centromeric DNA derived from the centromeric regions of a marker-free Asian elephant. Moreover, we present new information on the location and number of NORs in Asian and savanna elephants. We show that the metacentric marker was composed of heterochromatin with NORs at the terminal ends, originating most likely from the heterochromatic region of chromosome 27. In this context, we discuss the possible mechanism of marker formation. We also discuss the similarities between sSMCs and B chromosomes and whether the marker chromosome presented here could evolve into a B chromosome in the future.

Patterns and determinants of the global herbivorous mycobiome.

Despite their role in host nutrition, the anaerobic gut fungal (AGF) component of the herbivorous gut microbiome remains poorly characterized. Here, to examine global patterns and determinants of AGF diversity, we generate and analyze an amplicon dataset from 661 fecal samples from 34 mammalian species, 9 families, and 6 continents. We identify 56 novel genera, greatly expanding AGF diversity beyond current estimates (31 genera and candidate genera). Community structure analysis indicates that host phylogenetic affiliation, not domestication status and biogeography, shapes the community rather than. Fungal-host associations are stronger and more specific in hindgut fermenters than in foregut fermenters. Transcriptomics-enabled phylogenomic and molecular clock analyses of 52 strains from 14 genera indicate that most genera with preferences for hindgut hosts evolved earlier (44-58 Mya) than those with preferences for foregut hosts (22-32 Mya). Our results greatly expand the documented scope of AGF diversity and provide an ecologically and evolutionary-grounded model to explain the observed patterns of AGF diversity in extant animal hosts.

Development and validation of an LC-MS/MS method for determination of B vitamins and some its derivatives in whole blood.

Obligate symbiotic bacteria associated with the insects feeding exclusively on vertebrate blood are supposed to complement B vitamins presumably lacking in their diet. Recent genomic analyses revealed considerable differences in biosynthetic capacities across different symbionts, suggesting that levels of B vitamins may vary across different vertebrate hosts. However, a rigorous determination of B vitamins content in blood of various vertebrates has not yet been approached. A reliable analytical method focused on B vitamin complex in blood can provide valuable informative background and understanding of general principles of insect symbiosis. In this work, a chromatographic separation of eight B vitamins (thiamine, riboflavin, niacin, pantothenic acid, pyridoxine, biotin, folic acid, and cyanocobalamine), four B vitamin derivatives (niacinamide, pyridoxal-5-phosphate, 4-pyridoxic acid, and tetrahydrofolic acid), and 3 stable isotope labelled internal standards was developed. Detection was carried out using dual-pressure linear ion trap mass spectrometer in FullScan MS/MS and SIM mode. Except for vitamin B9 (tetrahydrofolic acid), the instrument quantitation limits of all analytes were ranging from 0.42 to 5.0 µg/L, correlation coefficients from 0.9997 to 1.0000, and QC coefficients from 0.53 to 3.2%. Optimization of whole blood sample preparation step was focused especially on evaluation of two types of protein-precipitation agents: trichloroacetic acid and zinc sulphate in methanol. The best results were obtained for zinc sulphate in methanol, but only nine analytes were successfully validated. Accuracy of the procedure using this protein-precipitating agent was ranging from 89 to 120%, precision from 0.5 to 13%, and process efficiency from 65 to 108%. The content of B vitamins in whole blood samples from human and various vertebrates is presented as an application example of this newly developed method.

The natural cytotoxicity receptor genes in the family Felidae.

Natural killer (NK) cells belong to the innate immune system. The germline-encoded natural killer cell receptors represent activating and inhibitory receptors regulating multiple NK cell activities. The natural cytotoxicity receptors (NCRs) are activating natural cytotoxicity triggering receptors 1, 2, and 3 (NKp46, NKp44, and NKp30), encoded by the genes NCR1, NCR2, and NCR3, respectively. NCRs may be expressed in different cell types engaged in mechanisms of innate and adaptive immunity. The family Felidae, comprising the domestic cat and a wide variety of free-ranging species represents a well-suited model for biomedical and evolutionary studies. We characterized the NCR1, NCR2, and NCR3 genes in a panel of felid species. We confirmed the presence of potentially functional genes NCR1, NCR2, and NCR3 in all species. All three genes are conserved within the family and are similar to other phylogenetically related mammalian families. The NCR1 and NCR2 phylogenetic trees based on both nucleotide and protein sequences corresponded to the current zoological taxonomy, with some exceptions suggesting effects of different selection pressures in some species. Highly conserved NCR3 sequences did not allow a robust phylogenetic analysis. Most interspecific differences both at the nucleotide and protein level were found in NCR2. Within species, the most polymorphic CDS was detected in NCR1. Selection analyses indicated the effects of purifying selection on individual amino acid sites in all three genes. In stray cats, a rather high intraspecific diversity was observed.

The primate gut mycobiome-bacteriome interface is impacted by environmental and subsistence factors.

The gut microbiome of primates is known to be influenced by both host genetic background and subsistence strategy. However, these inferences have been made mainly based on adaptations in bacterial composition - the bacteriome and have commonly overlooked the fungal fraction - the mycobiome. To further understand the factors that shape the gut mycobiome of primates and mycobiome-bacteriome interactions, we sequenced 16 S rRNA and ITS2 markers in fecal samples of four different nonhuman primate species and three human groups under different subsistence patterns (n = 149). The results show that gut mycobiome composition in primates is still largely unknown but highly plastic and weakly structured by primate phylogeny, compared with the bacteriome. We find significant gut mycobiome overlap between captive apes and human populations living under industrialized subsistence contexts; this is in contrast with contemporary hunter-gatherers and agriculturalists, who share more mycobiome traits with diverse wild-ranging nonhuman primates. In addition, mycobiome-bacteriome interactions were specific to each population, revealing that individual, lifestyle and intrinsic ecological factors affect structural correspondence, number, and kind of interactions between gut bacteria and fungi in primates. Our findings indicate a dominant effect of ecological niche, environmental factors, and diet over the phylogenetic background of the host, in shaping gut mycobiome composition and mycobiome-bacteriome interactions in primates.

Karyotypic relationships in Asiatic asses (kulan and kiang) as defined using horse chromosome arm-specific and region-specific probes.

Cross-species chromosome painting has been applied to most of the species making up the numerically small family Equidae. However, comparative mapping data were still lacking in Asiatic asses kulan (Equus hemionus kulan) and kiang (E. kiang). The set of horse arm-specific probes generated by laser microdissection was hybridized onto kulan (E. hemionus kulan) and kiang (E. kiang) chromosomes in order to establish a genome-wide chromosomal correspondence between these Asiatic asses and the horse. Moreover, region-specific probes were generated to determine fusion configuration and orientation of conserved syntenic blocks. The kulan karyotype (2n = 54) was ascertained to be almost identical to the previously investigated karyotype of onager E. h. onager (2n = 56). The only difference is in fusion/fission of chromosomes homologous to horse 2q/3q, which are involved in chromosome number polymorphism in many Equidae species. E. kiang karyotype differs from the karyotype of E. hemionus by two additional fusions 8q/15 and 7/25. Chromosomes equivalent to 2q and 3q are not fused in kiang individuals with 2n = 52. Several discrepancies in centromere positions among kulan, kiang and horse chromosomes have been described. Most of the chromosome fusions in Asiatic asses are of centromere-centromere type. Comparative chromosome painting in kiang completed the efforts to establish chromosomal homologies in all representatives of the family Equidae. Application of region-specific probes allows refinement comparative maps of Asiatic asses.

Polymorphism and selection in the major histocompatibility complex DRA and DQA genes in the family Equidae.

The major histocompatibility complex genes coding for antigen binding and presenting molecules are the most polymorphic genes in the vertebrate genome. We studied the DRA and DQA gene polymorphism of the family Equidae. In addition to 11 previously reported DRA and 24 DQA alleles, six new DRA sequences and 13 new DQA alleles were identified in the genus Equus. Phylogenetic analysis of both DRA and DQA sequences provided evidence for trans-species polymorphism in the family Equidae. The phylogenetic trees differed from species relationships defined by standard taxonomy of Equidae and from trees based on mitochondrial or neutral gene sequence data. Analysis of selection showed differences between the less variable DRA and more variable DQA genes. DRA alleles were more often shared by more species. The DQA sequences analysed showed strong amongst-species positive selection; the selected amino acid positions mostly corresponded to selected positions in rodent and human DQA genes.

Epizootological study on Toxoplasma gondii in zoo animals in the Czech Republic.

Toxoplasma gondii is protozoan parasite with ability of causing disease in wide-spectrum of animals; many species of animals in captivity died of clinical toxoplasmosis. The monitoring of T. gondii antibodies in zoo animals can be an important indicator of T. gondii circulation in zoo. The aim of this study was to examine sera of animals from eight Czech zoos by latex agglutination test with statistical evaluation and detect T. gondii DNA in stray cats and rodents captured in the zoos. T. gondii antibodies were detected in 33% of 1043 zoo animals without statistical difference between birds (27%, n = 74) and mammals (33%, n = 969). In birds, the chance to be infected with T. gondii was higher in Accipitriformes (71%) compared to Pelecaniformes (6%) (p < 0.0001). In mammals, the chance to be infected with T. gondii was higher in Carnivora (63%) compared to Cetarodactyla (30%), Perissodactyla (26%), Primates (28%) and Rodentia (13%) (p < 0.0001) and higher in Felidae (70%) compared to Bovidae (28%) and Equidae (28%) (p < 0.0001). Mammals with carnivore/scavenger way of feeding were in a higher risk of T. gondii infection compared to herbivores and omnivores (p < 0.0001). T. gondii DNA was detected in tissue of one stray cat while in none of 77 rodents caught in zoo. This study is the first report on toxoplasmosis in zoos from the Czech Republic including seroepidemiology and molecular detection.

Karyotype relationships among selected deer species and cattle revealed by bovine FISH probes.

The Cervidae family comprises more than fifty species divided into three subfamilies: Capreolinae, Cervinae and Hydropotinae. A characteristic attribute for the species included in this family is the great karyotype diversity, with the chromosomal numbers ranging from 2n = 6 observed in female Muntiacus muntjak vaginalis to 2n = 70 found in Mazama gouazoubira as a result of numerous Robertsonian and tandem fusions. This work reports chromosomal homologies between cattle (Bos taurus, 2n = 60) and nine cervid species using a combination of whole chromosome and region-specific paints and BAC clones derived from cattle. We show that despite the great diversity of karyotypes in the studied species, the number of conserved chromosomal segments detected by 29 cattle whole chromosome painting probes was 35 for all Cervidae samples. The detailed analysis of the X chromosomes revealed two different morphological types within Cervidae. The first one, present in the Capreolinae is a sub/metacentric X with the structure more similar to the bovine X. The second type found in Cervini and Muntiacini is an acrocentric X which shows rearrangements in the proximal part that have not yet been identified within Ruminantia. Moreover, we characterised four repetitive sequences organized in heterochromatic blocks on sex chromosomes of the reindeer (Rangifer tarandus). We show that these repeats gave no hybridization signals to the chromosomes of the closely related moose (Alces alces) and are therefore specific to the reindeer.

New species of Torque Teno miniviruses infecting gorillas and chimpanzees.

Anelloviridae family is comprised of small, non-enveloped viruses of various genome lengths, high sequence diversity, sharing the same genome organization. Infections and co-infections by different genotypes in humans are ubiquitous. Related viruses were described in number of mammalian hosts, but very limited data are available from the closest human relatives - great apes and non-human primates. Here we report the 100% prevalence determined by semi-nested PCR from fecal samples of 16 captive primate species. Only the Mandrillus sphinx, showed the prevalence only 8%. We describe three new species of gorillas׳ and four new species of chimpanzees׳ Betatorqueviruses and their co-infections in one individual. This study is also first report and analysis of nearly full length TTMV genomes infecting gorillas. Our attempts to sequence the complete genomes of anelloviruses from host feces invariably failed. Broader usage of blood /tissue material is necessary to understand the diversity and interspecies transmission of anelloviruses.