Transcriptome of embryonic and neonatal mouse cortex by high-throughput RNA sequencing

Prenatal progesterone treatment modulates fetal brain transcriptome and impacts adult offspring behavior in mice

Maternal exposure to elevated levels of steroid hormones during pregnancy is associated with the development of chronic conditions in offspring that manifest in adulthood. However, the effects of progesterone (P4) administration during early pregnancy on fetal development and subsequent offspring behavior remain poorly understood. In this study, we aimed to investigate the effects of P4 treatment during early pregnancy on the transcript abundance in the fetal brain and assess the behavioral consequences in the offspring during adolescence and adulthood. Using RNA-seq analysis, we examined the impact of P4 treatment on the fetal brain transcriptome in a dosage-dependent manner. Our results revealed differential regulation of genes involved in neurotransmitter transport, synaptic transmission, and transcriptional regulation. Specifically, we observed bidirectional regulation of transcription factors (TFs) by P4 at different doses, highlighting the critical role of these TFs in neurodevelopment. To assess behavioral outcomes, we conducted open field and elevated plus maze tests. Offspring treated with low-dose P4 (LP4) displayed increased exploratory behavior during both adolescence and adulthood. In contrast, the high-dose P4 (HP4) group exhibited impaired exploration and heightened anxiety-like behaviors compared to the control mice. Moreover, in a novel object recognition test, HP4-treated offspring demonstrated impaired object recognition memory during both developmental stages. Additionally, both LP4 and HP4 groups showed reduced social interaction in the three-chamber test. These results suggest that prenatal exposure to P4 exerts a notable influence on the expression of genes associated with neurodevelopment and may induce alterations in behavioral characteristics in progeny, highlighting the need to monitor progesterone levels during pregnancy for long-term impacts on fetal brain development and behavior.

New insights into the regulation mechanism of red claw crayfish (Cherax quadricarinatus) hepatopancreas under air exposure using transcriptome analysis

Red claw crayfish (Cherax quadricarinatus) is an important freshwater shrimp species worldwide with enormous economic value. Waterless transportation is an inherent feature of red claw crayfish transportation. However, the high mortality of red claw crayfish is a severe problem in the aquaculture of crayfish after waterless transportation. In this study, we investigated the responses of the hepatopancreas from the red claw crayfish undergoing air exposure stress and normal conditions on transcriptome levels. We used Illumina-based RNA sequencing (RNA-Seq) to perform a transcriptome analysis from the hepatopancreas of red claw crayfish challenged by air exposure. An average of 57,148,800 clean reads per library was obtained, and 33,567 unigenes could be predicted and classified according to their homology with matches in the National Center for Biotechnology Information (NCBI) non-redundant protein sequences (Nr), Gene Ontology (GO), a manually annotated and reviewed protein sequence database (Swiss-Prot), protein families (Pfam), Clusters of Orthologous Groups (COG) of proteins, and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. 690 and 3407 differentially expressed genes (DEGs) were identified between the two stress stages of the red claw crayfish. More DEGs were identified in 12 h, indicating that gene expressions were largely changed at 12 h. Some immune-related pathways and genes were identified according to KEGG and GO enrichment analysis. A total of 12 DEGs involved in immune response and trehalose mechanism were verified by quantitative real-time-polymerase chain reaction (qRT-PCR). The results indicated that the red claw crayfish might counteract the stress of air exposure at the transcriptomic level by increasing expression levels of antioxidant-, immune-, and trehalose metabolism-related genes. These transcriptome results from the hepatopancreas provide significant insights into the influence mechanism of air exposure to the trehalose mechanism and immune response in the red claw crayfish.

Integrated gene co-expression network analysis reveals unique developmental processes of Aurelia aurita

The typical life cycle of the moon jellyfish (Aurelia aurita) includes the planula, polyp, strobila, ephyra, and medusa developmental stages. These stages exhibit huge differences in both external morphology and internal physiological functions. However, the gene co-expression network involved in these post-embryonic developmental processes has not been studied yet. Here, based on 15 RNA sequencing samples covering all five stages of the A. aurita life cycle, we systematically analyzed the gene co-expression network and obtained 35 relevant modules. Furthermore, we identified the highly correlated modules and hub genes for each stage. These hub genes are implicated to play important roles in the developmental processes of A. aurita, which should help improve our understanding of the jellyfish life cycle.

Nutrient sensitive protein O-GlcNAcylation modulates the transcriptome through epigenetic mechanisms during embryonic neurogenesis

Protein O-GlcNAcylation is a dynamic, nutrient-sensitive mono-glycosylation deposited on numerous nucleo-cytoplasmic and mitochondrial proteins, including transcription factors, epigenetic regulators, and histones. However, the role of protein O-GlcNAcylation on epigenome regulation in response to nutrient perturbations during development is not well understood. Herein we recapitulated early human embryonic neurogenesis in cell culture and found that pharmacological up-regulation of O-GlcNAc levels during human embryonic stem cells' neuronal differentiation leads to up-regulation of key neurogenic transcription factor genes. This transcriptional de-repression is associated with reduced H3K27me3 and increased H3K4me3 levels on the promoters of these genes, perturbing promoter bivalency possibly through increased EZH2-Thr311 phosphorylation. Elevated O-GlcNAc levels also lead to increased Pol II-Ser5 phosphorylation and affect H2BS112O-GlcNAc and H2BK120Ub1 on promoters. Using an in vivo rat model of maternal hyperglycemia, we show similarly elevated O-GlcNAc levels and epigenetic dysregulations in the developing embryo brains because of hyperglycemia, whereas pharmacological inhibition of O-GlcNAc transferase (OGT) restored these molecular changes. Together, our results demonstrate O-GlcNAc mediated sensitivity of chromatin to nutrient status, and indicate how metabolic perturbations could affect gene expression during neurodevelopment.

The toggle switch model for gene expression change during the prenatal-to-postnatal transition in mammals

The prenatal-to-postnatal transition is a pivotal process in the life cycle whereby an organism shifts from responding to intrauterine cues to undergoing extrauterine stresses with many physiological adaptations. However, the molecular basis underlying the evolutionarily conserved physiological adaptations remains elusive. Here, we analyze the transcriptomes of seven organs across developmental time points from five mammalian species by constructing computational coexpression networks and report a developmental shift of gene expression at the perinatal stage. The low-to-high and high-to-low expressed genes tightly coalesce in the functional categories and gene regulatory pathways that implicate the physiological adaptions during the prenatal-to-postnatal transition, including lipid metabolism, circadian rhythm, immune response, cell cycle, and cell division. The low-to-high and high-to-low expressed genes around the perinatal stage tend to form the mutually inhibitory toggle switch gene pairs linking the gene regulatory networks in response to the environmental changes. We thus propose the toggle switch model for the developmental shift of gene expression as a mechanic framework to investigate how the physiological adaptations occur during the prenatal-to-postnatal transition.

Pathogenic Infections during Pregnancy and the Consequences for Fetal Brain Development

Pathogens comprised of viruses, bacteria, gut microbiome, and parasites are a leading cause of ever-emerging diseases in humans. Studying pathogens for their ability to cause diseases is a topic of critical discussion among scientists and pharmaceutical centers for effective drug development that diagnose, treat, and prevent infection-associated disorders. Pathogens impact health either directly by invading the host or by eliciting an acute inflammatory immune response. This paradigm of inflammatory immune responses is even more consequential in people who may be immunocompromised. In this regard, pregnancy offers an altered immunity scenario, which may allow the onset of severe diseases. Viruses, such as Influenza, HIV, and now SARS-CoV-2, associated with the COVID-19 pandemic, raise new concerns for maternal and fetal/neonatal health. Intrauterine bacterial and parasitic infections are also known to impact pregnancy outcomes and neonatal health. More importantly, viral and bacterial infections during pregnancy have been identified as a common contributor to fetal brain development defects. Infection-mediated inflammatory uterine immune milieu is thought to be the main trigger for causing poor fetal brain development, resulting in long-term cognitive impairments. The concept of in utero programming of childhood and adult disorders has revolutionized the field of neurodevelopment and its associated complications. Recent findings in mice and humans clearly support the idea that uterine immunity during pregnancy controls the health trajectory of the child and considerably influences the cognitive function and mental health. In this review, we focus on the in utero programming of autism spectrum disorders (ASD) and assess the effects of pathogens on the onset of ASD-like symptoms.

