Heavy metal pollution is an increasing global concern. Among heavy metals, mercury (Hg) is especially dangerous because of its massive release into the environment and high toxicity, especially for aquatic organisms. The molecular response mechanisms of algae to Hg exposure are mostly unknown. Here, we combine physiological, biochemical, and transcriptomic analysis to provide, for the first time, a comprehensive view on the pathways activated in Chromera velia in response to toxic levels of Hg. Production of hydrogen peroxide and superoxide anion, two reactive oxygen species (ROS), showed opposite patterns in response to Hg²⁺ while reactive nitrogen species (RNS) levels did not change. A deep RNA sequencing analysis generated a total of 307,738,790 high-quality reads assembled in 122,874 transcripts, representing 89,853 unigenes successfully annotated in databases. Detailed analysis of the differently expressed genes corroborate the biochemical results observed in ROS production and suggests novel putative molecular mechanisms in the algal response to Hg²⁺. Moreover, we indicated that important transcription factor (TF) families associated with stress responses differentially expressed in C. velia cultures under Hg stress. Our study presents the first in-depth transcriptomic analysis of C. velia, focusing on the expression of genes involved in different detoxification defense systems in response to heavy metal stress.

Mucosal melanomas (MM) are rare aggressive cancers in humans and one of the most common forms of oral cancers in dogs. Similar biological and histological features are shared between MM in both species making dogs a powerful model for comparative oncology studies of melanomas. Although exome sequencing recently identified recurrent coding mutations in canine MM, little is known about changes in non-coding gene expression and more particularly in canine long non-coding RNAs (lncRNAs), which are commonly dysregulated in human cancers. Here, we sampled a large cohort (n= 52) of canine normal/tumor oral MM from three predisposed breeds (poodles, Labrador retrievers and golden retrievers) and used deep transcriptome sequencing to identify more than 400 differentially expressed (DE) lncRNAs. We further prioritized candidate lncRNAs by comparative genomic analysis to pinpoint 26 dog-human conserved DE lncRNAs, including SOX21-AS, ZEB2-AS and CASC15 lncRNAs. Using unsupervised co-expression networks analysis with coding genes, we inferred potential functions of DE lncRNAs suggesting associations with cancer-related genes, cell cycle and carbohydrate metabolism GO terms. Finally, we exploited our multi-breed design to identify DE lncRNAs per breed. This study provides a unique transcriptomic resource for studying oral melanoma in dogs and highlights lncRNAs that may potentially be diagnostic or therapeutic targets for human and veterinary medicine.

Objective: Rett syndrome (RTT) is an X-linked neurodevelopmental disorder caused by mutations in MeCP2, a transcription factor. MeCP2 mutations cause abnormal expression of downstream genes and eventually lead to brain dysfunction. The role of MeCP2 in brain neural development remains unclear. To further elucidate this role, a MeCP2-null rat model was created with the CRISPR/cas9 system. Method: A MeCP2-cas9 vector was constructed and then microinjected into fertilized rat ova in vitro. Two mutations by CRISPR/cas9 were confirmed to cause deletions in exon 2 of MeCP2 via DNA sequencing, and protein expression was measured by Western blotting. Transcriptome data for three brain tissues, the cerebellum, cerebral cortex and hippocampus, were obtained via next-generation sequencing. Results: MeCP2-null rats were successfully obtained, and preliminary analysis showed that the MeCP2-null rats exhibited motor dysfunction and anxious and depressed behaviour. In addition, differentially expressed genes (DEGs) were identified in the three MeCP2-null brain tissues compared to wild-type rat brain tissues. In the rat cerebellum, 388 downregulated DEGs were mainly involved in the calcium ion signalling pathway and the PI3K-Akt signalling pathway. In the cerebral cortex, 386 upregulated DEGs were primarily involved in intracellular signal transduction, protein phosphorylation and the MAPK signalling pathway. In the hippocampus, the DEGs were related to the MAPK signalling pathway. Conclusion: We constructed a MeCP2-null rat model with unique features with CRISPR/cas9 technology to study RTT and analysed DEGs in three rat brain tissues to highlight potential targets for the development of new medicines.

