Colonialism in South Africa leaves a lasting legacy of reduced genetic diversity in Cape buffalo.

The iconic Cape buffalo has experienced several documented population declines in recent history. These declines have been largely attributed to the late 19th century rinderpest pandemic. However, the effect of the rinderpest pandemic on their genetic diversity remains contentious, and other factors that have potentially affected this diversity include environmental changes during the Pleistocene, range expansions and recent human activity. Motivated by this, we present analyses of whole genome sequencing data from 59 individuals from across the Cape buffalo range to assess present-day levels of genome-wide genetic diversity and what factors have influenced these levels. We found that the Cape buffalo has high average heterozygosity overall (0.40%), with the two southernmost populations having significantly lower heterozygosity levels (0.33% and 0.29%) on par with that of the domesticated water buffalo (0.29%). Interestingly, we found that these lower levels are probably due to recent inbreeding (average fraction of runs of homozygosity 23.7% and 19.9%) rather than factors further back in time during the Pleistocene. Moreover, detailed investigations of recent demographic history show that events across the past three centuries were the main drivers of the exceptional loss of genetic diversity in the southernmost populations, coincident with the onset of colonialism in the southern extreme of the Cape buffalo range. Hence, our results add to the growing body of studies suggesting that multiple recent human-mediated impacts during the colonial period caused massive losses of large mammal abundance in southern Africa.

Gaur genome reveals expansion of sperm odorant receptors in domesticated cattle.

BACKGROUND: The gaur (Bos gaurus) is the largest extant wild bovine species, native to South and Southeast Asia, with unique traits, and is listed as vulnerable by the International Union for Conservation of Nature (IUCN). RESULTS: We report the first gaur reference genome and identify three biological pathways including lysozyme activity, proton transmembrane transporter activity, and oxygen transport with significant changes in gene copy number in gaur compared to other mammals. These may reflect adaptation to challenges related to climate and nutrition. Comparative analyses with domesticated indicine (Bos indicus) and taurine (Bos taurus) cattle revealed genomic signatures of artificial selection, including the expansion of sperm odorant receptor genes in domesticated cattle, which may have important implications for understanding selection for male fertility. CONCLUSIONS: Apart from aiding dissection of economically important traits, the gaur genome will also provide the foundation to conserve the species.

Impact of Four Ovine TMEM154 Haplotypes on Ewes during Multiyear Lentivirus Exposure.

Polypeptide variation encoded by the ovine transmembrane protein 154 gene (TMEM154) is associated with susceptibility to ovine lentivirus, the causative agent of Ovine Progressive Pneumonia (OPP) and Visna/Maedi. Our aim was to compare the four most prevalent TMEM154 haplotypes on the incidence of infection and ewe productivity during natural multiyear virus exposure. Prospective cohort studies were designed to test gene action and estimate effects of TMEM154 haplotypes encoding distinctive variant residues: K35 ("1"), I70 ("2"), ancestral ("3"), and A4del/M44 ("4"). Exposure consisted of co-mingling infected ewes at a rate greater than 30% with serological status evaluated every four months. For ewes with one or two copies of the highly susceptible haplotypes "2" and "3", the infection prevalence steadily increased to nearly 100% at 55 months. Haplotypes "2" and "3" were equally susceptible and dominant to haplotype "1". A difference was not detected (p < 0.53) in the magnitude of effect with haplotype combinations of "1" and "4". The ewe infection prevalence with "1,1"; "1,4"; and "4,4" was 10% to 40% at 55 months. The latter suggested that two copies of the K35 amino acid substitution ("1") were as effective as a homozygous TMEM154 "knockout" with the frame-shift deletion mutation ("4") in reducing infection susceptibility. When considering ewe reproductive performance, a difference was not detected when comparing haplotypes "2", and "3" to each other, or "1" and "4" to each other. Our study indicated that ewes with two copies of the severely truncated versions of TMEM154 ("4,4") had normal lamb productivity. Without complete understanding of the natural function of TMEM154 our recommendations to producers interested in using TMEM154 selection to reduce their flock's genetic predisposition to OPP are encouraged to increase the frequency of TMEM154 haplotype K35 ("1") since it encodes a full-length protein with minimal difference to the ancestral polypeptide.

A Reference Genome Assembly of Simmental Cattle, Bos taurus taurus.

