Ancient and modern genomes unravel the evolutionary history of the rhinoceros family.

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

Epidemiology of Mycobacterium bovis infection in free-ranging rhinoceros in Kruger National Park, South Africa.

Mycobacterium bovis infection, which is a prominent cause of bovine tuberculosis, has been confirmed by mycobacterial culture in African rhinoceros species in Kruger National Park (KNP), South Africa. In this population-based study of the epidemiology of M. bovis in 437 African rhinoceros (Diceros bicornis, Ceratotherium simum), we report an estimated prevalence of 15.4% (95% CI: 10.4 to 21.0%), based on results from mycobacterial culture and an antigen-specific interferon gamma release assay from animals sampled between 2016 and 2020. A significant spatial cluster of cases was detected near the southwestern park border, although infection was widely distributed. Multivariable logistic regression models, including demographic and spatiotemporal variables, showed a significant, increasing probability of M. bovis infection in white rhinoceros based on increased numbers of African buffalo (Syncerus caffer) herds in the vicinity of the rhinoceros sampling location. Since African buffaloes are important maintenance hosts for M. bovis in KNP, spillover of infection from these hosts to white rhinoceros sharing the environment is suspected. There was also a significantly higher proportion of M. bovis infection in black rhinoceros in the early years of the study (2016­2018) than in 2019 and 2020, which coincided with periods of intense drought, although other temporal factors could be implicated. Species of rhinoceros, age, and sex were not identified as risk factors for M. bovis infection. These study findings provide a foundation for further epidemiological investigation of M. bovis, a multihost pathogen, in a complex ecosystem that includes susceptible species that are threatened and endangered.

Historic Sampling of a Vanishing Beast: Population Structure and Diversity in the Black Rhinoceros.

The black rhinoceros (Diceros bicornis L.) is a critically endangered species historically distributed across sub-Saharan Africa. Hunting and habitat disturbance have diminished both its numbers and distribution since the 19th century, but a poaching crisis in the late 20th century drove them to the brink of extinction. Genetic and genomic assessments can greatly increase our knowledge of the species and inform management strategies. However, when a species has been severely reduced, with the extirpation and artificial admixture of several populations, it is extremely challenging to obtain an accurate understanding of historic population structure and evolutionary history from extant samples. Therefore, we generated and analyzed whole genomes from 63 black rhinoceros museum specimens collected between 1775 and 1981. Results showed that the black rhinoceros could be genetically structured into six major historic populations (Central Africa, East Africa, Northwestern Africa, Northeastern Africa, Ruvuma, and Southern Africa) within which were nested four further subpopulations (Maasailand, southwestern, eastern rift, and northern rift), largely mirroring geography, with a punctuated north-south cline. However, we detected varying degrees of admixture among groups and found that several geographical barriers, most prominently the Zambezi River, drove population discontinuities. Genomic diversity was high in the middle of the range and decayed toward the periphery. This comprehensive historic portrait also allowed us to ascertain the ancestry of 20 resequenced genomes from extant populations. Lastly, using insights gained from this unique temporal data set, we suggest management strategies, some of which require urgent implementation, for the conservation of the remaining black rhinoceros diversity.

Mapping the follicle-specific regulation of extracellular vesicle-mediated microRNA transport in the southern white rhinoceros (Ceratotherium simum simum)†.

Efforts to implement effective assisted reproductive technologies (ARTs) for the conservation of the northern white rhinoceros (NWR; Ceratotherium simum cottoni) to prevent its forthcoming extinction, could be supported by research conducted on the closely related southern white rhinoceros (SWR; Ceratotherium simum simum). Within the follicle, extracellular vesicles (EVs) play a fundamental role in the bidirectional communication facilitating the crucial transport of regulatory molecules such as microRNAs (miRNAs) that control follicular growth and oocyte development. This study aimed to elucidate the dynamics of EV-miRNAs in stage-dependent follicular fluid (FF) during SWR ovarian antral follicle development. Three distinct follicular stages were identified based on diameter: Growing (G; 11-17 mm), Dominant (D; 18-29 mm), and Pre-ovulatory (P; 30-34 mm). Isolated EVs from the aspirated FF of segmented follicle stages were used to identify EV-miRNAs previously known via subsequent annotation to all equine (Equus caballus; eca), bovine (Bos taurus; bta), and human (Homo sapiens; hsa) miRNAs. A total of 417 miRNAs were detected, with 231 being mutually expressed across all three stages, including eca-miR-148a and bta-miR-451 as the top highly expressed miRNAs. Distinct expression dynamics in miRNA abundance were observed across the three follicular stages, including 31 differentially expressed miRNAs that target various pathways related to follicular growth and development, with 13 miRNAs commonly appearing amidst two different comparisons. In conclusion, this pioneering study provides a comprehensive understanding of the stage-specific expression dynamics of FF EV-miRNAs in the SWR. These findings provide insights that may lead to novel approaches in enhancing ARTs to catalyze rhinoceros conservation efforts.

