Search

Background: Insects evolved ingenious adaptations to use extraordinary food sources. Particularly, the diet of herbivores enriched with noxious plant secondary metabolites requires detoxification mechanisms. Sequestration, which involves the uptake, transfer, and concentration of occasionally modified phytochemicals into specialized tissues or hemolymph, is one of the most successful detoxification strategies found in most insect orders. Due to the ability of ATP-binding cassette (ABC) carriers to transport a wide range of molecules including phytochemicals and xenobiotics, it is highly likely that they play a role in this sequestration process. To shed light on the role of ABC proteins in sequestration, we describe an inventory of putative ABC transporters in various tissues in the sequestering juvenile poplar leaf beetle, Chrysomela populi. Results: In the transcriptome of C. populi, we predicted 65 ABC transporters. To link the proteins with a possible function, we performed comparative phylogenetic analyses with ABC transporters of other insects and of humans. While tissue-specific profiling of each ABC transporter subfamily suggests that ABCB, C and G influence the plant metabolite absorption in the gut, ABCC with 14 members is the preferred subfamily responsible for the excretion of these metabolites via Malpighian tubules. Moreover, salicin, which is sequestered from poplar plants, is translocated into the defensive glands for further deterrent production. In these glands and among all identified ABC transporters, an exceptionally high transcript level was observed only for Cpabc35 (Cpmrp). RNAi revealed the deficiency of other ABC pumps to compensate the function of CpABC35, demonstrating its key role during sequestration. Conclusion: We provide the first comprehensive phylogenetic study of the ABC family in a phytophagous beetle species. RNA-seq data from different larval tissues propose the importance of ABC pumps to achieve a homeostasis of plant-derived compounds and offer a basis for future analyses of their physiological function in sequestration processes. [ABSTRACT FROM AUTHOR]

There is growing evidence for nongenetic effects of maternal experience on offspring. For example, previous studies have shown that female threespined stickleback fish (Gasterosteus aculeatus) exposed to predation risk produce offspring with altered behavior, metabolism and stress physiology. Here, we investigate the effect of maternal exposure to predation risk on the embryonic transcriptome in sticklebacks. Using RNA-sequencing we compared genome-wide transcription in three day post-fertilization embryos of predator-exposed and control mothers. There were hundreds of differentially expressed transcripts between embryos of predator-exposed mothers and embryos of control mothers including several non-coding RNAs. Gene Ontology analysis revealed biological pathways involved in metabolism, epigenetic inheritance, and neural proliferation and differentiation that differed between treatments. Interestingly, predation risk is associated with an accelerated life history in many vertebrates, and several of the genes and biological pathways that were identified in this study suggest that maternal exposure to predation risk accelerates the timing of embryonic development. Consistent with this hypothesis, embryos of predator-exposed mothers were larger than embryos of control mothers. These findings point to some of the molecular mechanisms that might underlie maternal effects. [ABSTRACT FROM AUTHOR]

The current replication study confirmed significant differences in gene expression profiles of the cerebrum among the two commercial layer lines Lohmann Selected Leghorn (LSL) and Lohmann Brown (LB). Microarray analyses were performed for 30 LSL and another 30 LB laying hens kept in the small group housing system Eurovent German. A total of 14,103 microarray probe sets using customized Affymetrix ChiGene-1_0-st Arrays with 20,399 probe sets were differentially expressed among the two layer lines LSL and LB (FDR adjusted P-value <0.05). An at least 2-fold change in expression levels could be observed for 388 of these probe sets. In LSL, 214 of the 388 probe sets were down- and 174 were up-regulated and vice versa for the LB layer line. Among the 174 up-regulated probe sets in LSL, we identified 51 significantly enriched Gene ontology (GO) terms of the biological process category. A total of 63 enriched GO-terms could be identified for the 214 down-regulated probe sets of the layer line LSL. We identified nine genes significantly differentially expressed between the two layer lines in both microarray experiments. These genes play a crucial role in protection of neuronal cells from oxidative stress, bone mineral density and immune response among the two layer lines LSL and LB. Thus, the different regulation of these genes may significantly contribute to phenotypic trait differences among these layer lines. In conclusion, these novel findings provide a basis for further research to improve animal welfare in laying hens and these layer lines may be of general interest as an animal model. [ABSTRACT FROM AUTHOR]

Background: Originating from Primordial Germ Cells/gonocytes and developing via a precursor lesion called Carcinoma In Situ (CIS), Germ Cell Cancers (GCC) are the most common cancer in young men, subdivided in seminoma (SE) and non-seminoma (NS). During physiological germ cell formation/maturation, epigenetic processes guard homeostasis by regulating the accessibility of the DNA to facilitate transcription. Epigenetic deregulation through genetic and environmental parameters (i.e. genvironment) could disrupt embryonic germ cell development, resulting in delayed or blocked maturation. This potentially facilitates the formation of CIS and progression to invasive GCC. Therefore, determining the epigenetic and functional genomic landscape in GCC cell lines could provide insight into the pathophysiology and etiology of GCC and provide guidance for targeted functional experiments. Results: This study aims at identifying epigenetic footprints in SE and EC cell lines in genome-wide profiles by studying the interaction between gene expression, DNA CpG methylation and histone modifications, and their function in the pathophysiology and etiology of GCC. Two well characterized GCC-derived cell lines were compared, one representative for SE (TCam-2) and the other for EC (NCCIT). Data were acquired using the Illumina HumanHT-12-v4 (gene expression) and HumanMethylation450 BeadChip (methylation) microarrays as well as ChIP-sequencing (activating histone modifications (H3K4me3, H3K27ac)). Results indicate known germ cell markers not only to be differentiating between SE and NS at the expression level, but also in the epigenetic landscape. Conclusion: The overall similarity between TCam-2/NCCIT support an erased embryonic germ cell arrested in early gonadal development as common cell of origin although the exact developmental stage from which the tumor cells are derived might differ. Indeed, subtle difference in the (integrated) epigenetic and expression profiles indicate TCam-2 to exhibit a more germ cell-like profile, whereas NCCIT shows a more pluripotent phenotype. The results provide insight into the functional genome in GCC cell lines. [ABSTRACT FROM AUTHOR]

Genome-wide association studies (GWASs) may help to understand the effects of genetic polymorphisms on breast cancer (BC) progression and survival. However, they give only a focused view, which cannot capture the tremendous complexity of this disease. Therefore, we investigated data from a previously conducted GWAS on BC survival for enriched pathways by different enrichment analysis tools using the two main annotation databases Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The goal was to identify the functional categories (GO terms and KEGG pathways) that are consistently overrepresented in a statistically significant way in the list of genes generated from the single nucleotide polymorphism (SNP) data. The SNPs with allelic p-value cut-offs 0.005 and 0.01 were annotated to the genes by excluding or including a 20 kb up-and down-stream sequence of the genes and analyzed by six different tools. We identified eleven consistently enriched categories, the most significant ones relating to cell adhesion and calcium ion binding. Moreover, we investigated the similarity between our GWAS and the enrichment analyses of twelve published gene expression signatures for breast cancer prognosis. Five of them were commonly used and commercially available, five were based on different aspects of metastasis formation and two were developed from meta-analyses of published prognostic signatures. This comparison revealed similarities between our GWAS data and the general and the specific brain metastasis gene signatures as well as the Oncotype DX signature. As metastasis formation is a strong indicator of a patient’s prognosis, this result reflects the survival aspect of the conducted GWAS and supports cell adhesion and calcium signaling as important pathways in cancer progression. [ABSTRACT FROM AUTHOR]

Microsporidia comprise a phylum of over 1400 species of obligate intracellular pathogens that can infect almost all animals, but little is known about the host response to these parasites. Here we use the whole-animal host C. elegans to show an in vivo role for ubiquitin-mediated response to the microsporidian species Nematocida parisii, as well to the Orsay virus, another natural intracellular pathogen of C. elegans. We analyze gene expression of C. elegans in response to N. parisii, and find that it is similar to response to viral infection. Notably, we find an upregulation of SCF ubiquitin ligase components, such as the cullin ortholog cul-6, which we show is important for ubiquitin targeting of N. parisii cells in the intestine. We show that ubiquitylation components, the proteasome, and the autophagy pathway are all important for defense against N. parisii infection. We also find that SCF ligase components like cul-6 promote defense against viral infection, where they have a more robust role than against N. parisii infection. This difference may be due to suppression of the host ubiquitylation system by N. parisii: when N. parisii is crippled by anti-microsporidia drugs, the host can more effectively target pathogen cells for ubiquitylation. Intriguingly, inhibition of the ubiquitin-proteasome system (UPS) increases expression of infection-upregulated SCF ligase components, indicating that a trigger for transcriptional response to intracellular infection by N. parisii and virus may be perturbation of the UPS. Altogether, our results demonstrate an in vivo role for ubiquitin-mediated defense against microsporidian and viral infections in C. elegans. [ABSTRACT FROM AUTHOR]

