Draft genome of the red harvester ant Pogonomyrmex barbatus.

We report the draft genome sequence of the red harvester ant, Pogonomyrmex barbatus. The genome was sequenced using 454 pyrosequencing, and the current assembly and annotation were completed in less than 1 y. Analyses of conserved gene groups (more than 1,200 manually annotated genes to date) suggest a high-quality assembly and annotation comparable to recently sequenced insect genomes using Sanger sequencing. The red harvester ant is a model for studying reproductive division of labor, phenotypic plasticity, and sociogenomics. Although the genome of P. barbatus is similar to other sequenced hymenopterans (Apis mellifera and Nasonia vitripennis) in GC content and compositional organization, and possesses a complete CpG methylation toolkit, its predicted genomic CpG content differs markedly from the other hymenopterans. Gene networks involved in generating key differences between the queen and worker castes (e.g., wings and ovaries) show signatures of increased methylation and suggest that ants and bees may have independently co-opted the same gene regulatory mechanisms for reproductive division of labor. Gene family expansions (e.g., 344 functional odorant receptors) and pseudogene accumulation in chemoreception and P450 genes compared with A. mellifera and N. vitripennis are consistent with major life-history changes during the adaptive radiation of Pogonomyrmex spp., perhaps in parallel with the development of the North American deserts.

Revisiting sodium and water reabsorption with functional genomics tools.

PURPOSE OF REVIEW: The kidney plays an essential role in maintaining sodium and water balance, thereby controlling the volume and osmolarity of the extracellular body fluids, the blood volume and the blood pressure. The final adjustment of sodium and water reabsorption in the kidney takes place in cells of the distal part of the nephron in which a set of apical and basolateral transporters participate in vectorial sodium and water transport from the tubular lumen to the interstitium and, finally, to the general circulation. According to a current model, the activity and/or cell-surface expression of these transporters is/are under the control of a gene network composed of the hormonally regulated, as well as constitutively expressed, genes. It is proposed that this gene network may include new candidate genes for salt- and water-losing syndromes and for salt-sensitive hypertension. A new generation of functional genomics techniques have recently been applied to the characterization of this gene network. The purpose of this review is to summarize these studies and to discuss the potential of the different techniques for characterization of the renal transcriptome. RECENT FINDINGS: Recently, DNA microarrays and serial analysis of gene expression have been applied to characterize the kidney transcriptome in different in-vivo and in-vitro models. In these studies, a set of new interesting genes potentially involved in the regulation of sodium and water reabsorption by the kidney have been identified and are currently under detailed investigation. SUMMARY: Characterization of the kidney transcriptome is greatly expanding our knowledge of the gene networks involved in multiple kidney functions, including the maintenance of sodium and water homeostasis.

Differential miRNA and IncRNA expression in embryonic stem cells lacking the Notch1 receptor