Genomics research has relied principally on the establishment and curation of a reference genome for the species. However, it is increasingly recognized that a single reference genome cannot fully describe the extent of genetic variation within many widely distributed species. Pangenome representations are based on high-quality genome assemblies of multiple individuals and intended to represent the broadest possible diversity within a species. A Bovine Pangenome Consortium (BPC) has recently been established to begin assembling genomes from more than 600 recognized breeds of cattle, together with other related species to provide information on ancestral alleles and haplotypes. Previously reported de novo genome assemblies for Angus, Brahman, Hereford, and Highland breeds of cattle are part of the initial BPC effort. The present report describes a complete single haplotype assembly at chromosome-scale for a fullblood Simmental cow from an F1 bison-cattle hybrid fetus by trio binning. Simmental cattle, also known as Fleckvieh due to their red and white spots, originated in central Europe in the 1830s as a triple-purpose breed selected for draught, meat, and dairy production. There are over 50 million Simmental cattle in the world, known today for their fast growth and beef yields. This assembly (ARS_Simm1.0) is similar in length to the other bovine assemblies at 2.86 Gb, with a scaffold N50 of 102 Mb (max scaffold 156.8 Mb) and meets or exceeds the continuity of the best Bos taurus reference assemblies to date.

A Reference Genome Assembly of American Bison, Bison bison bison.

Bison are an icon of the American West and an ecologically, commercially, and culturally important species. Despite numbering in the hundreds of thousands today, conservation concerns remain for the species, including the impact on genetic diversity of a severe bottleneck around the turn of the 20th century and genetic introgression from domestic cattle. Genetic diversity and admixture are best evaluated at genome-wide scale, for which a high-quality reference is necessary. Here, we use trio binning of long reads from a bison-Simmental cattle (Bos taurus taurus) male F1 hybrid to sequence and assemble the genome of the American plains bison (Bison bison bison). The male haplotype genome is chromosome-scale, with a total length of 2.65 Gb across 775 scaffolds (839 contigs) and a scaffold N50 of 87.8 Mb. Our bison genome is ~13× more contiguous overall and ~3400× more contiguous at the contig level than the current bison reference genome. The bison genome sequence presented here (ARS-UCSC_bison1.0) will enable new research into the evolutionary history of this iconic megafauna species and provide a new tool for the management of bison populations in federal and commercial herds.

Evaluation of CD109, PCP4 and SEMA3D genes for their association with Ovine Johne's disease in Turkish sheep.

Johne's disease is a chronic, contagious, zoonotic disease that affects numerous species including livestock and sometimes humans. The disease is globally distributed in sheep populations and caused by Mycobacterium avium Subsp. paratuberculosis (MAP). A previous genome-wide association study identified single nucleotide polymorphism (SNP) markers associated with OJD serostatus in CD109, PCP4, and SEMA3D genes. Our aim was to evaluate the same markers for association with OJD seroprevalence in Turkish sheep in a retrospective matched case-control study. The serological status for OJD in 1801 sheep was determined for four native and four composite breeds from three research flocks. One hundred eleven matched case-control pairs were constructed according to breed type and age from 1750 comingled ewes reared in the same environment. A Single Nucleotide Primer Extension (SNuPE) assay was designed to genotype PCP4-Intron 1, PCP4-3'UTR, SEMA3D, CD109-intron 2 and CD109-intron 8 markers and a McNemar's test was performed on the matched pairs. An association with these five markers was not detected with the OJD serostatus in Turkish sheep (power of detection, 0.95; odds ratio >3; McNemar's p < .05). Thus, a wider search may be needed to identify any major underlying genetic risk factors for OJD in Turkish sheep.

Structure-Guided Identification of a Nonhuman Morbillivirus with Zoonotic Potential.

