Gene expression changes in the coccolithophore Emiliania huxleyi after 500 generations of selection to ocean acidification

Coccolithophores are unicellular marine algae that produce biogenic calcite scales and substantially contribute to marine primary production and carbon export to the deep ocean. Ongoing ocean acidification particularly impairs calcifying organisms, mostly resulting in decreased growth and calcification. Recent studies revealed that the immediate physiological response in the coccolithophore Emiliania huxleyi to ocean acidification may be partially compensated by evolutionary adaptation, yet the underlying molecular mechanisms are currently unknown. Here, we report on the expression levels of 10 candidate genes putatively relevant to pH regulation, carbon transport, calcification and photosynthesis in E. huxleyi populations short-term exposed to ocean acidification conditions after acclimation (physiological response) and after 500 generations of high CO2 adaptation (adaptive response). The physiological response revealed downregulation of candidate genes, well reflecting the concomitant decrease of growth and calcification. In the adaptive response, putative pH regulation and carbon transport genes were up-regulated, matching partial restoration of growth and calcification in high CO2-adapted populations. Adaptation to ocean acidification in E. huxleyi likely involved improved cellular pH regulation, presumably indirectly affecting calcification. Adaptive evolution may thus have the potential to partially restore cellular pH regulatory capacity and thereby mitigate adverse effects of ocean acidification.

Image Processing Challenges in the Creation of Spatiotemporal Gene Expression Atlases of Developing Embryos

To properly understand and model animal embryogenesis it is crucial to obtain detailed measurements, both in time and space, about their gene expression domains and cell dynamics. Such challenge has been confronted in recent years by a surge of atlases which integrate a statistically relevant number of different individuals to get robust, complete information about their spatiotemporal locations of gene patterns. This paper will discuss the fundamental image analysis strategies required to build such models and the most common problems found along the way. We also discuss the main challenges and future goals in the field.

Image Processing Challenges in the Creation of Spatiotemporal Gene Expression Atlases of Developing Embryos

To properly understand and model animal embryogenesis it is crucial to obtain detailed measurements, both in time and space, about their gene expression domains and cell dynamics. Such challenge has been confronted in recent years by a surge of atlases which integrate a statistically relevant number of different individuals to get robust, complete information about their spatiotemporal locations of gene patterns. This paper will discuss the fundamental image analysis strategies required to build such models and the most common problems found along the way. We also discuss the main challenges and future goals in the field.

An Automatic Quantification and Registration Strategy to Create a Gene Expression Atlas of Zebrafish Embryogenesis

In order to properly understand and model the gene regulatory networks in animals development, it is crucial to obtain detailed measurements, both in time and space, about their gene expression domains. In this paper, we propose a complete computational framework to fulfill this task and create a 3D Atlas of the early zebrafish embryogenesis annotated with both the cellular localizations and the level of expression of different genes at different developmental stages. The strategy to construct such an Atlas is described here with the expression pattern of 5 different genes at 6 hours of development post fertilization.

Gene Expression Profiling in Limb-Girdle Muscular Dystrophy 2A

Limb-girdle muscular dystrophy type 2A (LGMD2A) is a recessive genetic disorder caused by mutations in calpain 3 (CAPN3). Calpain 3 plays different roles in muscular cells, but little is known about its functions or in vivo substrates. The aim of this study was to identify the genes showing an altered expression in LGMD2A patients and the possible pathways they are implicated in. Ten muscle samples from LGMD2A patients with in which molecular diagnosis was ascertained were investigated using array technology to analyze gene expression profiling as compared to ten normal muscle samples. Upregulated genes were mostly those related to extracellular matrix (different collagens), cell adhesion (fibronectin), muscle development (myosins and melusin) and signal transduction. It is therefore suggested that different proteins located or participating in the costameric region are implicated in processes regulated by calpain 3 during skeletal muscle development. Genes participating in the ubiquitin proteasome degradation pathway were found to be deregulated in LGMD2A patients, suggesting that regulation of this pathway may be under the control of calpain 3 activity. As frizzled-related protein (FRZB) is upregulated in LGMD2A muscle samples, it could be hypothesized that β-catenin regulation is also altered at the Wnt signaling pathway, leading to an incorrect myogenesis. Conversely, expression of most transcription factor genes was downregulated (MYC, FOS and EGR1). Finally, the upregulation of IL-32 and immunoglobulin genes may induce the eosinophil chemoattraction explaining the inflammatory findings observed in presymptomatic stages. The obtained results try to shed some light on identification of novel therapeutic targets for limb-girdle muscular dystrophies

