An Atypical Parvovirus Drives Chronic Tubulointerstitial Nephropathy and Kidney Fibrosis.

The occurrence of a spontaneous nephropathy with intranuclear inclusions in laboratory mice has puzzled pathologists for over 4 decades, because its etiology remains elusive. The condition is more severe in immunodeficient animals, suggesting an infectious cause. Using metagenomics, we identify the causative agent as an atypical virus, termed "mouse kidney parvovirus" (MKPV), belonging to a divergent genus of Parvoviridae. MKPV was identified in animal facilities in Australia and North America, is transmitted via a fecal-oral or urinary-oral route, and is controlled by the adaptive immune system. Detailed analysis of the clinical course and histopathological features demonstrated a stepwise progression of pathology ranging from sporadic tubular inclusions to tubular degeneration and interstitial fibrosis and culminating in renal failure. In summary, we identify a widely distributed pathogen in laboratory mice and establish MKPV-induced nephropathy as a new tool for elucidating mechanisms of tubulointerstitial fibrosis that shares molecular features with chronic kidney disease in humans.

Orchestrated intron retention regulates normal granulocyte differentiation.

Intron retention (IR) is widely recognized as a consequence of mis-splicing that leads to failed excision of intronic sequences from pre-messenger RNAs. Our bioinformatic analyses of transcriptomic and proteomic data of normal white blood cell differentiation reveal IR as a physiological mechanism of gene expression control. IR regulates the expression of 86 functionally related genes, including those that determine the nuclear shape that is unique to granulocytes. Retention of introns in specific genes is associated with downregulation of splicing factors and higher GC content. IR, conserved between human and mouse, led to reduced mRNA and protein levels by triggering the nonsense-mediated decay (NMD) pathway. In contrast to the prevalent view that NMD is limited to mRNAs encoding aberrant proteins, our data establish that IR coupled with NMD is a conserved mechanism in normal granulopoiesis. Physiological IR may provide an energetically favorable level of dynamic gene expression control prior to sustained gene translation.

Modeling RNA-Binding Protein Specificity In Vivo by Precisely Registering Protein-RNA Crosslink Sites.

RNA-binding proteins (RBPs) regulate post-transcriptional gene expression by recognizing short and degenerate sequence motifs in their target transcripts, but precisely defining their binding specificity remains challenging. Crosslinking and immunoprecipitation (CLIP) allows for mapping of the exact protein-RNA crosslink sites, which frequently reside at specific positions in RBP motifs at single-nucleotide resolution. Here, we have developed a computational method, named mCross, to jointly model RBP binding specificity while precisely registering the crosslinking position in motif sites. We applied mCross to 112 RBPs using ENCODE eCLIP data and validated the reliability of the discovered motifs by genome-wide analysis of allelic binding sites. Our analyses revealed that the prototypical SR protein SRSF1 recognizes clusters of GGA half-sites in addition to its canonical GGAGGA motif. Therefore, SRSF1 regulates splicing of a much larger repertoire of transcripts than previously appreciated, including HNRNPD and HNRNPDL, which are involved in multivalent protein assemblies and phase separation.

Intron retention: importance, challenges, and opportunities.

Recent landmark discoveries have underpinned the physiological importance of intron retention (IR) across multiple domains of life and revealed an unexpected breath of functions in a large variety of biological processes. Despite significant progress in the field, some challenges remain. Once solved, opportunities will arise for discovering more functions of IR.

The multifaceted effects of YTHDC1-mediated nuclear m6A recognition.

N6-methyladenosine or m6A modification to mRNAs is now recognised as a key regulator of gene expression and protein translation. The fate of m6A-modified mRNAs is decoded by m6A readers, mostly found in the cytoplasm, except for the nuclear-localised YTHDC1. While earlier studies have implicated YTHDC1-m6A functions in alternative splicing and mRNA export, recent literature has expanded its close association to the chromatin-associated, noncoding and regulatory RNAs to fine-tune transcription and gene expression in cells. Here, we summarise current progress in the study of YTHDC1 function in cells, highlighting its multiple modes of action in regulating gene expression, and propose the formation of YTHDC1 nuclear condensates as a general mechanism that underlies its diverse functions in the nucleus.

