Process Gene Ontology Research Articles

Regulatory Effects of CsrA in Vibrio cholerae.

CsrA is a posttranscriptional global regulator in Vibrio cholerae Although CsrA is critical for V. cholerae survival within the mammalian host, the regulatory targets of CsrA remain mostly unknown. To identify pathways controlled by CsrA, RNA-seq transcriptome analysis was carried out by comparing the wild type and the csrA mutant grown to early exponential, mid-exponential, and stationary phases of growth. This enabled us to identify the global effects of CsrA-mediated regulation throughout the V. cholerae growth cycle. We found that CsrA regulates 22% of the V. cholerae transcriptome, with significant regulation within the gene ontology (GO) processes that involve amino acid transport and metabolism, central carbon metabolism, lipid metabolism, iron uptake, and flagellum-dependent motility. Through CsrA-RNA coimmunoprecipitation experiments, we found that CsrA binds to multiple mRNAs that encode regulatory proteins. These include transcripts encoding the major sigma factors RpoS and RpoE, which may explain how CsrA regulation affects such a large proportion of the V. cholerae transcriptome. Other direct targets include flrC, encoding a central regulator in flagellar gene expression, and aphA, encoding the virulence gene transcription factor AphA. We found that CsrA binds to the aphA mRNA both in vivo and in vitro, and CsrA significantly increases AphA protein synthesis. The increase in AphA was due to increased translation, not transcription, in the presence of CsrA, consistent with CsrA binding to the aphA transcript and enhancing its translation. CsrA is required for the virulence of V. cholerae and this study illustrates the central role of CsrA in virulence gene regulation.IMPORTANCEVibrio cholerae, a Gram-negative bacterium, is a natural inhabitant of the aqueous environment. However, once ingested, this bacterium can colonize the human host and cause the disease cholera. In order to successfully transition between its aqueous habitat and the human host, the bacterium must sense changes in its environment and rapidly alter gene expression. Global regulators, including CsrA, play an integral role in altering the expression of a large number of genes to promote adaptation and survival, which is required for intestinal colonization. We used transcriptomics and a directed CsrA-RNA coimmunoprecipitation to characterize the CsrA regulon and found that CsrA alters the expression of more than 800 transcripts in V. cholerae Processes regulated by CsrA include motility, the rugose phenotype, and virulence pathways. CsrA directly binds to the aphA transcript and positively regulates the production of the virulence regulator AphA. Thus, CsrA regulates multiple processes that have been linked to pathogenesis.

The Bromodomain Containing 8 (BRD8) transcriptional network in human lung epithelial cells

Mechanisms regulating gene expression in the airway epithelium underlie its response to the environment. A network of transcription factors (TFs) and architectural proteins, modulate chromatin accessibility and recruit activating or repressive signals. Bromodomain-containing proteins function as TFs or by engaging methyltransferase or acetyltransferase activity to induce chromatin modifications. Here we investigate the role of Bromodomain Containing 8 (BRD8) in coordinating lung epithelial function. Sites of BRD8 occupancy genome-wide were mapped in human lung epithelial cell lines (Calu-3 and 16HBE14o-). CCCTC-Binding Factor (CTCF) was identified as a predicted co-factor of BRD8, based upon motif over-representation under BRD8 ChIP-seq peaks. Following siRNA-mediated depletion of BRD8, differentially expressed genes with nearby peaks of BRD8 occupancy were subject to gene ontology process enrichment analysis. BRD8 targets are enriched for genes involved in the innate immune response and the cell cycle. Depletion of BRD8 increased the secretion of the antimicrobial peptide beta-defensin 1 and multiple chemokines, and reduced cell proliferation.

Lignin Synthesis Related Genes with Potential Significance in the Response of Upland Cotton to Fusarium Wilt Identified by Transcriptome Profiling