Maternal nematode infection upregulates expression of Th2/Treg and diapedesis related genes in the neonatal brain

Intestinal nematode infections common during pregnancy have recently been shown to have impacts that extend to their uninfected offspring including altered brain gene expression. If maternal immune signals reach the neonatal brain, they might alter neuroimmune development. We explored expression of genes associated with four distinct types of T cells (Th1, Th2, Th17, Treg) and with leukocyte transendothelial migration and endocytosis transport across the blood–brain barrier (BBB) in the postnatal brain of offspring of nematode-infected mice, through secondary analysis of a whole brain gene expression database. Th1/Th17 expression was lowered by maternal infection as evidenced by down-regulated expression of IL1β, Th1 receptors and related proteins, and of IL22 and several Th17 genes associated with immunopathology. In contrast, Th2/Treg related pathways were upregulated as shown by higher expression of IL4 and TGF-β family genes. Maternal infection also upregulated expression of pathways and integrin genes involved in transport of leukocytes in between endothelial cells but downregulated endosome vesicle formation related genes that are necessary for endocytosis of immunoglobulins across the BBB. Taken together, pup brain gene expression indicates that maternal nematode infection enhanced movement of leukocytes across the neonatal BBB and promoted a Th2/Treg environment that presumably minimizes the proinflammatory Th1 response in the pup brain.

RNA-Seq Transcriptome Analysis of the Liver and Brain of the Black Carp (Mylopharyngodon piceus) During Fasting

The black carp (Mylopharyngodon piceus) is an important carnivorous freshwater-cultured species. To understand the molecular basis underlying the response of black carp to fasting, we used RNA-Seq to analyze the liver and brain transcriptome of fasting fish. Annotation to the NCBI database identified 66,609 unigenes, of which 22,841 were classified into the Gene Ontology database and 15,925 were identified in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Comparative analysis of the expression profile between fasting and normal feeding fish revealed 13,737 differentially expressed genes (P < 0.05), of which 12,480 were found in liver tissue and 1257 were found in brain tissue. The KEGG pathway analysis showed significant differences in expression of genes involved in metabolic and immune pathways, such as the insulin signaling pathway, PI3K-Akt signaling pathway, cAMP signaling pathway, FoxO signaling pathway, AMPK signaling pathway, endocytosis, and apoptosis. Quantitative real-time PCR analysis confirmed that expression of the genes encoding the factors involved in those pathways differed between fasting and feeding fish. These results provide valuable information about the molecular response mechanism of black carp under fasting conditions.

Gene expression profiles underlying aggressive behavior in the prefrontal cortex of cattle

BACKGROUND

Aggressive behavior is an ancient and conserved trait, habitual for most animals in order to eat, protect themselves, compete for mating and defend their territories. Genetic factors have been shown to play an important role in the development of aggression both in animals and humans, displaying moderate to high heritability estimates. Although such types of behaviors have been studied in different animal models, the molecular architecture of aggressiveness remains poorly understood. This study compared gene expression profiles of 16 prefrontal cortex (PFC) samples from aggressive and non-aggressive cattle breeds: Lidia, selected for agonistic responses, and Wagyu, selected for tameness.

RESULTS

A total of 918 up-regulated and 278 down-regulated differentially expressed genes (DEG) were identified, representing above-chance overlap with genes previously identified in studies of aggression across species, as well as those implicated in recent human evolution. The functional interpretation of the up-regulated genes in the aggressive cohort revealed enrichment of pathways such as Alzheimer disease-presenilin, integrins and the ERK/MAPK signaling cascade, all implicated in the development of abnormal aggressive behaviors and neurophysiological disorders. Moreover, gonadotropins, are up-regulated as natural mechanisms enhancing aggression. Concomitantly, heterotrimeric G-protein pathways, associated with low reactivity mental states, and the GAD2 gene, a repressor of agonistic reactions associated with PFC activity, are down-regulated, promoting the development of the aggressive responses selected for in Lidia cattle. We also identified six upstream regulators, whose functional activity fits with the etiology of abnormal behavioral responses associated with aggression.

CONCLUSIONS

These transcriptional correlates of aggression, resulting, at least in part, from controlled artificial selection, can provide valuable insights into the complex architecture that underlies naturally developed agonistic behaviors. This analysis constitutes a first important step towards the identification of the genes and metabolic pathways that promote aggression in cattle and, providing a novel model species to disentangle the mechanisms underlying variability in aggressive behavior.

A reliable strategy for single-cell RNA sequencing analysis using cryoconserved primary cortical cells

BACKGROUND

The application of single-cell RNA sequencing (scRNASeq) represents a unique approach to identify hundreds to millions of cells in mammalian cortical multilayers at different stages of embryogenesis. ScRNASeq technology applied to neurological studies requires the use of fresh starting materials because standard cryopreservation methods do not guarantee high viability of cortical primary cells derived from dissected brain areas.

NEW METHOD

Here we set up and validate an innovative strategy to perform scRNASeq studies in cryopreserved primary cortical cells isolated from E15.5 mouse embryo. In order to freeze cortical primary cells, we have employed Neurostore, a medium able to guarantee high viability and cell composition of embryonic cortex after thawing.

COMPARISON WITH EXISTING METHODS

We showed for the first time the possibility to run scRNASeq experiments on primary cortical cells in an off-line set-up, ensuring cellular integrity and diversity.

RESULTS

By trypan blue assay and flow cytometry analysis, we found that Neurostore-cryopreserved cortical cells showed approximately 95 % of viability. Satisfactory RNA recovery and cDNA libraries were achieved. Transcriptome sequencing of 35,763 cryoconserved single cells yielded a robust data-set, identifying 25 cell clusters in three biological samples. Prevalence of peculiar neural populations before and after the cryopreservation-resuscitation procedure was verified by marker gene expression and immunofluorescence analysis.

CONCLUSIONS

Our findings support the evidence that frozen primary cortical cells can be successfully employed in scRNASeq experiments allowing an unprecedented flexibility in experimental procedures, such as sample preparation and subsequent processing steps performed in different locations.

Comprehensive transcriptional changes in the liver of Kanglang white minnow (Anabarilius grahami) in response to the infection of parasite Ichthyophthirius multifiliis

Simple Summary

Kanglang white minnow (KWM, Anabarilius grahami), is a typical “3E” (Endangered, Endemic and Economic) fish species in Yunnan-Guizhou Plateau. As one of the traditional “Four Famous Fishes” in Yunnan province, it has become the major local aquaculture species with increasing demand after the success of artificial breeding. However, this economically important fish is highly susceptible to the infection of a parasite ciliate, Ichthyophthirius multifiliis (Ich), during the practical procedure of artificial breeding. To examine the host immune responses to Ich, we divided the experimental fishes into three groups (including control, early-infected stage, and late-infected stage) for transcriptome sequencing to analyze the differentially expressed genes (DEGs) and immune response mechanisms.

Abstract

The notorious parasite Ichthyophthirius multifiliis (Ich) has been recorded worldwide in fish species and causes white spot disease, posing major threats and resulting in severe losses to international fish production. Extensively effective strategies for treating Ich are not available yet, and genetic mechanisms of hosts in response to the parasite are still largely unknown. In this study, we selected Kanglang white minnow (KWM, Anabarilius grahami) to examine its liver transcriptional changes after Ich infection, as white spot disease is one bottleneck problem in exploring this economically important species. We divided the experimental fishes into three groups (control, early-infected, and late-infected) to examine differentially expressed genes (DEGs). A total of 831 DEGs were identified and classified into 128 significantly enriched GO (Gene Ontology) terms and 71 significantly enriched KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. Most of these terms or pathways were functionally enriched in immunity, inflammatory response, and apoptosis, such as nucleotide-binding oligomerization domain-like (NOD-like) receptor signaling, tumor necrosis factor (TNF) signaling, interleukin-17 (IL-17) signaling, and apoptosis pathways. We also identified 178 putative antimicrobial peptides (AMPs) and AMP precursors based on our previously reported genome assembly of KWM, and revealed that the expressional patterns varied according to different types. In summary, our work reported the first comprehensive transcriptional changes in KWM in response to the exogenous infection of Ich, which would lay a solid foundation for in-depth studies on disease defense or resistant strains selection in this valuable fish.