Evolutionarily Young African Rhinoceros Gammaretroviruses.

High-throughput sequences were generated from DNA and cDNA from four Southern white rhinoceros (Ceratotherium simum simum) located in the Taronga Western Plain Zoo in Australia. Virome analysis identified reads that were similar to Mus caroli endogenous gammaretrovirus (McERV). Previous analysis of perissodactyl genomes did not recover gammaretroviruses. Our analysis, including the screening of the updated white rhinoceros (Ceratotherium simum) and black rhinoceros (Diceros bicornis) draft genomes identified high-copy orthologous gammaretroviral ERVs. Screening of Asian rhinoceros, extinct rhinoceros, domestic horse, and tapir genomes did not identify related gammaretroviral sequences in these species. The newly identified proviral sequences were designated SimumERV and DicerosERV for the white and black rhinoceros retroviruses, respectively. Two long terminal repeat (LTR) variants (LTR-A and LTR-B) were identified in the black rhinoceros, with different copy numbers associated with each (n = 101 and 373, respectively). Only the LTR-A lineage (n = 467) was found in the white rhinoceros. The African and Asian rhinoceros lineages diverged approximately 16 million years ago. Divergence age estimation of the identified proviruses suggests that the exogenous retroviral ancestor of the African rhinoceros ERVs colonized their genomes within the last 8 million years, a result consistent with the absence of these gammaretroviruses from Asian rhinoceros and other perissodactyls. The black rhinoceros germ line was colonized by two lineages of closely related retroviruses and white rhinoceros by one. Phylogenetic analysis indicates a close evolutionary relationship with ERVs of rodents including sympatric African rats, suggesting a possible African origin of the identified rhinoceros gammaretroviruses. IMPORTANCE Rhinoceros genomes were thought to be devoid of gammaretroviruses, as has been determined for other perissodactyls (horses, tapirs, and rhinoceros). While this may be true of most rhinoceros, the African white and black rhinoceros genomes have been colonized by evolutionarily young gammaretroviruses (SimumERV and DicerosERV for the white and black rhinoceros, respectively). These high-copy endogenous retroviruses (ERVs) may have expanded in multiple waves. The closest relative of SimumERV and DicerosERV is found in rodents, including African endemic species. Restriction of the ERVs to African rhinoceros suggests an African origin for the rhinoceros gammaretroviruses.

Estimation of the forces exerted on the limb long bones of a white rhinoceros (Ceratotherium simum) using musculoskeletal modelling and simulation.

Heavy animals incur large forces on their limb bones, due to the transmission of body weight and ground reaction forces, and the contractions of the various muscles of the limbs. This is particularly true for rhinoceroses, the heaviest extant animals capable of galloping. Several studies have examined their musculoskeletal system and the forces their bones incur, but no detailed quantification has ever been attempted. Such quantification could help understand better the link between form and function in giant land animals. Here we constructed three-dimensional musculoskeletal models of the forelimb and hindlimb of Ceratotherium simum, the heaviest extant rhino species, and used static optimisation (inverse) simulations to estimate the forces applied on the bones when standing at rest, including magnitudes and directions. Overall, unsurprisingly, the most active muscles were antigravity muscles, which generate moments opposing body weight (thereby incurring the ground reaction force), and thus keep the joints extended, avoiding joint collapse via flexion. Some muscles have an antigravity action around several joints, and thus were found to be highly active, likely specialised in body weight support (ulnaris lateralis; digital flexors). The humerus was subjected to the greatest amount of forces in terms of total magnitude; forces on the humerus furthermore came from a great variety of directions. The radius was mainly subject to high-magnitude compressive joint reaction forces, but to little muscular tension, whereas the opposite pattern was observed for the ulna. The femur had a pattern similar to that of the humerus, and the tibia's pattern was intermediate, being subject to great compression in its caudal side but to great tension in its cranial side (i.e. bending). The fibula was subject to by far the lowest force magnitude. Overall, the forces estimated were consistent with the documented morphofunctional adaptations of C. simum's long bones, which have larger insertion areas for several muscles and a greater robusticity overall than those of lighter rhinos, likely reflecting the intense forces we estimated here. Our estimates of muscle and bone (joint) loading regimes for this giant tetrapod improve the understanding of the links between form and function in supportive tissues and could be extended to other aspects of bone morphology, such as microanatomy.

