The relationship between extreme inter-individual variation in macrophage gene expression and genetic susceptibility to inflammatory bowel disease

The differentiation of resident intestinal macrophages from blood monocytes depends upon signals from the macrophage colony-stimulating factor receptor (CSF1R). Analysis of genome-wide association studies (GWAS) indicates that dysregulation of macrophage differentiation and response to microorganisms contributes to susceptibility to chronic inflammatory bowel disease (IBD). Here, we analyzed transcriptomic variation in monocyte-derived macrophages (MDM) from affected and unaffected sib pairs/trios from 22 IBD families and 6 healthy controls. Transcriptional network analysis of the data revealed no overall or inter-sib distinction between affected and unaffected individuals in basal gene expression or the temporal response to lipopolysaccharide (LPS). However, the basal or LPS-inducible expression of individual genes varied independently by as much as 100-fold between subjects. Extreme independent variation in the expression of pairs of HLA-associated transcripts (HLA-B/C, HLA-A/F and HLA-DRB1/DRB5) in macrophages was associated with HLA genotype. Correlation analysis indicated the downstream impacts of variation in the immediate early response to LPS. For example, variation in early expression of IL1B was significantly associated with local SNV genotype and with subsequent peak expression of target genes including IL23A, CXCL1, CXCL3, CXCL8 and NLRP3. Similarly, variation in early IFNB1 expression was correlated with subsequent expression of IFN target genes. Our results support the view that gene-specific dysregulation in macrophage adaptation to the intestinal milieu is associated with genetic susceptibility to IBD.

Reversible expansion of tissue macrophages in response to macrophage colony-stimulating factor (CSF1) transforms systemic lipid and carbohydrate metabolism

Macrophages regulate metabolic homeostasis in health and disease. Macrophage colony-stimulating factor (CSF1)-dependent macrophages contribute to homeostatic control of the size of the liver. This study aimed to determine the systemic metabolic consequences of elevating circulating CSF1. Acute administration of a CSF1-Fc fusion protein to mice led to monocytosis, increased resident tissue macrophages in the liver and all major organs, and liver growth. These effects were associated with increased hepatic glucose uptake and extensive mobilization of body fat. The impacts of CSF1 on macrophage abundance, liver size, and body composition were rapidly reversed to restore homeostasis. The effects of CSF1 on metabolism were independent of several known endocrine regulators and did not impact the physiological fasting response. Analysis using implantable telemetry in metabolic cages revealed progressively reduced body temperature and physical activity with no change in diurnal food intake. These results demonstrate the existence of a dynamic equilibrium between CSF1, the mononuclear phagocyte system, and control of liver-to-body weight ratio, which in turn controls systemic metabolic homeostasis. This novel macrophage regulatory axis has the potential to promote fat mobilization, without changes in appetence, which may have novel implications for managing metabolic syndrome.NEW & NOTEWORTHY CSF1 administration expands tissue macrophages, which transforms systemic metabolism. CSF1 drives fat mobilization and glucose uptake to support liver growth. The effects of CSF1 are independent of normal hormonal metabolic regulation. The effects of CSF1 are rapidly reversible, restoring homeostatic body composition. CSF1-dependent macrophages and liver size are coupled in a dynamic equilibrium.

Genetic models of fibrillinopathies

The fibrillinopathies represent a group of diseases in which the 10-12 nm extracellular microfibrils are disrupted by genetic variants in one of the genes encoding fibrillin molecules, large glycoproteins of the extracellular matrix. The best-known fibrillinopathy is Marfan syndrome, an autosomal dominant condition affecting the cardiovascular, ocular, skeletal, and other systems, with a prevalence of around 1 in 3,000 across all ethnic groups. It is caused by variants of the FBN1 gene, encoding fibrillin-1, which interacts with elastin to provide strength and elasticity to connective tissues. A number of mouse models have been created in an attempt to replicate the human phenotype, although all have limitations. There are also natural bovine models and engineered models in pig and rabbit. Variants in FBN2 encoding fibrillin-2 cause congenital contractural arachnodactyly and mouse models for this condition have also been produced. In most animals, including birds, reptiles, and amphibians, there is a third fibrillin, fibrillin-3 (FBN3 gene) for which the creation of models has been difficult as the gene is degenerate and nonfunctional in mice and rats. Other eukaryotes such as the nematode C. elegans and zebrafish D. rerio have a gene with some homology to fibrillins and models have been used to discover more about the function of this family of proteins. This review looks at the phenotype, inheritance, and relevance of the various animal models for the different fibrillinopathies.

Intraperitoneal transfer of wild-type bone marrow repopulates tissue macrophages in the Csf1r knockout rat without contributing to monocytopoiesis

Homozygous null mutation of the Csf1r gene (Csf1rko) in rats leads to the loss of most tissue macrophage populations and pleiotropic impacts on postnatal growth and organ maturation, leading to early mortality. The phenotype can be reversed by intraperitoneal transfer of WT BM cells (BMT) at weaning. Here, we used a Csf1r-mApple transgenic reporter to track the fate of donor-derived cells. Following BMT into Csf1rko recipients, mApple+ve cells restored IBA1+ tissue macrophage populations in every tissue. However, monocytes, neutrophils, and B cells in the BM, blood, and lymphoid tissues remained of recipient (mApple-ve ) origin. An mApple+ve cell population expanded in the peritoneal cavity and invaded locally in the mesentery, fat pads, omentum, and diaphragm. One week after BMT, distal organs contained foci of mApple+ve , IBA1-ve immature progenitors that appeared to proliferate, migrate, and differentiate locally. We conclude that rat BM contains progenitor cells that are able to restore, replace, and maintain all tissue macrophage populations in a Csf1rko rat directly without contributing to the BM progenitor or blood monocyte populations.

Defining infertility: a qualitative interview study of patients and physicians

PURPOSE

To investigate if infertility patients and physicians apply a traditional biomedical model of disease in their conceptualisation of infertility, examine any contradictions and conflicts in conceptualisations, and examine areas of concordance and discordance between physicians and patients.

METHODS

Semi-structured interviews were conducted with 20 infertility patients and 18 infertility physicians between September 2010 and April 2012. Interviews were analysed qualitatively to determine physician and patient conceptualisations of infertility, reactions to the definition of infertility as a disease, and potential benefits and concerns related to application of a disease label to the condition.

RESULTS

Most physicians (n = 14/18) and a minority of patients (n = 6/20) were supportive of defining infertility as a disease. Many of the patients who agreed with classifying infertility as a disease expressed that they had not personally defined it as such previously. Physicians (n = 14) and patients (n = 13) described potential benefits of a disease label, including increases in research funding, insurance coverage, and social acceptability. Some patients (n = 10) described potential stigma as a negative consequence. When describing appraisals of infertility, both physicians (n = 7) and patients (n = 8) invoked religious/spiritual concepts. The potential for religious/spiritual appraisal to contribute to stigmatising or de-stigmatising infertility was discussed.

CONCLUSION

Our findings contradict the assumption that infertility physicians and patients are fully supportive of defining infertility as a disease. While potential benefits of the disease label were recognised by both groups, caution against potential for stigmatisation and unsolicited invocation of religion/spirituality suggest a more holistic model may be appropriate.

Fibrillin-1 and asprosin, novel players in metabolic syndrome

