Bos grunniens and Bos mutus (Artiodactyla: Bovidae)

Cathemerality: a key temporal niche

Given the marked variation in abiotic and biotic conditions between day and night, many species specialise their physical activity to being diurnal or nocturnal, and it was long thought that these strategies were commonly fairly fixed and invariant. The term 'cathemeral', was coined in 1987, when Tattersall noted activity in a Madagascan primate during the hours of both daylight and darkness. Initially thought to be rare, cathemerality is now known to be a quite widespread form of time partitioning amongst arthropods, fish, birds, and mammals. Herein we provide a synthesis of present understanding of cathemeral behaviour, arguing that it should routinely be included alongside diurnal and nocturnal strategies in schemes that distinguish and categorise species across taxa according to temporal niche. This synthesis is particularly timely because (i) the study of animal activity patterns is being revolutionised by new and improved technologies; (ii) it is becoming apparent that cathemerality covers a diverse range of obligate to facultative forms, each with their own common sets of functional traits, geographic ranges and evolutionary history; (iii) daytime and nighttime activity likely plays an important but currently neglected role in temporal niche partitioning and ecosystem functioning; and (iv) cathemerality may have an important role in the ability of species to adapt to human-mediated pressures.

Ancestral alleles defined for 70 million cattle variants using a population-based likelihood ratio test

BACKGROUND

The study of ancestral alleles provides insights into the evolutionary history, selection, and genetic structures of a population. In cattle, ancestral alleles are widely used in genetic analyses, including the detection of signatures of selection, determination of breed ancestry, and identification of admixture. Having a comprehensive list of ancestral alleles is expected to improve the accuracy of these genetic analyses. However, the list of ancestral alleles in cattle, especially at the whole genome sequence level, is far from complete. In fact, the current largest list of ancestral alleles (~ 42 million) represents less than 28% of the total number of detected variants in cattle. To address this issue and develop a genomic resource for evolutionary studies, we determined ancestral alleles in cattle by comparing prior derived whole-genome sequence variants to an out-species group using a population-based likelihood ratio test.

RESULTS

Our study determined and makes available the largest list of ancestral alleles in cattle to date (70.1 million) and includes 2.3 million on the X chromosome. There was high concordance (97.6%) of the determined ancestral alleles with those from previous studies when only high-probability ancestral alleles were considered (29.8 million positions) and another 23.5 million high-confidence ancestral alleles were novel, expanding the available reference list to improve the accuracies of genetic analyses involving ancestral alleles. The high concordance of the results with previous studies implies that our approach using genomic sequence variants and a likelihood ratio test to determine ancestral alleles is appropriate.

CONCLUSIONS

Considering the high concordance of ancestral alleles across studies, the ancestral alleles determined in this study including those not previously listed, particularly those with high-probability estimates, may be used for further genetic analyses with reasonable accuracy. Our approach that used predetermined variants in species and the likelihood ratio test to determine ancestral alleles is applicable to other species for which sequence level genotypes are available.

Disruption of trait-environment relationships in African megafauna occurred in the middle Pleistocene

Mammalian megafauna have been critical to the functioning of Earth's biosphere for millions of years. However, since the Plio-Pleistocene, their biodiversity has declined concurrently with dramatic environmental change and hominin evolution. While these biodiversity declines are well-documented, their implications for the ecological function of megafaunal communities remain uncertain. Here, we adapt ecometric methods to evaluate whether the functional link between communities of herbivorous, eastern African megafauna and their environments (i.e., functional trait-environment relationships) was disrupted as biodiversity losses occurred over the past 7.4 Ma. Herbivore taxonomic and functional diversity began to decline during the Pliocene as open grassland habitats emerged, persisted, and expanded. In the mid-Pleistocene, grassland expansion intensified, and climates became more variable and arid. It was then that phylogenetic diversity declined, and the trait-environment relationships of herbivore communities shifted significantly. Our results divulge the varying implications of different losses in megafaunal biodiversity. Only the losses that occurred since the mid-Pleistocene were coincident with a disturbance to community ecological function. Prior diversity losses, conversely, occurred as the megafaunal species and trait pool narrowed towards those adapted to grassland environments.

First Insights into the Occurrence of Circular Single-Stranded DNA Genomes in Asian and African Cattle