Plants grown under fluctuating light impact plant developments compared with those grown under non-fluctuating light conditions. However, our knowledge on the underlying regulatory mechanisms is still quite limited, particularly from the transcriptional perspective. In order to investigate the influence of different light intensity distributions on tomato plant development, we designed three fluctuating light intensity distributions with the non-fluctuating light intensity as control and compared the transcriptional differences after five weeks of treatment. We found plant height and aerial/root weight were significantly reduced under all fluctuating light treatments. Transcriptome analysis revealed that the number of up and down regulated genes had a distinct distribution pattern between different treatments and control. The largest difference between the numbers of down and up regulated genes was found between treatment 1 and 3, reaching to a total of 416 genes. The number and type of the top 20 enriched pathways differed between treatments and control. The largest number of genes enriched was involved in the biosynthesis of secondary metabolites. These results provide insights into the transcriptional regulations of tomato under different light intensity distributions.

The rapid expansion of transcriptomics from increased affordability of next-generation sequencing (NGS) technologies generates rocketing amounts of gene expression data across biology and medicine, and notably in cancer research. Concomitantly, many bioinformatics tools were developed to streamline gene expression analysis and quantification. We tested the concordance of NGS RNA sequencing (RNA-seq) analysis outcomes between the two predominant programs for reads alignment, HISAT2 and STAR, and the two most popular programs for quantifying gene expression in NGS experiments, edgeR and DESeq2, using RNA-seq data from a series of breast cancer progression specimens, which include histologically confirmed normal, early neoplasia, ductal carcinoma in situ and infiltrating ductal carcinoma samples microdissected from formalin fixed, paraffin embedded (FFPE) breast tissue blocks. We identified significant differences in aligners’ performance: HISAT2 was prone to misalign reads to retrogene genomic loci, STAR generated more precise alignments, especially for early neoplasia samples. edgeR and DESeq2 produced similar lists of differentially expressed genes in stage comparisons, with edgeR producing more conservative, though shorter, lists of genes. Albeit, Gene Ontology (GO) enrichment analysis revealed no skewness in significant GO categories identified among differentially expressed genes by edgeR vs DESeq2. As transcriptome analysis of archived FFPE samples becomes a vanguard of precision medicine, identification and fine-tuning of bioinformatics tools becomes critical for clinical research. Our results indicate that STAR and edgeR are well-suited tools for differential gene expression analysis from FFPE samples.

Combined awareness about the power and limitations of bioinformatics and molecular biology enables advanced research based on high-throughput data. Despite an increasing demand for scientists with a combined background in both fields, the education in dry lab and wet lab is often separated. This work describes an example of integrated education with focus on genomics and transcriptomics. Participants learn computational and molecular biology methods in the same practical course. Peer-review is applied as a teaching method to foster cooperative learning of students with heterogeneous backgrounds. Evaluation results indicate acceptance and appreciation of this approach.

Viroids are smallest pathogen that consist of non-capsidated, single-stranded non-coding RNA replicons and exploits host factors for their replication and propagation. The severe stunting disease caused by Citrus bark cracking viroid (CBCVd) is a serious threat, which spread rapidly within hop gardens. In this study, we employed comprehensive transcriptome analyses to dissect host-viroid interactions and identify gene expression changes associated with disease development in hop. Our analysis revealed that CBCVd-infection resulted in the massive modulation of activity of over 2000 genes. Expression of genes associated with plant immune responses (protein kinase and mitogen-activated protein kinase), hypersensitive responses, phytohormone signaling pathways, photosynthesis, pigment metabolism, protein metabolism, sugar metabolism and modification and others were altered, which could be attributed to systemic symptom development upon CBCVd-infection in hop. In addition, genes encoding RNA-dependent RNA polymerase, pathogenesis-related protein, chitinase as well as those related to basal defense responses were up-regulated. The expression levels of several genes identified from RNA sequencing analysis were confirmed by qRT-PCR. Our systematic comprehensive CBCVd-responsive transcriptome analysis provides a better understanding and insights into complex viroid-hop plant interaction. This information will assist further in the development of future measures for the prevention of CBCVd spread in hop fields.

