Developmental roles of glomerular epithelial protein-1 in mice molar morphogenesis.

Glomerular epithelial protein-1 (Glepp1), a R3 subtype family of receptor-type protein tyrosine phosphatases, plays important role in the activation of Src family kinases and regulates cellular processes such as cell proliferation, differentiation, and apoptosis. In this study, we firstly examined the functional evaluation of Glepp1 in tooth development and morphogenesis. The precise expression level and developmental function of Glepp1 were examined by RT-qPCR, in situ hybridization, and loss and gain of functional study using a range of in vitro organ cultivation methods. Expression of Glepp1 was detected in the developing tooth germs in cap and bell stage of tooth development. Knocking down Glepp1 at E13 for 2 days showed the altered expression levels of tooth development-related signaling molecules, including Bmps, Dspp, Fgf4, Lef1, and Shh. Moreover, transient knock down of Glepp1 revealed alterations in cellular physiology, examined by the localization patterns of Ki67 and E-cadherin. Similarly, knocking down of Glepp1 showed disrupted enamel rod and interrod formation in 3-week renal transplanted teeth. In addition, due to attrition of odontoblastic layers, the expression signals of Dspp and the localization of NESTIN were almost not detected after knock down of Glepp1; however, their expressions were increased after Glepp1 overexpression. Thus, our results suggested that Glepp1 plays modulating roles during odontogenesis by regulating the expression levels of signaling molecules and cellular events to achieve the proper structural formation of hard tissue matrices in mice molar development.

Comparability of reference-based and reference-free transcriptome analysis approaches at the gene expression level.

BACKGROUND: Lately, high-throughput RNA sequencing has been extensively used to elucidate the transcriptome landscape and dynamics of cell types of different species. In particular, for most non-model organisms lacking complete reference genomes with high-quality annotation of genetic information, reference-free (RF) de novo transcriptome analyses, rather than reference-based (RB) approaches, are widely used, and RF analyses have substantially contributed toward understanding the mechanisms regulating key biological processes and functions. To date, numerous bioinformatics studies have been conducted for assessing the workflow, production rate, and completeness of transcriptome assemblies within and between RF and RB datasets. However, the degree of consistency and variability of results obtained by analyzing gene expression levels through these two different approaches have not been adequately documented. RESULTS: In the present study, we evaluated the differences in expression profiles obtained with RF and RB approaches and revealed that the former tends to be satisfactorily replaced by the latter with respect to transcriptome repertoires, as well as from a gene expression quantification perspective. In addition, we urge cautious interpretation of these findings. Several genes that are lowly expressed, have long coding sequences, or belong to large gene families must be validated carefully, whenever gene expression levels are calculated using the RF method. CONCLUSIONS: Our empirical results indicate important contributions toward addressing transcriptome-related biological questions in non-model organisms.

PIC-Me: paralogs and isoforms classifier based on machine-learning approaches.

BACKGROUND: Paralogs formed through gene duplication and isoforms formed through alternative splicing have been important processes for increasing protein diversity and maintaining cellular homeostasis. Despite their recognized importance and the advent of large-scale genomic and transcriptomic analyses, paradoxically, accurate annotations of all gene loci to allow the identification of paralogs and isoforms remain surprisingly incomplete. In particular, the global analysis of the transcriptome of a non-model organism for which there is no reference genome is especially challenging. RESULTS: To reliably discriminate between the paralogs and isoforms in RNA-seq data, we redefined the pre-existing sequence features (sequence similarity, inverse count of consecutive identical or non-identical blocks, and match-mismatch fraction) previously derived from full-length cDNAs and EST sequences and described newly discovered genomic and transcriptomic features (twilight zone of protein sequence alignment and expression level difference). In addition, the effectiveness and relevance of the proposed features were verified with two widely used support vector machine (SVM) and random forest (RF) models. From nine RNA-seq datasets, all AUC (area under the curve) scores of ROC (receiver operating characteristic) curves were over 0.9 in the RF model and significantly higher than those in the SVM model. CONCLUSIONS: In this study, using an RF model with five proposed RNA-seq features, we implemented our method called Paralogs and Isoforms Classifier based on Machine-learning approaches (PIC-Me) and showed that it outperformed an existing method. Finally, we envision that our tool will be a valuable computational resource for the genomics community to help with gene annotation and will aid in comparative transcriptomics and evolutionary genomics studies, especially those on non-model organisms.

