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

Target repositioning using multi-layer networks and machine learning: The case of prostate cancer

The discovery of novel therapeutic targets, defined as proteins which drugs can interact with to induce therapeutic benefits, typically represent the first and most important step of drug discovery. One solution for target discovery is target repositioning, a strategy which relies on the repurposing of known targets for new diseases, leading to new treatments, less side effects and potential drug synergies. Biological networks have emerged as powerful tools for integrating heterogeneous data and facilitating the prediction of biological or therapeutic properties. Consequently, they are widely employed to predict new therapeutic targets by characterizing potential candidates, often based on their interactions within a Protein-Protein Interaction (PPI) network, and their proximity to genes associated with the disease. However, over-reliance on PPI networks and the assumption that potential targets are necessarily near known genes can introduce biases that may limit the effectiveness of these methods. This study addresses these limitations in two ways. First, by exploiting a multi-layer network which incorporates additional information such as gene regulation, metabolite interactions, metabolic pathways, and several disease signatures such as Differentially Expressed Genes, mutated genes, Copy Number Alteration, and structural variants. Second, by extracting relevant features from the network using several approaches including proximity to disease-associated genes, but also unbiased approaches such as propagation-based methods, topological metrics, and module detection algorithms. Using prostate cancer as a case study, the best features were identified and utilized to train machine learning algorithms to predict 5 novel promising therapeutic targets for prostate cancer: IGF2R, C5AR, RAB7, SETD2 and NPBWR1.

Deciphering the molecular landscape of the FAM72 gene family: Implications for stem cell biology and cancer

Family with sequence similarity 72 (FAM72) is a protein-coding gene family located on chromosome 1 in humans, uniquely featuring four paralogs: FAM72A, FAM72B, FAM72C, and FAM72D. While FAM72's presence as a gene pair with the SLIT-ROBO Rho GTPase-activating protein 2 (SRGAP2) is intriguing, its functional roles, particularly in neural stem cells, remain incompletely understood. This review explores the distinct characteristics of FAM72, shedding light on its expression patterns, potential roles in cell cycle regulation, stem cell renewal and implications in neurogenesis and tumorigenesis.

Genome-Wide Association Study of Conformation Traits in Brazilian Holstein Cattle

The linear conformation of animals exerts an influence on health, reproduction, production, and welfare, in addition to longevity, which directly affects the profitability of milk-producing farms. The objectives of this study were (1) to perform genome-wide association studies (GWASs) of conformation traits, namely the Rump, Feet and Legs, Mammary System, Dairy Strength, and Final Classification traits, and (2) to identify genes and related pathways involved in physiological processes associated with conformation traits in Brazilian Holstein cattle. Phenotypic and genotypic data from 2339 Holstein animals distributed across the states of Rio Grande do Sul, Paraná, São Paulo, and Minas Gerais were used. The genotypic data were obtained with a 100 K SNP marker panel. The single-step genome-wide association study (ssGWAS) method was employed in the analyses. Genes close to a significant SNP were identified in an interval of 100 kb up- and downstream using the Ensembl database available in the BioMart tool. The DAVID database was used to identify the main metabolic pathways and the STRING program was employed to create the gene regulatory network. In total, 36 significant SNPs were found on 15 chromosomes; 27 of these SNPs were linked to genes that may influence the traits studied. Fourteen genes most closely related to the studied traits were identified, as well as four genes that showed interactions in important metabolic pathways such as myogenesis, adipogenesis, and angiogenesis. Among the total genes, four were associated with myogenesis (TMOD2, TMOD3, CCND2, and CTBP2), three with angiogenesis (FGF23, FGF1, and SCG3), and four with adipogenesis and body size and development (C5H12orf4, CCND2, EMILIN1, and FGF6). These results contribute to a better understanding of the biological mechanisms underlying phenotypic variability in conformation traits in Brazilian Holstein cattle.

Analysis of the heterogeneity and complexity of murine extraorbital lacrimal gland via single-cell RNA sequencing

PURPOSE

The lacrimal gland is essential for maintaining ocular surface health and avoiding external damage by secreting an aqueous layer of the tear film. However, a healthy lacrimal gland's inventory of cell types and heterogeneity remains understudied.

METHODS

Here, 10X Genome-based single-cell RNA sequencing was used to generate an unbiased classification of cellular diversity in the extraorbital lacrimal gland (ELG) of C57BL/6J mice. From 43,850 high-quality cells, we produced an atlas of cell heterogeneity and defined cell types using classic marker genes. The possible functions of these cells were analyzed through bioinformatics analysis. Additionally, the CellChat was employed for a preliminary analysis of the cell-cell communication network in the ELG.

RESULTS

Over 37 subclasses of cells were identified, including seven types of glandular epithelial cells, three types of fibroblasts, ten types of myeloid-derived immune cells, at least eleven types of lymphoid-derived immune cells, and five types of vascular-associated cell subsets. The cell-cell communication network analysis revealed that fibroblasts and immune cells play a pivotal role in the dense intercellular communication network within the mouse ELG.

CONCLUSIONS

This study provides a comprehensive transcriptome atlas and related database of the mouse ELG.

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.

Development and function of chicken XCR1+ conventional dendritic cells

Conventional dendritic cells (cDCs) are antigen-presenting cells (APCs) that play a central role in linking innate and adaptive immunity. cDCs have been well described in a number of different mammalian species, but remain poorly characterised in the chicken. In this study, we use previously described chicken cDC specific reagents, a novel gene-edited chicken line and single-cell RNA sequencing (scRNAseq) to characterise chicken splenic cDCs. In contrast to mammals, scRNAseq analysis indicates that the chicken spleen contains a single, chemokine receptor XCR1 expressing, cDC subset. By sexual maturity the XCR1+ cDC population is the most abundant mononuclear phagocyte cell subset in the chicken spleen. scRNAseq analysis revealed substantial heterogeneity within the chicken splenic XCR1+ cDC population. Immature MHC class II (MHCII)LOW XCR1+ cDCs expressed a range of viral resistance genes. Maturation to MHCIIHIGH XCR1+ cDCs was associated with reduced expression of anti-viral gene expression and increased expression of genes related to antigen presentation via the MHCII and cross-presentation pathways. To visualise and transiently ablate chicken XCR1+ cDCs in situ, we generated XCR1-iCaspase9-RFP chickens using a CRISPR-Cas9 knockin transgenesis approach to precisely edit the XCR1 locus, replacing the XCR1 coding region with genes for a fluorescent protein (TagRFP), and inducible Caspase 9. After inducible ablation, the chicken spleen is initially repopulated by immature CD1.1+ XCR1+ cDCs. XCR1+ cDCs are abundant in the splenic red pulp, in close association with CD8+ T-cells. Knockout of XCR1 prevented this clustering of cDCs with CD8+ T-cells. Taken together these data indicate a conserved role for chicken and mammalian XCR1+ cDCs in driving CD8+ T-cells responses.