Circular replicase-encoding single-stranded (CRESS) DNA viruses and other circular DNA agents are increasingly found in various samples and animals. A specific class of these agents-termed bovine meat and milk factors (BMMF)-has been supposed to act as a factor in indirect carcinogenesis in humans. Initial observations attributed the BMMF to European cattle breeds and foodstuffs produced thereof. In the present study, blood and fecal samples from African and Asian cattle were examined. BMMF molecules and genomoviruses were detected in all bovids under study. The majority (79%) of the 29 circular elements could be assigned to BMMF groups 1 and 2, whereas CRESS viruses of the family Genomoviridae accounted for the smaller part (21%). Two genomoviruses belong to the genus Gemykibivirus and one to the genus Gemykrogvirus. The remaining three might be considered as novel species within the genus Gemycircularvirus. The majority of all isolated molecules originated from fecal samples, whereas only three derived from blood. The results from this study expand our knowledge on the diversity and presence of circular DNA in different ruminants that serve for food production in many countries over the world.

Functional substitutability of native herbivores by livestock for soil carbon stock is mediated by microbial decomposers

Grazing by large mammalian herbivores impacts climate as it can favor the size and stability of a large carbon (C) pool in the soils of grazing ecosystems. As native herbivores in the world's grasslands, steppes, and savannas are progressively being displaced by livestock, it is important to ask whether livestock can emulate the functional roles of their native counterparts. While livestock and native herbivores can have remarkable similarity in their traits, they can differ greatly in their impacts on vegetation composition which can affect soil-C. It is uncertain how these similarities and differences impact soil-C via their influence on microbial decomposers. We test competing alternative hypotheses with a replicated, long-term, landscape-level, grazing-exclusion experiment to ask whether livestock in the Trans-Himalayan ecosystem of northern India can match decadal-scale (2005-2016) soil-C stocks under native herbivores. We evaluate multiple lines of evidence from 17 variables that influence soil-C (quantity and quality of C-input from plants, microbial biomass and metabolism, microbial community composition, eDNA, veterinary antibiotics in soil), and assess their inter-relationships. Livestock and native herbivores differed in their effects on several soil microbial processes. Microbial carbon use efficiency (CUE) was 19% lower in soils under livestock. Compared to native herbivores, areas used by livestock contained 1.5 kg C m^-2 less soil-C. Structural equation models showed that alongside the effects arising from plants, livestock alter soil microbial communities which is detrimental for CUE, and ultimately also for soil-C. Supporting evidence pointed toward a link between veterinary antibiotics used on livestock, microbial communities, and soil-C. Overcoming the challenges of sequestering antibiotics to minimize their potential impacts on climate, alongside microbial rewilding under livestock, may reconcile the conflicting demands from food-security and ecosystem services. Conservation of native herbivores and alternative management of livestock is crucial for soil-C stewardship to envision and achieve natural climate solutions.

Characterization of olfactory receptor repertoires provides insights into the high-altitude adaptation of the yak based on the chromosome-level genome

Olfaction in vertebrates plays pivotal parts in many aspects, such as localizing prey or food, mating behavior, avoiding predators, and social communication. Yak (Bos grunniens) is the only Bos species that can thrive in high-altitude areas. In view of the critical role of olfactory receptors (ORs) in the specific recognition of diverse stimuli, investigating the evolutionary dynamics of ORs in the yak means a lot. In this study, we used the chromosome-level genome of the yak to identify the ORs genes and discussed the effects of high altitude on the yak's olfaction by comparing the yak with other low-altitude living Bos species (Bos frontalis (gayal), Bos gaurus (gaur), Bos indicus (zebu) and Bos taurus (cattle)). The yak had 400 OR genes, including 264 functional genes, 16 partial genes and 120 OR pseudo genes. There were 387 OR genes mapped to yak 31 chromosomes, and chromosomes 13 and 8 had the most OR genes and functional OR genes. Among these five Bos species, yak had the least number of OR gene subfamilies, OR genes and functional OR genes, while the total number of OR genes in gayal (n = 784) was almost twice as many as that of yak, indicating that the olfaction of yak may be less developed. In addition, the phylogenetic relationships of the functional Bos OR genes were illustrated, which comprised 79 families and 466 subfamilies distributed in two classes (Class I and Class II). There were 76 OR gene subfamilies shared by these five Bos species and 17 OR gene subfamilies were unique to the yak. The potential odor specificity of 44 yak OR genes was identified through the similarity to human OR protein sequences. Remarkably, yak lacks β-ionone and Isovaleric acid(IVA)-related ORs, which may be related to the decline of high-altitude herbaceous plant diversity and underdeveloped yak sweat glands. The conserved motifs of OR genes were highly conserved in Bos species. These results provided a solid foundation for further studies on the molecular mechanisms of the yak's adaptation to the high-altitude environment in olfaction.

