Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gene expression networks”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Rapid Cis–Trans Coevolution Driven by a Novel Gene Retroposed from a Eukaryotic Conserved CCR4–NOT Component in Drosophila

Young, or newly evolved, genes arise ubiquitously across the tree of life, and they can rapidly acquire novel functions that influence a diverse array of biological processes. Previous work identified a young regulatory duplicate gene in Drosophila, Zeus that unexpectedly diverged rapidly from its parent, Caf40, an extremely conserved component in the CCR4–NOT machinery in post-transcriptional and post-translational regulation of eukaryotic cells, and took on roles in the male reproductive system. This neofunctionalization was accompanied by differential binding of the Zeus protein to loci throughout the Drosophila melanogaster genome. However, the way in which new DNA-binding proteins acquire and coevolve with their targets in the genome is not understood. Here, by comparing Zeus ChIP-Seq data from D. melanogaster and D. simulans to the ancestral Caf40 binding events from D. yakuba, a species that diverged before the duplication event, we found a dynamic pattern in which Zeus binding rapidly coevolved with a previously unknown DNA motif, which we term Caf40 and Zeus-Associated Motif (CAZAM), under the influence of positive selection. Interestingly, while both copies of Zeus acquired targets at male-biased and testis-specific genes, D. melanogaster and D. simulans proteins have specialized binding on different chromosomes, a pattern echoed in the evolution of the associated motif. Using CRISPR-Cas9-mediated gene knockout of Zeus and RNA-Seq, we found that Zeus regulated the expression of 661 differentially expressed genes (DEGs). Our results suggest that the evolution of young regulatory genes can be coupled to substantial rewiring of the transcriptional networks into which they integrate, even over short evolutionary timescales. Our results thus uncover dynamic genome-wide evolutionary processes associated with new genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Metabolic engineering of low-pH-tolerant non-model yeast, Issatchenkia orientalis , for production of citramalate

Methyl methacrylate (MMA) is an important petrochemical with many applications. However, its manufacture has a large environmental footprint. Combined biological and chemical synthesis (semisynthesis) may be a promising alternative to reduce both cost and environmental impact, but strains that can produce the MMA precursor (citramalate) at low pH are required. A non-conventional yeast, Issatchenkia orientalis, may prove ideal, as it can survive extremely low pH. Here, we demonstrate the engineering of I. orientalis for citramalate pro- duction. Using sequence similarity network analysis and subsequent DNA synthesis, we selected a more active citramalate synthase gene (cimA) variant for expression in I. orientalis. We then adapted a piggyBac transposon system for I. orientalis that allowed us to simultaneously explore the effects of different cimA gene copy numbers and integration locations. A batch fermentation showed the genome-integrated-cimA strains produced 2.0 g/L citramalate in 48 h and a yield of up to 7% mol citramalate/mol consumed glucose. These results demonstrate the potential of I. orientalis as a chassis for citramalate production.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Data for Metabolic Engineering of Low-pH-Tolerant Non-Model Yeast, Issatchenkia orientalis , for Production of Citramalate

Methyl methacrylate (MMA) is an important petrochemical with many applications. However, its manufacture has a large environmental footprint. Combined biological and chemical synthesis (semisynthesis) may be a promising alternative to reduce both cost and environmental impact, but strains that can produce the MMA precursor (citramalate) at low pH are required. A non-conventional yeast, Issatchenkia orientalis , may prove ideal, as it can survive extremely low pH. Here, we demonstrate the engineering of I. orientalis for citramalate production. Using sequence similarity network analysis and subsequent DNA synthesis, we selected a more active citramalate synthase gene (cimA) variant for expression in I. orientalis . We then adapted a piggyBac transposon system for I. orientalis that allowed us to simultaneously explore the effects of different cimA gene copy numbers and integration locations. A batch fermentation showed the genome-integrated-cimA strains produced 2.0 g/L citramalate in 48 h and a yield of up to 7% mol citramalate/mol consumed glucose. These results demonstrate the potential of I. orientalis as a chassis for citramalate production.

Conversion↗
DOE OSTI2021

Repulsive expansion dynamics in colony growth and gene expression