Fibrillin-1 is a major component of the extracellular microfibrils, where it interacts with other extracellular matrix proteins to provide elasticity to connective tissues, and regulates the bioavailability of TGFβ family members. A peptide consisting of the C-terminal 140 amino acids of fibrillin-1 has recently been identified as a glucogenic hormone, secreted from adipose tissue during fasting and targeting the liver to release glucose. This fragment, called asprosin, also signals in the hypothalamus to stimulate appetite. Asprosin levels are correlated with many of the pathologies indicative of metabolic syndrome, including insulin resistance and obesity. Previous studies and reviews have addressed the therapeutic potential of asprosin as a target in obesity, diabetes and related conditions without considering mechanisms underlying the relationship between generation of asprosin and expression of the much larger fibrillin-1 protein. Profibrillin-1 undergoes obligatory cleavage at the cell surface as part of its assembly into microfibrils, producing the asprosin peptide as well as mature fibrillin-1. Patterns of FBN1 mRNA expression are inconsistent with the necessity for regulated release of asprosin. The asprosin peptide may be protected from degradation in adipose tissue. We present evidence for an alternative possibility, that asprosin mRNA is generated independently from an internal promoter within the 3' end of the FBN1 gene, which would allow for regulation independent of fibrillin-synthesis and is more economical of cellular resources. The discovery of asprosin opened exciting possibilities for treatment of metabolic syndrome related conditions, but there is much to be understood before such therapies could be introduced into the clinic.

A kinase-dead Csf1r mutation associated with adult-onset leukoencephalopathy has a dominant inhibitory impact on CSF1R signalling

Amino acid substitutions in the kinase domain of the human CSF1R gene are associated with autosomal dominant adult-onset leukoencephalopathy with axonal spheroids and pigmented glia (ALSP). To model the human disease, we created a disease-associated mutation (pGlu631Lys; E631K) in the mouse Csf1r locus. Homozygous mutation (Csf1rE631K/E631K) phenocopied the Csf1r knockout, with prenatal mortality or severe postnatal growth retardation and hydrocephalus. Heterozygous mutation delayed the postnatal expansion of tissue macrophage populations in most organs. Bone marrow cells from Csf1rE631K/+mice were resistant to CSF1 stimulation in vitro, and Csf1rE631K/+ mice were unresponsive to administration of a CSF1-Fc fusion protein, which expanded tissue macrophage populations in controls. In the brain, microglial cell numbers and dendritic arborisation were reduced in Csf1rE631K/+ mice, as in patients with ALSP. The microglial phenotype is the opposite of microgliosis observed in Csf1r+/- mice. However, we found no evidence of brain pathology or impacts on motor function in aged Csf1rE631K/+ mice. We conclude that heterozygous disease-associated CSF1R mutations compromise CSF1R signalling. We speculate that leukoencephalopathy associated with dominant human CSF1R mutations requires an environmental trigger and/or epistatic interaction with common neurodegenerative disease-associated alleles.

Generation and network analysis of an RNA-seq transcriptional atlas for the rat

The laboratory rat is an important model for biomedical research. To generate a comprehensive rat transcriptomic atlas, we curated and downloaded 7700 rat RNA-seq datasets from public repositories, downsampled them to a common depth and quantified expression. Data from 585 rat tissues and cells, averaged from each BioProject, can be visualized and queried at http://biogps.org/ratatlas. Gene co-expression network (GCN) analysis revealed clusters of transcripts that were tissue or cell type restricted and contained transcription factors implicated in lineage determination. Other clusters were enriched for transcripts associated with biological processes. Many of these clusters overlap with previous data from analysis of other species, while some (e.g. expressed specifically in immune cells, retina/pineal gland, pituitary and germ cells) are unique to these data. GCN analysis on large subsets of the data related specifically to liver, nervous system, kidney, musculoskeletal system and cardiovascular system enabled deconvolution of cell type-specific signatures. The approach is extensible and the dataset can be used as a point of reference from which to analyse the transcriptomes of cell types and tissues that have not yet been sampled. Sets of strictly co-expressed transcripts provide a resource for critical interpretation of single-cell RNA-seq data.

Fragmentation of tissue-resident macrophages during isolation confounds analysis of single-cell preparations from mouse hematopoietic tissues

Mouse hematopoietic tissues contain abundant tissue-resident macrophages that support immunity, hematopoiesis, and bone homeostasis. A systematic strategy to characterize macrophage subsets in mouse bone marrow (BM), spleen, and lymph node unexpectedly reveals that macrophage surface marker staining emanates from membrane-bound subcellular remnants associated with unrelated cells. Intact macrophages are not present within these cell preparations. The macrophage remnant binding profile reflects interactions between macrophages and other cell types in vivo. Depletion of CD169⁺ macrophages in vivo eliminates F4/80⁺ remnant attachment. Remnant-restricted macrophage-specific membrane markers, cytoplasmic fluorescent reporters, and mRNA are all detected in non-macrophage cells including isolated stem and progenitor cells. Analysis of RNA sequencing (RNA-seq) data, including publicly available datasets, indicates that macrophage fragmentation is a general phenomenon that confounds bulk and single-cell analysis of disaggregated hematopoietic tissues. Hematopoietic tissue macrophage fragmentation undermines the accuracy of macrophage ex vivo molecular profiling and creates opportunity for misattribution of macrophage-expressed genes to non-macrophage cells.

Macrophages form erythropoietic niches and regulate iron homeostasis to adapt erythropoiesis in response to infections and inflammation

It has recently emerged that tissue-resident macrophages are key regulators of several stem cell niches orchestrating tissue formation during development, as well as postnatally, when they also organize the repair and regeneration of many tissues including the hemopoietic tissue. The fact that macrophages are also master regulators and effectors of innate immunity and inflammation allows them to coordinate hematopoietic response to infections, injuries, and inflammation. After recently reviewing the roles of phagocytes and macrophages in regulating normal and pathologic hematopoietic stem cell niches, we now focus on the key roles of macrophages in regulating erythropoiesis and iron homeostasis. We review herein the recent advances in understanding how macrophages at the center of erythroblastic islands form an erythropoietic niche that controls the terminal differentiation and maturation of erythroblasts into reticulocytes; how red pulp macrophages in the spleen control iron recycling and homeostasis; how these macrophages coordinate emergency erythropoiesis in response to blood loss, infections, and inflammation; and how persistent infections or inflammation can lead to anemia of inflammation via macrophages. Finally, we discuss the technical challenges associated with the molecular characterization of erythroid island macrophages and red pulp macrophages.

Role of macrophages and phagocytes in orchestrating normal and pathologic hematopoietic niches

The bone marrow (BM) contains a mosaic of niches specialized in supporting different maturity stages of hematopoietic stem and progenitor cells such as hematopoietic stem cells and myeloid, lymphoid, and erythroid progenitors. Recent advances in BM imaging and conditional gene knockout mice have revealed that niches are a complex network of cells of mesenchymal, endothelial, neuronal, and hematopoietic origins, together with local physicochemical parameters. Within these complex structures, phagocytes, such as neutrophils, macrophages, and dendritic cells, all of which are of hematopoietic origin, have been found to be important in regulating several niches in the BM, including hematopoietic stem cell niches, erythropoietic niches, and niches involved in endosteal bone formation. There is also increasing evidence that these macrophages have an important role in adapting hematopoiesis, erythropoiesis, and bone formation in response to inflammatory stressors and play a key part in maintaining the integrity and function of these. Likewise, there is also accumulating evidence that subsets of monocytes, macrophages, and other phagocytes contribute to the progression and response to treatment of several lymphoid malignancies such as multiple myeloma, Hodgkin lymphoma, and non-Hodgkin lymphoma, as well as lymphoblastic leukemia, and may also play a role in myelodysplastic syndrome and myeloproliferative neoplasms associated with Noonan syndrome and aplastic anemia. In this review, the potential functions of macrophages and other phagocytes in normal and pathologic niches are discussed, as are the challenges in studying BM and other tissue-resident macrophages at the molecular level.

CSF1R-dependent macrophages control postnatal somatic growth and organ maturation

