Rapidly evolving proteins can aid the identification of genes underlying phenotypic adaptation across taxa, but functional and structural elements of genes can also affect evolutionary rates. In plants, the 'edges' of exons, flanking intron junctions, are known to contain splice enhancers and to have a higher degree of conservation compared to the remainder of the coding region. However, the extent to which these regions may be masking indicators of positive selection or account for the relationship between dN/dS and other genomic parameters is unclear. We investigate the effects of exon edge conservation on the relationship of dN/dS to various sequence characteristics and gene expression parameters in the model plant Arabidopsis thaliana. We also obtain lineage-specific dN/dS estimates, making use of the recently sequenced genome of Thellungiella parvula, the second closest sequenced relative after the sister species Arabidopsis lyrata. Overall, we find that the effect of exon edge conservation, as well as the use of lineage-specific substitution estimates, upon dN/dS ratios partly explains the relationship between the rates of protein evolution and expression level. Furthermore, the removal of exon edges shifts dN/dS estimates upwards, increasing the proportion of genes potentially under adaptive selection. We conclude that lineage-specific substitutions and exon edge conservation have an important effect on dN/dS ratios and should be considered when assessing their relationship with other genomic parameters.

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

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.

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.

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.

Alternative splicing is widespread throughout eukaryotic genomes and greatly increases transcriptomic diversity. Many alternative isoforms have functional roles in developmental processes and are precisely temporally regulated. To facilitate the study of alternative splicing in a developmental context, we created MeDAS, a Metazoan Developmental Alternative Splicing database. MeDAS is an added-value resource that re-analyses publicly archived RNA-seq libraries to provide quantitative data on alternative splicing events as they vary across the time course of development. It has broad temporal and taxonomic scope and is intended to assist the user in identifying trends in alternative splicing throughout development. To create MeDAS, we re-analysed a curated set of 2232 Illumina polyA+ RNA-seq libraries that chart detailed time courses of embryonic and post-natal development across 18 species with a taxonomic range spanning the major metazoan lineages from Caenorhabditis elegans to human. MeDAS is freely available at https://das.chenlulab.com both as raw data tables and as an interactive browser allowing searches by species, tissue, or genomic feature (gene, transcript or exon ID and sequence). Results will provide details on alternative splicing events identified for the queried feature and can be visualised at the gene-, transcript- and exon-level as time courses of expression and inclusion levels, respectively.

Accurately identifying single-nucleotide polymorphisms (SNPs) from bacterial sequencing data is an essential requirement for using genomics to track transmission and predict important phenotypes such as antimicrobial resistance. However, most previous performance evaluations of SNP calling have been restricted to eukaryotic (human) data. Additionally, bacterial SNP calling requires choosing an appropriate reference genome to align reads to, which, together with the bioinformatic pipeline, affects the accuracy and completeness of a set of SNP calls obtained. This study evaluates the performance of 209 SNP-calling pipelines using a combination of simulated data from 254 strains of 10 clinically common bacteria and real data from environmentally sourced and genomically diverse isolates within the genera Citrobacter, Enterobacter, Escherichia, and Klebsiella.

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

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.

The laboratory rat is widely used as a model for human diseases. Many of these diseases involve monocytes and tissue macrophages in different states of activation. Whilst methods for in vitro differentiation of mouse macrophages from embryonic stem cells (ESC) and bone marrow (BM) are well established, these are lacking for the rat. The gene expression profiles of rat macrophages have also not been characterised to the same extent as mouse. We have established the methodology for production of rat ESC-derived macrophages and compared their gene expression profiles to macrophages obtained from the lung and peritoneal cavity and those differentiated from BM and blood monocytes. We determined the gene signature of Kupffer cells in the liver using rats deficient in macrophage colony stimulating factor receptor (CSF1R). We also examined the response of BM-derived macrophages to lipopolysaccharide (LPS). The results indicate that many, but not all, tissue-specific adaptations observed in mice are conserved in the rat. Importantly, we show that unlike mice, rat macrophages express the CSF1R ligand, colony stimulating factor 1 (CSF1).