ZNF714 Supports Pro-Oncogenic Features in Lung Cancer Cells

Despite the ongoing progress in diagnosis and treatments, cancer remains a threat to more than one-third of the human population. The emerging data indicate that many Krüppel-associated box zinc finger proteins (KRAB-ZNF) belonging to a large gene family may be involved in carcinogenesis. Our previous study identified Zinc Finger Protein 714 (ZNF714), a KRAB-ZNF gene of unknown function, as being commonly overexpressed in many tumors, pointing to its hypothetical oncogenic role. Here, we harnessed The Cancer Genome Atlas (TCGA)-centered databases and performed functional studies with transcriptomic and methylomic profiling to explore ZNF714 function in cancer. Our pan-cancer analyses confirmed frequent ZNF714 overexpression in multiple tumors, possibly due to regional amplification, promoter hypomethylation, and Nuclear Transcription Factor Y Subunit Beta (NFYB) signaling. We also showed that ZNF714 expression correlates with tumor immunosuppressive features. The in vitro studies indicated that ZNF714 expression positively associates with proliferation, migration, and invasion. The transcriptomic analysis of ZNF714 knocked-down cells demonstrated deregulation of cell adhesion, migration, proliferation, apoptosis, and differentiation. Importantly, we provided evidence that ZNF714 negatively regulates the expression of several known TSGs indirectly via promoter methylation. However, as ZNF714 did not show nuclear localization in our research model, the regulatory mechanisms exerted by ZNF714 require further investigation. In conclusion, our results reveal, for the first time, that ZNF714 may support pro-oncogenic features in lung cancer cells.

Single-cell chromatin accessibility profiling reveals a self-renewing muscle satellite cell state

A balance between self-renewal and differentiation is critical for the regenerative capacity of tissue-resident stem cells. In skeletal muscle, successful regeneration requires the orchestrated activation, proliferation, and differentiation of muscle satellite cells (MuSCs) that are normally quiescent. A subset of MuSCs undergoes self-renewal to replenish the stem cell pool, but the features that identify and define self-renewing MuSCs remain to be elucidated. Here, through single-cell chromatin accessibility analysis, we reveal the self-renewal versus differentiation trajectories of MuSCs over the course of regeneration in vivo. We identify Betaglycan as a unique marker of self-renewing MuSCs that can be purified and efficiently contributes to regeneration after transplantation. We also show that SMAD4 and downstream genes are genetically required for self-renewal in vivo by restricting differentiation. Our study unveils the identity and mechanisms of self-renewing MuSCs, while providing a key resource for comprehensive analysis of muscle regeneration.

A single-cell transcriptome atlas profiles early organogenesis in human embryos

The early window of human embryogenesis is largely a black box for developmental biologists. Here we probed the cellular diversity of 4-6 week human embryos when essentially all organs are just laid out. On the basis of over 180,000 single-cell transcriptomes, we generated a comprehensive atlas of 313 clusters in 18 developmental systems, which were annotated with a collection of ontology and markers from 157 publications. Together with spatial transcriptome on embryonic sections, we characterized the molecule and spatial architecture of previously unappreciated cell types. Combined with data from other vertebrates, the rich information shed light on spatial patterning of axes, systemic temporal regulation of developmental progression and potential human-specific regulation. Our study provides a compendium of early progenitor cells of human organs, which can serve as the root of lineage analysis in organogenesis.

Isogenic Cell Lines Derived from Specific Organ Metastases Exhibit Divergent Cytogenomic Aberrations

Aneuploidy, a deviation in chromosome numbers from the normal diploid set, is now recognized as a fundamental characteristic of all cancer types and is found in 70-90% of all solid tumors. The majority of aneuploidies are generated by chromosomal instability (CIN). CIN/aneuploidy is an independent prognostic marker of cancer survival and is a cause of drug resistance. Hence, ongoing research has been directed towards the development of therapeutics aimed at targeting CIN/aneuploidy. However, there are relatively limited reports on the evolution of CIN/aneuploidies within or across metastatic lesions. In this work, we built on our previous studies using a human xenograft model system of metastatic disease in mice that is based on isogenic cell lines derived from the primary tumor and specific metastatic organs (brain, liver, lung, and spine). As such, these studies were aimed at exploring distinctions and commonalities between the karyotypes; biological processes that have been implicated in CIN; single-nucleotide polymorphisms (SNPs); losses, gains, and amplifications of chromosomal regions; and gene mutation variants across these cell lines. Substantial amounts of inter- and intra-heterogeneity were found across karyotypes, along with distinctions between SNP frequencies across each chromosome of each metastatic cell line relative the primary tumor cell line. There were disconnects between chromosomal gains or amplifications and protein levels of the genes in those regions. However, commonalities across all cell lines provide opportunities to select biological processes as druggable targets that could have efficacy against the primary tumor, as well as metastases.

In vivo and in vitro matured bovine oocytes present a distinct pattern of single-cell gene expression

Oocyte gene expression is a well controlled event that promotes gamete competence to undergo maturation, fertilization, and to support early embryo development, directly affecting reproductive outcomes. Considering that in vivo controlled ovarian stimulation or in vitro maturation (IVM) for the acquisition of mature oocytes has distinct implications for gene expression, we sought to evaluate the effects of these procedures on the expression of competence-related genes in single-cell oocytes. Healthy Nelore cows of reproductive age were synchronized to harvest in vivo matured oocytes; ovaries from slaughtered animals were used to obtain cumulus-oocyte complexes that were in vitro matured. Single-cell gene expression was performed using TaqMan Low-Density Arrays and 42 genes were evaluated. In silico analysis of protein interactions and Gene Ontology (GO) analysis was performed. Reduced gene expression was observed for 24 targets in IVM oocytes when compared with those of in vivo matured oocytes (P < 0.05). Differences ranged from 1.5-fold to 4.8-fold higher in in vivo oocytes and the BMP15 (5.28), GDF9 (6.23), NOBOX (7.25), HSPA8 (7.85) and MSX1 (11.00) showed the greatest fold increases. The strongest score of functional interactions was observed between the CDC20 and CKS2, with the differentially expressed gene CDC20 being the main marker behind GO enrichment. IVM negatively affected the expression of important genes related to oocyte competency, and showed higher expression levels in in vivo matured oocytes. In vivo controlled ovarian stimulation may be a better strategy to achieve proper oocyte competence and increase the success of assisted reproductive technologies.

