Bayesian inference of the number of factors in gene-expression analysis: application to human virus challenge studies.

BACKGROUND: Nonparametric Bayesian techniques have been developed recently to extend the sophistication of factor models, allowing one to infer the number of appropriate factors from the observed data. We consider such techniques for sparse factor analysis, with application to gene-expression data from three virus challenge studies. Particular attention is placed on employing the Beta Process (BP), the Indian Buffet Process (IBP), and related sparseness-promoting techniques to infer a proper number of factors. The posterior density function on the model parameters is computed using Gibbs sampling and variational Bayesian (VB) analysis. RESULTS: Time-evolving gene-expression data are considered for respiratory syncytial virus (RSV), Rhino virus, and influenza, using blood samples from healthy human subjects. These data were acquired in three challenge studies, each executed after receiving institutional review board (IRB) approval from Duke University. Comparisons are made between several alternative means of per-forming nonparametric factor analysis on these data, with comparisons as well to sparse-PCA and Penalized Matrix Decomposition (PMD), closely related non-Bayesian approaches. CONCLUSIONS: Applying the Beta Process to the factor scores, or to the singular values of a pseudo-SVD construction, the proposed algorithms infer the number of factors in gene-expression data. For real data the "true" number of factors is unknown; in our simulations we consider a range of noise variances, and the proposed Bayesian models inferred the number of factors accurately relative to other methods in the literature, such as sparse-PCA and PMD. We have also identified a "pan-viral" factor of importance for each of the three viruses considered in this study. We have identified a set of genes associated with this pan-viral factor, of interest for early detection of such viruses based upon the host response, as quantified via gene-expression data.

Specialized motor-driven dusp1 expression in the song systems of multiple lineages of vocal learning birds.

Mechanisms for the evolution of convergent behavioral traits are largely unknown. Vocal learning is one such trait that evolved multiple times and is necessary in humans for the acquisition of spoken language. Among birds, vocal learning is evolved in songbirds, parrots, and hummingbirds. Each time similar forebrain song nuclei specialized for vocal learning and production have evolved. This finding led to the hypothesis that the behavioral and neuroanatomical convergences for vocal learning could be associated with molecular convergence. We previously found that the neural activity-induced gene dual specificity phosphatase 1 (dusp1) was up-regulated in non-vocal circuits, specifically in sensory-input neurons of the thalamus and telencephalon; however, dusp1 was not up-regulated in higher order sensory neurons or motor circuits. Here we show that song motor nuclei are an exception to this pattern. The song nuclei of species from all known vocal learning avian lineages showed motor-driven up-regulation of dusp1 expression induced by singing. There was no detectable motor-driven dusp1 expression throughout the rest of the forebrain after non-vocal motor performance. This pattern contrasts with expression of the commonly studied activity-induced gene egr1, which shows motor-driven expression in song nuclei induced by singing, but also motor-driven expression in adjacent brain regions after non-vocal motor behaviors. In the vocal non-learning avian species, we found no detectable vocalizing-driven dusp1 expression in the forebrain. These findings suggest that independent evolutions of neural systems for vocal learning were accompanied by selection for specialized motor-driven expression of the dusp1 gene in those circuits. This specialized expression of dusp1 could potentially lead to differential regulation of dusp1-modulated molecular cascades in vocal learning circuits.

Examination of Endogenous Rotund Expression and Function in Developing Drosophila Olfactory System Using CRISPR-Cas9-Mediated Protein Tagging.

The zinc-finger protein Rotund (Rn) plays a critical role in controlling the development of the fly olfactory system. However, little is known about its molecular function in vivo. Here, we added protein tags to the rn locus using CRISPR-Cas9 technology in Drosophila to investigate its subcellular localization and the genes that it regulates . We previously used a reporter construct to show that rn is expressed in a subset of olfactory receptor neuron (ORN) precursors and it is required for the diversification of ORN fates. Here, we show that tagged endogenous Rn protein is functional based on the analysis of ORN phenotypes. Using this method, we also mapped the expression pattern of the endogenous isoform-specific tags in vivo with increased precision. Comparison of the Rn expression pattern from this study with previously published results using GAL4 reporters showed that Rn is mainly present in early steps in antennal disc patterning, but not in pupal stages when ORNs are born. Finally, using chromatin immunoprecipitation, we showed a direct binding of Rotund to a previously identified regulatory element upstream of the bric-a-brac gene locus in the developing antennal disc.