Following the diagnosis of a paediatric disorder caused by an apparently de novo mutation, a recurrence risk of 1-2% is frequently quoted due to the possibility of parental germline mosaicism; but for any specific couple, this figure is usually incorrect. We present a systematic approach to providing individualized recurrence risk. By combining locus-specific sequencing of multiple tissues to detect occult mosaicism with long-read sequencing to determine the parent-of-origin of the mutation, we show that we can stratify the majority of couples into one of seven discrete categories associated with substantially different risks to future offspring. Among 58 families with a single affected offspring (representing 59 de novo mutations in 49 genes), the recurrence risk for 35 (59%) was decreased below 0.1%, but increased owing to parental mixed mosaicism for 5 (9%)-that could be quantified in semen for paternal cases (recurrence risks of 5.6-12.1%). Implementation of this strategy offers the prospect of driving a major transformation in the practice of genetic counselling.

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.

Disseminated cancer cells frequently lodge near vasculature in secondary organs. However, our understanding of the cellular crosstalk invoked at perivascular sites is still rudimentary. Here, we identify intercellular machinery governing formation of a pro-metastatic vascular niche during breast cancer colonization in the lung. We show that specific secreted factors, induced in metastasis-associated endothelial cells (ECs), promote metastasis in mice by enhancing stem cell properties and the viability of cancer cells. Perivascular macrophages, activated via tenascin C (TNC) stimulation of Toll-like receptor 4 (TLR4), were shown to be crucial in niche activation by secreting nitric oxide (NO) and tumor necrosis factor (TNF) to induce EC-mediated production of niche components. Notably, this mechanism was independent of vascular endothelial growth factor (VEGF), a key regulator of EC behavior and angiogenesis. However, targeting both macrophage-mediated vascular niche activation and VEGF-regulated angiogenesis resulted in added potency to curb lung metastasis in mice. Together, our findings provide mechanistic insights into the formation of vascular niches in metastasis.

Hypoxemia is a defining feature of acute respiratory distress syndrome (ARDS), an often-fatal complication of pulmonary or systemic inflammation, yet the resulting tissue hypoxia, and its impact on immune responses, is often neglected. In the present study, we have shown that ARDS patients were hypoxemic and monocytopenic within the first 48 h of ventilation. Monocytopenia was also observed in mouse models of hypoxic acute lung injury, in which hypoxemia drove the suppression of type I interferon signaling in the bone marrow. This impaired monopoiesis resulted in reduced accumulation of monocyte-derived macrophages and enhanced neutrophil-mediated inflammation in the lung. Administration of colony-stimulating factor 1 in mice with hypoxic lung injury rescued the monocytopenia, altered the phenotype of circulating monocytes, increased monocyte-derived macrophages in the lung and limited injury. Thus, tissue hypoxia altered the dynamics of the immune response to the detriment of the host and interventions to address the aberrant response offer new therapeutic strategies for ARDS.

Colony-stimulating factor-1 receptor (CSF1R) is a receptor tyrosine kinase that controls the differentiation and maintenance of most tissue-resident macrophages, and the inhibition of CSF1R has been suggested as a possible therapy for a range of human disorders. Herein, we present the synthesis, development, and structure-activity relationship of a series of highly selective pyrrolo[2,3-d]pyrimidines, showing subnanomolar enzymatic inhibition of this receptor and with excellent selectivity toward other kinases in the platelet-derived growth factor receptor (PDGFR) family. The crystal structure of the protein and 23 revealed that the binding conformation of the protein is DFG-out-like. The most promising compounds in this series were profiled for cellular potency and subjected to pharmacokinetic profiling and in vivo stability, indicating that this compound class could be relevant in a potential disease setting. Additionally, these compounds inhibited primarily the autoinhibited form of the receptor, contrasting the behavior of pexidartinib, which could explain the exquisite selectivity of these structures.