Kaistella rhinocerotis sp. nov., isolated from the faeces of rhinoceros and reclassification of Chryseobacterium faecale as Kaistella faecalis comb. nov.

Strain Ran72T, a novel Gram-stain-negative, obligately aerobic, non-motile, and rod-shaped bacterium, was isolated from the faeces of the rhinoceros species Ceratotherium simum. The novel bacterial strain grew optimally in Reasoner's 2A medium under the following conditions: 0 % (w/v) NaCl, pH 7.5, and 30 °C. Based on phylogenetic analysis using 16S rRNA gene sequencing, strain Ran72T was found to be most closely related to Chryseobacterium faecale F4T (98.4 %), Kaistella soli DKR-2T (98.0 %), and Kaistella haifensis H38T (97.4 %). A comprehensive genome-level comparison between strain Ran72T with C. faecale F4T, K. soli DKR-2T, and K. haifensis H38T revealed average nucleotide identity, digital DNA-DNA hybridization, and average amino acid identity values of ≤74.9, ≤19.3, and ≤78.7 %, respectively. The major fatty acids were anteiso-C15 : 0 (22.3 %), with MK-6 being the predominant respiratory quinone. The major polar lipids of strain Ran72T were phosphatidylethanolamine, four unidentified aminolipids, and two unidentified lipids. Based on our chemotaxonomic, genotypic, and phenotype characterizations, strain Ran72T was identified as representing a novel species in the genus Kaistella, for which the name Kaistella rhinocerotis sp. nov. is proposed, with the type strain Ran72T (=KACC 23136T=JCM 36038T). Based on the outcomes of our phylogenomic study, Chryseobacterium faecale should be reclassified under the genus Kaistella as Kaistella faecalis comb. nov.

Bacteroides rhinocerotis sp. nov., isolated from the fresh feces of rhinoceros in Beijing Zoo.

A Gram-negative strain, anaerobic, non-motile, non-spore-forming, rod-shaped bacterial strain named as NGMCC 1.200684 T was isolated from the fresh feces of rhinoceros in Beijing Zoo. Based on 16S rRNA gene sequences, phylogenetic analysis indicated that strain NGMCC 1.200684 T belonged to the genus Bacteroides and was most strongly related to the type strain of Bacteroides uniformis ATCC 8492 T (96.88%). The G + C content of the genomic DNA was determined to be 46.62%. Between strains NGMCC 1.200684 T and B. uniformis ATCC 8492 T, the average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) were 93.89 and 67.60%, respectively. Strain NGMCC 1.200684 T can produce acid from fermentation of several substrates, including glucose, mannitol, lactose, saccharose, maltose, salicin, xylose, cellobiose, mannose, raffinose, sorbitol, trehalose, D-galactose, and maltotriose. The major cellular fatty acids (> 10%) were identified as anteiso-C15:0, iso-C15:0, iso-C14:0, and iso-C17:0 3-OH. The polar lipid profiles of strain NGMCC 1.200684 T were determined to contain diphosphatidyl glycerol, phosphatidylglycerol, phosphatidylethanolamine, three unknown phospholipids, and two unknown amino-phospholipids. Based on phenotypic, phylogenetic, and chemotaxonomic characteristics, a novel species of the genus Bacteroides, Bacteroides rhinocerotis sp. nov. is proposed. The type strain is NGMCC 1.200684 T (= CGMCC 1.18013 T = JCM 35702 T).

Synthesis and Characterization of Curcumin-Chitosan Loaded Gold Nanoparticles by Oryctes rhinoceros' Chitin for Cosmeceutical Application.