The eyes are the sense organs through we view the world around us, and the cornea is the transparent layer which covers the outer, visible part of the eye. It is known that the gene expression of corneal endothelium depends on age. And the expressed genes of the endothelium of cornea for pediatric samples is different than that of the adult samples. The purpose or objective of this study was to characterize human corneal endothelial cell (HCEnC) gene expression and differential gene expression and to detect the expressed genes mapped to chromosomal loci associated with some corneal endothelial dystrophies. Both upregulated and downregulated genes were analyzed. For this purpose, total RNA was isolated from ex vivo corneal endothelium taken from six pediatric and five adult corneas. The complementary DNA was hybridized to the Affymetrix GeneChip 1.1 ST array. The data analysis was performed using the Enrichr software for both upregulated and downregulated genes. These are described in more detail in the results and discussion section. This study uses bioinformatics tools to identify and analyse gene sets present in the transcriptome of the corneal endothelium, and tries to find out and observe the relation of the aging effect on the corneal endothelium gene expression. The human subjects had participated voluntarily and informed consent was obtained from all before carrying out any testing procedure. Proper guidelines from the hospital ethical committee were also followed and no harmful chemicals were used on the participants. This study simply aims to raise some awareness of the given topic among the local people so that they are better able to take informed choices about their health in the near future, and also so that they seek medical help when necessary and have no inhibitions in doing so.

To understand the molecular mechanism associated with growth variability in bivalves, the Solexa/Illumina technology was employed to analyze the transcriptomic profiles of extreme growth rate differences (fast- VS. slow-growing individuals) in one full-sib family of the ark shell Scapharca subcrenata. De novo assembly of S. subcrenata transcriptome yielded 276,082,016 raw reads, which were assembled into 98,502 unique transcripts by Trinity strategy. A total of 6,357 differentially expressed genes (DEGs) were obtained between fast- and slow-growing individuals, with 580 up-regulated expression and 5777 down-regulated expression. Functional annotation revealed that the largest proportion of DEGs were classified to the large or small subunit ribosomal protein, all of which showed significantly lower expression levels in fast-growing group than those in slow-growing group. GO enrichment analysis identified the maximum of DEGs to biological process, followed by molecular function and cellular component. Most of the top enriched KEGG pathways were related to energy metabolism, protein synthesis and degradation. These findings reveal the link between gene expression and contrasting phenotypes in ark shells, which support that fast-growing individuals may be resulted from decreased energy requirements for metabolism maintenance, accompanying with greater efficiency of protein synthesis and degradation in bivalves.

The early generative phase of cauliflower (Brassica oleracea var. botrytis) curd ripening is sensitive to the water deficit. Mitochondrial responses under drought within Brassica genus are poorly understood. The main goal of this study was to investigate the mitochondrial biogenesis of three cauliflower cultivars varying with drought tolerance. Diverse quantitative changes (down-regulations mostly) in the mitochondrial proteome were assayed by 2D PAGE coupled with LC-MS/MS. Respiratory (e.g. CII, CIV and ATP synthase subunits), transporter (including diverse porin isoforms) and matrix multifunctional proteins (e.g. components of RNA editing machinery) appeared diversely affected in their abundance under two drought levels. Western immunoassays showed also cultivar-specific responses of selected mitochondrial proteins. Dehydrin-related tryptic peptides found in few 2D spots that appeared immunopositive with dehydrin-specific antisera highlighted the relevance of mitochondrial dehydrin-like proteins for the drought response. The level of selected messengers participating in drought response was also determined. We conclude that the mitochondrial biogenesis was strongly, but diversely affected in various cauliflower cultivars and associated with drought tolerance on the proteomic and functional levels. However, transcriptomic and proteomic regulations were largely uncoordinated due to the suggested altered availability of messengers for translation, mRNA/ribosome interactions and/or miRNA impact on transcript abundance and translation.