Identification and analysis of C17orf53 as a prognostic signature for hepatocellular carcinoma

C17orf53 is a novel gene for DNA synthesis and homologous recombination. However, the exact role of C17orf53 in hepatocellular carcinoma (HCC) remains unclear. In this study, we analyzed it using a set of public datasets. UALCAN, Human Protein Atlas (HPA), Kaplan‒Meier Plotter, Tumor Immune Estimation Resource (TIMER), cBioPortal, GEPIA, GeneMANIA, and LinkedOmics were used. Functional analysis was conducted in SK-Hep-1 cells by using small interfering RNA (siRNA). C17orf53 was highly expressed and predicted unfavorable survival in HCC patients. Moreover, it showed positive correlations with the abundance of B cells, macrophages and dendritic cells. In addition, we identified 126 genes that were positively correlated with C17orf53 and its coeffector minichromosome maintenance 8 (MCM8). These genes were mainly enriched in the cell cycle, DNA replication and Fanconi anemia pathways. Knockdown of C17orf53 significantly inhibited the proliferation of SK-Hep-1 cells and decreased the expression of MCM8, cyclin D1 and proliferating cell nuclear antigen (PCNA). Overall, C17orf53 is a novel prognostic signature for HCC.

Pan-cancer surveys indicate cell cycle-related roles of primate-specific genes in tumors and embryonic cerebrum

BACKGROUND

Despite having been extensively studied, it remains largely unclear why humans bear a particularly high risk of cancer. The antagonistic pleiotropy hypothesis predicts that primate-specific genes (PSGs) tend to promote tumorigenesis, while the molecular atavism hypothesis predicts that PSGs involved in tumors may represent recently derived duplicates of unicellular genes. However, these predictions have not been tested.

RESULTS

By taking advantage of pan-cancer genomic data, we find the upregulation of PSGs across 13 cancer types, which is facilitated by copy-number gain and promoter hypomethylation. Meta-analyses indicate that upregulated PSGs (uPSGs) tend to promote tumorigenesis and to play cell cycle-related roles. The cell cycle-related uPSGs predominantly represent derived duplicates of unicellular genes. We prioritize 15 uPSGs and perform an in-depth analysis of one unicellular gene-derived duplicate involved in the cell cycle, DDX11. Genome-wide screening data and knockdown experiments demonstrate that DDX11 is broadly essential across cancer cell lines. Importantly, non-neutral amino acid substitution patterns and increased expression indicate that DDX11 has been under positive selection. Finally, we find that cell cycle-related uPSGs are also preferentially upregulated in the highly proliferative embryonic cerebrum.

CONCLUSIONS

Consistent with the predictions of the atavism and antagonistic pleiotropy hypotheses, primate-specific genes, especially those PSGs derived from cell cycle-related genes that emerged in unicellular ancestors, contribute to the early proliferation of the human cerebrum at the cost of hitchhiking by similarly highly proliferative cancer cells.

Intestinal single-cell atlas reveals novel lymphocytes in pigs with similarities to human cells

Lymphocytes can heavily influence intestinal health, but resolving intestinal lymphocyte function is challenging as the intestine contains a vastly heterogeneous mixture of cells. Pigs are an advantageous biomedical model, but deeper understanding of intestinal lymphocytes is warranted to improve model utility. Twenty-six cell types were identified in the porcine ileum by single-cell RNA sequencing and further compared with cells in human and murine ileum. Though general consensus of cell subsets across species was revealed, some porcine-specific lymphocyte subsets were identified. Differential tissue dissection and in situ analyses conferred spatial context, revealing similar locations of lymphocyte subsets in Peyer's patches and epithelium in pig-to-human comparisons. Like humans, activated and effector lymphocytes were abundant in the ileum but not periphery of pigs, suggesting tissue-specific and/or activation-associated gene expression. Gene signatures for peripheral and ileal innate lymphoid cells newly discovered in pigs were defined and highlighted similarities to human innate lymphoid cells. Overall, we reveal novel lymphocyte subsets in pigs and highlight utility of pigs for intestinal research applications.

Sci-Seq of Human Fetal Salivary Tissue Introduces Human Transcriptional Paradigms and a Novel Cell Population

Multiple pathologies and non-pathological factors can disrupt the function of the non-regenerative human salivary gland including cancer and cancer therapeutics, autoimmune diseases, infections, pharmaceutical side effects, and traumatic injury. Despite the wide range of pathologies, no therapeutic or regenerative approaches exist to address salivary gland loss, likely due to significant gaps in our understanding of salivary gland development. Moreover, identifying the tissue of origin when diagnosing salivary carcinomas requires an understanding of human fetal development. Using computational tools, we identify developmental branchpoints, a novel stem cell-like population, and key signaling pathways in the human developing salivary glands by analyzing our human fetal single-cell sequencing data. Trajectory and transcriptional analysis suggest that the earliest progenitors yield excretory duct and myoepithelial cells and a transitional population that will yield later ductal cell types. Importantly, this single-cell analysis revealed a previously undescribed population of stem cell-like cells that are derived from SD and expresses high levels of genes associated with stem cell-like function. We have observed these rare cells, not in a single niche location but dispersed within the developing duct at later developmental stages. Our studies introduce new human-specific developmental paradigms for the salivary gland and lay the groundwork for the development of translational human therapeutics.

Single-cell analysis of endometriosis reveals a coordinated transcriptional programme driving immunotolerance and angiogenesis across eutopic and ectopic tissues

Endometriosis is characterized by the growth of endometrial-like tissue outside the uterus. It affects many women during their reproductive age, causing years of pelvic pain and potential infertility. Its pathophysiology remains largely unknown, which limits early diagnosis and treatment. We characterized peritoneal and ovarian lesions at single-cell transcriptome resolution and compared them to matched eutopic endometrium, unaffected endometrium and organoids derived from these tissues, generating data on over 122,000 cells across 14 individuals. We spatially localized many of the cell types using imaging mass cytometry. We identify a perivascular mural cell specific to the peritoneal lesions, with dual roles in angiogenesis promotion and immune cell trafficking. We define an immunotolerant peritoneal niche, fundamental differences in eutopic endometrium and between lesion microenvironments and an unreported progenitor-like epithelial cell subpopulation. Altogether, this study provides a holistic view of the endometriosis microenvironment that represents a comprehensive cell atlas of the disease in individuals undergoing hormonal treatment, providing essential information for future therapeutics and diagnostics.

