Insulin-like growth factor binding protein 2 (IGFBP2) promotes glioma development and progression

Overexpression of insulin-like growth factor binding protein 2 (IGFBP2) is associated with progression and poor survival in many types of human cancer (such as prostate, ovarian, adrenocortical, breast, colorectal carcinomas, leukemia, and high-grade gliomas). We therefore hypothesize that IGFBP2 is a key regulator of tumor progression. We tested our hypothesis in gliomas using the somatic gene transfer RCAS-tva mouse model system, which permits the introduction of specific genes into specific, cell lineages, in this case glial cells (RCAS: Replication competent avian sarcomavirus, tv-a: avian RCAS virus receptor). Mice are transgenic and harbor the tv-a receptor under the control of a glial-specific promoter and study genes are cloned into the RCAS vector for post-natal intracranial delivery. For these experiments, the study genes were IGFBP2, platelet-derived growth factor B (PDGFB), K-Ras, Akt, and IIp45 (invasion inhibitory protein 45 kDa; known to bind and block IGFBP2 activity), which were delivered separately and in combination. Our results show that PDGFB signaling leads exclusively to the formation of low-grade (WHO grade II) oligodendrogliomas. PDGFB delivered in combination with IGFBP2 results in the formation of anaplastic oligodendrogliomas (WHO grade III), which are characterized by increased cellularity, vascular proliferation, small regions of necrosis, increased mitotic activity, and increased activation of the Akt pathway. IIp45 injected in combination with PDGFB and IGFBP2 ablates IGFBP2-induced tumor progression, which results in formation of low-grade oligodendrogliomas, and an overall reduction in tumor incidence. K-Ras expression was required to form astrocytomas with either IGFBP2 or Akt, indicating the activation of two separate pathways is necessary for gliomagenesis. In ex vivo experiments, blockade of Akt by an inhibitor led to decreased viability of cells co-expressing IGFBP2 versus PDGFB expression alone. This study provides definitive evidence, for the first time, that: (1) IGFBP2 plays a role in activation of the Akt pathway, (2) IGFBP2 collaborates with K-Ras or PDGFB in the development and progression of two major types of glioma, and (3) IGFBP2-induced tumor progression can be ablated by IIp45 or by specific inhibition of the Akt pathway. ^

Establishment of a complex inflammatory and protumorigenic microenvironment in mouse pancreas through cyclooxygenase-2 overexpression

Chronic inflammation is an established risk factor in the pathogenesis of many cancers. Pancreatic ductal adenocarcinoma, a malignancy with a particularly dismal prognosis, is no exception. Cyclooxygenase-2, a key enzyme induced by tissue injury, has a critical role in the generation of bioactive lipids known as prostaglandins. COX-2 overexpression is a frequent finding in pancreatic cancer, chronic pancreatitis and pancreatic intraepithelial neoplasias. To explore mechanisms through which chronic inflammation establishes and maintains a protumorigenic environment, we designed a mouse model overexpressing COX-2 in pancreatic parenchyma (BK5.COX-2 mice). We discovered that constitutive expression of COX-2 has a number of important sequelae, including upregulation of additional eicosanoid-generating enzymes and proinflammatory cytokines. Many of these molecular alterations precede the onset of significant histopathological changes. Increased levels of prostaglandins E2, D2, and F2α, 5-, 12-, and 15-hydroxyeiosatetraenoic acid (HETEs) were documented in tumors and pancreata of younger transgenic mice. Using a TaqMan™ Mouse Immune Panel, we detected elevated mRNAs for a number of proinflammatory cytokines (e.g., TNFα, IL-1β, IL-6). ^ Histological examination revealed early changes in the pancreas with similarities to human chronic pancreatitis, including loss of acinar cells, appearance of metaplastic ducts, and increased deposition of stroma. As the lesions progress, features typical of dysplastic and neoplastic cells emerged within the metaplastic ductal complexes, including cellular and nuclear atypia, crowding of cells, and loss of normal tissue architecture. The amount of fibroinflammatory stroma increased considerably; numerous small vessels were evident. A number of immunocytes from both the myeloid and lymphoid lineages were identified in transgenic pancreata. Neutrophils were the earliest to infiltrate, followed shortly by macrophages and mast cells. B and T cells generally began to appear by 8–12 weeks, and organized aggregates of lymphoid cells were often found in advanced lesions. ^ We tested the efficacy of several chemopreventive agents in this model, including celecoxib, a COX-2 selective inhibitor, pentoxifylline, a cytokine inhibitor, curcumin, a polyphenol with antioxidant and anti-inflammatory properties, and GW2974, a dual EGFR/ErbB2 inhibitor. Effects on lesion development were modest in the GW2974 and pentoxifylline treated groups, but significant prevention effects were observed with curcumin and celecoxib. ^

