Mary E. B. McCulloch

A high resolution atlas of gene expression in the domestic sheep (Ovis aries)

Sheep are a key source of meat, milk and fibre for the global livestock sector, and an important biomedical model. Global analysis of gene expression across multiple tissues has aided genome annotation and supported functional annotation of mammalian genes. We present a large-scale RNA-Seq dataset representing all the major organ systems from adult sheep and from several juvenile, neonatal and prenatal developmental time points. The Ovis aries reference genome (Oar v3.1) includes 27,504 genes (20,921 protein coding), of which 25,350 (19,921 protein coding) had detectable expression in at least one tissue in the sheep gene expression atlas dataset. Network-based cluster analysis of this dataset grouped genes according to their expression pattern. The principle of ‘guilt by association’ was used to infer the function of uncharacterised genes from their co-expression with genes of known function. We describe the overall transcriptional signatures present in the sheep gene expression atlas and assign those signatures, where possible, to specific cell populations or pathways. The findings are related to innate immunity by focusing on clusters with an immune signature, and to the advantages of cross-breeding by examining the patterns of genes exhibiting the greatest expression differences between purebred and crossbred animals. This high-resolution gene expression atlas for sheep is, to our knowledge, the largest transcriptomic dataset from any livestock species to date. It provides a resource to improve the annotation of the current reference genome for sheep, presenting a model transcriptome for ruminants and insight into gene, cell and tissue function at multiple developmental stages.

Macrophage colony-stimulating factor (CSF1) controls monocyte production and maturation and the steady-state size of the liver in pigs

This study is based on extensive studies in the mouse of the role of CSF1 in monocyte-macrophage production and differentiation and the function of macrophages in the control of hepatocyte proliferation. We use a novel form of CSF1, an Fc fusion protein, to demonstrate that the findings in mice can be extended to large animals. We discuss the possible role for CSF1 in homeostatic control of the size of the liver.

Integration of quantitated expression estimates from polyA-selected and rRNA-depleted RNA-seq libraries

BackgroundThe availability of fast alignment-free algorithms has greatly reduced the computational burden of RNA-seq processing, especially for relatively poorly assembled genomes. Using these approaches, previous RNA-seq datasets could potentially be processed and integrated with newly sequenced libraries. Confounding factors in such integration include sequencing depth and methods of RNA extraction and selection. Different selection methods (typically, either polyA-selection or rRNA-depletion) omit different RNAs, resulting in different fractions of the transcriptome being sequenced. In particular, rRNA-depleted libraries sample a broader fraction of the transcriptome than polyA-selected libraries. This study aimed to develop a systematic means of accounting for library type that allows data from these two methods to be compared.ResultsThe method was developed by comparing two RNA-seq datasets from ovine macrophages, identical except for RNA selection method. Gene-level expression estimates were obtained using a two-part process centred on the high-speed transcript quantification tool Kallisto. Firstly, a set of reference transcripts was defined that constitute a standardised RNA space, with expression from both datasets quantified against it. Secondly, a simple ratio-based correction was applied to the rRNA-depleted estimates. The outcome is an almost perfect correlation between gene expression estimates, independent of library type and across the full range of levels of expression.ConclusionA combination of reference transcriptome filtering and a ratio-based correction can create equivalent expression profiles from both polyA-selected and rRNA-depleted libraries. This approach will allow meta-analysis and integration of existing RNA-seq data into transcriptional atlas projects.

The genomic architecture of mastitis resistance in dairy sheep