Individual mRNA expression profiles reveal the effects of specific microRNAs

BACKGROUND: MicroRNAs (miRNAs) are oligoribonucleotides with an important role in regulation of gene expression at the level of translation. Despite imperfect target complementarity, they can also significantly reduce mRNA levels. The validity of miRNA target gene predictions is difficult to assess at the protein level. We sought, therefore, to determine whether a general lowering of predicted target gene mRNA expression by endogenous miRNAs was detectable within microarray gene expression profiles. RESULTS: The target gene sets predicted for each miRNA were mapped onto known gene expression data from a range of tissues. Whether considering mean absolute target gene expression, rank sum tests or 'ranked ratios', many miRNAs with significantly reduced target gene expression corresponded to those known to be expressed in the cognate tissue. Expression levels of miRNAs with reduced target mRNA levels were higher than those of miRNAs with no detectable effect on mRNA expression. Analysis of microarray data gathered after artificial perturbation of expression of a specific miRNA confirmed the predicted increase or decrease in influence of the altered miRNA upon mRNA levels. Strongest associations were observed with targets predicted by TargetScan. CONCLUSION: We have demonstrated that the effect of a miRNA on its target mRNAs' levels can be measured within a single gene expression profile. This emphasizes the extent of this mode of regulation in vivo and confirms that many of the predicted miRNA-mRNA interactions are correct. The success of this approach has revealed the vast potential for extracting information about miRNA function from gene expression profiles.

Achieving hypoxia-inducible gene expression in tumors

Hypoxia is an inevitable feature of solid tumors and a common cause of treatment failure. Hypoxia acts as a trigger to genetic instability, apoptosis and possibly metastases. The adaptive response to cellular hypoxia involves the modulation of the synthesis of multiple proteins controlling processes such as glucose homeostasis, angiogenesis, vascular permeability and inflammation. The hypoxia responsive element (HRE) sequences isolated from oxygen-responsive genes have been shown to selectively induce gene expression in response to hypoxia when placed upstream of a promoter. The levels of induced gene expression were dependent on the number of HRE copies and the oxygen tension. Hypoxia-mediated cancer gene therapy strategies may represent a promising mean to significantly improve the efficacy of standard radiation therapy and chemotherapy approaches.

Chromatin immunoprecipitation (ChIP) methodology and readouts

The identification of direct nuclear hormone receptor gene targets provides clues to their contribution to both development and cancer progression. Until recently, the identification of such direct target genes has relied on a combination of expression analysis and in silico searches for consensus binding motifs in gene promoters. Consensus binding motifs for transcription factors are often defined using in vitro DNA binding strategies. Such in vitro strategies fail to account for the many factors that contribute significantly to target selection by transcription factors in cells beyond the recognition of these short consensus DNA sequences. These factors include DNA methylation, chromatin structure, posttranslational modifications of transcription factors, and the cooperative recruitment of transcription factor complexes. Chromatin immunoprecipitation (ChIP) provides a means of isolating transcription factor complexes in the context of endogenous chromatin, allowing the identification of direct transcription factor targets. ChIP can be combined with site-specific PCR for candidate binding sites or alternatively with cloning, genomic microarrays or more recently direct high throughput sequencing to identify novel genomic targets. The application of ChIP-based approaches has redefined consensus binding motifs for transcription factors and provided important insights into transcription factor biology.

An atlas of human gene expression from massively parallel signature sequencing (MPSS).