Characteristics of pulmonary microvascular structure in postnatal yaks

Yaks are typical plateau-adapted animals, however the microvascular changes and characteristics in their lungs after birth are still unclear. Pulmonary microvasculature characteristics and changes across age groups were analysed using morphological observation and molecular biology detection in yaks aged 1, 30 and 180 days old in addition to adults. Results: Our experiments demonstrated that yaks have fully developed pulmonary alveolar at birth but that interalveolar thickness increased with age. Immunofluorescence observations showed that microvessel density within the interalveolar septum in the yak gradually increased with age. In addition, transmission electron microscopy (TEM) results showed that the blood–air barrier of 1-day old and 30-days old yaks was significantly thicker than that observed at 180-days old and in adults (P < 0.05), which was caused by the thinning of the membrane of alveolar epithelial cells. Furthermore, Vegfa and Epas1 expression levels in 30-day old yaks were the highest in comparison to the other age groups (P < 0.05), whilst levels in adult yaks were the lowest (P < 0.05). The gradual increase in lung microvessel density can effectively satisfy the oxygen requirements of ageing yaks. In addition, these results suggest that the key period of yak lung development is from 30 to 180 days.

Fat Deposition in the Muscle of Female and Male Yak and the Correlation of Yak Meat Quality with Fat

This study aimed to explore the differences in fat deposition between female (FYs) and male yaks (MYs). Compared with MYs, the tenderness, L*, marbling, absolute content of fat, and most fatty acids (FAs) of longissimus dorsi (LD) in FYs were higher or better (p < 0.05), whereas the relative content of polyunsaturated fatty acids (PUFAs) and n-3 PUFAs were lower (p < 0.01). The absolute content of fat, C18:0, cis-C18:2, cis-C18:1, and C24:0 were positively correlated with L*_{45 min}, b*_{24 h}, tenderness, and marbling score of LD in FYs and MYs (p < 0.05), respectively. LPL, FATP2, ELOVL6, HADH, HACD, and PLINS genes play a crucial role in improving the marbling score and tenderness of yak meat. The results of gene expression and protein synthesis showed the effect of gender to FA biosynthesis, FA transport, lipolysis, and FA oxidation in the adipose tissue of yak was realized by the expressions of ME1, SCD, ACSL5, LPL, FABP1, PLIN4, and PLIN2 in peroxisome proliferators-activated receptor (PPAR) signaling. This study established a theoretical basis for the improvement of the meat quality of yak and molecular breeding.

The complete mitochondrial genome and phylogenetic analysis of Yanglong yak (Bos grunniens)

In this study, we assembled the complete mitochondrial genome of Yanglong yak (Bos grunniens) from Illumina sequencing reads. The mitochondrial genome is 16,323 bp long with an A + T-biased nucleotide composition, and encodes 13 protein-coding, 22 tRNA, and two rRNA genes along with a noncoding control region. In addition, its gene order is identical to those of the previously published mitochondrial genomes of its congeners. Phylogenetic analysis indicates that this breed is closely related to Datong yak, Pamir yak, Tianjun yak, polled yak, Seron yak, Sunnan yak, a series of Domestic Yak and wild yak, followed by Jinchuan yak and Gannan yak, and slightly far away from Bison and Bos taurus.

Functional traits of the world's late Quaternary large-bodied avian and mammalian herbivores