Spatial expansion of a population of cells can arise from growth of microorganisms, plant cells, and mammalian cells. It underlies normal or dysfunctional tissue development, and it can be exploited as the foundation for programming spatial patterns. This expansion is often driven by continuous growth and division of cells within a colony, which in turn pushes the peripheral cells outward. This process generates a repulsion velocity field at each location within the colony. Here we show that this process can be approximated as coarse-grained repulsive-expansion kinetics. This framework enables accurate and efficient simulation of growth and gene expression dynamics in radially symmetric colonies with homogenous z-directional distribution. It is robust even if cells are not spherical and vary in size. The simplicity of the resulting mathematical framework also greatly facilitates generation of mechanistic insights.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A network-enabled pipeline for gene discovery and validation in non-model plant species

Identifying key regulators of important genes in non-model crop species is challenging due to limited multi-omics resources. To address this, we introduce the network-enabled gene discovery pipeline NEEDLE, a user-friendly tool that systematically generates coexpression gene network modules, measures gene connectivity, and establishes network hierarchy to pinpoint key transcriptional regulators from dynamic transcriptome datasets. After validating its accuracy with two independent datasets, we applied NEEDLE to identify transcription factors (TFs) regulating the expression of cellulose synthase-like F6 ( CSLF6 ), a crucial cell wall biosynthetic gene, in Brachypodium and sorghum. Our analyses uncover regulators of CSLF6 and also shed light on the evolutionary conservation or divergence of gene regulatory elements among grass species. These results highlight NEEDLE’s capability to provide biologically relevant TF predictions and demonstrate its value for non-model plant species with dynamic transcriptome datasets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes

Effective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms. We investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables. Snap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Creb5 controls its own expression and directly induces the joint interzone regulatory program

Prior studies have indicated that the transcription factor Creb5 is expressed in the joint interzone, which contains the progenitors for all synovial joint tissues in both mouse and human embryos. In the absence of Creb5 function, most synovial joint interzones fail to form and the cartilage templates in the long bones remain fused. This earlier work did not clarify whether Creb5 initiates a cascade of signaling molecules, such as growth and differentiation factor 5 (Gdf5) and Wnt-family members, that in turn induce the formation of the joint interzone, or instead directly activates the expression of joint interzone markers. In the present study, an integrative analysis of the transcriptome, chromatin accessibility, and Creb5-occupancy in joint progenitors revealed that Creb5 directly binds to both its own two promoters and to the regulatory regions of Gdf5 and Sfrp2, each of whose expression in the joint interzone is Creb5-dependent. Functional enhancer analysis indicated that Creb5 binding sites in either the two Creb5 promoters, or in Gdf5 and Sfrp2 regulatory elements are necessary for these sequences to drive transgene expression in the developing synovial joints. While Creb5 directly drives Gdf5 and Sfrp2 expression in the inner joint interzone, Creb5 activates Barx1 expression specifically in the outer joint interzone. Our findings indicate that Creb5 initiates a regulatory network that both promotes the formation of synovial joints, and subsequently activates distinct transcriptional targets in the inner versus the outer regions of the joint interzone, thus regionalizing gene expression in the developing joint.

Zhang, Cheng-Hai↗
DOE OSTI2020

Network Modeling of Complex Data Sets

We demonstrate a selection of network and machine learning techniques useful in the analysis of complex datasets, including 2-way similarity networks, Markov clustering, enrichment statistical networks, FCROS differential analysis, and random forests. We demonstrate each of these techniques on the Populus trichocarpa gene expression atlas.

Jones, Piet C.↗
DOE OSTI2022

Host genotype controls ecological change in the leaf fungal microbiome