Homozygous mutation of the Csf1r locus (Csf1rko) in mice, rats and humans leads to multiple postnatal developmental abnormalities. To enable analysis of the mechanisms underlying the phenotypic impacts of Csf1r mutation, we bred a rat Csf1rko allele to the inbred dark agouti (DA) genetic background and to a Csf1r-mApple reporter transgene. The Csf1rko led to almost complete loss of embryonic macrophages and ablation of most adult tissue macrophage populations. We extended previous analysis of the Csf1rko phenotype to early postnatal development to reveal impacts on musculoskeletal development and proliferation and morphogenesis in multiple organs. Expression profiling of 3-week old wild-type (WT) and Csf1rko livers identified 2760 differentially expressed genes associated with the loss of macrophages, severe hypoplasia, delayed hepatocyte maturation, disrupted lipid metabolism and the IGF1/IGF binding protein system. Older Csf1rko rats developed severe hepatic steatosis. Consistent with the developmental delay in the liver Csf1rko rats had greatly-reduced circulating IGF1. Transfer of WT bone marrow (BM) cells at weaning without conditioning repopulated resident macrophages in all organs, including microglia in the brain, and reversed the mutant phenotypes enabling long term survival and fertility. WT BM transfer restored osteoclasts, eliminated osteopetrosis, restored bone marrow cellularity and architecture and reversed granulocytosis and B cell deficiency. Csf1rko rats had an elevated circulating CSF1 concentration which was rapidly reduced to WT levels following BM transfer. However, CD43hi non-classical monocytes, absent in the Csf1rko, were not rescued and bone marrow progenitors remained unresponsive to CSF1. The results demonstrate that the Csf1rko phenotype is autonomous to BM-derived cells and indicate that BM contains a progenitor of tissue macrophages distinct from hematopoietic stem cells. The model provides a unique system in which to define the pathways of development of resident tissue macrophages and their local and systemic roles in growth and organ maturation.

CRISPR-Cas9 Editing of Human Histone Deubiquitinase Gene USP16 in Human Monocytic Leukemia Cell Line THP-1

USP16 is a histone deubiquitinase which facilitates G2/M transition during the cell cycle, regulates DNA damage repair and contributes to inducible gene expression. We mutated the USP16 gene in a high differentiation clone of the acute monocytic leukemia cell line THP-1 using the CRISPR-Cas9 system and generated four homozygous knockout clones. All were able to proliferate and to differentiate in response to phorbol ester (PMA) treatment. One line was highly proliferative prior to PMA treatment and shut down proliferation upon differentiation, like wild type. Three clones showed sustained expression of the progenitor cell marker MYB, indicating that differentiation had not completely blocked proliferation in these clones. Network analysis of transcriptomic differences among wild type, heterozygotes and homozygotes showed clusters of genes that were up- or down-regulated after differentiation in all cell lines. Prior to PMA treatment, the homozygous clones had lower levels than wild type of genes relating to metabolism and mitochondria, including SRPRB, encoding an interaction partner of USP16. There was also apparent loss of interferon signaling. In contrast, a number of genes were up-regulated in the homozygous cells compared to wild type at baseline, including other deubiquitinases (USP12, BAP1, and MYSM1). However, three homozygotes failed to fully induce USP3 during differentiation. Other network clusters showed effects prior to or after differentiation in the homozygous clones. Thus the removal of USP16 affected the transcriptome of the cells, although all these lines were able to survive, which suggests that the functions attributed to USP16 may be redundant. Our analysis indicates that the leukemic line can adapt to the extreme selection pressure applied by the loss of USP16, and the harsh conditions of the gene editing and selection protocol, through different compensatory pathways. Similar selection pressures occur during the evolution of a cancer in vivo, and our results can be seen as a case study in leukemic cell adaptation. USP16 has been considered a target for cancer chemotherapy, but our results suggest that treatment would select for escape mutants that are resistant to USP16 inhibitors.

The Mononuclear Phagocyte System of the Rat

The laboratory rat continues to be the model of choice for many studies of physiology, behavior, and complex human diseases. Cells of the mononuclear phagocyte system (MPS; monocytes, macrophages, and dendritic cells) are abundant residents in every tissue in the body and regulate postnatal development, homeostasis, and innate and acquired immunity. Recruitment and proliferation of MPS cells is an essential component of both initiation and resolution of inflammation. The large majority of current knowledge of MPS biology is derived from studies of inbred mice, but advances in technology and resources have eliminated many of the advantages of the mouse as a model. In this article, we review the tools available and the current state of knowledge of development, homeostasis, regulation, and diversity within the MPS of the rat.

Microdeletion of 9q22.3: A patient with minimal deletion size associated with a severe phenotype

Basal cell nevus syndrome (also known as Gorlin Syndrome; MIM109400) is an autosomal dominant disorder characterized by recurrent pathological features such as basal cell carcinomas and odontogenic keratocysts as well as skeletal abnormalities. Most affected individuals have point mutations or small insertions or deletions within the PTCH1 gene on human chromosome 9, but there are some cases with more extensive deletion of the region, usually including the neighboring FANCC and/or ERCC6L2 genes. We report a 16‐year‐old patient with a deletion of approximately 400,000 bases which removes only PTCH1 and some non‐coding RNA genes but leaves FANCC and ERCC6L2 intact. In spite of the small amount of DNA for which he is haploid, his phenotype is more extreme than many individuals with longer deletions in the region. This includes early presentation with a large number of basal cell nevi and other skin lesions, multiple jaw keratocysts, and macrosomia. We found that the deletion was in the paternal chromosome, in common with other macrosomia cases. Using public databases, we have examined possible interactions between sequences within and outside the deletion and speculate that a regulatory relationship exists with flanking genes, which is unbalanced by the deletion, resulting in abnormal activation or repression of the target genes and hence the severity of the phenotype.

Clinical and Echocardiographic Findings in an Aged Population of Cavalier King Charles Spaniels

Myxomatous mitral valve disease (MMVD) is the most common cardiac disease in dogs. It varies from dogs without clinical signs to those developing left-sided congestive heart failure, leading to death. Cavalier King Charles Spaniels (CKCSs) are particularly susceptible to MMVD. We hypothesised that within the elderly CKCS population, there is a sub-cohort of MMVD-affected dogs that do not have cardiac remodelling. The objectives of the present study were (i) to determine the prevalence and the degree of cardiac remodelling associated with MMVD; and (ii) assess the effect of age, gender, and body weight on echocardiographic status in a population of aged CKCSs. A total of 126 CKCSs ≥ 8 years old were prospectively included. They all had a physical and echocardiographic examination. A systolic murmur was detected in 89% of dogs; the presence of clinical signs was reported in 19% of them; and echocardiographic evidence of MMVD was described in 100%. Despite the high prevalence, 44.4% of the dogs were clear of echocardiographic signs of cardiac remodelling. Age was significantly associated with the presence and severity of cardiac remodelling and mitral valve prolapse. Our results showed that a proportion of elderly CKCS with confirmed MMVD did not undergo advanced stages of this pathology.

Embryo blastulation and quality between days 5 and 6 of extended embryo culture

PURPOSE

This study aims to know what proportion of culture day 5 pre-blastocyst-stage embryos develop into blastocysts by culture day 6 and what patient and cycle characteristics are associated with delayed blastocyst formation.

METHODS

A retrospective observational cohort analysis was performed including a total of 9886 embryos from 1008 IVF cycles in 835 patients, who underwent treatment between January 1, 2016, and December 31, 2018. Autologous fresh in vitro fertilization (IVF) cycles at a single academic center were included in the analysis. Embryos were group-cultured using single-step culture media. Blastulation was defined as the presence of a new blastocyst. Usable blastulation was defined as the presence of a new good or excellent quality, expanded, hatching, or hatched blastocysts.