Morbilliviruses infect a broad range of mammalian hosts, including ruminants, carnivores, and humans. The recent eradication of rinderpest virus (RPV) and the active campaigns for eradication of the human-specific measles virus (MeV) have raised significant concerns that the remaining morbilliviruses may emerge in so-called vacated ecological niches. Seeking to assess the zoonotic potential of nonhuman morbilliviruses within human populations, we found that peste des petits ruminants virus (PPRV)-the small-ruminant morbillivirus-is restricted at the point of entry into human cells due to deficient interactions with human SLAMF1-the immune cell receptor for morbilliviruses. Using a structure-guided approach, we characterized a single amino acid change, mapping to the receptor-binding domain in the PPRV hemagglutinin (H) protein, which overcomes this restriction. The same mutation allowed escape from some cross-protective, human patient, anti-MeV antibodies, raising concerns that PPRV is a pathogen with zoonotic potential. Analysis of natural variation within human and ovine SLAMF1 also identified polymorphisms that could correlate with disease resistance. Finally, the mechanistic nature of the PPRV restriction was also investigated, identifying charge incompatibility and steric hindrance between PPRV H and human SLAMF1 proteins. Importantly, this research was performed entirely using surrogate virus entry assays, negating the requirement for in situ derivation of a human-tropic PPRV and illustrating alternative strategies for identifying gain-of-function mutations in viral pathogens.IMPORTANCE A significant proportion of viral pandemics occur following zoonotic transmission events, where animal-associated viruses jump species into human populations. In order to provide forewarnings of the emergence of these viruses, it is necessary to develop a better understanding of what determines virus host range, often at the genetic and structural levels. In this study, we demonstrated that the small-ruminant morbillivirus, a close relative of measles, is unable to use human receptors to enter cells; however, a change of a single amino acid in the virus is sufficient to overcome this restriction. This information will be important for monitoring this virus's evolution in the field. Of note, this study was undertaken in vitro, without generation of a fully infectious virus with this phenotype.

Worldwide patterns of ancestry, divergence, and admixture in domesticated cattle.

The domestication and development of cattle has considerably impacted human societies, but the histories of cattle breeds and populations have been poorly understood especially for African, Asian, and American breeds. Using genotypes from 43,043 autosomal single nucleotide polymorphism markers scored in 1,543 animals, we evaluate the population structure of 134 domesticated bovid breeds. Regardless of the analytical method or sample subset, the three major groups of Asian indicine, Eurasian taurine, and African taurine were consistently observed. Patterns of geographic dispersal resulting from co-migration with humans and exportation are recognizable in phylogenetic networks. All analytical methods reveal patterns of hybridization which occurred after divergence. Using 19 breeds, we map the cline of indicine introgression into Africa. We infer that African taurine possess a large portion of wild African auroch ancestry, causing their divergence from Eurasian taurine. We detect exportation patterns in Asia and identify a cline of Eurasian taurine/indicine hybridization in Asia. We also identify the influence of species other than Bos taurus taurus and B. t. indicus in the formation of Asian breeds. We detect the pronounced influence of Shorthorn cattle in the formation of European breeds. Iberian and Italian cattle possess introgression from African taurine. American Criollo cattle originate from Iberia, and not directly from Africa with African ancestry inherited via Iberian ancestors. Indicine introgression into American cattle occurred in the Americas, and not Europe. We argue that cattle migration, movement and trading followed by admixture have been important forces in shaping modern bovine genomic variation.

Genomic signatures of Mannheimia haemolytica that associate with the lungs of cattle with respiratory disease, an integrative conjugative element, and antibiotic resistance genes.

BACKGROUND: Mannheimia haemolytica typically resides in cattle as a commensal member of the upper respiratory tract microbiome. However, some strains can invade their lungs and cause respiratory disease and death, including those with multi-drug resistance. A nucleotide polymorphism typing system was developed for M. haemolytica from the genome sequences of 1133 North American isolates, and used to identify genetic differences between isolates from the lungs and upper respiratory tract of cattle with and without clinical signs of respiratory disease. RESULTS: A total of 26,081 nucleotide polymorphisms were characterized after quality control filtering of 48,403 putative polymorphisms. Phylogenetic analyses of nucleotide polymorphism genotypes split M. haemolytica into two major genotypes (1 and 2) that each were further divided into multiple subtypes. Multiple polymorphisms were identified with alleles that tagged genotypes 1 or 2, and their respective subtypes. Only genotype 2 M. haemolytica associated with the lungs of diseased cattle and the sequence of a particular integrative and conjugative element (ICE). Additionally, isolates belonging to one subtype of genotype 2 (2b), had the majority of antibiotic resistance genes detected in this study, which were assorted into seven combinations that ranged from 1 to 12 resistance genes. CONCLUSIONS: Typing of diverse M. haemolytica by nucleotide polymorphism genotypes successfully identified associations with diseased cattle lungs, ICE sequence, and antibiotic resistance genes. Management of cattle by their carriage of M. haemolytica could be an effective intervention strategy to reduce the prevalence of respiratory disease and supplemental needs for antibiotic treatments in North American herds.

