PIM kinases regulate early human Th17 cell differentiation

The serine/threonine-specific Moloney murine leukemia virus (PIM) kinase family (i.e., PIM1, PIM2, and PIM3) has been extensively studied in tumorigenesis. PIM kinases are downstream of several cytokine signaling pathways that drive immune-mediated diseases. Uncontrolled T helper 17 (Th17) cell activation has been associated with the pathogenesis of autoimmunity. However, the detailed molecular function of PIMs in human Th17 cell regulation has yet to be studied. In the present study, we comprehensively investigated how the three PIMs simultaneously alter transcriptional gene regulation during early human Th17 cell differentiation. By combining PIM triple knockdown with bulk and scRNA-seq approaches, we found that PIM deficiency promotes the early expression of key Th17-related genes while suppressing Th1-lineage genes. Further, PIMs modulate Th cell signaling, potentially via STAT1 and STAT3. Overall, our study highlights the inhibitory role of PIMs in human Th17 cell differentiation, thereby suggesting their association with autoimmune phenotypes.

Serum Insights: Leveraging the Power of miRNA Profiling as an Early Diagnostic Tool for Non-Small Cell Lung Cancer

Non-small cell lung cancer is the predominant form of lung cancer and is associated with a poor prognosis. MiRNAs implicated in cancer initiation and progression can be easily detected in liquid biopsy samples and have the potential to serve as non-invasive biomarkers. In this study, we employed next-generation sequencing to globally profile miRNAs in serum samples from 71 early-stage NSCLC patients and 47 non-cancerous pulmonary condition patients. Preliminary analysis of differentially expressed miRNAs revealed 28 upregulated miRNAs in NSCLC compared to the control group. Functional enrichment analyses unveiled their involvement in NSCLC signaling pathways. Subsequently, we developed a gradient-boosting decision tree classifier based on 2588 miRNAs, which demonstrated high accuracy (0.837), sensitivity (0.806), and specificity (0.859) in effectively distinguishing NSCLC from non-cancerous individuals. Shapley Additive exPlanations analysis improved the model metrics by identifying the top 15 miRNAs with the strongest discriminatory value, yielding an AUC of 0.96 ± 0.04, accuracy of 0.896, sensitivity of 0.884, and specificity of 0.903. Our study establishes the potential utility of a non-invasive serum miRNA signature as a supportive tool for early detection of NSCLC while also shedding light on dysregulated miRNAs in NSCLC biology. For enhanced credibility and understanding, further validation in an independent cohort of patients is warranted.

Cell-connectivity-guided trajectory inference from single-cell data

MOTIVATION

Single-cell RNA-sequencing enables cell-level investigation of cell differentiation, which can be modelled using trajectory inference methods. While tremendous effort has been put into designing these methods, inferring accurate trajectories automatically remains difficult. Therefore, the standard approach involves testing different trajectory inference methods and picking the trajectory giving the most biologically sensible model. As the default parameters are often suboptimal, their tuning requires methodological expertise.

RESULTS

We introduce Totem, an open-source, easy-to-use R package designed to facilitate inference of tree-shaped trajectories from single-cell data. Totem generates a large number of clustering results, estimates their topologies as minimum spanning trees, and uses them to measure the connectivity of the cells. Besides automatic selection of an appropriate trajectory, cell connectivity enables to visually pinpoint branching points and milestones relevant to the trajectory. Furthermore, testing different trajectories with Totem is fast, easy, and does not require in-depth methodological knowledge.

AVAILABILITY AND IMPLEMENTATION

Totem is available as an R package at https://github.com/elolab/Totem.

Spectrum and Density of Gamma and X-ray Induced Mutations in a Non-Model Rice Cultivar

Physical mutagens are a powerful tool used for genetic research and breeding for over eight decades. Yet, when compared to chemical mutagens, data sets on the effect of different mutagens and dosages on the spectrum and density of induced mutations remain lacking. To address this, we investigated the landscape of mutations induced by gamma and X-ray radiation in the most widely cultivated crop species: rice. A mutant population of a tropical upland rice, Oryza sativa L., was generated and propagated via self-fertilization for seven generations. Five dosages ranging from 75 Gy to 600 Gy in both X-ray and gamma-irradiated material were applied. In the process of a forward genetic screens, 11 unique rice mutant lines showing phenotypic variation were selected for mutation analysis via whole-genome sequencing. Thousands of candidate mutations were recovered in each mutant with single base substitutions being the most common, followed by small indels and structural variants. Higher dosages resulted in a higher accumulation of mutations in gamma-irradiated material, but not in X-ray-treated plants. The in vivo role of all annotated rice genes is yet to be directly investigated. The ability to induce a high density of single nucleotide and structural variants through mutagenesis will likely remain an important approach for functional genomics and breeding.