A breakthrough in cosmeceuticals by utilizing insects as major ingredients in cosmetic products is gaining popularity. Therefore, the interest in rare sources of ingredients, for instance, from the Oryctes rhinoceros beetle, can bring huge benefits in terms of turning pests into wealth. In this study, curcumin was chosen as the active ingredient loaded into chitosan-gold nanoparticles (CCG-NP). Curcumin is unstable and has poor absorption, a high rate of metabolism, and high sensitivity to light. These are all factors that contribute to the low bioavailability of any substance to reach the target cells. Therefore, chitosan extracted from O. rhinoceros could be used as a drug carrier to overcome these limitations. In order to overcome these limitations, CCG-NPs were synthesized and characterized. Chitosan was isolated from O. rhinoceros and CCG-NPs were successfully synthesized at 70 °C for 60 min under optimal conditions of a reactant ratio of 2:0.5 (0.5 mM HAuCl4: 0.1% curcumin). Characterizations of CCG-NP involved FTIR analysis, zeta potential, morphological properties determination by FE-SEM, particle size analysis, crystallinity study by XRD, and elemental analysis by EDX. The shape of the CCG-NP was round, its size was 128.27 d.nm, and the value of the zeta potential was 20.2 ± 3.81 mV. The IC50 value for cell viability is 58%, indicating a mild toxicity trait. To conclude, CCG-NP is a stable, spherical, nano-sized, non-toxic, and homogeneous solution.

Postoccupancy evaluation of staff, visitors, and three species of animals in a zoo setting.

Postoccupancy evaluation (POE) was used to assess newly constructed zoo exhibits from the perspective of three user groups: zoo staff, zoo visitors, and the animals. Staff represents a generally understudied user group in zoo-based POEs. We asked staff to rate the animals' space, the visitors' space, and the staff's space at previous and new exhibits. We also compared zoo visitors' ratings of the animals' behavior and environments, overall exhibit impressions, and the time visitors spent viewing previous and new exhibits. Lastly, we compared activity and space use of a Komodo dragon (Varanus komodoensis), two red pandas (Ailurus fulgens), and one rhinoceros (Rhinoceros unicornis) in their previous and new exhibits. Staff rated animal, visitor, and staff areas higher at the new exhibits compared to the previous exhibits. Visitors also rated several factors higher and spent more time at the new exhibits. The most naturalistic exhibit received the most favorable ratings in all categories and animal activity increased visitor stay time. We found that red pandas were less active in their new exhibit, and the Komodo dragon and rhino showed no difference in activity. The red pandas and the Komodo dragon used more available space in their new exhibits; however, we recommend using Electivity index to examine resource preferences for these species, whose enclosure use has been less studied than large mammals. We emphasize the importance of including staff in zoo-based POE, make other recommendations for future POE studies, and discuss various factors that could have influenced our results.

A high-quality de novo genome assembly based on nanopore sequencing of a wild-caught coconut rhinoceros beetle (Oryctes rhinoceros).

BACKGROUND: An optimal starting point for relating genome function to organismal biology is a high-quality nuclear genome assembly, and long-read sequencing is revolutionizing the production of this genomic resource in insects. Despite this, nuclear genome assemblies have been under-represented for agricultural insect pests, particularly from the order Coleoptera. Here we present a de novo genome assembly and structural annotation for the coconut rhinoceros beetle, Oryctes rhinoceros (Coleoptera: Scarabaeidae), based on Oxford Nanopore Technologies (ONT) long-read data generated from a wild-caught female, as well as the assembly process that also led to the recovery of the complete circular genome assemblies of the beetle's mitochondrial genome and that of the biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). As an invasive pest of palm trees, O. rhinoceros is undergoing an expansion in its range across the Pacific Islands, requiring new approaches to management that may include strategies facilitated by genome assembly and annotation. RESULTS: High-quality DNA isolated from an adult female was used to create four ONT libraries that were sequenced using four MinION flow cells, producing a total of 27.2 Gb of high-quality long-read sequences. We employed an iterative assembly process and polishing with one lane of high-accuracy Illumina reads, obtaining a final size of the assembly of 377.36 Mb that had high contiguity (fragment N50 length = 12 Mb) and accuracy, as evidenced by the exceptionally high completeness of the benchmarked set of conserved single-copy orthologous genes (BUSCO completeness = 99.1%). These quality metrics place our assembly ahead of the published Coleopteran genomes, including that of an insect model, the red flour beetle (Tribolium castaneum). The structural annotation of the nuclear genome assembly contained a highly-accurate set of 16,371 protein-coding genes, with only 2.8% missing BUSCOs, and the expected number of non-coding RNAs. The number and structure of paralogous genes in a gene family like Sigma GST is lower than in another scarab beetle (Onthophagus taurus), but higher than in the red flour beetle (Tribolium castaneum), which suggests expansion of this GST class in Scarabaeidae. The quality of our gene models was also confirmed with the correct placement of O. rhinoceros among other members of the rhinoceros beetles (subfamily Dynastinae) in a phylogeny based on the sequences of 95 protein-coding genes in 373 beetle species from all major lineages of Coleoptera. Finally, we provide a list of 30 candidate dsRNA targets whose orthologs have been experimentally validated as highly effective targets for RNAi-based control of several beetles. CONCLUSIONS: The genomic resources produced in this study form a foundation for further functional genetic research and management programs that may inform the control and surveillance of O. rhinoceros populations, and we demonstrate the efficacy of de novo genome assembly using long-read ONT data from a single field-caught insect.