We have used massively parallel signature sequencing (MPSS) to sample the transcriptomes of 32 normal human tissues to an unprecedented depth, thus documenting the patterns of expression of almost 20,000 genes with high sensitivity and specificity. The data confirm the widely held belief that differences in gene expression between cell and tissue types are largely determined by transcripts derived from a limited number of tissue-specific genes, rather than by combinations of more promiscuously expressed genes. Expression of a little more than half of all known human genes seems to account for both the common requirements and the specific functions of the tissues sampled. A classification of tissues based on patterns of gene expression largely reproduces classifications based on anatomical and biochemical properties. The unbiased sampling of the human transcriptome achieved by MPSS supports the idea that most human genes have been mapped, if not functionally characterized. This data set should prove useful for the identification of tissue-specific genes, for the study of global changes induced by pathological conditions, and for the definition of a minimal set of genes necessary for basic cell maintenance. The data are available on the Web at http://mpss.licr.org and http://sgb.lynxgen.com.

Downregulation of VAPB expression in motor neurons derived from induced pluripotent stem cells of ALS8 patients

Amyotrophic lateral sclerosis (ALS) is an incurable neuromuscular disease that leads to a profound loss of life quality and premature death. Around 10% of the cases are inherited and ALS8 is an autosomal dominant form of familial ALS caused by mutations in the vamp-associated protein B/C (VAPB) gene. The VAPB protein is involved in many cellular processes and it likely contributes to the pathogenesis of other forms of ALS besides ALS8. A number of successful drug tests in ALS animal models could not be translated to humans underscoring the need for novel approaches. The induced pluripotent stem cells (iPSC) technology brings new hope, since it can be used to model and investigate diseases in vitro. Here we present an additional tool to study ALS based on ALS8-iPSC. Fibroblasts from ALS8 patients and their non-carrier siblings were successfully reprogrammed to a pluripotent state and differentiated into motor neurons. We show for the first time that VAPB protein levels are reduced in ALS8-derived motor neurons but, in contrast to over-expression systems, cytoplasmic aggregates could not be identified. Our results suggest that optimal levels of VAPB may play a central role in the pathogenesis of ALS8, in agreement with the observed reduction of VAPB in sporadic ALS.

Molecular analysis of maltotriose active transport and fermentation by Saccharomyces cerevisiae reveals a determinant role for the AGT1 permease

Incomplete and/or sluggish maltotriose fermentation causes both quality and economic problems in the ale-brewing industry. Although it has been proposed previously that the sugar uptake must be responsible for these undesirable phenotypes, there have been conflicting reports on whether all the known alpha-glucoside transporters in Saccharomyces cerevisiae (MALx1, AGT1, and MPH2 and MPH3 transporters) allow efficient maltotriose utilization by yeast cells. We characterized the kinetics of yeast cell growth, sugar consumption, and ethanol production during maltose or maltotriose utilization by several S. cerevisiae yeast strains (both MAL constitutive and AM inducible) and by their isogenic counterparts with specific deletions of the AGT1 gene. Our results clearly showed that yeast strains carrying functional permeases encoded by the MAL21, MAL31, and/or MAL41 gene in their plasma membranes were unable to utilize maltotriose. While both high-and low-affinity transport activities were responsible for maltose uptake from the medium, in the case of maltotriose, the only low-affinity (K-m, 36 +/- 2 mM) transport activity was mediated by the AGT1 permease. In conclusion, the AGT1 transporter is required for efficient maltotriose fermentation by S. cerevisiae yeasts, highlighting the importance of this permease for breeding and/or selection programs aimed at improving sluggish maltotriose fermentations.

Development of a molecular test to determine the vitality status of Norway spruce (Picea abies) seedlings during frozen storage