Benchmarking methods for detecting differential states between conditions from multi-subject single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) enables researchers to quantify transcriptomes of thousands of cells simultaneously and study transcriptomic changes between cells. scRNA-seq datasets increasingly include multisubject, multicondition experiments to investigate cell-type-specific differential states (DS) between conditions. This can be performed by first identifying the cell types in all the subjects and then by performing a DS analysis between the conditions within each cell type. Naïve single-cell DS analysis methods that treat cells statistically independent are subject to false positives in the presence of variation between biological replicates, an issue known as the pseudoreplicate bias. While several methods have already been introduced to carry out the statistical testing in multisubject scRNA-seq analysis, comparisons that include all these methods are currently lacking. Here, we performed a comprehensive comparison of 18 methods for the identification of DS changes between conditions from multisubject scRNA-seq data. Our results suggest that the pseudobulk methods performed generally best. Both pseudobulks and mixed models that model the subjects as a random effect were superior compared with the naïve single-cell methods that do not model the subjects in any way. While the naïve models achieved higher sensitivity than the pseudobulk methods and the mixed models, they were subject to a high number of false positives. In addition, accounting for subjects through latent variable modeling did not improve the performance of the naïve methods.

COVID-19-specific transcriptomic signature detectable in blood across multiple cohorts

The coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is spreading across the world despite vast global vaccination efforts. Consequently, many studies have looked for potential human host factors and immune mechanisms associated with the disease. However, most studies have focused on comparing COVID-19 patients to healthy controls, while fewer have elucidated the specific host factors distinguishing COVID-19 from other infections. To discover genes specifically related to COVID-19, we reanalyzed transcriptome data from nine independent cohort studies, covering multiple infections, including COVID-19, influenza, seasonal coronaviruses, and bacterial pneumonia. The identified COVID-19-specific signature consisted of 149 genes, involving many signals previously associated with the disease, such as induction of a strong immunoglobulin response and hemostasis, as well as dysregulation of cell cycle-related processes. Additionally, potential new gene candidates related to COVID-19 were discovered. To facilitate exploration of the signature with respect to disease severity, disease progression, and different cell types, we also offer an online tool for easy visualization of the selected genes across multiple datasets at both bulk and single-cell levels.

Long Intergenic Noncoding RNA MIAT as a Regulator of Human Th17 Cell Differentiation

T helper 17 (Th17) cells protect against fungal and bacterial infections and are implicated in autoimmunity. Several long intergenic noncoding RNAs (lincRNA) are induced during Th17 differentiation, however, their contribution to Th17 differentiation is poorly understood. We aimed to characterize the function of the lincRNA Myocardial Infarction Associated Transcript (MIAT) during early human Th17 cell differentiation. We found MIAT to be upregulated early after induction of human Th17 cell differentiation along with an increase in the chromatin accessibility at the gene locus. STAT3, a key regulator of Th17 differentiation, directly bound to the MIAT promoter and induced its expression during the early stages of Th17 cell differentiation. MIAT resides in the nucleus and regulates the expression of several key Th17 genes, including IL17A, IL17F, CCR6 and CXCL13, possibly by altering the chromatin accessibility of key loci, including IL17A locus. Further, MIAT regulates the expression of protein kinase C alpha (PKCα), an upstream regulator of IL17A. A reanalysis of published single-cell RNA-seq data showed that MIAT was expressed in T cells from the synovium of RA patients. Our results demonstrate that MIAT contributes to human Th17 differentiation by upregulating several genes implicated in Th17 differentiation. High MIAT expression in T cells of RA patient synovia suggests a possible role of MIAT in Th17 mediated autoimmune pathologies.

A systematic comparison of FOSL1, FOSL2 and BATF-mediated transcriptional regulation during early human Th17 differentiation

