Background and Aims: Methionine has been proven to inhibit addictive behaviors of cocaine dependence. This study aimed to identify the potential mechanisms of MET relating to its inhibitory effects on cocaine induced cellular and behavioral changes. Methods: MRNA and miRNA high-throughput sequencing of the prefrontal cortex in a mouse model of cocaine conditioned place preference (CPP) combined with L-methionine was performed. Differentially expressed miRNAs (DE-miRNAs) and differentially expressed genes (DEGs) regulated by cocaine and inhibited by L-methionine were identified. DEGs were mapped to STRING database to construct a protein-protein interaction (PPI) network. Then, the identified DEGs were subjected to the DAVID webserver for functional annotation. Finally, miRNA-mRNA regulatory network and miRNA-mRNA-TF regulatory networks were established to screen key DE-miRNAs and coregulation network in Cytoscape. Results: Sequencing data analysis showed that L-methionine reversely regulated genes and miRNAs affected by cocaine. Pathways associated with drug addiction only enriched in CS-down with MC-up genes targeted by DE-miRNAs including GABAergic synapse, Glutamatergic synapse, Circadian entrainment, Axon guidance and Calcium signaling pathway. Drug addiction associated network was formed of 22 DEGs including calcium channel (Cacna1c, Cacna1e, Cacna1g and Cacng8), ephrin receptor genes (Ephb6 and Epha8) and ryanodine receptor genes (Ryr1 and Ryr2). Calcium channel gene network were identified as a core gene network modulated by L-methionine in response to cocaine dependence. Moreover, it was predicted that Grin1 and Fosb presented in TF-miRNA-mRNA coregulation network with a high degree of interaction as hub genes and interacted calcium channels. Conclusion: These identified key genes, miRNA and coregulation network demonstrated the efficacy of L-methionine in counteracting the effects of cocaine CPP. To a certain degree, it may provide some hints to better understand the underlying mechanism on L-methionine in response to cocaine abuse.

Radish (Raphanus sativus L.) is one of the most important vegetable crops worldwide. Taproot thickening represents a critical developmental period that determines yield and quality in radish life cycle. To isolate differentially expressed genes (DGEs) involved in radish taproot thickening process and explore the molecular mechanism underlying taproot development, three cDNA libraries from radish taproot collected at pre-cortex splitting stage (L1), cortex splitting stage (L2), and expanding stage (L3) were constructed and sequenced by RNA-Seq technology. More than seven million clean reads were obtained from the three libraries, from which 4,717,617 (L1, 65.35%), 4,809,588 (L2, 68.24%) and 4,973,745 (L3, 69.45%) reads were matched to the radish reference genes, respectively. A total of 85,939 transcripts were generated from three libraries, from which 10,450, 12,325, and 7392 differentially expressed transcripts (DETs) were detected in L1 vs. L2, L1 vs. L3, and L2 vs. L3 comparisons, respectively. Gene Ontology and pathway analysis showed that many DEGs, including EXPA9, Cyclin, CaM, Syntaxin, MADS-box, SAUR, and CalS were involved in cell events, cell wall modification, regulation of plant hormone levels, signal transduction and metabolisms, which may relate to taproot thickening. Furthermore, the integrated analysis of mRNA-miRNA revealed that 43 miRNAs and 92 genes formed 114 miRNA-target mRNA pairs were co-expressed, and three miRNA-target regulatory networks of taproot were constructed from different libraries. Finally, the expression patterns of 16 selected genes were confirmed using RT-qPCR analysis. A hypothetical model of genetic regulatory network associated with taproot thickening in radish was put forward. The taproot formation of radish is mainly attributed to cell differentiation, division and expansion, which are regulated and promoted by certain specific signal transduction pathways and metabolism processes. These results could provide new insights into the complex molecular mechanism underlying taproot thickening and facilitate genetic improvement of taproot in radish.