Embryonic stem (ES) cells-derived cardiomyocytes represent an attractive source of cells in cell replacement therapies for heart disease. However, controlled cardiogenic differentiation of ES cells requires a complete understanding of the complex molecular mechanisms regulating the differentiation process. We have previously shown that differentiation of ES cells into cardiomyocytes is favored by inactivation of the Notch 1 receptor pathway. In the present study, we therefore compared two ES cell lines, one with normal Notchl expression and one carrying deleted Notchl receptor alleles (Notchl-deleted ES cells) in order to identify genes responsible for the increased propensity of Notchl-deleted ES cells to produce cardiomyocytes. Using RNA-sequencing, we found approximately 300 coding and noncoding transcripts, which are differently expressed in undifferentiated Notchl-deleted ES cells. Since accumulating evidences indicate that long noncoding RNAs (IncRNAs) play important roles in ES cell pluripotency and differentiation, we focused our analysis on modulated IncRNAs. In particular, two IncRNAs, named here lnc 1230 and lnc 1335, are highly induced in the absence of Notchl receptor expression. These represent therefore prime candidates that could favor cardiogenic commitment in undifferentiated ES cells. Indeed, we demonstrate that forced expression of these two IncRNAs in wild-type ES cells result in a significant increase of the number of cardiac progenitor cells and cardiomyocytes in the differentiated progeny of these ES cells. Furthermore, we also identify several microRNAs that are differentially modulated in absence of Notchl expression. Among these are miR-142-5p and miR- 381-3p. Interestingly, both lncl230 and lncl335 are targets of these two microRNAs. Altogether, these data suggest that Notchl-dependent noncoding gene networks, implicating microRNAs and IncRNAs, control embryonic stem cell commitment into the mesodermal and cardiac lineages already at the undifferentiated state. - Les cardiomyocytes issus cellules souches embryonnaires sont une source très prometteuse pour les thérapies cellulaire de remplacement dans le cadre des maladies cardiaques. Cependant, l'utilisation de telles cellules requiert une compréhension poussée des mécanismes moléculaire régulant la différenciation. Nous avons par le passé démontré que la différenciation des cellules souches embryonnaires en cardiomyocytes est favorisée par l'inactivation de la voie d'activation intracellulaire dépendante du récepteur Notch 1. Nous avons donc comparé deux lignées de cellules souches embryonnaires, une présentant une voie d'activation Notchl normale et une chez laquelle les allèles codant pour le récepteur Notchl avaient été invalidés, de façon à identifier les gènes impliqués dans la capacité augmentée des cellules déficientes à produire des cardiomyocytes. En utilisant du séquençage d'ARN à haut débit, nous avons trouvé environ 300 gènes différemment exprimés dans les cellules déficientes pour Notchl. Par ailleurs, des évidences de plus en plus nombreuses suggèrent qu'une nouvelle classe de molécules appelée « long noncoding RNAs » joue un rôle prépondérant dans la maintenance de l'état non différencié et de la capacité de différenciation des cellules souches embryonnaires. Nous avons trouvé que plusieurs « long noncoding RNAs » étaient modulés en l'absence de Notchl, et en particulier deux molécules que nous avons appelées lncl230 et lncl335. Ces derniers représentent des candidats potentiels devant permettre de favoriser la production de cardiomyocytes. Nous avons en effet démontré que la surexpression de ces deux candidats dans des cellules souches embryonnaires résultait en une surproduction de cardiomyocytes. De plus, nous avons également identifié plusieurs microRNAs dont l'expression était modulée dans les cellules souches embryonnaires déficientes dans la voie Notchl. De façon intéressante, parmi ces microRNAs, le miR-142-5p et le miR-381-3p sont capables de cibler lncl230 and lncl335. Dans l'ensemble, ces résultats indiquent donc que des réseaux d'interaction dépendant de la voie d'activation Notch 1 et impliquant des ARNs non codant existent dans les cellules souches embryonnaires pour réguler leur différenciation en différent types cellulaires spécifiques.

Modifier Genes as Therapeutics: The Nuclear Hormone Receptor Rev Erb Alpha (Nr1d1) Rescues Nr2e3 Associated Retinal Disease

Nuclear hormone receptors play a major role in many important biological processes. Most nuclear hormone receptors are ubiquitously expressed and regulate processes such as metabolism, circadian function, and development. They function in these processes to maintain homeostasis through modulation of transcriptional gene networks. In this study we evaluate the effectiveness of a nuclear hormone receptor gene to modulate retinal degeneration and restore the integrity of the retina. Currently, there are no effective treatment options for retinal degenerative diseases leading to progressive and irreversible blindness. In this study we demonstrate that the nuclear hormone receptor gene Nr1d1 (Rev-Erba) rescues Nr2e3- associated retinal degeneration in the rd7 mouse, which lacks a functional Nr2e3 gene. Mutations in human NR2E3 are associated with several retinal degenerations including enhanced S cone syndrome and retinitis pigmentosa. The rd7 mouse, lacking Nr2e3, exhibits an increase in S cones and slow, progressive retinal degeneration. A traditional genetic mapping approach previously identified candidate modifier loci. Here, we demonstrate that in vivo delivery of the candidate modifier gene, Nr1d1 rescues Nr2e3 associated retinal degeneration. We observed clinical, histological, functional, and molecular restoration of the rd7 retina. Furthermore, we demonstrate that the mechanism of rescue at the molecular and functional level is through the re-regulation of key genes within the Nr2e3-directed transcriptional network. Together, these findings reveal the potency of nuclear receptors as modulators of disease and specifically of NR1D1 as a novel therapeutic for retinal degenerations.