Genetic subgroup of small ruminant lentiviruses that infects sheep homozygous for TMEM154 frameshift deletion mutation A4Δ53.

Small ruminant lentivirus (SRLV) infections of sheep are influenced by genetics on both the host and pathogen sides. Genetic variation in the ovine transmembrane 154 (TMEM154) gene associates with infection susceptibility, and distinct SRLV genetic subgroups infect sheep in association with their TMEM154 diplotypes. In this study, a novel SRLV subgroup was identified that naturally infected sheep with various TMEM154 diplotypes, including those homozygous for a rare frameshift mutation (A4 delta53), which is predicted to abolish TMEM154 protein function. Thus, these SRLVs may infect sheep that lack functional TMEM154, and may not be restricted by TMEM154 diplotypes in establishing infections.

Reduced lentivirus susceptibility in sheep with TMEM154 mutations.

Visna/Maedi, or ovine progressive pneumonia (OPP) as it is known in the United States, is an incurable slow-acting disease of sheep caused by persistent lentivirus infection. This disease affects multiple tissues, including those of the respiratory and central nervous systems. Our aim was to identify ovine genetic risk factors for lentivirus infection. Sixty-nine matched pairs of infected cases and uninfected controls were identified among 736 naturally exposed sheep older than five years of age. These pairs were used in a genome-wide association study with 50,614 markers. A single SNP was identified in the ovine transmembrane protein (TMEM154) that exceeded genome-wide significance (unadjusted p-value 3×10(-9)). Sanger sequencing of the ovine TMEM154 coding region identified six missense and two frameshift deletion mutations in the predicted signal peptide and extracellular domain. Two TMEM154 haplotypes encoding glutamate (E) at position 35 were associated with infection while a third haplotype with lysine (K) at position 35 was not. Haplotypes encoding full-length E35 isoforms were analyzed together as genetic risk factors in a multi-breed, matched case-control design, with 61 pairs of 4-year-old ewes. The odds of infection for ewes with one copy of a full-length TMEM154 E35 allele were 28 times greater than the odds for those without (p-value<0.0001, 95% CI 5-1,100). In a combined analysis of nine cohorts with 2,705 sheep from Nebraska, Idaho, and Iowa, the relative risk of infection was 2.85 times greater for sheep with a full-length TMEM154 E35 allele (p-value<0.0001, 95% CI 2.36-3.43). Although rare, some sheep were homozygous for TMEM154 deletion mutations and remained uninfected despite a lifetime of significant exposure. Together, these findings indicate that TMEM154 may play a central role in ovine lentivirus infection and removing sheep with the most susceptible genotypes may help eradicate OPP and protect flocks from reinfection.

Characterization of ovine Nectin-4, a novel peste des petits ruminants virus receptor.

Small ruminants infected with peste des petits ruminants virus exhibit lesions typical of epithelial infection and necrosis. However, the only established host receptor for this virus is the immune cell marker signaling lymphocyte activation molecule (SLAM). We have confirmed that the ovine Nectin-4 protein, when overexpressed in epithelial cells, permits efficient replication of PPRV. Furthermore, this gene was predominantly expressed in epithelial tissues and encoded by multiple haplotypes in sheep breeds from around the world.

African bushpigs exhibit porous species boundaries and appeared in Madagascar concurrently with human arrival.

Several African mammals exhibit a phylogeographic pattern where closely related taxa are split between West/Central and East/Southern Africa, but their evolutionary relationships and histories remain controversial. Bushpigs (Potamochoerus larvatus) and red river hogs (P. porcus) are recognised as separate species due to morphological distinctions, a perceived lack of interbreeding at contact, and putatively old divergence times, but historically, they were considered conspecific. Moreover, the presence of Malagasy bushpigs as the sole large terrestrial mammal shared with the African mainland raises intriguing questions about its origin and arrival in Madagascar. Analyses of 67 whole genomes revealed a genetic continuum between the two species, with putative signatures of historical gene flow, variable FST values, and a recent divergence time (<500,000 years). Thus, our study challenges key arguments for splitting Potamochoerus into two species and suggests their speciation might be incomplete. Our findings also indicate that Malagasy bushpigs diverged from southern African populations and underwent a limited bottleneck 1000-5000 years ago, concurrent with human arrival in Madagascar. These results shed light on the evolutionary history of an iconic and widespread African mammal and provide insight into the longstanding biogeographic puzzle surrounding the bushpig's presence in Madagascar.

