The sheep genome illuminates biology of the rumen and lipid metabolism

Genomic insights into ruminant evolution: from past to future prospects

Ruminants (Ruminantia) are among the most successful herbivorous mammals, exhibiting wide-ranging morphological and ecological characteristics (such as headgear and multichambered stomach) and including various key livestock species (e.g., cattle, buffalo, yak, sheep, and goat). Understanding their evolution is of great significance not only in scientific research but also in applications potential for human society. The rapid growth of genomic resources provides unprecedented opportunities to dissect the evolutionary histories and molecular mechanisms underlying the distinct characteristics of ruminants. Here we summarize our current understanding of the genetic, morphological, and ecological diversity of ruminants and provide prospects for future studies.

From our roots, we grow

Located at the head of the Indo-Burma biodiversity hotspot, the Kunming Institute of Zoology (KIZ), Chinese Academy of Sciences (CAS), serves as China's main center for research into the diverse animal and ecological resources of southwestern China, Eastern Himalayas, and Southeast Asia. As of October 2019, it has been 60 years since the inception of KIZ. Since 1959, strong roots have been laid down by generations of researchers, allowing KIZ to grow and evolve into a comprehensive research institution renowned for its remarkable achievements in evolutionary mechanisms of animal biodiversity, animal resources protection, and sustainable utilization. It is now recognized as "a major powerhouse in evolutionary biology research in China" and is"establishing itself in the world stage" (Overseas Experts Review Committee, organized by the Bureau of Development Planning, CAS, during international evaluation in 2014).To celebrate the 60th anniversary of KIZ and the 70th anniversary of CAS, Zoological Research presents this commemorative issue, composed primarily of contributions from KIZ researchers. In addition, it is our great honor to provide here a brief retrospective of the pioneering work undertaken by the earlier scientists at KIZ and recent achievements, which will hopefully serve to motivate and inspire present and future successors.

Towards the Complete Goat Pan-Genome by Recovering Missing Genomic Segments From the Reference Genome