Unveiling the transcriptome alteration of POMC neuron in diet-induced obesity

Loss of neuron homeostasis in the arcuate nucleus (ARC) is responsible for diet-induced-obesity (DIO). We previously reported that loss of Rb1 gene compromised the homeostasis of anorexigenic POMC neurons in ARC and induced obesity in mice. To evaluate the development of DIO, we propose to analyze the transcriptomic alteration of POMC neurons in mice following high fat diet (HFD) feeding. We isolated these neurons from established DIO mice and performed transcriptomic profiling using RNA-seq. In total, 1066 genes (628 upregulated and 438 downregulated) were identified as differentially expressed genes (DEGs). Pathway enrichment analysis with these DEGs further revealed that "cell cycle," "apoptosis," "chemokine signaling," and "sphingolipid metabolism" pathways were correlated with DIO development. Moreover, we validated that the pRb protein, a key regulator of "cell cycle pathway," was inactivated by phosphorylation in POMC neurons by HFD feeding. Importantly, the reversal of deregulated cell cycle by stereotaxic delivering of the unphosphorylated pRbΔP in ARC significantly meliorated the DIO. Collectively, our study provides insights into the mechanisms related to the loss of homeostasis of POMC neurons in DIO, and suggests pRb phosphorylation as a potential intervention target to treat DIO.

Prediction of Long Non-Coding RNAs Based on Deep Learning

With the rapid development of high-throughput sequencing technology, a large number of transcript sequences have been discovered, and how to identify long non-coding RNAs (lncRNAs) from transcripts is a challenging task. The identification and inclusion of lncRNAs not only can more clearly help us to understand life activities themselves, but can also help humans further explore and study the disease at the molecular level. At present, the detection of lncRNAs mainly includes two forms of calculation and experiment. Due to the limitations of bio sequencing technology and ineluctable errors in sequencing processes, the detection effect of these methods is not very satisfactory. In this paper, we constructed a deep-learning model to effectively distinguish lncRNAs from mRNAs. We used k-mer embedding vectors obtained through training the GloVe algorithm as input features and set up the deep learning framework to include a bidirectional long short-term memory model (BLSTM) layer and a convolutional neural network (CNN) layer with three additional hidden layers. By testing our model, we have found that it obtained the best values of 97.9%, 96.4% and 99.0% in F1score, accuracy and auROC, respectively, which showed better classification performance than the traditional PLEK, CNCI and CPC methods for identifying lncRNAs. We hope that our model will provide effective help in distinguishing mature mRNAs from lncRNAs, and become a potential tool to help humans understand and detect the diseases associated with lncRNAs.

Maternal Gastrointestinal Nematode Infection Up-regulates Expression of Genes Associated with Long-Term Potentiation in Perinatal Brains of Uninfected Developing Pups

Establishment of neural networks critical for memory and cognition begins during the perinatal period but studies on the impact of maternal infection are limited. Using a nematode parasite that remains in the maternal intestine, we tested our hypothesis that maternal infection during pregnancy and early lactation would alter perinatal brain gene expression, and that the anti-inflammatory nature of this parasite would promote synaptic plasticity and long-term potentiation. Brain gene expression was largely unaffected two days after birth, but in seven-day old pups, long-term potentiation and four related pathways essential for the development of synaptic plasticity, cognition and memory were up-regulated in pups of infected dams. Interestingly, our data suggest that a lowering of Th1 inflammatory processes may underscore the apparent beneficial impact of maternal intestinal infection on long-term potentiation.

Maternal immune activation impairs cognitive flexibility and alters transcription in frontal cortex

BACKGROUND

Epidemiological studies suggest that the risk of neurodevelopmental disorders such as autism spectrum disorder (ASD) and schizophrenia is increased by prenatal exposure to viral or bacterial infection during pregnancy. It is still unclear how activation of the maternal immune response interacts with underlying genetic factors to influence observed ASD phenotypes.

METHODS

The current study investigated how maternal immune activation (MIA) in mice impacts gene expression in the frontal cortex in adulthood, and how these molecular changes relate to deficits in cognitive flexibility and social behavior, and increases in repetitive behavior that are prevalent in ASD. Poly(I:C) (20 mg/kg) was administered to dams on E12.5 and offspring were tested for social approach behavior, repetitive grooming, and probabilistic reversal learning in adulthood (n = 8 vehicle; n = 9 Poly(I:C)). We employed next-generation high-throughput mRNA sequencing (RNA-seq) to comprehensively investigate the transcriptome profile in frontal cortex of adult offspring of Poly(I:C)-exposed dams.

RESULTS

Exposure to poly(I:C) during gestation impaired probabilistic reversal learning and decreased social approach in MIA offspring compared to controls. We found long-term effects of MIA on expression of 24 genes, including genes involved in glutamatergic neurotransmission, mTOR signaling and potassium ion channel activity. Correlations between gene expression and specific behavioral measures provided insight into genes that may be responsible for ASD-like behavioral alterations.

CONCLUSIONS

These findings suggest that MIA can lead to impairments in cognitive flexibility in mice similar to those exhibited in ASD individuals, and that these impairments are associated with altered gene expression in frontal cortex.

Liver transcriptome analysis of the Sparus macrocephlus in response to Vibrio parahaemolyticus infection

The black seabream (Sparus macrocephlus) is an economically pivotal aquaculture species cultured in China and Southeast Asian countries. To understand the molecular immune mechanisms underlying the response to Vibrio parahaemolyticus, a comparative gene transcription analysis were performed with utilized fresh livers of V. parahaemolyticus-immunized Sparus macrocephlus with a control group through RNA-Seq technology. A total of 256663 contigs were obtained after excluded the low-quality sequences and assembly. The average length of contigs collected from this research is 1066.93 bp. Furthermore, blast analysis indicates 30747 contigs were annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was employed to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. A total of 14470 genes were discovered in 303 KEGG pathways. RSEM and EdgeR were introduced to estimate 3841 genes significantly different expressed (False Discovery Rate<0.001) which includes 4072 up-regulated genes and 3771 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes were conducted to reveal the major immune-related pathways which refer to the toll-like receptor, complement, coagulation cascades, and chemokine signaling pathways. In addition, 92175 potential simple sequence repeats (SSRs) and 121912 candidate single nucleotide polymorphisms (SNPs) were detected and identified sequencely in the Sparus macrocephlus liver transcriptome. This research characterized a gene expression pattern for normal and the V. parahaemolyticus -immunized Sparus macrocephlus for the first time and not only sheds new light on the molecular mechanisms underlying the host-V. parahaemolyticus interaction but contribute to facilitate future studies on Sparus macrocephlus gene expression and functional genomics.

Neural stem cell heterogeneity in the mammalian forebrain

The brain was long considered an organ that underwent very little change after development. It is now well established that the mammalian central nervous system contains neural stem cells that generate progeny that are capable of making new neurons, astrocytes, and oligodendrocytes throughout life. The field has advanced rapidly as it strives to understand the basic biology of these precursor cells, and explore their potential to promote brain repair. The purpose of this review is to present current knowledge about the diversity of neural stem cells in vitro and in vivo, and highlight distinctions between neural stem cell populations, throughout development, and within the niche. A comprehensive understanding of neural stem cell heterogeneity will provide insights into the cellular and molecular regulation of neural development and lifelong neurogenesis, and will guide the development of novel strategies to promote regeneration and neural repair.

Identification of Intrinsic Axon Growth Modulators for Intact CNS Neurons after Injury

