The parvorder Rhynchophthirina with a single genus Haematomyzus is a small group of ectoparasites of unclear phylogenetic position, related to sucking and chewing lice. Previous screening based on the 16S rRNA gene indicated that Haematomyzus harbors a symbiotic bacterium whose DNA exhibits a strong shift in nucleotide composition typical of obligate mutualistic symbionts in insects. Within Phthiraptera, the smallest known genomes are found in the symbionts associated with sucking lice, which feed exclusively on mammal blood, compared to the generally larger genomes of the symbionts inhabiting chewing lice, which feed on skin derivatives. In this study, we investigate the genome characteristics of the symbiont associated with Haematomyzus elephantis. We sequenced and assembled the H. elephantis metagenome, extracted a genome draft of its symbiotic bacterium, and showed that the symbiont has a significantly reduced genome, which is with 0.39 Mbp the smallest genome among the symbionts known from Phthiraptera. Multigenic phylogenetic analysis places the symbiont into one of three clusters composed of long-branched symbionts from other insects. More specifically, it clusters together with symbionts from several other sucking lice and also with Wigglesworthia glossinidia, an obligate symbiont of tsetse flies. Consistent with the dramatic reduction of its genome, the H. elephantis symbiont lost many metabolic capacities. However, it retained functional pathways for four B vitamins, a trait typical for symbionts in blood-feeding insects. Considering genomic, metabolic, and phylogenetic characteristics, the new symbiont closely resembles those known from several sucking lice rather than chewing lice.IMPORTANCERhynchophthirina is a unique small group of permanent ectoparasites that is closely related to both sucking and chewing lice. These two groups of lice differ in their morphology, ecology, and feeding strategies. As a consequence of their different dietary sources, i.e., mammals' blood vs vertebrate skin derivatives, they also exhibit distinct patterns of symbiosis with obligate bacterial symbionts. While Rhynchophthirina shares certain traits with sucking and chewing lice, the nature of its obligate symbiotic bacterium and its metabolic role is not known. In this study, we assemble the genome of symbiotic bacterium from Haematomyzus elephantis (Rhynchophthirina), demonstrating its close similarity and phylogenetic proximity to several symbionts of sucking lice. The genome is highly reduced (representing the smallest genome among louse-associated symbionts) and exhibits a significant loss of metabolic pathways. However, similar to other sucking louse symbionts, it retains essential pathways for the synthesis of several B vitamins.

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.

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 face-masks 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.

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.

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.

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 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.

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.

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.

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.

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.

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 .

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.

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.

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.

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.

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.

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.

Although wild chimpanzees and other African great apes live in regions endemic for African sleeping sickness, very little is known about their trypanosome infections, mainly due to major difficulties in obtaining their blood samples. In present work, we established a diagnostic ITS1-based PCR assay that allows detection of the DNA of all four Trypanosoma brucei subspecies (Trypanosoma brucei brucei, Trypanosoma brucei rhodesiense, Trypanosoma brucei gambiense, and Trypanosoma brucei evansi) in feces of experimentally infected mice. Next, using this assay we revealed the presence of trypanosomes in the fecal samples of wild chimpanzees and this finding was further supported by results obtained using a set of primate tissue samples. Phylogenetic analysis of the ITS1 region showed that the majority of obtained sequences fell into the robust T. brucei group, providing strong evidence that these infections were caused by T. b. rhodesiense and/or T. b. gambiense. The optimized technique of trypanosome detection in feces will improve our knowledge about the epidemiology of trypanosomes in primates and possibly also other endangered mammals, from which blood and tissue samples cannot be obtained.Finally, we demonstrated that the mandrill serum was able to efficiently lyse T. b. brucei and T. b. rhodesiense, and to some extent T. b. gambiense, while the chimpanzee serum failed to lyse any of these subspecies.