It is broadly expected that next generation sequencing will ultimately generate a complete genome as is the latest goat reference genome (ARS1), which is considered to be one of the most continuous assemblies in livestock. However, the rich diversity of worldwide goat breeds indicates that a genome from one individual would be insufficient to represent the whole genomic contents of goats. By comparing nine de novo assemblies from seven sibling species of domestic goat with ARS1 and using resequencing and transcriptome data from goats for verification, we identified a total of 38.3 Mb sequences that were absent in ARS1. The pan-sequences contain genic fractions with considerable expression. Using the pan-genome (ARS1 together with the pan-sequences) as a reference genome, variation calling efficacy can be appreciably improved. A total of 56,657 spurious SNPs per individual were repressed and 24,414 novel SNPs per individual on average were recovered as a result of better reads mapping quality. The transcriptomic mapping rate was also increased by ∼1.15%. Our study demonstrated that comparing de novo assemblies from closely related species is an efficient and reliable strategy for finding missing sequences from the reference genome and could be applicable to other species. Pan-genome can serve as an improved reference genome in animals for a better exploration of the underlying genomic variations and could increase the probability of finding genotype-phenotype associations assessed by a comprehensive variation database containing much more differences between individuals. We have constructed a goat pan-genome web interface for data visualization (http://animal.nwsuaf.edu.cn/panGoat).

The genetic architecture of helminth-specific immune responses in a wild population of Soay sheep (Ovis aries)

Much of our knowledge of the drivers of immune variation, and how these responses vary over time, comes from humans, domesticated livestock or laboratory organisms. While the genetic basis of variation in immune responses have been investigated in these systems, there is a poor understanding of how genetic variation influences immunity in natural, untreated populations living in complex environments. Here, we examine the genetic architecture of variation in immune traits in the Soay sheep of St Kilda, an unmanaged population of sheep infected with strongyle gastrointestinal nematodes. We assayed IgA, IgE and IgG antibodies against the prevalent nematode Teladorsagia circumcincta in the blood plasma of > 3,000 sheep collected over 26 years. Antibody levels were significantly heritable (h² = 0.21 to 0.57) and highly stable over an individual’s lifespan. IgA levels were strongly associated with a region on chromosome 24 explaining 21.1% and 24.5% of heritable variation in lambs and adults, respectively. This region was adjacent to two candidate loci, Class II Major Histocompatibility Complex Transactivator (CIITA) and C-Type Lectin Domain Containing 16A (CLEC16A). Lamb IgA levels were also associated with the immunoglobulin heavy constant loci (IGH) complex, and adult IgE levels and lamb IgA and IgG levels were associated with the major histocompatibility complex (MHC). This study provides evidence of high heritability of a complex immunological trait under natural conditions and provides the first evidence from a genome-wide study that large effect genes located outside the MHC region exist for immune traits in the wild.

Understanding how immune responses vary in natural populations can give an insight into how infection affects the ability of hosts and parasites to survive and reproduce, and how this drives evolutionary and ecological dynamics. Yet, very little is known about how immune responses vary over an individual’s lifetime and how genes contribute to this variation under natural conditions. Our study investigates the genetic architecture of variation in three antibody types, IgA, IgE and IgG in a wild population of Soay sheep on the St Kilda archipelago in North-West Scotland. Using data collected over 26 years, we show that antibody levels have a heritable basis in lambs and adults and are stable over an individual’s lifetime. We also identify several genomic regions with large effects on immune responses. Our study offers the first insights into the genetic control of immunity in a wild population, which is essential to understand how immune profiles vary in challenging natural conditions and how natural selection maintains genetic variation in complex immune traits.

An evaluation of sequencing coverage and genotyping strategies to assess neutral and adaptive diversity

Whole genome sequences (WGS) greatly increase our ability to precisely infer population genetic parameters, demographic processes, and selection signatures. However, WGS may still be not affordable for a representative number of individuals/populations. In this context, our goal was to assess the efficiency of several SNP genotyping strategies by testing their ability to accurately estimate parameters describing neutral diversity and to detect signatures of selection. We analysed 110 WGS at 12× coverage for four different species, i.e., sheep, goats and their wild counterparts. From these data we generated 946 data sets corresponding to random panels of 1K to 5M variants, commercial SNP chips and exome capture, for sample sizes of five to 48 individuals. We also extracted low-coverage genome resequencing of 1×, 2× and 5× by randomly subsampling reads from the 12× resequencing data. Globally, 5K to 10K random variants were enough for an accurate estimation of genome diversity. Conversely, commercial panels and exome capture displayed strong ascertainment biases. Besides the characterization of neutral diversity, the detection of the signature of selection and the accurate estimation of linkage disequilibrium (LD) required high-density panels of at least 1M variants. Finally, genotype likelihoods increased the quality of variant calling from low coverage resequencing but proportions of incorrect genotypes remained substantial, especially for heterozygote sites. Whole genome resequencing coverage of at least 5× appeared to be necessary for accurate assessment of genomic variations. These results have implications for studies seeking to deploy low-density SNP collections or genome scans across genetically diverse populations/species showing similar genetic characteristics and patterns of LD decay for a wide variety of purposes.

A Population Genomics Analysis of the Native Irish Galway Sheep Breed

The Galway sheep population is the only native Irish sheep breed and this livestock genetic resource is currently categorised as 'at-risk'. In the present study, comparative population genomics analyses of Galway sheep and other sheep populations of European origin were used to investigate the microevolution and recent genetic history of the breed. These analyses support the hypothesis that British Leicester sheep were used in the formation of the Galway. When compared to conventional and endangered breeds, the Galway breed was intermediate in effective population size, genomic inbreeding and runs of homozygosity. This indicates that, although the Galway breed is declining, it is still relatively genetically diverse and that conservation and management plans informed by genomic information may aid its recovery. The Galway breed also exhibited distinct genomic signatures of artificial or natural selection when compared to other breeds, which highlighted candidate genes that may be involved in production and health traits.

Genetic Variability in the Loci of FABP4, PPARγ and SCD Genes of Sheep Breeds Raised for Different Purposes

The present study was undertaken to analyse the genetic variation in coding sequences, splicing sites and regulatory sequences of FABP4, PPARγ and SCD genes in five breeds of sheep raised in Poland with different purposes: meat (Suffolk, Ile de France), dairy/prolific (Olkuska sheep, Kołuda) and primitive breeds (Polish Mountain Sheep). To identify genetic variability HRM-PCR, sequencing and PCR-RFLP method were used. Determining the genetic structure of the mentioned genes revealed six SNPs: FABP4 – c.73 + 13A>T and c.348 + 85G>A, PPARγ – c.391 – 66C>T (c.481 – 66C>T) and c.529 + 27G>C (c.619 + 27G>C), SCD – c.*945G>A and c.*1116A>G. For the c.*1116A>G SNP, a potential association with specific sheep body type and breeding purpose was found. In turn, the substitution c.*945G>A located in the regulatory region of the 3'-UTR of the Ovis aries SCD gene was identified for the first time. Based on the relationship demonstrated between sheep body constitution and productive type (dairy, meat) and the polymorphism of the SCD gene, further research is needed. The correlation between c.*1116A>G polymorphism and growth rate, slaughter and carcass value as well as meat quality of lambs needs to be studied in more detail. Such studies may lead to more effective selection processes in sheep breeding in the future.

Elimination of Reference Mapping Bias Reveals Robust Immune Related Allele-Specific Expression in Crossbred Sheep

Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.

Whole Genome Sequencing of the Giant Grouper (Epinephelus lanceolatus) and High-Throughput Screening of Putative Antimicrobial Peptide Genes

Giant groupers, the largest grouper type in the world, are of economic importance in marine aquaculture for their rapid growth. At the same time, bacterial and viral diseases have become the main threats to the grouper industry. Here, we report a high-quality genome of a giant grouper sequenced by an Illumina HiSeq X-Ten and PacBio Bioscience Sequel platform. A total of 254 putative antimicrobial peptide (AMP) genes were identified, which can be divided into 34 classes according to the annotation of the Antimicrobial Peptides Database (APD3). Their locations in pseudochromosomes were also determined. Thrombin-, lectin-, and scolopendin-derived putative AMPs were the three largest parts. In addition, expressions of putative AMPs were measured by our transcriptome data. Two putative AMP genes (gapdh1 and gapdh2) were involved in glycolysis, which had extremely high expression levels in giant grouper muscle. As it has been reported that AMPs inhibit the growth of a broad spectrum of microbes and participate in regulating innate and adaptive immune responses, genome sequencing of this study provides a comprehensive cataloging of putative AMPs of groupers, supporting antimicrobial research and aquaculture therapy. These genomic resources will be beneficial to further molecular breeding of this economically important fish.

Gene expression analysis of ovine prepubertal testicular tissue vitrified with a novel cryodevice (E.Vit)

PURPOSE

Testicular tissue cryopreservation prior to gonadotoxic therapies is a method to preserve fertility in children. However, the technique still requires development, especially when the tissue is immature and rather susceptible to stress derived from in vitro manipulation. This study aimed to investigate the effects of vitrification with a new cryodevice (E.Vit) on cell membrane integrity and gene expression of prepubertal testicular tissue in the ovine model.

METHODS

Pieces of immature testicular tissue (1 mm³) were inserted into "E.Vit" devices and vitrified with a two-step protocol. After warming, tissues were cultured in vitro and cell membrane integrity was assessed after 0, 2, and 24 h by trypan blue exclusion test. Controls consisted of non-vitrified tissue analyzed after 0, 2, and 24 h in vitro culture (IVC). Expression of genes involved in transcriptional stress response (BAX, SOD1, CIRBP, HSP90AB1), cell proliferation (KIF11), and germ- (ZBDB16, TERT, POU5F1, KIT) and somatic- (AR, FSHR, STAR) cell specific markers was evaluated 2 and 24 h after warming.

RESULTS

Post-warming trypan blue staining showed the survival of most cells, although membrane integrity immediately after warming (66.00% ± 4.73) or after 2 h IVC (59.67% ± 4.18) was significantly lower than controls (C0h 89.67% ± 1.45). Extended post-warming IVC (24 h) caused an additional decrease to 31% ± 3.46 (P < 0.05). Germ- and somatic-cell specific markers showed the survival of both cell types after cryopreservation and IVC. All genes were affected by cryopreservation and/or IVC, and moderate stress conditions were indicated by transcriptional stress response.

CONCLUSIONS

Vitrification with the cryodevice E.Vit is a promising strategy to cryopreserve prepubertal testicular tissue.

Selection signature analysis reveals genes underlying sheep milking performance

Sheep milk is the most important feed resource for newborn lambs and an important food resource for humans. Sheep milk production and ingredients are influenced by genetic and environmental factors. In this study, we implemented selection signature analysis using Illumina Ovine SNP50 BeadChip data of 78 meat Lacaune and 103 milk Lacaune sheep, which have similar genetic backgrounds, from the Sheep HapMap project to identify candidate genes related to ovine milk traits. Since different methods can detect different variation types and complement each other, we used a haplotype-based method (hapFLK) to implement selection signature analysis. The results revealed six selection signature regions showing signs of being selected ( P < 0.001 ): chromosomes 1, 2, 3, 6, 13 and 18. In addition, 38 quantitative trait loci (QTLs) related to sheep milk performance were identified in selection signature regions, which contain 334 candidate genes. Of those, SUCNR1 (succinate receptor 1) and PPARGC1A (PPARG coactivator 1 alpha) may be the most significant genes that affect sheep milking performance, which supply a significant indication for future studies to investigate candidate genes that play an important role in milk production and quality.

Complete Inactivation of Sebum-Producing Genes Parallels the Loss of Sebaceous Glands in Cetacea

Genomes are dynamic biological units, with processes of gene duplication and loss triggering evolutionary novelty. The mammalian skin provides a remarkable case study on the occurrence of adaptive morphological innovations. Skin sebaceous glands (SGs), for instance, emerged in the ancestor of mammals serving pivotal roles, such as lubrication, waterproofing, immunity, and thermoregulation, through the secretion of sebum, a complex mixture of various neutral lipids such as triacylglycerol, free fatty acids, wax esters, cholesterol, and squalene. Remarkably, SGs are absent in a few mammalian lineages, including the iconic Cetacea. We investigated the evolution of the key molecular components responsible for skin sebum production: Dgat2l6, Awat1, Awat2, Elovl3, Mogat3, and Fabp9. We show that all analyzed genes have been rendered nonfunctional in Cetacea species (toothed and baleen whales). Transcriptomic analysis, including a novel skin transcriptome from blue whale, supports gene inactivation. The conserved mutational pattern found in most analyzed genes, indicates that pseudogenization events took place prior to the diversification of modern Cetacea lineages. Genome and skin transcriptome analysis of the common hippopotamus highlighted the convergent loss of a subset of sebum-producing genes, notably Awat1 and Mogat3. Partial loss profiles were also detected in non-Cetacea aquatic mammals, such as the Florida manatee, and in terrestrial mammals displaying specialized skin phenotypes such as the African elephant, white rhinoceros and pig. Our findings reveal a unique landscape of “gene vestiges” in the Cetacea sebum-producing compartment, with limited gene loss observed in other mammalian lineages: suggestive of specific adaptations or specializations of skin lipids.

Applications of the Soil, Plant and Rumen Microbiomes in Pastoral Agriculture

The production of dairy, meat, and fiber by ruminant animals relies on the biological processes occurring in soils, forage plants, and the animals' rumens. Each of these components has an associated microbiome, and these have traditionally been viewed as distinct ecosystems. However, these microbiomes operate under similar ecological principles and are connected via water, energy flows, and the carbon and nitrogen nutrient cycles. Here, we summarize the microbiome research that has been done in each of these three environments (soils, forage plants, animals' rumen) and investigate what additional benefits may be possible through understanding the interactions between the various microbiomes. The challenge for future research is to enhance microbiome function by appropriate matching of plant and animal genotypes with the environment to improve the output and environmental sustainability of pastoral agriculture.

