Stage-specific Gene Research Articles

Plasmodium falciparum Development from Gametocyte to Oocyst: Insight from Functional Studies.

Malaria elimination may never succeed without the implementation of transmission-blocking strategies. The transmission of Plasmodium spp. parasites from the human host to the mosquito vector depends on circulating gametocytes in the peripheral blood of the vertebrate host. Once ingested by the mosquito during blood meals, these sexual forms undergo a series of radical morphological and metabolic changes to survive and progress from the gut to the salivary glands, where they will be waiting to be injected into the vertebrate host. The design of effective transmission-blocking strategies requires a thorough understanding of all the mechanisms that drive the development of gametocytes, gametes, sexual reproduction, and subsequent differentiation within the mosquito. The drastic changes in Plasmodium falciparum shape and function throughout its life cycle rely on the tight regulation of stage-specific gene expression. This review outlines the mechanisms involved in Plasmodium falciparum sexual stage development in both the human and mosquito vector, and zygote to oocyst differentiation. Functional studies unravel mechanisms employed by P. falciparum to orchestrate the expression of stage-specific functional products required to succeed in its complex life cycle, thus providing us with potential targets for developing new therapeutics. These mechanisms are based on studies conducted with various Plasmodium species, including predominantly P. falciparum and the rodent malaria parasites P. berghei. However, the great potential of epigenetics, genomics, transcriptomics, proteomics, and functional genetic studies to improve the understanding of malaria as a disease remains partly untapped because of limitations in studies using human malaria parasites and field isolates.

Transcriptomic profiling of Parkinson's disease brains reveals disease stage specific gene expression changes

Parkinson´s disease (PD) is a progressive neurodegenerative disorder characterized by both motor and non-motor symptoms. Aggravation of symptoms is mirrored by accumulation of protein aggregates mainly composed by alpha-synuclein in different brain regions, called Lewy bodies (LB). Previous studies have identified several molecular mechanisms as autophagy and inflammation playing a role in PD pathogenesis. Increased insights into mechanisms involved in early disease stages and driving the progression of the LB pathology are required for the development of disease-modifying strategies. Here, we aimed to elucidate disease stage-specific transcriptomic changes in brain tissue of well-characterized PD and control donors. We collected frontal cortex samples from 84 donors and sequenced both the coding and non-coding RNAs. We categorized our samples into groups based on their degree of LB pathology aiming to recapitulate a central aspect of disease progression. Using an analytical pipeline that corrected for sex, age at death, RNA quality, cell composition and unknown sources of variation, we found major disease stage-specific transcriptomic changes. Gene expression changes were most pronounced in donors at the disease stage when microscopic LB changes first occur in the sampled brain region. Additionally, we identified disease stage-specific enrichment of brain specific pathways and immune mechanisms. On the contrary, we showed that mitochondrial mechanisms are enriched throughout the disease course. Our data-driven approach also suggests a role for several poorly characterized lncRNAs in disease development and progression of PD. Finally, by combining genetic and epigenetic information, we highlighted two genes (MAP4K4 and PHYHIP) as candidate genes for future functional studies. Together our results indicate that transcriptomic dysregulation and associated functional changes are highly disease stage-specific, which has major implications for the study of neurodegenerative disorders.

Rhoptry neck protein 4 plays important roles during Plasmodium sporozoite infection of the mammalian liver.

During invasion, Plasmodium parasites secrete proteins from rhoptry and microneme apical end organelles, which have crucial roles in attaching to and invading target cells. A sporozoite stage-specific gene silencing system revealed that rhoptry neck protein 2 (RON2), RON4, and RON5 are important for sporozoite invasion of mosquito salivary glands. Here, we further investigated the roles of RON4 during sporozoite infection of the liver in vivo. Following intravenous inoculation of RON4-knockdown sporozoites into mice, we demonstrated that sporozoite RON4 has multiple functions during sporozoite traversal of sinusoidal cells and infection of hepatocytes. In vitro infection experiments using a hepatoma cell line revealed that secreted RON4 is involved in sporozoite adhesion to hepatocytes and has an important role in the early steps of hepatocyte infection. In addition, in vitro motility assays indicated that RON4 is required for sporozoite attachment to the substrate and the onset of migration. These findings indicate that RON4 is crucial for sporozoite migration toward and invasion of hepatocytes via attachment ability and motility.IMPORTANCEMalarial parasite transmission to mammals is established when sporozoites are inoculated by mosquitoes and migrate through the bloodstream to infect hepatocytes. Many aspects of the molecular mechanisms underpinning migration and cellular invasion remain largely unelucidated. By applying a sporozoite stage-specific gene silencing system in the rodent malarial parasite, Plasmodium berghei, we demonstrated that rhoptry neck protein 4 (RON4) is crucial for sporozoite infection of the liver in vivo. Combined with in vitro investigations, it was revealed that RON4 functions during a crossing of the sinusoidal cell layer and invading hepatocytes, at an early stage of liver infection, by mediating the sporozoite capacity for adhesion and the onset of motility. Since RON4 is also expressed in Plasmodium merozoites and Toxoplasma tachyzoites, our findings contribute to understanding the conserved invasion mechanisms of Apicomplexa parasites.