A cell spot microarray method for high-throughput biology

High-throughput screening of cellular effects of RNA interference (RNAi) libraries is now being increasingly applied to explore the role of genes in specific cell biological processes and disease states. However, the technology is still limited to specialty laboratories, due to the requirements for robotic infrastructure, access to expensive reagent libraries, expertise in high-throughput screening assay development, standardization, data analysis and applications. In the future, alternative screening platforms will be required to expand functional large-scale experiments to include more RNAi constructs, allow combinatorial loss-of-function analyses (e.g. genegene or gene-drug interaction), gain-of-function screens, multi-parametric phenotypic readouts or comparative analysis of many different cell types. Such comprehensive perturbation of gene networks in cells will require a major increase in the flexibility of the screening platforms, throughput and reduction of costs. As an alternative for the conventional multi-well based high-throughput screening -platforms, here the development of a novel cell spot microarray method for production of high density siRNA reverse transfection arrays is described. The cell spot microarray platform is distinguished from the majority of other transfection cell microarray techniques by the spatially confined array layout that allow highly parallel screening of large-scale RNAi reagent libraries with assays otherwise difficult or not applicable to high-throughput screening. This study depicts the development of the cell spot microarray method along with biological application examples of high-content immunofluorescence and phenotype based cancer cell biological analyses focusing on the regulation of prostate cancer cell growth, maintenance of genomic integrity in breast cancer cells, and functional analysis of integrin protein-protein interactions in situ.

Coping with genetic diversity: the contribution of pathogen and human genomics to modern vaccinology

Vaccine development faces major difficulties partly because of genetic variation in both infectious organisms and humans. This causes antigenic variation in infectious agents and a high interindividual variability in the human response to the vaccine. The exponential growth of genome sequence information has induced a shift from conventional culture-based to genome-based vaccinology, and allows the tackling of challenges in vaccine development due to pathogen genetic variability. Additionally, recent advances in immunogenetics and genomics should help in the understanding of the influence of genetic factors on the interindividual and interpopulation variations in immune responses to vaccines, and could be useful for developing new vaccine strategies. Accumulating results provide evidence for the existence of a number of genes involved in protective immune responses that are induced either by natural infections or vaccines. Variation in immune responses could be viewed as the result of a perturbation of gene networks; this should help in understanding how a particular polymorphism or a combination thereof could affect protective immune responses. Here we will present: i) the first genome-based vaccines that served as proof of concept, and that provided new critical insights into vaccine development strategies; ii) an overview of genetic predisposition in infectious diseases and genetic control in responses to vaccines; iii) population genetic differences that are a rationale behind group-targeted vaccines; iv) an outlook for genetic control in infectious diseases, with special emphasis on the concept of molecular networks that will provide a structure to the huge amount of genomic data.

Human Stem Cell Cultures from Cleft Lip/Palate Patients Show Enrichment of Transcripts Involved in Extracellular Matrix Modeling By Comparison to Controls