Joining forces toward proactive elephant and rhinoceros conservation.

Proactive approaches that anticipate the long-term effects of current and future conservation threats could increase the effectiveness and efficiency of biodiversity conservation. However, such approaches can be obstructed by a lack of knowledge of habitat requirements for wildlife. To aggregate and assess the suitability of current information available on habitat requirements needed for proactive conservation, we conducted a systematic review of the literature on elephant and rhinoceros habitat requirements and synthesized data by combining a vote counting assessment with bibliometric and term maps. We contextualized these numeric and terminological results with a narrative review. We mapped current methods, results, terminology, and collaborations of 693 studies. Quantitative evidence for factors that influence the suitability of an area for elephants and rhinoceros was biased toward African savanna elephants and ecological variables. Less than one third of holistic approaches considered equal amounts of ecological and anthropogenic variables in their assessments. There was a general lack of quantitative evidence for direct proxies of anthropogenic variables that were expected to play an important role based on qualitative evidence and policy documents. However, there was evidence for a segregation in conceptual frameworks among countries and species and between science versus policy literature. There was also evidence of unused potential for collaborations among southern hemisphere researchers. Our results indicated that the success of proactive conservation interventions can be increased if ecological and anthropogenic dimensions are integrated into holistic habitat assessments and holistic carrying capacities and quantitative evidence for anthropogenic variables is improved. To avoid wasting limited resources, it is necessary to form inclusive collaborations within and across networks of researchers studying different species across regional and continental borders and in the science-policy realm.

Colaboración de Fuerzas hacia la Conservación Proactiva de Elefantes y Rinocerontes Resumen Los enfoques proactivos que anticipan los efectos a largo plazo de las amenazas a la conservación actuales y futuras podrían incrementar la efectividad y la eficiencia de la conservación de la biodiversidad. Sin embargo, dichos enfoques pueden ser frenados por una falta de conocimiento de los requerimientos de hábitat para la vida silvestre. Para sumar y evaluar la idoneidad de la información disponible sobre los requerimientos del hábitat necesitados para la conservación proactiva realizamos una revisión sistemática de la literatura sobre los requerimientos de hábitat de los elefantes y rinocerontes y sintetizamos los datos combinándolos en una evaluación de conteo de votos con mapas bibliométricos y de términos. Contextualizamos estos resultados numéricos y terminológicos con una revisión narrativa. Mapeamos los métodos, resultados, terminologías y colaboraciones actuales de 693 estudios. La evidencia cuantitativa de los factores que influyen sobre la idoneidad de un área para los elefantes y los rinocerontes estuvo sesgada hacia los elefantes de la sabana africana y las variables ecológicas. Menos de un tercio de las estrategias holísticas consideró cantidades iguales de variables ecológicas y antropogénicas en sus evaluaciones. Hubo una carencia generalizada de evidencias cuantitativas para los indicadores directos de las variables antropogénicas que se esperaba tendrían un papel importante con base en la evidencia cualitativa y los documentos de las políticas. Sin embargo, hubo evidencias de una segregación en el marco conceptual entre los países y las especies y entre la ciencia versus la literatura política. También hubo evidencias de un potencial sin explotar para las colaboraciones entre los investigadores del hemisferio sur. Nuestros resultados indicaron que el éxito de las intervenciones de conservación proactiva puede incrementarse si las dimensiones ecológicas y antropogénicas se integran a las evaluaciones holísticas del hábitat y si se mejoran las capacidades de carga holísticas y las evidencias cuantitativas de las variables antropogénicas. Para evitar gastar los recursos limitados, es necesario formar colaboraciones inclusivas dentro y a lo largo de las redes de investigadores que están estudiando a diferentes especies en las fronteras regionales y continentales y dentro del ámbito de la ciencia política.

Non-invasive assessment of fecal glucocorticoid, progesterone, and androgen metabolites and microbiome in free-ranging southern white rhinoceros (Ceratotherium simum simum) in South Africa.