Fusarium wilt, caused by the fungus Fusarium oxysporum Schlecht. f. sp. vasinfectum (Fov) is a destructive soil-borne cotton disease. To profile the genes and pathways responding to Fov infection, we compared transcriptomic responses before and after F. oxysporum inoculation in a highly resistant cotton cultivar, Yumian21, and a highly susceptible cultivar, Jimian11. Although the overall gene expression pattern was downregulated in both cultivars, the global gene expression in the resistant cultivar was stronger than that in the susceptible cultivar. In addition, the expressed genes of two cultivars mostly differed in “cellular process,” “single-organism process,” “metabolic process,” and “response to stimulus” functional groups in the biological process Gene Ontology category: the upregulated differentially expressed genes (DEG) were largely enriched in the resistant cultivar, while the downregulated DEGs were largely enriched in the susceptible cultivar. Phenylpropanoid biosynthesis and phenylalanine metabolism are the key metabolic pathways in cotton in response to Fov. We found that lignin plays a potential role in cotton resistance to Fov. Two coding genes, caffeic acid 3-O-methyltransferase and peroxidase2, as well as the two transcription factors MYB46 and MYB86, are possibly involved in the accumulation and synthesis of lignin. Furthermore, the result showed that the quantification of lignin could be potentially used as a selection tool to identify Fusarium wilt resistant cotton.

Identification of core miRNAs and regulatory pathways in breast cancer by integrated bioinformatics analysis.

Breast cancer (BC) ranks first among malignancies in the female population due to its complicated pathological progression and poor prognosis. Hence, the aim of the present study was to identify potential molecular prognostic biomarkers able to predict the prognosis of BC patients. We integrated two microRNA (miRNA) expression microarrays and three gene microarrays related to BC from the NCBI Gene Expression Comprehensive (GEO) database to screen for differentially expressed miRNAs and identify their regulatory networks. The Kaplan-Meier plotter online analysis tool was used to assess the overall survival value of miRNAs expression in BC patients. The LinkedOmics online tool was used to analyze genes correlated with miRNAs expression. To clarify the upstream regulation mechanism of genes, we used ChIP-Atlas to identify and screen for transcription factors and visually verify them using the Integrative Genomics Viewer. To further analyze the downstream regulatory mechanism of miRNA in BC, we verified differentially expressed genes (DEGs) correlated to miRNAs in three GEO gene microarrays and the gene set predicted by miRWalk. The open access Metascape program allowed analysis of Gene Ontology (GO) processes, KEGG pathways and GO enrichment was performed on the DEGs. To further identify hub genes, Cytoscape software and its plug-in were applied to construct protein-protein interaction networks. In the present study, several possible molecules and related pathways related to miR-483 were identified by bioinformatics analysis. These molecules and pathways might represent key mechanisms involved in BC progression and development. This work provides a novel view and insight in the pathogenesis, treatment and prognosis for BC.

Preliminary comparison of skin transcriptome from sheep with different wool fibre diameters

Context Wool is one of the most important animal fibres for the textile industry, and its diameter directly affects its economic value. However, the molecular mechanisms underlying wool fibre diameter (FD) in sheep have not been fully elucidated. Aims The aims of the work were to make an initial comparison of skin transcriptomes from sheep with wool of high and low FD, and to identify key genes affecting FD. Methods High-throughput RNA-Seq technology was employed to explore the skin transcriptome, using three sheep with fine wool (FD <21.0 μm) and three sheep with coarse wool (FD >27.0 μm). Key results We obtained 28 607 228 bp of clean sequence data, 78.9% (±3.8%) of which uniquely aligned to the reference genome across the six samples. In total, 19 914 mRNA transcripts were expressed (FPKM >0) in the six skin samples, among which were certain well-known genes involved in the skin–hair cycle, such as KRTAP7-1, KRT14, Wnt10b, Wnt2b, β-catenin and FGF5. Furthermore, 467 genes were significantly differentially expressed between the fine-wool and coarse-wool groups, including 21 genes with upregulated and 446 genes with downregulated expression in the sheep with lower FD. These differentially expressed genes were particularly enriched in the gene ontology processes related to lipid metabolism, skin development, differentiation and immune function (P < 0.05). The biological processes were involved in collagen catabolism, negative regulation of macromolecule metabolism, steroid hormone stimulation and lipid metabolism. A significant Kyoto Encyclopedia of Genes and Genomes pathway involving the metabolism of lipids and lipoproteins was also enriched, revealing that lipid metabolism might be one of the key factors affecting FD. The expression of these differentially expressed genes that were involved in the metabolism of lipids and lipoproteins pathway was verified by quantitative real-time PCR (qPCR). The correlation between the mRNA expression level from qPCR and RNA-Seq data was 0.999 (P < 0.001). Conclusions The 467 differentially expressed genes, especially those involved in lipid metabolism and immune function, may play key roles in wool follicle metabolism and the expression of wool FD. Implications This study provided valuable data for future studies aimed at elucidating the mechanisms that underlie wool follicle metabolism and wool FD. The work may also have implications for studies of the human hair follicle.