RESULTS

The mean blastulation rate between days 5 and 6 of extended embryo culture was 30.9%. The mean percentage of embryos developing into usable blastocyst-stage embryos was 19.8%. The factors associated with blastulation on day 6 included the total number of embryos and the number of pre-blastocysts on day 5, as well as the use of ICSI. Age, the number of total embryos, those remained in culture and pre-blastocysts, as well as the blastulation rate on day 5 were associated with usable blastulation.

CONCLUSION

It is important to know the usable blastocyst development rate between culture days 5 and 6 in order to adequately counsel patients debating whether to proceed with fresh ET on day 5 or forego ET with the expectation that embryos will be biopsied for PGT and/or cryopreserved on culture day 6. Our findings provide evidence to help guide patients in this difficult decision.

Analysis of homozygous and heterozygous Csf1r knockout in the rat as a model for understanding microglial function in brain development and the impacts of human CSF1R mutations

Mutations in the human CSF1R gene have been associated with dominant and recessive forms of neurodegenerative disease. Here we describe the impacts of Csf1r mutation in the rat on development of the brain. Diffusion imaging indicated small reductions in major fiber tracts that may be associated in part with ventricular enlargement. RNA-seq profiling revealed a set of 105 microglial markers depleted in all brain regions of the Csf1rko rats. There was no evidence of region or sex-specific expression of microglia-associated transcripts. Other than the microglial signature, Csf1rko had no effect on any neuronal or region-specific transcript cluster. Expression of markers of oligodendrocytes, astrocytes, dopaminergic neurons and Purkinje cells was minimally affected. However, there were defects in dendritic arborization of doublecortin-positive neurogenic precursors and expression of poly-sialylated neural cell adhesion molecule (PS-NCAM) in the dentate gyrus of the hippocampus. Heterozygous Csf1rko rats had no detectable brain phenotype. We conclude that most brain developmental processes occur normally in the absence of microglia and that CSF1R haploinsufficiency is unlikely to cause leukoencephalopathy.

Network analysis of transcriptomic diversity amongst resident tissue macrophages and dendritic cells in the mouse mononuclear phagocyte system

The mononuclear phagocyte system (MPS) is a family of cells including progenitors, circulating blood monocytes, resident tissue macrophages, and dendritic cells (DCs) present in every tissue in the body. To test the relationships between markers and transcriptomic diversity in the MPS, we collected from National Center for Biotechnology Information Gene Expression Omnibus (NCBI-GEO) a total of 466 quality RNA sequencing (RNA-seq) data sets generated from mouse MPS cells isolated from bone marrow, blood, and multiple tissues. The primary data were randomly downsized to a depth of 10 million reads and requantified. The resulting data set was clustered using the network analysis tool BioLayout. A sample-to-sample matrix revealed that MPS populations could be separated based upon tissue of origin. Cells identified as classical DC subsets, cDC1s and cDC2s, and lacking Fcgr1 (encoding the protein CD64) were contained within the MPS cluster, no more distinct than other MPS cells. A gene-to-gene correlation matrix identified large generic coexpression clusters associated with MPS maturation and innate immune function. Smaller coexpression gene clusters, including the transcription factors that drive them, showed higher expression within defined isolated cells, including monocytes, macrophages, and DCs isolated from specific tissues. They include a cluster containing Lyve1 that implies a function in endothelial cell (EC) homeostasis, a cluster of transcripts enriched in intestinal macrophages, and a generic lymphoid tissue cDC cluster associated with Ccr7. However, transcripts encoding Adgre1, Itgax, Itgam, Clec9a, Cd163, Mertk, Mrc1, Retnla, and H2-a/e (encoding class II major histocompatibility complex [MHC] proteins) and many other proposed macrophage subset and DC lineage markers each had idiosyncratic expression profiles. Coexpression of immediate early genes (for example, Egr1, Fos, Dusp1) and inflammatory cytokines and chemokines (tumour necrosis factor [Tnf], Il1b, Ccl3/4) indicated that all tissue disaggregation and separation protocols activate MPS cells. Tissue-specific expression clusters indicated that all cell isolation procedures also co-purify other unrelated cell types that may interact with MPS cells in vivo. Comparative analysis of RNA-seq and single-cell RNA-seq (scRNA-seq) data from the same lung cell populations indicated that MPS heterogeneity implied by global cluster analysis may be even greater at a single-cell level. This analysis highlights the power of large data sets to identify the diversity of MPS cellular phenotypes and the limited predictive value of surface markers to define lineages, functions, or subpopulations.

Comparative transcriptomic profiling of myxomatous mitral valve disease in the cavalier King Charles spaniel

BACKGROUND

Almost all elderly dogs develop myxomatous mitral valve disease by the end of their life, but the cavalier King Charles spaniel (CKCS) has a heightened susceptibility, frequently resulting in death at a young age and suggesting that there is a genetic component to the condition in this breed. Transcriptional profiling can reveal the impact of genetic variation through differences in gene expression levels. The aim of this study was to determine whether expression patterns were different in mitral valves showing myxomatous degeneration from CKCS dogs compared to valves from non-CKCS dogs.

RESULTS

Gene expression patterns in three groups of canine valves resulted in distinct separation of normal valves, diseased valves from CKCS and diseased valves from other breeds; the latter were more similar to the normal valves than were the valves from CKCS. Gene expression patterns in diseased valves from CKCS dogs were quite different from those in the valves from other dogs, both affected and normal. Patterns in all diseased valves (from CKCS and other breeds) were also somewhat different from normal non-diseased samples. Analysis of differentially expressed genes showed enrichment in GO terms relating to cardiac development and function and to calcium signalling canonical pathway in the genes down-regulated in the diseased valves from CKCS, compared to normal valves and to diseased valves from other breeds. F2 (prothrombin) (CKCS diseased valves compared to normal) and MEF2C pathway activation (CKCS diseased valves compared to non-CKCS diseased valves) had the strongest association with the gene changes. A large number of genes that were differentially expressed in the CKCS diseased valves compared with normal valves and diseased valves from other breeds were associated with cardiomyocytes including CASQ2, TNNI3 and RYR2.

CONCLUSION

Transcriptomic profiling identified gene expression changes in CKCS diseased valves that were not present in age and disease severity-matched non-CKCS valves. These genes are associated with cardiomyocytes, coagulation and extra-cellular matrix remodelling. Identification of genes that vary in the CKCS will allow exploration of genetic variation to understand the aetiology of the disease in this breed, and ultimately development of breeding strategies to eliminate this disease from the breed.

Expression of Calcification and Extracellular Matrix Genes in the Cardiovascular System of the Healthy Domestic Sheep (Ovis aries)

The maintenance of a healthy cardiovascular system requires expression of genes that contribute to essential biological activities and repression of those that are associated with functions likely to be detrimental to cardiovascular homeostasis. Vascular calcification is a major disruption to cardiovascular homeostasis, where tissues of the cardiovascular system undergo ectopic calcification and consequent dysfunction, but little is known about the expression of calcification genes in the healthy cardiovascular system. Large animal models are of increasing importance in cardiovascular disease research as they demonstrate more similar cardiovascular features (in terms of anatomy, physiology and size) to humans than do rodent species. We used RNA sequencing results from the sheep, which has been utilized extensively to examine calcification of prosthetic cardiac valves, to explore the transcriptome of the heart and cardiac valves in this large animal, in particular looking at expression of calcification and extracellular matrix genes. We then examined genes implicated in the process of vascular calcification in a wide array of cardiovascular tissues and across multiple developmental stages, using RT-qPCR. Our results demonstrate that there is a balance between genes that promote and those that suppress mineralization during development and across cardiovascular tissues. We show extensive expression of genes encoding proteins involved in formation and maintenance of the extracellular matrix in cardiovascular tissues, and high expression of hematopoietic genes in the cardiac valves. Our analysis will support future research into the functions of implicated genes in the development of valve calcification, and increase the utility of the sheep as a large animal model for understanding ectopic calcification in cardiovascular disease. This study provides a foundation to explore the transcriptome of the developing cardiovascular system and is a valuable resource for the fields of mammalian genomics and cardiovascular research.