Functional characterization of Cinnamate 4-hydroxylase gene family in soybean (Glycine max).

Cinnamate 4-hydroxylase (C4H) is the first key cytochrome P450 monooxygenase (P450) enzyme in the phenylpropanoid pathway. It belongs to the CYP73 family of P450 superfamily, and catalyzes the conversion of trans-cinnamic acid to p-coumaric acid. Since p-coumaric acid serves as the precursor for the synthesis of a wide variety of metabolites involved in plant development and stress resistance, alteration in the expression of soybean C4H genes is expected to affect the downstream metabolite levels, and its ability to respond to stress. In this study, we identified four C4H genes in the soybean genome that are distributed into both class I and class II CYP73 family. GmC4H2, GmC4H14 and GmC4H20 displayed tissue- and developmental stage-specific gene expression patterns with their transcript accumulation at the highest level in root tissues. GmC4H10 appears to be a pseudogene as its transcript was not detected in any soybean tissues. Furthermore, protein homology modelling revealed substrate docking only for GmC4H2, GmC4H14 and GmC4H20. To demonstrate the function of GmC4Hs, we modified a cloning vector for the heterologous expression of P450s in yeast, and used it for microsomal protein production and enzyme assay. Our results confirmed that GmC4H2, GmC4H14 and GmC4H20 contain the ability to hydroxylate trans-cinnamic acid with varying efficiencies.

Dynamic transcriptome profiles and novel markers in bovine spermatogenesis revealed by single-cell sequencing

Testicular development is an important biological process in male and requires interaction between the male germ cells and somatic cells. However, the mechanisms of testicular development in livestock, particularly in cattle, are poorly understood. Furthermore, cellular heterogeneity hinders the profiling of different cell types at different developmental stages. In this study, we first performed a single-cell transcriptomic study of the bovine testis development during puberty by using 10× genomics single-cell RNA sequencing (scRNA-seq). By collecting the scRNA-seq data from 11,083 cells from prepubertal and pubertal bovine testes, a high-resolution scRNA-seq atlas was described, identifying 9 somatic and 13 spermatogenic clusters. We also distinguished several stage-specific marker genes for bovine germ cells and somatic cells, such as GRAF2 and MORC1 for SSC (spermatogonial stem cells), HJURP and TCF19 for differentiating spermatogonia, ARSE for immature Sertoli, CLEC12B for mature Sertoli, LOC112441470 for Leydig. In conclusion, we have examined the transcription levels and constructed the single-cell developmental maps of germ cells and somatic cells during testicular development in Angus cattle. The datasets provided new insights into spermatogenesis and testicular somatic cell development in cattle.

The Clinical Utility of lncRNAs and Their Application as Molecular Biomarkers in Breast Cancer.

Given their tumor-specific and stage-specific gene expression, long non-coding RNAs (lncRNAs) have demonstrated to be potential molecular biomarkers for diagnosis, prognosis, and treatment response. Particularly, the lncRNAs DSCAM-AS1 and GATA3-AS1 serve as examples of this because of their high subtype-specific expression profile in luminal B-like breast cancer. This makes them candidates to use as molecular biomarkers in clinical practice. However, lncRNA studies in breast cancer are limited in sample size and are restricted to the determination of their biological function, which represents an obstacle for its inclusion as molecular biomarkers of clinical utility. Nevertheless, due to their expression specificity among diseases, such as cancer, and their stability in body fluids, lncRNAs are promising molecular biomarkers that could improve the reliability, sensitivity, and specificity of molecular techniques used in clinical diagnosis. The development of lncRNA-based diagnostics and lncRNA-based therapeutics will be useful in routine medical practice to improve patient clinical management and quality of life.