Prehistoric and recent extinctions of large-bodied terrestrial herbivores had significant and lasting impacts on Earth's ecosystems due to the loss of their distinct trait combinations. The world's surviving large-bodied avian and mammalian herbivores remain among the most threatened taxa. As such, a greater understanding of the ecological impacts of large herbivore losses is increasingly important. However, comprehensive and ecologically-relevant trait datasets for extinct and extant herbivores are lacking. Here, we present HerbiTraits, a comprehensive functional trait dataset for all late Quaternary terrestrial avian and mammalian herbivores ≥10 kg (545 species). HerbiTraits includes key traits that influence how herbivores interact with ecosystems, namely body mass, diet, fermentation type, habitat use, and limb morphology. Trait data were compiled from 557 sources and comprise the best available knowledge on late Quaternary large-bodied herbivores. HerbiTraits provides a tool for the analysis of herbivore functional diversity both past and present and its effects on Earth's ecosystems.

Reproduction in female wild cattle: Influence of seasonality on ARTs

Wild cattle species, often considered less alluring than certain conservation-dependent species, have not attracted the same level of interest as the charismatic megafauna from the general public, private or corporate donors, and other funding agencies. Currently, most wild cattle populations are vulnerable or threatened with extinction. The implementation of reproductive technologies to maintain genetically healthy cattle populations in situ and ex situ has been considered for more than 30 years. Protocols developed for domestic cattle breeds have been used with some success in various wild cattle species. However, inherent differences in the natural life history of these species makes extrapolation of domestic cattle protocols difficult, and in some cases, minimally effective. Reproductive seasonality, driven by either photoperiod or nutritional resource availability, has significant influence on the success of assisted reproductive technologies (ARTs). This review focuses on the physiological processes that differ in breeding (ovulatory) and non-breeding (anovulatory) seasons in female cattle, and the potential methods used to overcome these challenges. Techniques to be discussed within the context of seasonality include: estrus synchronization and ovulation induction, ovarian superstimulation, artificial insemination (AI), multiple ovulation embryo transfer (MOET), and ovum pick-up (OPU) with in vitro fertilization (IVF) and embryo transfer (ET).

Combined dental wear and cementum analyses in ungulates reveal the seasonality of Neanderthal occupations in Covalejos Cave (Northern Iberia)

We propose for the first time the use of the combination of two high-resolution techniques, dental wear (meso- and microwear) and dental cementum analyses, to gain a better understanding of Neanderthal subsistence strategies and occupational patterns. Dental wear analysis provides information not only on ungulate palaeodiet and palaeoenvironments but also on hunting time and seasons. Dental cementum analysis allows the accurate determination of the age and season at death of a prey. Our study has focused on the Cantabrian region and has applied both methods to investigate the Mousterian faunal assemblages in Covalejos Cave. Identification of the ungulate palaeodiet reveals information on the environmental conditions of the studied region. Moreover, it may facilitate observation on the evolution of both palaeodiet and palaeoenvironment throughout the site sequence. Results show a general stability in the palaeoenvironmental conditions and in the ungulate palaeodiet throughout the Mousterian sequence; this finding may be attributed to the role of the area as a climate refuge, and slight differences in levels 8, 7 and 4 suggest long- or short-term but repeated Neanderthal occupations at different seasons in the annual cycle.

Characterization of the mitochondrial genome of the Qilian yak (Bos grunniens) with a phylogenetic analysis of the family Bovidae (Artiodactyla)

Qilian yak (Bos grunniens) is a local breed of yak with high adaptation to the high-elevation, cold, and anoxic environments. In this study, its complete mitochondrial genome was assembled from high-throughput sequencing reads. The mitogenome is 16,324 bp long with an A + T-biased nucleotide composition and has the typical set of 37 animal mitochondrial genes and the noncoding control region. Phylogenetic analysis supported the monophyly of the subfamilies Aepycerotinae, Alcelaphinae, Bovinae, Caprinae, and Cephalophinae, confirmed the polyphyly of the family Antilopinae within the family Bovidae, and indicated the close relatedness between domestic yaks (Bos grunniens) and wild yaks (Bos mutus).

Phylogeny of anaerobic fungi (phylum Neocallimastigomycota), with contributions from yak in China