Large-scale ruminant genome sequencing provides insights into their evolution and distinct traits

The ruminants are one of the most successful mammalian lineages, exhibiting morphological and habitat diversity and containing several key livestock species. To better understand their evolution, we generated and analyzed de novo assembled genomes of 44 ruminant species, representing all six Ruminantia families. We used these genomes to create a time-calibrated phylogeny to resolve topological controversies, overcoming the challenges of incomplete lineage sorting. Population dynamic analyses show that population declines commenced between 100,000 and 50,000 years ago, which is concomitant with expansion in human populations. We also reveal genes and regulatory elements that possibly contribute to the evolution of the digestive system, cranial appendages, immune system, metabolism, body size, cursorial locomotion, and dentition of the ruminants.

A high-density BAC physical map covering the entire MHC region of addax antelope genome

BACKGROUND

The mammalian major histocompatibility complex (MHC) harbours clusters of genes associated with the immunological defence of animals against infectious pathogens. At present, no complete MHC physical map is available for any of the wild ruminant species in the world.

RESULTS

The high-density physical map is composed of two contigs of 47 overlapping bacterial artificial chromosome (BAC) clones, with an average of 115 Kb for each BAC, covering the entire addax MHC genome. The first contig has 40 overlapping BAC clones covering an approximately 2.9 Mb region of MHC class I, class III, and class IIa, and the second contig has 7 BAC clones covering an approximately 500 Kb genomic region that harbours MHC class IIb. The relative position of each BAC corresponding to the MHC sequence was determined by comparative mapping using PCR screening of the BAC library of 192,000 clones, and the order of BACs was determined by DNA fingerprinting. The overlaps of neighboring BACs were cross-verified by both BAC-end sequencing and co-amplification of identical PCR fragments within the overlapped region, with their identities further confirmed by DNA sequencing.

CONCLUSIONS

We report here the successful construction of a high-quality physical map for the addax MHC region using BACs and comparative mapping. The addax MHC physical map we constructed showed one gap of approximately 18 Mb formed by an ancient autosomal inversion that divided the MHC class II into IIa and IIb. The autosomal inversion provides compelling evidence that the MHC organizations in all of the ruminant species are relatively conserved.

The genome resources for conservation of Indo-Pacific humpback dolphin, Sousa chinensis

The Indo-Pacific humpback dolphin (Sousa chinensis), is a threatened marine mammal and belongs to the First Order of the National Key Protected Wild Aquatic Animals List in China. However, limited genomic information is available for studies of its population genetics and biological conservation. Here, we have assembled a genomic sequence of this species using a whole genome shotgun (WGS) sequencing strategy after a pilot low coverage genome survey. The total assembled genome size was 2.34 Gb: with a contig N50 of 67 kb and a scaffold N50 of 9 Mb (107.6-fold sequencing coverage). The S. chinensis genome contained 24,640 predicted protein-coding genes and had approximately 37% repeated sequences. The completeness of the genome assembly was evaluated by benchmarking universal single copy orthologous genes (BUSCOs): 94.3% of a total 4,104 expected mammalian genes were identified as complete, and 2.3% were identified as fragmented. This newly produced high-quality assembly and annotation of the genome will greatly promote the future studies of the genetic diversity, conservation and evolution.

High-density genotyping reveals signatures of selection related to acclimation and economically important traits in 15 local sheep breeds from Russia

Background

Domestication and centuries of selective breeding have changed genomes of sheep breeds to respond to environmental challenges and human needs. The genomes of local breeds, therefore, are valuable sources of genomic variants to be used to understand mechanisms of response to adaptation and artificial selection. As a step toward this we performed a high-density genotyping and comprehensive scans for signatures of selection in the genomes from 15 local sheep breeds reared across Russia.

Results

Results demonstrated that the genomes of Russian sheep breeds contain multiple regions under putative selection. More than 50% of these regions matched with intervals identified in previous scans for selective sweeps in sheep genomes. These regions contain well-known candidate genes related to morphology, adaptation, and domestication (e.g., KITLG, KIT, MITF, and MC1R), wool quality and quantity (e.g., DSG@, DSC@, and KRT@), growth and feed intake (e.g., HOXA@, HOXC@, LCORL, NCAPG, LAP3, and CCSER1), reproduction (e.g., CMTM6, HTRA1, GNAQ, UBQLN1, and IFT88), and milk-related traits (e.g., ABCG2, SPP1, ACSS1, and ACSS2). In addition, multiple genes that are putatively related to environmental adaptations were top-ranked in selected intervals (e.g., EGFR, HSPH1, NMUR1, EDNRB, PRL, TSHR, and ADAMTS5). Moreover, we observed that multiple key genes involved in human hereditary sensory and autonomic neuropathies, and genetic disorders accompanied with an inability to feel pain and environmental temperatures, were top-ranked in multiple or individual sheep breeds from Russia pointing to a possible mechanism of adaptation to harsh climatic conditions.

Conclusions

Our work represents the first comprehensive scan for signatures of selection in genomes of local sheep breeds from the Russian Federation of both European and Asian origins. We confirmed that the genomes of Russian sheep contain previously identified signatures of selection, demonstrating the robustness of our integrative approach. Multiple novel signatures of selection were found near genes which could be related to adaptation to the harsh environments of Russia. Our study forms a basis for future work on using Russian sheep genomes to spot specific genetic variants or haplotypes to be used in efforts on developing next-generation highly productive breeds, better suited to diverse Eurasian environments.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5537-0) contains supplementary material, which is available to authorized users.

Genome-wide association studies revealed candidate genes for tail fat deposition and body size in the Hulun Buir sheep

Fat-tailed sheep have a unique characteristic of depositing fat in their tails. In the present study, we conducted genome-wide association studies (GWAS) on traits related to tail fat deposition and body size in the Hulun Buir sheep. A total number of 300 individuals belonging to two fat-tailed lines of the Hulun Buir sheep breed genotyped with the Ovine Infinium HD SNP BeadChip were included in the current study. Two mixed models, one for continuous and one for binary phenotypic traits, were employed to analyse ten traits, that is, body length (BL), body height (BH), chest girth (CG), tail length (TL), tail width (TW), tail circumference (TC), carcass weight (CW), tail fat weight (TF), ratio of CW to TF (RCT) and tail type (TT). We identified 7, 6, 7, 2, 10 and 1 SNPs significantly associated with traits TF, CW, RCT, TW, TT and CG, respectively. Their associated genomic regions harboured 42 positional candidate genes. Out of them, 13 candidate genes including SMURF2, FBF1, DTNBP1, SETD7 and RBM11 have been associated with fat metabolism in sheep. The RBM11 gene has already been identified in a previous study on signatures of selection in this specific sheep population. Two more genes, that is, SMARCA5 and GAB1 were associated with body size in sheep. The present study has identified candidate genes that might be implicated in tail fat deposition and body size in sheep.