Genome-wide chromatin accessibility and gene expression profiling during flatfish metamorphosis

Metamorphosis is a widely studied post-embryonic process in which many tissues undergo dramatic modifications to adapt to the new adult lifestyle. Flatfishes represent a good example of metamorphosis in teleost fishes. During metamorphosis of flatfish, organ regression and neoformation occur, with one of the most notable changes being the migration of one of the eyes to the other side of the body. In order to create a useful and reliable tool to advance the molecular study of metamorphosis in flatfish, we generated a chromatin accessible atlas as well as gene expression profile during four developmental stages ranging from a phylotypic to a post-metamorphic stage. We identified 29,019 differentially accessible chromatin regions and 3,253 differentially expressed genes. We found stage-specific regulatory regions and gene expression profiles, supporting the quality of the results. Our work provides strongly reproducible data for further studies to elucidate the regulatory elements that ensure successful metamorphosis in flatfish species.

Plasmodium falciparum gametocytes display global chromatin remodelling during sexual differentiation

BackgroundThe protozoan malaria parasite Plasmodium falciparum has a complex life cycle during which it needs to differentiate into multiple morphologically distinct life forms. A key process for transmission of the disease is the development of male and female gametocytes in the human blood, yet the mechanisms determining sexual dimorphism in these haploid, genetically identical sexual precursor cells remain largely unknown. To understand the epigenetic program underlying the differentiation of male and female gametocytes, we separated the two sexual forms by flow cytometry and performed RNAseq as well as comprehensive ChIPseq profiling of several histone variants and modifications.ResultsWe show that in female gametocytes the chromatin landscape is globally remodelled with respect to genome-wide patterns and combinatorial usage of histone variants and histone modifications. We identified sex specific differences in heterochromatin distribution, implicating exported proteins and ncRNAs in sex determination. Specifically in female gametocytes, the histone variants H2A.Z/H2B.Z were highly enriched in H3K9me3-associated heterochromatin. H3K27ac occupancy correlated with stage-specific gene expression, but in contrast to asexual parasites this was unlinked to H3K4me3 co-occupancy at promoters in female gametocytes.ConclusionsCollectively, we defined novel combinatorial chromatin states differentially organising the genome in gametocytes and asexual parasites and unravelled fundamental, sex-specific differences in the epigenetic code. Our chromatin maps represent an important resource for future understanding of the mechanisms driving sexual differentiation in P. falciparum.

A comparative gene expression matrix in Apoe-deficient mice identifies unique and atherosclerotic disease stage-specific gene regulation patterns in monocytes and macrophages

Background and aimsAtherosclerosis is a systemic and chronic inflammatory disease propagated by monocytes and macrophages. Yet, our knowledge on how transcriptome of these cells evolves in time and space is limited. We aimed at characterizing gene expression changes in site-specific macrophages and in circulating monocytes during the course of atherosclerosis. MethodsWe utilized apolipoprotein E-deficient mice undergoing one- and six-month high cholesterol diet to model early and advanced atherosclerosis. Aortic macrophages, peritoneal macrophages, and circulating monocytes from each mouse were subjected to bulk RNA-sequencing (RNA-seq). We constructed a comparative directory that profiles lesion- and disease stage-specific transcriptomic regulation of the three cell types in atherosclerosis. Lastly, the regulation of one gene, Gpnmb, whose expression positively correlated with atheroma growth, was validated using single-cell RNA-seq (scRNA-seq) of atheroma plaque from murine and human. ResultsThe convergence of gene regulation between the three investigated cell types was surprisingly low. Overall 3245 differentially expressed genes were involved in the biological modulation of aortic macrophages, among which less than 1% were commonly regulated by the remote monocytes/macrophages. Aortic macrophages regulated gene expression most actively during atheroma initiation. Through complementary interrogation of murine and human scRNA-seq datasets, we showcased the practicality of our directory, using the selected gene, Gpnmb, whose expression in aortic macrophages, and a subset of foamy macrophages in particular, strongly correlated with disease advancement during atherosclerosis initiation and progression. ConclusionsOur study provides a unique toolset to explore gene regulation of macrophage-related biological processes in and outside the atheromatous plaque at early and advanced disease stages.

Dynamics of histone acetylation during human early embryogenesis