Computational methods to create and analyze a digital gene expression atlas of embryo development from microscopy images

Abstract The creation of atlases, or digital models where information from different subjects can be combined, is a field of increasing interest in biomedical imaging. When a single image does not contain enough information to appropriately describe the organism under study, it is then necessary to acquire images of several individuals, each of them containing complementary data with respect to the rest of the components in the cohort. This approach allows creating digital prototypes, ranging from anatomical atlases of human patients and organs, obtained for instance from Magnetic Resonance Imaging, to gene expression cartographies of embryo development, typically achieved from Light Microscopy. Within such context, in this PhD Thesis we propose, develop and validate new dedicated image processing methodologies that, based on image registration techniques, bring information from multiple individuals into alignment within a single digital atlas model. We also elaborate a dedicated software visualization platform to explore the resulting wealth of multi-dimensional data and novel analysis algo-rithms to automatically mine the generated resource in search of bio¬logical insights. In particular, this work focuses on gene expression data from developing zebrafish embryos imaged at the cellular resolution level with Two-Photon Laser Scanning Microscopy. Disposing of quantitative measurements relating multiple gene expressions to cell position and their evolution in time is a fundamental prerequisite to understand embryogenesis multi-scale processes. However, the number of gene expressions that can be simultaneously stained in one acquisition is limited due to optical and labeling constraints. These limitations motivate the implementation of atlasing strategies that can recreate a virtual gene expression multiplex. The developed computational tools have been tested in two different scenarios. The first one is the early zebrafish embryogenesis where the resulting atlas constitutes a link between the phenotype and the genotype at the cellular level. The second one is the late zebrafish brain where the resulting atlas allows studies relating gene expression to brain regionalization and neurogenesis. The proposed computational frameworks have been adapted to the requirements of both scenarios, such as the integration of partial views of the embryo into a whole embryo model with cellular resolution or the registration of anatom¬ical traits with deformable transformation models non-dependent on any specific labeling. The software implementation of the atlas generation tool (Match-IT) and the visualization platform (Atlas-IT) together with the gene expression atlas resources developed in this Thesis are to be made freely available to the scientific community. Lastly, a novel proof-of-concept experiment integrates for the first time 3D gene expression atlas resources with cell lineages extracted from live embryos, opening up the door to correlate genetic and cellular spatio-temporal dynamics. La creación de atlas, o modelos digitales, donde la información de distintos sujetos puede ser combinada, es un campo de creciente interés en imagen biomédica. Cuando una sola imagen no contiene suficientes datos como para describir apropiadamente el organismo objeto de estudio, se hace necesario adquirir imágenes de varios individuos, cada una de las cuales contiene información complementaria respecto al resto de componentes del grupo. De este modo, es posible crear prototipos digitales, que pueden ir desde atlas anatómicos de órganos y pacientes humanos, adquiridos por ejemplo mediante Resonancia Magnética, hasta cartografías de la expresión genética del desarrollo de embrionario, típicamente adquiridas mediante Microscopía Optica. Dentro de este contexto, en esta Tesis Doctoral se introducen, desarrollan y validan nuevos métodos de procesado de imagen que, basándose en técnicas de registro de imagen, son capaces de alinear imágenes y datos provenientes de múltiples individuos en un solo atlas digital. Además, se ha elaborado una plataforma de visualization específicamente diseñada para explorar la gran cantidad de datos, caracterizados por su multi-dimensionalidad, que resulta de estos métodos. Asimismo, se han propuesto novedosos algoritmos de análisis y minería de datos que permiten inspeccionar automáticamente los atlas generados en busca de conclusiones biológicas significativas. En particular, este trabajo se centra en datos de expresión genética del desarrollo embrionario del pez cebra, adquiridos mediante Microscopía dos fotones con resolución celular. Disponer de medidas cuantitativas que relacionen estas expresiones genéticas con las posiciones celulares y su evolución en el tiempo es un prerrequisito fundamental para comprender los procesos multi-escala característicos de la morfogénesis. Sin embargo, el número de expresiones genéticos que pueden ser simultáneamente etiquetados en una sola adquisición es reducido debido a limitaciones tanto ópticas como del etiquetado. Estas limitaciones requieren la implementación de estrategias de creación de atlas que puedan recrear un multiplexado virtual de expresiones genéticas. Las herramientas computacionales desarrolladas han sido validadas en dos escenarios distintos. El primer escenario es el desarrollo embrionario temprano del pez cebra, donde el atlas resultante permite constituir un vínculo, a nivel celular, entre el fenotipo y el genotipo de este organismo modelo. El segundo escenario corresponde a estadios tardíos del desarrollo del cerebro del pez cebra, donde el atlas resultante permite relacionar expresiones genéticas con la regionalización del cerebro y la formación de neuronas. La plataforma computacional desarrollada ha sido adaptada a los requisitos y retos planteados en ambos escenarios, como la integración, a resolución celular, de vistas parciales dentro de un modelo consistente en un embrión completo, o el alineamiento entre estructuras de referencia anatómica equivalentes, logrado mediante el uso de modelos de transformación deformables que no requieren ningún marcador específico. Está previsto poner a disposición de la comunidad científica tanto la herramienta de generación de atlas (Match-IT), como su plataforma de visualización (Atlas-IT), así como las bases de datos de expresión genética creadas a partir de estas herramientas. Por último, dentro de la presente Tesis Doctoral, se ha incluido una prueba conceptual innovadora que permite integrar los mencionados atlas de expresión genética tridimensionales dentro del linaje celular extraído de una adquisición in vivo de un embrión. Esta prueba conceptual abre la puerta a la posibilidad de correlar, por primera vez, las dinámicas espacio-temporales de genes y células.