Functional deficits persist after spinal cord injury (SCI) because axons in the adult mammalian central nervous system (CNS) fail to regenerate. However, modest levels of spontaneous functional recovery are typically observed after trauma and are thought to be mediated by the plasticity of intact circuitry. The mechanisms underlying intact circuit plasticity are not delineated. Here, we characterize the in vivo transcriptome of sprouting intact neurons from Ngr1 null mice after partial SCI. We identify the lysophosphatidic acid signaling modulators LPPR1 and LPAR1 as intrinsic axon growth modulators for intact corticospinal motor neurons after adjacent injury. Furthermore, in vivo LPAR1 inhibition or LPPR1 overexpression enhances sprouting of intact corticospinal tract axons and yields greater functional recovery after unilateral brainstem lesion in wild-type mice. Thus, the transcriptional profile of injury-induced sprouting of intact neurons reveals targets for therapeutic enhancement of axon growth initiation and new synapse formation.

Transcriptomic profile analysis of mouse neural tube development by RNA-Seq

The neural tube is the primordium of the central nervous system (CNS) in which its development is not entirely clear. Understanding the cellular and molecular basis of neural tube development could, therefore, provide vital clues to the mechanism of neural tube defects (NTDs). Here, we investigated the gene expression profiles of three different time points (embryonic day (E) 8.5, 9.5 and 10.5) of mouse neural tube by using RNA-seq approach. About 391 differentially expressed genes (DEGs) were screened during mouse neural tube development, including 45 DEGs involved in CNS development, among which Bmp2, Ascl1, Olig2, Lhx1, Wnt7b and Eomes might play the important roles. Of 45 DEGs, Foxp2, Eomes, Hoxb3, Gpr56, Hap1, Nkx2-1, Sez6l2, Wnt7b, Tbx20, Nfib, Cntn1 and Dcx had different isoforms, and the opposite expression pattern of different isoforms was observed for Gpr56, Nkx2-1 and Sez6l2. In addition, alternative splicing, such as mutually exclusive exon, retained intron, skipped exon and alternative 3' splice site was identified in 10 neural related differentially splicing genes, including Ngrn, Ddr1, Dctn1, Dnmt3b, Ect2, Map2, Mbnl1, Meis2, Vcan and App. Moreover, seven neural splicing factors, such as Nova1/2, nSR100/Srrm4, Elavl3/4, Celf3 and Rbfox1 were differentially expressed during mouse neural tube development. Interestingly, nine DEGs identified above were dysregulated in retinoic acid-induced NTDs model, indicating the possible important role of these genes in NTDs. Taken together, our study provides more comprehensive information on mouse neural tube development, which might provide new insights on NTDs occurrence. © 2017 IUBMB Life, 69(9):706-719, 2017.

Sex differences in the molecular signature of the developing mouse hippocampus

BACKGROUND

A variety of neurological disorders, including Alzheimer's disease, Parkinson's disease, major depressive disorder, dyslexia and autism, are differentially prevalent between females and males. To better understand the possible molecular basis for the sex-biased nature of neurological disorders, we used a developmental series of female and male mice at 1, 2, and 4 months of age to assess both mRNA and protein in the hippocampus with RNA-sequencing and mass-spectrometry, respectively.

RESULTS

The transcriptomic analysis identifies 2699 genes that are differentially expressed between animals of different ages. The bulk of these differentially expressed genes are changed in both sexes at one or more ages, but a total of 198 transcripts are differentially expressed between females and males at one or more ages. The number of transcripts that are differentially expressed between females and males is greater in adult animals than in younger animals. Additionally, we identify 69 transcripts that show complex and sex-specific patterns of temporal regulation through postnatal development, 8 of which are heat-shock proteins. We also find a modest correlation between levels of mRNA and protein in the mouse hippocampus (Rho = 0.53).

CONCLUSION

This study adds to the substantial body of evidence for transcriptomic regulation in the hippocampus during postnatal development. Additionally, this analysis reveals sex differences in the transcriptome of the developing mouse hippocampus, and further clarifies the need to include both female and male mice in longitudinal studies involving molecular changes in the hippocampus.

RNA sequencing analysis to capture the transcriptome landscape during skin ulceration syndrome progression in sea cucumber Apostichopus japonicus

Background

Sea cucumber Apostichopus japonicus is an important economic species in China, which is affected by various diseases; skin ulceration syndrome (SUS) is the most serious. In this study, we characterized the transcriptomes in A. japonicus challenged with Vibrio splendidus to elucidate the changes in gene expression throughout the three stages of SUS progression.

Results

RNA sequencing of 21 cDNA libraries from various tissues and developmental stages of SUS-affected A. japonicus yielded 553 million raw reads, of which 542 million high-quality reads were generated by deep-sequencing using the Illumina HiSeq™ 2000 platform. The reference transcriptome comprised a combination of the Illumina reads, 454 sequencing data and Sanger sequences obtained from the public database to generate 93,163 unigenes (average length, 1,052 bp; N50 = 1,575 bp); 33,860 were annotated. Transcriptome comparisons between healthy and SUS-affected A. japonicus revealed greater differences in gene expression profiles in the body walls (BW) than in the intestines (Int), respiratory trees (RT) and coelomocytes (C). Clustering of expression models revealed stable up-regulation as the main pattern occurring in the BW throughout the three stages of SUS progression. Significantly affected pathways were associated with signal transduction, immune system, cellular processes, development and metabolism. Ninety-two differentially expressed genes (DEGs) were divided into four functional categories: attachment/pathogen recognition (17), inflammatory reactions (38), oxidative stress response (7) and apoptosis (30). Using quantitative real-time PCR, twenty representative DEGs were selected to validate the sequencing results. The Pearson’s correlation coefficient (R) of the 20 DEGs ranged from 0.811 to 0.999, which confirmed the consistency and accuracy between these two approaches.

Conclusions

Dynamic changes in global gene expression occur during SUS progression in A. japonicus. Elucidation of these changes is important in clarifying the molecular mechanisms associated with the development of SUS in sea cucumber.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2810-3) contains supplementary material, which is available to authorized users.

Liver Transcriptome Analysis of the Large Yellow Croaker (Larimichthys crocea) during Fasting by Using RNA-Seq

The large yellow croaker (Larimichthys crocea) is an economically important fish species in Chinese mariculture industry. To understand the molecular basis underlying the response to fasting, Illumina HiSeq^TM 2000 was used to analyze the liver transcriptome of fasting large yellow croakers. A total of 54,933,550 clean reads were obtained and assembled into 110,364 contigs. Annotation to the NCBI database identified a total of 38,728 unigenes, of which 19,654 were classified into Gene Ontology and 22,683 were found in Kyoto Encyclopedia of Genes and Genomes (KEGG). Comparative analysis of the expression profiles between fasting fish and normal-feeding fish identified a total of 7,623 differentially expressed genes (P < 0.05), including 2,500 upregulated genes and 5,123 downregulated genes. Dramatic differences were observed in the genes involved in metabolic pathways such as fat digestion and absorption, citrate cycle, and glycolysis/gluconeogenesis, and the similar results were also found in the transcriptome of skeletal muscle. Further qPCR analysis confirmed that the genes encoding the factors involved in those pathways significantly changed in terms of expression levels. The results of the present study provide insights into the molecular mechanisms underlying the metabolic response of the large yellow croaker to fasting as well as identified areas that require further investigation.

Expression and localization of Inter-alpha Inhibitors in rodent brain

Inter-alpha Inhibitor Proteins (IAIPs) are a family of related serine protease inhibitors. IAIPs are important components of the systemic innate immune system. We have identified endogenous IAIPs in the central nervous system (CNS) of sheep during development and shown that treatment with IAIPs reduces neuronal cell death and improves behavioral outcomes in neonatal rats after hypoxic-ischemic brain injury. The presence of IAIPs in CNS along with their exogenous neuroprotective properties suggests that endogenous IAIPs could be part of the innate immune system in CNS. The purpose of this study was to characterize expression and localization of IAIPs in CNS. We examined cellular expressions of IAIPs in vitro in cultured cortical mouse neurons, in cultured rat neurons, microglia, and astrocytes, and in vivo on brain sections by immunohistochemistry from embryonic (E) day 18 mice and postnatal (P) day 10 rats. Cultured cortical mouse neurons expressed the light chain gene Ambp and heavy chain genes Itih-1, 2, 3, 4, and 5 mRNA transcripts and IAIP proteins. IAIP proteins were detected by immunohistochemistry in cultured cells as well as brain sections from E18 mice and P10 rats. Immunoreactivity was found in neurons, microglia, astrocytes and oligodendroglia in multiple brain regions including cortex and hippocampus, as well as within both the ependyma and choroid plexus. Our findings suggest that IAIPs are endogenous proteins expressed in a wide variety of cell types and regions both in vitro and in vivo in rodent CNS. We speculate that endogenous IAIPs may represent endogenous neuroprotective immunomodulatory proteins within the CNS.