β-adrenergic receptor (βAR)-mediated transactivation of epidermal growth factor receptor (EGFR) has been shown to promote cardioprotection in a mouse model of heart failure and we recently showed that this mechanism leads to enhanced cell survival in part via regulation of apoptotic transcript expression in isolated primary rat neonatal cardiomyocytes. Thus, we hypothesized that this process could regulate cardiac transcript expression in vivo. To comprehensively assess cardiac transcript alterations in response to acute βAR-dependent EGFR transactivation, we performed whole transcriptome analysis of hearts from C57BL/6 mice given i.p. injections of the βAR agonist isoproterenol in the presence or absence of the EGFR antagonist gefitinib for 1 hour. Total cardiac RNA from each treatment group underwent transcriptome analysis, revealing a substantial number of transcripts regulated by each treatment. Gefitinib alone significantly altered the expression of 405 transcripts, while isoproterenol either alone or in conjunction with gefitinib significantly altered 493 and 698 distinct transcripts, respectively. Further statistical analysis was performed, confirming 473 transcripts whose regulation by isoproterenol were significantly altered by gefitinib (isoproterenol-induced up/downregulation antagonized/promoted by gefinitib), including several known to be involved in the regulation of numerous processes including cell death and survival. Thus, βAR-dependent regulation of cardiac transcript expression in vivo can be modulated by the EGFR antagonist gefitinib. [ABSTRACT FROM AUTHOR]

Background: Insects evolved ingenious adaptations to use extraordinary food sources. Particularly, the diet of herbivores enriched with noxious plant secondary metabolites requires detoxification mechanisms. Sequestration, which involves the uptake, transfer, and concentration of occasionally modified phytochemicals into specialized tissues or hemolymph, is one of the most successful detoxification strategies found in most insect orders. Due to the ability of ATP-binding cassette (ABC) carriers to transport a wide range of molecules including phytochemicals and xenobiotics, it is highly likely that they play a role in this sequestration process. To shed light on the role of ABC proteins in sequestration, we describe an inventory of putative ABC transporters in various tissues in the sequestering juvenile poplar leaf beetle, Chrysomela populi. Results: In the transcriptome of C. populi, we predicted 65 ABC transporters. To link the proteins with a possible function, we performed comparative phylogenetic analyses with ABC transporters of other insects and of humans. While tissue-specific profiling of each ABC transporter subfamily suggests that ABCB, C and G influence the plant metabolite absorption in the gut, ABCC with 14 members is the preferred subfamily responsible for the excretion of these metabolites via Malpighian tubules. Moreover, salicin, which is sequestered from poplar plants, is translocated into the defensive glands for further deterrent production. In these glands and among all identified ABC transporters, an exceptionally high transcript level was observed only for Cpabc35 (Cpmrp). RNAi revealed the deficiency of other ABC pumps to compensate the function of CpABC35, demonstrating its key role during sequestration. Conclusion: We provide the first comprehensive phylogenetic study of the ABC family in a phytophagous beetle species. RNA-seq data from different larval tissues propose the importance of ABC pumps to achieve a homeostasis of plant-derived compounds and offer a basis for future analyses of their physiological function in sequestration processes. [ABSTRACT FROM AUTHOR]

There is growing evidence for nongenetic effects of maternal experience on offspring. For example, previous studies have shown that female threespined stickleback fish (Gasterosteus aculeatus) exposed to predation risk produce offspring with altered behavior, metabolism and stress physiology. Here, we investigate the effect of maternal exposure to predation risk on the embryonic transcriptome in sticklebacks. Using RNA-sequencing we compared genome-wide transcription in three day post-fertilization embryos of predator-exposed and control mothers. There were hundreds of differentially expressed transcripts between embryos of predator-exposed mothers and embryos of control mothers including several non-coding RNAs. Gene Ontology analysis revealed biological pathways involved in metabolism, epigenetic inheritance, and neural proliferation and differentiation that differed between treatments. Interestingly, predation risk is associated with an accelerated life history in many vertebrates, and several of the genes and biological pathways that were identified in this study suggest that maternal exposure to predation risk accelerates the timing of embryonic development. Consistent with this hypothesis, embryos of predator-exposed mothers were larger than embryos of control mothers. These findings point to some of the molecular mechanisms that might underlie maternal effects. [ABSTRACT FROM AUTHOR]

Background: Parasite biology, by its very nature, cannot be understood without integrating it with that of the host, nor can the host response be adequately explained without considering the activity of the parasite. However, due to experimental limitations, molecular studies of parasite-host systems have been predominantly one-sided investigations focusing on either of the partners involved. Here, we conducted a dual RNA-seq time course analysis of filarial worm parasite and host mosquito to better understand the parasite processes underlying development in and interaction with the host tissue, from the establishment of infection to the development of infective-stage larva. Methodology/Principal Findings: Using the Brugia malayi–Aedes aegypti system, we report parasite gene transcription dynamics, which exhibited a highly ordered developmental program consisting of a series of cyclical and state-transitioning temporal patterns. In addition, we contextualized these parasite data in relation to the concurrent dynamics of the host transcriptome. Comparative analyses using uninfected tissues and different host strains revealed the influence of parasite development on host gene transcription as well as the influence of the host environment on parasite gene transcription. We also critically evaluated the life-cycle transcriptome of B. malayi by comparing developmental stages in the mosquito relative to those in the mammalian host, providing insight into gene expression changes underpinning the mosquito-borne parasitic lifestyle of this heteroxenous parasite. Conclusions/Significance: The data presented herein provide the research community with information to design wet lab experiments and select candidates for future study to more fully dissect the whole set of molecular interactions of both organisms in this mosquito-filarial worm symbiotic relationship. Furthermore, characterization of the transcriptional program over the complete life cycle of the parasite, including stages within the mosquito, could help devise novel targets for control strategies. [ABSTRACT FROM AUTHOR]

It has long been recognized that oncogenic viruses often integrate close to common fragile sites. The papillomavirus E2 protein, in complex with BRD4, tethers the viral genome to host chromatin to ensure persistent replication. Here, we map these targets to a number of large regions of the human genome and name them Persistent E2 and BRD4-Broad Localized Enrichments of Chromatin or PEB-BLOCs. PEB-BLOCs frequently contain deletions, have increased rates of asynchronous DNA replication, and are associated with many known common fragile sites. Cell specific fragile sites were mapped in human C-33 cervical cells by FANCD2 ChIP-chip, confirming the association with PEB-BLOCs. HPV-infected cells amplify viral DNA in nuclear replication foci and we show that these form adjacent to PEB-BLOCs. We propose that HPV replication, which hijacks host DNA damage responses, occurs adjacent to highly susceptible fragile sites, greatly increasing the chances of integration here, as is found in HPV-associated cancers. [ABSTRACT FROM AUTHOR]

Although sex chromosome meiotic drive has been observed in a variety of species for over 50 years, the genes causing drive are only known in a few cases, and none of these cases cause distorted sex-ratios in nature. In stalk-eyed flies (Teleopsis dalmanni), driving X chromosomes are commonly found at frequencies approaching 30% in the wild, but the genetic basis of drive has remained elusive due to reduced recombination between driving and non-driving X chromosomes. Here, we used RNAseq to identify transcripts that are differentially expressed between males carrying either a driving X (XSR) or a standard X chromosome (XST), and found hundreds of these, the majority of which are X-linked. Drive-associated transcripts show increased levels of sequence divergence (dN/dS) compared to a control set, and are predominantly expressed either in testes or in the gonads of both sexes. Finally, we confirmed that XSR and XST are highly divergent by estimating sequence differentiation between the RNAseq pools. We found that X-linked transcripts were often strongly differentiated (whereas most autosomal transcripts were not), supporting the presence of a relatively large region of recombination suppression on XSR presumably caused by one or more inversions. We have identified a group of genes that are good candidates for further study into the causes and consequences of sex-chromosome drive, and demonstrated that meiotic drive has had a profound effect on sequence evolution and gene expression of X-linked genes in this species. [ABSTRACT FROM AUTHOR]

Mapping the polymorphisms responsible for variation in gene expression, known as Expression Quantitative Trait Loci (eQTL), is a common strategy for investigating the molecular basis of disease. Despite numerous eQTL studies, the relationship between the explanatory power of variants on gene expression versus their power to explain ultimate phenotypes remains to be clarified. We addressed this question using four naturally occurring Quantitative Trait Nucleotides (QTN) in three transcription factors that affect sporulation efficiency in wild strains of the yeast, Saccharomyces cerevisiae. We compared the ability of these QTN to explain the variation in both gene expression and sporulation efficiency. We find that the amount of gene expression variation explained by the sporulation QTN is not predictive of the amount of phenotypic variation explained. The QTN are responsible for 98% of the phenotypic variation in our strains but the median gene expression variation explained is only 49%. The alleles that are responsible for most of the variation in sporulation efficiency do not explain most of the variation in gene expression. The balance between the main effects and gene-gene interactions on gene expression variation is not the same as on sporulation efficiency. Finally, we show that nucleotide variants in the same transcription factor explain the expression variation of different sets of target genes depending on whether the variant alters the level or activity of the transcription factor. Our results suggest that a subset of gene expression changes may be more predictive of ultimate phenotypes than the number of genes affected or the total fraction of variation in gene expression variation explained by causative variants, and that the downstream phenotype is buffered against variation in the gene expression network. [ABSTRACT FROM AUTHOR]

