A Mini-Atlas of Gene Expression for the Domestic Goat (Capra hircus)

Goats (Capra hircus) are an economically important livestock species providing meat and milk across the globe. They are of particular importance in tropical agri-systems contributing to sustainable agriculture, alleviation of poverty, social cohesion, and utilisation of marginal grazing. There are excellent genetic and genomic resources available for goats, including a highly contiguous reference genome (ARS1). However, gene expression information is limited in comparison to other ruminants. To support functional annotation of the genome and comparative transcriptomics, we created a mini-atlas of gene expression for the domestic goat. RNA-Seq analysis of 17 transcriptionally rich tissues and 3 cell-types detected the majority (90%) of predicted protein-coding transcripts and assigned informative gene names to more than 1000 previously unannotated protein-coding genes in the current reference genome for goat (ARS1). Using network-based cluster analysis, we grouped genes according to their expression patterns and assigned those groups of coexpressed genes to specific cell populations or pathways. We describe clusters of genes expressed in the gastro-intestinal tract and provide the expression profiles across tissues of a subset of genes associated with functional traits. Comparative analysis of the goat atlas with the larger sheep gene expression atlas dataset revealed transcriptional similarities between macrophage associated signatures in the sheep and goats sampled in this study. The goat transcriptomic resource complements the large gene expression dataset we have generated for sheep and contributes to the available genomic resources for interpretation of the relationship between genotype and phenotype in small ruminants.

Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions

Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.

Developmental Stage-Specific Distribution of Macrophages in Mouse Mammary Gland

Mammary gland development begins in the embryo and continues throughout the reproductive life of female mammals. Tissue macrophages (Mϕs), dependent on signals from the Mϕ colony stimulating factor 1 receptor (CSF1R), have been shown to regulate the generation, regression and regeneration of this organ, which is central for mammalian offspring survival. However, the distribution of Mϕs in the pre- and post-natal mammary gland, as it undergoes distinct phases of development and regression, is unknown or has been inferred from immunostaining of thin tissue sections. Here, we used optical tissue clearing and 3-dimensional imaging of mammary tissue obtained from Csf1r-EGFP mice. Whilst tissue Mϕs were observed at all developmental phases, their abundance, morphology, localization and association with luminal and basal epithelial cells exhibited stage-specific differences. Furthermore, sexual dimorphism was observed at E14.5, when the male mammary bud is severed from the overlying epidermis. These findings provide new insights into the localization and possible functions of heterogeneous tissue Mϕ populations in mammogenesis.

Antigen Sampling CSF1R-Expressing Epithelial Cells Are the Functional Equivalents of Mammalian M Cells in the Avian Follicle-Associated Epithelium

The follicle-associated epithelium (FAE) is a specialized structure that samples luminal antigens and transports them into mucosa-associated lymphoid tissues (MALT). In mammals, transcytosis of antigens across the gut epithelium is performed by a subset of FAE cells known as M cells. Here we show that colony-stimulating factor 1 receptor (CSF1R) is expressed by a subset of cells in the avian bursa of Fabricius FAE. Expression was initially detected using a CSF1R-reporter transgene that also label subsets of bursal macrophages. Immunohistochemical detection using a specific monoclonal antibody confirmed abundant expression of CSF1R on the basolateral membrane of FAE cells. CSF1R-transgene expressing bursal FAE cells were enriched for expression of markers previously reported as putative M cell markers, including annexin A10 and CD44. They were further distinguished from a population of CSF1R-transgene negative epithelial cells within FAE by high apical F-actin expression and differential staining with the lectins jacalin, PHA-L and SNA. Bursal FAE cells that express the CSF1R-reporter transgene were responsible for the bulk of FAE transcytosis of labeled microparticles in the size range 0.02-0.1 μm. Unlike mammalian M cells, they did not readily take up larger bacterial sized microparticles (0.5 μm). Their role in uptake of bacteria was tested using Salmonella, which can enter via M cells in mammals. Labeled Salmonella enterica serovar Typhimurium entered bursal tissue via the FAE. Entry was partially dependent upon Type III secretion system-1. However, the majority of invading bacteria were localized to CSF1R-negative FAE cells and in resident phagocytes that express the phosphatidylserine receptor TIM4. CSF1R-expressing FAE cells in infected follicles showed evidence of cell death and shedding into the bursal lumen. In mammals, CSF1R expression in the gut is restricted to macrophages which only indirectly control M cell differentiation. The novel expression of CSF1R in birds suggests that these functional equivalents to mammalian M cells may have different ontological origins and their development and function are likely to be regulated by different growth factors.

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.

