The complete genome sequence of the African buffalo (Syncerus caffer)

Continent-wide genomic analysis of the African buffalo (Syncerus caffer)

The African buffalo (Syncerus caffer) is a wild bovid with a historical distribution across much of sub-Saharan Africa. Genomic analysis can provide insights into the evolutionary history of the species, and the key selective pressures shaping populations, including assessment of population level differentiation, population fragmentation, and population genetic structure. In this study we generated the highest quality de novo genome assembly (2.65 Gb, scaffold N50 69.17 Mb) of African buffalo to date, and sequenced a further 195 genomes from across the species distribution. Principal component and admixture analyses provided little support for the currently described four subspecies. Estimating Effective Migration Surfaces analysis suggested that geographical barriers have played a significant role in shaping gene flow and the population structure. Estimated effective population sizes indicated a substantial drop occurring in all populations 5-10,000 years ago, coinciding with the increase in human populations. Finally, signatures of selection were enriched for key genes associated with the immune response, suggesting infectious disease exert a substantial selective pressure upon the African buffalo. These findings have important implications for understanding bovid evolution, buffalo conservation and population management.

Research Progress and Applications of Bovine Genome in the Tribe Bovini

Various bovine species have been domesticated and bred for thousands of years, and they provide adequate animal-derived products, including meat, milk, and leather, to meet human requirements. Despite the review studies on economic traits in cattle, the genetic basis of traits has only been partially explained by phenotype and pedigree breeding methods, due to the complexity of genomic regulation during animal development and growth. With the advent of next-generation sequencing technology, genomics projects, such as the 1000 Bull Genomes Project, Functional Annotation of Animal Genomes project, and Bovine Pangenome Consortium, have advanced bovine genomic research. These large-scale genomics projects gave us a comprehensive concept, technology, and public resources. In this review, we summarize the genomics research progress of the main bovine species during the past decade, including cattle (Bos taurus), yak (Bos grunniens), water buffalo (Bubalus bubalis), zebu (Bos indicus), and gayal (Bos frontalis). We mainly discuss the development of genome sequencing and functional annotation, focusing on how genomic analysis reveals genetic variation and its impact on phenotypes in several bovine species.

The microscopic five of the big five: Managing zoonotic diseases within and beyond African wildlife protected areas

African protected areas strive to conserve the continent's great biodiversity with a targeted focus on the flagship 'Big Five' megafauna. Though often not considered, this biodiversity protection also extends to the lesser-known microbes and parasites that are maintained in these diverse ecosystems, often in a silent and endemically stable state. Climate and anthropogenic change, and associated diversity loss, however, are altering these dynamics leading to shifts in ecological interactions and pathogen spill over into new niches and hosts. As many African protected areas are bordered by game and livestock farms, as well as villages, they provide an ideal study system to assess infection dynamics at the human-livestock-wildlife interface. Here we review five zoonotic, multi-host diseases (bovine tuberculosis, brucellosis, Rift Valley fever, schistosomiasis and cryptosporidiosis)-the 'Microscopic Five'-and discuss the biotic and abiotic drivers of parasite transmission using the iconic Kruger National Park, South Africa, as a case study. We identify knowledge gaps regarding the impact of the 'Microscopic Five' on wildlife within parks and highlight the need for more empirical data, particularly for neglected (schistosomiasis) and newly emerging (cryptosporidiosis) diseases, as well as zoonotic disease risk from the rising bush meat trade and game farm industry. As protected areas strive to become further embedded in the socio-economic systems that surround them, providing benefits to local communities, One Health approaches can help maintain the ecological integrity of ecosystems, while protecting local communities and economies from the negative impacts of disease.

A draft genome of Drung cattle reveals clues to its chromosomal fusion and environmental adaptation

Drung cattle (Bos frontalis) have 58 chromosomes, differing from the Bos taurus 2n = 60 karyotype. To date, its origin and evolution history have not been proven conclusively, and the mechanisms of chromosome fusion and environmental adaptation have not been clearly elucidated. Here, we assembled a high integrity and good contiguity genome of Drung cattle with 13.7-fold contig N50 and 4.1-fold scaffold N50 improvements over the recently published Indian mithun assembly, respectively. Speciation time estimation and phylogenetic analysis showed that Drung cattle diverged from Bos taurus into an independent evolutionary clade. Sequence evidence of centromere regions provides clues to the breakpoints in BTA2 and BTA28 centromere satellites. We furthermore integrated a circulation and contraction-related biological process involving 43 evolutionary genes that participated in pathways associated with the evolution of the cardiovascular system. These findings may have important implications for understanding the molecular mechanisms of chromosome fusion, alpine valleys adaptability and cardiovascular function.