Influence of the MUC1 Cell Surface Mucin on Gastric Mucosal Gene Expression Profiles in Response to Helicobacter pylori Infection in Mice

The cell surface mucin MUC1 is an important host factor limiting Helicobacter pylori (H. pylori) pathogenesis in both humans and mice by providing a protective barrier and modulating mucosal epithelial and leukocyte responses. The aim of this study was to establish the time-course of molecular events in MUC1-modulated gene expression profiles in response to H. pylori infection in wild type (WT) and MUC1-deficient mice using microarray-determined mRNA expression, gene network analysis and Ingenuity Pathway Analysis (IPA). A time-course over the first 72 h of infection showed significantly higher mucosal loads of bacteria at 8 h of infection in Muc1^−/− mice compared with WT, confirming its importance in the early stages of infection (P = 0.0003). Microarray analysis revealed 266 differentially expressed genes at one or more time-points over 72 h in the gastric mucosa of Muc1^−/− mice compared with WT control using a threshold of 2-fold change. The SPINK1 pancreatic cancer canonical pathway was strongly inhibited in Muc1^−/− mice compared with WT at sham and 8 h infection (P = 6.08E-14 and P = 2.25 E-19, respectively) but potently activated at 24 and 72 h post-infection (P = 1.38E-22 and P = 5.87E-13, respectively). The changes in this pathway are reflective of higher expression of genes mediating digestion and absorption of lipids, carbohydrates, and proteins at sham and 8 h infection in the absence of MUC1, but that this transcriptional signature is highly down regulated as infection progresses in the absence of MUC1. Uninfected Muc1^−/− gastric tissue was highly enriched for expression of factors involved in lipid metabolism and 8 h infection further activated this network compared with WT. As infection progressed, a network of antimicrobial and anti-inflammatory response genes was more highly activated in Muc1^−/− than WT mice. Key target genes identified by time-course microarrays were independently validated using RT-qPCR. These results highlight the dynamic interplay between the host and H. pylori, and the role of MUC1 in host defense, and provide a general picture of changes in cellular gene expression modulated by MUC1 in a time-dependent manner in response to H. pylori infection.

Species-Specificity of Transcriptional Regulation and the Response to Lipopolysaccharide in Mammalian Macrophages

Mammalian macrophages differ in their basal gene expression profiles and response to the toll-like receptor 4 (TLR4) agonist, lipopolysaccharide (LPS). In human macrophages, LPS elicits a temporal cascade of transient gene expression including feed forward activators and feedback regulators that limit the response. Here we present a transcriptional network analysis of the response of sheep bone marrow-derived macrophages (BMDM) to LPS based upon RNA-seq at 0, 2, 4, 7, and 24 h post-stimulation. The analysis reveals a conserved transcription factor network with humans, and rapid induction of feedback regulators that constrain the response at every level. The gene expression profiles of sheep BMDM at 0 and 7 h post LPS addition were compared to similar data obtained from goat, cow, water buffalo, horse, pig, mouse and rat BMDM. This comparison was based upon identification of 8,200 genes annotated in all species and detected at >10TPM in at least one sample. Analysis of expression of transcription factors revealed a conserved transcriptional millieu associated with macrophage differentiation and LPS response. The largest co-expression clusters, including genes encoding cell surface receptors, endosome–lysosome components and secretory activity, were also expressed in all species and the combined dataset defines a macrophage functional transcriptome. All of the large animals differed from rodents in lacking inducible expression of genes involved in arginine metabolism and nitric oxide production. Instead, they expressed inducible transporters and enzymes of tryptophan and kynurenine metabolism. BMDM from all species expressed high levels of transcripts encoding transporters and enzymes involved in glutamine metabolism suggesting that glutamine is a major metabolic fuel. We identify and discuss transcripts that were uniquely expressed or regulated in rodents compared to large animals including ACOD1, CXC and CC chemokines, CD163, CLEC4E, CPM, CSF1, CSF2, CTSK, MARCO, MMP9, SLC2A3, SLC7A7, and SUCNR1. Conversely, the data confirm the conserved regulation of multiple transcripts for which there is limited functional data from mouse models and knockouts. The data provide a resource for functional annotation and interpretation of loci involved in susceptibility to infectious and inflammatory disease in humans and large animal species.

The Transcriptional Network That Controls Growth Arrest and Macrophage Differentiation in the Human Myeloid Leukemia Cell Line THP-1

The response of the human acute myeloid leukemia cell line THP-1 to phorbol esters has been widely studied to test candidate leukemia therapies and as a model of cell cycle arrest and monocyte-macrophage differentiation. Here we have employed Cap Analysis of Gene Expression (CAGE) to analyze a dense time course of transcriptional regulation in THP-1 cells treated with phorbol myristate acetate (PMA) over 96 h. PMA treatment greatly reduced the numbers of cells entering S phase and also blocked cells exiting G2/M. The PMA-treated cells became adherent and expression of mature macrophage-specific genes increased progressively over the duration of the time course. Within 1–2 h PMA induced known targets of tumor protein p53 (TP53), notably CDKN1A, followed by gradual down-regulation of cell-cycle associated genes. Also within the first 2 h, PMA induced immediate early genes including transcription factor genes encoding proteins implicated in macrophage differentiation (EGR2, JUN, MAFB) and down-regulated genes for transcription factors involved in immature myeloid cell proliferation (MYB, IRF8, GFI1). The dense time course revealed that the response to PMA was not linear and progressive. Rather, network-based clustering of the time course data highlighted a sequential cascade of transient up- and down-regulated expression of genes encoding feedback regulators, as well as transcription factors associated with macrophage differentiation and their inferred target genes. CAGE also identified known and candidate novel enhancers expressed in THP-1 cells and many novel inducible genes that currently lack functional annotation and/or had no previously known function in macrophages. The time course is available on the ZENBU platform allowing comparison to FANTOM4 and FANTOM5 data.

Disease Severity-Associated Gene Expression in Canine Myxomatous Mitral Valve Disease Is Dominated by TGFβ Signaling

Myxomatous mitral valve disease (MMVD) is the most common acquired canine cardiovascular disease and shares many similarities with human mitral valvulopathies. While transcriptomic datasets are available for the end-stage disease in both species, there is no information on how gene expression changes as the disease progresses, such that it cannot be stated with certainty if the changes seen in end-stage disease are casual or consequential. In contrast to humans, the disease in dogs can be more readily examined as it progresses, and this allows an opportunity for insight into disease pathogenesis relevant to both species. The aim of this study was to identify changes in valve gene expression as canine MMVD advances over an entire life-time, from normal (grade 0) to severely affected (grade 4), and differences in gene expression comparing normal and disease areas of the same valve. Transcriptomic profiling identified 1002 differentially expressed genes (DEGs) across all four disease grades when compared with normal valves with the greatest number of DEGs in grade 3 (673) and grade 4 (507). DEGs were associated with a large number of gene families, including genes encoding cytoskeletal filaments, peptidases, extra-cellular matrix (ECM) proteins, chemokines and integrins. Gene enrichment analysis identified significant grade-dependent changes in gene clustering, with clusters trending both up and down as disease progressed. Significant grade-dependent changes in hallmark disease gene expression intensity were identified, including ACTA2, HTR2B, MMP12, and CDKN2A. Gene Ontology terms were dominated by terms for ECM and inflammation with TGFβ1, TNF, IFGN identified as the top up-stream regulators in both whole and dissected diseased valve samples. These data show that while disease progression in MMVD is associated with increasing numbers of DEGs, TGFβ appears to be the dominant signaling pathway controlling pathogenesis irrespective of disease severity.