Personal exome and genome sequencing provides access to loss-of-function and rare deleterious alleles whose interpretation is expected to provide insight into individual disease burden. However, for each allele, accurate interpretation of its effect will depend on both its penetrance and the trait's expressivity. In this regard, an important factor that can modify the effect of a pathogenic coding allele is its level of expression; a factor which itself characteristically changes across tissues. To better inform the degree to which pathogenic alleles can be modified by expression level across multiple tissues, we have conducted exome, RNA and deep, targeted allele-specific expression (ASE) sequencing in ten tissues obtained from a single individual. By combining such data, we report the impact of rare and common loss-of-function variants on allelic expression exposing stronger allelic bias for rare stop-gain variants and informing the extent to which rare deleterious coding alleles are consistently expressed across tissues. This study demonstrates the potential importance of transcriptome data to the interpretation of pathogenic protein-coding variants. [ABSTRACT FROM AUTHOR]

Although patients with Glioblastoma multiforme (GBM) have grave prognosis, significant variability in patient outcome is observed. The objective of this study is to identify a molecular signature for GBM prognosis. We subjected 355 mRNA and microRNA expression profiles to elastic net-regulated Cox regression for identification of an integrated RNA signature for GBM prognosis. A prognostic index (PI) was generated for patient stratification. Survival comparison was conducted by Kaplan-Meier method and a general multivariate Cox regression procedure was applied to evaluate the independence of the PI. The abilities and efficiencies of signatures to predict GBM patient outcome was assessed and compared by the area under the curve (AUC) of the receiver-operator characteristic (ROC). An integrated RNA prognostic signature consisted by 4 protective mRNAs, 12 risky mRNAs, and 1 risky microRNA was identified. Decreased survival was associated with being in the high-risk group (hazard ratio = 2.864, P<0.0001). The prognostic value of the integrated signature was validated in five independent GBM expression datasets (n = 201, hazard ratio = 2.453, P<0.0001). The PI outperformed the known clinical factors, mRNA-only, and miRNA-only prognostic signatures for GBM prognosis (area under the ROC curve for the integrated RNA, mRNA-only, and miRNA-only signatures were 0.828, 0.742, and 0.757 at 3 years of overall survival, respectively, P<0.0001 by permutation test). We describe the first, to our knowledge, robust transcriptome-based integrated RNA signature that improves the current GBM prognosis based on clinical variables, mRNA-only, and miRNA-only signatures. [ABSTRACT FROM AUTHOR]

Background: Chronic inflammatory disorders have been increasing in incidence over the past decades following geographical patterns of industrialization. Fetal exposure to maternal inflammation may alter organ functions and the offspring's disease risk. We studied the development of genetically-driven ileitis and colitis in response to maternal inflammation using mouse models. Methods: Disease susceptible (TnfΔARE/+ and IL10−/−) and disease-free (Tnf+/+ and IL10−/+) offspring were raised in inflamed and non-inflamed dams. Ileal, caecal and colonic pathology was evaluated in the offspring at 8 or 12 weeks of age. Ly6G-positive cells in inflamed sections from the distal ileum and distal colon were analysed by immunofluorescence microscopy. Gene expression of pro-inflammatory cytokines was measured in whole tissue specimens by quantitative PCR. Microarray analyses were performed on laser microdissected intestinal epithelium. Caecal bacterial communities were assessed by Illumina sequencing of 16S rRNA amplicons. Results: Disease severity, the number of infiltrated neutrophils as well as Tnf and Il12p40 mRNA expression were independent of maternal inflammation in the offspring of mouse models for ileitis (TnfΔARE/+) and colitis (IL10−/−). Although TNF-driven maternal inflammation regulated 2,174 (wild type) and 3,345 (TnfΔARE/+) genes in the fetal epithelium, prenatal gene expression patterns were completely overwritten after birth. In addition, co-housing experiments revealed no change in phylogenetic diversity of the offspring's caecal microbiota in response to maternal inflammation. This is independent of the offspring's genotype before and after the onset of tissue pathology. Conclusions: Disease risk and activity in mouse models of chronic ileitis and colitis was independent of the fetal exposure to maternal inflammation. Likewise, maternal inflammation did not alter the diversity and composition of offspring's caecal microbiota, clearly demonstrating that changes of the gene expression program in the fetal gut epithelium were not relevant for the development of chronic inflammatory disorders in the gut. [ABSTRACT FROM AUTHOR]

Mannitol plays a crucial role in brown algae, acting as carbon storage, organic osmolytes and antioxidant. Transcriptomic analysis of Saccharina japonica revealed that the relative genes involved in the mannitol cycle are existent. Full-length sequence of mannitol-2-dehydrogenase (M2DH) gene was obtained, with one open reading frame of 2,007 bp which encodes 668 amino acids. Cis-regulatory elements for response to methyl jasmonic acid, light and drought existed in the 5′-upstream region. Phylogenetic analysis indicated that SjM2DH has an ancient prokaryotic origin, and is probably acquired by horizontal gene transfer event. Multiple alignment and spatial structure prediction displayed a series of conserved functional residues, motifs and domains, which favored that SjM2DH belongs to the polyol-specific long-chain dehydrogenases/reductase (PSLDR) family. Expressional profiles of SjM2DH in the juvenile sporophytes showed that it was influenced by saline, oxidative and desiccative factors. SjM2DH was over-expressed in Escherichia coli, and the cell-free extracts with recombinant SjM2DH displayed high activity on D-fructose reduction reaction. The analysis on SjM2DH gene structure and biochemical parameters reached a consensus that activity of SjM2DH is NADH-dependent and metal ion-independent. The characterization of SjM2DH showed that M2DH is a new member of PSLDR family and play an important role in mannitol metabolism in S. japonica. [ABSTRACT FROM AUTHOR]

Background: Salt stress interferes with plant growth and production. Plants have evolved a series of molecular and morphological adaptations to cope with this abiotic stress, and overexpression of salt response genes reportedly enhances the productivity of various crops. However, little is known about the salt responsive genes in the energy plant physic nut (Jatropha curcas L.). Thus, excavate salt responsive genes in this plant are informative in uncovering the molecular mechanisms for the salt response in physic nut. Methodology/Principal Findings: We applied next-generation Illumina sequencing technology to analyze global gene expression profiles of physic nut plants (roots and leaves) 2 hours, 2 days and 7 days after the onset of salt stress. A total of 1,504 and 1,115 genes were significantly up and down-regulated in roots and leaves, respectively, under salt stress condition. Gene ontology (GO) analysis of physiological process revealed that, in the physic nut, many “biological processes” were affected by salt stress, particular those categories belong to “metabolic process”, such as “primary metabolism process”, “cellular metabolism process” and “macromolecule metabolism process”. The gene expression profiles indicated that the associated genes were responsible for ABA and ethylene signaling, osmotic regulation, the reactive oxygen species scavenging system and the cell structure in physic nut. Conclusions/Significance: The major regulated genes detected in this transcriptomic data were related to trehalose synthesis and cell wall structure modification in roots, while related to raffinose synthesis and reactive oxygen scavenger in leaves. The current study shows a comprehensive gene expression profile of physic nut under salt stress. The differential expression genes detected in this study allows the underling the salt responsive mechanism in physic nut with the aim of improving its salt resistance in the future. [ABSTRACT FROM AUTHOR]

We investigated the role of tumor copy number (CN)–altered genome (CN-AG) in the carcinogenesis of cervical cancer (CC), especially its effect on gene expression, biological processes, and patient survival. Fifty-nine human papillomavirus 16 (HPV16)-positive CCs were investigated with microarrays–31 for mapping CN-AG and 55 for global gene expression, with 27 CCs in common. Five-year survival was investigated in 55 patients. Deletions and amplifications >2.5 Mb were defined as CN alterations. The %CN-AG varied from 0 to 32.2% (mean = 8.1±8.9). Tumors were classified as low (mean = 0.5±0.6, n = 11), medium (mean = 5.4±2.4, n = 10), or high (mean = 19.2±6.6, n = 10) CN. The highest %CN-AG was found in 3q, which contributed an average of 55% of all CN alterations. Genome-wide, only 5.3% of CN-altered genes were deregulated directly by gene dosage. In contrast, the rate in fully duplicated 3q was twice as high. Amplification of 3q explained 23.2% of deregulated genes in whole tumors (r2 = 0.232, p = 0.006; analysis of variance), including genes located in 3q and other chromosomes. A total of 862 genes were deregulated exclusively in high-CN tumors, but only 22.9% were CN altered. This suggests that the remaining genes are not deregulated directly by gene dosage, but by mechanisms induced in trans by CN-altered genes. Anaphase-promoting complex/cyclosome (APC/C)-dependent proteasome proteolysis, glycolysis, and apoptosis were upregulated, whereas cell adhesion and angiogenesis were downregulated exclusively in high-CN tumors. The high %CN-AG and upregulated gene expression profile of APC/C-dependent proteasome proteolysis were associated with poor patient survival (p<0.05, log-rank test). Along with glycolysis, they were linearly associated with FIGO stage (r>0.38, p<0.01, Spearman test). Therefore, inhibition of APC/C-dependent proteasome proteolysis and glycolysis could be useful for CC treatment. However, whether they are indispensable for tumor growth remains to be demonstrated. [ABSTRACT FROM AUTHOR]