Kouprey (Bos sauveli) genomes unveil polytomic origin of wild Asian Bos

The evolution of the genera Bos and Bison, and the nature of gene flow between wild and domestic species, is poorly understood, with genomic data of wild species being limited. We generated two genomes from the likely extinct kouprey (Bos sauveli) and analyzed them alongside other Bos and Bison genomes. We found that B. sauveli possessed genomic signatures characteristic of an independent species closely related to Bos javanicus and Bos gaurus. We found evidence for extensive incomplete lineage sorting across the three species, consistent with a polytomic diversification of the major ancestry in the group, potentially followed by secondary gene flow. Finally, we detected significant gene flow from an unsampled Asian Bos-like source into East Asian zebu cattle, demonstrating both that the full genomic diversity and evolutionary history of the Bos complex has yet to be elucidated and that museum specimens and ancient DNA are valuable resources to do so.

•

We generated two genomes from the likely extinct kouprey (Bos sauveli)

•

Extensive mt and nuclear-genome-wide incomplete lineage sorting across wild Asian Bos

•

Initial polytomic diversification of the wild Asian Bos—kouprey, banteng, and gaur

The genome of the lowland anoa (Bubalus depressicornis) illuminates the origin of river and swamp buffalo

Two types of domestic water buffalo are currently recognized: the river buffalo from the Indian subcontinent and Mediterranean countries and the swamp buffalo from China and Southeast Asia. To test the hypothesis of two separate species of water buffalo, we sequenced the genome of the lowland anoa, Bubalus depressicornis, which is a dwarf wild buffalo endemic to Sulawesi, and two genomes of swamp buffalo, and made comparisons with 12 additional genomes. Three genomic data sets were constructed to infer phylogenetic relationships: the mitochondrial genome (15,468 bp; maternal transmission), two concatenated Y-chromosomal genes, AMELY and DDX3Y (20,036 bp; paternal transmission), and a selection of 30 nuclear genes representing all cattle chromosomes (364,887 bp; biparental transmission). The comparisons between our 30 nuclear gene sequences obtained by read mapping and those directly extracted from Bos taurus and Bubalus bubalis genome assemblies show that the mapping approach revealed higher levels of heterozygosity at both nucleotide sites and indels (insertions and deletions) (0.09-0.15%), as well as several sequence errors (0.07%). Our phylogenetic and molecular dating analyses provide strong evidence that the lowland anoa, river buffalo, and swamp buffalo are three distinct taxa which separated rapidly from each other during the Pleistocene epoch. We therefore conclude that two species of domestic water buffalo should be distinguished: Bubalus bubalis for the river buffalo and Bubalus kerabau for the swamp buffalo. The new classification can have deep implications for understanding the evolution and selection of domesticated forms and for the conservation and management of wild buffalo populations in South and Southeast Asia.

Whole-Genome Sequencing and Characterization of Buffalo Genetic Resources: Recent Advances and Future Challenges

The buffalo was domesticated around 3000-6000 years ago and has substantial economic significance as a meat, dairy, and draught animal. The buffalo has remained underutilized in terms of the development of a well-annotated and assembled reference genome de novo. It is mandatory to explore the genetic architecture of a species to understand the biology that helps to manage its genetic variability, which is ultimately used for selective breeding and genomic selection. Morphological and molecular data have revealed that the swamp buffalo population has strong geographical genomic diversity with low gene flow but strong phenotypic consistency, while the river buffalo population has higher phenotypic diversity with a weak phylogeographic structure. The availability of recent high-quality reference genome and genotyping marker panels has invigorated many genome-based studies on evolutionary history, genetic diversity, functional elements, and performance traits. The increasing molecular knowledge syndicate with selective breeding should pave the way for genetic improvement in the climatic resilience, disease resistance, and production performance of water buffalo populations globally.

Transcriptome analysis of buffalo granulosa cells in three dimensional culture systems