Assembly of a parts list of the human mitotic cell cycle machinery

The set of proteins required for mitotic division remains poorly characterized. Here, an extensive series of correlation analyses of human and mouse transcriptomics data were performed to identify genes strongly and reproducibly associated with cells undergoing S/G2-M phases of the cell cycle. In so doing, 701 cell cycle-associated genes were defined and while it was shown that many are only expressed during these phases, the expression of others is also driven by alternative promoters. Of this list, 496 genes have known cell cycle functions, whereas 205 were assigned as putative cell cycle genes, 53 of which are functionally uncharacterized. Among these, 27 were screened for subcellular localization revealing many to be nuclear localized and at least three to be novel centrosomal proteins. Furthermore, 10 others inhibited cell proliferation upon siRNA knockdown. This study presents the first comprehensive list of human cell cycle proteins, identifying many new candidate proteins.

A Gene Expression Atlas of the Domestic Water Buffalo (Bubalus bubalis)

The domestic water buffalo (Bubalus bubalis) makes a major contribution to the global agricultural economy in the form of milk, meat, hides, and draught power. The global water buffalo population is predominantly found in Asia, and per head of population more people depend upon the buffalo than on any other livestock species. Despite its agricultural importance, there are comparatively fewer genomic and transcriptomic resources available for buffalo than for other livestock species. We have generated a large-scale gene expression atlas covering multiple tissue and cell types from all major organ systems collected from three breeds of riverine water buffalo (Mediterranean, Pandharpuri and Bhadawari) and used the network analysis tool Graphia Professional to identify clusters of genes with similar expression profiles. Alongside similar data, we and others have generated for ruminants as part of the Functional Annotation of Animal Genomes Consortium; this comprehensive transcriptome supports functional annotation and comparative analysis of the water buffalo genome.

Phenotypic impacts of CSF1R deficiencies in humans and model organisms

Mϕ proliferation, differentiation, and survival are controlled by signals from the Mϕ CSF receptor (CSF1R). Mono-allelic gain-of-function mutations in CSF1R in humans are associated with an autosomal-dominant leukodystrophy and bi-allelic loss-of-function mutations with recessive skeletal dysplasia, brain disorders, and developmental anomalies. Most of the phenotypes observed in these human disease states are also observed in mice and rats with loss-of-function mutations in Csf1r or in Csf1 encoding one of its two ligands. Studies in rodent models also highlight the importance of genetic background and likely epistatic interactions between Csf1r and other loci. The impacts of Csf1r mutations on the brain are usually attributed solely to direct impacts on microglial number and function. However, analysis of hypomorphic Csf1r mutants in mice and several other lines of evidence suggest that primary hydrocephalus and loss of the physiological functions of Mϕs in the periphery contribute to the development of brain pathology. In this review, we outline the evidence that CSF1R is expressed exclusively in mononuclear phagocytes and explore the mechanisms linking CSF1R mutations to pleiotropic impacts on postnatal growth and development.

Deletion of a Csf1r enhancer selectively impacts CSF1R expression and development of tissue macrophage populations

The proliferation, differentiation and survival of mononuclear phagocytes depend on signals from the receptor for macrophage colony-stimulating factor, CSF1R. The mammalian Csf1r locus contains a highly conserved super-enhancer, the fms-intronic regulatory element (FIRE). Here we show that genomic deletion of FIRE in mice selectively impacts CSF1R expression and tissue macrophage development in specific tissues. Deletion of FIRE ablates macrophage development from murine embryonic stem cells. Csf1rΔFIRE/ΔFIRE mice lack macrophages in the embryo, brain microglia and resident macrophages in the skin, kidney, heart and peritoneum. The homeostasis of other macrophage populations and monocytes is unaffected, but monocytes and their progenitors in bone marrow lack surface CSF1R. Finally, Csf1rΔFIRE/ΔFIRE mice are healthy and fertile without the growth, neurological or developmental abnormalities reported in Csf1r-/- rodents. Csf1rΔFIRE/ΔFIRE mice thus provide a model to explore the homeostatic, physiological and immunological functions of tissue-specific macrophage populations in adult animals.