Aberrant methylation of specific CpG sites at the promoter is widely responsible for genesis and development of various cancer types. Even though the microarray-based methylome analyzing techniques have contributed to the elucidation of the methylation change at the genome-wide level, the identification of key methylation markers or top regulatory networks appearing common in highly incident cancers through comparison analysis is still limited. In this study, we in silico performed the genome-wide methylation analysis on each 10 sets of normal and cancer pairs of five tissues: breast, colon, liver, lung, and stomach. The methylation array covers 27,578 CpG sites, corresponding to 14,495 genes, and significantly hypermethylated or hypomethylated genes in the cancer were collected (FDR adjusted p-value <0.05; methylation difference >0.3). Analysis of the dataset confirmed the methylation of previously known methylation markers and further identified novel methylation markers, such as GPX2, CLDN15, and KL. Cluster analysis using the methylome dataset resulted in a diagram with a bipartite mode distinguishing cancer cells from normal cells regardless of tissue types. The analysis further revealed that breast cancer was closest with lung cancer, whereas it was farthest from colon cancer. Pathway analysis identified that either the “cancer” related network or the “cancer” related bio-function appeared as the highest confidence in all the five cancers, whereas each cancer type represents its tissue-specific gene sets. Our results contribute toward understanding the essential abnormal epigenetic pathways involved in carcinogenesis. Further, the novel methylation markers could be applied to establish markers for cancer prognosis. [ABSTRACT FROM AUTHOR]

The edible mushroom Volvariella volvacea is an important crop in Southeast Asia and is predominantly harvested in the egg stage. One of the main factors that negatively affect its yield and value is the rapid transition from the egg to the elongation stage, which has a decreased commodity value and shelf life. To improve our understanding of the changes during stipe development and the transition from egg to elongation stage in particular, we analyzed gene transcription in stipe tissue of V. volvacea using 3′-tag based digital expression profiling. Stipe development turned out to be fairly complex with high numbers of expressed genes, and regulation of stage differences is mediated mainly by changes in expression levels of genes, rather than on/off modulation. Most explicit is the strong up-regulation of cell division from button to egg, and the very strong down-regulation hereof from egg to elongation, that continues in the maturation stage. Button and egg share cell division as means of growth, followed by a major developmental shift towards rapid stipe elongation based on cell extension as demonstrated by inactivation of cell division throughout elongation and maturation. Examination of regulatory genes up-regulated from egg to elongation identified three potential high upstream regulators for this switch. The new insights in stipe dynamics, together with a series of new target genes, will provide a sound base for further studies on the developmental mechanisms of mushroom stipes and the switch from egg to elongation in V. volvacea in particular. [ABSTRACT FROM AUTHOR]

Microarray analyses of gene knockouts have traditionally focused on the identification of genes whose mean expression is different in knockout and wild-type mice. However, recent work suggests that changes in the variability of gene expression can have important phenotypic consequences as well. Here, in an unbiased sample of publicly available microarray data on gene expression in various knockouts, highly significant differences from wild-type (either increases or decreases) are noted in the gene expression coefficients of variation (CVs) of virtually every knockout considered. Examination of the distribution of gene-by-gene CV differences indicates that these findings are not attributable to a few outlier genes, but rather to broadly increased or decreased CV in the various knockouts over all the (tens of thousands of) transcripts assayed. These global differences in variability may reflect either authentic biological effects of the knockouts or merely experimental inconsistencies. However, regardless of the underlying explanation, the variability differences are of importance as they will influence both the statistical detection of gene expression changes and, potentially, the knockout phenotype itself. [ABSTRACT FROM AUTHOR]

Regulation of transcription is one of the mechanisms involved in animal development, directing changes in patterning and cell fate specification. Large temporal data series, based on microarrays across the life cycle of the fly Drosophila melanogaster, revealed the existence of groups of genes which expression increases or decreases temporally correlated during the life cycle. These groups of genes are enriched in different biological functions. Here, instead of searching for temporal coincidence in gene expression using the entire genome expression data, we searched for temporal coincidence in gene expression only within predefined catalogues of functionally related genes and investigated whether a catalogue's expression profile can be used to generate larger catalogues, enriched in genes necessary for the same function. We analyzed the expression profiles from genes already associated with early neurodevelopment and late neurodifferentiation, at embryonic stages 16 and 17 of Drosophila life cycle. We hypothesized that during this interval we would find global downregulation of genes important for early neuronal development together with global upregulation of genes necessary for the final differentiation of neurons. Our results were consistent with this hypothesis. We then investigated if the expression profile of gene catalogues representing particular processes of neural development matched the temporal sequence along which these processes occur. The profiles of genes involved in patterning, neurogenesis, axogenesis or synaptic transmission matched the prediction, with largest transcript values at the time when the corresponding biological process takes place in the embryo. Furthermore, we obtained catalogues enriched in genes involved in temporally matching functions by performing a genome-wide systematic search for genes with their highest expression levels at the corresponding embryonic intervals. These findings imply the use of gene expression data in combination with known biological information to predict the involvement of functionally uncharacterized genes in particular biological events. [ABSTRACT FROM AUTHOR]

Defining the aggressiveness and growth rate of a malignant cell population is a key step in the clinical approach to treating tumor disease. The correct grading of breast cancer (BC) is a fundamental part in determining the appropriate treatment. Biological variables can make it difficult to elucidate the mechanisms underlying BC development. To identify potential markers that can be used for BC classification, we analyzed mRNAs expression profiles, gene copy numbers, microRNAs expression and their association with tumor grade in BC microarray-derived datasets. From mRNA expression results, we found that grade 2 BC is most likely a mixture of grade 1 and grade 3 that have been misclassified, being described by the gene signature of either grade 1 or grade 3. We assessed the potential of the new approach of integrating mRNA expression profile, copy number alterations, and microRNA expression levels to select a limited number of genomic BC biomarkers. The combination of mRNA profile analysis and copy number data with microRNA expression levels led to the identification of two gene signatures of 42 and 4 altered genes (FOXM1, KPNA4, H2AFV and DDX19A) respectively, the latter obtained through a meta-analytical procedure. The 42-based gene signature identifies 4 classes of up- or down-regulated microRNAs (17 microRNAs) and of their 17 target mRNA, and the 4-based genes signature identified 4 microRNAs (Hsa-miR-320d, Hsa-miR-139-5p, Hsa-miR-567 and Hsa-let-7c). These results are discussed from a biological point of view with respect to pathological features of BC. Our identified mRNAs and microRNAs were validated as prognostic factors of BC disease progression, and could potentially facilitate the implementation of assays for laboratory validation, due to their reduced number. [ABSTRACT FROM AUTHOR]

Sponges and bacteria have lived together in complex consortia for 700 million years. As filter feeders, sponges prey on bacteria. Nevertheless, some bacteria are associated with sponges in symbiotic relationships. To enable this association, sponges and bacteria are likely to have developed molecular communication systems. These may include molecules such as N-acyl-L-homoserine lactones, produced by Gram-negative bacteria also within sponges. In this study, we examined the role of N-3-oxododecanoyl-L-homoserine lactone (3-oxo-C12-HSL) on the expression of immune and apoptotic genes of the host sponge Suberites domuncula. This molecule seemed to inhibit the sponge innate immune system through a decrease of the expression of genes coding for proteins sensing the bacterial membrane: a Toll-Like Receptor and a Toll-like Receptor Associated Factor 6 and for an anti-bacterial perforin-like molecule. The expression of the pro-apoptotic caspase-like 3/7 gene decreased as well, whereas the level of mRNA of anti-apoptotic genes Bcl-2 Homolog Proteins did not change. Then, we demonstrated the differential expression of proteins in presence of this 3-oxo-C12-HSL using 3D sponge cell cultures. Proteins involved in the first steps of the endocytosis process were highlighted using the 2D electrophoresis protein separation and the MALDI-TOF/TOF protein characterization: α and β subunits of the lysosomal ATPase, a cognin, cofilins-related proteins and cytoskeleton proteins actin, α tubulin and α actinin. The genetic expression of some of these proteins was subsequently followed. We propose that the 3-oxo-C12-HSL may participate in the tolerance of the sponge apoptotic and immune systems towards the presence of bacteria. Besides, the sponge may sense the 3-oxo-C12-HSL as a molecular evidence of the bacterial presence and/or density in order to regulate the populations of symbiotic bacteria in the sponge. This study is the first report of a bacterial secreted molecule acting on sponge cells and regulating the symbiotic relationship. [ABSTRACT FROM AUTHOR]