Genomic Approaches Identify Novel Gene Associations with Out of Season Lambing in Sheep

Sheep are seasonally polyestrous, traditionally breeding when the day length shortens in the autumn. The changing photoperiod stimulates reproductive hormones through a series of chemical pathways, ultimately leading to cyclicity. Some breeds of sheep, such as the Polypay and Dorset, have been selected for reduced seasonality and can lamb year-round. Despite this selection, there is still variation within these breeds in the ability to lamb out of season. The identification of out of season lambing quantitative trait loci has the potential to improve genetic progress using genomic selection schemes. Association studies, fixation index (FST), and runs of homozygosity (ROH) were evaluated to identify regions of the genome that influence the ability of ewes to lamb out of season. All analyses used genotypic data from the Illumina Ovine HD beadchip. Genome-wide associations were tested both across breeds in 257 ewes and within the Dorset and Polypay breeds. FST was measured across breeds and between UK and US Dorsets to assess population differences. ROH were estimated in ewes to identify homozygous regions contributing to out of season lambing. Significant associations after multiple testing correction were found through these approaches, leading to the identification of several candidate genes for further study. Genes involved with eye development, reproductive hormones, and neuronal changes were identified as the most promising for influencing the ewe’s ability to lamb year-round. These candidate genes could be advantageous for selection for improved year-round lamb production and provide better insight into the complex regulation of seasonal reproduction.

Absence of a human ortholog of rodent Kupffer cell galactose-binding receptor encoded by the CLEC4f gene

The murine CLEC4f gene encodes the Kupffer cell receptor, a galactose-binding receptor containing a C-type carbohydrate-recognition domain. Orthologs have been identified in nearly 100 species. The receptors from rat and mouse have previously been characterized and data presented here show that functional CLEC4f protein is expressed in domestic cattle (Bos taurus). However, the human CLEC4f gene does not encode a functional receptor because a mutation in the splice acceptor site of the final exon prevents appropriate splicing and a missense mutation disrupts the sugar-binding site. Transcriptomic and PCR analysis of transcripts confirms the absence of a spliced transcript containing the final exon and only background levels of transcripts are detected in human tissues. These mutations are also present in the CLEC4f gene in Neanderthals. In contrast to humans, closely related species, including chimpanzees, do have CLEC4f genes that encode full-length receptors. Affinity chromatography and glycan array results demonstrate that the chimpanzee, bovine and murine proteins all bind to galactose, but they show preferences for different subsets of galactose-containing glycans. In non-human primates, the receptor is expressed in spleen rather than in liver. The results indicate that the CLEC4f protein probably has distinct functions in different species. Absence of the receptor precludes using it for targeting of glycoconjugates to cells in human liver. The fact that CLEC4f protein is expressed in spleen in non-human primates and the close evolutionary relationship of the CLEC4f protein to langerin (CD207) suggest that it may function in the immune system, possibly as a pathogen receptor.

Proposed U.S. regulation of gene-edited food animals is not fit for purpose

Dietary DNA is generally regarded as safe to consume, and is a routine ingredient of food obtained from any living organism. Millions of naturally-occurring DNA variations are observed when comparing the genomic sequence of any two healthy individuals of a given species. Breeders routinely select desired traits resulting from this DNA variation to develop new cultivars and varieties of food plants and animals. Regulatory agencies do not evaluate these new varieties prior to commercial release. Gene editing tools now allow plant and animal breeders to precisely introduce useful genetic variation into agricultural breeding programs. The U.S. Department of Agriculture (USDA) announced that it has no plans to place additional regulations on gene-edited plants that could otherwise have been developed through traditional breeding prior to commercialization. However, the U.S. Food and Drug Administration (FDA) has proposed mandatory premarket new animal drug regulatory evaluation for all food animals whose genomes have been intentionally altered using modern molecular technologies including gene editing technologies. This runs counter to U.S. biotechnology policy that regulatory oversight should be triggered by unreasonable risk, and not by the fact that an organism has been modified by a particular process or technique. Breeder intention is not associated with product risk. Harmonizing the regulations associated with gene editing in food species is imperative to allow both plant and animal breeders access to gene editing tools to introduce useful sustainability traits like disease resistance, climate adaptability, and food quality attributes into U.S. agricultural breeding programs.

Programmable Base Editing of the Sheep Genome Revealed No Genome-Wide Off-Target Mutations

Since its emergence, CRISPR/Cas9-mediated base editors (BEs) with cytosine deaminase activity have been used to precisely and efficiently introduce single-base mutations in genomes, including those of human cells, mice, and crop species. Most production traits in livestock are induced by point mutations, and genome editing using BEs without homology-directed repair of double-strand breaks can directly alter single nucleotides. The p.96R > C variant of Suppressor cytokine signaling 2 (SOCS2) has profound effects on body weight, body size, and milk production in sheep. In the present study, we successfully obtained lambs with defined point mutations resulting in a p.96R > C substitution in SOCS2 by the co-injection of BE3 mRNA and a single guide RNA (sgRNA) into sheep zygotes. The observed efficiency of the single nucleotide exchange in newborn animals was as high as 25%. Observations of body size and body weight in the edited group showed that gene modification contributes to enhanced growth traits in sheep. Moreover, targeted deep sequencing and unbiased family trio-based whole genome sequencing revealed undetectable off-target mutations in the edited animals. This study demonstrates the potential for the application of BE-mediated point mutations in large animals for the improvement of production traits in livestock species.

Comprehensive Transcriptional Profiling of the Gastrointestinal Tract of Ruminants from Birth to Adulthood Reveals Strong Developmental Stage Specific Gene Expression

One of the most significant physiological challenges to neonatal and juvenile ruminants is the development and establishment of the rumen. Using a subset of RNA-Seq data from our high-resolution atlas of gene expression in sheep (Ovis aries) we have provided the first comprehensive characterization of transcription of the entire gastrointestinal (GI) tract during the transition from pre-ruminant to ruminant. The dataset comprises 164 tissue samples from sheep at four different time points (birth, one week, 8 weeks and adult). Using network cluster analysis we illustrate how the complexity of the GI tract is reflected in tissue- and developmental stage-specific differences in gene expression. The most significant transcriptional differences between neonatal and adult sheep were observed in the rumen complex. Comparative analysis of gene expression in three GI tract tissues from age-matched sheep and goats revealed species-specific differences in genes involved in immunity and metabolism. This study improves our understanding of the transcriptomic mechanisms involved in the transition from pre-ruminant to ruminant by identifying key genes involved in immunity, microbe recognition and metabolism. The results form a basis for future studies linking gene expression with microbial colonization of the developing GI tract and provide a foundation to improve ruminant efficiency and productivity through identifying potential targets for novel therapeutics and gene editing.

Genome-wide association study of lung lesions and pleurisy in New Zealand lambs