Azotobacter vinelandii nitrogenase: “Kinetics of nif gene expression and insights into the roles of FdxN and NifQ in FeMo-co biosynthesis”

The Molybdenum-nitrogenase is responsible for most biological nitrogen fixation activity (BNF) in the biosphere. Due to its great agronomical importance, it has been the subject of profound genetic and biochemical studies. The Mo nitrogenase carries at its active site a unique iron-molybdenum cofactor (FeMoco) that consists of an inorganic 7 Fe, 1 Mo, 1 C, 9 S core coordinated to the organic acid homocitrate. Biosynthesis of FeMo-co occurs outside nitrogenase through a complex and highly regulated pathway involving proteins acting as molecular scaffolds, metallocluster carriers or enzymes that provide substrates in appropriate chemical forms. Specific expression regulatory factors tightly control the accumulation levels of all these other components. Insertion of FeMo-co into a P-cluster containing apo-NifDK polypeptide results in nitrogenase reconstitution. Investigation of FeMo-co biosynthesis has uncovered new radical chemistry reactions and new roles for Fe-S clusters in biology.

Growth rate and TRI5 gene expression profiles of Fusarium equiseti strains isolated from Spanish cereals cultivated on wheat and barley media at different environmental conditions

Fusarium equiseti is a toxigenic species that often contaminates ce real crops from diverse climatic regions such as Northern and Southern Europe. Previous results suggested the existence of two distinct populations within this species with differences in toxin pro file which largely corresponded to North and South Europe (Spain). In this work, growth rate profiles of 4 F. equiseti strains isolated from different cereals and distinct Spanish regions were determined on wheat and barley based media at a range of temperatures (15, 20, 25, 30, 35 and 40 °C) and water potentialregimens(−0.7,−2.8,−7.0,and −9.8MPa,correspondingto 0.99,0.98,0.95 and 0.93aw values).Growth was observed at all temperatures except at 40 °C, and at all the solute potential values except at−9.8 MPa when combined with 15 °C. Optimal growth was observed at 20– 30 °C and −0.7/−2.8 MPa. The effect of these factors on trichothecene biosynthesis was examined on a F. equiseti strain using a newly developed real time RT-PCR protocol to quantify TRI5 gene expression at 15, 25 and 35 °C and −0.7, −2.8, − 7.0 and −9.8 MPa on wheat and barley based media. Induction of TRI5 expression was detected between 25 and 35 °C and −0.7 and − 2.8 MPa, with maximum values at 35 °C and −2.8 MPa being higher in barley than in wheat medium. These results appeared to be consistent with a population well adapted to the present climatic conditions and predicted scenarios for Southern Europe and suggested some differences depending on the cereal considered. These are also discussed in relation to other Fusarium species co-occurring in cereals grown in this region and to their significance for prediction and control strategies of toxigenic risk in future scenarios of climate change for this region.