Linked Mutations at Adjacent Nucleotides Have Shaped Human Population Differentiation and Protein Evolution

Despite the fundamental importance of single nucleotide polymorphisms (SNPs) to human evolution, there are still large gaps in our understanding of the forces that shape their distribution across the genome. SNPs have been shown to not be distributed evenly, with directly adjacent SNPs found unusually frequently. Why this is the case is unclear. We illustrate how neighboring SNPs that cannot be explained by a single mutation event (that we term here sequential dinucleotide mutations [SDMs]) are driven by distinct processes to SNPs and multinucleotide polymorphisms (MNPs). By studying variation across populations, including a novel cohort of 1,358 Scottish genomes, we show that, SDMs are over twice as common as MNPs and like SNPs display distinct mutational spectra across populations. These biases are not only different to those observed among SNPs and MNPs but are also more divergent between human population groups. We show that the changes that make up SDMs are not independent and identify a distinct mutational profile, CA → CG → TG, that is observed an order of magnitude more often than expected from background SNP rates and the numbers of other SDMs involving the gain and deamination of CpG sites. Intriguingly particular pathways through the amino acid code appear to have been favored relative to that expected from intergenic SDM rates and the occurrences of coding SNPs, and in particular those that lead to the creation of single codon amino acids. We finally present evidence that epistatic selection has potentially disfavored sequential nonsynonymous changes in the human genome.

The Effect of Race Training on the Basal Gene Expression of Alveolar Macrophages Derived From Standardbred Racehorses

Mild-to-moderate equine asthma is prevalent in young racehorses, particularly early in their training period. Although the precise etiopathogenesis remains undetermined, it is possible that the susceptibility of this population might partly reflect an exercise-associated immune derangement at the level of the airway. We performed a genome-wide basal gene expression scan on alveolar macrophages (AMs) isolated from Standardbred racehorses before and after commencement of competition race training with a view to identifying any exercise-associated gene expression modulation consistent with functional alterations, which might reflect training-associated immunological derangement. Microarray technology was used to analyze the basal gene expression profiles of bronchoalveolar fluid-derived AMs, harvested from six systemically healthy Standardbred racehorses before (T0) and after (T1) entry into training. In addition, AM lipopolysaccharide (LPS)-induced TNF-α and IL-10 release at T0 and T1 was assessed. Although the data revealed significant interhorse heterogeneity in relation to the magnitude of individual gene expression at each timepoint, within each horse, several inflammatory-related genes [e.g., chemokine ligands, interferons, and nuclear factor kappa-light-chain-enhancer of activated B cells (NFKB)] declined in expression from T0 to T1. Entry into training did not significantly alter AM LPS-induced TNF-α or IL-10 release. The data support a direct effect of training on AM basal gene expression, particularly with respect to immune-related genes. The pattern of training-associated differential gene expression may indicate relative downregulation of inflammatory-related genes, consistent with an immunosuppressive effect of training and an increased susceptibility to opportunistic pathogens.

Measurement of serum macrophage migration inhibitory factor (MIF) and correlation with severity and pruritus scores in client owned dogs with atopic dermatitis

BACKGROUND

Atopic dermatitis (AD) is a common inflammatory skin disease of dogs. Macrophage migration inhibitory factor (MIF) initiates pro-inflammatory cytokine release in human AD and serum concentrations are correlated with disease severity.

HYPOTHESIS

Canine serum MIF concentrations are increased in dogs with AD and correlate with clinical lesion and pruritus scores.

ANIMALS

Client owned dogs (n = 49) diagnosed with AD and 17 healthy, unaffected control dogs were used for the study.

METHODS AND MATERIALS

A commercially available MIF ELISA was optimized for the dog and serum from clinical cases used. Information regarding treatment, Canine Atopic Dermatitis Extent and Severity Index, (CADESI-4) and pruritus Visual Analog Scale (pVAS) were recorded for each dog at the time of serum collection.