Transcriptome analysis of the Larimichthys crocea liver in response to Cryptocaryon irritans

The large yellow croaker (Larimichthys crocea) is an economically important marine fish cultured in China and East Asian countries and is facing a serious threat from Cryptocaryon irritans, which is a protozoan ectoparasite that infects most reared marine fish species. To understand the molecular immune mechanisms underlying the response to C. irritans, we first performed a comparative gene transcription analysis using livers from C. irritans-immunized L. croceas and from a control group through RNA-Seq technology. After the removal of low-quality sequences and assembly, 51360 contigs were obtained, with an average length of 1066.93 bp. Further, a blast analysis indicates that 30747 contigs can be annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was used to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. Moreover, 14470 genes were found in 303 KEGG pathways. We used RSEM and EdgeR to determine that 3841 genes were significantly differentially expressed (FDR < 0.001), including 2129 up-regulated genes and 1712 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes revealed major immune-related pathways, including the toll-like receptor, complement and coagulation cascades, and chemokine signaling pathways. In addition, 28748 potential simple sequence repeats (SSRs) were detected from 12776 transcripts, and 62992 candidate single nucleotide polymorphisms (SNPs) were identified in the L. croceas liver transcriptome. This study characterized a gene expression pattern for normal and C. irritans-immunized L. croceas for the first time and not only sheds new light on the molecular mechanisms underlying the host-C. irritans interaction but also facilitates future studies on L. croceas gene expression and functional genomics.

Suspension arrays based on nanoparticle-encoded microspheres for high-throughput multiplexed detection

Spectrometrically or optically encoded microsphere based suspension array technology (SAT) is applicable to the high-throughput, simultaneous detection of multiple analytes within a small, single sample volume. Thanks to the rapid development of nanotechnology, tremendous progress has been made in the multiplexed detecting capability, sensitivity, and photostability of suspension arrays. In this review, we first focus on the current stock of nanoparticle-based barcodes as well as the manufacturing technologies required for their production. We then move on to discuss all existing barcode-based bioanalysis patterns, including the various labels used in suspension arrays, label-free platforms, signal amplification methods, and fluorescence resonance energy transfer (FRET)-based platforms. We then introduce automatic platforms for suspension arrays that use superparamagnetic nanoparticle-based microspheres. Finally, we summarize the current challenges and their proposed solutions, which are centered on improving encoding capacities, alternative probe possibilities, nonspecificity suppression, directional immobilization, and "point of care" platforms. Throughout this review, we aim to provide a comprehensive guide for the design of suspension arrays, with the goal of improving their performance in areas such as multiplexing capacity, throughput, sensitivity, and cost effectiveness. We hope that our summary on the state-of-the-art development of these arrays, our commentary on future challenges, and some proposed avenues for further advances will help drive the development of suspension array technology and its related fields.

De novo assembly and characterization of the spleen transcriptome of common carp (Cyprinus carpio) using Illumina paired-end sequencing

Common carp (Cyprinus carpio) is one of the most important aquacultured species of the family Cyprinidae, and breeding this species for disease resistance is becoming more and more important. However, at the genome or transcriptome levels, study of the immunogenetics of disease resistance in the common carp is lacking. In this study, 60,316,906 and 75,200,328 paired-end clean reads were obtained from two cDNA libraries of the common carp spleen by Illumina paired-end sequencing technology. Totally, 130,293 unique transcript fragments (unigenes) were assembled, with an average length of 1400.57 bp. Approximately 105,612 (81.06%) unigenes could be annotated according to their homology with matches in the Nr, Nt, Swiss-Prot, COG, GO, or KEGG databases, and they were found to represent 46,747 non-redundant genes. Comparative analysis showed that 59.82% of the unigenes have significant similarity to zebrafish Refseq proteins. Gene expression comparison revealed that 10,432 and 6889 annotated unigenes were, respectively, up- and down-regulated with at least twofold changes between two developmental stages of the common carp spleen. Gene ontology and KEGG analysis were performed to classify all unigenes into functional categories for understanding gene functions and regulation pathways. In addition, 46,847 simple sequence repeats (SSRs) were detected from 35,618 unigenes, and a large number of single nucleotide polymorphism (SNP) and insertion/deletion (INDEL) sites were identified in the spleen transcriptome of common carp. This study has characterized the spleen transcriptome of the common carp for the first time, providing a valuable resource for a better understanding of the common carp immune system and defense mechanisms. This knowledge will also facilitate future functional studies on common carp immunogenetics that may eventually be applied in breeding programs.

Solexa-Sequencing Based Transcriptome Study of Plaice Skin Phenotype in Rex Rabbits (Oryctolagus cuniculus)

BACKGROUND

Fur is an important genetically-determined characteristic of domestic rabbits; rabbit furs are of great economic value. We used the Solexa sequencing technology to assess gene expression in skin tissues from full-sib Rex rabbits of different phenotypes in order to explore the molecular mechanisms associated with fur determination.

METHODOLOGY/PRINCIPAL FINDINGS

Transcriptome analysis included de novo assembly, gene function identification, and gene function classification and enrichment. We obtained 74,032,912 and 71,126,891 short reads of 100 nt, which were assembled into 377,618 unique sequences by Trinity strategy (N50=680 nt). Based on BLAST results with known proteins, 50,228 sequences were identified at a cut-off E-value ≥ 10-5. Using Blast to Gene Ontology (GO), Clusters of Orthologous Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG), we obtained several genes with important protein functions. A total of 308 differentially expressed genes were obtained by transcriptome analysis of plaice and un-plaice phenotype animals; 209 additional differentially expressed genes were not found in any database. These genes included 49 that were only expressed in plaice skin rabbits. The novel genes may play important roles during skin growth and development. In addition, 99 known differentially expressed genes were assigned to PI3K-Akt signaling, focal adhesion, and ECM-receptor interactin, among others. Growth factors play a role in skin growth and development by regulating these signaling pathways. We confirmed the altered expression levels of seven target genes by qRT-PCR. And chosen a key gene for SNP to found the differentially between plaice and un-plaice phenotypes rabbit.

CONCLUSIONS/SIGNIFICANCE

The rabbit transcriptome profiling data provide new insights in understanding the molecular mechanisms underlying rabbit skin growth and development.

Identification and expression patterns of novel long non-coding RNAs in neural progenitors of the developing mammalian cortex

Long non-coding (lnc)RNAs play key roles in many biological processes. Elucidating the function of lncRNAs in cell type specification during organ development requires knowledge about their expression in individual progenitor types rather than in whole tissues. To achieve this during cortical development, we used a dual-reporter mouse line to isolate coexisting proliferating neural stem cells, differentiating neurogenic progenitors and newborn neurons and assessed the expression of lncRNAs by paired-end, high-throughput sequencing. We identified 379 genomic loci encoding novel lncRNAs and performed a comprehensive assessment of cell-specific expression patterns for all, annotated and novel, lncRNAs described to date. Our study provides a powerful new resource for studying these elusive transcripts during stem cell commitment and neurogenesis.

Tox: a multifunctional transcription factor and novel regulator of mammalian corticogenesis

Major efforts are invested to characterize the factors controlling the proliferation of neural stem cells. During mammalian corticogenesis, our group has identified a small pool of genes that are transiently downregulated in the switch of neural stem cells to neurogenic division and reinduced in newborn neurons. Among these switch genes, we found Tox, a transcription factor with hitherto uncharacterized roles in the nervous system. Here, we investigated the role of Tox in corticogenesis by characterizing its expression at the tissue, cellular and temporal level. We found that Tox is regulated by calcineurin/Nfat signalling. Moreover, we combined DNA adenine methyltransferase identification (DamID) with deep sequencing to characterize the chromatin binding properties of Tox including its motif and downstream transcriptional targets including Sox2, Tbr2, Prox1 and other key factors. Finally, we manipulated Tox in the developing brain and validated its multiple roles in promoting neural stem cell proliferation and neurite outgrowth of newborn neurons. Our data provide a valuable resource to study the role of Tox in other tissues and highlight a novel key player in brain development.