SpliceWiz: interactive analysis and visualization of alternative splicing in R.

Alternative splicing (AS) is a crucial mechanism for regulating gene expression and isoform diversity in eukaryotes. However, the analysis and visualization of AS events from RNA sequencing data remains challenging. Most tools require a certain level of computer literacy and the available means of visualizing AS events, such as coverage and sashimi plots, have limitations and can be misleading. To address these issues, we present SpliceWiz, an R package with an interactive Shiny interface that allows easy and efficient AS analysis and visualization at scale. A novel normalization algorithm is implemented to aggregate splicing levels within sample groups, thereby allowing group differences in splicing levels to be accurately visualized. The tool also offers downstream gene ontology enrichment analysis, highlighting ASEs belonging to functional pathways of interest. SpliceWiz is optimized for speed and efficiency and introduces a new file format for coverage data storage that is more efficient than BigWig. Alignment files are processed orders of magnitude faster than other R-based AS analysis tools and on par with command-line tools. Overall, SpliceWiz streamlines AS analysis, enabling reliable identification of functionally relevant AS events for further characterization. SpliceWiz is a Bioconductor package and is also available on GitHub (https://github.com/alexchwong/SpliceWiz).

Dynamic changes in RNA m6A and 5 hmC influence gene expression programs during macrophage differentiation and polarisation.

RNA modifications are essential for the establishment of cellular identity. Although increasing evidence indicates that RNA modifications regulate the innate immune response, their role in monocyte-to-macrophage differentiation and polarisation is unclear. While m6A has been widely studied, other RNA modifications, including 5 hmC, remain poorly characterised. We profiled m6A and 5 hmC epitranscriptomes, transcriptomes, translatomes and proteomes of monocytes and macrophages at rest and pro- and anti-inflammatory states. Transcriptome-wide mapping of m6A and 5 hmC reveals enrichment of m6A and/or 5 hmC on specific categories of transcripts essential for macrophage differentiation. Our analyses indicate that m6A and 5 hmC modifications are present in transcripts with critical functions in pro- and anti-inflammatory macrophages. Notably, we also discover the co-occurrence of m6A and 5 hmC on alternatively-spliced isoforms and/or opposing ends of the untranslated regions (UTR) of mRNAs with key roles in macrophage biology. In specific examples, RNA 5 hmC controls the decay of transcripts independently of m6A. This study provides (i) a comprehensive dataset to interrogate the role of RNA modifications in a plastic system (ii) a resource for exploring different layers of gene expression regulation in the context of human monocyte-to-macrophage differentiation and polarisation, (iii) new insights into RNA modifications as central regulators of effector cells in innate immunity.

Ultrafast two-photon fluorescence imaging of cerebral blood circulation in the mouse brain in vivo.

Characterizing blood flow dynamics in vivo is critical to understanding the function of the vascular network under physiological and pathological conditions. Existing methods for hemodynamic imaging have insufficient spatial and temporal resolution to monitor blood flow at the cellular level in large blood vessels. By using an ultrafast line-scanning module based on free-space angular chirped enhanced delay, we achieved two-photon fluorescence imaging of cortical blood flow at 1,000 two-dimensional (2D) frames and 1,000,000 one-dimensional line scans per second in the awake mouse. This orders-of-magnitude increase in temporal resolution allowed us to measure cerebral blood flow at up to 49 mm/s and observe pulsatile blood flow at harmonics of heart rate. Directly visualizing red blood cell (RBC) flow through vessels down to >800 µm in depth, we characterized cortical layer­dependent flow velocity distributions of capillaries, obtained radial velocity profiles and kilohertz 2D velocity mapping of multifile blood flow, and performed RBC flux measurements from penetrating blood vessels.

The m6A-epitranscriptome in brain plasticity, learning and memory.

Activity-dependent gene expression and protein translation underlie the ability of neurons to dynamically adjust their synaptic strength in response to sensory experience and during learning. The emerging field of epitranscriptomics (RNA modifications) has rapidly shifted our views on the mechanisms that regulate gene expression. Among hundreds of biochemical modifications on RNA, N6-methyladenosine (m6A) is the most abundant reversible mRNA modification in the brain. Its dynamic nature and ability to regulate all aspects of mRNA processing have positioned m6A as an important and versatile regulator of nervous system functions, including neuronal plasticity, learning and memory. In this review, we summarise recent experimental evidence that supports the role of m6A signalling in learning and memory, as well as providing an overview of the underlying molecular mechanisms in neurons. We also discuss the consequences of perturbed m6A signalling and/or its regulatory networks which are increasingly being linked to various cognitive disorders in humans.

Association of clonal hematopoiesis and mosaic chromosomal alterations with solid malignancy incidence and mortality.

BACKGROUND: Understanding the impact of clonal hematopoiesis of indeterminate potential (CHIP) and mosaic chromosomal alterations (mCAs) on solid tumor risk and mortality can shed light on novel cancer pathways. METHODS: The authors analyzed whole genome sequencing data from the Trans-Omics for Precision Medicine Women's Health Initiative study (n = 10,866). They investigated the presence of CHIP and mCA and their association with the development and mortality of breast, lung, and colorectal cancers. RESULTS: CHIP was associated with higher risk of breast (hazard ratio [HR], 1.30; 95% confidence interval [CI], 1.03-1.64; p = .02) but not colorectal (p = .77) or lung cancer (p = .32). CHIP carriers who developed colorectal cancer also had a greater risk for advanced-stage (p = .01), but this was not seen in breast or lung cancer. CHIP was associated with increased colorectal cancer mortality both with (HR, 3.99; 95% CI, 2.41-6.62; p < .001) and without adjustment (HR, 2.50; 95% CI, 1.32-4.72; p = .004) for advanced-stage and a borderline higher breast cancer mortality (HR, 1.53; 95% CI, 0.98-2.41; p = .06). Conversely, mCA (cell fraction [CF] >3%) did not correlate with cancer risk. With higher CFs (mCA >5%), autosomal mCA was associated with increased breast cancer risk (HR, 1.39; 95% CI, 1.06-1.83; p = .01). There was no association of mCA (>3%) with breast, colorectal, or lung mortality except higher colon cancer mortality (HR, 2.19; 95% CI, 1.11-4.3; p = .02) with mCA >5%. CONCLUSIONS: CHIP and mCA (CF >5%) were associated with higher breast cancer risk and colorectal cancer mortality individually. These data could inform on novel pathways that impact cancer risk and lead to better risk stratification.

Epigenetic activation of SOX11 is associated with recurrence and progression of ductal carcinoma in situ to invasive breast cancer.

BACKGROUND: Risk of recurrence and progression of ductal carcinoma in situ (DCIS) to invasive cancer remains uncertain, emphasizing the need for developing predictive biomarkers of aggressive DCIS. METHODS: Human cell lines and mouse models of disease progression were analyzed for candidate risk predictive biomarkers identified and validated in two independent DCIS cohorts. RESULTS: RNA profiling of normal mammary and DCIS tissues (n = 48) revealed that elevated SOX11 expression correlates with MKI67, EZH2, and DCIS recurrence score. The 21T human cell line model of DCIS progression to invasive cancer and two mouse models developing mammary intraepithelial neoplasia confirmed the findings. AKT activation correlated with chromatin accessibility and EZH2 enrichment upregulating SOX11 expression. AKT and HER2 inhibitors decreased SOX11 expression along with diminished mammosphere formation. SOX11 was upregulated in HER2+ and basal-like subtypes (P < 0.001). Longitudinal DCIS cohort (n = 194) revealed shorter recurrence-free survival in SOX11+ than SOX11- patients (P = 0.0056 in all DCIS; P < 0.0001 in HER2+ subtype) associated with increased risk of ipsilateral breast event/IBE (HR = 1.9, 95%CI = 1.2-2.9; P = 0.003). DISCUSSION: Epigenetic activation of SOX11 drives recurrence of DCIS and progression to invasive cancer, suggesting SOX11 as a predictive biomarker of IBE.

Overexpression of Estrogen Receptor α in Mammary Glands of Aging Mice Is Associated with a Proliferative Risk Signature and Generation of Estrogen Receptor α-Positive Mammary Adenocarcinomas.

Age is a risk factor for human estrogen receptor-positive breast cancer, with highest prevalence following menopause. While transcriptome risk profiling is available for human breast cancers, it is not yet developed for prognostication for primary or secondary breast cancer development utilizing at-risk breast tissue. Both estrogen receptor α (ER) and aromatase overexpression have been linked to human breast cancer. Herein, conditional genetically engineered mouse models of estrogen receptor 1 (Esr1) and cytochrome P450 family 19 subfamily A member 1 (CYP19A1) were used to show that induction of Esr1 overexpression just before or with reproductive senescence and maintained through age 30 months resulted in significantly higher prevalence of estrogen receptor-positive adenocarcinomas than CYP19A1 overexpression. All adenocarcinomas tested showed high percentages of ER+ cells. Mammary cancer development was preceded by a persistent proliferative transcriptome risk signature initiated within 1 week of transgene induction that showed parallels to the Prosigna/Prediction Analysis of Microarray 50 human prognostic signature for early-stage human ER+ breast cancer. CYP19A1 mice also developed ER+ mammary cancers, but histology was more divided between adenocarcinoma and adenosquamous, with one ER- adenocarcinoma. Results demonstrate that, like humans, generation of ER+ adenocarcinoma in mice was facilitated by aging mice past the age of reproductive senescence. Esr1 overexpression was associated with a proliferative estrogen pathway-linked signature that preceded appearance of ER+ mammary adenocarcinomas.

Variation in microbial feature perception in the Rutaceae family with immune receptor conservation in citrus.

Although much is known about the responses of model plants to microbial features, we still lack an understanding of the extent of variation in immune perception across members of a plant family. In this work, we analyzed immune responses in Citrus and wild relatives, surveying 86 Rutaceae genotypes with differing leaf morphologies and disease resistances. We found that responses to microbial features vary both within and between members. Species in 2 subtribes, the Balsamocitrinae and Clauseninae, can recognize flagellin (flg22), cold shock protein (csp22), and chitin, including 1 feature from Candidatus Liberibacter species (csp22CLas), the bacterium associated with Huanglongbing. We investigated differences at the receptor level for the flagellin receptor FLAGELLIN SENSING 2 (FLS2) and the chitin receptor LYSIN MOTIF RECEPTOR KINASE 5 (LYK5) in citrus genotypes. We characterized 2 genetically linked FLS2 homologs from "Frost Lisbon" lemon (Citrus ×limon, responsive) and "Washington navel" orange (Citrus ×aurantium, nonresponsive). Surprisingly, FLS2 homologs from responsive and nonresponsive genotypes were expressed in Citrus and functional when transferred to a heterologous system. "Washington navel" orange weakly responded to chitin, whereas "Tango" mandarin (C. ×aurantium) exhibited a robust response. LYK5 alleles were identical or nearly identical between the 2 genotypes and complemented the Arabidopsis (Arabidopsis thaliana) lyk4/lyk5-2 mutant with respect to chitin perception. Collectively, our data indicate that differences in chitin and flg22 perception in these citrus genotypes are not the results of sequence polymorphisms at the receptor level. These findings shed light on the diversity of perception of microbial features and highlight genotypes capable of recognizing polymorphic pathogen features.

CYP eicosanoid pathway mediates colon cancer-promoting effects of dietary linoleic acid.

Human and animal studies support that consuming a high level of linoleic acid (LA, 18:2ω-6), an essential fatty acid and key component of the human diet, increases the risk of colon cancer. However, results from human studies have been inconsistent, making it challenging to establish dietary recommendations for optimal LA intake. Given the importance of LA in the human diet, it is crucial to better understand the molecular mechanisms underlying its potential colon cancer-promoting effects. Using LC-MS/MS-based targeted lipidomics, we find that the cytochrome P450 (CYP) monooxygenase pathway is a major pathway for LA metabolism in vivo. Furthermore, CYP monooxygenase is required for the colon cancer-promoting effects of LA, since the LA-rich diet fails to exacerbate colon cancer in CYP monooxygenase-deficient mice. Finally, CYP monooxygenase mediates the pro-cancer effects of LA by converting LA to epoxy octadecenoic acids (EpOMEs), which have potent effects on promoting colon tumorigenesis via gut microbiota-dependent mechanisms. Overall, these results support that CYP monooxygenase-mediated conversion of LA to EpOMEs plays a crucial role in the health effects of LA, establishing a unique mechanistic link between dietary fatty acid intake and cancer risk. These results could help in developing more effective dietary guidelines for optimal LA intake and identifying subpopulations that may be especially vulnerable to LA's negative effects.

Overexpression of VIRMA confers vulnerability to breast cancers via the m6A-dependent regulation of unfolded protein response.

Virilizer-like m6A methyltransferase-associated protein (VIRMA) maintains the stability of the m6A writer complex. Although VIRMA is critical for RNA m6A deposition, the impact of aberrant VIRMA expression in human diseases remains unclear. We show that VIRMA is amplified and overexpressed in 15-20% of breast cancers. Of the two known VIRMA isoforms, the nuclear-enriched full-length but not the cytoplasmic-localised N-terminal VIRMA promotes m6A-dependent breast tumourigenesis in vitro and in vivo. Mechanistically, we reveal that VIRMA overexpression upregulates the m6A-modified long non-coding RNA, NEAT1, which contributes to breast cancer cell growth. We also show that VIRMA overexpression enriches m6A on transcripts that regulate the unfolded protein response (UPR) pathway but does not promote their translation to activate the UPR under optimal growth conditions. Under stressful conditions that are often present in tumour microenvironments, VIRMA-overexpressing cells display enhanced UPR and increased susceptibility to death. Our study identifies oncogenic VIRMA overexpression as a vulnerability that may be exploited for cancer therapy.

Increased chromatin accessibility facilitates intron retention in specific cell differentiation states.

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

Clinical and radiographic outcomes of shoulder hemiarthroplasty for patients with glenoid medialization.

BACKGROUND: Glenoid bone loss in shoulder arthroplasty is a difficult problem that is prone to complications because of challenges with achieving glenoid component fixation and stability. The purpose of this study was to evaluate the outcomes of primary shoulder hemiarthroplasty for patients with severe glenoid medialization precluding placement of a glenoid component. METHODS: This was a retrospective case series evaluating patients who underwent shoulder hemiarthroplasty for severe glenoid erosion and medialization between 2010 and 2020. Patients were evaluated via chart review and phone survey to determine if there were any reoperations at final follow-up and to obtain Single Assessment Numeric Evaluation (SANE), American Shoulder and Elbow Surgeons Standardized Shoulder Assessment Form (ASES), and Simple Shoulder Test (SST) scores. Preoperative and postoperative radiographs were evaluated and compared to determine glenoid morphology, version, medialization, acromiohumeral distance, and humeral offset. Final postoperative films were also evaluated for anterosuperior migration and signs of mechanical failure, including loosening or periprosthetic fracture. RESULTS: Overall, there were 28 patients during this period who underwent shoulder hemiarthroplasty for severe glenoid medialization. Eight patients were deceased at the time of the study, 2 were unable to complete surveys because of dementia, and 7 were lost to follow-up. The final cohort included 11 shoulders and 11 patients with a mean age of 71 ± 7.1 years and mean follow-up of 6.7 years (range 1.6-13.0 years). Mean postoperative SANE, ASES, and SST scores were 80.6 ± 17.6, 71.5 ± 29.3, and 7.6 ± 2.0, respectively. There were no reoperations or revision surgeries at final follow-up. Radiographic evaluation demonstrated severe glenoid medialization and decreased lateral humeral offset, which was unchanged postoperatively. There were 2 patients with signs of anterosuperior migration at final radiographic follow-up but no signs of implant failure. CONCLUSION: Shoulder hemiarthroplasty for severe medial glenoid bone loss provides modest clinical outcomes and low rates of reoperation at mid- to long-term follow-up and is an option worth considering in cases where placement of a glenoid component is challenging because of deficient bone stock and high risk for complications.

Chromosomal imbalances detected via RNA-sequencing in 28 cancers.

MOTIVATION: RNA-sequencing (RNA-seq) of tumor tissue is typically only used to measure gene expression. Here, we present a statistical approach that leverages existing RNA-seq data to also detect somatic copy number alterations (SCNAs), a pervasive phenomenon in human cancers, without a need to sequence the corresponding DNA. RESULTS: We present an analysis of 4942 participant samples from 28 cancers in The Cancer Genome Atlas (TCGA), demonstrating robust detection of SCNAs from RNA-seq. Using genotype imputation and haplotype information, our RNA-based method had a median sensitivity of 85% to detect SCNAs defined by DNA analysis, at high specificity (∼95%). As an example of translational potential, we successfully replicated SCNA features associated with breast cancer subtypes. Our results credential haplotype-based inference based on RNA-seq to detect SCNAs in clinical and population-based settings. AVAILABILITY AND IMPLEMENTATION: The analyses presented use the data publicly available from TCGA Research Network (http://cancergenome.nih.gov/). See Methods for details regarding data downloads. hapLOHseq software is freely available under The MIT license and can be downloaded from http://scheet.org/software.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MLH1-methylated endometrial cancer under 60 years of age as the "sentinel" cancer in female carriers of high-risk constitutional MLH1 epimutation.

OBJECTIVE: Universal screening of endometrial carcinoma (EC) for mismatch repair deficiency (MMRd) and Lynch syndrome uses presence of MLH1 methylation to omit common sporadic cases from follow-up germline testing. However, this overlooks rare cases with high-risk constitutional MLH1 methylation (epimutation), a poorly-recognized mechanism that predisposes to Lynch-type cancers with MLH1 methylation. We aimed to determine the role and frequency of constitutional MLH1 methylation among EC cases with MMRd, MLH1-methylated tumors. METHODS: We screened blood for constitutional MLH1 methylation using pyrosequencing and real-time methylation-specific PCR in patients with MMRd, MLH1-methylated EC ascertained from (i) cancer clinics (n = 4, <60 years), and (ii) two population-based cohorts; "Columbus-area" (n = 68, all ages) and "Ohio Colorectal Cancer Prevention Initiative (OCCPI)" (n = 24, <60 years). RESULTS: Constitutional MLH1 methylation was identified in three out of four patients diagnosed between 36 and 59 years from cancer clinics. Two had mono-/hemiallelic epimutation (∼50% alleles methylated). One with multiple primaries had low-level mosaicism in normal tissues and somatic "second-hits" affecting the unmethylated allele in all tumors, demonstrating causation. In the population-based cohorts, all 68 cases from the Columbus-area cohort were negative and low-level mosaic constitutional MLH1 methylation was identified in one patient aged 36 years out of 24 from the OCCPI cohort, representing one of six (∼17%) patients <50 years and one of 45 patients (∼2%) <60 years in the combined cohorts. EC was the first/dual-first cancer in three patients with underlying constitutional MLH1 methylation. CONCLUSIONS: A correct diagnosis at first presentation of cancer is important as it will significantly alter clinical management. Screening for constitutional MLH1 methylation is warranted in patients with early-onset EC or synchronous/metachronous tumors (any age) displaying MLH1 methylation.

Dynamic intron retention modulates gene expression in the monocytic differentiation pathway.

Monocytes play a crucial role in maintaining homeostasis and mediating a successful innate immune response. They also act as central players in diverse pathological conditions, thus making them an attractive therapeutic target. Within the bone marrow, monocytes arise from a committed precursor termed Common Monocyte Progenitor (cMoP). However, molecular mechanisms that regulate the differentiation of cMoP to various monocytic subsets remain unclear. Herein, we purified murine myeloid precursors for deep poly-A-enriched RNA sequencing to understand the role of alternative splicing in the development and differentiation of monocytes under homeostasis. Our analyses revealed intron retention to be the major alternative splicing mechanism involved in the monocyte differentiation cascade, especially in the differentiation of Ly6Chi monocytes to Ly6Clo monocytes. Furthermore, we found that the intron retention of key genes involved in the differentiation of murine Ly6Chi to Ly6Clo monocytes was also conserved in humans. Our data highlight the unique role of intron retention in the regulation of the monocytic differentiation pathway.