Nonsyndromic cleft lip and palate (NSCL/P) is a complex disease resulting from failure of fusion of facial primordia, a complex developmental process that includes the epithelial-mesenchymal transition (EMT). Detection of differential gene transcription between NSCL/P patients and control individuals offers an interesting alternative for investigating pathways involved in disease manifestation. Here we compared the transcriptome of 6 dental pulp stem cell (DPSC) cultures from NSCL/P patients and 6 controls. Eighty-seven differentially expressed genes (DEGs) were identified. The most significant putative gene network comprised 13 out of 87 DEGs of which 8 encode extracellular proteins: ACAN, COL4A1, COL4A2, GDF15, IGF2, MMP1, MMP3 and PDGFa. Through clustering analyses we also observed that MMP3, ACAN, COL4A1 and COL4A2 exhibit co-regulated expression. Interestingly, it is known that MMP3 cleavages a wide range of extracellular proteins, including the collagens IV, V, IX, X, proteoglycans, fibronectin and laminin. It is also capable of activating other MMPs. Moreover, MMP3 had previously been associated with NSCL/P. The same general pattern was observed in a further sample, confirming involvement of synchronized gene expression patterns which differed between NSCL/P patients and controls. These results show the robustness of our methodology for the detection of differentially expressed genes using the RankProd method. In conclusion, DPSCs from NSCL/P patients exhibit gene expression signatures involving genes associated with mechanisms of extracellular matrix modeling and palate EMT processes which differ from those observed in controls. This comparative approach should lead to a more rapid identification of gene networks predisposing to this complex malformation syndrome than conventional gene mapping technologies.

Expression of genes related to apoptosis, cell cycle and signaling pathways are independent of TP53 status in urinary bladder cancer cells

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

Genomic signatures of evolutionary transitions from solitary to group living

The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

Apoptosis rate and transcriptional response of pancreatic islets exposed to the PPAR gamma agonist Pioglitazone

To explore the molecular pathways underlying thiazolidinediones effects on pancreatic islets in conditions mimicking normo- and hyperglycemia, apoptosis rate and transcriptional response to Pioglitazone at both physiological and supraphysiological glucose concentrations were evaluated. Adult rat islets were cultured at physiological (5.6 mM) and supraphysiological (23 mM) glucose concentrations in presence of 10 μM Pioglitazone or vehicle. RNA expression profiling was evaluated with the PancChip 13k cDNA microarray after 24-h, and expression results for some selected genes were validated by qRT-PCR. The effects of Pioglitazone were investigated regarding apoptosis rate after 24-, 48- and 72-h. At 5.6 mM glucose, 101 genes were modulated by Pioglitazone, while 1,235 genes were affected at 23 mM glucose. Gene networks related to lipid metabolism were identified as altered by Pioglitazone at both glucose concentrations. At 23 mM glucose, cell cycle and cell death pathways were significantly regulated as well. At 5.6 mM glucose, Pioglitazone elicited a transient reduction in islets apoptosis rate while at 23 mM, Bcl2 expression was reduced and apoptosis rate was increased by Pioglitazone. Our data demonstrate that the effect of Pioglitazone on gene expression profile and apoptosis rate depends on the glucose concentration. The modulation of genes related to cell death and the increased apoptosis rate observed at supraphysiological glucose concentration raise concerns about Pioglitazone’s direct effects in conditions of hyperglycemia and reinforce the necessity of additional studies designed to evaluate TZDs effects on the preservation of β-cell function in situations where glucotoxicity might be more relevant than lipotoxicity.

Development of synthetic molecular circuits for the control of cell systems