Pneumonia is an important issue for sheep production, leading to reduced growth rate and a predisposition to pleurisy. The objective of this study was to identify loci associated with pneumonic lesions and pleurisy in New Zealand progeny test lambs. The lungs from 3,572 progeny-test lambs were scored for presence and severity of pneumonic lesions and pleurisy at slaughter. Animals were genotyped using the Illumina Ovine Infinium HD SNP BeadChip (606,006 markers). The heritability of lung lesion score and pleurisy were calculated using the genomic relationship matrix, and genome-wide association analyses were conducted using EMMAX and haplotype trend regression. At slaughter, 35% of lambs had pneumonic lesions, with 9% showing lesions on more than half of any individual lobe. The number of lambs recorded as having pleurisy by the processing plants was 9%. Heritability estimates for pneumonic lesions and pleurisy scores adjusted for heteroscedasticity (CPSa and PLEURa) were 0.16 (± 0.03) and 0.05 (± 0.02), respectively. Five single-nucleotide polymorphisms (SNPs) were significantly associated with pneumonic lesions at the genome-wide level, and additional 37 SNPs were suggestively significant. Four SNPs were significantly associated with pleurisy, with an additional 11 SNPs reaching the suggestive level of significance. There were no regions that overlapped between the 2 traits. Multiple SNPs were in regions that contained genes involved in either the DNA damage response or the innate immune response, including several that had previously been reported to have associations with respiratory disease. Both EMMAX and HTR analyses of pleurisy data showed a significant peak on chromosome 2, located downstream from the transcription factor SP3. SP3 activates or suppresses the expression of numerous genes, including several genes with known functions in the immune system. This study identified several SNPs associated with genes involved in both the innate immune response and the response to DNA damage that are associated with pneumonic lesions and pleurisy in lambs at slaughter. Additionally, the identification in sheep of several SNPs within genes that have previously been associated with the respiratory system in cattle, pigs, rats, and mice indicates that there may be common pathways that underlie the response to invasion by respiratory pathogens in multiple species.

The genetic legacy of 50 years of desert bighorn sheep translocations

Conservation biologists have increasingly used translocations to mitigate population declines and restore locally extirpated populations. Genetic data can guide the selection of source populations for translocations and help evaluate restoration success. Bighorn sheep (Ovis canadensis) are a managed big game species that suffered widespread population extirpations across western North America throughout the early 1900s. Subsequent translocation programs have successfully re-established many formally extirpated bighorn herds, but most of these programs pre-date genetically informed management practices. The state of Nevada presents a particularly well-documented case of decline followed by restoration of extirpated herds. Desert bighorn sheep (O. c. nelsoni) populations declined to less than 3,000 individuals restricted to remnant herds in the Mojave Desert and a few locations in the Great Basin Desert. Beginning in 1968, the Nevada Department of Wildlife translocated ~2,000 individuals from remnant populations to restore previously extirpated areas, possibly establishing herds with mixed ancestries. Here, we examined genetic diversity and structure among remnant herds and the genetic consequences of translocation from these herds using a genotyping-by-sequencing approach to genotype 17,095 loci in 303 desert bighorn sheep. We found a signal of population genetic structure among remnant Mojave Desert populations, even across geographically proximate mountain ranges. Further, we found evidence of a genetically distinct, potential relict herd from a previously hypothesized Great Basin lineage of desert bighorn sheep. The genetic structure of source herds was clearly reflected in translocated populations. In most cases, herds retained genetic evidence of multiple translocation events and subsequent admixture when founded from multiple remnant source herds. Our results add to a growing literature on how population genomic data can be used to guide and monitor restoration programs.

Dosage Compensation and Gene Expression of the X Chromosome in Sheep

Ohno's hypothesis predicts that the expression of the single X chromosome in males needs compensatory upregulation to balance its dosage with that of the diploid autosomes. Additionally, X chromosome inactivation ensures that quadruple expression of the two X chromosomes is avoided in females. These mechanisms have been actively studied in mice and humans but lag behind in domestic species. Using RNA sequencing data, we analyzed the X chromosome upregulation in sheep fetal tissues from day 135 of gestation under control, over or restricted maternal diets (100%, 140% and 60% of National Research Council Total Digestible Nutrients), and in conceptuses, juvenile, and adult somatic tissues. By computing the mean expression ratio of all X-linked genes to all autosomal genes (X:A), we found that all samples displayed some levels of X chromosome upregulation. The degrees of X upregulation were not significant (P-value = 0.74) between ovine females and males in the same somatic tissues. Brain, however, displayed complete X upregulation. Interestingly, the male and female reproduction-related tissues exhibited divergent X dosage upregulation. Moreover, expression upregulation of the X chromosome in fetal tissues was not affected by maternal diets. Maternal nutrition, however, did change expression levels of several X-linked genes, such as sex determination genes SOX3 and NR0B1 In summary, our results showed that X chromosome upregulation occurred in nearly all sheep somatic tissues analyzed, thus support Ohno's hypothesis in a new species. However, the levels of upregulation differed by different subgroups of genes such as those that are house-keeping and "dosage-sensitive".

Genome-Wide Variation, Candidate Regions and Genes Associated With Fat Deposition and Tail Morphology in Ethiopian Indigenous Sheep

Variations in body weight and in the distribution of body fat are associated with feed availability, thermoregulation, and energy reserve. Ethiopia is characterized by distinct agro-ecological and human ethnic farmer diversity of ancient origin, which have impacted on the variation of its indigenous livestock. Here, we investigate autosomal genome-wide profiles of 11 Ethiopian indigenous sheep populations using the Illumina Ovine 50 K SNP BeadChip assay. Sheep from the Caribbean, Europe, Middle East, China, and western, northern and southern Africa were included to address globally, the genetic variation and history of Ethiopian populations. Population relationship and structure analysis separated Ethiopian indigenous fat-tail sheep from their North African and Middle Eastern counterparts. It indicates two main genetic backgrounds and supports two distinct genetic histories for African fat-tail sheep. Within Ethiopian sheep, our results show that the short fat-tail sheep do not represent a monophyletic group. Four genetic backgrounds are present in Ethiopian indigenous sheep but at different proportions among the fat-rump and the long fat-tail sheep from western and southern Ethiopia. The Ethiopian fat-rump sheep share a genetic background with Sudanese thin-tail sheep. Genome-wide selection signature analysis identified eight putative candidate regions spanning genes influencing growth traits and fat deposition (NPR2, HINT2, SPAG8, INSR), development of limbs and skeleton, and tail formation (ALX4, HOXB13, BMP4), embryonic development of tendons, bones and cartilages (EYA2, SULF2), regulation of body temperature (TRPM8), body weight and height variation (DIS3L2), control of lipogenesis and intracellular transport of long-chain fatty acids (FABP3), the occurrence and morphology of horns (RXFP2), and response to heat stress (DNAJC18). Our findings suggest that Ethiopian fat-tail sheep represent a uniquely admixed but distinct genepool that presents an important resource for understanding the genetic control of skeletal growth, fat metabolism and associated physiological processes.

Whole-genome sequences of 89 Chinese sheep suggest role of RXFP2 in the development of unique horn phenotype as response to semi-feralization

Background

Animal domestication has been extensively studied, but the process of feralization remains poorly understood.

Results

Here, we performed whole-genome sequencing of 99 sheep and identified a primary genetic divergence between 2 heterogeneous populations in the Tibetan Plateau, including 1 semi-feral lineage. Selective sweep and candidate gene analysis revealed local adaptations of these sheep associated with sensory perception, muscle strength, eating habit, mating process, and aggressive behavior. In particular, a horn-related gene, RXFP2, showed signs of rapid evolution specifically in the semi-feral breeds. A unique haplotype and repressed horn-related tissue expression of RXFP2 were correlated with higher horn length, as well as spiral and horizontally extended horn shape.

Conclusions