Identification of Differentially Expressed Genes Associated with Extracellular Matrix Degradation and Inflammatory Regulation in Calcific Tendinopathy Using RNA Sequencing.

Calcific tendinopathy (CT), developed due to calcium hydroxyapatite deposition in the rotator cuff tendon, mostly affects women in their 40 s and 50 s and causes severe shoulder pain. However, the molecular basis of its pathogenesis and appropriate treatment methods are largely unknown. In this study, we identified 202 differentially expressed genes (DEGs) between calcific and adjacent normal tendon tissues of rotator cuff using RNA sequencing-based transcriptome analysis. The DEGs were highly enriched in extracellular matrix (ECM) degradation and inflammation-related processes. Further, matrix metalloproteinase 9 (MMP9) and matrix metalloproteinase 13 (MMP13), two of the enzymes associated with ECM degradation, were found to be highly upregulated 25.85- and 19.40-fold, respectively, in the calcific tendon tissues compared to the adjacent normal tendon tissues. Histopathological analyses indicated collagen degradation and macrophage infiltration at the sites of calcific deposit in the rotator cuff tendon. Our study acts as a foundation that may help in better understanding of the pathogenesis associated with CT, and thus in better management of the disease.

Accelerated temperature and humidity testing of 2D SnS2 thin films made via four-inch-wafer-scale atomic layer deposition.

Tin disulfide (SnS2) has emerged as a promising two-dimensional (2D) material due to its excellent electrical and optical properties. However, research into 2D SnS2 has mainly focused on its synthesis procedures and applications; its stability to humidity and temperature has yet to be studied. In this work, 2D SnS2 thin films were grown by atomic layer deposition (ALD) and characterized by various tools, such as x-ray diffraction, Raman analysis, and transmission electron spectroscopy. Characterization reveals that ALD-grown SnS2 thin films are a high-quality 2D material. After characterization, a four-inch-wafer-scale uniformity test was performed by Raman analysis. Owing to the quality, large-area growth enabled by the ALD process, 98.72% uniformity was obtained. Finally, we calculated the thermodynamic equations for possible reactions between SnS2 and H2O to theoretically presurmise the oxidation of SnS2 during accelerated humidity and temperature testing. After the accelerated humidity and temperature test, x-ray diffraction, Raman analysis, and Auger electron spectroscopy were performed to check whether SnS2 was oxidized or not. Our data revealed that 2D SnS2 thin films were stable at humid conditions.

COVID-19 Patients Upregulate Toll-like Receptor 4-mediated Inflammatory Signaling That Mimics Bacterial Sepsis.

BACKGROUND: Observational studies of the ongoing coronavirus disease 2019 (COVID-19) outbreak suggest that a 'cytokine storm' is involved in the pathogenesis of severe illness. However, the molecular mechanisms underlying the altered pathological inflammation in COVID-19 are largely unknown. We report here that toll-like receptor (TLR) 4-mediated inflammatory signaling molecules are upregulated in peripheral blood mononuclear cells (PBMCs) from COVID-19 patients, compared with healthy controls (HC). METHODS: A total of 48 subjects including 28 COVID-19 patients (8 severe/critical vs. 20 mild/moderate cases) admitted to Chungnam National University Hospital, and age/sex-matched 20 HC were enrolled in this study. PBMCs from the subjects were processed for nCounter Human Immunology gene expression assay to analyze the immune related transcriptome profiles. Recombinant proteins of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) were used to stimulate the PBMCs and monocyte-derived macrophages, and real-time polymerase chain reaction was performed to quantify the mRNA expressions of the pro-inflammatory cytokines/chemokines. RESULTS: Among the most highly increased inflammatory mediators in severe/critically ill patients, S100A9, an alarmin and TLR4 ligand, was found as a noteworthy biomarker, because it inversely correlated with the serum albumin levels. We also observed that recombinant S2 and nucleocapsid proteins of SARS-CoV-2 significantly increased pro-inflammatory cytokines/chemokines and S100A9 in human primary PBMCs. CONCLUSION: These data support a link between TLR4 signaling and pathological inflammation during COVID-19 and contribute to develop therapeutic approaches through targeting TLR4-mediated inflammation.