The phylum Neocallimastigomycota contains eight genera (about 20 species) of strictly anaerobic fungi. The evolutionary relationships of these genera are uncertain due to insufficient sequence data to infer their phylogenies. Based on morphology and molecular phylogeny, thirteen isolates obtained from yak faeces and rumen digesta in China were assigned to Neocallimastix frontalis (nine isolates), Orpinomyces joyonii (two isolates) and Caecomyces sp. (two isolates), respectively. The phylogenetic relationships of the eight genera were evaluated using complete ITS and partial LSU sequences, compared to the ITS1 region which has been widely used in this phylum in the past. Five monophyletic lineages corresponding to six of the eight genera were statistically supported. Isolates of Caecomyces and Cyllamyces were present in a single lineage and could not be separated properly. Members of Neocallimastigomycota with uniflagellate zoospores represented by Piromyces were polyphyletic. The Piromyces-like genus Oontomyces was consistently closely related to the traditional Anaeromyces, and separated the latter genus into two clades. The phylogenetic position of the Piromyces-like genus Buwchfawromyces remained unresolved. Orpinomyces and Neocallimastix, sharing polyflagellate zoospores, were supported as sister genera in the LSU phylogeny. Apparently ITS, specifically ITS1 alone, is not a good marker to resolve the generic affinities of the studied fungi. The LSU sequences are easier to align and appear to work well to resolve generic relationships. This study provides a comparative phylogenetic revision of Neocallimastigomycota isolates known from culture and sequence data.

Understanding habitat selection of the Vulnerable wild yak Bos mutus on the Tibetan Plateau

We tested a series of hypotheses on drivers of habitat selection by the Vulnerable wild yak Bos mutus, combining distribution-wide sighting data with species distribution modelling approaches. The results indicate that climatic conditions are of paramount importance in shaping the wild yak's distribution on the Tibetan Plateau. Habitat selection patterns were seasonal, with yaks appearing to select areas closer to villages during the vegetation-growing season. Unexpectedly, our index of forage quantity had a limited effect in determining the distribution of the species. Overall, our results suggest that expected changes in climate for this region could have a significant impact on habitat availability for wild yaks, and we call for more attention to be focused on the unique wildlife in this ecosystem.

Legacies of past exploitation and climate affect mammalian sexes differently on the roof of the world - the case of wild yaks

In polar environments, a lack of empirical knowledge about biodiversity prompts reliance on species distribution models to predict future change, yet these ignore the role of biotic interactions including the role of long past human exploitation. To explore how mammals of extreme elevation respond to glacial recession and past harvest, we combined our fieldwork with remote sensing and used analyses of ~60 expeditions from 1850–1925 to represent baseline conditions for wildlife before heavy exploitation on the Tibetan Plateau. Focusing on endangered wild yaks (Bos mutus), we document female changes in habitat use across time whereupon they increasingly relied on steeper post-glacial terrain, and currently have a 20x greater dependence on winter snow patches than males. Our twin findings—that the sexes of a cold-adapted species respond differently to modern climate forcing and long-past exploitation—indicate that effective conservation planning will require knowledge of the interplay between past and future if we will assure persistence of the region's biodiversity.

Complete mitochondrial genome of the Yakow (Bos primigenius taurus × Bos grunniens) in China

The Chinese Yakow, Bos primigenius taurus × Bos grunniens, is a large and commercially important hybrid in family Bovidae. We first determined and annotated its complete mitochondrial genome. The mitogenome is 16,322 bp in length, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes and a control region. As in other mammals, most mitochondrial genes are encoded on the heavy strand, except for ND6 and eight tRNA genes, which are encoded on the light strand. Its overall base composition is A: 33.7%, T: 27.3%, C: 25.8% and G: 13.2%. The complete mitogenome of Yakow (B. p. taurus × B. grunniens) could provide an important data to further explore the taxonomic status of Yakow (B. p. taurus × B. grunniens) in B. grunniens and enrich the genetic information or evolutionary history of the Bovidae.

Complete mitochondrial genome of the a rare subspecies of genus Bos, Tianzhu white yak from Tibetan area in China

The population of domestic yak, Tianzhu white yak, from Tibetan area in China is considered as a rare Bos grunniens species. We first determined and annotated its complete mitochondrial genome. The mitogenome is 16,319 bp in length, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes and a control region. As in other mammals, most mitochondrial genes are encoded on the heavy strand, except for ND6 and eight tRNA genes, which are encoded on the light strand. Its overall base composition is A: 33.7%, T: 27.2%, C: 25.8% and G: 13.2%. The complete mitogenome of the new subspecies of Bos grunniens could provide an important data to further explore the taxonomic status of the subspecies.