Th17 cells are essential for protection against extracellular pathogens, but their aberrant activity can cause autoimmunity. Molecular mechanisms that dictate Th17 cell-differentiation have been extensively studied using mouse models. However, species-specific differences underscore the need to validate these findings in human. Here, we characterized the human-specific roles of three AP-1 transcription factors, FOSL1, FOSL2 and BATF, during early stages of Th17 differentiation. Our results demonstrate that FOSL1 and FOSL2 co-repress Th17 fate-specification, whereas BATF promotes the Th17 lineage. Strikingly, FOSL1 was found to play different roles in human and mouse. Genome-wide binding analysis indicated that FOSL1, FOSL2 and BATF share occupancy over regulatory regions of genes involved in Th17 lineage commitment. These AP-1 factors also share their protein interacting partners, which suggests mechanisms for their functional interplay. Our study further reveals that the genomic binding sites of FOSL1, FOSL2 and BATF harbour hundreds of autoimmune disease-linked SNPs. We show that many of these SNPs alter the ability of these transcription factors to bind DNA. Our findings thus provide critical insights into AP-1-mediated regulation of human Th17-fate and associated pathologies.

Early DNA methylation changes in children developing beta cell autoimmunity at a young age

AIMS/HYPOTHESIS

Type 1 diabetes is a chronic autoimmune disease of complex aetiology, including a potential role for epigenetic regulation. Previous epigenomic studies focused mainly on clinically diagnosed individuals. The aim of the study was to assess early DNA methylation changes associated with type 1 diabetes already before the diagnosis or even before the appearance of autoantibodies.

METHODS

Reduced representation bisulphite sequencing (RRBS) was applied to study DNA methylation in purified CD4⁺ T cell, CD8⁺ T cell and CD4^-CD8^- cell fractions of 226 peripheral blood mononuclear cell samples longitudinally collected from seven type 1 diabetes-specific autoantibody-positive individuals and control individuals matched for age, sex, HLA risk and place of birth. We also explored correlations between DNA methylation and gene expression using RNA sequencing data from the same samples. Technical validation of RRBS results was performed using pyrosequencing.

RESULTS

We identified 79, 56 and 45 differentially methylated regions in CD4⁺ T cells, CD8⁺ T cells and CD4^-CD8^- cell fractions, respectively, between type 1 diabetes-specific autoantibody-positive individuals and control participants. The analysis of pre-seroconversion samples identified DNA methylation signatures at the very early stage of disease, including differential methylation at the promoter of IRF5 in CD4⁺ T cells. Further, we validated RRBS results using pyrosequencing at the following CpG sites: chr19:18118304 in the promoter of ARRDC2; chr21:47307815 in the intron of PCBP3; and chr14:81128398 in the intergenic region near TRAF3 in CD4⁺ T cells.

CONCLUSIONS/INTERPRETATION

These preliminary results provide novel insights into cell type-specific differential epigenetic regulation of genes, which may contribute to type 1 diabetes pathogenesis at the very early stage of disease development. Should these findings be validated, they may serve as a potential signature useful for disease prediction and management.

scShaper: an ensemble method for fast and accurate linear trajectory inference from single-cell RNA-seq data

MOTIVATION

Computational models are needed to infer a representation of the cells, i.e. a trajectory, from single-cell RNA-sequencing data that model cell differentiation during a dynamic process. Although many trajectory inference methods exist, their performance varies greatly depending on the dataset and hence there is a need to establish more accurate, better generalizable methods.

RESULTS

We introduce scShaper, a new trajectory inference method that enables accurate linear trajectory inference. The ensemble approach of scShaper generates a continuous smooth pseudotime based on a set of discrete pseudotimes. We demonstrate that scShaper is able to infer accurate trajectories for a variety of trigonometric trajectories, including many for which the commonly used principal curves method fails. A comprehensive benchmarking with state-of-the-art methods revealed that scShaper achieved superior accuracy of the cell ordering and, in particular, the differentially expressed genes. Moreover, scShaper is a fast method with few hyperparameters, making it a promising alternative to the principal curves method for linear pseudotemporal ordering.

AVAILABILITY AND IMPLEMENTATION

scShaper is available as an R package at https://github.com/elolab/scshaper. The test data are available at https://doi.org/10.5281/zenodo.5734488.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Differential ATAC-seq and ChIP-seq peak detection using ROTS

Changes in cellular chromatin states fine-tune transcriptional output and ultimately lead to phenotypic changes. Here we propose a novel application of our reproducibility-optimized test statistics (ROTS) to detect differential chromatin states (ATAC-seq) or differential chromatin modification states (ChIP-seq) between conditions. We compare the performance of ROTS to existing and widely used methods for ATAC-seq and ChIP-seq data using both synthetic and real datasets. Our results show that ROTS outperformed other commonly used methods when analyzing ATAC-seq data. ROTS also displayed the most accurate detection of small differences when modeling with synthetic data. We observed that two-step methods that require the use of a separate peak caller often more accurately called enrichment borders, whereas one-step methods without a separate peak calling step were more versatile in calling sub-peaks. The top ranked differential regions detected by the methods had marked correlation with transcriptional differences of the closest genes. Overall, our study provides evidence that ROTS is a useful addition to the available differential peak detection methods to study chromatin and performs especially well when applied to study differential chromatin states in ATAC-seq data.