XIPHOPHORUS ENZYME GENE EXPRESSION: TISSUE SPECIFICITY, GENETIC CONTROL, AND STABILITY IN CULTURED CELLS

Five permanent cell lines were developed from Xiphophorus maculatus, X. helleri, and their hybrids using three tissue sources, including adults and embryos of different stages. To evaluate cell line gene expression for retention of either tissue-of-origin-specific or ontogenetic stage-specific characters, the activity distribution of 44 enzyme loci was determined in 11 X. maculatus tissues, and the developmental genetics of 17 enzyme loci was charted in X. helleri and in helleri x maculatus hybrids using starch gel electrophoresis. In the process, eight new loci were discovered and characterized for Xiphophorus.^ No Xiphophorus cell line showed retention of tissue-of-origin-specific or ontogenetic stage-specific enzyme gene expressional traits. Instead, gene expression was similar among the cell lines. One enzyme, adenosine deaminase (ADA) was an exception. Two adult-origin cell lines expressed ADA, whereas, three cell lines derived independently from embryos did not. ADA('-) expression of Xiphophorus embryo-derived cell lines may represent retention of an embryonic gene expressional trait. In one cell line (T(,3)) derived from 13 pooled interspecific hybrid (F(,2)) embryos, shifts with time were observed at enzyme loci polymorphic between the two species. This suggested shifts in ratios of cells of different genotypes in the population rather than unstable gene expression in one dominant cell type.^ Verification of this hypothesis was attempted by cloning the culture--seeking clones having different genetic signatures. The large number of loci electrophoretically polymorphic between the two species and whose alleles segregated independently into the 13 progeny from which this culture originated almost guaranteed the presence of different genetic signatures (lineages) in T(,3).^ Seven lineages of cells were found within T(,3), each expressing genotypes at some loci not characteristic of the expression of the culture-as-a-whole, supporting the hypothesis tested. Quantitative studies of ADA expression in the whole culture (ADA('-)) and in clones of these seven lineages suggested the predominance in T(,3) of ADA deficient cell lineages, although moderate to high ADA output clones also occurred. Thus, T(,3) has the potential to shift phenotypes from ADA('-) to ADA('+) by simply changing proportions of its constituent cell types, demonstrating that such shifts can occur in any cell culture containing cells of mixed expressional characteristics.^

Molecular and evolutionary genetics of the X -linked visual pigment genes in humans and New World monkeys

Normal humans have one red and at least one green visual pigment genes. These genes are tightly linked as tandem repeats on the X chromosome and each of them has six exons. There is only one X-linked visual pigment gene in New World monkeys (NWMs) but the locus has three polymorphic alleles encoding red, yellow and green visual pigments, respectively. The spectral properties of the squirrel monkey and the marmoset (both NWMs) have been studied and partial sequences of the three alleles are available. To study the evolutionary history of these X-linked opsin genes in humans and NWMs, coding and intron sequences of the three squirrel monkey alleles and the three marmoset alleles were amplified by PCR followed by subcloning and sequencing. Introns 2 and 4 of the human red and green pigment genes were also sequenced. The results obtained are as follows: (1) The sequences of introns 2 and 4 of the human red and green opsin genes are significantly more similar between the two genes than are coding sequences, contrary to the usual situation where coding regions are better conserved in evolution than are introns. The high similarities in the two introns are probably due to recent gene conversion events during evolution of the human lineage. (2) Phylogenetic analysis of both intron and exon sequences indicates that the phylogenetic tree of the available primate opsin genes is the same as the species tree. The two human genes were derived from a gene duplication event after the divergence of the human and NWM lineages. The three alleles in each of the two NWM species diverged after the split of the two NWMs but have persisted in the population for at least 5 million years. (3) Allelic gene conversion might have occurred between the three squirrel monkey alleles. (4) A model of additive effect of hydroxyl-bearing amino acids on spectral tuning is proposed by treating some unknown variables as groups. Under the assumption that some residues have no effect, it is found that at least five amino acid residues, at positions 178 (3 nm), 180 (5 nm), 230 ($-$4 nm), 277 (9 nm) and 285 (13 nm), have linear spectral tuning effects. (5) Adaptive evolution of the opsin genes to different spectral peaks was observed at four residues that are important for spectral tuning. ^