Increased poaching in northern South Africa has necessitated relocation of large numbers of southern white rhinoceros (Ceratotherium simum simum) to the Eastern Cape Province. The climate and grassland ecology of this province differ from that of northern South Africa which may impact the health of this species. This assessment of fecal steroid levels and microbiome in 10 free-ranging southern white rhinoceros in the Eastern Cape will provide insights into white rhinoceros physiology in this biome. Fecal steroid metabolites were analyzed using enzyme immunoassay (EIA) and ultra-performance convergence chromatography tandem mass spectrometry (UPC2-MS/MS). Fecal microbial composition was assessed via next generation sequencing. EIAs with antibodies raised against progesterone (P4; mouse monoclonal - CL425 clone), testosterone (T; rabbit polyclonal), corticosterone (B; sheep polyclonal) were utilized. Pregnant females had large quantities of fecal progesterone metabolites (FPMs) detected by CL425 EIA. Pregnant females also had native P4 and 11α-hydroxydihydroprogesterone (11αOHDHP4; 4-pregnen-11α-ol-3,20-dione) detected by UPC2-MS/MS but these concentrations were 1000-fold less than the concentrations of FPMs detected by the CL425 EIA. By contrast, non-pregnant females had FPM concentrations detected by CL425 EIA which were similar to native P4 and 11αOHDHP4 concentrations detected by UPC2-MS/MS. Mean fecal androgen metabolite (FAM) concentrations detected by the T EIA were similar between males and females. 11-ketoandrostenedione (11KA4) detected by UPC2-MS/MS was higher in females than males. However, there was no difference between males and females in the concentration of fecal glucocorticoid metabolites (FGMs) detected by the B EIA. Bacteroidia, followed by Clostridia, was the most abundant classes of fecal microbes. The unfiltered microbiome of females was more diverse than that of males. The core fecal microbiome of young rhinoceros had a higher observed species richness (Shannon diversity index, and Simpson diversity index) than that of old rhinoceros. In the alpha male, immobilization was associated with an increase in FGMs detected by 11-deoxycortisol (S) detected by UPC2-MS/MS coupled with decreased abundance of Spirochaetia. We detected substantially different FAM and FPM concentrations from those previously reported for both captive and wild white rhinoceros. Comparison of our UPC2-MS/MS and EIA results underscores the fact that most EIAs are highly cross reactive for many steroid metabolites. Our data also demonstrates a distinct effect of stress not only on FGMs but also on the fecal microbiome. This is the first non-invasive assessment of fecal steroid metabolites by UPC2-MS/MS and the fecal microbiome in wild white rhinoceros.

Identification and Characterization of Genes for the Allo A Lectins in Japanese Rhinoceros Beetle (Trypoxylus dichotomus [Allomyrina dichotoma]).

The Japanese rhinoceros beetle (Trypoxylus dichotomus [Allomyrina dichotoma]) produces the lectins allo A-I and allo A-II, which have strong N-acetyllactosamine (Galß1-4GlcNAc)-binding activity. It has been suggested that the two lectins are formed from three subunits (α, ß, and γ), with allo A-I comprising α and γ subunits and allo A-II comprising ß and γ subunits. Here, we determined the cDNA sequences of these subunits using both conventional polymerase chain reaction (PCR)-cloning-sequencing and transcriptome-sequencing analyses. For the α and ß subunits, one gene (locus) for each was predicted, whereas for the γ subunit, two types of cDNA sequences were obtained, which we named γ1 and γ2. These two types probably have distinct loci. Average nucleotide sequence identities among the subunits ranged from 87.6% (between α and γ1) to 92.6% (between γ1 and γ2), suggesting that they form a gene family. Although no homology was found between the sequences of allo A and other known lectin proteins in a protein database search, some unknown proteins containing the DUF3421 domain were identified. Those DUF3421 domain-encoding proteins are upregulated in the insect larval midgut. Thus, we infer that allo A genes also play an important role in larvae and that their lectin activity may have been obtained collaterally.

Phenylbutazone pharmacokinetics in southern white rhinoceros (Ceratotherium simum simum) after oral administration.