Modeling trophic resource availability for the first human settlers of Europe: the case of Atapuerca TD6

Food resource availability strongly influences the survival opportunities of all organisms. The effect of animal food resource availability on the survival and dispersal of hominin populations is hotly debated. In this article, we present a mathematical model that provides estimations of the maximum and minimum available resources for secondary consumers in a palaeocommunity. This model provides insights into the intensity of competition and the available niche space for hominins in Europe during the early Galerian (1.2-0.8 Ma). Published data from the Atapuerca TD6 assemblage were used in combination with the model to investigate trophic dynamics and resource availability for a Homo antecessor population 800,000 years ago. The effect on our results of the possible presence at Atapuerca of some large carnivores not recorded in the fossil assemblage is also evaluated. Results indicate the existence of a rich ecosystem at Atapuerca at the end of the Early Pleistocene. Secondary production was abundant enough to maintain a hunter-gatherer population and a rich carnivore guild more diverse than that recorded in the TD6 assemblage. Based on these results, the practice of cannibalism by H. antecessor cannot be explained by a long-term scarcity of resources. High food availability at TD6 implies a low to moderate level of competition for resources between carnivores and humans. According to this interpretation, an empty niche for a highly carnivorous omnivore existed in Europe during the early Galerian, and it was successfully exploited by Homo.

Combining multiple autosomal introns for studying shallow phylogeny and taxonomy of Laurasiatherian mammals: Application to the tribe Bovini (Cetartiodactyla, Bovidae)

Mitochondrial sequences are widely used for species identification and for studying phylogenetic relationships among closely related species or populations of the same species. However, many studies of mammals have shown that the maternal history of the mitochondrial genome can be discordant with the true evolutionary history of the taxa. In such cases, the analyses of multiple nuclear genes can be more powerful for deciphering interspecific relationships. Here, we designed primers for amplifying 13 new exon-primed intron-crossing (EPIC) autosomal loci for studying shallow phylogeny and taxonomy of Laurasiatherian mammals. Three criteria were used for the selection of the markers: gene orthology, a PCR product length between 600 and 1200 nucleotides, and different chromosomal locations in the bovine genome. Positive PCRs were obtained from different species representing the orders Carnivora, Cetartiodactyla, Chiroptera, Perissodactyla and Pholidota. The newly developed markers were analyzed in a phylogenetic study of the tribe Bovini (the group containing domestic and wild cattle, bison, yak, African buffalo, Asian buffalo, and saola) based on 17 taxa and 18 nuclear genes, representing a total alignment of 13,095 nucleotides. The phylogenetic results were compared to those obtained from analyses of the complete mitochondrial genome and Y chromosomal genes. Our analyses support a basal divergence of the saola (Pseudoryx) and a sister-group relationship between yak and bison. These results contrast with recent molecular studies but are in better agreement with morphology. The comparison of pairwise nucleotide distances shows that our nuDNA dataset provides a good signal for identifying taxonomic levels, such as species, genera, subtribes, tribes and subfamilies, whereas the mtDNA genome fails because of mtDNA introgression and higher levels of homoplasy. Accordingly, we conclude that the genus Bison should be regarded as a synonym of Bos, with the European bison relegated to a subspecies rank within Bos bison. We compared our molecular dating estimates to the fossil record in order to propose a biogeographic scenario for the evolution of Bovini during the Neogene.

Out of Tibet: Pliocene woolly rhino suggests high-plateau origin of Ice Age megaherbivores

Ice Age megafauna have long been known to be associated with global cooling during the Pleistocene, and their adaptations to cold environments, such as large body size, long hair, and snow-sweeping structures, are best exemplified by the woolly mammoths and woolly rhinos. These traits were assumed to have evolved as a response to the ice sheet expansion. We report a new Pliocene mammal assemblage from a high-altitude basin in the western Himalayas, including a primitive woolly rhino. These new Tibetan fossils suggest that some megaherbivores first evolved in Tibet before the beginning of the Ice Age. The cold winters in high Tibet served as a habituation ground for the megaherbivores, which became preadapted for the Ice Age, successfully expanding to the Eurasian mammoth steppe.