Cloning and molecular analysis of paraxis: A {\it trans\/}-acting factor required for somite maturation

During vertebrate embryogenesis, cells from the paraxial mesoderm coalesce in a rostral-to-caudal progression to form the somites. Subsequent compartmentalization of the somites yields the sclerotome, myotome and dermatome, which give rise to the axial skeleton, axial musculature, and dermis, respectively. Recently, we cloned a novel basic-Helix-Loop-Helix (bHLH) protein, called scleraxis, which is expressed in the sclerotome, in mesenchymal precursors of bone and cartilage, and in connective tissues. This dissertation focuses on the cloning, expression and functional analysis of a bHLH protein termed paraxis, which is nearly identical to scleraxis within the bHLH region but diverges in both its amino and carboxyl termini. During the process of mouse embryogenesis, paraxis transcripts are first detected at about day 7.5 post coitum within the primitive mesoderm lying posterior to the head and heart primordia. Subsequently, paraxis expression progresses caudally through the paraxial mesoderm, immediately preceding somite formation. Paraxis is expressed at high levels in newly formed somites before the first detectable expression of the myogenic bHLH genes, and as the somite becomes compartmentalized, paraxis becomes downregulated within the myotome.^ To determine the function of paraxis during mammalian embryogenesis, mice were generated with a null mutation in the paraxis locus. Paraxis null mice survived until birth, but exhibited severe foreshortening along the anteroposterior axis due to the absence of vertebrae caudal to the midthoracic region. The phenotype also included axial skeletal defects, particularly shortened bifurcated ribs which were detached from the vertebral column, fused vertebrae and extensive truncation and disorganization caudal to the hindlimbs. Mutant neonates also lacked normal levels of trunk muscle and exhibited defects in the dermis as well as the stratification of the epidermis. Analysis of paraxis -/- mutant embryos has revealed a failure of the somites to both properly epithelialize and compartmentalize, resulting in defects in somite-derived cell lineages. These results suggest that paraxis is an essential component of the genetic pathway regulating somitogenesis. ^

Functional analysis of the sea urchin {\it orthodenticle\/} -related gene, {\it SPOTX\/}, in aboral ectoderm-specific gene activation and aboral ectoderm cell fate specification

An important question in developmental biology is how embryonic cell types are derived from a fertilized egg. To address this question, this thesis investigates the mechanisms by which the aboral ectoderm-specific Spec2a gene is spatially and temporally regulated during sea urchin embryogenesis. The Spec2a gene of the sea urchin Strongylocentratus purpuratus has served as a valuable maker to understand the basis of lineage-specific gene activation and the role of transcription factors in cell fate specification. The hypothesis is that transcription factors responsible for cell type-specific gene activation are key components in the initial cell specification step. The Spec2a gene, which encodes a small cytosolic calcium-binding protein, is expressed exclusively in aboral ectoderm cell lineages. The 1516-bp control region of the Spec2a gene contains a 188-bp enhancer element required for temporal activation and aboral ectoderm/mesenchyme cell expression, while an unidentified element upstream of the enhancer represses expression in mesenchyme cells. Using an enhancer activation assay, combined with site-directed mutagenesis, I showed that three TAATCC/T sites within the enhancer are responsible for enhancer activity. Mutagenizing these sites and a fourth one just upstream abolished all activity from the Spec2a control region. A 77-bp DNA fragment from the Spec2a enhancer containing two of the TAATCC/T sites is sufficient for aboral ectoderm/mesenchyme cell expression. A cDNA encoding SpOtx, an orthodenticle-related protein, was cloned from S. purpuratus and shown to bind with high affinity to the TAATCC/T sequences within the Spec2a control region. SpOtx transcripts were found initially in all cells of the cleaving embryo, but they gradually became restricted to oral ectoderm and endoderm cells, suggesting that SpOtx might play a role in the initial temporal activation of the Spec2a gene and most likely has additional functions in the developing embryo. To reveal the broader biological functions of SpOtx, I injected SpOtx mRNA into living sea urchin eggs to determine what effects overexpressing the SpOtx protein might have on embryo development. SpOtx mRNA-injected embryos displayed dramatic alterations in development. Instead of developing into pluteus larvae with 15 different cell types, uniform epithelia balls were formed. These balls consisted of a thin layer of squamous cells with short cilia highly reminiscent of aboral ectoderm. Immunohistochemical staining and RT-PCR demonstrated that the SpOtx-injected embryoids expressed aboral ectoderm markers uniformly, but showed very weak or no expression of markers for non-aboral ectoderm cell types. These data strongly suggested that overexpression of SpOtx redirected the normal fate of non-aboral ectoderm cells to that of aboral ectoderm. These results show that SpOtx is involved in aboral ectoderm differentiation by activating aboral ectoderm-specific genes and that modulating its expression can lead to changes in cell fate. ^