Comparative Analysis of Systemic and Tumor Microenvironment Proteomes From Children With B-Cell Acute Lymphocytic Leukemia at Diagnosis and After Induction Treatment.

Among the childhood diseases, B-cell acute lymphocytic leukemia (B-ALL) is the most frequent type of cancer. Despite recent advances concerning disease treatment, cytotoxic chemotherapy remains the first line of treatment in several countries, and the modifications induced by such drugs in the organism are still poorly understood. In this context, the present study provided a comparative high-throughput proteomic analysis of the cumulative changes induced by chemotherapeutic drugs used in the induction phase of B-ALL treatment in both peripheral blood (PB) and bone marrow compartment (BM) samples. To reach this goal, PB and BM plasma samples were comparatively analyzed by using label-free proteomics at two endpoints: at diagnosis (D0) and the end of the cumulative induction phase treatment (D28). Proteomic data was available via ProteomeXchange with identifier PXD021584. The resulting differentially expressed proteins were explored by bioinformatics approaches aiming to identify the main gene ontology processes, pathways, and transcription factors altered by chemotherapy, as well as to understand B-ALL biology in each compartment at D0. At D0, PB was characterized as a pro-inflammatory environment, with the involvement of several downregulated coagulation proteins as KNG, plasmin, and plasminogen. D28 was characterized predominantly by immune response-related processes and the super expression of the transcription factor IRF3 and transthyretin. RUNX1 was pointed out as a common transcription factor found in both D0 and D28. We chose to validate the proteins transthyretin and interferon-gamma (IFN-γ) by commercial kits and expressed the results as PB/BM ratios. Transthyretin ratio was augmented after induction chemotherapy, while IFN-γ was reduced at the end of the treatment. Considering that most of these proteins were not yet described in B-ALL literature, these findings added to understanding disease biology at diagnosis and highlighted a possible role for transthyretin and IFN-γ as mechanisms related to disease resolution.

High levels of glucose alter Physcomitrella patens metabolism and trigger a differential proteomic response.

Sugars act not only as substrates for plant metabolism, but also have a pivotal role in signaling pathways. Glucose signaling has been widely studied in the vascular plant Arabidopsis thaliana, but it has remained unexplored in non-vascular species such as Physcomitrella patens. To investigate P. patens response to high glucose treatment, we explored the dynamic changes in metabolism and protein population by applying a metabolomic fingerprint analysis (DIESI-MS), carbohydrate and chlorophyll quantification, Fv/Fm determination and label-free untargeted proteomics. Glucose feeding causes specific changes in P. patens metabolomic fingerprint, carbohydrate contents and protein accumulation, which is clearly different from those of osmotically induced responses. The maximal rate of PSII was not affected although chlorophyll decreased in both treatments. The biological process, cellular component, and molecular function gene ontology (GO) classifications of the differentially expressed proteins indicate the translation process is the most represented category in response to glucose, followed by photosynthesis, cellular response to oxidative stress and protein refolding. Importantly, although several proteins have high fold changes, these proteins have no predicted identity. The most significant discovery of our study at the proteome level is that high glucose increase abundance of proteins related to the translation process, which was not previously evidenced in non-vascular plants, indicating that regulation by glucose at the translational level is a partially conserved response in both plant lineages. To our knowledge, this is the first time that metabolome fingerprint and proteomic analyses are performed after a high sugar treatment in non-vascular plants. These findings unravel evolutionarily shared and differential responses between vascular and non-vascular plants.

A thorough analysis of the contribution of experimental, derived and sequence-based predicted protein-protein interactions for functional annotation of proteins.