BackgroundMastitis is the most prevalent disease in dairy sheep with major economic, hygienic and welfare implications. The disease persists in all dairy sheep production systems despite the implementation of improved management practises. Selective breeding for enhanced mastitis resistance may provide the means to further control the disease. In the present study, we investigated the genetic architecture of four mastitis traits in dairy sheep. Individual animal records for clinical mastitis occurrence and three mastitis indicator traits (milk somatic cell count, total viable bacterial count in milk and the California mastitis test) were collected monthly throughout lactation for 609 ewes of the Greek Chios breed. All animals were genotyped with a custom-made 960-single nucleotide polymorphism (SNP) DNA array based on markers located in quantitative trait loci (QTL) regions for mastitis resistance previously detected in three other distinct dairy sheep populations.ResultsHeritable variation and strong positive genetic correlations were estimated for clinical mastitis occurrence and the three mastitis indicator traits. SNP markers significantly associated with these mastitis traits were confirmed on chromosomes 2, 3, 5, 16 and 19. We identified pathways, molecular interaction networks and functional gene clusters for mastitis resistance. Candidate genes within the detected regions were identified based upon analysis of an ovine transcriptional atlas and transcriptome data derived from milk somatic cells. Relevant candidate genes implicated in innate immunity included SOCS2, CTLA4, C6, C7, C9, PTGER4, DAB2, CARD6, OSMR, PLXNC1, IDH1, ICOS, FYB, and LYFR.ConclusionsThe results confirmed the presence of animal genetic variability in mastitis resistance and identified genomic regions associated with specific mastitis traits in the Chios sheep. The conserved genetic architecture of mastitis resistance between distinct dairy sheep breeds suggests that across-breed selection programmes would be feasible.

Species-Specific Transcriptional Regulation of Genes Involved in Nitric Oxide Production and Arginine Metabolism in Macrophages

Activated mouse macrophages metabolize arginine via NO synthase (NOS2) to produce NO as an antimicrobial effector. Published gene expression datasets provide little support for the activation of this pathway in human macrophages. Generation of NO requires the coordinated regulation of multiple genes. We have generated RNA-sequencing data from bone marrow–derived macrophages from representative rodent (rat), monogastric (pig and horse), and ruminant (sheep, goat, cattle, and water buffalo) species, and analyzed the expression of genes involved in arginine metabolism in response to stimulation with LPS. In rats, as in mice, LPS strongly induced Nos2, the arginine transporter Slc7a2, arginase 1 (Arg1), GTP cyclohydrolase (Gch1), and argininosuccinate synthase (Ass1). None of these responses was conserved across species. Only cattle and water buffalo showed substantial NOS2 induction. The species studied also differed in expression and regulation of arginase (ARG2, rather than ARG1), and amino acid transporters. Variation between species was associated with rapid promoter evolution. Differential induction of NOS2 and ARG2 between the ruminant species was associated with insertions of the Bov-A2 retrotransposon in the promoter region. Bov-A2 was shown to possess LPS-inducible enhancer activity in transfected RAW264.7 macrophages. Consistent with a function in innate immunity, NO production and arginine metabolism vary greatly between species and differences may contribute to pathogen host restriction.

ADGRE1 (EMR1, F4/80) Is a Rapidly-Evolving Gene Expressed in Mammalian Monocyte-Macrophages

The F4/80 antigen, encoded by the Adgre1 locus, has been widely-used as a monocyte-macrophage marker in mice, but its value as a macrophage marker in other species is unclear, and has even been questioned. ADGRE1 is a seven transmembrane G protein-coupled receptor with an extracellular domain containing repeated Epidermal Growth Factor (EGF)-like calcium binding domains. Using a new monoclonal antibody, we demonstrated that ADGRE1 is a myeloid differentiation marker in pigs, absent from progenitors in bone marrow, highly-expressed in mature granulocytes, monocytes, and tissue macrophages and induced by macrophage colony-stimulating factor (CSF1) treatment in vivo. Based upon these observations, we utilized RNA-Seq to assess the expression of ADGRE1 mRNA in bone marrow or monocyte-derived macrophages (MDM) and alveolar macrophages from 8 mammalian species including pig, human, rat, sheep, goat, cow, water buffalo, and horse. ADGRE1 mRNA was expressed by macrophages in each species, with inter-species variation both in expression level and response to lipopolysaccharide (LPS) stimulation. Analysis of the RNA-Seq data also revealed additional exons in several species compared to current Ensembl annotations. The ruminant species and horses appear to encode a complete duplication of the 7 EGF-like domains. In every species, Sashimi plots revealed evidence of exon skipping of the EGF-like domains, which are highly-variable between species and polymorphic in humans. Consistent with these expression patterns, key elements of the promoter and a putative enhancer are also conserved across all species. The rapid evolution of this molecule and related ADGRE family members suggests immune selection and a role in pathogen recognition.