The sequencing of multiple genomes of the same plant species has revealed polymorphic gene and exon loss. Genes associated with disease resistance are overrepresented among those showing structural variations, suggesting an adaptive role for gene and exon presence-absence variation (PAV). To shed light on the possible functional relevance of polymorphic coding region loss and the mechanisms driving this process, we characterized genes that have lost entire exons or their whole coding regions in 17 fully sequenced Arabidopsis thaliana accessions. We found that although a significant enrichment in genes associated with certain functional categories is observed, PAV events are largely restricted to genes with signatures of reduced essentiality: PAV genes tend to be newer additions to the genome, tissue specific, and lowly expressed. In addition, PAV genes are located in regions of lower gene density and higher transposable element density. Partial coding region PAV events were associated with only a marginal reduction in gene expression level in the affected accession and occurred in genes with higher levels of alternative splicing in the Col-0 accession. Together, these results suggest that although adaptive scenarios cannot be ruled out, PAV events can be explained without invoking them.

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.

The origins and functions of kidney macrophages in the adult have been explored, but their roles during development remain largely unknown. Here we characterise macrophage arrival, localisation, heterogeneity, and functions during kidney organogenesis. Using genetic approaches to ablate macrophages, we identify a role for macrophages in nephron progenitor cell clearance as mouse kidney development begins. Throughout renal organogenesis, most kidney macrophages are perivascular and express F4/80 and CD206. These macrophages are enriched for mRNAs linked to developmental processes, such as blood vessel morphogenesis. Using antibody-mediated macrophage-depletion, we show macrophages support vascular anastomoses in cultured kidney explants. We also characterise a subpopulation of galectin-3+ (Gal3+) myeloid cells within the developing kidney. Our findings may stimulate research into macrophage-based therapies for renal developmental abnormalities and have implications for the generation of bioengineered kidney tissues.

CSF1 is the primary growth factor controlling macrophage numbers, but whether expression of the CSF1 receptor differs between discrete populations of mononuclear phagocytes remains unclear. We have generated a Csf1r-mApple transgenic fluorescent reporter mouse that, in combination with lineage tracing, Alexa Fluor 647-labeled CSF1-Fc and CSF1, and a modified ΔCsf1-enhanced cyan fluorescent protein (ECFP) transgene that lacks a 150 bp segment of the distal promoter, we have used to dissect the differentiation and CSF1 responsiveness of mononuclear phagocyte populations in situ. Consistent with previous Csf1r-driven reporter lines, Csf1r-mApple was expressed in blood monocytes and at higher levels in tissue macrophages, and was readily detectable in whole mounts or with multiphoton microscopy. In the liver and peritoneal cavity, uptake of labeled CSF1 largely reflected transgene expression, with greater receptor activity in mature macrophages than monocytes and tissue-specific expression in conventional dendritic cells. However, CSF1 uptake also differed between subsets of monocytes and discrete populations of tissue macrophages, which in macrophages correlated with their level of dependence on CSF1 receptor signaling for survival rather than degree of transgene expression. A double ΔCsf1r-ECFP-Csf1r-mApple transgenic mouse distinguished subpopulations of microglia in the brain, and permitted imaging of interstitial macrophages distinct from alveolar macrophages, and pulmonary monocytes and conventional dendritic cells. The Csf1r-mApple mice and fluorescently labeled CSF1 will be valuable resources for the study of macrophage and CSF1 biology, which are compatible with existing EGFP-based reporter lines.

Colony-stimulating factor 1 (CSF1) controls the growth and differentiation of macrophages.CSF1R signaling has been implicated in the maintenance of the intestinal stem cell niche and differentiation of Paneth cells, but evidence of expression of CSF1R within the crypt is equivocal. Here we show that CSF1R-dependent macrophages influence intestinal epithelial differentiation and homeostasis. In the intestinal lamina propria CSF1R mRNA expression is restricted to macrophages which are intimately associated with the crypt epithelium, and is undetectable in Paneth cells. Macrophage ablation following CSF1R blockade affects Paneth cell differentiation and leads to a reduction of Lgr5+ intestinal stem cells. The disturbances to the crypt caused by macrophage depletion adversely affect the subsequent differentiation of intestinal epithelial cell lineages. Goblet cell density is enhanced, whereas the development of M cells in Peyer's patches is impeded. We suggest that modification of the phenotype or abundance of macrophages in the gut wall alters the development of the intestinal epithelium and the ability to sample gut antigens.