It remains poorly understood about the regulation of gene and transposon transcription during human early embryogenesis. Here, we report that broad H3K27ac domains are genome-widely distributed in human 2-cell and 4-cell embryos and transit into typical peaks in the 8-cell embryos. The broad H3K27ac domains in early embryos before zygotic genome activation (ZGA) are also observed in mouse. It suggests that broad H3K27ac domains play conserved functions before ZGA in mammals. Intriguingly, a large portion of broad H3K27ac domains overlap with broad H3K4me3 domains. Further investigation reveals that histone deacetylases are required for the removal or transition of broad H3K27ac domains and ZGA. After ZGA, the number of typical H3K27ac peaks is dynamic, which is associated with the stage-specific gene expression. Furthermore, P300 is important for the establishment of H3K27ac peaks and the expression of associated genes in early embryos after ZGA. Our data also indicate that H3K27ac marks active transposons in early embryos. Interestingly, H3K27ac and H3K18ac signals rather than H3K9ac signals are enriched at ERVK elements in mouse embryos after ZGA. It suggests that different types of histone acetylations exert distinct roles in the activation of transposons. In summary, H3K27ac modification undergoes extensive reprogramming during early embryo development in mammals, which is associated with the expression of genes and transposons.

Re-analysis of publicly available methylomes using signal detection yields new information

Cytosine methylation is an epigenetic mark that participates in regulation of gene expression and chromatin stability in plants. Advancements in whole genome sequencing technologies have enabled investigation of methylome dynamics under different conditions. However, the computational methods for analyzing bisulfite sequence data have not been unified. Contention remains in the correlation of differentially methylated positions with the investigated treatment and exclusion of noise, inherent to these stochastic datasets. The prevalent approaches apply Fisher’s exact test, logistic, or beta regression, followed by an arbitrary cut-off for differences in methylation levels. A different strategy, the MethylIT pipeline, utilizes signal detection to determine cut-off based on a fitted generalized gamma probability distribution of methylation divergence. Re-analysis of publicly available BS-seq data from two epigenetic studies in Arabidopsis and applying MethylIT revealed additional, previously unreported results. Methylome repatterning in response to phosphate starvation was confirmed to be tissue-specific and included phosphate assimilation genes in addition to sulfate metabolism genes not implicated in the original study. During seed germination plants undergo major methylome reprogramming and use of MethylIT allowed us to identify stage-specific gene networks. We surmise from these comparative studies that robust methylome experiments must account for data stochasticity to achieve meaningful functional analyses.

Multiscale 3D genome reorganization during skeletal muscle stem cell lineage progression and aging.

Little is known about three-dimensional (3D) genome organization in skeletal muscle stem cells [also called satellite cells (SCs)]. Here, we comprehensively map the 3D genome topology reorganization during mouse SC lineage progression. Specifically, rewiring at the compartment level is most pronounced when SCs become activated. Marked loss in topologically associating domain (TAD) border insulation and chromatin looping also occurs during early activation process. Meanwhile, TADs can form TAD clusters and super-enhancer-containing TAD clusters orchestrate stage-specific gene expression. Furthermore, we uncover that transcription factor PAX7 is pivotal in enhancer-promoter (E-P) loop formation. We also identify cis-regulatory elements that are crucial for local chromatin organization at Pax7 locus and Pax7 expression. Lastly, we unveil that geriatric SC displays a prominent gain in long-range contacts and loss of TAD border insulation. Together, our results uncover that 3D chromatin extensively reorganizes at multiple architectural levels and underpins the transcriptome remodeling during SC lineage development and SC aging.

Transcriptome and morphological analyses of double flower formation in Dianthus chinensis

The double flower developmental process is regulated via a complex transcriptional regulatory network. To understand this highly dynamic and complex developmental process of Dianthus spp., we performed a comparative analysis of floral morphology and transcriptome dynamics in simple flowers and double flowers. We found that the primordium of double flowers of ‘X’ was larger in size compared to that of simple flowers of ‘L’ in Dianthus chinensis. RNA-seq and Weighted Gene Co-expression Network Analysis (WGCNA) during flower development, identified stage-specific gene network modules. Expression analysis by RNA-seq indicated that a group of genes related to floral meristem identity, primordia position and polarity were highly expressed in double flowers genotypes compared to simple flowers genotypes, suggesting their roles in double-petal formation. A total of 21 DEGs related to petal number were identified between simple and double flowers. The experiments of in situ hybridization revealed that DcaAP2L, DcaLFY and DcaUFO genes were expressed in the intra-sepal boundary and petal boundary. We proposed a potential transcriptional regulatory network for simple and double flower development. This study provides novel insights into the molecular mechanism underlying double flower formation in Dianthus spp.