Conventionally used reference genes are not outstanding for normalization of gene expression in human cancer research.

BACKGROUND: The selection of reference genes is essential for quantifying gene expression. Theoretically they should be expressed stably and not regulated by experimental or pathological conditions. However, identification and validation of reference genes for human cancer research are still being regarded as a critical point, because cancerous tissues often represent genetic instability and heterogeneity. Recent pan-cancer studies have demonstrated the importance of the appropriate selection of reference genes for use as internal controls for the normalization of gene expression; however, no stably expressed, consensus reference genes valid for a range of different human cancers have yet been identified. RESULTS: In the present study, we used large-scale cancer gene expression datasets from The Cancer Genome Atlas (TCGA) database, which contains 10,028 (9,364 cancerous and 664 normal) samples from 32 different cancer types, to confirm that the expression of the most commonly used reference genes is not consistent across a range of cancer types. Furthermore, we identified 38 novel candidate reference genes for the normalization of gene expression, independent of cancer type. These genes were found to be highly expressed and highly connected to relevant gene networks, and to be enriched in transcription-translation regulation processes. The expression stability of the newly identified reference genes across 29 cancerous and matched normal tissues were validated via quantitative reverse transcription PCR (RT-qPCR). CONCLUSIONS: We reveal that most commonly used reference genes in current cancer studies cannot be appropriate to serve as representative control genes for quantifying cancer-related gene expression levels, and propose in this study three potential reference genes (HNRNPL, PCBP1, and RER1) to be the most stably expressed across various cancerous and normal human tissues.

Identification and sequencing of the gene encoding DNA methyltransferase 3 (DNMT3) from sea cucumber, Apostichopus japonicus.

The sea cucumber Apostichopus japonicus is well known as a traditional tonic food and as a commercially important cultured aquatic species. This species produces saponins, and has a high potential to cope with environmental stress, such as aestivation, organ regeneration, and wound healing. Recently, several studies have shown that cellular reprogramming and the physiological responses of the sea cucumber to environmental changes, including aestivation, are potentially mediated by epigenetic DNA methylation. The DNA methyltransferase (DNMT)1 and DNMT3 genes are independent participants in the maintenance and de novo methylation of specific sequences. Sea urchin (Strongylocentrotus purpuratus) and starfish (Asterina pectinifera), which belong to the same phylum as A. japonicus, have both DNMT1 and DNMT3 genes. However, it was previously reported that DNMT1 is present, but DNMT3 is absent, in A. japonicus. In the present study, we sequenced the full-length cDNA of the A. japonicus DNMT3 gene. The newly sequenced DNMT3 gene comprises three major conserved domains (Pro-Trp-Trp-Pro (PWWP), plant homeodomain (PHD), and S-adenosylmethionine-dependent methyltransferase (AdoMet-MTase)), indicating that the DNMT3 possibly has de novo DNA methylation catalytic activity. Gene structure and phylogenetic analysis showed that sea cucumber DNMT3 is evolutionarily conserved in the Echinodermata. Next, we demonstrated the conservation of DNMT3 gene expression in sea cucumber and starfish belong to same phylum, echinoderm. Using reverse transcription-polymerase chain reaction, sea cucumber DNMT3 mRNA was detected in testis tissue, but not in other tissues tested, including the respiratory tree, muscle, tentacle, intestine, and ovary. This is inconsistent with previous reports, which showed the expression of DNMT3 in ovary, but not in testis of the starfish A. pectinifera, indicating the tissue- and species-specific expression of DNMT3 gene. Although further studies are needed to clarify the epigenetic regulatory mechanisms of DNMT3 and its application to the aquaculture industry, our findings may provide insights into the sea cucumber biology.