Physical interaction between two proteins is strong evidence that the proteins are involved in the same biological process, making Protein-Protein Interaction (PPI) networks a valuable data resource for predicting the cellular functions of proteins. However, PPI networks are largely incomplete for non-model species. Here, we tested to what extent these incomplete networks are still useful for genome-wide function prediction. We used two network-based classifiers to predict Biological Process Gene Ontology terms from protein interaction data in four species: Saccharomyces cerevisiae, Escherichia coli, Arabidopsis thaliana and Solanum lycopersicum (tomato). The classifiers had reasonable performance in the well-studied yeast, but performed poorly in the other species. We showed that this poor performance can be considerably improved by adding edges predicted from various data sources, such as text mining, and that associations from the STRING database are more useful than interactions predicted by a neural network from sequence-based features.

Transcriptome analysis for molecular landscaping of genes controlling diterpene andrographolide biosynthesis in Andrographis paniculata (Burm . f.) Nees.

Kalmegh [Andrographis paniculata (Burm. f.) Nees.] is one of the essential medicinal plants due to an important terpenoid, i.e. andrographolide possesses immense therapeutic and pharmacological uses. The experiment was performed to elucidate the expression of candidate genes associated with andrographolide biosynthesis. Based on results obtained in chromatography for andrographolide content analysis of six genotypes, two contrast genotypes, i.e. IC-520361 (maximum andrographolide content—2.33%) and Anand Local (lowest andrographolide content—1.01%) were selected for the transcriptome analysis. A total of 1.04 Gb of raw data were produced using MiSeq Illumina platform, in which IC 520361 generated 645 million base pairs sequence along with 4,524,251 raw reads and Anand Local produced 419 million base pairs sequence along with 3,021,316 raw reads. The combined assembly of high quality reads generated for both the samples had 33,247,454 bp of total assembled bases and 38,292 of transcripts. The GC percent of assembled transcripts was 44.79%, an average read length was 800 bp and N50 value was 1186 bp. Species-specific distribution using BLAST X (Nr), showed the highest Blast hits with Sesamum indicum. Out of 23,346 transcripts, 87% of transcripts annotated in UniProt KB (Universal Protein Resource KnowledgeBase) database and only 0.21% of transcripts were annotated in TAIR (The Arabidopsis Information Resources). Biological processes gene ontology classified based on Blast2GO showed, out of 6853 transcripts, 1370 of transcripts were represented by terpenoid biosynthetic pathway, which involved in secondary metabolite andrographolide biosynthesis. The heat map showed 1016 transcripts were differentially expressed between two kalmegh genotypes, in which nine important differentially expressed transcripts related to MEP (2C methyl-d-erythritol 4-phosphate) and MVA (Mevalonic acid) andrographolide biosynthesis pathways such as, geranyl diphosphate synthase small subunit, Isopentenyl-diphosphate delta-isomerase i-like, 4, 13-hydroxy-3-methylglutaryl-coenzyme a reductase etc. were upregulated in IC 520361 as compared to Anand Local, which were validated through RT-qPCR. The highest expression of gene 13-hydroxy-3-methylglutaryl-coenzyme a reductase (HMGR) was reported, which is responsible for accumulation of andrographolide in leaf. This comparative transcriptome analysis confirmed the expression level of genes were higher in accession IC 520361 as compare to Anand Local related to andrographolide biosynthesis pathways i.e. MEP and MVA. These up-regulated genes could be over-expressed to enhance the andrographolide content using genetic engineering of these metabolic pathways. It will also give an idea to the breeder for development of molecular markers for direct screening of the genotypes.Electronic supplementary materialThe online version of this article (10.1007/s13205-020-02511-y) contains supplementary material, which is available to authorized users.

Rare Functional Variants Associated with Antidepressant Remission in Mexican-Americans: Short title: Antidepressant remission and pharmacogenetics in Mexican-Americans

Rare genetic functional variants can contribute to 30-40% of functional variability in genes relevant to drug action. Therefore, we investigated the role of rare functional variants in antidepressant response. Mexican-American individuals meeting the Diagnostic and Statistical Manual-IV criteria for major depressive disorder (MDD) participated in a prospective randomized, double-blind study with desipramine or fluoxetine. The rare variant analysis was performed using whole-exome genotyping data. Network and pathway analyses were carried out with the list of significant genes. The Kernel-Based Adaptive Cluster method identified functional rare variants in 35 genes significantly associated with treatment remission (False discovery rate, FDR <0.01). Pathway analysis of these genes supports the involvement of the following gene ontology processes: olfactory/sensory transduction, regulation of response to cytokine stimulus, and meiotic cell cycleprocess. Our study did not have a placebo arm. We were not able to use antidepressant blood level as a covariate. Our study is based on a small sample size of only 65 Mexican-American individuals. Further studies using larger cohorts are warranted. Our data identified several rare functional variants in antidepressant drug response in MDD patients. These have the potential to serve as genetic markers for predicting drug response. ClinicalTrials.gov NCT00265291.