Pecan (Carya illinoinensis), as a popular nut tree, is widely planted in China in recent years. Grafting is an important technique for its cultivation. For a successful grafting, graft union development generally involves the formation of callus and vascular bundles at the graft union. To explore the molecular mechanism of graft union development, we applied high through-put RNA sequencing to investigate transcriptomic profiles of graft union at four time points (0d, 8d, 15d, and 30d) during pecan grafting process. We identified a total of 12,180 differentially expressed genes. In addition, we found that the content of auxin, cytokinin and gibberellin were accumulated at the graft unions during the grafting process. Correspondingly, genes involved in those hormone signaling were found to be differentially expressed. Interestingly, we found that most genes associated with cell division were up-regulated at callus formative stages, while genes related to cell elongation, secondary cell wall deposition, and programmed cell death were generally up-regulated at vascular bundle formative stages. In the meantime, genes responsible for reactive oxygen species were highly up-regulated across the graft union developmental process. These results will aid in our understanding of successful grafting in the future.

Nine crossbred finishing barrows randomly assigned to 3 dietary treatments were used to investigate the effects of dietary lysine on muscle growth related metabolic and signaling pathways. Muscle samples were collected from the longissimus dorsi of individual pigs after feeding the lysine-deficient, lysine-adequate, or lysine-excess diet for 5 weeks, and the total RNA was extracted afterwards. Affymetrix Porcine Gene 1.0 ST Array was used to quantify the expression levels of 19,211 genes. A total of 674 transcripts were differentially expressed (P ≤ 0.05); 60 out of 131 transcripts (P ≤ 0.01) were annotated in the NetAffx database. Ingenuity pathway analysis showed that dietary lysine deficiency may lead to (1) increased muscle protein degradation via the ubiquitination pathway as indicated by the up-regulated DNAJA1, HSP90AB1 and UBE2B mRNA, (2) reduced muscle protein synthesis via the up-regulated RND3 and ZIC1 mRNA, (3) increased serine and glycine synthesis via the up-regulated PHGDH and PSPH mRNA, and (4) increased lipid accumulation via the up-regulated ME1, SCD, and CIDEC mRNA. Dietary lysine excess may lead to (1) decreased muscle protein degradation via the down-regulated DNAJA1, HSP90AA1, HSPH1, and UBE2D3 mRNA, and (2) reduced lipid biosynthesis via the down-regulated CFD and ME1 mRNA.

A high intake of carbohydrates, associated with obesity, is one of the major causes of fatty liver disease in humans. This study investigated how high carbohydrate intake induces fatty liver disease in Blunt snout bream (Megalobrama amblycephala). Blunt snout bream were fed a high-carbohydrate diet (HCBD) for 60 days. Their growth indices were evaluated, and the transcriptomes, metabolites, biochemistry, and histology of their blood and livers were analyzed. The final weight, weight gain, specific growth rate, and feed conversion ratio were all higher in the HCBD group than in the control group, but not significantly so (P > 0.05). The hepatosomatic index (HSI) differed significantly in the two groups (P < 0.05), and the metabolomics results showed that a high carbohydrate intake induced significant increases in plasma α/β-glucose, succinate, and tyrosine, which could increase hepatic glycogen and triglyceride. Low levels of betaine were also found in the livers of the HCBD group. The histology and blood biochemistry results suggested abnormal liver, with excessive lipid accumulation and liver damage. A transcriptome analysis and quantitative reverse transcription–PCR (RT–qPCR) indicated that the expression of the factors INSR, IRS, PI3K, PDK, AKT, ACC, IL6, AP1, ChREBP-MLX, PEPCK, and FBP in the insulin signaling pathway was significantly upregulated and that of SOCS3, GSK3β, and AMPK significantly downregulated in the HCBD. This pattern is associated with the nonalcoholic fatty liver disease (NAFLD) pathway. This study extends our understanding of how high carbohydrate causes increased fat deposition in the liver, enhanced glycolysis (α/β-glucose) in the plasma, and reduced betaine in the liver. This leads to activation of hepatocyte insulin resistance and lipogenesis by regulating the expression of genes related to fatty liver disease.

Ectopic expression of the MYB transcription factor of AmROSEA1 from Antirrhinum majus has been reported to change anthocyanin and other metabolites in several species. In this study, we found that overexpression of AmRosea1 significantly improved the tolerance of transgenic rice to drought and salinity stress. Transcriptome analysis revealed that a considerable amount of stress-related genes were affected by exogenous AmRosea1 during both drought and salinity stress treatments. These affected genes are involved in stress signal transduction, the hormone signal pathway, ion homeostasis and the enzymes that remove peroxides. This work suggests that the AmRosea1 gene is a potential candidate for genetic engineering of crops.