Variant- and vaccination-specific alternative splicing profiles in SARS-CoV-2 infections.

The COVID-19 pandemic, driven by the SARS-CoV-2 virus and its variants, highlights the important role of understanding host-viral molecular interactions influencing infection outcomes. Alternative splicing post-infection can impact both host responses and viral replication. We analyzed RNA splicing patterns in immune cells across various SARS-CoV-2 variants, considering immunization status. Using a dataset of 190 RNA-seq samples from our prior studies, we observed a substantial deactivation of alternative splicing and RNA splicing-related genes in COVID-19 patients. The alterations varied significantly depending on the infecting variant and immunization history. Notably, Alpha or Beta-infected patients differed from controls, while Omicron-infected patients displayed a splicing profile closer to controls. Particularly, vaccinated Omicron-infected individuals showed a distinct dynamic in alternative splicing patterns not widely shared among other groups. Our findings underscore the intricate interplay between SARS-CoV-2 variants, vaccination-induced immunity, and alternative splicing, emphasizing the need for further investigations to deepen understanding and guide therapeutic development.

The Loss of an Orphan Nuclear Receptor NR2E3 Augments Wnt/ß-catenin Signaling via Epigenetic Dysregulation that Enhances Sp1-ß catenin-p300 Interactions in Hepatocellular Carcinoma.

The orphan nuclear receptor NR2E3 (Nuclear receptor subfamily 2 group E, Member 3) is an epigenetic player that modulates chromatin accessibility to activate p53 during liver injury. Nonetheless, a precise tumor suppressive and epigenetic role of NR2E3 in hepatocellular carcinoma (HCC) development remains unclear. HCC patients expressing low NR2E3 exhibit unfavorable clinical outcomes, aligning with heightened activation of the Wnt/ß-catenin signaling pathway. The murine HCC models utilizing NR2E3 knockout mice consistently exhibits accelerated liver tumor formation accompanied by enhanced activation of Wnt/ß-catenin signaling pathway and inactivation of p53 signaling. At cellular level, the loss of NR2E3 increases the acquisition of aggressive cancer cell phenotype and tumorigenicity and upregulates key genes in the WNT/ß-catenin pathway with increased chromatin accessibility. This event is mediated through increased formation of active transcription complex involving Sp1, ß-catenin, and p300, a histone acetyltransferase, on the promoters of target genes. These findings demonstrate that the loss of NR2E3 activates Wnt/ß-catenin signaling at cellular and organism levels and this dysregulation is associated with aggressive HCC development and poor clinical outcomes. In summary, NR2E3 is a novel tumor suppressor with a significant prognostic value, maintaining epigenetic homeostasis to suppress the Wnt/ß-catenin signaling pathway that promotes HCC development.

Variant- and Vaccination-Specific Alternative Splicing Profiles in SARS-CoV-2 Infections.

The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, and its subsequent variants has underscored the importance of understanding the host-viral molecular interactions to devise effective therapeutic strategies. A significant aspect of these interactions is the role of alternative splicing in modulating host responses and viral replication mechanisms. Our study sought to delineate the patterns of alternative splicing of RNAs from immune cells across different SARS-CoV-2 variants and vaccination statuses, utilizing a robust dataset of 190 RNA-seq samples from our previous studies, encompassing an average of 212 million reads per sample. We identified a dynamic alteration in alternative splicing and genes related to RNA splicing were highly deactivated in COVID-19 patients and showed variant- and vaccination-specific expression profiles. Overall, Omicron-infected patients exhibited a gene expression profile akin to healthy controls, unlike the Alpha or Beta variants. However, significantly, we found identified a subset of infected individuals, most pronounced in vaccinated patients infected with Omicron variant, that exhibited a specific dynamic in their alternative splicing patterns that was not widely shared amongst the other groups. Our findings underscore the complex interplay between SARS-CoV-2 variants, vaccination-induced immune responses, and alternative splicing, emphasizing the necessity for further investigations into these molecular cross-talks to foster deeper understanding and guide strategic therapeutic development.