Southern white rhinoceros (Ceratotherium simum simum) frequently develop painful conditions, such as traumatic injuries or osteoarthritis, necessitating the administration of pain-relieving medications. One of the preferred treatments is the nonsteroidal anti-inflammatory drug phenylbutazone because of the availability of oral formulations and the familiarity of its use in horses. For the main study, a single oral dose of phenylbutazone at 2 mg/kg was administered to healthy adult rhinoceros (n = 33) housed at six North American zoological institutions. Each rhinoceros had up to four blood samples collected under voluntary behavioural restraint at up to four predetermined time points (0, 1, 1.5, 2, 3, 4, 6, 8, 10, 24, 30 and 48 h). Drug analysis was performed by high-performance liquid chromatography. The population pharmacokinetic parameters were calculated with nonlinear mixed-effects modelling, and analysis showed a peak concentration (CMAX ) of 3.8 µg/ml at 1.8 h and an elimination half-life of 9 h. The concentrations achieved were similar to what has been reported for horses and were within the half maximal effective concentration for horses for at least 10 h. A multi-dose trial in five rhinoceros receiving 2 mg/kg orally once daily for five days found mild accumulation at a predicted factor of 1.2. This study represents the first pharmacokinetic data of phenylbutazone in any rhinoceros species.

Effects of dehorning on population productivity in four Namibia sub-populations of black rhinoceros (Diceros bicornis bicornis).

The black rhinoceros (rhino) (Diceros bicornis) is critically endangered, with poaching being one of several threats to the species' survival. Many reserves across several countries, including Namibia, South Africa, and Zimbabwe, now dehorn their rhinos in an attempt to reduce poaching. Historical data collected by the Namibian Ministry of Environment, Forestry, and Tourism and Save the Rhino Trust were used to investigate whether dehorning has an effect on age at first reproduction (AFR), inter-calving interval (ICI), birth sex ratios, calf survival, cause of death, and lifespan. These were assessed in four Namibian sub-populations (hereafter referred to as 'populations') of black rhino (denoted A, B, C, and D) which have undergone varying levels of dehorning. No significant difference was found in any of the variables between dehorned and horned individuals. Population was a significant predictor of AFR (LRT = 7.433, p = 0.024) and ICI (LRT = 13.281, p = 0.004), although pairwise comparisons only found populations A and B to be significantly different (AFR: z = -2.736, p = 0.017, ICI: z = -3.705, p = 0.001). Additionally, a significantly higher number of males than females were born in population D (p = 0.021, CI = 0.544, 0.960). The main cause of death across all individuals was poaching, although there was no significant difference in the proportion of rhinos poached between dehorned and horned individuals (X 2 = 0.638, p = 0.424, n = 265). No evidence was found to suggest that dehorning has any effect on AFR, ICI, birth sex ratios, calf survival, or lifespan in the black rhino, which is reassuring in the continued use of dehorning as an anti-poaching technique in this species.

Interspecific Gene Flow and the Evolution of Specialization in Black and White Rhinoceros.

Africa's black (Diceros bicornis) and white (Ceratotherium simum) rhinoceros are closely related sister-taxa that evolved highly divergent obligate browsing and grazing feeding strategies. Although their precursor species Diceros praecox and Ceratotherium mauritanicum appear in the fossil record ∼5.2 Ma, by 4 Ma both were still mixed feeders, and were even spatiotemporally sympatric at several Pliocene sites in what is today Africa's Rift Valley. Here, we ask whether or not D. praecox and C. mauritanicum were reproductively isolated when they came into Pliocene secondary contact. We sequenced and de novo assembled the first annotated black rhinoceros reference genome and compared it with available genomes of other black and white rhinoceros. We show that ancestral gene flow between D. praecox and C. mauritanicum ceased sometime between 3.3 and 4.1 Ma, despite conventional methods for the detection of gene flow from whole genome data returning false positive signatures of recent interspecific migration due to incomplete lineage sorting. We propose that ongoing Pliocene genetic exchange, for up to 2 My after initial divergence, could have potentially hindered the development of obligate feeding strategies until both species were fully reproductively isolated, but that the more severe and shifting paleoclimate of the early Pleistocene was likely the ultimate driver of ecological specialization in African rhinoceros.

Transcription Profile and Genomic Variations of Oryctes Rhinoceros Nudivirus in Coconut Rhinoceros Beetles.