A digital framework to build, visualize and analyze a gene expression atlas with cellular resolution in zebrafish early embryogenesis

A gene expression atlas is an essential resource to quantify and understand the multiscale processes of embryogenesis in time and space. The automated reconstruction of a prototypic 4D atlas for vertebrate early embryos, using multicolor fluorescence in situ hybridization with nuclear counterstain, requires dedicated computational strategies. To this goal, we designed an original methodological framework implemented in a software tool called Match-IT. With only minimal human supervision, our system is able to gather gene expression patterns observed in different analyzed embryos with phenotypic variability and map them onto a series of common 3D templates over time, creating a 4D atlas. This framework was used to construct an atlas composed of 6 gene expression templates from a cohort of zebrafish early embryos spanning 6 developmental stages from 4 to 6.3 hpf (hours post fertilization). They included 53 specimens, 181,415 detected cell nuclei and the segmentation of 98 gene expression patterns observed in 3D for 9 different genes. In addition, an interactive visualization software, Atlas-IT, was developed to inspect, supervise and analyze the atlas. Match-IT and Atlas-IT, including user manuals, representative datasets and video tutorials, are publicly and freely available online. We also propose computational methods and tools for the quantitative assessment of the gene expression templates at the cellular scale, with the identification, visualization and analysis of coexpression patterns, synexpression groups and their dynamics through developmental stages.

Regulation of T cell receptor (TCR) β gene expression by CD3 complex signaling in immature thymocytes: Implications for TCRβ allelic exclusion

During αβ thymocyte development, clonotype-independent CD3 complexes are expressed at the cell surface before the pre-T cell receptor (TCR). Signaling through clonotype-independent CD3 complexes is required for expression of rearranged TCRβ genes. On expression of a TCRβ polypeptide chain, the pre-TCR is assembled, and TCRβ locus allelic exclusion is established. We investigated the putative contribution of clonotype-independent CD3 complex signaling to TCRβ locus allelic exclusion in mice single-deficient or double-deficient for CD3ζ/η and/or p56lck. These mice display defects in the expression of endogenous TCRβ genes in immature thymocytes, proportional to the severity of CD3 complex malfunction. Exclusion of endogenous TCRβ VDJ (variable, diversity, joining) rearrangements by a functional TCRβ transgene was severely compromised in the single-deficient and double-deficient mutant mice. In contrast to wild-type mice, most of the CD25+ double-negative (DN) thymocytes of the mutant mice failed to express the TCRβ transgene, suggesting defective expression of the TCRβ transgene similar to endogenous TCRβ genes. In the mutant mice, a proportion of CD25+ DN thymocytes that failed to express the transgene expressed endogenous TCRβ polypeptide chains. Many double-positive cells of the mutant mice coexpressed endogenous and transgenic TCRβ chains or more than one endogenous TCRβ chain. The data suggest that signaling through clonotype-independent CD3 complexes may contribute to allelic exclusion of the TCRβ locus by inducing the expression of rearranged TCRβ genes in CD25+ DN thymocytes.

Induction of ARF tumor suppressor gene expression and cell cycle arrest by transcription factor DMP1

Expression of the DMP1 transcription factor, a cyclin D-binding Myb-like protein, induces growth arrest in mouse embryo fibroblast strains but is devoid of antiproliferative activity in primary diploid fibroblasts that lack the ARF tumor suppressor gene. DMP1 binds to a single canonical recognition site in the ARF promoter to activate gene expression, and in turn, p19ARF synthesis causes p53-dependent cell cycle arrest. Unlike genes such as Myc, adenovirus E1A, and E2F-1, which, when overexpressed, activate the ARF-p53 pathway and trigger apoptosis, DMP1, like ARF itself, does not induce programmed cell death. Therefore, apart from its recently recognized role in protecting cells from potentially oncogenic signals, ARF can be induced in response to antiproliferative stimuli that do not obligatorily lead to apoptosis.

Loss of the circadian clock-associated protein 1 in Arabidopsis results in altered clock-regulated gene expression