Semi-feralization has an extensive impact on diverse phenotypic traits of sheep. By acquiring features like those of their wild ancestors, semi-feral sheep were able to regain fitness while in frequent contact with wild surroundings and rare human interventions. This study provides a new insight into the evolution of domestic animals when human interventions are no longer dominant.

Whole-genome resequencing reveals signatures of selection and timing of duck domestication

Background

The genetic basis of animal domestication remains poorly understood, and systems with substantial phenotypic differences between wild and domestic populations are useful for elucidating the genetic basis of adaptation to new environments as well as the genetic basis of rapid phenotypic change. Here, we sequenced the whole genome of 78 individual ducks, from two wild and seven domesticated populations, with an average sequencing depth of 6.42X per individual.

Results

Our population and demographic analyses indicate a complex history of domestication, with early selection for separate meat and egg lineages. Genomic comparison of wild to domesticated populations suggests that genes that affect brain and neuronal development have undergone strong positive selection during domestication. Our F_ST analysis also indicates that the duck white plumage is the result of selection at the melanogenesis-associated transcription factor locus.

Conclusions

Our results advance the understanding of animal domestication and selection for complex phenotypic traits.

Improving pregnancy outcomes in humans through studies in sheep

Experimental studies that are relevant to human pregnancy rely on the selection of appropriate animal models as an important element in experimental design. Consideration of the strengths and weaknesses of any animal model of human disease is fundamental to effective and meaningful translation of preclinical research. Studies in sheep have made significant contributions to our understanding of the normal and abnormal development of the fetus. As a model of human pregnancy, studies in sheep have enabled scientists and clinicians to answer questions about the etiology and treatment of poor maternal, placental, and fetal health and to provide an evidence base for translation of interventions to the clinic. The aim of this review is to highlight the advances in perinatal human medicine that have been achieved following translation of research using the pregnant sheep and fetus.

Contrasting Patterns of Genomic Diversity Reveal Accelerated Genetic Drift but Reduced Directional Selection on X-Chromosome in Wild and Domestic Sheep Species

Analyses of genomic diversity along the X chromosome and of its correlation with autosomal diversity can facilitate understanding of evolutionary forces in shaping sex-linked genomic architecture. Strong selective sweeps and accelerated genetic drift on the X-chromosome have been inferred in primates and other model species, but no such insight has yet been gained in domestic animals compared with their wild relatives. Here, we analyzed X-chromosome variability in a large ovine data set, including a BeadChip array for 943 ewes from the world’s sheep populations and 110 whole genomes of wild and domestic sheep. Analyzing whole-genome sequences, we observed a substantially reduced X-to-autosome diversity ratio (∼0.6) compared with the value expected under a neutral model (0.75). In particular, one large X-linked segment (43.05–79.25 Mb) was found to show extremely low diversity, most likely due to a high density of coding genes, featuring highly conserved regions. In general, we observed higher nucleotide diversity on the autosomes, but a flat diversity gradient in X-linked segments, as a function of increasing distance from the nearest genes, leading to a decreased X: autosome (X/A) diversity ratio and contrasting to the positive correlation detected in primates and other model animals. Our evidence suggests that accelerated genetic drift but reduced directional selection on X chromosome, as well as sex-biased demographic events, explain low X-chromosome diversity in sheep species. The distinct patterns of X-linked and X/A diversity we observed between Middle Eastern and non-Middle Eastern sheep populations can be explained by multiple migrations, selection, and admixture during the domestic sheep’s recent postdomestication demographic expansion, coupled with natural selection for adaptation to new environments. In addition, we identify important novel genes involved in abnormal behavioral phenotypes, metabolism, and immunity, under selection on the sheep X-chromosome.

Functional Annotation of Animal Genomes (FAANG): Current Achievements and Roadmap

Functional annotation of genomes is a prerequisite for contemporary basic and applied genomic research, yet farmed animal genomics is deficient in such annotation. To address this, the FAANG (Functional Annotation of Animal Genomes) Consortium is producing genome-wide data sets on RNA expression, DNA methylation, and chromatin modification, as well as chromatin accessibility and interactions. In addition to informing our understanding of genome function, including comparative approaches to elucidate constrained sequence or epigenetic elements, these annotation maps will improve the precision and sensitivity of genomic selection for animal improvement. A scientific community-driven effort has already created a coordinated data collection and analysis enterprise crucial for the success of this global effort. Although it is early in this continuing process, functional data have already been produced and application to genetic improvement reported. The functional annotation delivered by the FAANG initiative will add value and utility to the greatly improved genome sequences being established for domesticated animal species.

Poor feed efficiency in sheep is associated with several structural abnormalities in the community metabolic network of their ruminal microbes

Ruminant animals have a symbiotic relationship with the microorganisms in their rumens. In this relationship, rumen microbes efficiently degrade complex plant-derived compounds into smaller digestible compounds, a process that is very likely associated with host animal feed efficiency. The resulting simpler metabolites can then be absorbed by the host and converted into other compounds by host enzymes. We used a microbial community metabolic network inferred from shotgun metagenomics data to assess how this metabolic system differs between animals that are able to turn ingested feedstuffs into body mass with high efficiency and those that are not. We conducted shotgun sequencing of microbial DNA from the rumen contents of 16 sheep that differed in their residual feed intake (RFI), a measure of feed efficiency. Metagenomic reads from each sheep were mapped onto a database-derived microbial metabolic network, which was linked to the sheep metabolic network by interface metabolites (metabolites transferred from microbes to host). No single enzyme was identified as being significantly different in abundance between the low and high RFI animals (P > 0.05, Wilcoxon test). However, when we analyzed the metabolic network as a whole, we found several differences between efficient and inefficient animals. Microbes from low RFI (efficient) animals use a suite of enzymes closer in network space to the host's reactions than those of the high RFI (inefficient) animals. Similarly, low RFI animals have microbial metabolic networks that, on average, contain reactions using shorter carbon chains than do those of high RFI animals, potentially allowing the host animals to extract metabolites more efficiently. Finally, the efficient animals possess community networks with greater Shannon diversity among their enzymes than do inefficient ones. Thus, our system approach to the ruminal microbiome identified differences attributable to feed efficiency in the structure of the microbes' community metabolic network that were undetected at the level of individual microbial taxa or reactions.

Draft genome of the Marco Polo Sheep (Ovis ammon polii)

Background

The Marco Polo Sheep (Ovis ammon polii), a subspecies of argali (Ovis ammon) that is distributed mainly in the Pamir Mountains, provides a mammalian model to study high altitude adaptation mechanisms. Due to over-hunting and subsistence poaching, as well as competition with livestock and habitat loss, O. ammon has been categorized as an endangered species on several lists. It can have fertile offspring with sheep. Hence, a high-quality reference genome of the Marco Polo Sheep will be very helpful in conservation genetics and even in exploiting useful genes in sheep breeding.

Findings

A total of 1022.43 Gb of raw reads resulting from whole-genome sequencing of a Marco Polo Sheep were generated using an Illumina HiSeq2000 platform. The final genome assembly (2.71 Gb) has an N50 contig size of 30.7 Kb and a scaffold N50 of 5.49 Mb. The repeat sequences identified account for 46.72% of the genome, and 20 336 protein-coding genes were predicted from the masked genome. Phylogenetic analysis indicated a close relationship between Marco Polo Sheep and the domesticated sheep, and the time of their divergence was approximately 2.36 million years ago. We identified 271 expanded gene families and 168 putative positively selected genes in the Marco Polo Sheep lineage.

Conclusions

We provide the first genome sequence and gene annotation for the Marco Polo Sheep. The availability of these resources will be of value in the future conservation of this endangered large mammal, for research into high altitude adaptation mechanisms, for reconstructing the evolutionary history of the Caprinae, and for the future conservation of the Marco Polo Sheep.