ILoReg: a tool for high-resolution cell population identification from single-cell RNA-seq data

MOTIVATION

Single-cell RNA-seq allows researchers to identify cell populations based on unsupervised clustering of the transcriptome. However, subpopulations can have only subtle transcriptomic differences and the high dimensionality of the data makes their identification challenging.

RESULTS

We introduce ILoReg, an R package implementing a new cell population identification method that improves identification of cell populations with subtle differences through a probabilistic feature extraction step that is applied before clustering and visualization. The feature extraction is performed using a novel machine learning algorithm, called iterative clustering projection (ICP), that uses logistic regression and clustering similarity comparison to iteratively cluster data. Remarkably, ICP also manages to integrate feature selection with the clustering through L1-regularization, enabling the identification of genes that are differentially expressed between cell populations. By combining solutions of multiple ICP runs into a single consensus solution, ILoReg creates a representation that enables investigating cell populations with a high resolution. In particular, we show that the visualization of ILoReg allows segregation of immune and pancreatic cell populations in a more pronounced manner compared with current state-of-the-art methods.

AVAILABILITY AND IMPLEMENTATION

ILoReg is available as an R package at https://bioconductor.org/packages/ILoReg.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Reproducibility-optimized detection of differential DNA methylation

Aim: DNA methylation is a key epigenetic mechanism regulating gene expression. Identifying differentially methylated regions is integral to DNA methylation analysis and there is a need for robust tools reliably detecting regions with significant differences in their methylation status. Materials & methods: We present here a reproducibility-optimized test statistic (ROTS) for detection of differential DNA methylation from high-throughput sequencing or array-based data. Results: Using both simulated and real data, we demonstrate the ability of ROTS to identify differential methylation between sample groups. Conclusion: Compared with state-of-the-art methods, ROTS shows competitive sensitivity and specificity in detecting consistently differentially methylated regions.

Leptin and fractalkine: novel subcutaneous cytokines in burn injury

Burn injury is a pathology underpinned by progressive and aberrant inflammation. It is a major clinical challenge to survival and quality of life. Although the complex local and disseminating pathological processes of a burn injury ultimately stem from local tissue damage, to date relatively few studies have attempted to characterise the local inflammatory mediator profile. Here, cytokine content and associated transcriptional changes were measured in rat skin for three hours immediately following induction of a scald-type (60°C, 2 min) burn injury model. Leptin (P=0.0002) and fractalkine (P=0.0478) concentrations were significantly elevated post-burn above pre-burn and control site values, coinciding with the development of burn site oedema and differential expression of leptin mRNA (P=0.0004). Further, gene sequencing enrichment analysis indicated cytokine-cytokine receptor interaction (P=1.45×10-6). Subsequent behavioural studies demonstrated that, following subcutaneous injection into the dorsum of the paw, both leptin and fractalkine induced mechanical allodynia, heat hyperalgesia and the recruitment of macrophages. This is the first report of leptin elevation specifically at the burn site, and the first report of fractalkine elevation in any tissue post-burn which, together with the functional findings, calls for exploration of the influence of these cytokines on pain, inflammation and burn wound progression. In addition, targeting these signalling molecules represents a therapeutic potential as early formative mediators of these pathological processes.

Evaluation of Transcriptomic Regulations behind Metabolic Syndrome in Obese and Lean Subjects

Multiple mechanisms have been suggested to confer to the pathophysiology of metabolic syndrome (MetS), however despite great interest from the scientific community, the exact contribution of each of MetS risk factors still remains unclear. The present study aimed to investigate molecular signatures in peripheral blood of individuals affected by MetS and different degrees of obesity. Metabolic health of 1204 individuals from 1000PLUS cohort was assessed, and 32 subjects were recruited to four study groups: MetS lean, MetS obese, “healthy obese”, and healthy lean. Whole-blood transcriptome next generation sequencing with functional data analysis were carried out. MetS obese and MetS lean study participants showed the upregulation of genes involved in inflammation and coagulation processes: granulocyte adhesion and diapedesis (p < 0.0001, p = 0.0063), prothrombin activation pathway (p = 0.0032, p = 0.0091), coagulation system (p = 0.0010, p = 0.0155). The results for “healthy obese” indicate enrichment in molecules associated with protein synthesis (p < 0.0001), mitochondrial dysfunction (p < 0.0001), and oxidative phosphorylation (p < 0.0001). Our results suggest that MetS is related to the state of inflammation and vascular system changes independent of excess body weight. Furthermore, “healthy obese”, despite not fulfilling the criteria for MetS diagnosis, seems to display an intermediate state with a lower degree of metabolic abnormalities, before they proceed to a full blown MetS.