Multi-omics analysis suggests enhanced epileptogenesis in the Cornu Ammonis 3 of the pilocarpine model of mesial temporal lobe epilepsy.

Mesial temporal lobe epilepsy (MTLE) is a chronic neurological disorder characterized by the occurrence of seizures, and histopathological abnormalities in the mesial temporal lobe structures, mainly hippocampal sclerosis (HS). We used a multi-omics approach to determine the profile of transcript and protein expression in the dorsal and ventral hippocampal dentate gyrus (DG) and Cornu Ammonis 3 (CA3) in an animal model of MTLE induced by pilocarpine. We performed label-free proteomics and RNAseq from laser-microdissected tissue isolated from pilocarpine-induced Wistar rats. We divided the DG and CA3 into dorsal and ventral areas and analyzed them separately. We performed a data integration analysis and evaluated enriched signaling pathways, as well as the integrated networks generated based on the gene ontology processes. Our results indicate differences in the transcriptomic and proteomic profiles among the DG and the CA3 subfields of the hippocampus. Moreover, our data suggest that epileptogenesis is enhanced in the CA3 region when compared to the DG, with most abnormalities in transcript and protein levels occurring in the CA3. Furthermore, our results show that the epileptogenesis in the pilocarpine model involves predominantly abnormal regulation of excitatory neuronal mechanisms mediated by N-methyl D-aspartate (NMDA) receptors, changes in the serotonin signaling, and neuronal activity controlled by calcium/calmodulin-dependent protein kinase (CaMK) regulation and leucine-rich repeat kinase 2 (LRRK2)/WNT signaling pathways.

Bio-activating ultrafine grain titanium: RNA sequencing reveals enhanced mechano-activation of osteoconduction on nanostructured substrates.

Titanium is essentially absent from biological systems yet reliably integrates into bone. To achieve osseointegration, titanium must activate biological processes without entering cells, defining it as a bio-activating material. Nanostructuring bulk titanium reduces grain size, increases strength, and improves other quantifiable physical properties, including cytocompatibility. The biological processes activated by increasing grain boundary availability were detected with total RNA-sequencing in mouse pre-osteoblasts grown for 72 hours on nanometrically smooth substrates of either coarse grain or nanostructured ultrafine grain titanium. The average grain boundary length under cells on the conventional coarse grain substrates is 273.0 μm, compared to 70,881.5 μm for cells adhered to the nanostructured ultrafine grain substrates; a 260-fold difference. Cells on both substrates exhibit similar expression profiles for genes whose products are critical for mechanosensation and transduction of cues that trigger osteoconduction. Biological process Gene Ontology term enrichment analysis of differentially expressed genes reveals that cell cycle, chromatin modification, telomere maintenance, and RNA metabolism processes are upregulated on ultrafine grain titanium. Processes related to immune response, including apoptosis, are downregulated. Tumor-suppressor genes are upregulated while tumor-promoting genes are downregulated. Upregulation of genes involved in chromatin remodeling and downregulation of genes under the control of the peripheral circadian clock implicate both processes in the transduction of mechanosensory information. Non-coding RNAs may also play a role in the response. Merging transcriptomics with well-established mechanobiology principles generates a unified model to explain the bio-activating properties of titanium. The modulation of processes is accomplished through chromatin remodeling in which the nucleus responds like a rheostat to grain boundary concentration. This convergence of biological and materials science reveals a pathway toward understanding the biotic-abiotic interface and will inform the development of effective bio-activating and bio-inactivating materials.

Protein Expression Analysis of an In Vitro Murine Model of Prostate Cancer Progression: Towards Identification of High-Potential Therapeutic Targets.