Leaf fungal microbiomes can be fundamental drivers of host plant success, as they contain pathogens that devastate crop plants and taxa that enhance nutrient uptake, discourage herbivory, and antagonize pathogens. We measured leaf fungal diversity with amplicon sequencing across an entire growing season in a diversity panel of switchgrass ( Panicum virgatum ). We also sampled a replicated subset of genotypes across 3 additional sites to compare the importance of time, space, ecology, and genetics. We found a strong successional pattern in the microbiome shaped both by host genetics and environmental factors. Further, we used genome-wide association (GWA) mapping and RNA sequencing to show that 3 cysteine-rich receptor-like kinases (crRLKs) were linked to a genetic locus associated with microbiome structure. We confirmed GWAS results in an independent set of genotypes for both the internal transcribed spacer (ITS) and large subunit (LSU) ribosomal DNA markers. Fungal pathogens were central to microbial covariance networks, and genotypes susceptible to pathogens differed in their expression of the 3 crRLKs, suggesting that host immune genes are a principal means of controlling the entire leaf microbiome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Identification of Novel Genomic Associations and Gene Candidates for Grain Starch Content in Sorghum

Starch accumulated in the endosperm of cereal grains as reserve energy for germination serves as a staple in human and animal nutrition. Unraveling genetic control for starch metabolism is important for breeding grains with high starch content. In this study, we used a sorghum association panel with 389 individuals and 141,557 single nucleotide polymorphisms (SNPs) to fit linear mixed models (LMM) for identifying genomic regions and potential candidate genes associated with starch content. Three associated genomic regions, one in chromosome (chr) 1 and two novel associations in chr-8, were identified using combination of LMM and Bayesian sparse LMM. All significant SNPs were located within protein coding genes, with SNPs ∼ 52 Mb of chr-8 encoding a Casperian strip membrane protein (CASP)-like protein (Sobic.008G111500) and a heat shock protein (HSP) 90 (Sobic.008G111600) that were highly expressed in reproductive tissues including within the embryo and endosperm. The HSP90 is a potential hub gene with gene network of 75 high-confidence first interactors that is enriched for five biochemical pathways including protein processing. The first interactors of HSP90 also showed high transcript abundance in reproductive tissues. The candidates of this study are likely involved in intricate metabolic pathways and represent candidate gene targets for source-sink activities and drought and heat stress tolerance during grain filling.

genomics↗
DOE OSTI2021

Characterization of the Ghd8 Flowering Time Gene in a Mini-Core Collection of Miscanthus sinensis

The optimal flowering time for bioenergy crop Miscanthus is essential for environmental adaptability and biomass accumulation. However, little is known about how genes controlling flowering in other grasses contribute to flowering regulation in Miscanthus. Here, we report on the sequence characterization and gene expression of Miscanthus sinensisGhd8, a transcription factor encoding a HAP3/NF-YB DNA-binding domain, which has been identified as a major quantitative trait locus in rice, with pleiotropic effects on grain yield, heading date and plant height. In M. sinensis, we identified two homoeologous loci, MsiGhd8A located on chromosome 13 and MsiGhd8B on chromosome 7, with one on each of this paleo-allotetraploid species’ subgenomes. A total of 46 alleles and 28 predicted protein sequence types were identified in 12 wild-collected accessions. Several variants of MsiGhd8 showed a geographic and latitudinal distribution. Quantitative real-time PCR revealed that MsiGhd8 expressed under both long days and short days, and MsiGhd8B showed a significantly higher expression than MsiGhd8A. The comparison between flowering time and gene expression indicated that MsiGhd8B affected flowering time in response to day length for some accessions. This study provides insight into the conserved function of Ghd8 in the Poaceae, and is an important initial step in elucidating the flowering regulatory network of Miscanthus.

geographic distribution↗
NASA NTRS2011

Stress-Driven Selection of Novel Phenotypes