Transcription factor-driven coordination of cell cycle exit and lineage-specification in vivo during granulocytic differentiation : In memoriam Professor Niels Borregaard

Differentiation of multipotent stem cells into mature cells is fundamental for development and homeostasis of mammalian tissues, and requires the coordinated induction of lineage-specific transcriptional programs and cell cycle withdrawal. To understand the underlying regulatory mechanisms of this fundamental process, we investigated how the tissue-specific transcription factors, CEBPA and CEBPE, coordinate cell cycle exit and lineage-specification in vivo during granulocytic differentiation. We demonstrate that CEBPA promotes lineage-specification by launching an enhancer-primed differentiation program and direct activation of CEBPE expression. Subsequently, CEBPE confers promoter-driven cell cycle exit by sequential repression of MYC target gene expression at the G1/S transition and E2F-meditated G2/M gene expression, as well as by the up-regulation of Cdk1/2/4 inhibitors. Following cell cycle exit, CEBPE unleashes the CEBPA-primed differentiation program to generate mature granulocytes. These findings highlight how tissue-specific transcription factors coordinate cell cycle exit with differentiation through the use of distinct gene regulatory elements.

Identifying Immune Cell Infiltration and Effective Diagnostic Biomarkers in Lung Adenocarcinoma by Comprehensive Bioinformatics Analysis and In Vitro Study

Family with sequence similarity 72B (FAM72B) has been characterized in the regulation of neuronal development. Nevertheless, the prognostic value of FAM72B expression and its function in the immune microenvironment of lung adenocarcinoma (LUAD) currently remains elusive. In this study, by adopting bioinformatics methodology and experimental verification, we found that FAM72B was upregulated in lung cancer tissues and cell lines, and a higher FAM72B level predicted an unfavorable clinical outcome in LUAD patients. The knockdown of FAM72B significantly inhibited cell proliferation, cell migration, and induced cell apoptosis in LUAD. The receiver operating characteristic curve suggested that FAM72B had a high predictive accuracy for the outcomes of LUAD. Kyoto Encyclopedia of Genes and Genomes and Gene Set Enrichment Analyses confirmed that FAM72B-related genes were involved in cell proliferation and immune-response signaling pathway. Moreover, upregulated FAM72B expression was significantly associated with immune cell infiltration in the LUAD tumor microenvironment. Meanwhile, a potential ceRNA network was constructed to identify the lncRNA-AL360270.2/TMPO-AS1/AC125807.2/has-let-7a/7b/7c/7e/7f/FAM72B regulatory axis that regulates FAM72B overexpression in LUAD and is associated with a poor prognosis. We also confirmed that AL360270.2, TMPO-AS1, and AC125807.2 were significantly upregulated in LUAD cell lines than in human bronchial epithelial cells. In conclusion, FAM72B may serve as a novel biomarker in predicting the clinical prognosis and immune status for lung adenocarcinoma.

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.

Modeling Progression of Single Cell Populations Through the Cell Cycle as a Sequence of Switches

Cell cycle is a biological process underlying the existence and propagation of life in time and space. It has been an object for mathematical modeling for long, with several alternative mechanistic modeling principles suggested, describing in more or less details the known molecular mechanisms. Recently, cell cycle has been investigated at single cell level in snapshots of unsynchronized cell populations, exploiting the new methods for transcriptomic and proteomic molecular profiling. This raises a need for simplified semi-phenomenological cell cycle models, in order to formalize the processes underlying the cell cycle, at a higher abstracted level. Here we suggest a modeling framework, recapitulating the most important properties of the cell cycle as a limit trajectory of a dynamical process characterized by several internal states with switches between them. In the simplest form, this leads to a limit cycle trajectory, composed by linear segments in logarithmic coordinates describing some extensive (depending on system size) cell properties. We prove a theorem connecting the effective embedding dimensionality of the cell cycle trajectory with the number of its linear segments. We also develop a simplified kinetic model with piecewise-constant kinetic rates describing the dynamics of lumps of genes involved in S-phase and G2/M phases. We show how the developed cell cycle models can be applied to analyze the available single cell datasets and simulate certain properties of the observed cell cycle trajectories. Based on our model, we can predict with good accuracy the cell line doubling time from the length of cell cycle trajectory.

WNT/beta-catenin signalling interrupts a senescence-induction cascade in human mesenchymal stem cells that restricts their expansion

Senescence, the irreversible cell cycle arrest of damaged cells, is accompanied by a deleterious pro-inflammatory senescence-associated secretory phenotype (SASP). Senescence and the SASP are major factors in aging, cancer, and degenerative diseases, and interfere with the expansion of adult cells in vitro, yet little is known about how to counteract their induction and deleterious effects. Paracrine signals are increasingly recognized as important senescence triggers and understanding their regulation and mode of action may provide novel opportunities to reduce senescence-induced inflammation and improve cell-based therapies. Here, we show that the signalling protein WNT3A counteracts the induction of paracrine senescence in cultured human adult mesenchymal stem cells (MSCs). We find that entry into senescence in a small subpopulation of MSCs triggers a secretome that causes a feed-forward signalling cascade that with increasing speed induces healthy cells into senescence. WNT signals interrupt this cascade by repressing cytokines that mediate this induction of senescence. Inhibition of those mediators by interference with NF-κB or interleukin 6 signalling reduced paracrine senescence in absence of WNT3A and promoted the expansion of MSCs. Our work reveals how WNT signals can antagonize senescence and has relevance not only for expansion of adult cells but can also provide new insights into senescence-associated inflammatory and degenerative diseases.

Comparative transcriptional profiling of the early host response to infection by typhoidal and non-typhoidal Salmonella serovars in human intestinal organoids