Background: Prostate cancer (PC) is the most frequently diagnosed cancer among men worldwide. The poor prognosis of PC is largely due to late diagnosis of the disease when it has progressed to advanced stages marked by androgen-independence. We interrogated proteomic signatures that embody the transition of PC from an androgen-dependent (AD) to an androgen-independent (AI) state. Methods: We have previously established AD and AI murine PC cell lines, PLum-AD and PLum-AI, respectively, which recapitulate primary and progressive PC at phenotypic and subcellular levels. We statistically surveyed global protein expression profiles in these cell lines. Differential profiles were functionally interrogated by pathways and protein–protein interaction network analyses. Results: Protein expression pattern analysis revealed a total of 683 proteins, among which 99 were significantly differentially altered in PLum-AI cells as compared to PLum-AD cells (45 increased and 54 decreased). Principal component analysis (PCA) revealed that the two different cell lines clearly separated apart, indicating a significant proteome expression difference between them. Four of the proteins (vimentin, catalase, EpCAM, and caspase 3) that were differentially expressed in PLum-AI cells compared to PLum-AD cells were subjected to biochemical validation by Western blotting. Biological process gene ontology (GO) analysis of the differentially expressed proteins demonstrated enrichment of biological functions and pathways in PLum-AI cells that are central to PI3 kinase and androgen receptor pathways. Besides, other relevant biological processes that are enriched in PLum-AI cells included cell adhesion and cell migration processes, cell and DNA damage, apoptosis, and cell cycle regulation. Conclusions: Our protein expression analysis of a murine in vitro model of PC progression identified differential protein spots that denote this progression and that comprise high-potential targets for early treatment of PC with a personalized patient-specific approach. Efforts are underway to functionally assess the potential roles of these proteins as therapeutic targets for PC progression.

Transcriptome Profiling Associated with Carcass Quality of Loin Muscles in Crossbred Pigs.

Simple SummaryCarcass quality traits, such as lean depth and loin depth, are of extreme economic importance for the swine industry. This study aimed to identify the gene expression pattern related to carcass quality in crossbred pigs ((Landrace × Yorkshire) × Duroc). In total, 20 crossbred pigs were used in this study and divided into two groups based on the loin muscle quality grade. Total RNA samples extracted from the loin muscles of both groups were submitted for RNA-seq. Differentially expressed gene (DEG) analysis of the two groups revealed 282 up-regulated and 189 down-regulated genes (p ≤ 0.01), linked to tissue development, striated muscle tissue development, tissue morphogenesis, and lipid metabolic process gene ontology (GO) terms. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis highlighted genes related to the calcium signaling pathway, apelin signaling pathway, and the mammalian target of rapamycin (mTOR) signaling pathway. We constructed an expressed gene catalog, which may serve as a resource for genomic studies focused on uncovering the molecular mechanisms underlying carcass quality in crossbred pigs.Carcass quality traits, such as lean depth and loin depth, are of extreme economic importance for the swine industry. This study aimed to identify the gene expression pattern related to carcass quality in crossbred pigs ((Landrace × Yorkshire) × Duroc). In total, 20 crossbred pigs were used in this study and they were divided into two groups (class I grade, n = 10; class II grade, n = 10) based on the carcass grades. Total RNA samples extracted from the loin muscles of both groups were submitted for RNA-seq. The quality assessment of the sequencing reads resulted in 25,458 unigenes and found 12,795 candidate coding unigenes with homology to other species after annotation. Differentially expressed gene (DEG) analysis of the two groups revealed 282 up-regulated and 189 down-regulated genes (p ≤ 0.01), linked to tissue development, striated muscle tissue development, tissue morphogenesis, and lipid metabolic process gene ontology (GO) terms. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis highlighted genes related to the calcium signaling pathway, melanogenesis, the sphingolipid signaling pathway, the apelin signaling pathway, and the mTOR signaling pathway. We constructed an expressed gene profile, which may serve as a resource for genomic studies focused on uncovering the molecular mechanisms underlying carcass quality in crossbred pigs.

Effect of Nutritional Restriction on the Hair Follicles Development and Skin Transcriptome of Chinese Merino Sheep.