Functional Annotation of the Transcriptome of the Pig, Sus scrofa, Based Upon Network Analysis of an RNAseq Transcriptional Atlas

The domestic pig (Sus scrofa) is both an economically important livestock species and a model for biomedical research. Two highly contiguous pig reference genomes have recently been released. To support functional annotation of the pig genomes and comparative analysis with large human transcriptomic data sets, we aimed to create a pig gene expression atlas. To achieve this objective, we extended a previous approach developed for the chicken. We downloaded RNAseq data sets from public repositories, down-sampled to a common depth, and quantified expression against a reference transcriptome using the mRNA quantitation tool, Kallisto. We then used the network analysis tool Graphia to identify clusters of transcripts that were coexpressed across the merged data set. Consistent with the principle of guilt-by-association, we identified coexpression clusters that were highly tissue or cell-type restricted and contained transcription factors that have previously been implicated in lineage determination. Other clusters were enriched for transcripts associated with biological processes, such as the cell cycle and oxidative phosphorylation. The same approach was used to identify coexpression clusters within RNAseq data from multiple individual liver and brain samples, highlighting cell type, process, and region-specific gene expression. Evidence of conserved expression can add confidence to assignment of orthology between pig and human genes. Many transcripts currently identified as novel genes with ENSSSCG or LOC IDs were found to be coexpressed with annotated neighbouring transcripts in the same orientation, indicating they may be products of the same transcriptional unit. The meta-analytic approach to utilising public RNAseq data is extendable to include new data sets and new species and provides a framework to support the Functional Annotation of Animals Genomes (FAANG) initiative.

Compound heterozygous mutations in FBN1 in a large family with Marfan syndrome

Background

Marfan syndrome (MFS) is a dominant monogenic disorder caused by mutations in fibrillin 1 (FBN1). Rarely, compound heterozygosity for FBN1 mutations has been described.

Methods

A large kindred with MFS was assessed clinically over decades, and genetically using exome and/or Sanger sequencing.

Results

A previously identified FBN1 missense variant (p.Tyr754Cys) was confirmed in all subjects with MFS. An additional variant (p.Met2273Thr), previously associated with incomplete MFS, was identified in three siblings. These three compound heterozygous individuals had aortic dilatation at early age (all <30 years): one also had cerebral and ocular aneurysms; and one, who had undergone surgical repair aged 18 years, died from aortic dissection at 31 years. In contrast, their heterozygous father (p.Tyr754Cys) with MFS died at 57 years (myocardial infarction) without requiring surgical intervention and one heterozygous (p.Tyr754Cys) sibling has aortic dilatation presenting >40 years but not requiring surgical intervention. Another heterozygous (p.Tyr754Cys) sibling did require aortic root repair (28 years). The heterozygous (p.Met2273Thr) mother had aortic dilatation diagnosed at age 68 years but has not required surgical repair.

Conclusion

Although compound heterozygosity or homozygosity is rare in MFS, it should be considered when there is an unusually severe phenotype in a subset of family members.

Visualization and analysis of RNA-Seq assembly graphs

RNA-Seq is a powerful transcriptome profiling technology enabling transcript discovery and quantification. Whilst most commonly used for gene-level quantification, the data can be used for the analysis of transcript isoforms. However, when the underlying transcript assemblies are complex, current visualization approaches can be limiting, with splicing events a challenge to interpret. Here, we report on the development of a graph-based visualization method as a complementary approach to understanding transcript diversity from short-read RNA-Seq data. Following the mapping of reads to a reference genome, a read-to-read comparison is performed on all reads mapping to a given gene, producing a weighted similarity matrix between reads. This is used to produce an RNA assembly graph, where nodes represent reads and edges similarity scores between them. The resulting graphs are visualized in 3D space to better appreciate their sometimes large and complex topology, with other information being overlaid on to nodes, e.g. transcript models. Here we demonstrate the utility of this approach, including the unusual structure of these graphs and how they can be used to identify issues in assembly, repetitive sequences within transcripts and splice variants. We believe this approach has the potential to significantly improve our understanding of transcript complexity.

Analysis of the Progeny of Sibling Matings Reveals Regulatory Variation Impacting the Transcriptome of Immune Cells in Commercial Chickens

There is increasing recognition that the underlying genetic variation contributing to complex traits influences transcriptional regulation and can be detected at a population level as expression quantitative trait loci. At the level of an individual, allelic variation in transcriptional regulation of individual genes can be detected by measuring allele-specific expression in RNAseq data. We reasoned that extreme variants in gene expression could be identified by analysis of inbred progeny with shared grandparents. Commercial chickens have been intensively selected for production traits. Selection is associated with large blocks of linkage disequilibrium with considerable potential for co-selection of closely linked "hitch-hiker alleles" affecting traits unrelated to the feature being selected, such as immune function, with potential impact on the productivity and welfare of the animals. To test this hypothesis that there is extreme allelic variation in immune-associated genes we sequenced a founder population of commercial broiler and layer birds. These birds clearly segregated genetically based upon breed type. Each genome contained numerous candidate null mutations, protein-coding variants predicted to be deleterious and extensive non-coding polymorphism. We mated selected broiler-layer pairs then generated cohorts of F2 birds by sibling mating of the F1 generation. Despite the predicted prevalence of deleterious coding variation in the genomic sequence of the founders, clear detrimental impacts of inbreeding on survival and post-hatch development were detected in only one F2 sibship of 15. There was no effect on circulating leukocyte populations in hatchlings. In selected F2 sibships we performed RNAseq analysis of the spleen and isolated bone marrow-derived macrophages (with and without lipopolysaccharide stimulation). The results confirm the predicted emergence of very large differences in expression of individual genes and sets of genes. Network analysis of the results identified clusters of co-expressed genes that vary between individuals and suggested the existence of trans-acting variation in the expression in macrophages of the interferon response factor family that distinguishes the parental broiler and layer birds and influences the global response to lipopolysaccharide. This study shows that the impact of inbreeding on immune cell gene expression can be substantial at the transcriptional level, and potentially opens a route to accelerate selection using specific alleles known to be associated with desirable expression levels.

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.

Lysine demethylases KDM6A and UTY: The X and Y of histone demethylation

Histone demethylases remove transcriptional repressive marks from histones in the nucleus. KDM6A (also known as UTX) is a lysine demethylase which acts on the trimethylated lysine at position 27 in histone 3. The KDM6A gene is located on the X chromosome but escapes X inactivation even though it is not located in the pseudoautosomal region. There is a homologue of KDM6A on the Y chromosome, known as UTY. UTY was thought to have lost its demethylase activity and to represent a non-functional remnant of the ancestral KDM6A gene. However, results with knockout mice suggest that the gene is expressed and the protein performs some function within the cell. Female mice with homozygous deletion of Kdm6a do not survive, but hemizygous males are viable, attributed to the presence of the Uty gene. KDM6A is mutated in the human condition Kabuki syndrome type 2 (OMIM 300867) and in many cases of cancer. The amino acid sequence of KDM6A has been conserved across animal phyla, although it is only found on the X chromosome in eutherian mammals. In this review, we reanalyse existing data from various sources (protein sequence comparison, evolutionary genetics, transcription factor binding and gene expression analysis) to determine the function, expression and evolution of KDM6A and UTY and show that UTY has a functional role similar to KDM6A in metabolism and development.