Network medicine-based epistasis detection in complex diseases: ready for quantum computing.

Most heritable diseases are polygenic. To comprehend the underlying genetic architecture, it is crucial to discover the clinically relevant epistatic interactions (EIs) between genomic single nucleotide polymorphisms (SNPs)1-3. Existing statistical computational methods for EI detection are mostly limited to pairs of SNPs due to the combinatorial explosion of higher-order EIs. With NeEDL (network-based epistasis detection via local search), we leverage network medicine to inform the selection of EIs that are an order of magnitude more statistically significant compared to existing tools and consist, on average, of five SNPs. We further show that this computationally demanding task can be substantially accelerated once quantum computing hardware becomes available. We apply NeEDL to eight different diseases and discover genes (affected by EIs of SNPs) that are partly known to affect the disease, additionally, these results are reproducible across independent cohorts. EIs for these eight diseases can be interactively explored in the Epistasis Disease Atlas (https://epistasis-disease-atlas.com). In summary, NeEDL is the first application that demonstrates the potential of seamlessly integrated quantum computing techniques to accelerate biomedical research. Our network medicine approach detects higher-order EIs with unprecedented statistical and biological evidence, yielding unique insights into polygenic diseases and providing a basis for the development of improved risk scores and combination therapies.

Head transcriptome profiling of glossiphoniid leech (Helobdella austinensis) reveals clues about proboscis development.

Cephalization refers to the evolutionary trend towards the concentration of neural tissues, sensory organs, mouth and associated structures at the front end of bilaterian animals. Comprehensive studies on gene expression related to the anterior formation in invertebrate models are currently lacking. In this study, we performed de novo transcriptional profiling on a proboscis-bearing leech (Helobdella austinensis) to identify differentially expressed genes (DEGs) in the anterior versus other parts of the body, in particular to find clues as to the development of the proboscis. Between the head and the body, 132 head-specific DEGs were identified, of which we chose 11 to investigate their developmental function during embryogenesis. Analysis of the spatial expression of these genes using in situ hybridization showed that they were characteristically expressed in the anterior region of the developing embryo, including the proboscis. Our results provide information on the genes related to head formation and insights into the function of proboscis-related genes during organogenesis with the potential roles of genes not yet characterized.

Dysregulated thrombospondin 1 and miRNA-29a-3p in severe COVID-19.

Although nearly a fifth of symptomatic COVID-19 patients suffers from severe pulmonary inflammation, the mechanism of developing severe illness is not yet fully understood. To identify significantly altered genes in severe COVID-19, we generated messenger RNA and micro-RNA profiling data of peripheral blood mononuclear cells (PBMCs) from five COVID-19 patients (2 severe and 3 mild patients) and three healthy controls (HC). For further evaluation, two publicly available RNA-Seq datasets (GSE157103 and GSE152418) and one single-cell RNA-Seq dataset (GSE174072) were employed. Based on RNA-Seq datasets, thrombospondin 1 (THBS1) and interleukin-17 receptor A (IL17RA) were significantly upregulated in severe COVID-19 patients' blood. From single-cell RNA-sequencing data, IL17RA level is increased in monocytes and neutrophils, whereas THBS1 level is mainly increased in the platelets. Moreover, we identified three differentially expressed microRNAs in severe COVID-19 using micro-RNA sequencings. Intriguingly, hsa-miR-29a-3p significantly downregulated in severe COVID-19 was predicted to bind the 3'-untranslated regions of both IL17RA and THBS1 mRNAs. Further validation analysis of our cohort (8 HC, 7 severe and 8 mild patients) showed that THBS1, but not IL17RA, was significantly upregulated, whereas hsa-miR-29a-3p was downregulated, in PBMCs from severe patients. These findings strongly suggest that dysregulated expression of THBS1, IL17RA, and hsa-miR-29a-3p involves severe COVID-19.

ppGpp signaling plays a critical role in virulence of Acinetobacter baumannii.

Acinetobacter baumannii, a major nosocomial pathogen, survives in diverse hospital environments, and its multidrug resistance is a major concern. The ppGpp-dependent stringent response mediates the reprogramming of genes with diverse functions in several bacteria. We investigated whether ppGpp is involved in A. baumannii's pathogenesis by examining biofilm formation, surface motility, adhesion, invasion, and mouse infection studies. Transcriptome analysis of early stationary phase cultures revealed 498 differentially-expressed genes (≥ 2-fold change) in a ppGpp-deficient A. baumannii strain; 220 and 278 genes were up and downregulated, respectively. Csu operon expression, important in pilus biosynthesis during early biofilm formation, was significantly reduced in the ppGpp-deficient strain. Our findings suggest that ppGpp signaling influences A. baumannii biofilm formation, surface motility, adherence, and virulence. We showed the association between ppGpp and pathogenicity in A. baumannii for the first time; ppGpp may be a novel antivirulence target in A. baumannii.

MiR-144-3p is associated with pathological inflammation in patients infected with Mycobacteroides abscessus.

Infection with rapidly growing nontuberculous mycobacteria is emerging as a global health issue; however, key host factors remain elusive. Here, we investigated the characteristic immune profiles of peripheral blood mononuclear cells (PBMCs) from patients infected with Mycobacteroides abscessus subsp. abscessus (Mabc) and M. abscessus subsp. massiliense (Mmass). Using an integrated analysis of global mRNA and microRNA expression profiles, we found that several inflammatory cytokines/chemokines [interleukin (IL)-1ß, IL-6, C-X-C motif chemokine ligand 2, and C-C motif chemokine ligand 2] and miR-144-3p were significantly upregulated in PBMCs from patients compared with those from healthy controls (HCs). Notably, there was a strong correlation between the expression levels of miR-144-3p and proinflammatory cytokines/chemokines. Similarly, upregulated expression of miR-144-3p and proinflammatory cytokines/chemokines was found in macrophages and lungs from mice after infection with Mabc and Mmass. We showed that the expression of negative regulators of inflammation (SARM1 and TNIP3) was significantly downregulated in PBMCs from the patients, although they were not putative targets of miR-144-3p. Furthermore, overexpression of miR-144-3p led to a marked increase in proinflammatory cytokines/chemokines and promoted bacterial growth in macrophages. Together, our results highlight the importance of miR-144-3p linking to pathological inflammation during M. abscessus infection.

Lipocalin2 Induced by Bacterial Flagellin Protects Mice against Cyclophosphamide Mediated Neutropenic Sepsis.

Neutropenic sepsis is a fatal consequence of chemotherapy, and septic complications are the principal cause of mortality. Chemotherapy-induced neutropenia leads to the formation of microscopic ulcers in the gastrointestinal epithelium that function as a portal of entry for intraluminal bacteria, which translocate across the intestinal mucosal barrier and gain access to systemic sites, causing septicemia. A cyclophosphamide-induced mouse model was developed to mimic the pathophysiologic sequence of events that occurs in patients with neutropenic sepsis. The TLR5 agonist bacterial flagellin derived from Vibrio vulnificus extended the survival of cyclophosphamide-treated mice by reducing the bacterial load in internal organs. The protective effect of flagellin was mediated by the antimicrobial protein lipocalin 2 (Lcn2), which is induced by TLR5-NF-κB activation in hepatocytes. Lcn2 sequestered iron from infecting bacteria, particularly siderophore enterobactin-dependent members of the Enterobacteriaceae family, thereby limiting their proliferation. Lcn2 should be considered for the treatment of neutropenic sepsis and gastrointestinal damage during chemotherapy to prevent or minimize the adverse effects of cancer chemotherapy.

RBM47-regulated alternative splicing of TJP1 promotes actin stress fiber assembly during epithelial-to-mesenchymal transition.

Morphological and functional changes in cells during the epithelial-mesenchymal transition (EMT) process are known to be regulated by alternative splicing. However, only a few splicing factors involved in EMT have been reported and their underlying mechanisms remain largely unknown. Here, we showed that an isoform of tight junction protein 1 (TJP1) lacking exon 20 (TJP1-α-) is predominantly expressed in tumor tissues and in A549 cells during transforming growth factor-ß (TGF-ß)-induced EMT. RBM47 promoted the inclusion of exon 20 of TJP1, the alternative exon encoding the α-domain, by which RBM47 recognizes to (U)GCAUG in the downstream intronic region of exon 20. We also found that the first RNA recognition motif (RRM) domain of RBM47 is critical in the regulation of alternative splicing and its recognition to pre-mRNA of TJP1. Furthermore, we demonstrated that the TJP1-α- isoform enhances the assembly of actin stress fibers, thereby promoting cellular migration in a wound healing assay. Our results suggest the regulatory mechanism for the alternative splicing of TJP1 pre-mRNA by RBM47 during EMT, providing a basis for studies related to the modulation of EMT via alternative splicing.

Characterisation of genes differentially expressed in macrophages by virulent and attenuated Mycobacterium tuberculosis through RNA-Seq analysis.

Tuberculosis (TB) remains a global healthcare issue. Understanding the host-pathogen interactions in TB is vital to develop strategies and therapeutic tools for the control of Mycobacterium tuberculosis (Mtb). In this study, transcriptome analyses of macrophages infected with either the virulent Mtb strain H37Rv (Rv) or the avirulent Mtb strain H37Ra (Ra) were carried out and 750 differentially expressed genes (DEGs) were identified. As expected, the DEGs were mainly involved in the induction of innate immune responses against mycobacterial infections. Among the DEGs, solute carrier family 7 member 2 (Slc7a2) was more strongly expressed in Ra-infected macrophages. Induction of SLC7A2 was important for macrophages to control the intracellular survival of Mtb. Our results imply that SLC7A2 plays an important role in macrophages during Mtb infection. Our findings could prove useful for the development of new therapeutic strategies to control TB infection.

The developmental transcriptome atlas of the spoon worm Urechis unicinctus (Echiurida: Annelida).

Background: Echiurida is one of the most intriguing major subgroups of annelida because, unlike most other annelids, echiurids lack metameric body segmentation as adults. For this reason, transcriptome analyses from various developmental stages of echiurid species can be of substantial value for understanding precise expression levels and the complex regulatory networks during early and larval development. Results: A total of 914 million raw RNA-Seq reads were produced from 14 developmental stages of Urechis unicinctus and were de novo assembled into contigs spanning 63,928,225 bp with an N50 length of 2700 bp. The resulting comprehensive transcriptome database of the early developmental stages of U. unicinctus consists of 20,305 representative functional protein-coding transcripts. Approximately 66% of unigenes were assigned to superphylum-level taxa, including Lophotrochozoa (40%). The completeness of the transcriptome assembly was assessed using benchmarking universal single-copy orthologs; 75.7% of the single-copy orthologs were presented in our transcriptome database. We observed 3 distinct patterns of global transcriptome profiles from 14 developmental stages and identified 12,705 genes that showed dynamic regulation patterns during the differentiation and maturation of U. unicinctus cells. Conclusions: We present the first large-scale developmental transcriptome dataset of U. unicinctus and provide a general overview of the dynamics of global gene expression changes during its early developmental stages. The analysis of time-course gene expression data is a first step toward understanding the complex developmental gene regulatory networks in U. unicinctus and will furnish a valuable resource for analyzing the functions of gene repertoires in various developmental phases.