Strawberry is often subjected to cold stress in temperate regions when insulation measures are not strictly applied in protected cultivation. Cold stress adversely influences plant growth and development by triggering a massive change to the transcriptome. To provide the potential strategies in improving strawberry cold tolerance and give a glimpse into the understanding of the complex cold signaling pathways in plants, this study identified attractive candidate genes and revealed diverse regulatory networks that responded to cold stress in strawberry (Fragaria×ananassa) by a transcriptomic analysis. Totally, there were 2397 differentially expressed genes (DEGs) under cold stress treatment (T1) vs. normal treatment (CK). Of these, 1180 DEGs were upregulated, while 1217 DEGs were downregulated. Functional enrichment analysis showed that DEGs were significantly (adjusted P value < 0.05) overrepresented in six pathways including plant hormone signal transduction, flavonoid biosynthesis, mitogen-activated protein kinase (MAPK) signaling, starch and sucrose metabolism, circadian rhythm, and alpha-linolenic acid metabolism. The cold signaling initiated expression of downstream cold-responsive (COR) genes with cis-acting element ABRE or CRT/DRE in the ABA-independent or ABA-dependent pathway to impel plant defense against the stress. Strikingly, GIGANTEA (gene id 101308922), two-component response regulator-like PRR95 (gene id 101295449), and ethylene-responsive transcription factor ERF105-like (gene id 101295082) were dramatically induced under low-temperature treatment, indicating that they played an important role in response to cold stress in strawberry.

Methionine has been proven to inhibit addictive behaviors of cocaine dependence. However, the mechanism of methionine response to cocaine CPP is unknown. Recent evidence highlights piRNAs to regulate genes via a miRNA-like mechanism. Here, next-generation sequencing is used to study mechanism on methionine response to drug-induced behaviors though piRNA.

MicroRNAs (miRNAs) are small noncoding RNAs that play pivotal roles in plant growth, development and stress response. Chromium (Cr) is one of common environmental contaminants possessing potential health hazards to living organisms. To date, little is known about the regulatory roles of miRNAs in response to Cr stress in radish. To systematically identify Cr-responsive miRNAs and their targets in radish, two sRNA libraries derived from Cr-free (CK) and Cr-treated (Cr200) roots were constructed. With Solexa sequencing, 81 known and 72 novel miRNAs were identified, from which 54 known and 16 novel miRNAs were significantly differentially expressed under Cr stress. Several target genes for Cr-responsive miRNAs encode different transcription factor (TF) families, including SPLs, MYBs, ERFs and bZIPs, might regulate corresponding HM-related transcriptional processes in plants. Notably, a few key responsive enzymes or proteins, including HMA, YSL1 and ABC transporter protein were involved in Cr uptake and homeostasis process. Furthermore, the expression patterns of some Cr-responsive miRNAs and their targets were validated by RT-qPCR. This study represents the first characterization of Cr-responsive miRNAs and their targets in radish. The outcomes of this study could provide novel insights into miRNA-mediated regulatory mechanisms underlying plant response to Cr stress in root vegetable crops.

MicroRNAs (miRNAs) play vital regulatory roles in plant growth and development. The phase transition from vegetative growth to flowering is crucial in the life cycle of plants. To date, miRNA-mediated flowering regulatory networks remain largely unexplored in radish. In this study, two small RNA libraries from radish leaves at vegetative and reproductive stages were constructed and sequenced by Solexa sequencing. A total of 94 known miRNAs representing 21 conserved and 13 non-conserved miRNA families, and 44 potential novel miRNAs, were identified from the two libraries. In addition, 42 known and 17 novel miRNAs were significantly differentially expressed and identified as bolting-related miRNAs. RT-qPCR analysis revealed that some miRNAs exhibited tissue- or developmental stage-specific expression patterns. Moreover, 154 target transcripts were identified for 50 bolting-related miRNAs, which were predominately involved in plant development, signal transduction and transcriptional regulation. Based on the characterization of bolting-related miRNAs and their target genes, a putative schematic model of miRNA-mediated bolting and flowering regulatory network was proposed. These results could provide insights into bolting and flowering regulatory networks in radish, and facilitate dissecting the molecular mechanisms underlying bolting and flowering time regulation in vegetable crops.