Estrogen Signaling Drives Ciliogenesis in Human Endometrial Organoids

The human endometrium is the inner lining of the uterus consisting of stromal and epithelial (secretory and ciliated) cells. It undergoes a hormonally regulated monthly cycle of growth, differentiation, and desquamation. However, how these cyclic changes control the balance between secretory and ciliated cells remains unclear. Here, we established endometrial organoids to investigate the estrogen (E2)-driven control of cell fate decisions in human endometrial epithelium. We demonstrate that they preserve the structure, expression patterns, secretory properties, and E2 responsiveness of their tissue of origin. Next, we show that the induction of ciliated cells is orchestrated by the coordinated action of E2 and NOTCH signaling. Although E2 is the primary driver, inhibition of NOTCH signaling provides a permissive environment. However, inhibition of NOTCH alone is not sufficient to trigger ciliogenesis. Overall, we provide insights into endometrial biology and propose endometrial organoids as a robust and powerful model for studying ciliogenesis in vitro.

Microfluidic nutrient gradient-based three-dimensional chondrocyte culture-on-a-chip as an in vitro equine arthritis model

In this work, we describe a microfluidic three-dimensional (3D) chondrocyte culture mimicking in vivo articular chondrocyte morphology, cell distribution, metabolism, and gene expression. This has been accomplished by establishing a physiologic nutrient diffusion gradient across the simulated matrix, while geometric design constraints of the microchambers drive native-like cellular behavior. Primary equine chondrocytes remained viable for the extended culture time of 3 weeks and maintained the low metabolic activity and high Sox9, aggrecan, and Col2 expression typical of articular chondrocytes. Our microfluidic 3D chondrocyte microtissues were further exposed to inflammatory cytokines to establish an animal-free, in vitro osteoarthritis model. Results of our study indicate that our microtissue model emulates the basic characteristics of native cartilage and responds to biochemical injury, thus providing a new foundation for exploration of osteoarthritis pathophysiology in both human and veterinary patients.

Transcriptional Alterations in the Trigeminal Ganglia, Nucleus and Peripheral Blood Mononuclear Cells in a Rat Orofacial Pain Model

Orofacial pain and headache disorders are among the most debilitating pain conditions. While the pathophysiological basis of these disorders may be diverse, it is generally accepted that a common mechanism behind the arising pain is the sensitization of extra- and intracranial trigeminal primary afferents. In the present study we investigated gene expression changes in the trigeminal ganglia (TRG), trigeminal nucleus caudalis (TNC) and peripheral blood mononuclear cells (PBMC) evoked by Complete Freund's Adjuvant (CFA)-induced orofacial inflammation in rats, as a model of trigeminal sensitization. Microarray analysis revealed 512 differentially expressed genes between the ipsi- and contralateral TRG samples 7 days after CFA injection. Time-dependent expression changes of G-protein coupled receptor 39 (Gpr39), kisspeptin-1 receptor (Kiss1r), kisspeptin (Kiss1), as well as synaptic plasticity-associated Lkaaear1 (Lkr) and Neurod2 mRNA were described on the basis of qPCR results. The greatest alterations were observed on day 3 ipsilaterally, when orofacial mechanical allodynia reached its maximum. This corresponded well with patterns of neuronal (Fosb), microglia (Iba1), and astrocyte (Gfap) activation markers in both TRG and TNC, and interestingly also in PBMCs. This is the first description of up- and downregulated genes both in primary and secondary sensory neurones of the trigeminovascular system that might play important roles in neuroinflammatory activation mechanisms. We are the first to show transcriptomic alterations in the PBMCs that are similar to the neuronal changes. These results open new perspectives and initiate further investigations in the research of trigeminal pain disorders.

Linear ubiquitination by LUBEL has a role in Drosophila heat stress response