A 5-Methylcytosine Site of Growth Differentiation Factor 9 (GDF9) Gene Affects Its Tissue-Specific Expression in Sheep

Simple Summary

Growth differentiation factor 9 (GDF9) is an important gene for ovine fertility. GDF9 is highly expressed in the ovary as opposed to other tissues, but the reason for this is unknown. Our study found this can be caused by the methylation level of the promoter CpG island mC-4 site. This finding contributes to the understanding of the regulatory mechanism of GDF9 gene in reproduction.

Abstract

Growth differentiation factor 9 (GDF9) plays an important role in the early folliculogenesis of sheep. This study investigated the mRNA expression of ovine GDF9 in different tissues by real-time PCR. GDF9 exhibits significantly higher levels of expression (p < 0.01) in the ovary, relative to other tissues, indicating that its expression is tissue specific. To explore the regulatory mechanism of this tissue-specific expression, the methylation level of one CpG island (−1453 to −1854) of GDF9 promoter in ovary and heart was determined. In this region (−1987 to −1750), only the mC-4 site was present in the Sp4 binding site showed differential methylation between the heart and ovary; with increased (p < 0.01) methylation being observed in the heart. Additionally, the methylation level was negatively correlated with GDF9 mRNA expression (R = −0.75, p = 0.012), indicating that the methylation of this site plays an important role in transcriptional regulation of GDF9. The methylation effect of the mC-4 site was confirmed by using dual-luciferase. Site-directed mutation (methylation) of mC-4 site significantly reduced (p < 0.05) basal transcriptional activity of GDF9 promoter in oocytes. These results imply that methylation of GDF9 promoter CpG island mC-4 site may affect the binding of the Sp4 transcription factor to the GDF9 promoter region in sheep, thereby regulating GDF9 expression and resulting in a tissue-specific expression.

Chromosome-level assembly reveals extensive rearrangement in saker falcon and budgerigar, but not ostrich, genomes

BACKGROUND

The number of de novo genome sequence assemblies is increasing exponentially; however, relatively few contain one scaffold/contig per chromosome. Such assemblies are essential for studies of genotype-to-phenotype association, gross genomic evolution, and speciation. Inter-species differences can arise from chromosomal changes fixed during evolution, and we previously hypothesized that a higher fraction of elements under negative selection contributed to avian-specific phenotypes and avian genome organization stability. The objective of this study is to generate chromosome-level assemblies of three avian species (saker falcon, budgerigar, and ostrich) previously reported as karyotypically rearranged compared to most birds. We also test the hypothesis that the density of conserved non-coding elements is associated with the positions of evolutionary breakpoint regions.

RESULTS

We used reference-assisted chromosome assembly, PCR, and lab-based molecular approaches, to generate chromosome-level assemblies of the three species. We mapped inter- and intrachromosomal changes from the avian ancestor, finding no interchromosomal rearrangements in the ostrich genome, despite it being previously described as chromosomally rearranged. We found that the average density of conserved non-coding elements in evolutionary breakpoint regions is significantly reduced. Fission evolutionary breakpoint regions have the lowest conserved non-coding element density, and intrachromomosomal evolutionary breakpoint regions have the highest.

CONCLUSIONS

The tools used here can generate inexpensive, efficient chromosome-level assemblies, with > 80% assigned to chromosomes, which is comparable to genomes assembled using high-density physical or genetic mapping. Moreover, conserved non-coding elements are important factors in defining where rearrangements, especially interchromosomal, are fixed during evolution without deleterious effects.

Phylogenetic analysis of the cytochrome P450 (CYP450) nucleotide sequences of the horse and predicted CYP450s of the white rhinoceros (Ceratotherium simum) and other mammalian species

Background

The plight of the white rhinoceros (Ceratotherium simum) and the increasing need of treatment options for injured poaching victims led to the necessity to expand the knowledge on applicable drugs in this endangered species. With very little information available on drug pharmacokinetics in rhino, veterinarians have to rely on information generated from other species. The horse being a closely related species, has served as the model for dose extrapolations. However, from recent research on enrofloxacin and carprofen, the white rhino showed considerable differences in the pharmacokinetic properties of these drugs in comparison to the horse. While the reason for the differences is unknown, a likely cause may be a difference in present cytochrome P450 (CYP450), which may result in the rhino being genetically deficient in certain enzyme families.

Methods

For this paper we assess the degree of similarity of the CYP genome sequences across the different species, using BLAT (BLAST-like alignment tool) for the alignment of the nucleotide sequences of the equine CYP450 with potential homologous nucleotide sequences of the published database from white rhinos and other mammalian species (cow, pig, dog, sheep, elephant, mouse and human).

Results

The white rhino nucleotide sequences were 90.74% identical to the equine sequences. This was higher than the degree of similarity between any of the other evaluated species sequences. While no specific CYP family were found to be deficient in the published rhino genome, the horse genome contained additional genetic sequence for a larger number of iso-enzymes that were not present in the rhino.

Discussion

In pharmacokinetic study, it is well known that absence of a metabolic enzyme will result in constraints in drug metabolism and drug elimination. While this was our speculation, comparison to the horse and other mammalian species indicate that all the described CYP genes required for metabolism are present within the rhino genome. These results leave functional differences in enzyme activity and a lack of isoenzymes as the likely reason for the constraint in drug metabolism. Despite a more than 90% similarity of the equine and rhino gene sequences, seemingly small differences can have major effects on drug metabolism. Thus, in spite of the close anatomical relationship, the rhino should not simply be treated like a big horse.

Longitudinal In Vivo Monitoring of the CNS Demonstrates the Efficacy of Gene Therapy in a Sheep Model of CLN5 Batten Disease

Neuronal ceroid lipofuscinoses (NCLs; Batten disease) are neurodegenerative lysosomal storage diseases predominantly affecting children. Single administration of brain-directed lentiviral or recombinant single-stranded adeno-associated virus 9 (ssAAV9) vectors expressing ovine CLN5 into six pre-clinically affected sheep with a naturally occurring CLN5 NCL resulted in long-term disease attenuation. Treatment efficacy was demonstrated by non-invasive longitudinal in vivo monitoring developed to align with assessments used in human medicine. The treated sheep retained neurological and cognitive function, and one ssAAV9-treated animal has been retained and is now 57 months old, almost triple the lifespan of untreated CLN5-affected sheep. The onset of visual deficits was much delayed. Computed tomography and MRI showed that brain structures and volumes remained stable. Because gene therapy in humans is more likely to begin after clinical diagnosis, self-complementary AAV9-CLN5 was injected into the brain ventricles of four 7-month-old affected sheep already showing early clinical signs in a second trial. This also halted disease progression beyond their natural lifespan. These findings demonstrate the efficacy of CLN5 gene therapy, using three different vector platforms, in a large animal model and, thus, the prognosis for human translation.

How to estimate kinship

The concept of kinship permeates many domains of fundamental and applied biology ranging from social evolution to conservation science to quantitative and human genetics. Until recently, pedigrees were the gold standard to infer kinship, but the advent of next‐generation sequencing and the availability of dense genetic markers in many species make it a good time to (re)evaluate the usefulness of genetic markers in this context. Using three published data sets where both pedigrees and markers are available, we evaluate two common and a new genetic estimator of kinship. We show discrepancies between pedigree values and marker estimates of kinship and explore via simulations the possible reasons for these. We find these discrepancies are attributable to two main sources: pedigree errors and heterogeneity in the origin of founders. We also show that our new marker‐based kinship estimator has very good statistical properties and behaviour and is particularly well suited for situations where the source population is of small size, as will often be the case in conservation biology, and where high levels of kinship are expected, as is typical in social evolution studies.