Little is known about plant circadian oscillators, in spite of how important they are to sessile plants, which require accurate timekeepers that enable the plants to respond to their environment. Previously, we identified a circadian clock-associated (CCA1) gene that encodes an Myb-related protein that is associated with phytochrome control and circadian regulation in plants. To understand the role CCA1 plays in phytochrome and circadian regulation, we have isolated an Arabidopsis line with a T DNA insertion that results in the loss of CCA1 RNA, of CCA1 protein, and of an Lhcb-promoter binding activity. This mutation affects the circadian expression of all four clock-controlled genes that we examined. The results show that, despite their similarity, CCA1 and LHY are only partially redundant. The lack of CCA1 also affects the phytochrome regulation of gene expression, suggesting that CCA1 has an additional role in a signal transduction pathway from light, possibly acting at the point of integration between phytochrome and the clock. Our results indicate that CCA1 is an important clock-associated protein involved in circadian regulation of gene expression.

Cooperation between the activin and Wnt pathways in the spatial control of organizer gene expression

The normal expression pattern of the Wnt responsive homeobox gene Siamois is restricted to the dorso-vegetal region of the Xenopus embryo. Because the Wnt signaling pathway (via β-catenin) is active on the entire dorsal side of the early embryo, we have asked why Siamois expression is not seen in the dorsal ectoderm. Only Wnt signaling, via activation of β-catenin, can induce directly Siamois, and signaling via the SMAD1 (BMP2/4) or SMAD2 (activin/Vg-1) pathways cannot. We now directly show that the SMAD2 pathway can cooperate with the Wnt pathway to induce expression of Siamois much more strongly than the Wnt pathway alone, in normal embryos. We demonstrate the significance of this cooperation in normal embryos by blocking the SMAD2 signaling pathway with a dominant negative activin receptor. The activin dominant negative receptor blocks this cooperative effect and reduces the expression of Siamois by threefold in early embryos. Furthermore, we find that this cooperative relationship between the SMAD2 and Wnt pathways is reciprocal. Thus, in normal embryos, the Wnt pathway can enhance induction, by the SMAD 2 pathway, of the organizer genes Gsc and Chd but not the pan-mesodermal marker genes Xbra and Eomes. We conclude that the Wnt and SMAD2 signaling pathways cooperate to induce the expression of Spemann-organizer specific genes and so help to localize their spatial expression.

Hypotension and inflammatory cytokine gene expression triggered by factor Xa–nitric oxide signaling

The signaling pathway initiated by factor Xa on vascular endothelial cells was investigated. Factor Xa stimulated a 5- to 10-fold increased release of nitric oxide (NO) in a dose-dependent reaction (0.1–2.5 μg/ml) unaffected by the thrombin inhibitor hirudin but abolished by active site inhibitors, tick anticoagulant peptide, or Glu-Gly-Arg-chloromethyl ketone. In contrast, the homologous clotting protease factor IXa or another endothelial cell ligand, fibrinogen, was ineffective. A factor Xa inter-epidermal growth factor synthetic peptide L83FTRKL88(G) blocking ligand binding to effector cell protease receptor-1 inhibited NO release by factor Xa in a dose-dependent manner, whereas a control scrambled peptide KFTGRLL was ineffective. Catalytically active factor Xa induced hypotension in rats and vasorelaxation in the isolated rat mesentery, which was blocked by the NO synthase inhibitor l-NG-nitroarginine methyl ester (l-NAME) but not by d-NAME. Factor Xa/NO signaling also produced a dose-dependent endothelial cell release of interleukin 6 (range 0.55–3.1 ng/ml) in a reaction inhibited by l-NAME and by the inter-epidermal growth factor peptide Leu83–Leu88 but unaffected by hirudin. Maximal induction of interleukin 6 mRNA required a brief, 30-min stimulation with factor Xa, unaffected by subsequent addition of tissue factor pathway inhibitor. These data suggest that factor Xa-induced NO release modulates endothelial cell-dependent vasorelaxation and cytokine gene expression. This pathway requiring factor Xa binding to effector cell protease receptor-1 and a secondary step of ligand-dependent proteolysis may preserve an anti-thrombotic phenotype of endothelium but also trigger acute phase responses during activation of coagulation in vivo.