Functions of GCN5 and P /CAF acetyltransferases in mammalian development

Histone acetyltransferases are important chromatin modifiers that function as transcriptional co-activators. The identification of the transcriptional regulator GCN5 as the first nuclear histone acetyltransferase in yeast directly linked chromatin remodeling to transcriptional regulation. Although emerging evidence suggests that acetyltransferases participate in multiple cellular processes, their roles in mammalian development remain undefined. In this study, I have cloned and characterized the mouse homolog of GCN5 and a closely related protein P/CAF that interacts with p300/CBP. In contrast to yeast GCN5, but similar to P/CAF, mouse GCN5 possesses an additional N-terminal domain that confers the ability to acetylate nucleosomal histones. GCN5 and P/CAF exhibit identical substrate specificity and both interact with p300/CBP. Interestingly, expression levels of GCN5 and P/CAF display a complementary pattern in mouse embryos and in adult tissues, suggesting that they have distinct tissue or developmental stage specific roles. To define the in vivo function of GCN5 and P/CAF, I have generated mice that are nullizygous for GCN5 or P/CAF. P/CAF null mice are viable and fertile with no gross morphological defects, indicating that P/CAF is dispensable for development and p300/CBP function in vivo. In contrast, mice lacking GCN5 die between 10.5–11 days of gestation. GCN5 null mice are severely retarded but have anterior ectopic outgrowth. Molecular marker analyses reveal that early mesoderm is formed in GCN5 null mice but further differentiation into distinct mesodermal lineages is perturbed. While presomitic mesoderm and chodamesoderm are missing in GCN5 mutant mice, extraembryonic tissues and lateral mesoderm are unaffected. This is consistent with our finding that GCN5 expression is absent in the heart and extraembryonic tissues but is uniform throughout the rest of the embryo. Remarkably, GCN5 mutant mice exhibit an unusually high incidence of apoptosis in the embryonic ectoderm and mesoderm. Finally, mice doubly null for GCN5 and P/CAF die much earlier than mice harboring the GCN5 mutation alone, suggesting that P/CAF and GCN5 share some overlapping function during embryogenesis. This work is the first study to show that specific acetyltransferase is important for cell survival as well as mesoderm differentiation or maintenance during early mammalian development. ^

Isolation and characterization of a novel human Mix -like homeobox gene MIXL

Molecular events involved in specification of early hematopoietic system are not well known. In Xenopus, a paired-box homeodomain family (Mix.1–4) has been implicated in this process. Although Mix-like homeobox genes have been isolated from zebrafish (bon), chicken (CMIX) and mice (MmI/MIXL1), isolation of a human Mix-like gene has remained elusive. ^ We have recently isolated and characterized a novel human Mix-like homeobox gene with a predicted open reading frame of 232 amino acids designated the Mix.1 homeobox (Xenopus laevis)-like gene (MIXL). The overall identity of this novel protein to CMIX and MmI/MIXL1 is 41% and 69%, respectively. However, the identity in the homeodomain is 66% to that of Xenopus Mix.1, 79% to that of CMIX, and 94% to that of MmI/MIXL1. In normal hematopoiesis, MIXL expression appears to be restricted immature B and T lymphoid cells. Several acute leukemic cell lines of B, T and myeloid lineages express MIXL suggesting a survival/block in differentiation advantage. Furthermore, Xenopus animal cap assay revealed that MIXL could induce expression of the α-globin gene, suggesting a functional conservation of the homeodomain. ^ Biochemical analysis revealed that MIXL proteins are phosphorylated at multiple sites. Immunoprecipitation and immunoblotting confirmed that MIXL is tyrosine phosphorylated. Mutational analysis determined that Tyr20 appears to be the site for phosphorylation. However, deletion analysis preliminarily showed that the proline-rich domain appears not to be necessary for tyrosine phosphorylation. The novel finding will help us make a deeper understanding of the regulation on homeodomain proteins by rarely reported tyrosine phosphorylation. ^ Taken together, isolation of the MIXL gene is the first step toward understanding novel regulatory circuits in early hematopoietic differentiation and malignant transformation. ^