Synthetic biology is a young field of applicative research aiming to design and build up artificial biological devices, useful for human applications. How synthetic biology emerged in past years and how the development of the Registry of Standard Biological Parts aimed to introduce one practical starting solution to apply the basics of engineering to molecular biology is presented in chapter 1 in the thesis The same chapter recalls how biological parts can make up a genetic program, the molecular cloning tecnique useful for this purpose, and an overview of the mathematical modeling adopted to describe gene circuit behavior. Although the design of gene circuits has become feasible the increasing complexity of gene networks asks for a rational approach to design gene circuits. A bottom-up approach was proposed, suggesting that the behavior of a complicated system can be predicted from the features of its parts. The option to use modular parts in large-scale networks will be facilitated by a detailed and shared characterization of their functional properties. Such a prediction, requires well-characterized mathematical models of the parts and of how they behave when assembled together. In chapter 2, the feasibility of the bottom-up approach in the design of a synthetic program in Escherichia coli bacterial cells is described. The rational design of gene networks is however far from being established. The synthetic biology approach can used the mathematical formalism to identify biological information not assessable with experimental measurements. In this context, chapter 3 describes the design of a synthetic sensor for identifying molecules of interest inside eukaryotic cells. The Registry of Standard parts collects standard and modular biological parts. To spread the use of BioBricks the iGEM competition was started. The ICM Laboratory, where Francesca Ceroni completed her Ph.D, partecipated with teams of students and Chapter 4 summarizes the projects developed.

Crosstalk between B lymphocytes, microbiota and the intestinal epithelium governs immunity versus metabolism in the gut

Using a systems biology approach, we discovered and dissected a three-way interaction between the immune system, the intestinal epithelium and the microbiota. We found that, in the absence of B cells, or of IgA, and in the presence of the microbiota, the intestinal epithelium launches its own protective mechanisms, upregulating interferon-inducible immune response pathways and simultaneously repressing Gata4-related metabolic functions. This shift in intestinal function leads to lipid malabsorption and decreased deposition of body fat. Network analysis revealed the presence of two interconnected epithelial-cell gene networks, one governing lipid metabolism and another regulating immunity, that were inversely expressed. Gene expression patterns in gut biopsies from individuals with common variable immunodeficiency or with HIV infection and intestinal malabsorption were very similar to those of the B cell-deficient mice, providing a possible explanation for a longstanding enigmatic association between immunodeficiency and defective lipid absorption in humans.

IDENTIFICATION OF GENES CONTROLLING BIOLOGICAL PROCESSES AND PATHWAYS THROUGH STATISTICAL ANALYSIS AND NETWORK RECONSTRUCTION

The developmental processes and functions of an organism are controlled by the genes and the proteins that are derived from these genes. The identification of key genes and the reconstruction of gene networks can provide a model to help us understand the regulatory mechanisms for the initiation and progression of biological processes or functional abnormalities (e.g. diseases) in living organisms. In this dissertation, I have developed statistical methods to identify the genes and transcription factors (TFs) involved in biological processes, constructed their regulatory networks, and also evaluated some existing association methods to find robust methods for coexpression analyses. Two kinds of data sets were used for this work: genotype data and gene expression microarray data. On the basis of these data sets, this dissertation has two major parts, together forming six chapters. The first part deals with developing association methods for rare variants using genotype data (chapter 4 and 5). The second part deals with developing and/or evaluating statistical methods to identify genes and TFs involved in biological processes, and construction of their regulatory networks using gene expression data (chapter 2, 3, and 6). For the first part, I have developed two methods to find the groupwise association of rare variants with given diseases or traits. The first method is based on kernel machine learning and can be applied to both quantitative as well as qualitative traits. Simulation results showed that the proposed method has improved power over the existing weighted sum method (WS) in most settings. The second method uses multiple phenotypes to select a few top significant genes. It then finds the association of each gene with each phenotype while controlling the population stratification by adjusting the data for ancestry using principal components. This method was applied to GAW 17 data and was able to find several disease risk genes. For the second part, I have worked on three problems. First problem involved evaluation of eight gene association methods. A very comprehensive comparison of these methods with further analysis clearly demonstrates the distinct and common performance of these eight gene association methods. For the second problem, an algorithm named the bottom-up graphical Gaussian model was developed to identify the TFs that regulate pathway genes and reconstruct their hierarchical regulatory networks. This algorithm has produced very significant results and it is the first report to produce such hierarchical networks for these pathways. The third problem dealt with developing another algorithm called the top-down graphical Gaussian model that identifies the network governed by a specific TF. The network produced by the algorithm is proven to be of very high accuracy.