Salmonella enterica represents over 2500 serovars associated with a wide-ranging spectrum of disease; from self-limiting gastroenteritis to invasive infections caused by non-typhoidal serovars (NTS) and typhoidal serovars, respectively. Host factors strongly influence infection outcome as malnourished or immunocompromised individuals can develop invasive infections from NTS, however, comparative analyses of serovar-specific host responses have been constrained by reliance on limited model systems. Here we used human intestinal organoids (HIOs), a three-dimensional “gut-like” in vitro system derived from human embryonic stem cells, to elucidate similarities and differences in host responses to NTS and typhoidal serovars. HIOs discriminated between the two most prevalent NTS, Salmonella enterica serovar Typhimurium (STM) and Salmonella enterica serovar Enteritidis (SE), and typhoidal serovar Salmonella enterica serovar Typhi (ST) in epithelial cell invasion, replication and transcriptional responses. Pro-inflammatory signaling and cytokine output was reduced in ST-infected HIOs compared to NTS infections, consistent with early stages of NTS and typhoidal diseases. While we predicted that ST would induce a distinct transcriptional profile from the NTS strains, more nuanced expression profiles emerged. Notably, pathways involved in cell cycle, metabolism and mitochondrial functions were downregulated in STM-infected HIOs and upregulated in SE-infected HIOs. These results correlated with suppression of cellular proliferation and induction of host cell death in STM-infected HIOs and in contrast, elevated levels of reactive oxygen species production in SE-infected HIOs. Collectively, these results suggest that the HIO model is well suited to reveal host transcriptional programming specific to infection by individual Salmonella serovars, and that individual NTS may provoke unique host epithelial responses during intestinal stages of infection.

Salmonella enterica is the major causative agent of bacterial infections associated with contaminated food and water. Salmonella enterica consists of over 2500 serovars of which Typhimurium (STM), Enteritidis (SE) and Typhi (ST) are the three major serovars with medical relevance to humans. These serovars elicit distinctive immune responses and cause different diseases in humans, including self-limiting diarrhea, gastroenteritis and typhoid fever. Differences in the human host response to these serovars are likely to be a major contributing factor to distinct disease outcomes but are not well characterized, possibly due to the limitations of human-derived physiological infection models. Distinct from immortalized epithelial cell culture models, human intestinal organoids (HIOs) are three-dimensional structures derived from embryonic stem cells that differentiate into intestinal mesenchymal and epithelial cells, mirroring key organizational aspects of the intestine. In this study, we used HIOs to monitor transcriptional changes during early stages of STM, SE and ST infection. Our comparative analysis showed that HIO inflammatory responses are the dominant response in all infections, but ST infection induces the weakest upregulation of inflammatory mediators relative to the other serovars. In addition, we identified several cellular processes, including cell cycle and mitochondrial functions, that were inversely regulated between STM and SE infection despite these serovars causing similar localized intestinal infection in humans. Our findings reinforce HIOs as an emerging model system to study Salmonella serovar infection and define global host transcriptional response profiles as a foundation for understanding human infection outcomes.

Ferric ammonium citrate (FAC)-induced inhibition of osteoblast proliferation/differentiation and its reversal by soybean-derived peptides (SDP)

Ferric citrate has been used to treat hyperphosphatemia, a prevalent symptom in patients with chronic kidney disease while ferric ammonium citrate (FAC), a more dissolvable format, is widely used as food additive. However, excess iron is associated with osteoporosis. Dietary soybean products have been shown to prevent the progression of osteoporosis. In this study, a group of peptides, referred as P3, was identified from the enzymolysis of soybean protein isolates, and its biological functions were investigated. The results showed that MC3T3-E1 cell cycle progression from G0/G1 to S phase was accelerated by P3 treatment. MC3T3-E1 cell proliferation was enhanced by P3 via ERK1/2 activation. Importantly, P3 treatment abolished the antiproliferative effect of FAC on MC3T3-E1 cell. In addition, P3 treatment increased the expression of ALP, COL-1, OCN, consequently promoting the differentiation and mineralization of MC3T3-E1 cells via activation of p38 MAPK pathway. Consequently, P3 treatment was able to reverse the inhibitory effect of FAC on osteoblasts differentiation and mineralization. Our findings suggest P3, as a dietary supplement, has a potential therapeutic function to attenuate the adverse effects of FAC on bone metabolism and to prevent osteoporosis progression.

Reference Transcriptomes of Porcine Peripheral Immune Cells Created Through Bulk and Single-Cell RNA Sequencing

Pigs are a valuable human biomedical model and an important protein source supporting global food security. The transcriptomes of peripheral blood immune cells in pigs were defined at the bulk cell-type and single cell levels. First, eight cell types were isolated in bulk from peripheral blood mononuclear cells (PBMCs) by cell sorting, representing Myeloid, NK cells and specific populations of T and B-cells. Transcriptomes for each bulk population of cells were generated by RNA-seq with 10,974 expressed genes detected. Pairwise comparisons between cell types revealed specific expression, while enrichment analysis identified 1,885 to 3,591 significantly enriched genes across all 8 cell types. Gene Ontology analysis for the top 25% of significantly enriched genes (SEG) showed high enrichment of biological processes related to the nature of each cell type. Comparison of gene expression indicated highly significant correlations between pig cells and corresponding human PBMC bulk RNA-seq data available in Haemopedia. Second, higher resolution of distinct cell populations was obtained by single-cell RNA-sequencing (scRNA-seq) of PBMC. Seven PBMC samples were partitioned and sequenced that produced 28,810 single cell transcriptomes distributed across 36 clusters and classified into 13 general cell types including plasmacytoid dendritic cells (DC), conventional DCs, monocytes, B-cell, conventional CD4 and CD8 αβ T-cells, NK cells, and γδ T-cells. Signature gene sets from the human Haemopedia data were assessed for relative enrichment in genes expressed in pig cells and integration of pig scRNA-seq with a public human scRNA-seq dataset provided further validation for similarity between human and pig data. The sorted porcine bulk RNAseq dataset informed classification of scRNA-seq PBMC populations; specifically, an integration of the datasets showed that the pig bulk RNAseq data helped define the CD4CD8 double-positive T-cell populations in the scRNA-seq data. Overall, the data provides deep and well-validated transcriptomic data from sorted PBMC populations and the first single-cell transcriptomic data for porcine PBMCs. This resource will be invaluable for annotation of pig genes controlling immunogenetic traits as part of the porcine Functional Annotation of Animal Genomes (FAANG) project, as well as further study of, and development of new reagents for, porcine immunology.

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.

Cyclin D1 controls development of cerebellar granule cell progenitors through phosphorylation and stabilization of ATOH1