Hypothetical proteins comprise roughly half of the predicted gene complement of Toxoplasma gondii and Plasmodium falciparum and represent the largest class of uniquely functioning proteins in these parasites. Following the idea that functional relationships can be informed by the timing of gene expression, we devised a strategy to identify the core set of apicomplexan cell division cycling genes with important roles in parasite division, which includes many uncharacterized proteins. We assembled an expanded list of orthologs from the T. gondii and P. falciparum genome sequences (2781 putative orthologs), compared their mRNA profiles during synchronous replication, and sorted the resulting set of dual cell cycle regulated orthologs (744 total) into protein pairs conserved across many eukaryotic families versus those unique to the Apicomplexa. The analysis identified more than 100 ortholog gene pairs with unknown function in T. gondii and P. falciparum that displayed co-conserved mRNA abundance, dynamics of cyclical expression and similar peak timing that spanned the complete division cycle in each parasite. The unknown cyclical mRNAs encoded a diverse set of proteins with a wide range of mass and showed a remarkable conservation in the internal organization of ordered versus disordered structural domains. A representative sample of cyclical unknown genes (16 total) was epitope tagged in T. gondii tachyzoites yielding the discovery of new protein constituents of the parasite inner membrane complex, key mitotic structures and invasion organelles. These results demonstrate the utility of using gene expression timing and dynamic profile to identify proteins with unique roles in Apicomplexa biology. [ABSTRACT FROM AUTHOR]

The paper mulberry is one of the multifunctional tree species in agroforestry systems and is also commonly utilized in traditional medicine in China and other Asian countries. However, little is known about its molecular genetics, which hinders research on and exploitation of this valuable resource. To discern the correlation between gene expression and the essential properties of the paper mulberry, we performed a transcriptomics analysis, assembling a total of 37,725 unigenes from 54,638,676 reads generated by RNA-seq. Among these, 22,692 unigenes showed greater than 60% similarity with genes from other species. The lengths of 13,566 annotated unigenes were longer than 1,000 bp. Functional clustering analysis with COG (Cluster of Orthologous Groups) revealed that 17,184 unigenes are primarily involved in transcription, translation, signal transduction, carbohydrate metabolism, secondary metabolism, and energy metabolism. GO (Gene Ontology) annotation suggests enrichment of genes encoding antioxidant activity, transporter activity, biosynthesis, metabolism and stress response, with a total of 30,659 unigenes falling in these categories. KEGG (Kyoto Encyclopedia of Genes and Genomes) metabolic pathway analysis showed that 7,199 unigenes are associated with 119 metabolic pathways. In addition to the basic metabolism, these genes are enriched for plant pathogen interaction, flavonoid metabolism and other secondary metabolic processes. Furthermore, differences in the transcriptomes of leaf, stem and root tissues were analyzed and 7,233 specifically expressed unigenes were identified. This global expression analysis provided novel insights about the molecular mechanisms of the biosynthesis of flavonoid, lignin and cellulose, as well as on the response to biotic and abiotic stresses including the remediation of contaminated soil by the paper mulberry. [ABSTRACT FROM AUTHOR]

Background: Crassostrea gigas accumulates diarrheic shellfish toxins (DSP) associated to Prorocentrum lima of which Okadaic acid (OA) causes specific inhibitions of serine and threonine phosphatases 1 and 2A. Its toxic effects have been extensively reported in bivalve mollusks at cellular and physiological levels, but genomic approaches have been scarcely studied. Methodology/Principal Findings: Acute and sub-chronic exposure effects of P. lima were investigated on farmed juvenile C. gigas (3–5 mm). The Pacific oysters were fed with three dinoflagellate concentrations: 0.3, 3, and 30×103 cells mL−1 along with a nontoxic control diet of Isochrysis galbana. The effects of P. lima on C. gigas were followed by analyzing expression levels of a total of four genes, three involved in cell cycle regulation and one in immune response by polymerase chain reaction and real time quantitative PCR, where changes in time and cell concentration were found. The highest expression levels were found in oysters fed 3×103 cells mL−1 at 168 h for the cycle regulator p21 protein (9 fold), chromatin assembly factor 1 p55 subunit (8 fold), elongation factor 2 (2 fold), and lipopolysaccharide/β-1, 3 glucan binding protein (13 fold above base line). Additionally, the transcript level of all the genes decreased in oysters fed wich the mixed diet 30×103 cells mL−1 of dinoflagellate after 72 h and was lowest in the chromatin assembly factor 1 p55 subunit (0.9 fold below baseline). Conclusions: On C. gigas the whole cell ingestion of P lima caused a clear mRNA modulation expression of the genes involved in cell cycle regulation and immune system. Over-expression could be related to DNA damage, disturbances in cell cycle continuity, probably a genotoxic effect, as well as an activation of its innate immune system as first line of defense. [ABSTRACT FROM AUTHOR]

In transcriptomics research, design for experimentation by carefully considering biological, technological, practical and statistical aspects is very important, because the experimental design space is essentially limitless. Usually, the ranges of variable biological parameters of the design space are based on common practices and in turn on phenotypic endpoints. However, specific sub-cellular processes might only be partially reflected by phenotypic endpoints or outside the associated parameter range. Here, we provide a generic protocol for range finding in design for transcriptomics experimentation based on small-scale gene-expression experiments to help in the search for the right location in the design space by analyzing the activity of already known genes of relevant molecular mechanisms. Two examples illustrate the applicability: in-vitro UV-C exposure of mouse embryonic fibroblasts and in-vivo UV-B exposure of mouse skin. Our pragmatic approach is based on: framing a specific biological question and associated gene-set, performing a wide-ranged experiment without replication, eliminating potentially non-relevant genes, and determining the experimental ‘sweet spot’ by gene-set enrichment plus dose-response correlation analysis. Examination of many cellular processes that are related to UV response, such as DNA repair and cell-cycle arrest, revealed that basically each cellular (sub-) process is active at its own specific spot(s) in the experimental design space. Hence, the use of range finding, based on an affordable protocol like this, enables researchers to conveniently identify the ‘sweet spot’ for their cellular process of interest in an experimental design space and might have far-reaching implications for experimental standardization. [ABSTRACT FROM AUTHOR]

The RpoS/σS sigma subunit of RNA polymerase (RNAP) controls a global adaptive response that allows many Gram-negative bacteria to survive starvation and various stresses. σS also contributes to biofilm formation and virulence of the food-borne pathogen Salmonella enterica serovar Typhimurium (S. Typhimurium). In this study, we used directional RNA-sequencing and complementary assays to explore the σS-dependent transcriptome of S. Typhimurium during late stationary phase in rich medium. This study confirms the large regulatory scope of σS and provides insights into the physiological functions of σS in Salmonella. Extensive regulation by σS of genes involved in metabolism and membrane composition, and down-regulation of the respiratory chain functions, were important features of the σS effects on gene transcription that might confer fitness advantages to bacterial cells and/or populations under starving conditions. As an example, we show that arginine catabolism confers a competitive fitness advantage in stationary phase. This study also provides a firm basis for future studies to address molecular mechanisms of indirect regulation of gene expression by σS. Importantly, the σS-controlled downstream network includes small RNAs that might endow σS with post-transcriptional regulatory functions. Of these, four (RyhB-1/RyhB-2, SdsR, SraL) were known to be controlled by σS and deletion of the sdsR locus had a competitive fitness cost in stationary phase. The σS-dependent control of seven additional sRNAs was confirmed in Northern experiments. These findings will inspire future studies to investigate molecular mechanisms and the physiological impact of post-transcriptional regulation by σS. [ABSTRACT FROM AUTHOR]

Clostridium difficile is the main agent responsible for hospital acquired antibiotic associated diarrhoea. In recent years, epidemic strains have emerged causing more severe infections. Whilst C. difficile has two major virulence factors, toxins TcdA and TcdB, it is generally accepted that other virulence components of the bacterium contribute to disease. Previously, it has been suggested that flagella expression from pathogenic bacteria might be implicated in virulence. In a recent study, we observed an increased mortality in a gnotobiotic mouse model when animals were colonized with an isogenic fliC mutant constructed in the PCR-ribotype 027 (B1/NAP1) strain R20291, while animals survived when colonized by the parental strain or after colonization by other high-toxin-producing C. difficile strains. To understand the reasons for this increased virulence, we compared the global gene expression profiles between the fliC-R20291 mutant and its parental strain using an in vitro and in vivo transcriptomic approach. The latter made use of the gnotobiotic mouse model. Interestingly, in the fliC mutant, we observed considerable up-regulation of genes involved in mobility, membrane transport systems (PTS, ABC transporters), carbon metabolism, known virulence factors and sporulation. A smaller but significant up-regulation of genes involved in cell growth, fermentation, metabolism, stress and antibiotic resistance was also apparent. All of these genes may be associated with the increased virulence of the fliC-R20921 mutant. We confirmed that the fliC mutation is solely responsible for the observed changes in gene expression in the mutant strain since expression profiles were restored to that of the wild-type strain in the fliC-complemented strain. Thus, the absence of FliC is directly or indirectly involved in the high mortality observed in the fliC mutant infected animals. Therefore, we provide the first evidence that when the major structural component of the flagellum is neutralized, deregulation of gene expression can occur during infection. [ABSTRACT FROM AUTHOR]