The HOIP ubiquitin E3 ligase generates linear ubiquitin chains by forming a complex with HOIL‐1L and SHARPIN in mammals. Here, we provide the first evidence of linear ubiquitination induced by a HOIP orthologue in Drosophila. We identify Drosophila CG11321, which we named Linear Ubiquitin E3 ligase (LUBEL), and find that it catalyzes linear ubiquitination in vitro. We detect endogenous linear ubiquitin chain‐derived peptides by mass spectrometry in Drosophila Schneider 2 cells and adult flies. Furthermore, using CRISPR/Cas9 technology, we establish linear ubiquitination‐defective flies by mutating residues essential for the catalytic activity of LUBEL. Linear ubiquitination signals accumulate upon heat shock in flies. Interestingly, flies with LUBEL mutations display reduced survival and climbing defects upon heat shock, which is also observed upon specific LUBEL depletion in muscle. Thus, LUBEL is involved in the heat response by controlling linear ubiquitination in flies.

IRF4 Deficiency Leads to Altered BCR Signalling Revealed by Enhanced PI3K Pathway, Decreased SHIP Expression and Defected Cytoskeletal Responses

The graded expression of transcription factor interferon regulatory factor 4 (IRF4) regulates B cell development and is critical for plasma cell differentiation. However, the mechanisms, by which IRF4 elicits its crucial tasks, are largely unknown. To characterize the molecular targets of IRF4 in B cells, we established an IRF4-deficient DT40 B cell line. We found that in the absence of IRF4, the expression of several molecules involved in BCR signalling was altered. For example, the expression of B cell adaptor for PI3K (BCAP) was upregulated, whereas the SHIP (SH2-containing Inositol 5?-Phosphatase) expression was downregulated. These molecular unbalances were accompanied by increased BCR-induced calcium signalling, attenuated B cell linker protein (BLNK) and ERK activity and enhanced activity of PI3K/protein kinase B (Akt) pathway. Further, the IRF4-deficient cells showed dramatically diminished cytoskeletal responses to anti-IgM cross-linking. Our results show that IRF4 has an important role in the regulation of BCR signalling and help to shed light on the molecular mechanisms of B cell development and germinal centre response.

Retene causes multifunctional transcriptomic changes in the heart of rainbow trout (Oncorhynchus mykiss) embryos

Fish are particularly sensitive to aryl hydrocarbon receptor (AhR)-mediated developmental toxicity. The molecular mechanisms behind these adverse effects have remained largely unresolved in salmonids, and for AhR-agonistic polycyclic aromatic hydrocarbons (PAHs). This study explored the cardiac transcriptome of rainbow trout (Oncorhynchus mykiss) eleuteroembryos exposed to retene, an AhR-agonistic PAH. The embryos were exposed to retene (nominal concentration 32 μg/L) and control, their hearts were collected before, at and after the onset of the visible signs of developmental toxicity, and transcriptomic changes were studied by microarray analysis. Retene up- or down-regulated 122 genes. The largest Gene Ontology groups were signal transduction, transcription, apoptosis, cell growth, cytoskeleton, cell adhesion/mobility, cardiovascular development, xenobiotic metabolism, protein metabolism, lipid metabolism and transport, and amino acid metabolism. Together these findings suggest that retene affects multiple signaling cascades in the heart of rainbow trout embryos, and potentially disturbs processes related to cardiovascular development and function.

Differential Promoter Methylation of Macrophage Genes Is Associated With Impaired Vascular Growth in Ischemic Muscles of Hyperlipidemic and Type 2 Diabetic Mice: Genome-Wide Promoter Methylation Study

RATIONALE

Hyperlipidemia and type 2 diabetes mellitus (T2DM) severely impair adaptive vascular growth responses in ischemic muscles. This is largely attributed to dysregulated gene expression, although details of the changes are unknown.

OBJECTIVE

To define the role of promoter methylation in adaptive vascular growth in hyperlipidemia (LDLR(-/-)ApoB(100/100)) and T2DM (IGF-II/LDLR(-/-)ApoB(100/100)) mouse models of hindlimb ischemia.

METHODS AND RESULTS