(Table S1) Molecular analyses and morphological comparisons of deep-sea clams

Large vesicomyid clams are common inhabitants of sulphidic deep-sea habitats such as hydrothermal vents, hydrocarbon seeps and whale-falls. Yet, the species- and genus-level taxonomy of these diverse clams has been unstable due to insufficiencies in sampling and absence of detailed taxonomic studies that would consistently compare molecular and morphological characters. To clarify uncertainties about species-level assignments, we examined DNA sequences from mitochondrial cytochrome-c-oxidase subunit I (COI) in conjunction with morphological characters. New and published COI sequences were used to create a molecular database for 44 unique evolutionary lineages corresponding to species. Overall, the congruence between molecular and morphological characters was good. Several discrepancies due to synonymous species designations were recognized, and acceptable species names were rectified with published COI sequences in cases where morphological specimens were available. We identified seven species with trans-Pacific distributions, and two species with Indo-Pacific distributions. Presently, 27 species have only been documented from one region, which might reflect limited ranges, or insufficient geographical sampling. Vesicomyids exhibit the greatest species diversity along the northwest Pacific ridge systems and in the eastern Pacific, along the western America margin, where depth zonation typically results in segregation of closely related species. The broad distributions of several vesicomyid species suggest that their required chemosynthetic habitats might be more common than previously recognized and occur along most continental margins.

Geochemistry of ash and pumice of ODP Hole 132-810C

Air-fall volcanic ash and pumice were recovered from 22 intervals in upper Miocene-Pleistocene nannofossil oozes cored in Hole 810C on Shatsky Rise, northwest Pacific. Shatsky Rise is near the eastern limit of ash falls produced by explosive volcanism in arc systems in northern Japan and the Kuriles, more than 1600 km away. Electron probe analyses establish that the ash beds and pumice pebbles are andesitic to rhyolitic in composition, and belong to both tholeiitic and high-alumina lineages similar to tephra from Japanese volcanoes. High-speed winds in the polar-front and subtropical jets are evidently what propelled the ash for such a distance. The pumice arrived by flotation, driven from the same directions by winds, waves, and currents. It is not ice-rafted debris from the north. One thick pumice bed probably was deposited when a large pumice mat passed over Shatsky Rise. Far more abundant ash occurs in sediments cored at DSDP Sites 578 through 580, about 500 km west of Shatsky Rise. Most of the ash and pumice at Shatsky Rise can be correlated with specific ash beds at 1, 2, or all 3 of these sites by interpolating to precisely determined magnetic reversal sequences in the cores. Most of the correlations are to thick ash layers (5.7 +/- 3.0 cm) at one or more sites. These must represent extremely large eruptions that spread ash over very wide areas. Whereas several of the thicker correlative ashes fell from elongate east-trending plumes directed from central Japan, the majority of them - dating from about 2 Ma - came from the North Honshu and Kurile arc systems to the northwest. This direction probably was in response to both long-term and seasonal fluctuations in the location and velocity of the polar-front jet, and to more vigorous winter storm fronts originating over glaciated Siberia.

Geochemical data from pore water and solids from cores GS08-153-91GC12 and GS08-153-53-GC6, Arctic Mid-Ocean Ridge

Microbial communities and their associated metabolic activity in marine sediments have a profound impact on global biogeochemical cycles. Their composition and structure are attributed to geochemical and physical factors, but finding direct correlations has remained a challenge. Here we show a significant statistical relationship between variation in geochemical composition and prokaryotic community structure within deep-sea sediments. We obtained comprehensive geochemical data from two gravity cores near the hydrothermal vent field Loki's Castle at the Arctic Mid-Ocean Ridge, in the Norwegian-Greenland Sea. Geochemical properties in the rift valley sediments exhibited strong centimeter-scale stratigraphic variability. Microbial populations were profiled by pyrosequencing from 15 sediment horizons (59,364 16S rRNA gene tags), quantitatively assessed by qPCR, and phylogenetically analyzed. Although the same taxa were generally present in all samples, their relative abundances varied substantially among horizons and fluctuated between Bacteria- and Archaea-dominated communities. By independently summarizing covariance structures of the relative abundance data and geochemical data, using principal components analysis, we found a significant correlation between changes in geochemical composition and changes in community structure. Differences in organic carbon and mineralogy shaped the relative abundance of microbial taxa. We used correlations to build hypotheses about energy metabolisms, particularly of the Deep Sea Archaeal Group, specific Deltaproteobacteria, and sediment lineages of potentially anaerobic Marine Group I Archaea. We demonstrate that total prokaryotic community structure can be directly correlated to geochemistry within these sediments, thus enhancing our understanding of biogeochemical cycling and our ability to predict metabolisms of uncultured microbes in deep-sea sediments.