Background: Triticum monococcum (2n) is a close ancestor of T. urartu, the A-genome progenitor of cultivated hexaploid wheat, and is therefore a useful model for the study of components regulating photomorphogenesis in diploid wheat. In order to develop genetic and genomic resources for such a study, we constructed genome-wide transcriptomes of two Triticum monococcum subspecies, the wild winter wheat T. monococcum ssp. aegilopoides (accession G3116) and the domesticated spring wheat T. monococcum ssp. monococcum (accession DV92) by generating de novo assemblies of RNA-Seq data derived from both etiolated and green seedlings. Principal Findings: The de novo transcriptome assemblies of DV92 and G3116 represent 120,911 and 117,969 transcripts, respectively. We successfully mapped ∼90% of these transcripts from each accession to barley and ∼95% of the transcripts to T. urartu genomes. However, only ∼77% transcripts mapped to the annotated barley genes and ∼85% transcripts mapped to the annotated T. urartu genes. Differential gene expression analyses revealed 22% more light up-regulated and 35% more light down-regulated transcripts in the G3116 transcriptome compared to DV92. The DV92 and G3116 mRNA sequence reads aligned against the reference barley genome led to the identification of ∼500,000 single nucleotide polymorphism (SNP) and ∼22,000 simple sequence repeat (SSR) sites. Conclusions: De novo transcriptome assemblies of two accessions of the diploid wheat T. monococcum provide new empirical transcriptome references for improving Triticeae genome annotations, and insights into transcriptional programming during photomorphogenesis. The SNP and SSR sites identified in our analysis provide additional resources for the development of molecular markers. [ABSTRACT FROM AUTHOR]

Animals are frequently used as model systems for determination of safety and efficacy in pharmaceutical research and development. However, significant quantitative and qualitative differences exist between humans and the animal models used in research. This is as a result of genetic variation between human and the laboratory animal. Therefore the development of a system that would allow the assessment of all molecular differences between species after drug exposure would have a significant impact on drug evaluation for toxicity and efficacy. Here we describe a cross-species microarray methodology that identifies and selects orthologous probes after cross-species sequence comparison to develop an orthologous cross-species gene expression analysis tool. The assumptions made by the use of this orthologous gene expression strategy for cross-species extrapolation is that; conserved changes in gene expression equate to conserved pharmacodynamic endpoints. This assumption is supported by the fact that evolution and selection have maintained the structure and function of many biochemical pathways over time, resulting in the conservation of many important processes. We demonstrate this cross-species methodology by investigating species specific differences of the peroxisome proliferator-activator receptor (PPAR) α response in rat and human. [ABSTRACT FROM AUTHOR]

Arsenic is a ubiquitous contaminant and a toxic metalloid which presents two main redox states in nature: arsenite [AsIII] and arsenate [AsV]. Arsenic resistance in Synechocystis sp. strain PCC 6803 is mediated by the arsBHC operon and two additional arsenate reductases encoded by the arsI1 and arsI2 genes. Here we describe the genome-wide responses to the presence of arsenate and arsenite in wild type and mutants in the arsenic resistance system. Both forms of arsenic produced similar responses in the wild type strain, including induction of several stress related genes and repression of energy generation processes. These responses were transient in the wild type strain but maintained in time in an arsB mutant strain, which lacks the arsenite transporter. In contrast, the responses observed in a strain lacking all arsenate reductases were somewhat different and included lower induction of genes involved in metal homeostasis and Fe-S cluster biogenesis, suggesting that these two processes are targeted by arsenite in the wild type strain. Finally, analysis of the arsR mutant strain revealed that ArsR seems to only control 5 genes in the genome. Furthermore, the arsR mutant strain exhibited hypersentivity to nickel, copper and cadmium and this phenotype was suppressed by mutation in arsB but not in arsC gene suggesting that overexpression of arsB is detrimental in the presence of these metals in the media. [ABSTRACT FROM AUTHOR]

Transcriptome analysis based on a de novo assembly of next generation RNA sequences is now performed routinely in many laboratories. The generated results, including contig sequences, quantification figures, functional annotations and variation discovery outputs are usually bulky and quite diverse. This article presents a user oriented storage and visualisation environment permitting to explore the data in a top-down manner, going from general graphical views to all possible details. The software package is based on biomart, easy to install and populate with local data. The software package is available under the GNU General Public License (GPL) at http://bioinfo.genotoul.fr/RNAbrowse. [ABSTRACT FROM AUTHOR]

Maize (Zea mays) is the most widely grown crop species in the world and a classical model organism for plant research. The completion of a high-quality reference genome sequence and the advent of high-throughput sequencing have greatly empowered re-sequencing studies in maize. In this study, plants of maize inbred line B73 descended from two different sets of seed material grown for several generations either in the field or in the greenhouse were found to show a different growth phenotype and ionome under phosphate starvation conditions and moreover a different responsiveness towards mycorrhizal fungi of the species Glomus intraradices (syn: Rhizophagus irregularis). Whole genome re-sequencing of individuals from both sets and comparison to the B73 reference sequence revealed three cryptic introgressions on chromosomes 1, 5 and 10 in the line grown in the greenhouse summing up to a total of 5,257 single-nucleotide polymorphisms (SNPs). Transcriptome sequencing of three individuals from each set lent further support to the location of the introgression intervals and confirmed them to be fixed in all sequenced individuals. Moreover, we identified >120 genes differentially expressed between the two B73 lines. We thus have found a nearly-isogenic line (NIL) of maize inbred line B73 that is characterized by an altered growth phenotype under phosphate starvation conditions and an improved responsiveness towards symbiosis with mycorrhizal fungi. Through next-generation sequencing of the genomes and transcriptomes we were able to delineate exact introgression intervals. Putative de novo mutations appeared approximately uniformly distributed along the ten maize chromosomes mainly representing G:C -> A:T transitions. The plant material described in this study will be a valuable tool both for functional studies of genes differentially expressed in both B73 lines and for research on growth behavior especially in response to symbiosis between maize and mycorrhizal fungi. [ABSTRACT FROM AUTHOR]

Growth factor independence 1b (GFI1B) is a DNA binding repressor of transcription with vital functions in hematopoiesis. Gfi1b-null embryos die at midgestation very likely due to defects in erythro- and megakaryopoiesis. To analyze the full functionality of Gfi1b, we used conditionally deficient mice that harbor floxed Gfi1b alleles and inducible (Mx-Cre, Cre-ERT) or erythroid specific (EpoR-Cre) Cre expressing transgenes. In contrast to the germline knockout, EpoR-Cre mediated erythroid specific ablation of Gfi1b allows full gestation, but causes perinatal lethality with very few mice surviving to adulthood. Both the embryonic deletion of Gfi1b by EpoR-Cre and the deletion in adult mice by Mx-Cre or Cre-ERT leads to reduced numbers of erythroid precursors, perturbed and delayed erythroid maturation, anemia and extramedullary erythropoiesis. Global expression analyses showed that the Hba-x, Hbb-bh1 and Hbb-y embryonic globin genes were upregulated in Gfi1b deficient TER119+ fetal liver cells over the gestation period from day 12.5–17.5 p.c. and an increased level of Hbb-bh1 and Hbb-y embryonic globin gene expression was even maintained in adult Gfi1b deficient mice. While the expression of Bcl11a, a regulator of embryonic globin expression was not affected by Gfi1b deficiency, the expression of Gata1 was reduced and the expression of Sox6, also involved in globin switch, was almost entirely lost when Gfi1b was absent. These findings establish Gfi1b as a regulator of embryonic globin expression and embryonic and adult erythroid maturation. [ABSTRACT FROM AUTHOR]

Background: Meat quality depends on skeletal muscle structure and metabolic properties. While most studies carried on pigs focus on the Longissimus muscle (LM) for fresh meat consumption, Semimembranosus (SM) is also of interest because of its importance for cooked ham production. Even if both muscles are classified as glycolytic muscles, they exhibit dissimilar myofiber composition and metabolic characteristics. The comparison of LM and SM transcriptome profiles undertaken in this study may thus clarify the biological events underlying their phenotypic differences which might influence several meat quality traits. Methodology/Principal Findings: Muscular transcriptome analyses were performed using a custom pig muscle microarray: the 15 K Genmascqchip. A total of 3823 genes were differentially expressed between the two muscles (Benjamini-Hochberg adjusted P value ≤0.05), out of which 1690 and 2133 were overrepresented in LM and SM respectively. The microarray data were validated using the expression level of seven differentially expressed genes quantified by real-time RT-PCR. A set of 1047 differentially expressed genes with a muscle fold change ratio above 1.5 was used for functional characterization. Functional annotation emphasized five main clusters associated to transcriptome muscle differences. These five clusters were related to energy metabolism, cell cycle, gene expression, anatomical structure development and signal transduction/immune response. Conclusions/Significance: This study revealed strong transcriptome differences between LM and SM. These results suggest that skeletal muscle discrepancies might arise essentially from different post-natal myogenic activities. [ABSTRACT FROM AUTHOR]

