RAR-related orphan receptors A (RORA) and B (RORB) and voltage-gated sodium channel type 1 (SCN1A) genes play critical roles in the regulation of the circadian clock. Evidence has shown an association of RORA and RORB polymorphisms with susceptibility to autism and depression. Hence, we tested the association of RORA rs12912233, rs16943429, rs880626, rs2290430, and rs12900948; RORB rs1157358, rs7022435, rs3750420, and rs3903529; and SCN1A rs3812718 with epilepsy risk in the Malaysians. DNA was genotyped in 1789 subjects (39% epilepsy patients) by using MassARRAY (Sequenom). Significant association was obtained for rs12912233 in Malaysian Chinese (p=0.003). Interaction between rs12912233-rs880626 and rs3812718 was associated with the epilepsy risk in the subjects overall (p=0.001). Results show that RORA rs12912233 alone might be a possible risk variant for epilepsy in Malaysian Chinese, but that, together with RORA rs880626 and SCN1A rs3812718, this polymorphism may have a synergistic effect in the epilepsy risk in Malaysians.
Cocoa polyphenol (CP), due to their biological actions, may be supplementary treatments for adipose tissue-fat gain. However, the molecular mechanism of CPs is still ambiguous. This study investigated the hypothesis that CP treatment modulates expressing of lipid metabolism genes in mesenteric white adipose tissue (MES-WAT). Sprague-Dawley (SD) rats were fed a low-fat (LF) or high-fat (HF) diet for 12 weeks. Thereafter, HFD rats (n = 10/group) were treated at a dose of 600 mg/kg bw/day CPs (HFD + CPs) for 4 weeks. DNA microarray analysis resulted in 753 genes of the 13,008 genes expressed. Bioinformatics tools showed CP treatment significantly decreased gene expression levels for lipogenic enzymes, while increased the mRNA levels responsible for lipolysis enzymes. CP administration differentially regulates gene expression involved in lipid metabolism in MES-WAT. These data unveil a new insight into the molecular mechanisms underlying the pharmacological effect of CPs on obesity biomarkers in obese rats.
Long non-coding RNAs (lncRNAs) have been shown to possess a wide range of functions in both cellular and developmental processes including cancers. Although some of the lncRNAs have been implicated in the regulation of the immune response, the exact function of the large majority of lncRNAs still remains unknown. In this study, we characterized the lncRNAs in human primary monocytes, an essential component of the innate immune system. We performed RNA sequencing of monocytes from four individuals and combined our data with eleven other publicly available datasets. Our analysis led to identification of ~8000 lncRNAs of which >1000 have not been previously reported in monocytes. PCR-based validation of a subset of the identified novel long intergenic noncoding RNAs (lincRNAs) revealed distinct expression patterns. Our study provides a landscape of lncRNAs in monocytes, which could facilitate future experimental studies to characterize the functions of these molecules in the innate immune system.
Natural antisense transcripts (NATs) are involved in cellular development and regulatory processes. Multiple NATs at the Sox4 gene locus are spatiotemporally regulated throughout murine cerebral corticogenesis. In the study, we evaluated the potential functional role of Sox4 NATs at Sox4 gene locus. We demonstrated Sox4 sense and NATs formed dsRNA aggregates in the cytoplasm of brain cells. Over expression of Sox4 NATs in NIH/3T3 cells generally did not alter the level of Sox4 mRNA expression or protein translation. Upregulation of a Sox4 NAT known as Sox4ot1 led to the production of a novel small RNA, Sox4_sir3. Its biogenesis is Dicer1-dependent and has characteristics resemble piRNA. Expression of Sox4_sir3 was observed in the marginal and germinative zones of the developing and postnatal brains suggesting a potential role in regulating neurogenesis. We proposed that Sox4 sense-NATs serve as Dicer1-dependent templates to produce a novel endo-siRNA- or piRNA-like Sox4_sir3.
Red algae are a major source of marine sulfated galactans. In this study, orthologs and inparalogs from seven red algae were analyzed and compared with the aim to discover differences in algal galactan biosynthesis and related pathways of these algae. Red algal orthologs for putative carbohydrate sulfotransferases were found to be prevalent in Porphyridium purpureum, Florideophytes and Bangiophytes, while red algal orthologs for putative chondroitin sulfate synthases, sulfurylases, and porphyranases /carrageenases were found exclusively in Florideophytes and Bangiophytes. The acquirement of these genes could have happened after the divergence from Cyanidiales red algae. Cyanidiales red algae were found to have more number and types of putative sulfate permeases, suggesting that these genes could have been acquired in adaptation to the environmental stresses and biogeochemistry of respective habitats. The findings of this study shed lights on the evolution of different homeostasis mechanisms by the early and late diverging red algal orders.
Pieris rapae is a serious pest of brassicas worldwide. We performed de novo assembly of P. rapae transcriptome by next-generation sequencing and assembled approximately 65,727,422 clean paired-end reads into 32,118 unigenes, of which 13,585 were mapped to 255 pathways in the KEGG database. A total of 6173 novel transcripts were identified from reads directly mapped to P. rapae genome. Additionally, 1490 SSRs, 301,377 SNPs, and 29,284 InDels were identified as potential molecular markers to explore polymorphism within P. rapae populations. We screened and mapped 36 transcripts related to OBP, CSP, SNMP, PBAN, and OR. We analyzed the expression profiles of 7 selected genes involved in pheromone transport and degradation by quantitative real-time PCR; these genes are sex-specific and differentially expressed in the developmental stages. Overall, the comprehensive transcriptome resources described in this study could help understand and identify molecular targets particularly reproduction-related genes for developing effective P. rapae management tools.
The genus Meridianimaribacter is one of the least-studied genera within Cytophaga-Flavobacteria. To date, no genomic analysis of Meridianimaribacter has been reported. In this study, Meridianimaribacter sp. strain CL38, a lignocellulose degrading halophile was isolated from mangrove soil. The genome of strain CL38 was sequenced and analyzed. The assembled genome contains 17 contigs with 3.33 Mbp, a GC content of 33.13% and a total of 2982 genes predicted. Lignocellulose degrading enzymes such as cellulases (GH3, 5, 9, 16, 74 and 144), xylanases (GH43 and CE4) and mannanases (GH5, 26, 27 and 130) are encoded in the genome. Furthermore, strain CL38 demonstrated its ability to decompose empty fruit bunch, a lignocellulosic waste residue arising from palm oil industry. The genome information coupled with experimental studies confirmed the ability of strain CL38 to degrade lignocellulosic biomass. Therefore, Meridianimaribacter sp. strain CL38, with its halotolerance, could be useful for seawater based lignocellulosic biorefining.
The free-living flatworm Macrostoma lignano (M. lignano) is an emerging model organism for aging and regeneration research. Long intergenic non-coding RNAs (lincRNAs) have important roles in many biological processes such as aging, stem cell maintenance and differentiation. However, to date, there is no systematic identification of lincRNAs in M. lignano. By using public RNA-seq data, we identified a total of 2547 lincRNA transcripts in M. lignano genome. We discovered that M. lignano lincRNAs shared many characteristics with other species such as shorter in length, lower GC content, and lower in expression compared to protein-coding genes. Unlike protein-coding genes, M. lignano lincRNAs showed higher tendency to be expressed in temporal and region-specific fashion. Additionally, co-expression network analysis and functional enrichment suggest that M. lignano lincRNAs have potential roles in regeneration. This study will provide important resources and pave the way for investigations on non-coding genes involved in aging and regeneration.
In this study, we first identified male-specific SNP markers using restriction site-associated DNA sequencing, and further developed a PCR-based sex identification technique for Charybdis feriatus. A total of 296.96 million clean reads were obtained, with 114.95 and 182.01 million from females and males. After assembly and alignment, 10 SNP markers were identified being heterozygous in males but homozygous in females. Five markers were further confirmed to be male-specific in a large number of individuals. Moreover, two male-specific sense primers and a common antisense primer were designed, using which, a PCR-based genetic sex identification method was successfully developed and used to identify the sex of 103 individuals, with a result of 49 females and 54 males. The presence of male-specific SNP markers suggests an XX/XY sex determination system for C. feriatus. These findings should be helpful for better understanding sex determination mechanism, and drafting artificial breeding program in crustaceans.
Differences in expression of potential virulence and survival genes were associated with B. pseudomallei colony morphology variants. Microarray was used to investigate B. pseudomallei transcriptome alterations among the wild type and small colony variant (SCV) pre- and post-exposed to A549 cells. SCV pre- and post-exposed have lower metabolic requirements and consume lesser energy than the wild type pre- and post-exposed to A549. However, both the wild type and SCV limit their metabolic activities post- infection of A549 cells and this is indicated by the down-regulation of genes implicated in the metabolism of amino acids, carbohydrate, lipid, and other amino acids. Many well-known virulence and survival factors, including T3SS, fimbriae, capsular polysaccharides and stress response were up-regulated in both the wild type and SCV pre- and post-exposed to A549 cells. Microarray analysis demonstrated essential differences in bacterial response associated with virulence and survival pre- and post-exposed to A549 cells.
Infection by the rhizocephalan parasite Sacculina beauforti can have detrimental effects on mud crab Scylla olivacea. However, the molecular changes that occur during rhizocephalan infection are poorly understood. Due to the disruption in the reproductive system after infection, the gonadal transcriptomic profiles of non-infected and infected Scylla olivacea were compared. A total of 686 and 843 unigenes were differentially expressed between non-infected and infected males, and females, respectively. The number of DEGs increased after infection. By comparing shared DEGs of non-infected and infected individuals, potential immune- and reproduction-related of host, and immune- and metabolism-related genes of parasite are highlighted. The only shared KEGG pathway between non-infected and infected individuals was the ribosome pathway. In summary, findings in this study provide new insights into the host-parasite relationship of rhizocephalan parasites and their crustacean hosts.
The role of microbiota in gut-brain communication has led to the development of probiotics promoting brain health. Here we report a genomic study of a Lactobacillus fermentum PS150 and its patented bioactive protein, elongation factor Tu (EF-Tu), which is associated with cognitive improvement in rats. The L. fermentum PS150 circular chromosome is 2,238,401 bp and it consists of 2281 genes. Chromosome comparisons with other L. fermentum strains highlighted a cluster of glycosyltransferases as potential candidate probiotic factors besides EF-Tu. Molecular evolutionary analyses on EF-Tu genes (tuf) in 235 bacteria species revealed one to three copies of the gene per genome. Seven tuf pseudogenes were found and three species only possessed pseudogenes, which is an unprecedented finding. Protein variability analysis of EF-Tu showed five highly variable residues (40 K, 41G, 42 L, 44 K, and 46E) on the protein surface, which warrant further investigation regarding their potential roles as binding sites.
Differential gene and transcript expression pattern of human primary monocytes from healthy young subjects were profiled under different sequencing depths (50M, 100M, and 200M reads). The raw data consisted of 1.3 billion reads generated from RNA sequencing (RNA-Seq) experiments. A total of 17,657 genes and 75,392 transcripts were obtained at sequencing depth of 200M. Total splice junction reads showed an even more significant increase. Comparative analysis of the expression patterns of immune-related genes revealed a total of 217 differentially expressed (DE) protein-coding genes and 50 DE novel transcripts, in which 40 DE protein-coding genes were related to the immune system. At higher sequencing depth, more genes, known and novel transcripts were identified and larger proportion of reads were allowed to map across splice junctions. The results also showed that increase in sequencing depth has no effect on the sequence alignment.
High-depth next generation sequencing data provide valuable insights into the number and distribution of RNA editing events. Here, we report the RNA editing events at cellular level of human primary monocyte using high-depth whole genomic and transcriptomic sequencing data. We identified over a ten thousand putative RNA editing sites and 69% of the sites were A-to-I editing sites. The sites enriched in repetitive sequences and intronic regions. High-depth sequencing datasets revealed that 90% of the canonical sites were edited at lower frequencies (<0.7). Single and multiple human monocytes and brain tissues samples were analyzed through genome sequence independent approach. The later approach was observed to identify more editing sites. Monocytes was observed to contain more C-to-U editing sites compared to brain tissues. Our results establish comparable pipeline that can address current limitations as well as demonstrate the potential for highly sensitive detection of RNA editing events in single cell type.
Gene expression data are expected to make a great contribution in the producing of efficient cancer diagnosis and prognosis. Gene expression data are coded by large measured genes, and only of a few number of them carry precious information for different classes of samples. Recently, several researchers proposed gene selection methods based on metaheuristic algorithms for analysing and interpreting gene expression data. However, due to large number of selected genes with limited number of patient's samples and complex interaction between genes, many gene selection methods experienced challenges in order to approach the most relevant and reliable genes. Hence, in this paper, a hybrid filter/wrapper, called rMRMR-MBA is proposed for gene selection problem. In this method, robust Minimum Redundancy Maximum Relevancy (rMRMR) as filter to select the most promising genes and an modified bat algorithm (MBA) as search engine in wrapper approach is proposed to identify a small set of informative genes. The performance of the proposed method has been evaluated using ten gene expression datasets. For performance evaluation, MBA is evaluated by studying the convergence behaviour of MBA with and without TRIZ optimisation operators. For comparative evaluation, the results of the proposed rMRMR-MBA were compared against ten state-of-arts methods using the same datasets. The comparative study demonstrates that the proposed method produced better results in terms of classification accuracy and number of selected genes in two out of ten datasets and competitive results on the remaining datasets. In a nutshell, the proposed method is able to produce very promising results with high classification accuracy which can be considered a promising contribution for gene selection domain.
Exposing the skin to solar UV radiation induces cascades of signaling pathways and biological alterations such as redox imbalance, suppression of antioxidant genes and programmed cell death. Therefore, the aim of this study was to use RNA-Seq to unravel the effects of UV radiation on Normal Human Adult Fibroblast cells (NHDF). Cells were exposed to UV (20 mJ/cm2 for 3 mins) and incubated for 24 h. Total mRNA from the cells generated libraries of 72,080,648 and 40,750,939 raw reads from UV-treated and control cells respectively. Of the differentially expressed genes (DEGs) produced 2,007 were up-regulated and 2,791 were down-regulated (fold change ≥2, p
PIWI-interacting RNAs (piRNAs) are abundantly found in germ cells and involved in gametogenesis and gonadal development. Information on the regulatory roles of piRNAs in crustacean reproduction, however, is scarce. Thus, we identified gonadal piRNAs of mud crab Scylla paramamosain. Of the 115,491 novel piRNAs, 596 were differentially expressed. Subsequently, 389,887 potential piRNA-target genes were predicted. The expression of 4 piRNAs and 9 genes with high piRNA interactions were validated with the inclusion of additional immature specimens, including LRP2 that is involved in growth and reproduction, MDN1 in ribosome biogenesis pathway and gametogenesis, and PRKDC, a DNA repair gene involved in gonadal differentiation and maturation. KEGG analysis further revealed the involvement of predicted piRNA target genes in gametogenesis- and reproduction-related pathways. Our findings provide baseline information of mud crab piRNAs and their differential expression between testes and ovaries suggests that piRNAs play an essential role in regulating gametogenesis and gonadal development.
Bacterial dormancy plays a crucial role in maintaining the functioning and diversity of microbial communities in natural environments. However, the metabolic regulations of the dormancy of bacteria in natural habitats, especially marine habitats, have remained largely unknown. A marine bacterium, Microbulbifer aggregans CCB-MM1 exhibits rod-to-coccus cell shape change during the dormant state. Therefore, to clarify the metabolic regulation of the dormancy, differential gene expression analysis based on RNA-Seq was performed between rod- (vegetative), intermediate, and coccus-shaped cells (dormancy). The RNA-Seq data revealed that one of two distinct electron transfer chains was upregulated in the dormancy. Dissimilatory sulfite reductase and soluble hydrogenase were also highly upregulated in the dormancy. In addition, induction of the dormancy of MM1 in the absence of MgSO4 was slower than that in the presence of MgSO4. These results indicate that the sulfate-reducing pathway plays an important role in entering the dormancy of MM1.
Agrobacterium tumefaciens is considered a prominent phytopathogen, though most isolates are nonpathogenic. Agrobacteria can inhabit plant tissues interacting with other microorganisms. Yeasts are likewise part of these communities. We analyzed the quorum sensing (QS) systems of A. tumefaciens strain 6N2, and its relevance for the interaction with the yeast Meyerozyma guilliermondii, both sugarcane endophytes. We show that strain 6N2 is nonpathogenic, produces OHC8-HSL, OHC10-HSL, OC12-HSL and OHC12-HSL as QS signals, and possesses a complex QS architecture, with one truncated, two complete systems, and three additional QS-signal receptors. A proteomic approach showed differences in QS-regulated proteins between pure (64 proteins) and dual (33 proteins) cultures. Seven proteins were consistently regulated by quorum sensing in pure and dual cultures. M. guilliermondii proteins influenced by QS activity were also evaluated. Several up- and down- regulated proteins differed depending on the bacterial QS. These results show the QS regulation in the bacteria-yeast interactions.
Agar and agarose have wide applications in food and pharmaceutical industries. Knowledge on the genome of red seaweeds that produce them is still lacking. To fill the gap in genome analyses of these red algae, we have sequenced the nuclear and organellar genomes of an agarophyte, Gracilaria changii. The partial nuclear genome sequence of G. changii has a total length of 35.8Mb with 10,912 predicted protein coding sequences. Only 39.4% predicted proteins were found to have significant matches to protein sequences in SwissProt. The chloroplast genome of G. changii is 183,855bp with a total of 201 open reading frames (ORFs), 29 tRNAs and 3 rRNAs predicted. Five genes: ssrA, leuC and leuD CP76_p173 (orf139) and pbsA were absent in the chloroplast genome of G. changii. The genome information is valuable in accelerating functional studies of individual genes and resolving evolutionary relationship of red seaweeds.