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

Ere, a Family of Short Interspersed Elements in the Genomes of Odd-Toed Ungulates (Perissodactyla)

Short Interspersed Elements (SINEs) are eukaryotic retrotransposons transcribed by RNA polymerase III (pol III). Many mammalian SINEs (T+ SINEs) contain a polyadenylation signal (AATAAA), a pol III transcription terminator, and an A-rich tail in their 3'-end. The RNAs of such SINEs have the capacity for AAUAAA-dependent polyadenylation, which is unique to pol III-generated transcripts. The structure, evolution, and polyadenylation of the Ere SINE of ungulates (horses, rhinos, and tapirs) were investigated in this study. A bioinformatics analysis revealed the presence of up to ~4 × 105 Ere copies in representatives of all three families. These copies can be classified into two large subfamilies, EreA and EreB, the former distinguished by an additional 60 bp sequence. The 3'-end of numerous EreA and all EreB copies exhibit a 50 bp sequence designated as a terminal domain (TD). The Ere family can be further subdivided into subfamilies EreA_0TD, EreA_1TD, EreB_1TD, and EreB_2TD, depending on the presence and number of terminal domains (TDs). Only EreA_0TD copies can be assigned to T+ SINEs as they contain the AATAAA signal and the TCTTT transcription terminator. The analysis of young Ere copies identified by comparison with related perissodactyl genomes revealed that EreA_0TD and, to a much lesser extent, EreB_2TD have retained retrotranspositional activity in the recent evolution of equids and rhinoceroses. The targeted mutagenesis and transfection of HeLa cells were used to identify sequences in equine EreA_0TD that are critical for the polyadenylation of its pol III transcripts. In addition to AATAAA and the transcription terminator, two sites in the 3' half of EreA, termed the β and τ signals, were found to be essential for this process. The evolution of Ere, with a particular focus on the emergence of T+ SINEs, as well as the polyadenylation signals are discussed in comparison with other T+ SINEs.

Genetic load and viability of a future restored northern white rhino population

As biodiversity loss outpaces recovery, conservationists are increasingly turning to novel tools for preventing extinction, including cloning and in vitro gametogenesis of biobanked cells. However, restoration of populations can be hindered by low genetic diversity and deleterious genetic load. The persistence of the northern white rhino (Ceratotherium simum cottoni) now depends on the cryopreserved cells of 12 individuals. These banked genomes have higher genetic diversity than southern white rhinos (C. s. simum), a sister subspecies that successfully recovered from a severe bottleneck, but the potential impact of genetic load is unknown. We estimated how demographic history has shaped genome-wide genetic load in nine northern and 13 southern white rhinos. The bottleneck left southern white rhinos with more fixed and homozygous deleterious alleles and longer runs of homozygosity, whereas northern white rhinos retained more deleterious alleles masked in heterozygosity. To gauge the impact of genetic load on the fitness of a northern white rhino population restored from biobanked cells, we simulated recovery using fitness of southern white rhinos as a benchmark for a viable population. Unlike traditional restoration, cell-derived founders can be reintroduced in subsequent generations to boost lost genetic diversity and relieve inbreeding. In simulations with repeated reintroduction of founders into a restored population, the fitness cost of genetic load remained lower than that borne by southern white rhinos. Without reintroductions, rapid growth of the restored population (>20-30% per generation) would be needed to maintain comparable fitness. Our results suggest that inbreeding depression from genetic load is not necessarily a barrier to recovery of the northern white rhino and demonstrate how restoration from biobanked cells relieves some constraints of conventional restoration from a limited founder pool. Established conservation methods that protect healthy populations will remain paramount, but emerging technologies hold promise to bolster these tools to combat the extinction crisis.

ONT read assembly of the black rhino genome

OBJECTIVES

The black rhinoceros (Diceros bicornis) is an endangered mammal for which a captive breeding program is part of the conservation effort. Black rhinos in zoo's often suffer from chronic infections and heamochromatosis. Furthermore, breeding is hampered by low male fertility. To aid a research project studying these topics, we sequenced and assembled the genome of a captive male black rhino using ONT sequencing data only.