During development, neural progenitors are in proliferative and immature states; however, the molecular machinery that cooperatively controls both states remains elusive. Here, we report that cyclin D1 (CCND1) directly regulates both proliferative and immature states of cerebellar granule cell progenitors (GCPs). CCND1 not only accelerates cell cycle but also upregulates ATOH1 protein, an essential transcription factor that maintains GCPs in an immature state. In cooperation with CDK4, CCND1 directly phosphorylates S309 of ATOH1, which inhibits additional phosphorylation at S328 and consequently prevents S328 phosphorylation-dependent ATOH1 degradation. Additionally, PROX1 downregulates Ccnd1 expression by histone deacetylation of Ccnd1 promoter in GCPs, leading to cell cycle exit and differentiation. Moreover, WNT signaling upregulates PROX1 expression in GCPs. These findings suggest that WNT-PROX1-CCND1-ATOH1 signaling cascade cooperatively controls proliferative and immature states of GCPs. We revealed that the expression and phosphorylation levels of these molecules dynamically change during cerebellar development, which are suggested to determine appropriate differentiation rates from GCPs to GCs at distinct developmental stages. This study contributes to understanding the regulatory mechanism of GCPs as well as neural progenitors.

Somatostatin-expressing parafacial neurons are CO2/H+ sensitive and regulate baseline breathing

Glutamatergic neurons in the retrotrapezoid nucleus (RTN) function as respiratory chemoreceptors by regulating breathing in response to tissue CO2/H+. The RTN and greater parafacial region may also function as a chemosensing network composed of CO2/H+-sensitive excitatory and inhibitory synaptic interactions. In the context of disease, we showed that loss of inhibitory neural activity in a mouse model of Dravet syndrome disinhibited RTN chemoreceptors and destabilized breathing (Kuo et al., 2019). Despite this, contributions of parafacial inhibitory neurons to control of breathing are unknown, and synaptic properties of RTN neurons have not been characterized. Here, we show the parafacial region contains a limited diversity of inhibitory neurons including somatostatin (Sst)-, parvalbumin (Pvalb)-, and cholecystokinin (Cck)-expressing neurons. Of these, Sst-expressing interneurons appear uniquely inhibited by CO2/H+. We also show RTN chemoreceptors receive inhibitory input that is withdrawn in a CO2/H+-dependent manner, and chemogenetic suppression of Sst+ parafacial neurons, but not Pvalb+ or Cck+ neurons, increases baseline breathing. These results suggest Sst-expressing parafacial neurons contribute to RTN chemoreception and respiratory activity.

Induced pluripotent stem cells from subjects with Lesch-Nyhan disease

Lesch-Nyhan disease (LND) is an inherited disorder caused by pathogenic variants in the HPRT1 gene, which encodes the purine recycling enzyme hypoxanthine-guanine phosphoribosyltransferase (HGprt). We generated 6 induced pluripotent stem cell (iPSC) lines from 3 individuals with LND, along with 6 control lines from 3 normal individuals. All 12 lines had the characteristics of pluripotent stem cells, as assessed by immunostaining for pluripotency markers, expression of pluripotency genes, and differentiation into the 3 primary germ cell layers. Gene expression profiling with RNAseq demonstrated significant heterogeneity among the lines. Despite this heterogeneity, several anticipated abnormalities were readily detectable across all LND lines, including reduced HPRT1 mRNA. Several unexpected abnormalities were also consistently detectable across the LND lines, including decreases in FAR2P1 and increases in RNF39. Shotgun proteomics also demonstrated several expected abnormalities in the LND lines, such as absence of HGprt protein. The proteomics study also revealed several unexpected abnormalities across the LND lines, including increases in GNAO1 decreases in NSE4A. There was a good but partial correlation between abnormalities revealed by the RNAseq and proteomics methods. Finally, functional studies demonstrated LND lines had no HGprt enzyme activity and resistance to the toxic pro-drug 6-thioguanine. Intracellular purines in the LND lines were normal, but they did not recycle hypoxanthine. These cells provide a novel resource to reveal insights into the relevance of heterogeneity among iPSC lines and applications for modeling LND.

FAM72, Glioblastoma Multiforme (GBM) and Beyond

Simple Summary

Glioblastoma multiforme (GBM) is a serious and aggressive cancer disease that has not allowed scientists to rest for decades. In this review, we consider the new gene pair |-SRGAP2–FAM72-| and discuss its role in the cell cycle and the possibility of defining new therapeutic approaches for the treatment of GBM and other cancers via this gene pair |-SRGAP2–FAM72-|.

Abstract

Neural stem cells (NSCs) offer great potential for regenerative medicine due to their excellent ability to differentiate into various specialized cell types of the brain. In the central nervous system (CNS), NSC renewal and differentiation are under strict control by the regulation of the pivotal SLIT-ROBO Rho GTPase activating protein 2 (SRGAP2)—Family with sequence similarity 72 (FAM72) master gene (i.e., |-SRGAP2–FAM72-|) via a divergent gene transcription activation mechanism. If the gene transcription control unit (i.e., the intergenic region of the two sub-gene units, SRGAP2 and FAM72) gets out of control, NSCs may transform into cancer stem cells and generate brain tumor cells responsible for brain cancer such as glioblastoma multiforme (GBM). Here, we discuss the surveillance of this |-SRGAP2–FAM72-| master gene and its role in GBM, and also in light of FAM72 for diagnosing various types of cancers outside of the CNS.

Transcriptomic Analysis of Rat Macrophages

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

HEDGEHOG/GLI Modulates the PRR11-SKA2 Bidirectional Transcription Unit in Lung Squamous Cell Carcinomas

We previously demonstrated that proline-rich protein 11 (PRR11) and spindle and kinetochore associated 2 (SKA2) constituted a head-to-head gene pair driven by a prototypical bidirectional promoter. This gene pair synergistically promoted the development of non-small cell lung cancer. However, the signaling pathways leading to the ectopic expression of this gene pair remains obscure. In the present study, we first analyzed the lung squamous cell carcinoma (LSCC) relevant RNA sequencing data from The Cancer Genome Atlas (TCGA) database using the correlation analysis of gene expression and gene set enrichment analysis (GSEA), which revealed that the PRR11-SKA2 correlated gene list highly resembled the Hedgehog (Hh) pathway activation-related gene set. Subsequently, GLI1/2 inhibitor GANT-61 or GLI1/2-siRNA inhibited the Hh pathway of LSCC cells, concomitantly decreasing the expression levels of PRR11 and SKA2. Furthermore, the mRNA expression profile of LSCC cells treated with GANT-61 was detected using RNA sequencing, displaying 397 differentially expressed genes (203 upregulated genes and 194 downregulated genes). Out of them, one gene set, including BIRC5, NCAPG, CCNB2, and BUB1, was involved in cell division and interacted with both PRR11 and SKA2. These genes were verified as the downregulated genes via RT-PCR and their high expression significantly correlated with the shorter overall survival of LSCC patients. Taken together, our results indicate that GLI1/2 mediates the expression of the PRR11-SKA2-centric gene set that serves as an unfavorable prognostic indicator for LSCC patients, potentializing new combinatorial diagnostic and therapeutic strategies in LSCC.