Transcriptome and expression profiling analysis of the hemocytes reveals a large number of immune-related genes in mud crab Scylla paramamosain during Vibrio parahaemolyticus infection

BACKGROUND

Mud crab Scylla paramamosain is an economically important marine species in China. However, frequent outbreaks of infectious diseases caused by marine bacteria, such as Vibrio parahaemolyticus, result in great economic losses.

METHODOLOGY/PRINCIPAL FINDINGS

Comparative de novo transcriptome analysis of S. paramamosain infected with V. parahaemolyticus was carried out to investigate the molecular mechanisms underlying the immune response to pathogenic bacteria by using the Illumina paired-end sequencing platform. A total of 52,934,042 clean reads from the hemocytes of V. parahaemolyticus-infected mud crabs and controls were obtained and assembled into 186,193 contigs. 59,120 unigenes were identified from 81,709 consensus sequences of mud crabs and 48,934 unigenes were matched proteins in the Nr or Swissprot databases. Among these, 10,566 unigenes belong to 3 categories of Gene Ontology, 25,349 to 30 categories of KEGG, and 15,191 to 25 categories of COG database, covering almost all functional categories. By using the Solexa/Illumina's DGE platform, 1213 differentially expressed genes (P<0.05), including 538 significantly up-regulated and 675 down-regulated, were detected in V. parahaemolyticus-infected crabs as compared to that in the controls. Transcript levels of randomly-chosen genes were further measured by quantitative real-time PCR to confirm the expression profiles. Many differentially expressed genes are involved in various immune processes, including stimulation of the Toll pathway, Immune Deficiency (IMD) pathway, Ras-regulated endocytosis, and proPO-activating system.

CONCLUSIONS/SIGNIFICANCE

Analysis of the expression profile of crabs under infection provides invaluable new data for biological research in S. paramamosain, such as the identification of novel genes in the hemocytes during V. parahaemolyticus infection. These results will facilitate our comprehensive understanding of the mechanisms involved in the immune response to bacterial infection and will be helpful for diseases prevention in crab aquaculture.

De novo characterization of the spleen transcriptome of the large yellow croaker (Pseudosciaena crocea) and analysis of the immune relevant genes and pathways involved in the antiviral response

The large yellow croaker (Pseudosciaena crocea) is an economically important marine fish in China. To understand the molecular basis for antiviral defense in this species, we used Illumia paired-end sequencing to characterize the spleen transcriptome of polyriboinosinic:polyribocytidylic acid [poly(I:C)]-induced large yellow croakers. The library produced 56,355,728 reads and assembled into 108,237 contigs. As a result, 15,192 unigenes were found from this transcriptome. Gene ontology analysis showed that 4,759 genes were involved in three major functional categories: biological process, cellular component, and molecular function. We further ascertained that numerous consensus sequences were homologous to known immune-relevant genes. Kyoto Encyclopedia of Genes and Genomes orthology mapping annotated 5,389 unigenes and identified numerous immune-relevant pathways. These immune-relevant genes and pathways revealed major antiviral immunity effectors, including but not limited to: pattern recognition receptors, adaptors and signal transducers, the interferons and interferon-stimulated genes, inflammatory cytokines and receptors, complement components, and B-cell and T-cell antigen activation molecules. Moreover, the partial genes of Toll-like receptor signaling pathway, RIG-I-like receptors signaling pathway, Janus kinase-Signal Transducer and Activator of Transcription (JAK-STAT) signaling pathway, and T-cell receptor (TCR) signaling pathway were found to be changed after poly(I:C) induction by real-time polymerase chain reaction (PCR) analysis, suggesting that these signaling pathways may be regulated by poly(I:C), a viral mimic. Overall, the antivirus-related genes and signaling pathways that were identified in response to poly(I:C) challenge provide valuable leads for further investigation of the antiviral defense mechanism in the large yellow croaker.

voom: Precision weights unlock linear model analysis tools for RNA-seq read counts

New normal linear modeling strategies are presented for analyzing read counts from RNA-seq experiments. The voom method estimates the mean-variance relationship of the log-counts, generates a precision weight for each observation and enters these into the limma empirical Bayes analysis pipeline. This opens access for RNA-seq analysts to a large body of methodology developed for microarrays. Simulation studies show that voom performs as well or better than count-based RNA-seq methods even when the data are generated according to the assumptions of the earlier methods. Two case studies illustrate the use of linear modeling and gene set testing methods.

Transcriptomic changes in brain development

The transcriptome changes hugely during development of the brain. Whole genes, alternate exons, and single base pair changes related to RNA editing all show differences between embryonic and mature brain. Collectively, these changes control proteomic diversity as the brain develops. Additionally, there are many changes in noncoding RNAs (miRNA and lncRNA) that interact with mRNA to influence the overall transcriptional landscape. Here, we will discuss what is known about such changes in brain development, particularly focusing on high-throughput approaches and how those can be used to infer mechanisms by which gene expression is controlled in the brain as it matures.

Transcriptome sequencing during mouse brain development identifies long non-coding RNAs functionally involved in neurogenic commitment

Transcriptome analysis of somatic stem cells and their progeny is fundamental to identify new factors controlling proliferation versus differentiation during tissue formation. Here, we generated a combinatorial, fluorescent reporter mouse line to isolate proliferating neural stem cells, differentiating progenitors and newborn neurons that coexist as intermingled cell populations during brain development. Transcriptome sequencing revealed numerous novel long non-coding (lnc)RNAs and uncharacterized protein-coding transcripts identifying the signature of neurogenic commitment. Importantly, most lncRNAs overlapped neurogenic genes and shared with them a nearly identical expression pattern suggesting that lncRNAs control corticogenesis by tuning the expression of nearby cell fate determinants. We assessed the power of our approach by manipulating lncRNAs and protein-coding transcripts with no function in corticogenesis reported to date. This led to several evident phenotypes in neurogenic commitment and neuronal survival, indicating that our study provides a remarkably high number of uncharacterized transcripts with hitherto unsuspected roles in brain development. Finally, we focussed on one lncRNA, Miat, whose manipulation was found to trigger pleiotropic effects on brain development and aberrant splicing of Wnt7b. Hence, our study suggests that lncRNA-mediated alternative splicing of cell fate determinants controls stem-cell commitment during neurogenesis.

DNA methylation markers in the postnatal developing rat brain

In spite of intense interest in how altered epigenetic processes including DNA methylation may contribute to psychiatric and neurodevelopmental disorders, there is a limited understanding of how methylation processes change during early postnatal brain development. The present study used in situ hybridization to assess mRNA expression for the three major DNA methyltranserases (DNMTs)--DNMT1, DNMT3a and DNMT3b--in the developing rat brain at seven developmental timepoints: postnatal days (P) 1, 4, 7, 10, 14, 21, and 75. We also assessed 5-methylcytosine levels (an indicator of global DNA methylation) in selected brain regions during the first three postnatal weeks. DNMT1, DNMT3a and DNMT3b mRNAs are widely expressed throughout the adult and postnatal developing rat brain. Overall, DNMT mRNA levels reached their highest point in the first week of life and gradually decreased over the first three postnatal weeks within the hippocampus, amygdala, striatum, cingulate and lateral septum. Global DNA methylation levels did not follow this developmental pattern; methylation levels gradually increased over the first three postnatal weeks in the hippocampus, and remained stable in the developing amygdala and prefrontal cortex. Our results contribute to a growing understanding of how DNA methylation markers unfold in the developing brain, and highlight how these developmental processes may differ within distinct brain regions.

Long non-coding RNA identification over mouse brain development by integrative modeling of chromatin and genomic features