Introgression and disruption of migration routes have shaped the genetic integrity of wildebeest populations.

The blue wildebeest (Connochaetes taurinus) is a keystone species in savanna ecosystems from southern to eastern Africa, and is well known for its spectacular migrations and locally extreme abundance. In contrast, the black wildebeest (C. gnou) is endemic to southern Africa, barely escaped extinction in the 1900s and is feared to be in danger of genetic swamping from the blue wildebeest. Despite the ecological importance of the wildebeest, there is a lack of understanding of how its unique migratory ecology has affected its gene flow, genetic structure and phylogeography. Here, we analyze whole genomes from 121 blue and 22 black wildebeest across the genus' range. We find discrete genetic structure consistent with the morphologically defined subspecies. Unexpectedly, our analyses reveal no signs of recent interspecific admixture, but rather a late Pleistocene introgression of black wildebeest into the southern blue wildebeest populations. Finally, we find that migratory blue wildebeest populations exhibit a combination of long-range panmixia, higher genetic diversity and lower inbreeding levels compared to neighboring populations whose migration has recently been disrupted. These findings provide crucial insights into the evolutionary history of the wildebeest, and tangible genetic evidence for the negative effects of anthropogenic activities on highly migratory ungulates.

Small ruminant lentivirus genetic subgroups associate with sheep TMEM154 genotypes.

Small ruminant lentiviruses (SRLVs) are prevalent in North American sheep and a major cause of production losses for the U.S. sheep industry. Sheep susceptibility to SRLV infection is influenced by genetic variation within the ovine transmembrane 154 gene (TMEM154). Animals with either of two distinct TMEM154 haplotypes that both encode glutamate at position 35 of the protein (E35) are at greater risk of SRLV infection than those homozygous with a lysine (K35) haplotype. Prior to this study, it was unknown if TMEM154 associations with infection are influenced by SRLV genetic subgroups. Accordingly, our goals were to characterize SRLVs naturally infecting sheep from a diverse U.S. Midwestern flock and test them for associations with TMEM154 E35K genotypes. Two regions of the SRLV genome were targeted for proviral amplification, cloning, sequence analysis, and association testing with TMEM154 E35K genotypes: gag and the transmembrane region of env. Independent analyses of gag and env sequences showed that they clustered in two subgroups (1 and 2), they were distinct from SRLV subtypes originating from Europe, and that subgroup 1 associated with hemizygous and homozygous TMEM154 K35 genotypes and subgroup 2 with hemi- and homozygous E35 genotypes (gag p<0.001, env p=0.01). These results indicate that SRLVs in the U.S. have adapted to infect sheep with specific TMEM154 E35K genotypes. Consequently, both host and SRLV genotypes affect the relative risk of SRLV infection in sheep.

First gene-edited calf with reduced susceptibility to a major viral pathogen.

Bovine viral diarrhea virus (BVDV) is one of the most important viruses affecting the health and well-being of bovine species throughout the world. Here, we used CRISPR-mediated homology-directed repair and somatic cell nuclear transfer to produce a live calf with a six amino acid substitution in the BVDV binding domain of bovine CD46. The result was a gene-edited calf with dramatically reduced susceptibility to infection as measured by reduced clinical signs and the lack of viral infection in white blood cells. The edited calf has no off-target edits and appears normal and healthy at 20 months of age without obvious adverse effects from the on-target edit. This precision bred, proof-of-concept animal provides the first evidence that intentional genome alterations in the CD46 gene may reduce the burden of BVDV-associated diseases in cattle and is consistent with our stepwise, in vitro and ex vivo experiments with cell lines and matched fetal clones.

Resolving the evolution of extant and extinct ruminants with high-throughput phylogenomics.