Global Gene Expression of Cultured Human Dermal Fibroblasts: Focus on Cell Cycle and Proliferation Status in Improving the Condition of Face Skin

Chronological skin ageing is an inevitable physiological process that results in thin and sagging skin, fine wrinkles, and gradual dermal atrophy. The main therapeutic approaches to soft tissue augmentation involve using dermal fillers, where natural fillers, such as autologous fibroblasts, are involved in generating dermal matrix proteins. The aim of this study was to determine the global transcriptome profile of three passages of dermal autologous fibroblasts from a male volunteer, focusing on the processes of the cell cycle and cell proliferation status to estimate the optimal passage of the tested cells with respect to their reimplantation. We performed K-means clustering and validation of the expression of the selected mRNA by qRT-PCR. Ten genes were selected (ANLN, BUB1, CDC20, CCNA2, DLGAP5, MKI67, PLK1, PRC1, SPAG5, and TPX2) from the top five processes annotated to cluster 5. Detailed microarray analysis of the fibroblast genes indicated that the cell population of the third passage exhibited the highest number of upregulated genes involved in the cell cycle and cell proliferation. In all cases, the results of qRT-PCR confirmed the differences in expression of the selected mRNAs between fibroblasts from the primary culture (C0) and from the first (C1), second (C2), and third (C3) cell passage. Our results thus suggest that these cells might be useful for increasing fibroblast numbers after reimplantation into a recipient's skin, and the method used in this study seems to be an excellent tool for autologous transplantation allowing the rejuvenation of aging skin.

Functional reciprocity of proteins involved in mitosis and endocytosis

Mitosis and endocytosis are two fundamental cellular processes essential for maintaining a eukaryotic life. Mitosis partitions duplicated chromatin enveloped in the nuclear membrane into two new cells, whereas endocytosis takes in extracellular substances through membrane invagination. These two processes are spatiotemporally separated and seemingly unrelated. However, recent studies have uncovered that endocytic proteins have moonlighting functions in mitosis, and mitotic complexes manifest additional roles in endocytosis. In this review, we summarize important proteins or protein complexes that participate in both processes, compare their mechanism of action, and discuss the rationale behind this multifunctionality. We also speculate on the possible origin of the functional reciprocity from an evolutionary perspective.

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.

Impact of genetic polymorphisms in kinetochore and spindle assembly genes on chromosomal aberration frequency in healthy humans

Genomic instability is a characteristic of a majority of human malignancies. Chromosomal instability is a common form of genomic instability that can be caused by defects in mitotic checkpoint genes. Chromosomal aberrations in peripheral blood are also indicative of genotoxic exposure and potential cancer risk. We evaluated associations between inherited genetic variants in 33 mitotic checkpoint genes and the frequency of chromosomal aberrations (CAs) in the presence and absence of environmental genotoxic exposure. Associations with both chromosome and chromatid type of aberrations were evaluated in two cohorts of healthy individuals, namely an exposed and a reference group consisting of 607 and 866 individuals, respectively. Binary logistic and linear regression analyses were performed for the association studies. Bonferroni-corrected significant p-value was 5 × 10^-4 for 99 tests based on the number of analyzed genes and phenotypes. In the reference group the most prominent associations were found with variants in CCNB1, a master regulator of mitosis, and in genes involved in kinetochore function, including CENPH and TEX14, whereas in the exposed group the main association was found with variants in TTK, also an important gene in kinetochore function. How the identified variants may affect the fidelity of mitotic checkpoint remains to be investigated, however, the present study suggests that genetic variation may partly explain interindividual variation in the formation of CAs.

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.

FRET-Based Sorting of Live Cells Reveals Shifted Balance between PLK1 and CDK1 Activities During Checkpoint Recovery

Cells recovering from the G2/M DNA damage checkpoint rely more on Aurora A-PLK1 signaling than cells progressing through an unperturbed G2 phase, but the reason for this discrepancy is not known. Here, we devised a method based on a FRET reporter for PLK1 activity to sort cells in distinct populations within G2 phase. We employed mass spectroscopy to characterize changes in protein levels through an unperturbed G2 phase and validated that ATAD2 levels decrease in a proteasome-dependent manner. Comparing unperturbed cells with cells recovering from DNA damage, we note that at similar PLK1 activities, recovering cells contain higher levels of Cyclin B1 and increased phosphorylation of CDK1 targets. The increased Cyclin B1 levels are due to continuous Cyclin B1 production during a DNA damage response and are sustained until mitosis. Whereas partial inhibition of PLK1 suppresses mitotic entry more efficiently when cells recover from a checkpoint, partial inhibition of CDK1 suppresses mitotic entry more efficiently in unperturbed cells. Our findings provide a resource for proteome changes during G2 phase, show that the mitotic entry network is rewired during a DNA damage response, and suggest that the bottleneck for mitotic entry shifts from CDK1 to PLK1 after DNA damage.

C17orf53 is identified as a novel gene involved in inter-strand crosslink repair

Ataxia Telangiectasia and Rad3-Related kinase (ATR) is a master regulator of genome maintenance, and participates in DNA replication and various DNA repair pathways. In a genome-wide screen for ATR-dependent fitness genes, we identified a previously uncharacterized gene, C17orf53, whose loss led to hypersensitivity to ATR inhibition. C17orf53 is conserved in vertebrates and is required for efficient cell proliferation. Loss of C17orf53 slowed down DNA replication and led to pronounced interstrand crosslink (ICL) repair defect. We showed that C17orf53 is a ssDNA- and RPA-binding protein and both characteristics are important for its functions in the cell. In addition, using multiple omics methods, we found that C17orf53 works with MCM8/9 to promote cell survival in response to ICL lesions. Taken together, our data suggest that C17orf53 is a novel component involved in ICL repair pathway.

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.

The transcriptional signature associated with human motile cilia