In silico prediction of genomic long non-coding RNAs (lncRNAs) is prerequisite to the construction and elucidation of non-coding regulatory network. Chromatin modifications marked by chromatin regulators are important epigenetic features, which can be captured by prevailing high-throughput approaches such as ChIP sequencing. We demonstrate that the accuracy of lncRNA predictions can be greatly improved when incorporating high-throughput chromatin modifications over mouse embryonic stem differentiation toward adult Cerebellum by logistic regression with LASSO regularization. The discriminating features include H3K9me3, H3K27ac, H3K4me1, open reading frames and several repeat elements. Importantly, chromatin information is suggested to be complementary to genomic sequence information, highlighting the importance of an integrated model. Applying integrated model, we obtain a list of putative lncRNAs based on uncharacterized fragments from transcriptome assembly. We demonstrate that the putative lncRNAs have regulatory roles in vicinity of known gene loci by expression and Gene Ontology enrichment analysis. We also show that the lncRNA expression specificity can be efficiently modeled by the chromatin data with same developmental stage. The study not only supports the biological hypothesis that chromatin can regulate expression of tissue-specific or developmental stage-specific lncRNAs but also reveals the discriminating features between lncRNA and coding genes, which would guide further lncRNA identifications and characterizations.

Identification and characterization of long non-coding RNAs related to mouse embryonic brain development from available transcriptomic data

Long non-coding RNAs (lncRNAs) as a key group of non-coding RNAs have gained widely attention. Though lncRNAs have been functionally annotated and systematic explored in higher mammals, few are under systematical identification and annotation. Owing to the expression specificity, known lncRNAs expressed in embryonic brain tissues remain still limited. Considering a large number of lncRNAs are only transcribed in brain tissues, studies of lncRNAs in developmental brain are therefore of special interest. Here, publicly available RNA-sequencing (RNA-seq) data in embryonic brain are integrated to identify thousands of embryonic brain lncRNAs by a customized pipeline. A significant proportion of novel transcripts have not been annotated by available genomic resources. The putative embryonic brain lncRNAs are shorter in length, less spliced and show less conservation than known genes. The expression of putative lncRNAs is in one tenth on average of known coding genes, while comparable with known lncRNAs. From chromatin data, putative embryonic brain lncRNAs are associated with active chromatin marks, comparable with known lncRNAs. Embryonic brain expressed lncRNAs are also indicated to have expression though not evident in adult brain. Gene Ontology analysis of putative embryonic brain lncRNAs suggests that they are associated with brain development. The putative lncRNAs are shown to be related to possible cis-regulatory roles in imprinting even themselves are deemed to be imprinted lncRNAs. Re-analysis of one knockdown data suggests that four regulators are associated with lncRNAs. Taken together, the identification and systematic analysis of putative lncRNAs would provide novel insights into uncharacterized mouse non-coding regions and the relationships with mammalian embryonic brain development.

Skin transcriptome profiles associated with coat color in sheep

Background

Previous molecular genetic studies of physiology and pigmentation of sheep skin have focused primarily on a limited number of genes and proteins. To identify additional genes that may play important roles in coat color regulation, Illumina sequencing technology was used to catalog global gene expression profiles in skin of sheep with white versus black coat color.

Results

There were 90,006 and 74,533 unigenes assembled from the reads obtained from white and black sheep skin, respectively. Genes encoding for the ribosomal proteins and keratin associated proteins were most highly expressed. A total of 2,235 known genes were differentially expressed in black versus white sheep skin, with 479 genes up-regulated and 1,756 genes down-regulated. A total of 845 novel genes were differentially expressed in black versus white sheep skin, consisting of 107 genes which were up-regulated (including 2 highly expressed genes exclusively expressed in black sheep skin) and 738 genes that were down-regulated. There was also a total of 49 known coat color genes expressed in sheep skin, from which 13 genes showed higher expression in black sheep skin. Many of these up-regulated genes, such as DCT, MATP, TYR and TYRP1, are members of the components of melanosomes and their precursor ontology category.

Conclusion

The white and black sheep skin transcriptome profiles obtained provide a valuable resource for future research to understand the network of gene expression controlling skin physiology and melanogenesis in sheep.

Post-transcriptional regulatory elements and spatiotemporal specification of neocortical stem cells and projection neurons

The mature neocortex is a unique six-layered mammalian brain region. It is composed of morphologically and functionally distinct subpopulations of primary projection neurons that form complex circuits across the central nervous system. The precisely-timed generation of projection neurons from neural stem cells governs their differentiation, postmitotic specification, and signaling, and is critical for cognitive and sensorimotor ability. Developmental perturbations to the birthdate, location, and connectivity of neocortical neurons are observed in neurological and psychiatric disorders. These facts are highlighting the importance of the precise spatiotemporal development of the neocortex regulated by intricate transcriptional, but also complex post-transcriptional events. Indeed, mRNA transcripts undergo many post-transcriptional regulatory steps before the production of functional proteins, which specify neocortical neural stem cells and subpopulations of neocortical neurons. Therefore, particular attention is paid to the differential post-transcriptional regulation of key transcripts by RNA-binding proteins, including splicing, localization, stability, and translation. We also present a transcriptome screen of candidate molecules associated with post-transcriptional mRNA processing that are differentially expressed at key developmental time points across neocortical prenatal neurogenesis.

Deep mRNA sequencing analysis to capture the transcriptome landscape of zebrafish embryos and larvae

Transcriptome analysis is a powerful tool to obtain large amount genome-scale gene expression profiles. Despite its extensive usage to diverse biological problems in the last decade, transcriptomic researches approaching the zebrafish embryonic development have been very limited. Several recent studies have made great progress in this direction, yet the large gap still exists, especially regarding to the transcriptome dynamics of embryonic stages from early gastrulation onwards. Here, we present a comprehensive analysis about the transcriptomes of 9 different stages covering 7 major periods (cleavage, blastula, gastrula, segmentation, pharyngula, hatching and early larval stage) in zebrafish development, by recruiting the RNA-sequencing technology. We detected the expression for at least 24,065 genes in at least one of the 9 stages. We identified 16,130 genes that were significantly differentially expressed between stages and were subsequently classified into six clusters. Each revealed gene cluster had distinct expression patterns and characteristic functional pathways, providing a framework for the understanding of the developmental transcriptome dynamics. Over 4000 genes were identified as preferentially expressed in one of the stages, which could be of high relevance to stage-specific developmental and molecular events. Among the 68 transcription factor families active during development, most had enhanced average expression levels and thus might be crucial for embryogenesis, whereas the inactivation of the other families was likely required by the activation of the zygotic genome. We discussed our RNA-seq data together with previous findings about the Wnt signaling pathway and some other genes with known functions, to show how our data could be used to advance our understanding about these developmental functional elements. Our study provides ample information for further study about the molecular and cellular mechanisms underlying vertebrate development.

Transcriptome profiling of Chironomus kiinensis under phenol stress using Solexa sequencing technology

Phenol is a major pollutant in aquatic ecosystems due to its chemical stability, water solubility and environmental mobility. To date, little is known about the molecular modifications of invertebrates under phenol stress. In the present study, we used Solexa sequencing technology to investigate the transcriptome and differentially expressed genes (DEGs) of midges (Chironomus kiinensis) in response to phenol stress. A total of 51,518,972 and 51,150,832 clean reads in the phenol-treated and control libraries, respectively, were obtained and assembled into 51,014 non-redundant (Nr) consensus sequences. A total of 6,032 unigenes were classified by Gene Ontology (GO), and 18,366 unigenes were categorized into 238 Kyoto Encyclopedia of Genes and Genomes (KEGG) categories. These genes included representatives from almost all functional categories. A total of 10,724 differentially expressed genes (P value <0.05) were detected in a comparative analysis of the expression profiles between phenol-treated and control C. kiinensis including 8,390 upregulated and 2,334 downregulated genes. The expression levels of 20 differentially expressed genes were confirmed by real-time RT-PCR, and the trends in gene expression that were observed matched the Solexa expression profiles, although the magnitude of the variations was different. Through pathway enrichment analysis, significantly enriched pathways were identified for the DEGs, including metabolic pathways, aryl hydrocarbon receptor (AhR), pancreatic secretion and neuroactive ligand-receptor interaction pathways, which may be associated with the phenol responses of C. kiinensis. Using Solexa sequencing technology, we identified several groups of key candidate genes as well as important biological pathways involved in the molecular modifications of chironomids under phenol stress.

Transcriptome map of mouse isochores in embryonic and neonatal cortex