A Stage-Specific Gene Expression Framework for Promoting Genetic Algorithm

Genetic algorithm (GA) has become very popular for tackling computationally expensive numerical optimization problems. To optimize their performance and convergence rate, we propose and demonstrate an embedded stage-specific gene expression framework (SGEF) with a two-tier coding scheme and a filter operator based on variational auto-encoder. For the first time, the structure of variational auto-encoder is used to map the parent population of a genetic algorithm from the uniformly distributed high dimensions to the Gaussian distributed low dimensions. The design of the filter operator which can maintain the trade-off balance between exploitation and exploration is also analyzed in detail. The experimental results demonstrated the good performance of the proposed SGEF in promoting both single-objective optimization genetic algorithms and multiple-objective optimization genetic algorithms.

Cross-Talk between Transcriptome Analysis and Dynamic Changes of Carbohydrates Identifies Stage-Specific Genes during the Flower Bud Differentiation Process of Chinese Cherry (Prunus pseudocerasus L.)

Flower bud differentiation is crucial to reproductive success in plants. In the present study, RNA-Seq and nutrients quantification were used to identify the stage-specific genes for flower bud differentiation with buds which characterize the marked change during flower bud formation from a widely grown Chinese cherry (Prunus pseudocerasus L.) cultivar 'Manaohong'. A KEGG enrichment analysis revealed that the sugar metabolism pathways dynamically changed. The gradually decreasing trend in the contents of total sugar, soluble sugar and protein implies that the differentiation was an energy-consuming process. Changes in the contents of D-glucose and sorbitol were conformed with the gene expression trends of bglX and SORD, respectively, which at least partially reflects a key role of the two substances in the transition from physiological to morphological differentiation. Further, the WRKY and SBP families were also significantly differentially expressed during the vegetative-to-reproductive transition. In addition, floral meristem identity genes, e.g., AP1, AP3, PI, AGL6, SEP1, LFY, and UFO demonstrate involvement in the specification of the petal and stamen primordia, and FPF1 might promote the onset of morphological differentiation. Conclusively, the available evidence justifies the involvement of sugar metabolism in the flower bud differentiation of Chinese cherry, and the uncovered candidate genes are beneficial to further elucidate flower bud differentiation in cherries.

PfHMGB2 has a role in malaria parasite mosquito infection.

Differentiation of asexually replicating parasites into gametocytes is critical for successful completion of the sexual phase of the malaria parasite life cycle. Gametes generated from gametocytes fuse to form a zygote which differentiates into ookinetes and oocysts. The sporozoites are formed inside oocysts which migrate to the salivary glands for next cycle of human infection. These morphologically and functionally distinct stages require stage-specific gene expression via specific transcriptional regulators. The capacity of high mobility group box (HMGB) proteins to interact with DNA in a sequence independent manner enables them to regulate higher order chromosome organization and regulation of gene expression. Plasmodium falciparum HMGB2 (PfHMGB2) shows a typical L- shaped predicted structure which is similar to mammalian HMG box proteins and shows very high protein sequence similarity to PyHMGB2 and PbHMGB2. Functional characterization of PfHMGB2 by gene deletion (Pfhmgb2¯) showed that knockout parasites develop normally as asexual stages and undergo gametocytogenesis. Transmission experiments revealed that Pfhmgb2¯ can infect mosquitoes and develop as oocyst stages. However, transmission was reduced compared to wild type (WT) parasites and as a consequence, the salivary gland sporozoites were reduced in number. In summary, we demonstrate that PfHMGB2 has no role in asexual growth and a modest role in sexual phase development and parasite transmission to the mosquito.

The dynamic dysregulated network identifies stage-specific markers during lung adenocarcinoma malignant progression and metastasis

