Effects of lithium and valproic acid on gene expression and phenotypic markers in an NT2 neurosphere model of neural development

Mood stabilising drugs such as lithium (LiCl) and valproic acid (VPA) are the first line agents for treating conditions such as Bipolar disorder and Epilepsy. However, these drugs have potential developmental effects that are not fully understood. This study explores the use of a simple human neurosphere-based in vitro model to characterise the pharmacological and toxicological effects of LiCl and VPA using gene expression changes linked to phenotypic alterations in cells. Treatment with VPA and LiCl resulted in the differential expression of 331 and 164 genes respectively. In the subset of VPA targeted genes, 114 were downregulated whilst 217 genes were upregulated. In the subset of LiCl targeted genes, 73 were downregulated and 91 were upregulated. Gene ontology (GO) term enrichment analysis was used to highlight the most relevant GO terms associated with a given gene list following toxin exposure. In addition, in order to phenotypically anchor the gene expression data, changes in the heterogeneity of cell subtype populations and cell cycle phase were monitored using flow cytometry. Whilst LiCl exposure did not significantly alter the proportion of cells expressing markers for stem cells/undifferentiated cells (Oct4, SSEA4), neurons (Neurofilament M), astrocytes (GFAP) or cell cycle phase, the drug caused a 1.4-fold increase in total cell number. In contrast, exposure to VPA resulted in significant upregulation of Oct4, SSEA, Neurofilament M and GFAP with significant decreases in both G2/M phase cells and cell number. This neurosphere model might provide the basis of a human-based cellular approach for the regulatory exploration of developmental impact of potential toxic chemicals.

POPLAR GENE EXPRESSION DATA ANALYSIS PIPELINES