Oryctes rhinoceros nudivirus (OrNV) is a double-stranded DNA (dsDNA) virus which has been used as a biocontrol agent to suppress the coconut rhinoceros beetle (Oryctes rhinoceros) in Southeast Asia and the Pacific Islands. A new wave of O. rhinoceros incursions in Oceania is thought to be related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles. In this study, chronically infected beetles were collected from Philippines, Fiji, Papua New Guinea (PNG), and the Solomon Islands (SI). RNA sequencing (RNA-seq) was performed to investigate the global viral gene expression profiles and for comparative genomic analysis of structural variations. Maximum likelihood phylogenic analysis indicated that OrNV strains from the SI and Philippines are closely related, while OrNV strains from PNG and Fiji formed a distinct adjacent clade. We detected several polymorphic sites with a frequency higher than 35% in 892 positions of the viral genome. Nonsynonymous mutations were detected in several hypothetical proteins and 15 nudivirus core genes, such as gp034, lef-8, lef-4, and vp91 We found limited evidence of variation in viral gene expression among geographic populations. Only a few genes, such as gp01, gp022, and gp107, were differentially expressed among different strains. Additionally, small RNA sequencing from the SI population suggested that OrNV is targeted by the host RNA interference (RNAi) response with abundant 21-nucleotide small RNAs. Some of these genomic changes are specific to the geographic population and could be related to particular phenotypic characteristics of the strain, such as viral pathogenicity or transmissibility, and this requires further investigation.IMPORTANCE Oryctes rhinoceros nudivirus has been an effective biocontrol agent against the coconut rhinoceros beetle in Southeast Asia and the Pacific Islands for decades. The recent outbreak of these beetles in many South Pacific islands has had a significant impact on livelihoods in the region. It has been suggested that the resurgence and spread of the pest are related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles. We examined viral genomic and transcriptional variations in chronically infected beetles from different geographical populations. A high number of polymorphic sites among several geographical strains of OrNV were identified, but potentially only a few of these variations in the genome are involved in functional changes and can potentially alter the typical function. These findings provide valuable resources for future studies to improve our understanding of the OrNV genetic variations in different geographic regions and their potential link to virus pathogenicity.

Historical population declines prompted significant genomic erosion in the northern and southern white rhinoceros (Ceratotherium simum).

Large vertebrates are extremely sensitive to anthropogenic pressure, and their populations are declining fast. The white rhinoceros (Ceratotherium simum) is a paradigmatic case: this African megaherbivore has suffered a remarkable decline in the last 150 years due to human activities. Its subspecies, the northern (NWR) and the southern white rhinoceros (SWR), however, underwent opposite fates: the NWR vanished quickly, while the SWR recovered after the severe decline. Such demographic events are predicted to have an erosive effect at the genomic level, linked to the extirpation of diversity, and increased genetic drift and inbreeding. However, there is little empirical data available to directly reconstruct the subtleties of such processes in light of distinct demographic histories. Therefore, we generated a whole-genome, temporal data set consisting of 52 resequenced white rhinoceros genomes, representing both subspecies at two time windows: before and during/after the bottleneck. Our data reveal previously unknown population structure within both subspecies, as well as quantifiable genomic erosion. Genome-wide heterozygosity decreased significantly by 10% in the NWR and 36% in the SWR, and inbreeding coefficients rose significantly by 11% and 39%, respectively. Despite the remarkable loss of genomic diversity and recent inbreeding it suffered, the only surviving subspecies, the SWR, does not show a significant accumulation of genetic load compared to its historical counterpart. Our data provide empirical support for predictions about the genomic consequences of shrinking populations, and our findings have the potential to inform the conservation efforts of the remaining white rhinoceroses.

Can Biological Control Overcome the Threat From Newly Invasive Coconut Rhinoceros Beetle Populations (Coleoptera: Scarabaeidae)? A Review.

The coconut rhinoceros beetle (CRB: Oryctes rhinoceros Linnaeus) is one of the most damaging pests to coconut and oil palms in Asia and the Pacific Islands. Adults bore into the crown and damage developing fronds, which affects tree development and yield. The insect is native to South and Southeast Asia and was inadvertently introduced into the Pacific in 1909. It has since spread to several Pacific island nations and territories, causing significant economic impact on these important coconut and palm-growing regions. In the 1950s and 1960s, an international biological control effort was initiated to search for and release natural enemy species. Release of the Oryctes rhinoceros nudivirus Huger (OrNV) and the species complex of Metarhizium Sorokin (Hypocreales: Clavicipitaceae) was successful in controlling CRB in its invaded range. Recently a new biotype of the beetle, known as CRB-G, has spread into the Pacific Islands causing unprecedented levels of damage due to the failure of previously successful biological control agents (BCAs) to suppress this biotype. The re-emergence of CRB as a serious pest warrants a rigorous re-evaluation of potential BCAs and a new search for effective natural enemies if necessary. In this article, we review literature on CRB to 1) analyze past introductions of BCAs and their effectiveness; 2) identify potentially important natural enemies and their geographical origins; and 3) assess possible approaches for utilization of BCAs against the new wave of CRB invasion. Research gaps and directions deserving future attention are highlighted and a strategy for renovation of biological controls for CRB suggested.