RESULTS

Dogs with AD which had not received steroid therapy and those treated with oclacitinib had significantly elevated serum MIF concentrations compared to controls. Concentrations of MIF were not significantly different in AD dogs receiving steroids compared to controls. There was no significant correlation between MIF concentrations and clinical scores (CADESI-4 or pVAS).

CONCLUSIONS AND CLINICAL IMPORTANCE

Serum MIF concentrations are increased in dogs with AD and MIF might be a target for therapy.

Characterization of Subpopulations of Chicken Mononuclear Phagocytes That Express TIM4 and CSF1R

The phosphatidylserine receptor TIM4, encoded by TIMD4, mediates the phagocytic uptake of apoptotic cells. We applied anti-chicken TIM4 mAbs in combination with CSF1R reporter transgenes to dissect the function of TIM4 in the chick (Gallus gallus). During development in ovo, TIM4 was present on the large majority of macrophages, but expression became more heterogeneous posthatch. Blood monocytes expressed KUL01, class II MHC, and CSF1R-mApple uniformly. Around 50% of monocytes were positive for surface TIM4. They also expressed many other monocyte-specific transcripts at a higher level than TIM4^- monocytes. In liver, highly phagocytic TIM4^hi cells shared many transcripts with mammalian Kupffer cells and were associated with uptake of apoptotic cells. Although they expressed CSF1R mRNA, Kupffer cells did not express the CSF1R-mApple transgene, suggesting that additional CSF1R transcriptional regulatory elements are required by these cells. By contrast, CSF1R-mApple was detected in liver TIM4^lo and TIM4^- cells, which were not phagocytic and were more abundant than Kupffer cells. These cells expressed CSF1R alongside high levels of FLT3, MHCII, XCR1, and other markers associated with conventional dendritic cells in mice. In bursa, TIM4 was present on the cell surface of two populations. Like Kupffer cells, bursal TIM4^hi phagocytes coexpressed many receptors involved in apoptotic cell recognition. TIM4^lo cells appear to be a subpopulation of bursal B cells. In overview, TIM4 is associated with phagocytes that eliminate apoptotic cells in the chick. In the liver, TIM4 and CSF1R reporters distinguished Kupffer cells from an abundant population of dendritic cell-like cells.

A chicken bioreactor for efficient production of functional cytokines

BACKGROUND

The global market for protein drugs has the highest compound annual growth rate of any pharmaceutical class but their availability, especially outside of the US market, is compromised by the high cost of manufacture and validation compared to traditional chemical drugs. Improvements in transgenic technologies allow valuable proteins to be produced by genetically-modified animals; several therapeutic proteins from such animal bioreactors are already on the market after successful clinical trials and regulatory approval. Chickens have lagged behind mammals in bioreactor development, despite a number of potential advantages, due to the historic difficulty in producing transgenic birds, but the production of therapeutic proteins in egg white of transgenic chickens would substantially lower costs across the entire production cycle compared to traditional cell culture-based production systems. This could lead to more affordable treatments and wider markets, including in developing countries and for animal health applications.

RESULTS

Here we report the efficient generation of new transgenic chicken lines to optimize protein production in eggs. As proof-of-concept, we describe the expression, purification and functional characterization of three pharmaceutical proteins, the human cytokine interferon α2a and two species-specific Fc fusions of the cytokine CSF1.

CONCLUSION

Our work optimizes and validates a transgenic chicken system for the cost-effective production of pure, high quality, biologically active protein for therapeutics and other applications.

The Mononuclear Phagocyte System: The Relationship between Monocytes and Macrophages

The mononuclear phagocyte system (MPS) is defined as a cell lineage in which committed marrow progenitors give rise to blood monocytes and tissue macrophages. Here, we discuss the concept of self-proscribed macrophage territories and homeostatic regulation of tissue macrophage abundance through growth factor availability. Recent studies have questioned the validity of the MPS model and argued that tissue-resident macrophages are a separate lineage seeded during development and maintained by self-renewal. We address this issue; discuss the limitations of inbred mouse models of monocyte-macrophage homeostasis; and summarize the evidence suggesting that during postnatal life, monocytes can replace resident macrophages in all major organs and adopt their tissue-specific gene expression. We conclude that the MPS remains a valid and accurate framework for understanding macrophage development and homeostasis.