In boreal forest regions, a great portion of forest tree seedlings are stored indoors in late autumn to prevent seedlings from outdoor winter damage. For seedlings to be able to survive in storage it is crucial that they store well and can cope with the dark and cold storage environment. The aim of this study was to search for genes that can determine the vitality status of Norway spruce (Picea abies (L.) Karst.) seedlings during frozen storage. Furthermore, the sensitivity of the ColdNSure (TM) test, a gene activity test that predicts storability was assessed. The storability of seedlings was tested biweekly by evaluating damage with the gene activity test and the electrolyte leakage test after freezing seedlings to -25 A degrees C (the SELdiff-25 method). In parallel, seedlings were frozen stored at -3 A degrees C. According to both methods, seedlings were considered storable from week 41. This also corresponded to the post storage results determined at the end of the storage period. In order to identify vitality indicators, Next Generation Sequencing (NGS) was performed on bud samples collected during storage. Comparing physiological post storage data to gene analysis data revealed numerous vitality related genes. To validate the results, a second trial was performed. In this trial, gene activity was better in predicting seedling storability than the conventional freezing test; this indicates a high sensitivity level of this molecular assay. For multiple indicators a clear switch between damaged and vital seedlings was observed. A collection of indicators will be used in the future development of a commercial vitality test.

A custom-built insulin resistance gene chip

Objectives/Aim—Microarray (gene chip) technology offers a powerful new tool for analyzing the expression of large numbers of genes in many experimental samples. The aim of this study was to design, construct, and use a gene chip to measure the expression levels of key genes in metabolic pathways related to insulin resistance.
Methods—We selected genes that were implicated in the development of insulin resistance, including genes involved in insulin signaling; glucose uptake, oxidation, and storage; fat uptake, oxidation, and storage; cytoskeletal components; and transcription factors. The key regulatory genes in the pathways were identified, along with other recently identified candidate genes such as calpain-10. A total of 242 selected genes (including 32 internal control elements) were sequence-verified, purified, and arrayed on aldehyde-coated slides.
Results—Where more than 1 clone containing the gene of interest was available, we chose those containing the genes in the 5' orientation and an insert size of around 1.5 kb. Of the 262 clones purchased, 56 (21%) were found to contain sequences other than those expected. In addition, 2 (1%) did not grow under standard conditions and were assumed to be nonviable. In these cases, alternate clones containing the gene of interest were chosen as described above. The current version of the Insulin Resistance Gene Chip contains 210 genes of interest, plus 48 control elements. A full list of the genes is available at http://www.hbs.deakin.edu.au/mru/research/gene_chip_tech/genechip_three.htm/.
Conclusions—The human Insulin Resistance Gene Chip that we have constructed will be a very useful tool for investigating variation in the expression of genes relevant to insulin resistance under various experimental conditions. Initially, the gene chip will be used in studies such as exercise interventions, fasting, euglycemic-hyperinsulinemic clamps, and administration of antidiabetic agents

Zinc and DHA have opposing effects on the expression levels of histones H3 and H4 in human neuronal cells

Zn and DHA have putative neuroprotective effects and these two essential nutrients are known to interact biochemically. We aimed to identify novel protein candidates that are differentially expressed in human neuronal cell line M17 in response to Zn and DHA that would explain the molecular basis of this interaction. Two-dimensional gel electrophoresis and MS were applied to identify major protein expression changes in the protein lysates of human Ml7 neuronal cells that had been grown in the presence and absence of Zn and DHA. Proteomic findings were further investigated using Western immunoblot and real-time PCR analyses. Four protein spots, which had significant differential expression, were identified and selected for in-gel trypsin digestion followed by matrix-assisted laser desorption ionisation MS analysis. The resultant peptide mass fingerprint for each spot allowed their respective identities to be deduced. Two human histone variants H3 and H4 were identified. Both H3 and H4 were downregulated by Zn in the absence of DHA (Zn effect) and upregulated by DHA (DHA effect) in the presence of Zn (physiological condition). These proteomic findings were further supported by Western immunoblot and real-time PCR analyses using H3- and H4-specific monoclonal antibodies and oligonucleotide primers, respectively. We propose that dietary Zn and DHA cause a global effect on gene expression, which is mediated by histones. Such novel information provides possible clues to the molecular basis of neuroprotection by Zn and DHA that may contribute to the future treatment, prevention and management of neurodegenerative diseases such as Alzheimer's disease.

Molecular detection of enterotoxins E, G, H and I in Staphylococcus aureus and coagulase-negative staphylococci isolated from clinical samples of newborns in Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)