(Table 1) Number of alleles per locus and heterozygosity of Arctic charr (Salvelinus alpinus) populations in Iceland

The ecological theory of adaptive radiation predicts that the evolution of phenotypic diversity within species is generated by divergent natural selection arising from different environments and competition between species. Genetic connectivity among populations is likely also to have an important role in both the origin and maintenance of adaptive genetic diversity. Our goal was to evaluate the potential roles of genetic connectivity and natural selection in the maintenance of adaptive phenotypic differences among morphs of Arctic charr, Salvelinus alpinus, in Iceland. At a large spatial scale, we tested the predictive power of geographic structure and phenotypic variation for patterns of neutral genetic variation among populations throughout Iceland. At a smaller scale, we evaluated the genetic differentiation between two morphs in Lake Thingvallavatn relative to historically explicit, coalescent-based null models of the evolutionary history of these lineages. At the large spatial scale, populations are highly differentiated, but weakly structured, both geographically and with respect to patterns of phenotypic variation. At the intralacustrine scale, we observe modest genetic differentiation between two morphs, but this level of differentiation is nonetheless consistent with strong reproductive isolation throughout the Holocene. Rather than a result of the homogenizing effect of gene flow in a system at migration-drift equilibrium, the modest level of genetic differentiation could equally be a result of slow neutral divergence by drift in large populations. We conclude that contemporary and recent patterns of restricted gene flow have been highly conducive to the evolution and maintenance of adaptive genetic variation in Icelandic Arctic charr.

The effect of repeated exposure to parasite genotypes of the eyefluke Diplostomum pseudospathaceum on infection rate and immunization in three-spined sticklebacks

Adaptive immunity in vertebrates can confer increased resistance against invading pathogens upon re-infection. But how specific parasite genotypes affect the transition from innate to adaptive immunity is poorly understood. Here, we investigated the effects of homologous and heterologous exposures of genetically distinct parasite lineages of the eye fluke Diplostomum pseudospathaceum on gene expression patterns of adaptive immunity in sticklebacks (Gasterosteus aculeatus). We showed that observable differences were largely attributable to final exposures and that there is no transcription pattern characteristic for a general response to repeated infections with D. pseudospathaceum. Final exposure did not unify expression patterns of heterologous pre-exposed fish. Interestingly, heterologous final exposures showed similarities between different treatment groups subjected to homologous pre-exposure. The observed pattern was supported by parasite infection rates and suggests that host immunization was optimized towards an adaptive immune response that favored effectiveness against parasite diversity over specificity.

Lupinus mariae-josphi, a new lupin endemic of soils with active lime and high pH in South Eastern Spain, is nodulated by a new bacterial lineage within Bradyrhizobium genus

Lupinus mariae-josephi is a recently described species (Pascual, 2004) able to grow in soils with high pH and active lime content in the Valencia province (Spain). L. mariae-josephi endosymbionts are extremely slowgrowing bacteria with genetic and symbiotic characteristics that differentiate them from Bradyrhizobium strains nodulating Lupinus spp. native of the Iberian Peninsula and adapted to grow in acid soils. Cross-inoculation experiments revealed that all the endosymbiotic isolates from L. mariae-josephi tested are legume-host selective and are unable to nodulate species such as L. angustifolius, and L. luteus. In contrast, Bradyrhizobium strains from Lupinus spp. tested were able to nodulate L. mariae-josephi, although the nodules fixed nitrogen inefficiently. Phylogenetic analysis was performed with housekeeping genes (rrn, glnII, recA, atpD) and nodulation gene nodC. Housekeeping gene phylogeny revealed that L. mariae-josephi rhizobia form a strongly supported monophyletic group within Bradyrhizobium genus. This cluster also includes B. jicamae and certain strains of B. elkanii. Contrarily, isolates from other Lupinus spp. native of the Iberian Peninsula were grouped mainly within B. canariense and two B. japonicum lineages. Phylogenetic analysis of L. mariae-josephi isolates based on the nodC symbiotic gene defined a solid clade close to isolates from Algerian Retama spp. and to fast-growing rhizobia.

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.