A process has been developed that can confer novel properties, such as metal resistance, to a host bacterium. This same process can also be used to produce RNAs and peptides that have novel properties, such as the ability to bind particular compounds. It is inherent in the method that the peptide or RNA will behave as expected in the target organism. Plasmid-born mini-gene libraries coding for either a population of combinatorial peptides or stable, artificial RNAs carrying random inserts are produced. These libraries, which have no bias towards any biological function, are used to transform the organism of interest and to serve as an initial source of genetic variation for stress-driven evolution. The transformed bacteria are propagated under selective pressure in order to obtain variants with the desired properties. The process is highly distinct from in vitro methods because the variants are selected in the context of the cell while it is experiencing stress. Hence, the selected peptide or RNA will, by definition, work as expected in the target cell as the cell adapts to its presence during the selection process. Once the novel gene, which produces the sought phenotype, is obtained, it can be transferred to the main genome to increase the genetic stability in the organism. Alternatively, the cell line can be used to produce novel RNAs or peptides with selectable properties in large quantity for separate purposes. The system allows for easy, large-scale purification of the RNAs or peptide products. The process has been reduced to practice by imposing sub-inhibitory concentrations of NiCl2 on cells of the bacterium Escherichia coli that were transformed separately with the peptide library and RNA library. The evolved resistant clones were isolated, and sequences of the selected mini-gene variants were established. Clones resistant to NiCl2 were found to carry identical plasmid variants with a functional mini-gene that specifically conferred significant nickel tolerance on the host cells. Sequencing of the selected mini-gene revealed a propensity of the encoded peptide to bind transient metal ions. Expression of the mini-gene markedly improved growth parameters of the evolved clones at sub-inhibitory concentrations of NiCl2 while being slightly detrimental in the absence of stress. Similar results have been obtained with the RNA libraries. Overall, the results demonstrate a very natural outcome of the selection experiments in which the mini-genes were expected to be either successfully integrated into bacterial genetic networks, or rejected depending upon their effect on host fitness. This described approach can be useful as a laboratory model to study the dynamics of bacterial adaptive evolution on the molecular level. It can also provide a strategy for screening expressed DNA libraries in search of novel genes with desirable properties.

Fox, George E.↗
DOE OSTI2022

Kingdom-Wide Analysis of Fungal Protein-Coding and tRNA Genes Reveals Conserved Patterns of Adaptive Evolution

Protein-coding genes evolved codon usage bias due to the combined but uneven effects of adaptive and nonadaptive influences. Studies in model fungi agree on codon usage bias as an adaptation for fine-tuning gene expression levels; however, such knowledge is lacking for most other fungi. Our comparative genomics analysis of over 450 species supports codon usage and transfer RNAs (tRNAs) as coadapted for translation speed and this is most likely a realization of convergent evolution. Rather than drift, phylogenetic reconstruction inferred adaptive radiation as the best explanation for the variation of interspecific codon usage bias. Although the phylogenetic signals for individual codon and tRNAs frequencies are lower than expected by genetic drift, we found remarkable conservation of highly expressed genes being codon optimized for translation by the most abundant tRNAs, especially by inosine-modified tRNAs. As an application, we present a sequence-to-expression neural network that uses codons to reliably predict highly expressed transcripts. The kingdom Fungi, with over a million species, includes many key players in various ecosystems and good targets for biotechnology. Collectively, our results have implications for better understanding the evolutionary success of fungi, as well as informing the biosynthetic manipulation of fungal genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Metabolome and RNA-seq Analysis of Responses to Nitrogen Deprivation and Resupply in Tea Plant (Camellia sinensis) Roots

Nitrogen (N) is an important contributor in regulating plant growth and development as well as secondary metabolites synthesis, so as to promote the formation of tea quality and flavor. Theanine, polyphenols, and caffeine are important secondary metabolites in tea plant. In this study, the responses of Camellia sinensis roots to N deprivation and resupply were investigated by metabolome and RNA-seq analysis. N deficiency induced content increase for most amino acids (AAs) and reduction for the remaining AAs, polyphenols, and caffeine. After N recovery, the decreased AAs and polyphenols showed a varying degree of recovery in content, but caffeine did not. Meanwhile, theanine increased in content, but its related synthetic genes were down-regulated, probably due to coordination of the whole N starvation regulatory network. Flavonoids-related pathways were relatively active following N stress according to KEGG enrichment analysis. Gene co-expression analysis revealed TCS2 , AMT1;1 , TAT2 , TS , and GOGAT as key genes, and TFs like MYB, bHLH, and NAC were also actively involved in N stress responses in C. sinensis roots. These findings facilitate the understanding of the molecular mechanism of N regulation in tea roots and provide genetic reference for improving N use efficiency in tea plant.

Xu, Wenluan↗
NASA NTRSDate not supplied