The Pecorans (higher ruminants) are believed to have rapidly speciated in the Mid-Eocene, resulting in five distinct extant families: Antilocapridae, Giraffidae, Moschidae, Cervidae, and Bovidae. Due to the rapid radiation, the Pecoran phylogeny has proven difficult to resolve, and 11 of the 15 possible rooted phylogenies describing ancestral relationships among the Antilocapridae, Giraffidae, Cervidae, and Bovidae have each been argued as representations of the true phylogeny. Here we demonstrate that a genome-wide single nucleotide polymorphism (SNP) genotyping platform designed for one species can be used to genotype ancient DNA from an extinct species and DNA from species diverged up to 29 million years ago and that the produced genotypes can be used to resolve the phylogeny for this rapidly radiated infraorder. We used a high-throughput assay with 54,693 SNP loci developed for Bos taurus taurus to rapidly genotype 678 individuals representing 61 Pecoran species. We produced a highly resolved phylogeny for this diverse group based upon 40,843 genome-wide SNP, which is five times as many informative characters as have previously been analyzed. We also establish a method to amplify and screen genomic information from extinct species, and place Bison priscus within the Bovidae. The quality of genotype calls and the placement of samples within a well-supported phylogeny may provide an important test for validating the fidelity and integrity of ancient samples. Finally, we constructed a phylogenomic network to accurately describe the relationships between 48 cattle breeds and facilitate inferences concerning the history of domestication and breed formation.

Structural variant-based pangenome construction has low sensitivity to variability of haplotype-resolved bovine assemblies.

Advantages of pangenomes over linear reference assemblies for genome research have recently been established. However, potential effects of sequence platform and assembly approach, or of combining assemblies created by different approaches, on pangenome construction have not been investigated. Here we generate haplotype-resolved assemblies from the offspring of three bovine trios representing increasing levels of heterozygosity that each demonstrate a substantial improvement in contiguity, completeness, and accuracy over the current Bos taurus reference genome. Diploid coverage as low as 20x for HiFi or 60x for ONT is sufficient to produce two haplotype-resolved assemblies meeting standards set by the Vertebrate Genomes Project. Structural variant-based pangenomes created from the haplotype-resolved assemblies demonstrate significant consensus regardless of sequence platform, assembler algorithm, or coverage. Inspecting pangenome topologies identifies 90 thousand structural variants including 931 overlapping with coding sequences; this approach reveals variants affecting QRICH2, PRDM9, HSPA1A, TAS2R46, and GC that have potential to affect phenotype.

An improved ovine reference genome assembly to facilitate in-depth functional annotation of the sheep genome.

BACKGROUND: The domestic sheep (Ovis aries) is an important agricultural species raised for meat, wool, and milk across the world. A high-quality reference genome for this species enhances the ability to discover genetic mechanisms influencing biological traits. Furthermore, a high-quality reference genome allows for precise functional annotation of gene regulatory elements. The rapid advances in genome assembly algorithms and emergence of sequencing technologies with increasingly long reads provide the opportunity for an improved de novo assembly of the sheep reference genome. FINDINGS: Short-read Illumina (55× coverage), long-read Pacific Biosciences (75× coverage), and Hi-C data from this ewe retrieved from public databases were combined with an additional 50× coverage of Oxford Nanopore data and assembled with canu v1.9. The assembled contigs were scaffolded using Hi-C data with Salsa v2.2, gaps filled with PBsuitev15.8.24, and polished with Nanopolish v0.12.5. After duplicate contig removal with PurgeDups v1.0.1, chromosomes were oriented and polished with 2 rounds of a pipeline that consisted of freebayes v1.3.1 to call variants, Merfin to validate them, and BCFtools to generate the consensus fasta. The ARS-UI_Ramb_v2.0 assembly is 2.63 Gb in length and has improved continuity (contig NG50 of 43.18 Mb), with a 19- and 38-fold decrease in the number of scaffolds compared with Oar_rambouillet_v1.0 and Oar_v4.0. ARS-UI_Ramb_v2.0 has greater per-base accuracy and fewer insertions and deletions identified from mapped RNA sequence than previous assemblies. CONCLUSIONS: The ARS-UI_Ramb_v2.0 assembly is a substantial improvement in contiguity that will optimize the functional annotation of the sheep genome and facilitate improved mapping accuracy of genetic variant and expression data for traits in sheep.