Thyroid-carcinoma (THCA) is the most common malignancy with an increasing incidence. Recent evidence has emphasized the role of microRNA (miRNA) in THCA. However, knowledge concerning the roles of miRNAs in THCA is still limited. We therefore use a miRNA-target gene differential regulatory network (MGDRN) to identify key miRNAs and characterize their synergistic regulation in THCA. Both miRNA-target gene interactions from multiple databases and negative expression correlations between miRNA-target genes were used to characterize the interactions. Then, two regulatory networks involving normal and tumor conditions were constructed, respectively. The MGDRN was finally constructed using different interactions between the above two regulatory networks. By analyzing topological features of the MGDRN, four miRNAs (hsa-mir-152-3p, hsa-mir-148a, hsa-mir-130b and hsa-mir-15b) are identified as key miRNAs in THCA. Over-expression of mir-152-3p inhibited proliferation and colony formation of TPC-1 cells. Furthermore, mir-152-3p negatively regulated ERBB3 by binding to the 3'-UTR of ERBB3, and down-regulation of ERBB3 by small interfering (si)RNAs inhibited proliferation and colony formation of TPC-1 cells, indicating that mir-152-3p acted as an anti-tumor miRNA by negatively regulating ERBB3. Finally, two synergistically dysregulated modules were identified which may contribute to the initiation and progression of THCA. Overall, the results provided a better understanding of the molecular basis of THCA, and suggested novel treatment strategies for this cancer.