Analyzing large-scale gene expression data is a labor-intensive and time-consuming process. To make data analysis easier, we developed a set of pipelines for rapid processing and analysis poplar gene expression data for knowledge discovery. Of all pipelines developed, differentially expressed genes (DEGs) pipeline is the one designed to identify biologically important genes that are differentially expressed in one of multiple time points for conditions. Pathway analysis pipeline was designed to identify the differentially expression metabolic pathways. Protein domain enrichment pipeline can identify the enriched protein domains present in the DEGs. Finally, Gene Ontology (GO) enrichment analysis pipeline was developed to identify the enriched GO terms in the DEGs. Our pipeline tools can analyze both microarray gene data and high-throughput gene data. These two types of data are obtained by two different technologies. A microarray technology is to measure gene expression levels via microarray chips, a collection of microscopic DNA spots attached to a solid (glass) surface, whereas high throughput sequencing, also called as the next-generation sequencing, is a new technology to measure gene expression levels by directly sequencing mRNAs, and obtaining each mRNA’s copy numbers in cells or tissues. We also developed a web portal (http://sys.bio.mtu.edu/) to make all pipelines available to public to facilitate users to analyze their gene expression data. In addition to the analyses mentioned above, it can also perform GO hierarchy analysis, i.e. construct GO trees using a list of GO terms as an input.

TaNF-YC11, one of the light-upregulated NF-YC members in Triticum aestivum, is co-regulated with photosynthesis-related genes

NF-Y is a heterotrimeric transcription factor complex. Each of the NF-Y subunits (NF-YA, NF-YB and NF-YC) in plants is encoded by multiple genes. Quantitative RT-PCR analysis revealed that five wheat NF-YC members (TaNF-YC5, 8, 9, 11 & 12) were upregulated by light in both the leaf and seedling shoot. Co-expression analysis of Affymetrix wheat genome array datasets revealed that transcript levels of a large number of genes were consistently correlated with those of the TaNF-YC11 and TaNF-YC8 genes in 3-4 separate Affymetrix array datasets. TaNF-YC11-correlated transcripts were significantly enriched with the Gene Ontology term photosynthesis. Sequence analysis in the promoters of TaNF-YC11-correlated genes revealed the presence of putative NF-Y complex binding sites (CCAAT motifs). Quantitative RT-PCR analysis of a subset of potential TaNF-YC11 target genes showed that ten out of the thirteen genes were also light-upregulated in both the leaf and seedling shoot and had significantly correlated expression profiles with TaNF-YC11. The potential target genes for TaNF-YC11 include subunit members from all four thylakoid membrane bound complexes required for the conversion of solar energy into chemical energy and rate limiting enzymes in the Calvin cycle. These data indicate that TaNF-YC11 is potentially involved in regulation of photosynthesis-related genes.

TaNF-YB3 is involved in the regulation of photosynthesis genes in Triticum aestivum

Nuclear Factor Y (NF-Y) transcription factor is a heterotrimer comprised of three subunits: NF-YA, NF-YB and NF-YC. Each of the three subunits in plants is encoded by multiple genes with differential expression profiles, implying the functional specialisation of NF-Y subunit members in plants. In this study, we investigated the roles of NF-YB members in the light-mediated regulation of photosynthesis genes. We identified two NF-YB members from Triticum aestivum (TaNF-YB3 & 7) which were markedly upregulated by light in the leaves and seedling shoots using quantitative RT-PCR. A genome-wide coexpression analysis of multiple Affymetrix Wheat Genome Array datasets revealed that TaNF-YB3-coexpressed transcripts were highly enriched with the Gene Ontology term photosynthesis. Transgenic wheat lines constitutively overexpressing TaNF-YB3 had a significant increase in the leaf chlorophyll content, photosynthesis rate and early growth rate. Quantitative RT-PCR analysis showed that the expression levels of a number of TaNF-YB3-coexpressed transcripts were elevated in the transgenic wheat lines. The mRNA level of TaGluTR encoding glutamyl-tRNA reductase, which catalyses the rate limiting step of the chlorophyll biosynthesis pathway, was significantly increased in the leaves of the transgenic wheat. Significant increases in the expression level in the transgenic plant leaves were also observed for four photosynthetic apparatus genes encoding chlorophyll a/b-binding proteins (Lhca4 and Lhcb4) and photosystem I reaction center subunits (subunit K and subunit N), as well as for a gene coding for chloroplast ATP synthase  subunit. These results indicate that TaNF-YB3 is involved in the positive regulation of a number of photosynthesis genes in wheat.

Massively parallel sequencing and analysis of expressed sequence tags in a successful invasive plant

Background Invasive species pose a significant threat to global economies, agriculture and biodiversity. Despite progress towards understanding the ecological factors associated with plant invasions, limited genomic resources have made it difficult to elucidate the evolutionary and genetic factors responsible for invasiveness. This study presents the first expressed sequence tag (EST) collection for Senecio madagascariensis, a globally invasive plant species. Methods We used pyrosequencing of one normalized and two subtractive libraries, derived from one native and one invasive population, to generate an EST collection. ESTs were assembled into contigs, annotated by BLAST comparison with the NCBI non-redundant protein database and assigned gene ontology (GO) terms from the Plant GO Slim ontologies. Key Results Assembly of the 221 746 sequence reads resulted in 12 442 contigs. Over 50 % (6183) of 12 442 contigs showed significant homology to proteins in the NCBI database, representing approx. 4800 independent transcripts. The molecular transducer GO term was significantly over-represented in the native (South African) subtractive library compared with the invasive (Australian) library. Based on NCBI BLAST hits and literature searches, 40 % of the molecular transducer genes identified in the South African subtractive library are likely to be involved in response to biotic stimuli, such as fungal, bacterial and viral pathogens. Conclusions This EST collection is the first representation of the S. madagascariensis transcriptome and provides an important resource for the discovery of candidate genes associated with plant invasiveness. The over-representation of molecular transducer genes associated with defence responses in the native subtractive library provides preliminary support for aspects of the enemy release and evolution of increased competitive ability hypotheses in this successful invasive. This study highlights the contribution of next-generation sequencing to better understanding the molecular mechanisms underlying ecological hypotheses that are important in successful plant invasions.

An ensemble method for predicting subnuclear localizations from primary protein structures

Background Predicting protein subnuclear localization is a challenging problem. Some previous works based on non-sequence information including Gene Ontology annotations and kernel fusion have respective limitations. The aim of this work is twofold: one is to propose a novel individual feature extraction method; another is to develop an ensemble method to improve prediction performance using comprehensive information represented in the form of high dimensional feature vector obtained by 11 feature extraction methods. Methodology/Principal Findings A novel two-stage multiclass support vector machine is proposed to predict protein subnuclear localizations. It only considers those feature extraction methods based on amino acid classifications and physicochemical properties. In order to speed up our system, an automatic search method for the kernel parameter is used. The prediction performance of our method is evaluated on four datasets: Lei dataset, multi-localization dataset, SNL9 dataset and a new independent dataset. The overall accuracy of prediction for 6 localizations on Lei dataset is 75.2% and that for 9 localizations on SNL9 dataset is 72.1% in the leave-one-out cross validation, 71.7% for the multi-localization dataset and 69.8% for the new independent dataset, respectively. Comparisons with those existing methods show that our method performs better for both single-localization and multi-localization proteins and achieves more balanced sensitivities and specificities on large-size and small-size subcellular localizations. The overall accuracy improvements are 4.0% and 4.7% for single-localization proteins and 6.5% for multi-localization proteins. The reliability and stability of our classification model are further confirmed by permutation analysis. Conclusions It can be concluded that our method is effective and valuable for predicting protein subnuclear localizations. A web server has been designed to implement the proposed method. It is freely available at http://bioinformatics.awowshop.com/snlpr​ed_page.php.

Urinary biomolecular indicators of exercise-induced over exertion injury

Poor health and injury represent major obstacles to the future economic security of Australia. The national economic cost of work-related injury is estimated at $57.5 billion p/a. Since exposure to high physical demands is a major risk factor for musculoskeletal injury, monitoring and managing such physical activity levels in workers is a potentially important injury prevention strategy. Current injury monitoring practices are inadequate for the provision of clinically valuable information about the tissue specific responses to physical exertion. Injury of various soft tissue structures can manifest over time through accumulation of micro-trauma. Such micro-trauma has a propensity to increase the risk of acute injuries to soft-tissue structures such as muscle or tendon. As such, the capacity to monitor biomarkers that result from the disruption of these tissues offers a means of assisting the pre-emptive management of subclinical injury prior to acute failure or for evaluation of recovery processes. Here we have adopted an in-vivo exercise induced muscle damage model allowing the application of laboratory controlled conditions to assist in uncovering biochemical indicators associated with soft-tissue trauma and recovery. Importantly, urine was utilised as the diagnostic medium since it is non-invasive to collect, more acceptable to workers and less costly to employers. Moreover, it is our hypothesis that exercise induced tissue degradation products enter the circulation and are subsequently filtered by the kidney and pass through to the urine. To test this hypothesis a range of metabolomic and proteomic discovery-phase techniques were used, along with targeted approaches. Several small molecules relating to tissue damage were identified along with a series of skeletal muscle-specific protein fragments resulting from exercise induced soft-tissue damage. Each of the potential biomolecular markers appeared to be temporally present within urine. Moreover, the regulation of abundance seemed to be associated with functional recovery following the injury. This discovery may have important clinical applications for monitoring of a variety of inflammatory myopathies as well as novel applications in monitoring of the musculoskeletal health status of workers, professional athletes and/or military personnel to reduce the onset of potentially debilitating musculoskeletal injuries within these professions.

Modeling the cumulative genetic risk for multiple sclerosis from genome-wide association data

Background: Multiple sclerosis (MS) is the most common cause of chronic neurologic disability beginning in early to middle adult life. Results from recent genome-wide association studies (GWAS) have substantially lengthened the list of disease loci and provide convincing evidence supporting a multifactorial and polygenic model of inheritance. Nevertheless, the knowledge of MS genetics remains incomplete, with many risk alleles still to be revealed. Methods: We used a discovery GWAS dataset (8,844 samples, 2,124 cases and 6,720 controls) and a multi-step logistic regression protocol to identify novel genetic associations. The emerging genetic profile included 350 independent markers and was used to calculate and estimate the cumulative genetic risk in an independent validation dataset (3,606 samples). Analysis of covariance (ANCOVA) was implemented to compare clinical characteristics of individuals with various degrees of genetic risk. Gene ontology and pathway enrichment analysis was done using the DAVID functional annotation tool, the GO Tree Machine, and the Pathway-Express profiling tool. Results: In the discovery dataset, the median cumulative genetic risk (P-Hat) was 0.903 and 0.007 in the case and control groups, respectively, together with 79.9% classification sensitivity and 95.8% specificity. The identified profile shows a significant enrichment of genes involved in the immune response, cell adhesion, cell communication/ signaling, nervous system development, and neuronal signaling, including ionotropic glutamate receptors, which have been implicated in the pathological mechanism driving neurodegeneration. In the validation dataset, the median cumulative genetic risk was 0.59 and 0.32 in the case and control groups, respectively, with classification sensitivity 62.3% and specificity 75.9%. No differences in disease progression or T2-lesion volumes were observed among four levels of predicted genetic risk groups (high, medium, low, misclassified). On the other hand, a significant difference (F = 2.75, P = 0.04) was detected for age of disease onset between the affected misclassified as controls (mean = 36 years) and the other three groups (high, 33.5 years; medium, 33.4 years; low, 33.1 years). Conclusions: The results are consistent with the polygenic model of inheritance. The cumulative genetic risk established using currently available genome-wide association data provides important insights into disease heterogeneity and completeness of current knowledge in MS genetics.

Expression Profiles of Mitochondrial Genes in the Frontal Cortex and the Caudate Nucleus of Developing Humans and Mice Selectively Bred for High and Low Fear

A growing body of evidence suggests that mitochondrial function may be important in brain development and psychiatric disorders. However, detailed expression profiles of those genes in human brain development and fear-related behavior remain unclear. Using microarray data available from the public domain and the Gene Ontology analysis, we identified the genes and the functional categories associated with chronological age in the prefrontal cortex (PFC) and the caudate nucleus (CN) of psychiatrically normal humans ranging in age from birth to 50 years. Among those, we found that a substantial number of genes in the PFC (115) and the CN (117) are associated with the GO term: mitochondrion (FDR qv <0.05). A greater number of the genes in the PFC (91%) than the genes in the CN (62%) showed a linear increase in expression during postnatal development. Using quantitative PCR, we validated the developmental expression pattern of four genes including monoamine oxidase B (MAOB), NADH dehydrogenase flavoprotein (NDUFV1), mitochondrial uncoupling protein 5 (SLC25A14) and tubulin beta-3 chain (TUBB3). In mice, overall developmental expression pattern of MAOB, SLC25A14 and TUBB3 in the PFC were comparable to the pattern observed in humans (p<0.05). However, mice selectively bred for high fear did not exhibit normal developmental changes of MAOB and TUBB3. These findings suggest that the genes associated with mitochondrial function in the PFC play a significant role in brain development and fear-related behavior.

The Anadara trapezia transcriptome : a resource for molluscan physiological genomics

In this study we undertook deep sequencing of the blood cockle, Anadara trapezia, transcriptome to generate genomic resources for future functional genomics analyses. Over 27 million high quality paired end reads were assembled into 75 024 contigs. Of these contigs, 29 013 (38.7%) received significant BLASTx hits and gene ontology (GO) terms were assigned to 13 718 of these sequences. This resourcewill facilitate physiological genomic studies to test the gene expression response of A. trapezia to various environmental stresses.

Sexual selection in true fruit flies (Diptera: Tephritidae) : transcriptome and experimental evidences for phytochemicals increasing male competitive ability

In male tephritid fruit flies of the genus Bactrocera, feeding on secondary plant compounds (sensu lato male lures = methyl eugenol, raspberry ketone and zingerone) increases male mating success. Ingested male lures alter the male pheromonal blend, normally making it more attractive to females and this is considered the primary mechanism for the enhanced mating success. However, the male lures raspberry ketone and zingerone are known, across a diverse range of other organisms, to be involved in increasing energy metabolism. If this also occurs in Bactrocera, then this may represent an additional benefit to males as courtship is metabolically expensive and lure feeding may increase a fly's short-term energy. We tested this hypothesis by performing comparative RNA-seq analysis between zingerone-fed and unfed males of Bactrocera tryoni. We also carried out behavioural assays with zingerone- and cuelure-fed males to test whether they became more active. RNA-seq analysis revealed, in zingerone-fed flies, up-regulation of 3183 genes with homologues transcripts to those known to regulate intermale aggression, pheromone synthesis, mating and accessory gland proteins, along with significant enrichment of several energy metabolic pathways and gene ontology terms. Behavioural assays show significant increases in locomotor activity, weight reduction and successful mating after mounting; all direct/indirect measures of increased activity. These results suggest that feeding on lures leads to complex physiological changes, which result in more competitive males. These results do not negate the pheromone effect, but do strongly suggest that the phytochemical-induced sexual selection is governed by both female preference and male competitive mechanisms.

Molecular profiling of exudate from chronic ulcerated wounds

Non-healing wounds represent a significant burden to healthcare systems and societies worldwide. Current best practice treatments of chronic wounds can require patients to undergo extensive periods of therapy without any positive outcome. This consumes substantial healthcare resources and severely impacts patient quality of life. At present, there are no measures to predict a patient's response to best practice care. The hypothesis of this thesis was that biochemical markers could be found within the wound fluid of chronic ulcers and these markers could predict the healing outcome of an ulcer undergoing best practice care. Discovery phase proteomic and mass spectrometry techniques were utilised to determine novel proteins that correlated with the healing outcome of ulcers. These candidate biomarkers could be developed into simple dip-stick tools for use in clinical practice. This would aid clinicians in the choice of effective wound management strategies to address hard-to-heal wounds.

The tumour-promoting receptor tyrosine kinase, EphB4, regulates expression of Integrin-β8 in prostate cancer cells

Background The EphB4 receptor tyrosine kinase is overexpressed in many cancers including prostate cancer. The molecular mechanisms by which this ephrin receptor influences cancer progression are complex as there are tumor-promoting ligand-independent mechanisms in place as well as ligand-dependent tumor suppressive pathways. Methods We employed transient knockdown of EPHB4 in prostate cancer cells, coupled with gene microarray analysis, to identify genes that were regulated by EPHB4 and may represent linked tumor-promoting factors. We validated target genes using qRT-PCR and employed functional assays to determine their role in prostate cancer migration and invasion. Results We discovered that over 500 genes were deregulated upon EPHB4 siRNA knockdown, with integrin β8 (ITGB8) being the top hit (29-fold down-regulated compared to negative non-silencing siRNA). Gene ontology analysis found that the process of cell adhesion was highly deregulated and two other integrin genes, ITGA3 and ITGA10, were also differentially expressed. In parallel, we also discovered that over-expression of EPHB4 led to a concomitant increase in ITGB8 expression. In silico analysis of a prostate cancer progression microarray publically available in the Oncomine database showed that both EPHB4 and ITGB8 are highly expressed in prostatic intraepithelial neoplasia, the precursor to prostate cancer. Knockdown of ITGB8 in PC-3 and 22Rv1 prostate cancer cells in vitro resulted in significant reduction of cell migration and invasion. Conclusions These results reveal that EphB4 regulates integrin β8 expression and that integrin β8 plays a hitherto unrecognized role in the motility of prostate cancer cells and thus targeting integrin β8 may be a new treatment strategy for prostate cancer.

Analysis, characterisation and expression of gill-expressed carbonic anhydrase genes in the freshwater crayfish Cherax quadricarinatus

Changes in water quality parameters such as pH and salinity can have a significant effect on productivity of aquaculture species. Similarly, relative osmotic pressure influences various physiological processes and regulates expression of a number of osmoregulatory genes. Among those, carbonic anhydrase (CA) plays a key role in systemic acid–base balance and ion regulation. Redclaw crayfish (Cherax quadricarinatus) are unique in their ability to thrive in environments with naturally varied pH levels, suggesting unique adaptation to pH stress. To date, however, no studies have focused on identification and characterisation of CA or other osmoregulatory genes in C. quadricarinatus. Here, we analysed the redclaw gill transcriptome and characterized CA genes along with a number of other key osmoregulatory genes that were identified in the transcriptome. We also examined patterns of gene expression of these CA genes when exposed to three pH treatments. In total, 72,382,710 paired end Illumina reads were assembled into 36,128 contigs with an average length of 800 bp. Approximately 37% of contigs received significant BLAST hits and 22% were assigned gene ontology terms. Three full length CA isoforms; cytoplasmic CA (ChqCAc), glycosyl-phosphatidylinositol-linked CA (ChqCAg), and β-CA (ChqCA-beta) as well as two partial CA gene sequences were identified. Both partial CA genes showed high similarity to ChqCAg and appeared to be duplicated from the ChqCAg. Full length coding sequences of Na+/K+-ATPase, V-type H+-ATPase, sarcoplasmic Ca+-ATPase, arginine kinase, calreticulin and Cl− channel protein 2 were also identified. Only the ChqCAc gene showed significant differences in expression across the three pH treatments. These data provide valuable information on the gill expressed CA genes and their expression patterns in freshwater crayfish. Overall our data suggest an important role for the ChqCAc gene in response to changes in pH and in systemic acid–base balance in freshwater crayfish.