Cilia are complex microtubule-based organelles essential to a range of processes associated with embryogenesis and tissue homeostasis. Mutations in components of these organelles or those involved in their assembly may result in a diverse set of diseases collectively known as ciliopathies. Accordingly, many cilia-associated proteins have been described, while those distinguishing cilia subtypes are poorly defined. Here we set out to define genes associated with motile cilia in humans based on their transcriptional signature. To define the signature, we performed network deconvolution of transcriptomics data derived from tissues possessing motile ciliated cell populations. For each tissue, genes coexpressed with the motile cilia-associated transcriptional factor, FOXJ1, were identified. The consensus across tissues provided a transcriptional signature of 248 genes. To validate these, we examined the literature, databases (CilDB, CentrosomeDB, CiliaCarta and SysCilia), single cell RNA-Seq data, and the localisation of mRNA and proteins in motile ciliated cells. In the case of six poorly characterised signature genes, we performed new localisation experiments on ARMC3, EFCAB6, FAM183A, MYCBPAP, RIBC2 and VWA3A. In summary, we report a set of motile cilia-associated genes that helps shape our understanding of these complex cellular organelles.

Chromosomal instability in untreated primary prostate cancer as an indicator of metastatic potential

Background

Metastatic prostate cancer (PC) is highly lethal. The ability to identify primary tumors capable of dissemination is an unmet need in the quest to understand lethal biology and improve patient outcomes. Previous studies have linked chromosomal instability (CIN), which generates aneuploidy following chromosomal missegregation during mitosis, to PC progression. Evidence of CIN includes broad copy number alterations (CNAs) spanning > 300 base pairs of DNA, which may also be measured via RNA expression signatures associated with CNA frequency. Signatures of CIN in metastatic PC, however, have not been interrogated or well defined. We examined a published 70-gene CIN signature (CIN70) in untreated and castration-resistant prostate cancer (CRPC) cohorts from The Cancer Genome Atlas (TCGA) and previously published reports. We also performed transcriptome and CNA analysis in a unique cohort of untreated primary tumors collected from diagnostic prostate needle biopsies (PNBX) of localized (M0) and metastatic (M1) cases to determine if CIN was linked to clinical stage and outcome.

Methods

PNBX were collected from 99 patients treated in the VA Greater Los Angeles (GLA-VA) Healthcare System between 2000 and 2016. Total RNA was extracted from high-grade cancer areas in PNBX cores, followed by RNA sequencing and/or copy number analysis using OncoScan. Multivariate logistic regression analyses permitted calculation of odds ratios for CIN status (high versus low) in an expanded GLA-VA PNBX cohort (n = 121).

Results

The CIN70 signature was significantly enriched in primary tumors and CRPC metastases from M1 PC cases. An intersection of gene signatures comprised of differentially expressed genes (DEGs) generated through comparison of M1 versus M0 PNBX and primary CRPC tumors versus metastases revealed a 157-gene “metastasis” signature that was further distilled to 7-genes (PC-CIN) regulating centrosomes, chromosomal segregation, and mitotic spindle assembly. High PC-CIN scores correlated with CRPC, PC-death and all-cause mortality in the expanded GLA-VA PNBX cohort. Interestingly, approximately 1/3 of M1 PNBX cases exhibited low CIN, illuminating differential pathways of lethal PC progression.

Conclusions

Measuring CIN in PNBX by transcriptome profiling is feasible, and the PC-CIN signature may identify patients with a high risk of lethal progression at the time of diagnosis.

The hydroxymethylome of multiple myeloma identifies FAM72D as a 1q21 marker linked to proliferation

Cell identity relies on the cross-talk between genetics and epigenetics and their impact on gene expression. Oxidation of 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC) is the first step of an active DNA demethylation process occurring mainly at enhancers and gene bodies and, as such, participates in processes governing cell identity in normal and pathological conditions. Although genetic alterations are well documented in multiple myeloma (MM), epigenetic alterations associated with this disease have not yet been thoroughly analyzed. To gain insight into the biology of MM, genome-wide 5hmC profiles were obtained and showed that regions enriched in this modified base overlap with MM enhancers and super enhancers and are close to highly expressed genes. Through the definition of a MM-specific 5hmC signature, we identified FAM72D as a poor prognostic gene located on 1q21, a region amplified in high risk myeloma. We further uncovered that FAM72D functions as part of the FOXM1 transcription factor network controlling cell proliferation and survival and we evidenced an increased sensitivity of cells expressing high levels of FOXM1 and FAM72 to epigenetic drugs targeting histone deacetylases and DNA methyltransferases.

Polypyrimidine tract-binding proteins are essential for B cell development

Polypyrimidine tract-binding protein 1 (PTBP1) is a RNA-binding protein (RBP) expressed throughout B cell development. Deletion of Ptbp1 in mouse pro-B cells results in upregulation of PTBP2 and normal B cell development. We show that PTBP2 compensates for PTBP1 in B cell ontogeny as deletion of both Ptbp1 and Ptbp2 results in a complete block at the pro-B cell stage and a lack of mature B cells. In pro-B cells PTBP1 ensures precise synchronisation of the activity of cyclin dependent kinases at distinct stages of the cell cycle, suppresses S-phase entry and promotes progression into mitosis. PTBP1 controls mRNA abundance and alternative splicing of important cell cycle regulators including CYCLIN-D2, c-MYC, p107 and CDC25B. Our results reveal a previously unrecognised mechanism mediated by a RBP that is essential for B cell ontogeny and integrates transcriptional and post-translational determinants of progression through the cell cycle.

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.

Molecular Signature of Subtypes of Non-Small-Cell Lung Cancer by Large-Scale Transcriptional Profiling: Identification of Key Modules and Genes by Weighted Gene Co-Expression Network Analysis (WGCNA)

Non-small-cell lung cancer (NSCLC) represents a heterogeneous group of malignancies consisting essentially of adenocarcinoma (ADC) and squamous cell carcinoma (SCC). Although the diagnosis and treatment of ADC and SCC have been greatly improved in recent decades, there is still an urgent need to identify accurate transcriptome profile associated with the histological subtypes of NSCLC. The present study aims to identify the key dysregulated pathways and genes involved in the development of lung ADC and SCC and to relate them with the clinical traits. The transcriptional changes between tumour and normal lung tissues were investigated by RNA-seq. Gene ontology (GO), canonical pathways analysis with the prediction of upstream regulators, and weighted gene co-expression network analysis (WGCNA) to identify co-expressed modules and hub genes were used to explore the biological functions of the identified dysregulated genes. It was indicated that specific gene signatures differed significantly between ADC and SCC related to the distinct pathways. Of identified modules, four and two modules were the most related to clinical features in ADC and SCC, respectively. CTLA4, MZB1, NIP7, and BUB1B in ADC, as well as GNG11 and CCNB2 in SCC, are novel top hub genes in modules associated with tumour size, SUVmax, and recurrence-free survival. Our research provides a more effective understanding of the importance of biological pathways and the relationships between major genes in NSCLC in the perspective of searching for new molecular targets.

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.