Whole-exome and whole-genome sequencing have revealed millions of somatic mutations associated with different human cancers, and the vast majority of them are located outside of coding sequences, making it challenging to directly interpret their functional effects. With the rapid advances in high-throughput sequencing technologies, genome-scale long-range chromatin interactions were detected, and distal target genes of regulatory elements were determined using three-dimensional (3D) chromatin looping. Herein, we present OncoBase (http://www.oncobase.biols.ac.cn/), an integrated database for annotating 81 385 242 somatic mutations in 68 cancer types from more than 120 cancer projects by exploring their roles in distal interactions between target genes and regulatory elements. OncoBase integrates local chromatin signatures, 3D chromatin interactions in different cell types and reconstruction of enhancer-target networks using state-of-the-art algorithms. It employs informative visualization tools to display the integrated local and 3D chromatin signatures and effects of somatic mutations on regulatory elements. Enhancer-promoter interactions estimated from chromatin interactions are integrated into a network diffusion system that quantitatively prioritizes somatic mutations and target genes from a large pool. Thus, OncoBase is a useful resource for the functional annotation of regulatory noncoding regions and systematically benchmarking the regulatory effects of embedded noncoding somatic mutations in human carcinogenesis.

Preeclampsia (PE) confers a significant risk for subsequent diagnosis with autism spectrum disorder (ASD), with the mechanisms underlying this observation being largely unknown. To identify molecular networks affected by both PE and ASD, we conducted a large-scale literature data mining and a gene set enrichment analysis (GSEA), followed by an expression mega-analysis in 13 independently profiled ASD datasets. Sets of genes implicated in ASD and in PE significantly overlap (156 common genes; p = 3.14E-67), with many biological pathways shared (94 pathways; p < 1.00E-21). A set of PE-driven molecular triggers possibly contributing to worsening the risk of subsequent ASD was identified, possibly representing a regulatory shift toward greater vulnerability to the development of ASD. Mega-analysis of expression highlighted RPS4Y1, an inhibitor of STAT3 that is expressed in a sexually dimorphic manner, as a contributor to both PE and ASD, which should be evaluated as a possible contributor to male predominance in ASD. A set of PE-driven molecular triggers may shift the developing brain toward a greater risk of ASD. One of these triggers, chromosome Y encoded gene RPS4Y1, an inhibitor of STAT3 signaling, warrants evaluation as a possible contributor to male predominance in ASD.

Broodiness in laying hens results in atrophy of the ovary and consequently decreases productivity. However, the regulatory mechanisms that drive ovary development remain elusive. Thus, we collected atrophic ovaries (AO) from 380-day-old broody chickens (BC) and normal ovaries (NO) from even-aged egg-laying hens (EH) for RNA sequencing. We identified 3,480 protein-coding transcripts that were differentially expressed (DE), including 1,719 that were down-regulated and 1,761 that were up-regulated in AO. There were 959 lncRNA transcripts that were DE, including 56 that were down-regulated and 903 that were up-regulated. Among the116 miRNAs that were DE, 79 were down-regulated and 37 were up-regulated in AO. Numerous DE protein-coding transcripts and target genes for miRNAs/lncRNAs were significantly enriched in reproductive processes, cell proliferation, and apoptosis pathways. A miRNA-intersection gene-pathway network was constructed by considering target relationships and correlation of the expression levels between ovary development-related genes and miRNAs. We also constructed a competing endogenous RNA (ceRNA) network by integrating competing relationships between protein-coding genes and lncRNA transcripts, and identified several lncRNA transcripts predicted to regulate the CASP6, CYP1B1, GADD45, MMP2, and SMAS2 genes. In conclusion, we discovered protein-coding genes, miRNAs, and lncRNA transcripts that are candidate regulators of ovary development in broody chickens.

We attempted to investigate the mechanism and susceptibility genes for diabetes in Han and Kazak ethnic individuals.The abdominal omental adipose tissues were obtained from diabetic cases and healthy controls. The gene expression and methylation data were produced for Kazak and Han individuals, respectively, and analyzed by bioinformatics analysis.We obtained 921 differentially expressed genes (DEGs) in Han group and 1772 in Kazak group. DEGs in Han group were significantly related with type 2 diabetes mellitus, and biosynthesis of amino acids, while the DEGs specific to Kazak patients were significantly enriched in metabolism-related pathways such as carbon metabolism, propanoate metabolism, and 2-oxocarboxylic acid metabolism. Major facilitator superfamily domain containing 1 (MFSD1) was found to be a methylation associated gene at hypermethylation site of cg16289538 in Han group. Rho guanine nucleotide exchange factor 1 (ARHGEF1) was the susceptible gene corresponding to the methylation sites of cg18800192 and cg00759295 in Kazak group. ARHGEF1 was also a node in protein-protein interaction network and significantly enriched in hsa04270: vascular smooth muscle contraction pathways.The molecular mechanism of diabetes may be different in Han and Kazak patients. MFSD1 and ARHGEF1 may be the diabetes susceptible genes.

Radish (Raphanus sativus L.) is an important worldwide root vegetable crop with high nutrient values and is adversely affected by non-essential heavy metals including chromium (Cr). Little is known about the molecular mechanism underlying Cr stress response in radish. In this study, RNA-Seq technique was employed to identify differentially expressed genes (DEGs) under Cr stress. Based on de novo transcriptome assembly, there were 30,676 unigenes representing 60,881 transcripts isolated from radish root under Cr stress. Differential gene analysis revealed that 2985 uingenes were significantly differentially expressed between Cr-free (CK) and Cr-treated (Cr600) libraries, among which 1424 were up-regulated and 1561 down-regulated. Gene ontology (GO) analysis revealed that these DEGs were mainly involved in primary metabolic process, response to abiotic stimulus, cellular metabolic process and small molecule metabolic process. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis showed that the DEGs were mainly involved in protein processing in endoplasmic reticulum, starch and sucrose metabolism, amino acid metabolism, glutathione metabolism, drug and xenobiotics by cytochrome P450 metabolism. RT-qPCR analysis showed that the expression patterns of 12 randomly selected DEGs were highly accordant with the results from RNA-seq. Furthermore, many candidate genes including signaling protein kinases, transcription factors and metal transporters, chelate compound biosynthesis and antioxidant system, were involved in defense and detoxification mechanisms of Cr stress response regulatory networks. These results would provide novel insight into molecular mechanism underlying plant responsiveness to Cr stress and facilitate further genetic manipulation on Cr uptake and accumulation in radish.

Rheumatoid arthritis (RA) is a chronic autoimmune disease associated with an increased risk of death, but its underlying mechanisms are not fully understood. Circular RNAs (circRNAs) have recently been implicated in various biological processes. The aim of this study was to investigate the key circRNAs related to RA.

Mapping expression quantitative trait loci (eQTL) of targeted genes represents a powerful and widely adopted approach to identify putative regulatory variants. Linking regulation differences to specific genes might assist in the identification of networks and interactions. The objective of this study is to identify eQTL underlying expression of four gene families encoding isoflavone synthetic enzymes involved in the phenylpropanoid pathway, which are phenylalanine ammonia-lyase (PAL; EC 4.3.1.5), chalcone synthase (CHS; EC 2.3.1.74), 2-hydroxyisoflavanone synthase (IFS; EC1.14.13.136) and flavanone 3-hydroxylase (F3H; EC 1.14.11.9). A population of 130 recombinant inbred lines (F5:11), derived from a cross between soybean cultivar 'Zhongdou 27' (high isoflavone) and 'Jiunong 20' (low isoflavone), and a total of 194 simple sequence repeat (SSR) markers were used in this study. Overlapped loci of eQTLs and phenotypic QTLs (pQTLs) were analyzed to identify the potential candidate genes underlying the accumulation of isoflavone in soybean seed.

The present study is to screen lymph node metastasis-related microRNAs (miRNAs) in lung adenocarcinoma (LUAD) and uncover their underlying mechanisms.

The fungus Aspergillus flavus (A. flavus) is a serious threat to maize (Zea mays) production worldwide. It causes considerable yield and economic losses, and poses a health risk to humans and livestock due to the high toxicity of aflatoxin. However, key genes and regulatory networks conferring maize resistance to A. flavus are not clear, especially at the early stage of infection. Here, we performed a comprehensive transcriptome analysis of two maize inbred lines with contrasting resistance to A. flavus infection.

The appropriate timing of bolting and flowering is pivotal for reproductive success in Brassicaceae crops including radish (Raphanus sativus L.). Although several flowering regulatory pathways had been described in some plant species, no study on genetic networks of bolting and flowering regulation was performed in radish. In this study, to generate dataset of radish unigene sequences for large-scale gene discovery and functional pathway identification, a cDNA library from mixed radish leaves at different developmental stages was subjected to high-throughput RNA sequencing (RNA-seq).

MicroRNAs (miRNAs) are endogenous non-coding small RNAs that play vital regulatory roles in plant growth, development, and environmental stress responses. Cadmium (Cd) is a non-essential heavy metal that is highly toxic to living organisms. To date, a number of conserved and non-conserved miRNAs have been identified to be involved in response to Cd stress in some plant species. However, the miRNA-mediated gene regulatory networks responsive to Cd stress in radish (Raphanus sativus L.) remain largely unexplored. To dissect Cd-responsive miRNAs and their targets systematically at the global level, two small RNA libraries were constructed from Cd-treated and Cd-free roots of radish seedlings. Using Solexa sequencing technology, 93 conserved and 16 non-conserved miRNAs (representing 26 miRNA families) and 28 novel miRNAs (representing 22 miRNA families) were identified. In all, 15 known and eight novel miRNA families were significantly differently regulated under Cd stress. The expression patterns of a set of Cd-responsive miRNAs were validated by quantitative real-time PCR. Based on the radish mRNA transcriptome, 18 and 71 targets for novel and known miRNA families, respectively, were identified by the degradome sequencing approach. Furthermore, a few target transcripts including phytochelatin synthase 1 (PCS1), iron transporter protein, and ABC transporter protein were involved in plant response to Cd stress. This study represents the first transcriptome-based analysis of miRNAs and their targets responsive to Cd stress in radish roots. These findings could provide valuable information for functional characterization of miRNAs and their targets in regulatory networks responsive to Cd stress in radish.

Vibrio coralliilyticus is known to cause coral diseases, especially under environmental perturbation, but its impact on coral physiology and underpinning mechanism is poorly understood. In the present study, we investigated cytological, immunological, and metatranscriptomic responses of the scleractinian coral Pocillopora damicornis to V. coralliilyticus infection. The density and chlorophyll content of symbiotic zooxanthellae decreased significantly at 12 and 24 h after Vibrio challenge. The activities of antioxidant enzymes such as superoxide dismutase and catalase, nitric oxide synthase, phenoloxidase (PO), and the activation level of caspase3 all rose significantly in P. damicornis after Vibrio challenge. In the metatranscriptomic analysis, we found 10 significantly upregulated genes in the symbionts at 24 h after the challenge, which were mostly involved in the metabolism of nucleic acid and polysaccharide, and 133 significantly down-regulated symbiont genes, which were mainly related to amino acid catabolism and transport. Meanwhile, 1432 significantly upregulated coral genes were revealed, highly overrepresented in GO terms that are mostly related to the regulation of immune response, the regulation of cytokine production, and innate immune response. Furthermore, at 24 h after Vibrio challenge, 890 coral genes were significantly downregulated, highly overrepresented in four GO terms implicated in defense response. These results in concert suggest that V. coralliilyticus infection triggered the innate immune response including the redox, PO, and apoptosis systems, but repressed the response of the complement system in the scleractinian coral P. damicornis, accompanied by symbiont density decrease and symbiosis collapse through disordering the metabolism of the symbionts. These findings shed light on the molecular regulatory processes underlying bleaching and degradation of P. damicornis resulting from the infection of V. coralliilyticus.

Solute carrier family 38 (SLC38s) transporters play important roles in amino acid transportation and signaling transduction. However, their genetic alterations and biological roles in tumors are still largely unclear. This study aimed to elucidate the genetic signatures of SLC38s transporters and their implications in esophageal squamous cell carcinoma (ESCC).