Several studies on adult tissues agree on the presence of a positive effect of the genomic and genic base composition on mammalian gene expression. Recent literature supports the idea that during developmental processes GC-poor genomic regions are preferentially implicated. We investigate the relationship between the compositional properties of the isochores and of the genes with their respective expression activity during developmental processes. Using RNA-seq data from two distinct developmental stages of the mouse cortex, embryonic day 18 (E18) and postnatal day 7 (P7), we established for the first time a developmental-related transcriptome map of the mouse isochores. Additionally, for each stage we estimated the correlation between isochores' GC level and their expression activity, and the genes' expression patterns for each isochore family. Our analyses add evidence supporting the idea that during development GC-poor isochores are preferentially implicated, and confirm the positive effect of genes' GC level on their expression activity.

Sox11 modulates neocortical development by regulating the proliferation and neuronal differentiation of cortical intermediate precursors

Neural precursor cells play important roles in the neocortical development, but the mechanisms of neural progenitor proliferation, neuronal differentiation, and migration, as well as patterning are still unclear. Sox11, one of SoxC family members, has been reported to be essential for embryonic and adult neurogenesis. But there is no report about the roles of Sox11 in corticogenesis. In order to investigate Sox11 function during cortical development, loss of function experiment was performed in this study. Knockdown of Sox11 by Sox11 siRNA constructs resulted in a diminished neuronal differentiation, but enhanced proliferation of intermediate progenitors. Accompanied with the high expression of Sox11 in the postmitotic neurons, but low expression of Sox11 in the dividing neural progenitors, all the observations indicate that Sox11 induces neuronal differentiation during the neocortical development.

RNA-seq liver transcriptome analysis reveals an activated MHC-I pathway and an inhibited MHC-II pathway at the early stage of vaccine immunization in zebrafish

Background

Zebrafish (Danio rerio) is a prominent vertebrate model of human development and pathogenic disease and has recently been utilized to study teleost immune responses to infectious agents threatening the aquaculture industry. In this work, to clarify the host immune mechanisms underlying the protective effects of a putative vaccine and improve its immunogenicity in the future efforts, high-throughput RNA sequencing technology was used to investigate the immunization-related gene expression patterns of zebrafish immunized with Edwardsiella tarda live attenuated vaccine.

Results

Average reads of 18.13 million and 14.27 million were obtained from livers of zebrafish immunized with phosphate buffered saline (mock) and E. tarda vaccine (WED), respectively. The reads were annotated with the Ensembl zebrafish database before differential expressed genes sequencing (DESeq) comparative analysis, which identified 4565 significantly differentially expressed genes (2186 up-regulated and 2379 down-regulated in WED; p<0.05). Among those, functional classifications were found in the Gene Ontology database for 3891 and in the Kyoto Encyclopedia of Genes and Genomes database for 3467. Several pathways involved in acute phase response, complement activation, immune/defense response, and antigen processing and presentation were remarkably affected at the early stage of WED immunization. Further qPCR analysis confirmed that the genes encoding the factors involved in major histocompatibility complex (MHC)-I processing pathway were up-regulated, while those involved in MHC-II pathway were down-regulated.

Conclusion

These data provided insights into the molecular mechanisms underlying zebrafish immune response to WED immunization and might aid future studies to develop a highly immunogenic vaccine against gram-negative bacteria in teleosts.

Transcriptome analysis of microRNAs in developing cerebral cortex of rat

BACKGROUND

The morphogenesis of the cerebral cortex depends on the precise control of gene expression during development. Small non-coding RNAs, including microRNAs and other groups of small RNAs, play profound roles in various physiological and pathological processes via their regulation of gene expression. A systematic analysis of the expression profile of small non-coding RNAs in developing cortical tissues is important for clarifying the gene regulation networks mediating key developmental events during cortical morphogenesis.

RESULTS

Global profiling of the small RNA transcriptome was carried out in rat cerebral cortex from E10 till P28 using next-generation sequencing technique. We found an extraordinary degree of developmental stage-specific expression of a large group of microRNAs. A group of novel microRNAs with functional hints were identified, and brain-enriched expression and Dicer-dependent production of high-abundant novel microRNAs were validated. Profound editing of known microRNAs at "seed" sequence and flanking sequence was observed, with much higher editing events detected at late postnatal stages than embryonic stages, suggesting the necessity of microRNA editing for the fine tuning of gene expression during the formation of complicated synaptic connections at postnatal stages.

CONCLUSION

Our analysis reveals extensive regulation of microRNAs during cortical development. The dataset described here will be a valuable resource for clarifying new regulatory mechanisms for cortical development and diseases and will greatly contribute to our understanding of the divergence, modification, and function of microRNAs.

A comparative, developmental, and clinical perspective of neurobehavioral sexual dimorphisms

Women and men differ in a wide variety of behavioral traits and in their vulnerability to developing certain mental disorders. This review endeavors to explore how recent preclinical and clinical research findings have enhanced our understanding of the factors that underlie these disparities. We start with a brief overview of some of the important genetic, molecular, and hormonal determinants that contribute to the process of sexual differentiation. We then discuss the importance of animal models in studying the mechanisms responsible for sex differences in neuropsychiatric disorders (e.g., drug dependence) - with a special emphasis on experimental models based on the neurodevelopmental and "three hits" hypotheses. Next, we describe the most common brain phenotypes observed in vivo with magnetic resonance imaging. We discuss the challenges in interpreting these phenotypes vis-à-vis the underlying neurobiology and revisit the known sex differences in brain structure from birth, through adolescence, and into adulthood. This is followed by a presentation of pertinent clinical and epidemiological data that point to important sex differences in the prevalence, course, and expression of psychopathologies such as schizophrenia, and mood disorders including major depression and posttraumatic stress disorder. Recent evidence implies that mood disorders and psychosis share some common genetic predispositions and neurobiological bases. Therefore, modern research is emphasizing dimensional representation of mental disorders and conceptualization of schizophrenia and major depression as a continuum of cognitive deficits and neurobiological abnormalities. Herein, we examine available evidence on cerebral sexual dimorphism to verify if sex differences vary quantitatively and/or qualitatively along the psychoses-depression continuum. Finally, sex differences in the prevalence of posttraumatic disorder and drug abuse have been described, and we consider the genomic and molecular data supporting these differences.

Regulation of gene expression during early neuronal differentiation: evidence for patterns conserved across neuron populations and vertebrate classes

Analysis of transcription factor function during neurogenesis has provided a huge amount of data on the generation and specification of diverse neuron populations in the central and peripheral nervous systems of vertebrates. However, an understanding of the induction of key neuron functions including electrical information processing and synaptic transmission lags seriously behind. Whereas pan-neuronal markers such as neurofilaments, neuron-specific tubulin and RNA-binding proteins have often been included in developmental analysis, the molecular players underlying electrical activity and transmitter release have been neglected in studies addressing gene expression during neuronal induction. Here, I summarize the evidence for a distinct accumulation pattern of mRNAs for synaptic proteins, a pattern that is delayed compared with pan-neuronal gene expression during neurogenesis. The conservation of this pattern across diverse avian and mammalian neuron populations suggests a common mechanism for the regulation of various sets of neuronal genes during initial neuronal differentiation. The co-regulation of genes coding for synaptic proteins from embryonic to postnatal development indicates that the expression of the players required for synaptic transmission shares common regulatory features. For the ion channels involved in neuronal electrical activity, such as voltage-gated sodium channels, the situation is less clear because of the lack of comparative studies early during neurogenesis. Transcription factors have been characterized that regulate the expression of synaptic proteins in vitro and in vivo. They currently do not explain the co-regulation of these genes across different neuron populations. The neuron-restrictive silencing factor NRSF/REST targets a large gene set, but not all of the genes coding for pan-neuronal, synaptic and ion channel proteins. The discrepancy between NRSF/REST loss-of-function and silencer-to-activator-switch studies leaves the full functional implications of this factor open. Together with microRNAs, splicing regulators, chromatin remodellers and an increasing list of transcriptional regulators, the factor is embedded in feedback circuits with the potential to orchestrate neuronal differentiation. The precise regulation of the coordinated expression of proteins underlying key neuronal functions by these circuits during neuronal induction is a major emerging topic.