Hanging drop (HD) three-dimensional (3D) culture model for buffalo granulosa cells (GC) was reported to mimic the preovulatory stage of ovarian follicles in our previous study. To further verify its reliability, the present study attempted a comparative transcriptome profile of buffalo GC freshly isolated from ovarian follicles (<8 mm diameter) (FC) and their cultures in normal culture dish (ND or 2D), polyHEMA coated dish (PH) and HD culture systems (3D). Out of 223 significantly (-log2 fold change: >3; p < .0005; false discovery rate [FDR]: <0.1) differentially expressed genes (SDEGs) among different culture systems, 137 were found unannotated, and 94, 29, and 66 were exclusively expressed in FC, PH, and HD, respectively. However, on eliminating the fixed points of p values and FDR from the entire raw data, only 11 genes related to long noncoding RNA, 12 genes related to luteinization, and 3 genes related to follicular maturation were exclusively expressed in FC, PH, and HD culture systems, respectively. The quantitative real time-PCR validation and the next generation sequencing data had more than 90% correlation. Bioinformatics analyses of the exclusively expressed SDEG revealed that the freshly aspirated GCs were a true representative of GCs from small follicles (<8 mm diameter), the GC spheroids under PH maintained mitochondrial function, and those cultured in HD system for 6 days simulated the inflammatory milieu required for ovulation. Therefore, the comparative transcriptome profile also reinforced that HD culture system is better in vitro culture method than the other methods analyzed in this study for buffalo GC.

High diversity, inbreeding and a dynamic Pleistocene demographic history revealed by African buffalo genomes

Genomes retain records of demographic changes and evolutionary forces that shape species and populations. Remnant populations of African buffalo (Syncerus caffer) in South Africa, with varied histories, provide an opportunity to investigate signatures left in their genomes by past events, both recent and ancient. Here, we produce 40 low coverage (7.14×) genome sequences of Cape buffalo (S. c. caffer) from four protected areas in South Africa. Genome-wide heterozygosity was the highest for any mammal for which these data are available, while differences in individual inbreeding coefficients reflected the severity of historical bottlenecks and current census sizes in each population. PSMC analysis revealed multiple changes in N_e between approximately one million and 20 thousand years ago, corresponding to paleoclimatic changes and Cape buffalo colonisation of southern Africa. The results of this study have implications for buffalo management and conservation, particularly in the context of the predicted increase in aridity and temperature in southern Africa over the next century as a result of climate change.

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species - taurine cattle, zebu and yak - and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.

Sequencing and assembly of the Egyptian buffalo genome

Water buffalo (Bubalus bubalis) is an important source of meat and milk in countries with relatively warm weather. Compared to the cattle genome, a little has been done to reveal its genome structure and genomic traits. This is due to the complications stemming from the large genome size, the complexity of the genome, and the high repetitive content. In this paper, we introduce a high-quality draft assembly of the Egyptian water buffalo genome. The Egyptian breed is used as a dual purpose animal (milk/meat). It is distinguished by its adaptability to the local environment, quality of feed changes, as well as its high resistance to diseases. The genome assembly of the Egyptian water buffalo has been achieved using a reference-based assembly workflow. Our workflow significantly reduced the computational complexity of the assembly process, and improved the assembly quality by integrating different public resources. We also compared our assembly to the currently available draft assemblies of water buffalo breeds. A total of 21,128 genes were identified in the produced assembly. A list of milk virgin-related genes; milk pregnancy-related genes; milk lactation-related genes; milk involution-related genes; and milk mastitis-related genes were identified in the assembly. Our results will significantly contribute to a better understanding of the genetics of the Egyptian water buffalo which will eventually support the ongoing breeding efforts and facilitate the future discovery of genes responsible for complex processes of dairy, meat production and disease resistance among other significant traits.

Understanding divergent domestication traits from the whole-genome sequencing of swamp- and river-buffalo populations

Domesticated buffaloes have been integral to rice-paddy agro-ecosystems for millennia, yet relatively little is known about the buffalo genomics. Here, we sequenced and assembled reference genomes for both swamp and river buffaloes and we re-sequenced 230 individuals (132 swamp buffaloes and 98 river buffaloes) sampled from across Asia and Europe. Beyond the many actionable insights that our study revealed about the domestication, basic physiology and breeding of buffalo, we made the striking discovery that the divergent domestication traits between swamp and river buffaloes can be explained with recent selections of genes on social behavior, digestion metabolism, strengths and milk production.

Risk alleles for tuberculosis infection associate with reduced immune reactivity in a wild mammalian host