Brain metastasis occurs in approximately 30% of patients with lung adenocarcinoma (LUAD) and is closely associated with poor prognosis, recurrence, and death. However, dynamic gene regulation and molecular mechanism driving LUAD progression remain poorly understood. In this study, we performed a comprehensive single-cell transcriptome analysis using data from normal, early stage, advanced stage, and brain metastasis LUAD. Our single-cell-level analysis reveals the cellular composition heterogeneity at different stages during LUAD progression. We identified stage-specific risk genes that could contribute to LUAD progression and metastasis by reprogramming immune-related and metabolic-related functions. We constructed an early advanced metastatic dysregulated network and revealed the dynamic changes in gene regulations during LUAD progression. We identified 6 early advanced (HLA-DRB1, HLA-DQB1, SFTPB, SFTPC, PLA2G1B, and FOLR1), 8 advanced metastasis (RPS15, RPS11, RPL13A, RPS24, HLA-DRB5, LYPLA1, KCNJ15, and PSMA3), and 2 common risk genes in different stages (SFTPD and HLA-DRA) as prognostic markers in LUAD. Particularly, decreased expression of HLA-DRA, HLA-DRB1, HLA-DQB1, and HLA-DRB5 refer poor prognosis in LUAD by controlling antigen processing and presentation and Tcell activation. Increased expression of PSMA3 and LYPLA1 refer poor prognosis by reprogramming fatty acid metabolism and RNA catabolic process. Our findings will help further understanding the pathobiology of brain metastases in LUAD.

A Multi-Cohort Gene Expression Classifier (ALLCatchR) Identifies B-Precursor ALL Subtypes and Their Developmental Trajectories across Age Groups

The molecular landscape of B cell precursor acute lymphoblastic leukemia (BCP-ALL) describes >20 disease subtypes defined by initiating genomic alterations and corresponding gene expression signatures. Some subtypes represent phenocopies with shared gene expression signatures but distinct or absent genomic drivers. The majority of BCP-ALL subtypes have been acknowledged as diagnostic entities by the currently revised classifications of hematologic malignancies (WHO-HAEM5 / International Consensus Classification). Thus, unified diagnostic approaches are needed to systematically identify subtypes in individual BCP-ALL cases in clinical routine. Gene expression signatures represent the shared functional equivalent of heterogeneous genomic driver alterations. Therefore, we established a gene expression-based classifier to identify BCP-ALL subtypes and related clinical and biological phenotypes. We analyzed a total of n=3477 BCP-ALL gene expression profiles obtained from transcriptome sequencing of adult (n=1253) and pediatric patients (n=2224) including publicly available (n=2074; Gu Z et al. Nat. Genet. 2019; Schmidt B. et al. Blood Adv. 2022) and own data sets (n=1403) from adult and pediatric European study groups (GMALL, AIEOP-BFM, CELL). Raw read counts from 6 cohorts representing heterogeneous library preps (i.e., stranded/unstranded; mRNA/total RNA), sequencing platforms/depths, and gene count quantifications were included. All samples had been previously allocated to molecular BCP-ALL subtypes based on integration of independent genomic profiling, gene fusion calling and gene expression. For classification, data with defined subgroup allocation were randomly split into training (n=1831; 70%) and testing sets (n=785; 30%) and two hold-out studies (n=305). To define gene expression signatures of 21 BCP-ALL subtypes, we used training set samples and extracted specific gene sets employing different feature selection algorithms as well as ranking of genes in individual subtypes. Based on these definitions we established a BCP-ALL subtype classifier (ALLCatchR), which integrates machine learning models (support vector machines) and nearest-neighbors associations of single sample gene set enrichment analyses (Fig. A). ALLCatchR performed molecular subtype allocation for the testing set at high sensitivity (0.934±0.165) and specificity (0.997±0.005) on average across subtypes. A similar performance was achieved in the two hold-out cohorts despite completely independent data structures which have not been shown to the classifier (sensitivity: 0.946±0.097 / specificity: 0.996±0.007). Overall accuracies of 0.952 and 0.92 for testing and hold-out-data sets indicate that ALLCatchR reliably identifies molecular BCP-ALL subtypes in mixed adult / pediatric cohorts based on gene expression profiles alone (Fig. A, below). Among n=44/1090 (4%) samples of the combined hold-out/testing set with ambiguous subtype allocation (deemed 'unclassified'), ALLCatchR provides candidate subtype allocations (sensitivity: 0.671±0.393 / specificity: 0.978±0.036) for confirmation based on genomic driver alterations. In addition, we modified our classifier to predict baseline variables such as patient's sex (accuracy: 0.991), immunophenotype (accuracy: 0.935) or blast count (R=0.633). To define the corresponding B cell differentiation stage of origin, we compared BCP-ALL subtypes to differentiation stage-specific gene sets extracted from 7 FACS-sorted lymphopoietic subsets from healthy bone marrow donors (n=4). We observed defined BCP-ALL subtype-specific enrichment patterns which can be grouped along 4 major differentiation stages (Pro-B to Pre-B-II-Large, Fig. 1B) with a high degree of similarity in independent adult and pediatric cohorts. Thus, B cell developmental trajectories underlying BCP-ALL subtypes are dissected at high resolution by gene expression profiling and differentiation stages may favor the selection of specific leukemogenic drivers. ALLCatchR is a gene expression classifier for 21 established BCP-ALL molecular subtypes, which performs at independently validated cohorts with an accuracy of >90% and allows imputation of clinical baseline variables as well as underlying B cell developmental trajectories. It provides a reliable basis for systematic diagnostic approaches in pediatric and adult BCP-ALL patients. Figure 1View largeDownload PPTFigure 1View largeDownload PPT Close modal