Unilateral hindlimb ischemia was induced by ligating femoral artery. Perfusion was assessed using ultrasound, and capillary and arteriole parameters were assessed using immunohistochemistry. Genome-wide methylated DNA sequencing was performed with DNA isolated from ischemic muscle, tissue macrophages (Mϕs), and endothelial cells. Compared with the controls, hyperlipidemia and T2DM mice showed impaired perfusion recovery, which was associated with impaired angiogenesis and arteriogenesis. Genome-wide proximal promoter DNA methylation analysis suggested differential patterns of methylation in Mϕ genes in ischemic muscles. Classically activated M1-Mϕ gene promoters, including Cfb, Serping1, and Tnfsf15, were significantly hypomethylated, whereas alternatively activated M2-Mϕ gene promoters, including Nrp1, Cxcr4, Plxnd1, Arg1, Cdk18, and Fes, were significantly hypermethylated in Mϕs isolated from hyperlipidemia and T2DM ischemic muscles compared with controls. These results combined with mRNA expression and immunohistochemistry showed the predominance of proinflammatory M1-Mϕs, compared with anti-inflammatory and proangiogenic M2-Mϕs in hyperlipidemia and T2DM ischemic muscles.

CONCLUSIONS

We found significant promoter hypomethylation of genes typical for proinflammatory M1-Mϕs and hypermethylation of anti-inflammatory, proangiogenic M2-Mϕ genes in hyperlipidemia and T2DM ischemic muscles. Epigenetic alterations modify Mϕ phenotype toward proinflammatory M1 as opposed to anti-inflammatory, proangiogenic, and tissue repair M2 phenotype, which may contribute to the impaired adaptive vascular growth under these pathological conditions.

Whole transcriptome characterization of the effects of dehydration and rehydration on Cladonia rangiferina, the grey reindeer lichen

BACKGROUND

Lichens are symbiotic organisms with a fungal and an algal or a cyanobacterial partner. Lichens inhabit some of the harshest climates on earth and most lichen species are desiccation-tolerant. Lichen desiccation-tolerance has been studied at the biochemical level and through proteomics, but the underlying molecular genetic mechanisms remain largely unexplored. The objective of our study was to examine the effects of dehydration and rehydration on the gene expression of Cladonia rangiferina.

RESULTS

Samples of C. rangiferina were collected at several time points during both the dehydration and rehydration process and the gene expression intensities were measured using a custom DNA microarray. Several genes, which were differentially expressed in one or more time points, were identified. The microarray results were validated using qRT-PCR analysis. Enrichment analysis of differentially expressed transcripts was also performed to identify the Gene Ontology terms most associated with the rehydration and dehydration process.

CONCLUSIONS

Our data identify differential expression patterns for hundreds of genes that are modulated during dehydration and rehydration in Cladonia rangiferina. These dehydration and rehydration events clearly differ from each other at the molecular level and the largest changes to gene expression are observed within minutes following rehydration. Distinct changes are observed during the earliest stage of rehydration and the mechanisms not appear to be shared with the later stages of wetting or with drying. Several of the most differentially expressed genes are similar to genes identified in previous studies that have investigated the molecular mechanisms of other desiccation-tolerant organisms. We present here the first microarray experiment for any lichen species and have for the first time studied the genetic mechanisms behind lichen desiccation-tolerance at the whole transcriptome level.

SAP30L (Sin3A-associated protein 30-like) is involved in regulation of cardiac development and hematopoiesis in zebrafish embryos

The Sin3A-associated proteins SAP30 and SAP30L share 70% sequence identity and are part of the multiprotein Sin3A corepressor complex. They participate in gene repression events by linking members of the complex and stabilizing interactions among the protein members as well as between proteins and DNA. While most organisms have both SAP30 and SAP30L, the zebrafish is exceptional because it only has SAP30L. Here we demonstrate that SAP30L is expressed ubiquitously in embryonic and adult zebrafish tissues. Knockdown of SAP30L using morpholino-mediated technology resulted in a morphant phenotype manifesting as cardiac insufficiency and defective hemoglobinization of red blood cells. A microarray analysis of gene expression in SAP30L morphant embryos revealed changes in the expression of genes involved in regulation of transcription, TGF-beta signaling, Wnt-family transcription factors, and nuclear genes encoding mitochondrial proteins. The expression of the heart-specific nkx2.5 gene was markedly down-regulated in SAP30L morphants, and the cardiac phenotype could be partially rescued by nkx2.5 mRNA. In addition, changes were detected in the expression of genes known to be important in hemoglobin synthesis and erythropoiesis. Our results demonstrate that SAP30L regulates several transcriptional pathways in zebrafish embryos and is involved in the development of cardiac and hematopoietic systems.

Characterization of a transcriptome from a non-model organism, Cladonia rangiferina, the grey reindeer lichen, using high-throughput next generation sequencing and EST sequence data

BACKGROUND