Cardiovascular Responses to Simulated Spaceflight: Molecular Signatures and Surrogate Outputs to Measure CVD Risk

During extended space missions beyond low Earth orbit, astronauts will encounter prolonged periods of weightlessness and low dose space radiation. Previous studies have shown that exposure to small doses of high LET radiation (< 50 cGy) can lead to both short-term and long-term alterations in heart function, structure and underlying molecular mechanisms. In this study, we aim to identify the molecular signature associated with the cardiovascular response to simulated galactic cosmic radiation (5-ion GCR) alone or in combination with simulated weightlessness at time intervals relevant to mission length and recovery. Additionally, we aim to determine whether sex impacts cardiovascular responses to these spaceflight factors. Our overarching goal is to enhance our understanding of the cardiovascular risks associated with extended space missions and the clinical endpoints they suggest. We hypothesize that exposure to simulated space radiation leads to enduring alterations in the transcriptome, redox signaling and cytokine environment of cardiovascular tissue, some which have known links with reduced cardiovascular performance, aging, and increased risk of cardiovascular disease (CVD). Furthermore, we posit that simulated space radiation exposure in combination with simulated microgravity exacerbates cardiovascular deficits compared to single factor exposure. Female and male C57BL/6J mice, aged 23-24 weeks, were exposed to a single dose of 5, 15, or 50 cGy of 5-ion GCR, or sham-treated (0 cGy). Euthanasia was performed at 14 days and ~4 months post-irradiation. Hearts, aorta and blood plasma were collected shortly thereafter. RNA-sequencing of left ventricles at ~4 months post-GCR exposure revealed sex differences in the heart transcriptome with a few genes showing radiation-dependent changes in expression levels. Notably, some of the differentially expressed genes in 15 and 50 cGy GCR groups are known to play roles in the development of CVD. Analysis of protein levels of a subset of inflammatory cytokines in the heart indicated sex differences but no differences between sham and 50 cGy groups. Results also showed correlations among differentially expressed genes and a subset of inflammatory cytokines, with some correlations altered by GCR exposure. These findings suggest that GCR exposure can modify protein and gene networks linked to inflammation and CVD progression. In the aorta, telomere lengths were comparable across treatment groups sexes. Mitochondrial copy number is a biomarker for mitochondrial function with decreased copy numbers associated with cardiometabolic disease traits. Mitochondrial copy numbers of aorta also showed no sex nor dose differences. In a second study, mice underwent one week of simulated microgravity by hindlimb unloading (HU) and then exposed to a single dose of 15 cGy of 5-ion GCR. HU was conducted for an additional two weeks following GCR exposure. Single factor exposure groups (HU or GCR only) also were included in the study. Euthanasia was then performed and the same tissues were collected. Protein levels of select inflammatory cytokines in the heart showed sex-dependent differences in expression. In the aorta, telomere lengths and mitochondrial copy number also showed sex differences. In summary, our results indicate differences between sexes in biomarkers related to cardiovascular health. Exposure to 5-ion GCR or HU, alone or in combination, did not result in changes in most of the cardiovascular biomarkers that were examined. However, in the heart, simulated space radiation at doses of 15 and 50 cGy led to long-term alterations in the expression levels of a small group of genes known to be associated with the progression of CVD. The long-term transcriptomic changes resulting from exposure to simulated space radiation should be carefully investigated to mitigate adverse cardiovascular events during and after deep space missions. Our results also highlight the importance of sex-specific strategies in monitoring and maintaining cardiovascular health during and after deep space missions.

cardiovascular↗
DOE OSTI2023

The Mycobacterium tuberculosis protein O -phosphorylation landscape