Retinoic acid receptor beta and angiopoietin-like protein 1 are involved in the regulation of human androgen biosynthesis

Androgens are essential for sexual development and reproduction. However, androgen regulation in health and disease is poorly understood. We showed that human adrenocortical H295R cells grown under starvation conditions acquire a hyperandrogenic steroid profile with changes in steroid metabolizing enzymes HSD3B2 and CYP17A1 essential for androgen production. Here we studied the regulatory mechanisms underlying androgen production in starved H295R cells. Microarray expression profiling of normal versus starved H295R cells revealed fourteen differentially expressed genes; HSD3B2, HSD3B1, CYP21A2, RARB, ASS1, CFI, ASCL1 and ENC1 play a role in steroid and energy metabolism and ANGPTL1, PLK2, DUSP6, DUSP10 and FREM2 are involved in signal transduction. We discovered two new gene networks around RARB and ANGPTL1, and show how they regulate androgen biosynthesis. Transcription factor RARB stimulated the promoters of genes involved in androgen production (StAR, CYP17A1 and HSD3B2) and enhanced androstenedione production. For HSD3B2 regulation RARB worked in cooperation with Nur77. Secretory protein ANGPTL1 modulated CYP17A1 and DUSP6 expression by inducing ERK1/2 phosphorylation. By contrast, our studies revealed no evidence for hormones or cell cycle involvement in regulating androgen biosynthesis. In summary, these studies establish a firm role for RARB and ANGPTL1 in the regulation of androgen production in H295R cells.

Reactive oxygen species-dependent mechanisms of N-(4-hydroxyphenyl)retinamide (4HPR)-induced apoptosis in human carcinoma cells

4HPR is a synthetic retinoid that has shown chemopreventive and therapeutic efficacy against premalignant and malignant lesions including oral leukoplakia, ovarian and breast cancer, and neuroblastoma. 4HPR induces apoptosis in various cancer cells and production of reactive oxygen species (ROS) has been suggested as a possible cause underlying these effects. However, the mechanisms governing these effects by 4HPR are not fully elucidated. In this study, we explored the mechanisms of 4HPR-induced ROS increase and apoptosis in human cancer cells. ^ First, we identified genes modulated by 4HPR using oligonucleotide gene expression arrays and found that they fall into specific functional canonical pathways and gene networks using Ingenuity Pathways Analysis®. Further analysis has shown that 4HPR induced up-regulation of Endoplasmic Reticulum (ER)-related genes such as Heat shock proteins 70 and 90 and the transcriptional factor, GADD153. These findings were validated using quantitative real-time PCR. ^ Second, we found that 4HPR induced extensive ER stress evidenced by dilation of the ER and endoribonuclease-mediated splicing and activation of the transcriptional factor, XBP-1. In addition, 4HPR induced the up-regulation of various ER stress-related genes and their protein products, as well as cleavage and activation of the ER specific Caspase-4. Concomitantly with XBP-1 splicing, all of these effects were dependent on ROS generation by 4HPR. Furthermore, chemical inhibition and RNA interference studies revealed a novel pro-apoptotic role for HSP70/A1A in 4HPR-mediated apoptosis. ^ Third, we observed rapid activation of the small GTPase Rac by 4HPR which was upstream of ROS generation. Inhibition of Rac activity or silencing of its expression by RNA interference inhibited ROS generation and apoptosis induction by 4HPR. siRNA targeting PAK1 and expression of a dominant negative Rac, decreased 4HPR-mediated ROS generation, while expression of a constitutive active Rac increased basal and 4HPR-induced ROS generation and PARP cleavage. Furthermore, metastatic cancer cells exhibited higher Rac activation, ROS generation, and cell growth inhibition due to 4HPR exposure compared to their primary cancer cell counterparts. ^ These findings provide novel insights into 4HPR-mediated ROS generation and apoptosis induction and support the use of ROS inducing agents such as 4HPR against metastatic cancer cells. ^