Integrating biological processes across scales remains a central challenge in disease ecology. Genetic variation drives differences in host immune responses, which, along with environmental factors, generates temporal and spatial infection patterns in natural populations that epidemiologists seek to predict and control. However, genetics and immunology are typically studied in model systems, whereas population-level patterns of infection status and susceptibility are uniquely observable in nature. Despite obvious causal connections, organizational scales from genes to host outcomes to population patterns are rarely linked explicitly. Here we identify two loci near genes involved in macrophage (phagocyte) activation and pathogen degradation that additively increase risk of bovine tuberculosis infection by up to ninefold in wild African buffalo. Furthermore, we observe genotype-specific variation in IL-12 production indicative of variation in macrophage activation. Here, we provide measurable differences in infection resistance at multiple scales by characterizing the genetic and inflammatory variation driving patterns of infection in a wild mammal.

Effect of Dioxins in Milk on the 3D Cultured Primary Buffalo Hepatocyte Model System

Cow and human milk have been reported to contain dioxins ranging from 0.023 to 26.46 and 0.88 to 19 pg/g of fat, respectively. However, the toxic effects of the dioxins in the milk in this range of concentrations were not explored. Therefore, considering the outbred livestock tissues as better models than inbred laboratory animals, the present study targeted to study the effect of dioxins present in the milk on three-dimensionally (3D) cultured buffalo primary hepatocyte spheroids. The spheroids were treated with a model dioxin, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), directly and also through milk fat at different concentrations (i.e, 0.02-20 pg/mL) for 24 h. Among the liver-cell-specific (ALB, HNF4α, and AFP) genes, a similar ALB and upregulated HNF4α expression at all treatments indicated the functional and transcriptionally active hepatocyte spheroids. Supportingly, no significant difference in the antiapoptotic gene expression between the treatments of milk fat and milk fat containing dioxins indicated the survivability of the spheroids during dioxin treatments. Among the selected TCDD responsive (CYP1A1, CYP1A2, AHR, CYP1B1, and TIPARP) genes, a nonsignificant increasing trend of the CYP1A1 expression was observed from 0.2 to 10 pg/mL of TCDD concentration through milk fat. This pattern was similar to the reported insensitive response of human primary hepatocytes toward dioxins than that of rat primary hepatocytes. This may indicate that the buffalo hepatocyte spheroids could be better models than rats for TCDD hepatotoxic studies. Further, TCDD in the milk in the range of 0.02-20 pg/mL concentration may not be very hepatotoxic.

Draft genome of the river water buffalo

Water buffalo (Bubalus bubalis), a large-sized member of the Bovidae family, is considered as an important livestock species throughout Southeast Asia. In order to better understand the molecular basis of buffalo improvement and breeding, we sequenced and assembled the genome (2n=50) of a river buffalo species Bubalus bubalis from Bangladesh. Its genome size is 2.77 Gb, with a contig N50 of 25 kb and the scaffold N50 of 6.9 Mbp. Based on the assembled genome, we annotated 24,613 genes for future functional genomics studies. Phylogenetic tree analysis of cattle and water buffalo lineages showed that they diverged about 5.8-9.8 million years ago. Our findings provide an insight into the water buffalo genome which will contribute in further research on buffalo such as molecular breeding, understanding complex traits, conservation, and biodiversity.

Context-dependent costs and benefits of tuberculosis resistance traits in a wild mammalian host

Disease acts as a powerful driver of evolution in natural host populations, yet individuals in a population often vary in their susceptibility to infection. Energetic trade-offs between immune and reproductive investment lead to the evolution of distinct life history strategies, driven by the relative fitness costs and benefits of resisting infection. However, examples quantifying the cost of resistance outside of the laboratory are rare. Here, we observe two distinct forms of resistance to bovine tuberculosis (bTB), an important zoonotic pathogen, in a free-ranging African buffalo (Syncerus caffer) population. We characterize these phenotypes as "infection resistance," in which hosts delay or prevent infection, and "proliferation resistance," in which the host limits the spread of lesions caused by the pathogen after infection has occurred. We found weak evidence that infection resistance to bTB may be heritable in this buffalo population (h 2 = 0.10) and comes at the cost of reduced body condition and marginally reduced survival once infected, but also associates with an overall higher reproductive rate. Infection-resistant animals thus appear to follow a "fast" pace-of-life syndrome, in that they reproduce more quickly but die upon infection. In contrast, proliferation resistance had no apparent costs and was associated with measures of positive host health-such as having a higher body condition and reproductive rate. This study quantifies striking phenotypic variation in pathogen resistance and provides evidence for a link between life history variation and a disease resistance trait in a wild mammalian host population.