Arginine to Glutamine Variant in Olfactomedin Like 3 (OLFML3) Is a Candidate for Severe Goniodysgenesis and Glaucoma in the Border Collie Dog Breed

Goniodysgenesis is a developmental abnormality of the anterior chamber of the eye. It is generally considered to be congenital in dogs (Canis lupus familiaris), and has been associated with glaucoma and blindness. Goniodysgenesis and early-onset glaucoma initially emerged in Border Collies in Australia in the late 1990s and have subsequently been found in this breed in Europe and the USA. The objective of the present study was to determine the genetic basis of goniodysgenesis in Border Collies. Clinical diagnosis was based on results of examinations by veterinary ophthalmologists of affected and unaffected dogs from eleven different countries. Genotyping using the Illumina high density canine single nucleotide variant genotyping chip was used to identify a candidate genetic region. There was a highly significant peak of association over chromosome 17, with a p-value of 2 × 10-13 Expression profiles and evolutionary conservation of candidate genes were assessed using public databases. Whole genome sequences of three dogs with glaucoma, three severely affected by goniodysgenesis and three unaffected dogs identified a missense variant in the olfactomedin like 3 (OLFML3) gene in all six affected animals. This was homozygous for the risk allele in all nine cases with glaucoma and 12 of 14 other severely affected animals. Of 67 reportedly unaffected animals, only one was homozygous for this variant (offspring of parents both with goniodysgenesis who were also homozygous for the variant). Analysis of pedigree information was consistent with an autosomal recessive mode of inheritance for severe goniodysgenesis (potentially leading to glaucoma) in this breed. The identification of a candidate genetic region and putative causative variant will aid breeders to reduce the frequency of goniodysgenesis and the risk of glaucoma in the Border Collie population.

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.

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.

Pleiotropic Impacts of Macrophage and Microglial Deficiency on Development in Rats with Targeted Mutation of the Csf1r Locus

We have produced Csf1r-deficient rats by homologous recombination in embryonic stem cells. Consistent with the role of Csf1r in macrophage differentiation, there was a loss of peripheral blood monocytes, microglia in the brain, epidermal Langerhans cells, splenic marginal zone macrophages, bone-associated macrophages and osteoclasts, and peritoneal macrophages. Macrophages of splenic red pulp, liver, lung, and gut were less affected. The pleiotropic impacts of the loss of macrophages on development of multiple organ systems in rats were distinct from those reported in mice. Csf1r^-/- rats survived well into adulthood with postnatal growth retardation, distinct skeletal and bone marrow abnormalities, infertility, and loss of visceral adipose tissue. Gene expression analysis in spleen revealed selective loss of transcripts associated with the marginal zone and, in brain regions, the loss of known and candidate novel microglia-associated transcripts. Despite the complete absence of microglia, there was little overt phenotype in brain, aside from reduced myelination and increased expression of dopamine receptor-associated transcripts in striatum. The results highlight the redundant and nonredundant functions of CSF1R signaling and of macrophages in development, organogenesis, and homeostasis.

Environmentally enriched pigs have transcriptional profiles consistent with neuroprotective effects and reduced microglial activity

Environmental enrichment (EE) is widely used to study the effects of external factors on brain development, function and health in rodent models, but very little is known of the effects of EE on the brain in a large animal model such as the pig. Twenty-four young pigs (aged 5 weeks at start of study, 1:1 male: female ratio) were housed in environmentally enriched (EE) pens and provided with additional enrichment stimulation (a bag filled with straw) once daily. Litter, weight and sex matched controls n= (24) were housed in barren (B) conditions. Behaviour was recorded on alternate days from study day 10. After 21 days, RNA-sequencing of the frontal cortex of male piglets culled one hour after the enrichment stimulation, but not those at 4 h after stimulation, showed upregulation of genes involved in neuronal activity and synaptic plasticity in the EE compared to the B condition. This result is mirrored in the behavioural response to the stimulation which showed a peak in activity around the 1 h time-point. By contrast, EE piglets displayed a signature consistent with a relative decrease in microglial activity compared to those in the B condition. These results confirm those from rodents, suggesting that EE may also confer neuronal health benefits in large mammal models, through a potential relative reduction in neuroinflammatory process and increase in neuroprotection driven by an enrichment-induced increase in behavioural activity.

Exploiting novel valve interstitial cell lines to study calcific aortic valve disease

Calcific aortic valve disease (CAVD) involves progressive valve leaflet thickening and severe calcification, impairing leaflet motion. The in vitro calcification of primary rat, human, porcine and bovine aortic valve interstitial cells (VICs) is commonly employed to investigate CAVD mechanisms. However, to date, no published studies have utilised cell lines to investigate this process. The present study has therefore generated and evaluated the calcification potential of immortalized cell lines derived from sheep and rat VICs. Immortalised sheep (SAVIC) and rat (RAVIC) cell lines were produced by transduction with a recombinant lentivirus encoding the Simian virus (SV40) large and small T antigens (sheep), or large T antigen only (rat), which expressed markers of VICs (vimentin and α‑smooth muscle actin). Calcification was induced in the presence of calcium (Ca; 2.7 mM) in SAVICs (1.9 fold; P<0.001) and RAVICs (4.6 fold; P<0.01). Furthermore, a synergistic effect of calcium and phosphate was observed (2.7 mM Ca/2.0 mM Pi) on VIC calcification in the two cell lines (P<0.001). Analysis of SAVICs revealed significant increases in the mRNA expression of two key genes associated with vascular calcification in cells cultured under calcifying conditions, runt related transcription factor‑2 (RUNX2;1.3 fold; P<0.05 in 4.5 mM Ca) and sodium‑dependent phosphate transporter‑1 (PiT1; 1.2 fold; P<0.05 in 5.4 mM Ca). A concomitant decrease in the expression of the calcification inhibitor matrix Gla protein (MGP) was noted at 3.6 mM Ca (1.3 fold; P<0.01). Assessment of RAVICs revealed alterations in Runx2, Pit1 and Mgp mRNA expression levels (P<0.01). Furthermore, a significant reduction in calcification was observed in SAVICs following treatment with established calcification inhibitors, pyrophosphate (1.8 fold; P<0.01) and etidronate (3.2 fold; P<0.01). Overall, the present study demonstrated that the use of immortalised sheep and rat VIC cell lines is a convenient and cost effective system to investigate CAVD in vitro, and will make a useful contribution to increasing current understanding of the pathophysiological process.

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.

Comparative Transcriptomic Profiling and Gene Expression for Myxomatous Mitral Valve Disease in the Dog and Human

Myxomatous mitral valve disease is the single most important mitral valve disease in both dogs and humans. In the case of the dog it is ubiquitous, such that all aged dogs will have some evidence of the disease, and for humans it is known as Barlow's disease and affects up to 3% of the population, with an expected increase in prevalence as the population ages. Disease in the two species show many similarities and while both have the classic myxomatous degeneration only in humans is there extensive fibrosis. This dual pathology of the human disease markedly affects the valve transcriptome and the difference between the dog and human is dominated by changes in genes associated with fibrosis. This review will briefly examine the comparative valve pathology and then, in more detail, the transcriptomic profiling and gene expression reported so far for both species.

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

BACKGROUND

The 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.

RESULTS

The 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.

CONCLUSION

A 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.

Incidence rates and risk factor analyses for owner reported vomiting and diarrhoea in Labrador Retrievers - findings from the Dogslife Cohort