Photoperiod and thyroid hormone regulate expression of l-dopachrome tautomerase (Dct), a melanocyte stem-cell marker, in tanycytes of the ovine hypothalamus

The pars tuberalis (PT) of the pituitary is central to the control of seasonal breeding. In mammals, the PT translates the photoperiodic message carried by melatonin into an endocrine thyroid-stimulating hormone output, which controls local thyroid hormone (TH) signalling in tanycytes of the neighbouring hypothalamus. In the present study, we identify l-dopachrome tautomerase (Dct) as a novel marker of ovine tanycytes and show that Dct displays marked seasonal variations in expression, with higher levels during spring and summer. This seasonal profile is photoperiod-dependent because an acute exposure to long days induces Dct expression. In addition, we find that TH also modulates Dct expression. DCT functions as an enzyme in the melanin synthesis pathway within skin melanocytes, whereas expression in other tissues is comparatively low. We demonstrate that both Tyr and Tyrp1, which are enzymes that intervene upstream and downstream of Dct in the melanin synthesis pathway, respectively, are expressed at very low levels in the ovine hypothalamus. This suggests that Dct in tanycytes may not be involved in melanin synthesis. We speculate that DCT function is linked to its protective role towards oxidative stress and/or its function in the control of neural progenitor cell proliferation.

Across-Experiment Transcriptomics of Sheep Rumen Identifies Expression of Lipid/Oxo-Acid Metabolism and Muscle Cell Junction Genes Associated With Variation in Methane-Related Phenotypes

Ruminants are significant contributors to the livestock generated component of the greenhouse gas, methane (CH₄). The CH₄ is primarily produced by the rumen microbes. Although the composition of the diet and animal intake amount have the largest effect on CH₄ production and yield (CH₄ production/dry matter intake, DMI), the host also influences CH₄ yield. Shorter rumen feed mean retention time (MRT) is associated with higher dry matter intake and lower CH₄ yield, but the molecular mechanism(s) by which the host affects CH₄ production remain unclear. We integrated rumen wall transcriptome data and CH₄ phenotypes from two independent experiments conducted with sheep in Australia (AUS, n = 62) and New Zealand (NZ, n = 24). The inclusion of the AUS data validated the previously identified clusters and gene sets representing rumen epithelial, metabolic and muscular functions. In addition, the expression of the cell cycle genes as a group was consistently positively correlated with acetate and butyrate concentrations (p < 0.05, based on AUS and NZ data together). The expression of a group of metabolic genes showed positive correlations in both AUS and NZ datasets with CH₄ production (p < 0.05) and yield (p < 0.01). These genes encode key enzymes in the ketone body synthesis pathway and included members of the poorly characterized aldo-keto reductase 1C (AKR1C) family. Several AKR1C family genes appear to have ruminant specific evolution patterns, supporting their specialized roles in the ruminants. Combining differential gene expression in the rumen wall muscle of the shortest and longest MRT AUS animals (no data available for the NZ animals) with correlation and network analysis, we identified a set of rumen muscle genes involved in cell junctions as potential regulators of MRT, presumably by influencing contraction rates of the smooth muscle component of the rumen wall. Higher rumen expression of these genes, including SYNPO (synaptopodin, p < 0.01) and NEXN (nexilin, p < 0.05), was associated with lower CH₄ yield in both AUS and NZ datasets. Unlike the metabolic genes, the variations in the expression of which may reflect the availability of rumen metabolites, the muscle genes are currently our best candidates for causal genes that influence CH₄ yield.

Global genomic diversity and conservation priorities for domestic animals are associated with the economies of their regions of origin

Domestic animals play a key role in human survival and the development of civilization. However, the genetic resources of domestic animals are facing an alarming rate of erosion due to socioeconomic changes, economic globalization and financial constraints. In this study, through genome-wide SNP analysis, we estimated the heterozygosity, inbreeding coefficient, effective population size, and runs of homozygosity to identify the breeds facing the risk of extinction for sheep and cattle across the world. In particular, we quantified the contribution of 97 sheep breeds and 53 cattle breeds to genomic diversity (within-breed, between-breed and total) and prioritized the breeds for conservation. Additionally, we compared the average values of genomic diversity between breeds from regions (or countries) in different economic categories (underdeveloped, developing and developed), and found that breeds in developed regions exhibit significantly higher levels of total genomic diversity than those in underdeveloped and developing regions. Altogether, our results suggested that conservation priority should be given to breeds in developed regions to secure the future genomic diversity hotspots of domestic animal resources.

Comparative transcriptomics of 5 high-altitude vertebrates and their low-altitude relatives

Background

Species living at high altitude are subject to strong selective pressures due to inhospitable environments (e.g., hypoxia, low temperature, high solar radiation, and lack of biological production), making these species valuable models for comparative analyses of local adaptation. Studies that have examined high-altitude adaptation have identified a vast array of rapidly evolving genes that characterize the dramatic phenotypic changes in high-altitude animals. However, how high-altitude environment shapes gene expression programs remains largely unknown.

Findings

We generated a total of 910 Gb of high-quality RNA-seq data for 180 samples derived from 6 tissues of 5 agriculturally important high-altitude vertebrates (Tibetan chicken, Tibetan pig, Tibetan sheep, Tibetan goat, and yak) and their cross-fertile relatives living in geographically neighboring low-altitude regions. Of these, ∼75% reads could be aligned to their respective reference genomes, and on average ∼60% of annotated protein coding genes in each organism showed FPKM expression values greater than 0.5. We observed a general concordance in topological relationships between the nucleotide alignments and gene expression-based trees. Tissue and species accounted for markedly more variance than altitude based on either the expression or the alternative splicing patterns. Cross-species clustering analyses showed a tissue-dominated pattern of gene expression and a species-dominated pattern for alternative splicing. We also identified numerous differentially expressed genes that could potentially be involved in phenotypic divergence shaped by high-altitude adaptation.

Conclusions

These data serve as a valuable resource for examining the convergence and divergence of gene expression changes between species as they adapt or acclimatize to high-altitude environments.

Whole-genome resequencing of Ujumqin sheep to investigate the determinants of the multi-vertebral trait

The Ujumqin sheep is one of the most profitable breeds in China, with unique multi-vertebral characteristics. We performed high-throughput genome resequencing of five multi-vertebral and three non-multi-vertebral sheep in an Ujumqin population. We identified the genomic regions that correlated with the germplasm characteristics to establish the cause of the "multi-vertebral" phenotype in this breed. Sequencing generated a total of 314 952 000 000 bp of raw data. The alignment rate of all the samples was between 98.53% and 99.11%, and the mean depth of coverage relative to the reference genome was between 11.58× and 14.92×. After comparing the differences between the two groups, we identified 21 homozygous single nucleotide polymorphisms (SNPs) in the mutant exons of 14 genes. Nineteen loci of 10 genes contained nonsynonymous mutations, while two loci contained synonymous mutations. Resequencing revealed homozygous mutations comprised of 44 indels located within exons of 19 genes. These indels included 37 frameshift mutations, 6 non-frameshift mutations, and 1 stopgain single nucleotide variation (SNV). Finally, comparisons of genotypic variations revealed 17 genes with homozygous mutations in their coding regions, 5 of which have previously been associated with vertebral development and the remaining 12 genes were newly identified in this study.