Epithelial ovarian cancer (EOC) is the most common gynecologic malignancy. To identify the micro-ribonucleic acids (miRNAs) expression profile in EOC tissues that may serve as a novel diagnostic biomarker for EOC detection, the expression of 1722 miRNAs from 15 normal ovarian tissue samples and 48 ovarian cancer samples was profiled by using a quantitative real-time polymerase chain reaction (qRT-PCR) assay. A ten-microRNA signature (hsa-miR-1271-5p, hsa-miR-574-3p, hsa-miR-182-5p, hsa-miR-183-5p, hsa-miR-96-5p, hsa-miR-15b-5p, hsa-miR-182-3p, hsa-miR-141-5p, hsa-miR-130b-5p, and hsa-miR-135b-3p) was identified to be able to distinguish human ovarian cancer tissues from normal tissues with 97% sensitivity and 92% specificity. Two miRNA clusters of miR183-96-183 (miR-96-5p, and miR-182, miR183) and miR200 (miR-141-5p, miR200a, b, c and miR429) are significantly up-regulated in ovarian cancer tissue samples compared to those of normal tissue samples, suggesting theses miRNAs may be involved in ovarian cancer development. [ABSTRACT FROM AUTHOR]

The plant cell wall performs a number of essential functions including providing shape to many different cell types and serving as a defense against potential pathogens. The net pattern mutation creates breaks in the seed coat of soybean (Glycine max) because of ruptured cell walls. Using RNA-Seq, we examined the seed coat transcriptome from three stages of immature seed development in two pairs of isolines with normal or defective seed coat phenotypes due to the net pattern. The genome-wide comparative study of the transcript profiles of these isolines revealed 364 differentially expressed genes in common between the two varieties that were further divided into different broad functional categories. Genes related to cell wall processes accounted for 19% of the differentially expressed genes in the middle developmental stage of 100–200 mg seed weight. Within this class, the cell wall proline-rich and glycine-rich protein genes were highly differentially expressed in both genetic backgrounds. Other genes that showed significant expression changes in each of the isoline pairs at the 100–200 mg seed weight stage were xylem serine proteinase, fasciclin-related genes, auxin and stress response related genes, TRANSPARENT TESTA 1 (TT1) and other transcription factors. The mutant appears to shift the timing of either the increase or decrease in the levels of some of the transcripts. The analysis of these data sets reveals the physiological changes that the seed coat undergoes during the formation of the breaks in the cell wall. [ABSTRACT FROM AUTHOR]

Olfactory perception is mediated by a large array of olfactory receptor genes. The human genome contains 851 olfactory receptor gene loci. More than 50% of the loci are annotated as nonfunctional due to frame-disrupting mutations. Furthermore haplotypic missense alleles can be nonfunctional resulting from substitution of key amino acids governing protein folding or interactions with signal transduction components. Beyond their role in odor recognition, functional olfactory receptors are also required for a proper targeting of olfactory neuron axons to their corresponding glomeruli in the olfactory bulb. Therefore, we anticipate that profiling of olfactory receptor gene expression in whole human olfactory mucosa and analysis in the human population of their expression should provide an opportunity to select the frequently expressed and potentially functional olfactory receptors in view of a systematic deorphanization. To address this issue, we designed a TaqMan Low Density Array (Applied Biosystems), containing probes for 356 predicted human olfactory receptor loci to investigate their expression in whole human olfactory mucosa tissues from 26 individuals (13 women, 13 men; aged from 39 to 81 years, with an average of 67±11 years for women and 63±12 years for men). Total RNA isolation, DNase treatment, RNA integrity evaluation and reverse transcription were performed for these 26 samples. Then 384 targeted genes (including endogenous control genes and reference genes specifically expressed in olfactory epithelium for normalization purpose) were analyzed using the same real-time reverse transcription PCR platform. On average, the expression of 273 human olfactory receptor genes was observed in the 26 selected whole human olfactory mucosa analyzed, of which 90 were expressed in all 26 individuals. Most of the olfactory receptors deorphanized to date on the basis of sensitivity to known odorant molecules, which are described in the literature, were found in the expressed olfactory receptors gene set. [ABSTRACT FROM AUTHOR]

Background: There have been conflicting reports regarding the function of miR-20a in a variety of cancer types and we previously found it to be dysregulated in sporadic versus familial papillary thyroid cancer. In this study, we studied the expression of miR-20a in normal, benign and malignant thyroid samples, and its effect on thyroid cancer cells in vitro and in vivo. Methodology/Principal Findings: The expression of miR-20a in normal, benign and malignant thyroid tissue was determined by quantitative RT-PCR. Thyroid cancer cells were transfected with miR-20a and the effect on cellular proliferation, tumor spheroid formation, and invasion was evaluated. Target genes of miR-20 were determined by genome-wide mRNA expression analysis with miR-20a overexpression in thyroid cancer cells and target prediction database. Target genes were validated by quantitative PCR and immunoblotting, and luciferase assays. MiR-20a expression was significantly higher in anaplastic thyroid cancer than in differentiated thyroid cancer, and benign and normal thyroid tissues. MiR-20a significantly inhibited thyroid cancer cell proliferation in vitro (p<0.01) and in vivo (p<0.01), tumor spheroid formation (p<0.05) and invasion (p<0.05) in multiple thyroid cancer cell lines. We found that LIMK1 was a target of miR-20a in thyroid cancer cell lines and direct knockdown of LIMK1 recapitulated the effect of miR-20a in thyroid cancer cells. Conclusions/Significance: To our knowledge, this is the first study to demonstrate that miR-20a plays a role as a tumor suppressor in thyroid cancer cells and targets LIMK1. Our findings suggest the upregulated expression of miR-20a in anaplastic thyroid cancer counteracts thyroid cancer progression and may have therapeutic potential. [ABSTRACT FROM AUTHOR]

Statins exert atheroprotective effects through the induction of specific transcriptional factors in multiple organs. In endothelial cells, statin-dependent atheroprotective gene up-regulation is mediated by Kruppel-like factor (KLF) family transcription factors. To dissect the mechanism of gene regulation, we sought to determine molecular targets by performing microarray analyses of human umbilical vein endothelial cells (HUVECs) treated with pitavastatin, and KLF4 was determined to be the most highly induced gene. In addition, it was revealed that the atheroprotective genes induced with pitavastatin, such as nitric oxide synthase 3 (NOS3) and thrombomodulin (THBD), were suppressed by KLF4 knockdown. Myocyte enhancer factor-2 (MEF2) family activation is reported to be involved in pitavastatin-dependent KLF4 induction. We focused on MEF2C among the MEF2 family members and identified a novel functional MEF2C binding site 148 kb upstream of the KLF4 gene by chromatin immunoprecipitation along with deep sequencing (ChIP-seq) followed by luciferase assay. By applying whole genome and quantitative chromatin conformation analysis {chromatin interaction analysis with paired end tag sequencing (ChIA-PET), and real time chromosome conformation capture (3C) assay}, we observed that the MEF2C-bound enhancer and transcription start site (TSS) of KLF4 came into closer spatial proximity by pitavastatin treatment. 3D-Fluorescence in situ hybridization (FISH) imaging supported the conformational change in individual cells. Taken together, dynamic chromatin conformation change was shown to mediate pitavastatin-responsive gene induction in endothelial cells. [ABSTRACT FROM AUTHOR]

Purpose: The aim of the present study was to examine the biological differences between seminomas with occult and clinically apparent metastases at the time of diagnosis of the primary tumor to gain insight into the biology of these tumors and facilitate the identification of novel predictors of seminoma metastasis. Materials and Methods: Total RNA including small RNAs was isolated from testicular tumors of patients with pure seminoma presenting with lymphogenic metastasis (n = 5, clinical stage IIb/c) and occult metastasis (n = 5, clinical stage I). The regulation of biological processes was examined (1) throughout the mRNA transcriptome (whole genome microarrays, 8×60 K Array, Agilent with 4 samples/group) and (2) the miRNA transcriptome employing small RNA next generation sequencing (SOLID, Life Technologies with 5 samples/group). Protein coding genes (mRNAs) and small RNAs showing a significant (≥2-fold) difference between the groups were identified. Finally (3), we examined 95 candidate miRNAs in 36 apparent metastasized and another 5 occult metastasized seminoma using logistic regression analysis. Results: Among 19,596 genes, on average 12,894 mRNAs appeared expressed (65.8%, SD+/−2.4; range, 62.0–69.3%) and 16.99×106/13.94×106 small RNA reads were identified for apparent/occult metastasized seminoma. These reads on average convert into 9,901/9,675 small RNAs including 422/404 mature microRNAs. None of these mRNAs/small RNAs met our selection criteria for candidate genes. From 95 candidate miRNAs 44 appeared expressed, with 3 of them showing weak but significant (p = 0.05) differences among both groups. Conclusions: Occult and apparent metastasized seminomas are biologically almost indistinguishable and probably represent no separate tumor entities. These findings may simplify future research on seminoma metastasis. [ABSTRACT FROM AUTHOR]