Dogslife collects data directly from owners of Labrador Retrievers across the UK including information regarding signs of illness irrespective of whether the signs precipitated a veterinary visit. In December 2015, the cohort comprised 6084 dogs aged up to six years and their owners had made 2687 and 2601 reports of diarrhoea and vomiting respectively. The co-occurrence of vomiting and diarrhoea with other signs was described and the frequencies and durations of the two signs were examined with reference to veterinary visitation. Age-specific illness rates were described and Cox Proportional Hazards models were used to estimate risk factors. Just 37% of diarrhoea reports were associated with a veterinary visit and the proportion was even lower for vomiting at 28%; indicating that studies of veterinary practice data miss the majority of signs of gastrointestinal upset. In terms of frequency and duration, diarrhoea typically needed to last two days before the dog would be taken to the vet but if the dog vomited at least every six hours, the owner would be more likely to take the dog to the vet after one day. The illness rates of both signs peaked when the dogs were aged between three and six months. There was also a seasonal pattern to the incidents with the lowest hazards for both in May. Diarrhoea incidents peaked in August-September each year but, while vomiting appeared to be higher in September, it peaked in February. Having another dog in the household was associated with a lower hazard for both vomiting and diarrhoea but having a cat was only associated with a reduced hazard of vomiting. In addition to the distinct seasonal patterns of reporting, there were clear differences in the geographic risks for the two signs. The hazard of diarrhoea was positively associated with human population density within Great Britain (according to home post code) whereas no significant geographical association was found with vomiting. This study is particularly relevant for dog owners because it highlights the wealth of gastrointestinal illnesses in dogs that are dealt with by owners but never seen by veterinarians. The risk factor analyses make use of owner-reported demographic information, highlighting the differences between vomiting and diarrhoea. The analyses give rise to the possibility that the presence of other pets in households may affect rates of illness and indicate new avenues for investigations of these distinct, and oft-suffered conditions.

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

Macrophage colony-stimulating factor (CSF1) is an essential growth and differentiation factor for cells of the macrophage lineage. To explore the role of CSF1 in steady-state control of monocyte production and differentiation and tissue repair, we previously developed a bioactive protein with a longer half-life in circulation by fusing pig CSF1 with the Fc region of pig IgG1a. CSF1-Fc administration to pigs expanded progenitor pools in the marrow and selectively increased monocyte numbers and their expression of the maturation marker CD163. There was a rapid increase in the size of the liver, and extensive proliferation of hepatocytes associated with increased macrophage infiltration. Despite the large influx of macrophages, there was no evidence of liver injury and no increase in circulating liver enzymes. Microarray expression profiling of livers identified increased expression of macrophage markers, i.e., cytokines such as TNF, IL1, and IL6 known to influence hepatocyte proliferation, alongside cell cycle genes. The analysis also revealed selective enrichment of genes associated with portal, as opposed to centrilobular regions, as seen in hepatic regeneration. Combined with earlier data from the mouse, this study supports the existence of a CSF1-dependent feedback loop, linking macrophages of the liver with bone marrow and blood monocytes, to mediate homeostatic control of the size of the liver. The results also provide evidence of safety and efficacy for possible clinical applications of CSF1-Fc.

Expression of FBN1 during adipogenesis: Relevance to the lipodystrophy phenotype in Marfan syndrome and related conditions

Fibrillin-1 is a large glycoprotein encoded by the FBN1 gene in humans. It provides strength and elasticity to connective tissues and is involved in regulating the bioavailability of the growth factor TGFβ. Mutations in FBN1 may be associated with depleted or abnormal adipose tissue, seen in some patients with Marfan syndrome and lipodystrophies. As this lack of adipose tissue does not result in high morbidity or mortality, it is generally under-appreciated, but is a cause of psychosocial problems particularly to young patients. We examined the role of fibrillin-1 in adipogenesis. In inbred mouse strains we found significant variation in the level of expression in the Fbn1 gene that correlated with variation in several measures of body fat, suggesting that mouse fibrillin-1 is associated with the level of fat tissue. Furthermore, we found that FBN1 mRNA was up-regulated in the adipose tissue of obese women compared to non-obese, and associated with an increase in adipocyte size. We used human mesenchymal stem cells differentiated in culture to adipocytes to show that fibrillin-1 declines after the initiation of differentiation. Gene expression results from a similar experiment (available through the FANTOM5 project) revealed that the decline in fibrillin-1 protein was paralleled by a decline in FBN1 mRNA. Examination of the FBN1 gene showed that the region commonly affected in FBN1-associated lipodystrophy is highly conserved both across the three human fibrillin genes and across genes encoding fibrillin-1 in vertebrates. These results suggest that fibrillin-1 is involved as the undifferentiated mesenchymal stem cells transition to adipogenesis but then declines as the developing adipocytes take on their final phenotype. Since the C-terminal peptide of fibrillin-1 is a glucogenic hormone, individuals with low fibrillin-1 (for example with FBN1 mutations associated with lipodystrophy) may fail to differentiate adipocytes and/or to accumulate adipocyte lipids, although this still needs to be shown experimentally.

Datasets of genes coexpressed with FBN1 in mouse adipose tissue and during human adipogenesis

This article contains data related to the research article entitled “Expression of FBN1 during adipogenesis: relevance to the lipodystrophy phenotype in Marfan syndrome and related conditions” [1]. The article concerns the expression of FBN1, the gene encoding the extracellular matrix protein fibrillin-1, during adipogenesis in vitro and in relation to adipose tissue in vivo. The encoded protein has recently been shown to produce a short glucogenic peptide hormone, (Romere et al., 2016) [2], and this gene is therefore a key gene for regulating blood glucose levels. FBN1 and coexpressed genes were examined in mouse strains and in human cells undergoing adipogenesis. The data show the genes that were coexpressed with FBN1, including genes coding for other connective tissue proteins and the proteases that modify them and for the transcription factors that control their expression. Data analysed were derived from datasets available in the public domain and the analysis highlights the utility of such datasets for ongoing analysis and hence reduction in the use of experimental animals.

A Deletion in the Canine POMC Gene Is Associated with Weight and Appetite in Obesity-Prone Labrador Retriever Dogs

Sequencing of candidate genes for obesity in Labrador retriever dogs identified a 14 bp deletion in pro-opiomelanocortin (POMC) with an allele frequency of 12%. The deletion disrupts the β-MSH and β-endorphin coding sequences and is associated with body weight (per allele effect of 0.33 SD), adiposity, and greater food motivation. Among other dog breeds, the deletion was only found in the closely related flat-coat retriever (FCR), where it is similarly associated with body weight and food motivation. The mutation is significantly more common in Labrador retrievers selected to become assistance dogs than pets. In conclusion, the deletion in POMC is a significant modifier of weight and appetite in Labrador retrievers and FCRs and may influence other behavioral traits.

Large animal models of cardiovascular disease

The human cardiovascular system is a complex arrangement of specialized structures with distinct functions. The molecular landscape, including the genome, transcriptome and proteome, is pivotal to the biological complexity of both normal and abnormal mammalian processes. Despite our advancing knowledge and understanding of cardiovascular disease (CVD) through the principal use of rodent models, this continues to be an increasing issue in today's world. For instance, as the ageing population increases, so does the incidence of heart valve dysfunction. This may be because of changes in molecular composition and structure of the extracellular matrix, or from the pathological process of vascular calcification in which bone-formation related factors cause ectopic mineralization. However, significant differences between mice and men exist in terms of cardiovascular anatomy, physiology and pathology. In contrast, large animal models can show considerably greater similarity to humans. Furthermore, precise and efficient genome editing techniques enable the generation of tailored models for translational research. These novel systems provide a huge potential for large animal models to investigate the regulatory factors and molecular pathways that contribute to CVD in vivo. In turn, this will help bridge the gap between basic science and clinical applications by facilitating the refinement of therapies for cardiovascular disease.