DATA DESCRIPTION

This work produced over 100 Gb whole genome sequencing reads from whole blood. These were assembled into a 2.47 Gb draft genome consisting of 834 contigs with an N50 of 29.53 Mb. The genome annotation was lifted over from an available genome annotation for black rhino, which resulted in the retrieval of over 99% of gene features. This new genome assembly will be a valuable resource in for conservation genetic research in this species.

Metabolomic profiling implicates mitochondrial and immune dysfunction in disease syndromes of the critically endangered black rhinoceros (Diceros bicornis)

The critically endangered black rhinoceros (Diceros bicornis; black rhino) experiences extinction threats from poaching in-situ. The ex-situ population, which serves as a genetic reservoir against impending extinction threats, experiences its own threats to survival related to several disease syndromes not typically observed among their wild counterparts. We performed an untargeted metabolomic analysis of serum from 30 ex-situ housed black rhinos (Eastern black rhino, EBR, n = 14 animals; Southern black rhino, SBR, n = 16 animals) and analyzed differences in metabolite profiles between subspecies, sex, and health status (healthy n = 13 vs. diseased n = 14). Of the 636 metabolites detected, several were differentially (fold change > 1.5; p < 0.05) expressed between EBR vs. SBR (40 metabolites), female vs. male (36 metabolites), and healthy vs. diseased (22 metabolites). Results suggest dysregulation of propanoate, amino acid metabolism, and bile acid biosynthesis in the subspecies and sex comparisons. Assessment of healthy versus diseased rhinos indicates involvement of arachidonic acid metabolism, bile acid biosynthesis, and the pentose phosphate pathway in animals exhibiting inflammatory disease syndromes. This study represents the first systematic characterization of the circulating serum metabolome in the black rhinoceros. Findings further implicate mitochondrial and immune dysfunction as key contributors for the diverse disease syndromes reported in ex-situ managed black rhinos.

A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea

The impact of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae and Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including (i) phylogenomics, (ii) the distribution of shared derived alleles and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggest that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species' gene pools.

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.

High-Throughput Sequencing Reveals a Dynamic Bacterial Linkage between the Captive White Rhinoceros and Its Environment

Soil is an essential part of the animal habitat and has a large diversity of microbiota, while the animal body was colonized by a complex bacterial community; so far, the relationship between the animal host microbial community and the soil microbial ecosystem remains largely unknown. In this study, 15 white rhinoceros from three different captive grounds were selected and the bacterial community of the gut, skin, and environment of these rhinoceros were analyzed by 16S rRNA sequencing technology. Our results showed that Firmicutes and Bacteroidota were the predominant phyla in the gut microbiome, whereas skin and environment samples share similar microbiome profiles and are dominated by the phyla of Actinobacteriota, Chloroflexi, and Proteobacteria. Although the bacterial composition of the gut differs from that of the skin and environment, the Venn diagrams showed that there were 22 phyla and 186 genera shared by all the gut, skin, and environmental microbes in white rhinoceroses. Further cooccurrence network analysis indicated a bacterial linkage based on a complex interaction was established by the bacterial communities from the three different niches. In addition, beta diversity and bacterial composition analysis showed that both the captive ground and host ages induced shifts in the microbial composition of white rhinoceroses, which suggested that the bacterial linkage between the captive white rhinoceros and its environment is dynamic. Overall, our data contribute to a better understanding of the bacterial community of the captive white rhinoceros, especially for the relationship between the environment and animal bacterial communities. IMPORTANCE The white rhinoceros is one of the world's most endangered mammals. The microbial population plays a key role in animal health and welfare; however, studies regarding the microbial communities of the white rhinoceros are relatively limited. As the white rhinoceros has a common behavior of mud baths and thus is in direct contact with the environment, a relationship between the animal microbial community and the soil microbial ecosystem appears possible, but it remains unclear. Here, we described the characteristics and interaction of bacterial communities of the white rhinoceros in three different niches, including gut, skin, and environment. We also analyzed the effects of captive ground and age on the composition of the bacterial community. Our findings highlighted the relationship among the three niches and may have important implications for the conservation and management of this threatened species.