Sex separation unveils the functional plasticity of the vomeronasal organ in rabbits

Chemosensory cues are vital for social and sexual behaviours and are primarily detected and processed by the vomeronasal system (VNS), whose plastic capacity has been investigated in mice. However, studying chemosensory plasticity outside of laboratory conditions may give a more realistic picture of how the VNS adapts to a changing environment. Rabbits are a well-described model of chemocommunication since the discovery of the rabbit mammary pheromone and their vomeronasal organ (VNO) transcriptome was recently characterised, a first step to further study plasticity-mediated transcriptional changes. In this study, we assessed the plastic capacity of the rabbit male and female VNO under sex-separation vs. sex-combined scenarios, including adults and juveniles, to determine whether the rabbit VNO is plastic and, if so, whether such plasticity is already established at early stages of life. First, we characterised the number of differentially expressed genes (DEGs) between the VNO of rabbit male and female under sex-separation and compared it to sex-combined individuals, both in adults and juveniles, finding that differences between male and female were larger in a sex-separated scenario. Secondly, we analysed the number of DEGs between sex-separated and sex-combined scenarios, both in males and females. In adults, both sexes showed a high number of DEGs while in juveniles only females showed differences. Additionally, the vomeronasal receptor genes were strikingly downregulated in sex-separated adult females, whereas in juveniles upregulation was shown for the same condition, suggesting a role of VRs in puberty onset. Finally, we described the environment-modulated plastic capacity of genes involved in reproduction, immunity and VNO functional activity, including G-protein coupled receptors. Our results show that sex-separation induces sex- and stage-specific gene expression differences in the VNO of male and female rabbit, both in adults and juveniles. These results bring out for the first time the plastic capacity of the rabbit VNO, supporting its functional adaptation to specifically respond to a continuous changing environment. Finally, species-specific differences and individual variability should always be considered in VNO studies and overall chemocommunication research.

Stage specific gene expression of folate mediated one-carbon metabolism enzymes and transporters in buffalo oocytes and pre-implantation embryos

DNA synthesis and methylations are crucial during pre-implantation embryonic development, and are mediated by one-carbon metabolism of folates. Folates, transported into the cells via folate receptors (FOLR1 and FOLR2) and carriers (SLC19A1), are metabolized by various enzymes involved in folate-methionine cycle. However, the variations in temporal expression of folate transporters and folate-methionine cycle enzymes during pre-implantation embryo development is obscure. Thus, the present study aimed to investigate the differential expression of the genes for folate transporters and folate-methionine cycle enzymes. We also examined the expression of folate transport proteins in different pre-implantation development stages. Immature buffalo oocytes were matured in maturation medium followed by in vitro fertilization and culture at standard culture conditions. The temporal pattern of gene expression in buffalo, when compared to previous studies, indicated an inter-specific variation. The transcripts of some enzymes and folate transporters were significantly upregulated after zygotic genome activation. The transcripts as well as proteins for FOLR1, FOLR2 and SLC19A1 were present in oocytes and all the pre-implantation embryo stages. FOLR1 was present in the nuclei of different stages of developing embryos but not in the metaphase (MII) oocytes. As a result, the present study advocates the existence of active folate transport in buffalo oocytes and pre-implantation embryos. The data provided by the analysis of differential gene expression of folate transporters and metabolic enzymes would likely contribute to a better understanding of the role of folates in embryo development as well as advancements in assisted reproductive technologies.