Simple SummaryThe morphogenesis hair follicles begins in the embryonic stage. The first follicles formed are the primary wool follicles (PF), followed by secondary wool follicles (SF), and then secondary-derived follicles that branch from the SF in Merino sheep. The morphogenesis and development of SF in sheep determine wool yield and quality. At present, the mechanism of SF development and branching is still unclear. In this study, hair follicle morphogenesis and development were inhibited through nutritional restrictions during pregnancy, and transcriptome sequencing was performed to screen for differentially expressed genes related to the SF development in the fetus. Our findings are helpful to understand the mechanisms of SF development and branching of Merino sheep.The high concentration of secondary branched wool follicles is a distinctive feature of the Merino sheep. At present, the molecular control of the development and branching of secondary wool follicles (SF) remains elusive. To reveal the potential genes associated with the development of hair follicles, we investigated the characteristics of prenatal and postnatal development of wool follicles, and the transcriptional expression profile in fetuses/lambs from dams under either maternal maintenance or sub-maintenance (75% maintenance) nutrition. The density of SF and the ratio of SF to primary wool follicles (PF) were reduced (p < 0.05) in fetuses from day 105 to 135 of gestation under sub-maintenance nutrition. Differentially expressed genes were enriched in the binding, single-organism process, cellular process, cell and cell part Gene Ontology (GO) functional categories and metabolism, apoptosis, and ribosome pathways. Four candidate genes, SFRP4, PITX1, BAMBI, and KRT16, which were involved in secondary wool follicles branching and development, were identified. Our results indicate that nutritional intervention imposed on pregnant ewes by short-term sub-maintenance nutrition could provide a strategy for the study of wool follicle development. Overall insight into the global gene expression associated with SF development can be used to investigate the underlying mechanisms of SF branching in Merino sheep.

Old Mice Have Less Transcriptional Activation But Similar Periosteal Cell Proliferation Compared to Young-Adult Mice in Response to in vivo Mechanical Loading.

Mechanical loading is a potent strategy to induce bone formation, but with aging, the bone formation response to the same mechanical stimulus diminishes. Our main objectives were to (i) discover the potential transcriptional differences and (ii) compare the periosteal cell proliferation between tibias of young-adult and old mice in response to strain-matched mechanical loading. First, to discover potential age-related transcriptional differences, we performed RNA sequencing (RNA-seq) to compare the loading responses between tibias of young-adult (5-month) and old (22-month) C57BL/6N female mice following 1, 3, or 5 days of axial loading (loaded versus non-loaded). Compared to young-adult mice, old mice had less transcriptional activation following loading at each time point, as measured by the number of differentially expressed genes (DEGs) and the fold-changes of the DEGs. Old mice engaged fewer pathways and gene ontology (GO) processes, showing less activation of processes related to proliferation and differentiation. In tibias of young-adult mice, we observed prominent Wnt signaling, extracellular matrix (ECM), and neuronal responses, which were diminished with aging. Additionally, we identified several targets that may be effective in restoring the mechanoresponsiveness of aged bone, including nerve growth factor (NGF), Notum, prostaglandin signaling, Nell-1, and the AP-1 family. Second, to directly test the extent to which periosteal cell proliferation was diminished in old mice, we used bromodeoxyuridine (BrdU) in a separate cohort of mice to label cells that divided during the 5-day loading interval. Young-adult and old mice had an average of 15.5 and 16.7 BrdU+ surface cells/mm, respectively, suggesting that impaired proliferation in the first 5 days of loading does not explain the diminished bone formation response with aging. We conclude that old mice have diminished transcriptional activation following mechanical loading, but periosteal proliferation in the first 5 days of loading does not differ between tibias of young-adult and old mice. © 2020 American Society for Bone and Mineral Research.

Gene Expression in Patient-Derived Neural Progenitors Implicates WNT5A Signaling in the Etiology of Schizophrenia