The Contribution of Digital Sequence Information to Conservation Biology: A Southern African Perspective

Many recent contributions have made a compelling case that genetic diversity is not adequately reflected in international frameworks and policies, as well as in local governmental processes implementing such frameworks. Using digital sequence information (DSI) and other publicly available data is supported to assess genetic diversity, toward formulation of practical actions for long-term conservation of biodiversity, with the particular goal of maintaining ecological and evolutionary processes. Given the inclusion of specific goals and targets regarding DSI in the latest draft of the Global Biodiversity Framework negotiated at the 15^th Conference of the Parties (COP15) in Montreal in December 2022 and the crucial decisions on access and benefit sharing to DSI that will be taken in the coming months and future COP meetings, a southern African perspective on how and why open access to DSI is essential for the conservation of intraspecific biodiversity (genetic diversity and structure) across country borders is provided.

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.

Genomic insights into the evolutionary relationships and demographic history of kiwi

Kiwi are a unique and emblematic group of birds endemic to New Zealand. Deep-time evolutionary relationships among the five extant kiwi species have been difficult to resolve, in part due to the absence of pre-Quaternary fossils to inform speciation events. Here, we utilise single representative nuclear genomes of all five extant kiwi species (great spotted kiwi, little spotted kiwi, Okarito brown kiwi, North Island brown kiwi, and southern brown kiwi) and investigate their evolutionary histories with phylogenomic, genetic diversity, and deep-time (past million years) demographic analyses. We uncover relatively low levels of gene-tree phylogenetic discordance across the genomes, suggesting clear distinction between species. However, we also find indications of post-divergence gene flow, concordant with recent reports of interspecific hybrids. The four species for which unbiased levels of genetic diversity could be calculated, due to the availability of reference assemblies (all species except the southern brown kiwi), show relatively low levels of genetic diversity, which we suggest reflects a combination of older environmental as well as more recent anthropogenic influence. In addition, we suggest hypotheses regarding the impact of known past environmental events, such as volcanic eruptions and glacial periods, on the similarities and differences observed in the demographic histories of the five kiwi species over the past million years.

Evolutionary and Ecological Drivers of Local Adaptation and Speciation in a North American Avian Species Complex

Throughout the speciation process, genomic divergence can be differentially impacted by selective pressures, as well as gene flow and genetic drift. Disentangling the effects of these evolutionary mechanisms remains challenging, especially for non-model organisms. Accounting for complex evolutionary histories and contemporary population structure often requires sufficient sample sizes, for which the expense of full genomes remains prohibitive. Here, we demonstrate the utility of partial-genome sequence data for range-wide samples to shed light into the divergence process of two closely related ducks, the Mexican duck (Anas diazi) and mallard (A. platyrhynchos). We determine the role of selective and neutral processes during speciation of Mexican ducks by integrating evolutionary and demographic modelling with genotype-environment and genotype-phenotype association testing. First, evolutionary models and demographic analyses support the hypothesis that Mexican ducks originally diverged ~300,000 years ago in a climate refugia arising during a glacial period in in a southwestern North America, and that subsequent environmental selective pressures played a key role in divergence. Mexican ducks then showed cyclical demographic patterns that likely reflected repeated range expansions and contractions, along with bouts of gene flow with mallards during glacial cycles. Finally, we provide evidence that sexual selection acted on several phenotypic traits as a co-evolutionary process, facilitating the development of reproductive barriers that initially arose due to strong ecological selection. More broadly, this work reveals that the genomic and phenotypic patterns observed across species complexes are the result of myriad factors that contribute in dynamic ways to the evolutionary trajectories of a lineage.

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.

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.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.