Traditional methods for discovery and development of new drugs can be a very time-consuming and expensive process because it includes several stages such as compound identification, pre-clinical and clinical trials before the drug is approved by the US Food and Drug Administration (FDA). Therefore, drug repurposing, namely using currently FDA-approved drugs as therapeutics for other diseases than what they are originally prescribed for, is emerging to be a faster and more cost-effective alternative to current drug discovery methods. In this paper, we have described a three-step in silico protocol for analyzing transcriptomics data using online databases and bioinformatics tools for identifying potentially repurposable drugs. The efficacy of this protocol was evaluated by comparing its predictions with the findings of two case studies of recently reported repurposed drugs: HIV treating drug Zidovudine for the treatment of Dry Age-Related Macular Degeneration and the antidepressant Imipramine for Small-Cell Lung Carcinoma. The proposed protocol successfully identified the published findings, thus demonstrating the efficacy of this method. In addition, it also yielded several novel predictions that have not yet been published, including the finding that Imipramine could potentially treat Severe Acute Respiratory Syndrome (SARS), a disease that currently does not have any treatment or vaccine. Since this in-silico protocol is simple to use and does not require advanced computer skills, we believe any motivated participant with access to these databases and tools would be able to apply it to large datasets to identify other potentially repurposable drugs in the future.

Pardosa pseudoannulata is one of the most common wandering spiders in agricultural fields and a potentially good bioindicator for heavy metal contamination. However, little is known about the mechanism by which spiders respond to heavy metals at the molecular level. In this study, high-throughput transcriptome sequencing was employed to characterize the de novo transcriptome of the spiders and to identify differentially expressed genes (DEGs) after cadmium exposure. We obtained 60,489 assembled unigenes, 18,773 of which were annotated in the public databases. Ultimately, 3450 cDNA simple sequence repeats were identified and validated as potential molecular markers in the unigenes. A total of 2939, 2491 and 3759 DEGs were detected among the three libraries of two Cd-treated groups and the control. Functional enrichment analysis revealed that metabolism processes and digestive system function were predominately enriched in response to Cd stress. At the cellular and molecular levels, significantly enriched pathways in lysosomes and phagosomes as well as replication, recombination and repair demonstrated that oxidative damage resulted from Cd exposure. Based on the selected DEGs, certain critical genes involved in defence and detoxification were analysed. These results may elucidate the molecular mechanism underlying spiders' responses to heavy metal stress.

Lily tepals have a short lifespan. Once the tepals senesce, the ornamental value of the flower is lost. Some cultivars have attractive purple ovaries and fruits which greatly enhance the ornamental value of Asiatic hybrid lilies. However, little is known about the molecular mechanisms of anthocyanin biosynthesis in Asiatic hybrid lily ovaries. To investigate the transcriptional network that governs purple ovary coloration in Asiatic hybrid lilies, we obtained transcriptome data from green ovaries (S1) and purple ovaries (S2) of Asiatic ‘Tiny Padhye’. Comparative transcriptome analysis revealed 4228 differentially expressed genes. Differential expression analysis revealed that nine unigenes including four CHS genes, one CHI gene, one F3H gene, one F3′H gene, one DFR gene, and one UFGT gene were significantly up-regulated in purple ovaries. One MYB gene, LhMYB12-Lat, was identified as a key transcription factor determining the distribution of anthocyanins in Asiatic hybrid lily ovaries. Further qPCR results showed unigenes related to anthocyanin biosynthesis were highly expressed in purple ovaries of Asiatic ‘Tiny Padhye’ at stages 2 and 3, while they showed an extremely low level of expression in ovaries of Asiatic ‘Yellow Pixels’ during all developmental stages. In addition, shading treatment significantly decreased pigment accumulation by suppressing the expression of several unigenes related to anthocyanin biosynthesis in ovaries of Asiatic ‘Tiny Padhye’. These results could further our understanding of the molecular mechanisms of anthocyanin biosynthesis in Asiatic hybrid lily ovaries.