Chinese medicine is a complex system guided by traditional Chinese medicine (TCM) theories, which has proven to be especially effective in treating chronic and complex diseases. However, the underlying modes of action (MOA) are not always systematically investigated. Herein, a systematic study was designed to elucidate the multi-compound, multi-target and multi-pathway MOA of a Chinese medicine, QiShenYiQi (QSYQ), on myocardial infarction. QSYQ is composed of Astragalus membranaceus (Huangqi), Salvia miltiorrhiza (Danshen), Panax notoginseng (Sanqi), and Dalbergia odorifera (Jiangxiang). Male Sprague Dawley rat model of myocardial infarction were administered QSYQ intragastrically for 7 days while the control group was not treated. The differentially expressed genes (DEGs) were identified from myocardial infarction rat model treated with QSYQ, followed by constructing a cardiovascular disease (CVD)-related multilevel compound-target-pathway network connecting main compounds to those DEGs supported by literature evidences and the pathways that are functionally enriched in ArrayTrack. 55 potential targets of QSYQ were identified, of which 14 were confirmed in CVD-related literatures with experimental supporting evidences. Furthermore, three sesquiterpene components of QSYQ, Trans-nerolidol, (3S,6S,7R)-3,7,11-trimethyl-3,6-epoxy-1,10-dodecadien-7-ol and (3S,6R,7R)-3,7,11-trimethyl-3,6-epoxy-1,10-dodecadien-7-ol from Dalbergia odorifera T. Chen, were validated experimentally in this study. Their anti-inflammatory effects and potential targets including extracellular signal-regulated kinase-1/2, peroxisome proliferator-activated receptor-gamma and heme oxygenase-1 were identified. Finally, through a three-level compound-target-pathway network with experimental analysis, our study depicts a complex MOA of QSYQ on myocardial infarction. [ABSTRACT FROM AUTHOR]

Introduction: Animal models for the study of sepsis are being increasingly scrutinized, despite their essential role for early translational research. In particular, recent studies have suggested that at the level of the leukocyte transcriptome, murine models of burns, trauma and endotoxemia markedly differ from their human equivalents, and are only weakly similar amongst themselves. We compared the plasma cytokine and leukocyte transcriptome responses between two different low-lethality murine models of polymicrobial intra-abdominal sepsis. Methods: Six to ten week male C57BL/6j mice underwent either the ‘gold standard’ cecal ligation and puncture (CLP) model of intra-abdominal sepsis or administration of a cecal slurry (CS), where cecal contents are injected intraperitoneally. Surviving mice were euthanized at two hours, one or three days after sepsis. Results: The murine leukocyte transcriptomic response to the CLP and CS models of sepsis was surprisingly dissimilar at two hours, one, and three days after sepsis. The Pearson correlation coefficient for the maximum change in expression for the entire leukocyte transcriptome that changed significantly over time (n = 19,071) was R = 0.54 (R2 = 0.297). The CS model resulted in greater magnitude of early inflammatory gene expression changes in response to sepsis with associated increased production of inflammatory chemokines and cytokines. Two hours after sepsis, CLP had more significant expression of genes associated with IL-10 signaling pathways, whereas CS had greater expression of genes related to CD28, apoptosis, IL-1 and T-cell receptor signaling. By three days, the changes in gene expression in both sepsis models were returning to baseline in surviving animals. Conclusion: These analyses reveal that the murine blood leukocyte response to sepsis is highly dependent on which model of intra-abdominal sepsis is employed, despite their similar lethality. It may be difficult to extrapolate findings from one murine model to another, let alone to human sepsis. [ABSTRACT FROM AUTHOR]

Background: Gene expression studies of peripheral blood mononuclear cells from patients with systemic lupus erythematosus (SLE) have demonstrated a type I interferon signature and increased expression of inflammatory cytokine genes. Studies of patients with Aicardi Goutières syndrome, commonly cited as a single gene model for SLE, have suggested that accumulation of non-coding RNAs may drive some of the pathologic gene expression, however, no RNA sequencing studies of SLE patients have been performed. This study was designed to define altered expression of coding and non-coding RNAs and to detect globally altered RNA processing in SLE. Methods: Purified monocytes from eight healthy age/gender matched controls and nine SLE patients (with low-moderate disease activity and lack of biologic drug use or immune suppressive treatment) were studied using RNA-seq. Quantitative RT-PCR was used to validate findings. Serum levels of endotoxin were measured by ELISA. Results: We found that SLE patients had diminished expression of most endogenous retroviruses and small nucleolar RNAs, but exhibited increased expression of pri-miRNAs. Splicing patterns and polyadenylation were significantly altered. In addition, SLE monocytes expressed novel transcripts, an effect that was replicated by LPS treatment of control monocytes. We further identified increased circulating endotoxin in SLE patients. Conclusions: Monocytes from SLE patients exhibit globally dysregulated gene expression. The transcriptome is not simply altered by the transcriptional activation of a set of genes, but is qualitatively different in SLE. The identification of novel loci, inducible by LPS, suggests that chronic microbial translocation could contribute to the immunologic dysregulation in SLE, a new potential disease mechanism. [ABSTRACT FROM AUTHOR]

Pulmonary surfactant is required for lung function at birth and throughout postnatal life. Defects in the surfactant system are associated with common pulmonary disorders including neonatal respiratory distress syndrome and acute respiratory distress syndrome in children and adults. Lipogenesis is essential for the synthesis of pulmonary surfactant by type II epithelial cells lining the alveoli. This study sought to identify the role of pulmonary epithelial SREBP, a transcriptional regulator of cellular lipid homeostasis, during a critical time period of perinatal lung maturation in the mouse. Genome wide mRNA expression profiling of lung tissue from transgenic mice with epithelial-specific deletions of Scap (ScapΔ/Δ, resulting in inactivation of SREBP signaling) or Insig1 and Insig2 (Insig1/2Δ/Δ, resulting in activation of SREBP signaling) was assessed. Differentially expressed genes responding to SREBP perturbations were identified and subjected to functional enrichment analysis, pathway mapping and literature mining to predict upstream regulators and transcriptional networks regulating surfactant lipid homeostasis. Through comprehensive data analysis and integration, time dependent effects of epithelial SCAP/INSIG/SREBP deletion and defined SCAP/INSIG/SREBP-associated genes, bioprocesses and downstream pathways were identified. SREBP signaling influences epithelial development, cell death and cell proliferation at E17.5, while primarily influencing surfactant physiology, lipid/sterol synthesis, and phospholipid transport after birth. SREBP signaling integrated with the Wnt/β-catenin and glucocorticoid receptor signaling pathways during perinatal lung maturation. SREBP regulates perinatal lung lipogenesis and maturation through multiple mechanisms by interactions with distinct sets of regulatory partners. [ABSTRACT FROM AUTHOR]

The edible mushroom Volvariella volvacea is an important crop in Southeast Asia and is predominantly harvested in the egg stage. One of the main factors that negatively affect its yield and value is the rapid transition from the egg to the elongation stage, which has a decreased commodity value and shelf life. To improve our understanding of the changes during stipe development and the transition from egg to elongation stage in particular, we analyzed gene transcription in stipe tissue of V. volvacea using 3′-tag based digital expression profiling. Stipe development turned out to be fairly complex with high numbers of expressed genes, and regulation of stage differences is mediated mainly by changes in expression levels of genes, rather than on/off modulation. Most explicit is the strong up-regulation of cell division from button to egg, and the very strong down-regulation hereof from egg to elongation, that continues in the maturation stage. Button and egg share cell division as means of growth, followed by a major developmental shift towards rapid stipe elongation based on cell extension as demonstrated by inactivation of cell division throughout elongation and maturation. Examination of regulatory genes up-regulated from egg to elongation identified three potential high upstream regulators for this switch. The new insights in stipe dynamics, together with a series of new target genes, will provide a sound base for further studies on the developmental mechanisms of mushroom stipes and the switch from egg to elongation in V. volvacea in particular. [ABSTRACT FROM AUTHOR]

Regulation of transcription is one of the mechanisms involved in animal development, directing changes in patterning and cell fate specification. Large temporal data series, based on microarrays across the life cycle of the fly Drosophila melanogaster, revealed the existence of groups of genes which expression increases or decreases temporally correlated during the life cycle. These groups of genes are enriched in different biological functions. Here, instead of searching for temporal coincidence in gene expression using the entire genome expression data, we searched for temporal coincidence in gene expression only within predefined catalogues of functionally related genes and investigated whether a catalogue's expression profile can be used to generate larger catalogues, enriched in genes necessary for the same function. We analyzed the expression profiles from genes already associated with early neurodevelopment and late neurodifferentiation, at embryonic stages 16 and 17 of Drosophila life cycle. We hypothesized that during this interval we would find global downregulation of genes important for early neuronal development together with global upregulation of genes necessary for the final differentiation of neurons. Our results were consistent with this hypothesis. We then investigated if the expression profile of gene catalogues representing particular processes of neural development matched the temporal sequence along which these processes occur. The profiles of genes involved in patterning, neurogenesis, axogenesis or synaptic transmission matched the prediction, with largest transcript values at the time when the corresponding biological process takes place in the embryo. Furthermore, we obtained catalogues enriched in genes involved in temporally matching functions by performing a genome-wide systematic search for genes with their highest expression levels at the corresponding embryonic intervals. These findings imply the use of gene expression data in combination with known biological information to predict the involvement of functionally uncharacterized genes in particular biological events. [ABSTRACT FROM AUTHOR]