Dissecting the Genomic Architecture of Resistance to Eimeria maxima Parasitism in the Chicken

Coccidiosis in poultry, caused by protozoan parasites of the genus Eimeria, is an intestinal disease with substantial economic impact. With the use of anticoccidial drugs under public and political pressure, and the comparatively higher cost of live-attenuated vaccines, an attractive complementary strategy for control is to breed chickens with increased resistance to Eimeria parasitism. Prior infection with Eimeria maxima leads to complete immunity against challenge with homologous strains, but only partial resistance to challenge with antigenically diverse heterologous strains. We investigate the genetic architecture of avian resistance to E. maxima primary infection and heterologous strain secondary challenge using White Leghorn populations of derived inbred lines, C.B12 and 15I, known to differ in susceptibility to the parasite. An intercross population was infected with E. maxima Houghton (H) strain, followed 3 weeks later by E. maxima Weybridge (W) strain challenge, while a backcross population received a single E. maxima W infection. The phenotypes measured were parasite replication (counting fecal oocyst output or qPCR for parasite numbers in intestinal tissue), intestinal lesion score (gross pathology, scale 0-4), and for the backcross only, serum interleukin-10 (IL-10) levels. Birds were genotyped using a high density genome-wide DNA array (600K, Affymetrix). Genome-wide association study located associations on chromosomes 1, 2, 3, and 5 following primary infection in the backcross population, and a suggestive association on chromosome 1 following heterologous E. maxima W challenge in the intercross population. This mapped several megabases away from the quantitative trait locus (QTL) linked to the backcross primary W strain infection, suggesting different underlying mechanisms for the primary- and heterologous secondary- responses. Underlying pathways for those genes located in the respective QTL for resistance to primary infection and protection against heterologous challenge were related mainly to immune response, with IL-10 signaling in the backcross primary infection being the most significant. Additionally, the identified markers associated with IL-10 levels exhibited significant additive genetic variance. We suggest this is a phenotype of interest to the outcome of challenge, being scalable in live birds and negating the requirement for single-bird cages, fecal oocyst counts, or slaughter for sampling (qPCR).

Phenotypic and genetic variation in the response of chickens to Eimeria tenella induced coccidiosis

BACKGROUND

Coccidiosis is a major contributor to losses in poultry production. With emerging constraints on the use of in-feed prophylactic anticoccidial drugs and the relatively high costs of effective vaccines, there are commercial incentives to breed chickens with greater resistance to this important production disease. To identify phenotypic biomarkers that are associated with the production impacts of coccidiosis, and to assess their covariance and heritability, 942 Cobb500 commercial broilers were subjected to a defined challenge with Eimeria tenella (Houghton). Three traits were measured: weight gain (WG) during the period of infection, caecal lesion score (CLS) post mortem, and the level of a serum biomarker of intestinal inflammation, i.e. circulating interleukin 10 (IL-10), measured at the height of the infection.

RESULTS

Phenotypic analysis of the challenged chicken cohort revealed a significant positive correlation between CLS and IL-10, with significant negative correlations of both these traits with WG. Eigenanalysis of phenotypic covariances between measured traits revealed three distinct eigenvectors. Trait weightings of the first eigenvector, (EV1, eigenvalue = 59%), were biologically interpreted as representing a response of birds that were susceptible to infection, with low WG, high CLS and high IL-10. Similarly, the second eigenvector represented infection resilience/resistance (EV2, 22%; high WG, low CLS and high IL-10), and the third eigenvector tolerance (EV3, 19%; high WG, high CLS and low IL-10), respectively. Genome-wide association studies (GWAS) identified two SNPs that were associated with WG at the suggestive level.

CONCLUSIONS

Eigenanalysis separated the phenotypic impact of a defined challenge with E. tenella on WG, caecal inflammation/pathology, and production of IL-10 into three major eigenvectors, indicating that the susceptibility-resistance axis is not a single continuous quantitative trait. The SNPs identified by the GWAS for body weight were located in close proximity to two genes that are involved in innate immunity (FAM96B and RRAD).