Lichens are symbiotic organisms that have a remarkable ability to survive in some of the most extreme terrestrial climates on earth. Lichens can endure frequent desiccation and wetting cycles and are able to survive in a dehydrated molecular dormant state for decades at a time. Genetic resources have been established in lichen species for the study of molecular systematics and their taxonomic classification. No lichen species have been characterised yet using genomics and the molecular mechanisms underlying the lichen symbiosis and the fundamentals of desiccation tolerance remain undescribed. We report the characterisation of a transcriptome of the grey reindeer lichen, Cladonia rangiferina, using high-throughput next-generation transcriptome sequencing and traditional Sanger EST sequencing data.

RESULTS

Altogether 243,729 high quality sequence reads were de novo assembled into 16,204 contigs and 49,587 singletons. The genome of origin for the sequences produced was predicted using Eclat with sequences derived from the axenically grown symbiotic partners used as training sequences for the classification model. 62.8% of the sequences were classified as being of fungal origin while the remaining 37.2% were predicted as being of algal origin. The assembled sequences were annotated by BLASTX comparison against a non-redundant protein sequence database with 34.4% of the sequences having a BLAST match. 29.3% of the sequences had a Gene Ontology term match and 27.9% of the sequences had a domain or structural match following an InterPro search. 60 KEGG pathways with more than 10 associated sequences were identified.

CONCLUSIONS

Our results present a first transcriptome sequencing and de novo assembly for a lichen species and describe the ongoing molecular processes and the most active pathways in C. rangiferina. This brings a meaningful contribution to publicly available lichen sequence information. These data provide a first glimpse into the molecular nature of the lichen symbiosis and characterise the transcriptional space of this remarkable organism. These data will also enable further studies aimed at deciphering the genetic mechanisms behind lichen desiccation tolerance.

Alternate pathways for Bcl6-mediated regulation of B cell to plasma cell differentiation

The transcription factor Bcl6 regulates germinal center formation and differentiation of B cells into high-affinity antibody-producing plasma cells. The direct double-negative regulatory circuit between Bcl6 and Blimp-1 is well established. We now reveal alternative mechanisms for Bcl6-mediated regulation of B-cell differentiation to plasma cells and show with DT40 cells that Bcl6 directly promotes the expression of Bach2, a known suppressor of Blimp-1. Moreover, Bcl6 suppresses Blimp-1 expression through direct binding to the IRF4 gene, as well as by promoting the expression of MITF, a known suppressor of IRF4. We also provide evidence that Bcl6 is needed for the expression of AID and UNG, the indispensable proteins for somatic hypermutation and class-switch recombination, and UNG appears to be a direct Bcl6 target. Our findings reveal a complex regulatory network in which Bcl6 acts as a key element dictating the transition of DT40 B cells to plasma cells.

Optimization and comparison of different methods for RNA isolation for cDNA library construction from the reindeer lichen Cladonia rangiferina

BACKGROUND

The reindeer lichen is the product of a mutualistic relationship between a fungus and an algae. Lichen demonstrate a remarkable capacity to tolerate dehydration. This tolerance is driven by a variety of biochemical processes and the accumulation of specific secondary metabolites that may be of relevance to the pharmaceutical, biotechnology and agriculture industries. These protective metabolites hinder in vitro enzymatic reactions required in cDNA synthesis. Along with the low concentrations of RNA present within lichen tissues, the process of creating a cDNA library is technically challenging.

FINDINGS

An evaluation of existing commercial and published protocols for RNA extraction from plant or fungal tissues has been performed and experimental conditions have been optimised to balance the need for the highest quality total ribonucleotides and the constraints of budget, time and human resources.

CONCLUSION

We present a protocol that balances inexpensive RNA extraction methods with commercial RNA clean-up kits to yield sufficient RNA for cDNA library construction. Evaluation of the protocol and the construction of, and sampling from, a cDNA library is used to demonstrate the suitability of the RNA extraction method for expressed sequence tag production.

Rapid detection and isolation of DNA-binding compounds from Streptomyces xanthochromogenes

The culture medium of Streptomyces xanthochromogenes JH903 was found to show selective activity against DNA-repair-deficient Escherichia coli CM871 strain. In this report we describe a simple method to locate and isolate DNA binding compounds from the fermentation broth. The method is based on the retention of DNA-reacting compounds in cellulose complexed with DNA, and purification of these compounds with thin-layer chromatography. Screening of microbial metabolites from chloroform extracts of fermentation media resulted in detection of five genotoxic fractions.