Bacterial phosphosignaling has long been synonymous with the histidine kinases of the two component systems (TCSs), but many bacteria, including Mycobacterium tuberculosis (Mtb), also code for Ser/Thr protein kinases (STPKs). STPKs are the main phosphosignaling enzymes in eukaryotes, but the full extent of phosphorylation on protein Ser/Thr and Tyr (O-phosphorylation) in bacteria is unclear. Here, we explored the global signaling capacity of the STPKs in Mtb. We generated a panel of STPK loss- and gain-of-function strains and measured the resulting O-phosphorylation and transcriptional changes. This deep phosphoproteome shows that O-phosphorylation in Mtb is a vastly underexplored protein modification that affects >70% of the proteome. The substrate-kinase interactions of STPKs show an extensive interface with the transcriptional machinery, resulting in STPK regulation of gene expression of over 30% of Mtb genes. Finally, Mtb O-phosphorylation gives rise to an expansive, distributed, and cooperative network of a complexity that has previously only been associated with eukaryotic phosphosignaling networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Identification of a diguanylate cyclase expressed in the presence of plants and its application for discovering candidate gene products involved in plant colonization by Pantoea sp. YR343

Microbial colonization of plant roots is a highly complex process that requires the coordination and regulation of many gene networks, yet the identities and functions of many of these gene products have yet to be discovered. Pantoea sp. YR343, a gamma-proteobacterium isolated from the rhizosphere of Populus deltoides , forms robust biofilms along the root surfaces of Populus and possesses plant growth-promoting characteristics. In this work, we identified three diguanylate cyclases in the plant-associated microbe Pantoea sp. YR343 that are expressed in the presence of plant roots. One of these diguanylate cyclases, DGC2884, localizes to discrete sites in the cells and its overexpression results in reduced motility and increased EPS production and biofilm formation. We performed a genetic screen by expressing this diguanylate cyclase from an inducible promoter in order to identify candidate gene products that may be involved in root colonization by Pantoea sp. YR343. Further, we demonstrate the importance of other domains in DGC2884 to its activity, which in combination with the genes identified by transposon mutagenesis, may yield insights into the mechanisms of plant association as well as the activity and regulation of homologous enzymes in medically and agriculturally relevant microbes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Global transcriptome profile of the developmental principles of in vitro iPSC-to-motor neuron differentiation

Abstract Background Human induced pluripotent stem cells (iPSC) have opened new avenues for regenerative medicine. Consequently, iPSC-derived motor neurons have emerged as potentially viable therapies for spinal cord injuries and neurodegenerative disorders including Amyotrophic Lateral Sclerosis. However, direct clinical application of iPSC bears in itself the risk of tumorigenesis and other unforeseeable genetic or epigenetic abnormalities. Results Employing RNA-seq technology, we identified and characterized gene regulatory networks triggered by in vitro chemical reprogramming of iPSC into cells with the molecular features of motor neurons (MNs) whose function in vivo is to innervate effector organs. We present meta-transcriptome signatures of 5 cell types: iPSCs, neural stem cells, motor neuron progenitors, early motor neurons, and mature motor neurons. In strict response to the chemical stimuli, along the MN differentiation axis we observed temporal downregulation of tumor growth factor-β signaling pathway and consistent activation of sonic hedgehog, Wnt/β-catenin, and Notch signaling. Together with gene networks defining neuronal differentiation (neurogenin 2, microtubule-associated protein 2, Pax6, and neuropilin-1), we observed steady accumulation of motor neuron-specific regulatory genes, including Islet-1 and homeobox protein HB9. Interestingly, transcriptome profiling of the differentiation process showed that Ca 2+ signaling through cAMP and LPC was downregulated during the conversion of the iPSC to neural stem cells and key regulatory gene activity of the pathway remained inhibited until later stages of motor neuron formation. Pathways shaping the neuronal development and function were well-represented in the early motor neuron cells including, neuroactive ligand-receptor interactions, axon guidance, and the cholinergic synapse formation. A notable hallmark of our in vitro motor neuron maturation in monoculture was the activation of genes encoding G-coupled muscarinic acetylcholine receptors and downregulation of the ionotropic nicotinic acetylcholine receptors expression. We observed the formation of functional neuronal networks as spontaneous oscillations in the extracellular action potentials recorded on multi-electrode array chip after 20 days of differentiation. Conclusions Detailed transcriptome profile of each developmental step from iPSC to motor neuron driven by chemical induction provides the guidelines to novel therapeutic approaches in the re-construction efforts of muscle innervation.

59 BASIC BIOLOGICAL SCIENCES↗