BackgroundGenome-wide association studies of schizophrenia have demonstrated that variations in noncoding regions are responsible for most of the common variation heritability of the disease. It is hypothesized that these risk variants alter gene expression. Therefore, studying alterations in gene expression in schizophrenia may provide a direct approach to understanding the etiology of the disease. In this study we use cultured neural progenitor cells derived from olfactory neuroepithelium (CNON cells) as a genetically unaltered cellular model to elucidate the neurodevelopmental aspects of schizophrenia. MethodsWe performed a gene expression study using RNA sequencing of CNON cells from 111 control subjects and 144 individuals with schizophrenia. Differentially expressed genes were identified with DESeq2 software, using covariates to correct for sex, age, library batches, and 1 surrogate variable component. ResultsA total of 80 genes were differentially expressed (false discovery rate < 10%), showing enrichment in cell migration, cell adhesion, developmental process, synapse assembly, cell proliferation, and related Gene Ontology categories. Cadherin and Wnt signaling pathways were positive in overrepresentation test, and, in addition, many genes were specifically involved in WNT5A signaling. The differentially expressed genes were modestly, but significantly, enriched in the genes overlapping single nucleotide polymorphisms with genome-wide significant association from the Psychiatric Genomics Consortium genome-wide association study of schizophrenia. We also found substantial overlap with genes associated with other psychiatric disorders or brain development, enrichment in the same Gene Ontology categories as genes with mutations de novo in schizophrenia, and studies of induced pluripotent stem cell–derived neural progenitor cells. ConclusionsCNON cells are a good model of the neurodevelopmental aspects of schizophrenia and can be used to elucidate the etiology of the disorder.

Integrated analysis of a ceRNA network reveals potential prognostic lncRNAs in gastric cancer.

Long noncoding RNAs (lncRNAs) have important biological functions as competing endogenous RNAs (ceRNAs) in tumors, yet the functions and regulatory mechanisms of lncRNA‐related ceRNAs in gastric cancer have not been fully elucidated. In this study, we constructed a lncRNA‐miRNA‐mRNA ceRNA network and identified potential lncRNA biomarkers in gastric cancer. Basing on the RNA profiles downloaded from The Cancer Genome Atlas (TCGA) platform, the gastric cancer‐specific differentially expressed lncRNAs, miRNAs, and mRNAs were screened for constructing a ceRNA network using bioinformatic tools. The enrichment analysis of the biological processes in Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes pathways was performed on the ceRNA‐related DEmRNAs. According to the modularization of protein‐protein interaction (PPI) network, we extracted a ceRNA subnetwork and analyzed the correlation between the expression of the lncRNAs involved and specific clinical features of patients. Next, the expression of highly up‐regulated in liver cancer (HULC) and RP11‐314B1.2 showed significant changes in several pathological processes involved in gastric cancer, and nine lncRNAs were found to be correlated with the overall survival of patients with gastric cancer. Through the univariate and multivariate Cox regression analyses, two lncRNAs (LINC00106 and RP11‐999E24.3) were identified and utilized to establish a risk score model for assessing the prognosis of patients. The analysis results were also partially verified using quantitative real‐time PCR. The findings from this study indicate that HULC, RP11‐314B1.2, LINC00106, and RP11‐999E24.3 could be considered as potential therapeutic targets or prognostic biomarkers in gastric cancer, and provide a new perspective for cancer pathogenesis research.

Genome Analysis Reveals a Synergistic Mechanism of Ursodeoxycholic Acid and Jasminoidin in Mice Brain Repair After Ischemia/Reperfusion: Crosstalk Among Muti-Pathways.

Studies have shown that combination drug therapy which corresponding treatment involves multiple genes and targets is more effective against cerebral ischemia. To identify the synergistic mechanism of ursodeoxycholic acid and jasminoidin based on differential pathway network, which protect against brain ischemia-reperfusion injury. Totally 115 mice with focal cerebral ischemia-reperfusion injury were allocated into five groups: sham, vehicle, ursodeoxycholic acid (UA), jasminoidin (JA), and JA and UA combination group (JU). The differentially expressed genes identified by microarray which consisted of 11,644 complementary DNAs were loaded to the GeneGo MetaCore™ software to analyze the enriched pathways and processes among different groups. Of the top 10 pathways and process networks, 5, 6, and 3 overlapping pathways as well as 5, 3, and 4 overlapping process networks were observed between UA and JA, UA and JU, and JA and JU, respectively. Of these, three pathways and three process networks overlapped across the three groups. Interestingly, four representative pathways and six process networks were only noted in the JU group. Gene Ontology process analysis showed 2 processes were shared by all three treatment groups in the top 10 processes. The UA and JA combination resulted in synergistic effects through affecting multi-signal transduction pathways, different locations in the same pathway, and the new signaling pathway emerged in drug combination group, those together may enhance the treatment of cerebral ischemia-reperfusion injury through promoting neural cell apoptosis, decreasing calcium levels, inhibiting inflammation, and protecting neurons.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

BackgroundSalmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.ResultsAll selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.ConclusionsOur approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.