Genome assembly and transcriptome resource for river buffalo, Bubalus bubalis (2n = 50)

Water buffalo is a globally important species for agriculture and local economies. A de novo assembled, well-annotated reference sequence for the water buffalo is an important prerequisite for studying the biology of this species, and is necessary to manage genetic diversity and to use modern breeding and genomic selection techniques. However, no such genome assembly has been previously reported. There are 2 species of domestic water buffalo, the river (2 n = 50) and the swamp (2 n = 48) buffalo. Here we describe a draft quality reference sequence for the river buffalo created from Illumina GA and Roche 454 short read sequences using the MaSuRCA assembler. The assembled sequence is 2.83 Gb, consisting of 366 983 scaffolds with a scaffold N50 of 1.41 Mb and contig N50 of 21 398 bp. Annotation of the genome was supported by transcriptome data from 30 tissues and identified 21 711 predicted protein coding genes. Searches for complete mammalian BUSCO gene groups found 98.6% of curated single copy orthologs present among predicted genes, which suggests a high level of completeness of the genome. The annotated sequence is available from NCBI at accession GCA_000471725.1.

The La and related RNA-binding proteins (LARPs): structures, functions, and evolving perspectives

La was first identified as a polypeptide component of ribonucleic protein complexes targeted by antibodies in autoimmune patients and is now known to be a eukaryote cell-ubiquitous protein. Structure and function studies have shown that La binds to a common terminal motif, UUU-3'-OH, of nascent RNA polymerase III (RNAP III) transcripts and protects them from exonucleolytic decay. For precursor-tRNAs, the most diverse and abundant of these transcripts, La also functions as an RNA chaperone that helps to prevent their misfolding. Related to this, we review evidence that suggests that La and its link to RNAP III were significant in the great expansions of the tRNAomes that occurred in eukaryotes. Four families of La-related proteins (LARPs) emerged during eukaryotic evolution with specialized functions. We provide an overview of the high-resolution structural biology of La and LARPs. LARP7 family members most closely resemble La but function with a single RNAP III nuclear transcript, 7SK, or telomerase RNA. A cytoplasmic isoform of La protein as well as LARPs 6, 4, and 1 function in mRNA metabolism and translation in distinct but similar ways, sometimes with the poly(A)-binding protein, and in some cases by direct binding to poly(A)-RNA. New structures of LARP domains, some complexed with RNA, provide novel insights into the functional versatility of these proteins. We also consider LARPs in relation to ancestral La protein and potential retention of links to specific RNA-related pathways. One such link may be tRNA surveillance and codon usage by LARP-associated mRNAs. WIREs RNA 2017, 8:e1430. doi: 10.1002/wrna.1430 For further resources related to this article, please visit the WIREs website.

Factors That Shape Eukaryotic tRNAomes: Processing, Modification and Anticodon-Codon Use

Transfer RNAs (tRNAs) contain sequence diversity beyond their anticodons and the large variety of nucleotide modifications found in all kingdoms of life. Some modifications stabilize structure and fit in the ribosome whereas those to the anticodon loop modulate messenger RNA (mRNA) decoding activity more directly. The identities of tRNAs with some universal anticodon loop modifications vary among distant and parallel species, likely to accommodate fine tuning for their translation systems. This plasticity in positions 34 (wobble) and 37 is reflected in codon use bias. Here, we review convergent evidence that suggest that expansion of the eukaryotic tRNAome was supported by its dedicated RNA polymerase III transcription system and coupling to the precursor-tRNA chaperone, La protein. We also review aspects of eukaryotic tRNAome evolution involving G34/A34 anticodon-sparing, relation to A34 modification to inosine, biased codon use and regulatory information in the redundancy (synonymous) component of the genetic code. We then review interdependent anticodon loop modifications involving position 37 in eukaryotes. This includes the eukaryote-specific tRNA modification, 3-methylcytidine-32 (m³C₃₂) and the responsible gene, TRM140 and homologs which were duplicated and subspecialized for isoacceptor-specific substrates and dependence on i⁶A₃₇ or t⁶A₃₇. The genetics of tRNA function is relevant to health directly and as disease modifiers.