Implication du récepteur aux cannabinoïdes CB2 dans le développement de la voie rétinothalamique

Durant le développement du système visuel, les cellules ganglionnaires de la rétine (CGRs) envoient des axones qui seront influencés par divers signaux guidant leur cône de croissance, permettant ainsi la navigation des axones vers leurs cibles terminales. Les endocannabinoïdes, des dérivés lipidiques activant les récepteurs aux cannabinoides (CB1 et CB2), sont présents de manière importante au cours du développement. Nous avons démontré que le récepteur CB2 est exprimé à différents points du tractus visuel durant le développement du hamster. L’injection d’agonistes et d’agonistes inverses pour le récepteur CB2 a modifié l’aire du cône de croissance et le nombre de filopodes présents à sa surface. De plus, l’injection d’un gradient d’agoniste du récepteur CB2 produit la répulsion du cône de croissance tandis qu’un analogue de l’AMPc (db-AMPc) produit son attraction. Les effets du récepteur CB2 sur le cône de croissance sont produits en modulant l’activité de la protéine kinase A(PKA), influençant la présence à la membrane cellulaire d’un récepteur à la nétrine-1 nommé Deleted in Colorectal Cancer (DCC). Notamment, pour que le récepteur CB2 puisse moduler le guidage du cône de croissance, la présence fonctionnelle du récepteur DCC est essentielle.. Suite à une injection intra-occulaire d’un agoniste inverse du récepteur CB2, nous avons remarqué une augmentation de la longueur des branches collatérales des axones rétiniens au niveau du LTN (noyau lateral terminal). Nous avons également remarqué une diminution de la ségrégation des projections ganglionnaires au niveau du dLGN, le noyau genouillé lateral dorsal, chez les animaux transgéniques cnr2-/-, ayant le gène codant pour le récepteur CB2 inactif. Nos données suggèrent l’implication des endocannabinoïdes et de leur récepteur CB2 dans la modulation des processus de navigation axonale et de ségrégation lors du développement du système visuel.

Le rôle des récepteurs aux cannabinoïdes CB1 et CB2 dans le guidage axonal

Au cours du développement, les axones des cellules ganglionnaires de la rétine (CGRs) voyagent sur de longues distances pour établir des connexions avec leurs cellules cibles. La navigation des cônes de croissance est guidée par différentes molécules chimiotropiques présentes dans leur environnement. Les endocannabinoïdes (eCB) sont d’importants neuromodulateurs qui régulent de manière rétrograde la fonction de nombreuses synapses du cerveau. Ils agissent principalement par le biais de leurs récepteurs liés à une protéine Gi/o CB1 (CB1R) et CB2 (CB2R). La présence des eCBs durant le stade fœtal et la période postnatale suggère leur implication dans des événements régulant le développement du système nerveux. Cette thèse confirme l’expression des récepteurs aux cannabinoïdes CB1 et CB2 ainsi que l’enzyme dégradant les eCBs lors du développement embryonnaire et perinatal des CGRs et de la voie rétinothalamique in vivo. La manipulation pharmacologique de l’activité de CB1R et CB2R réorganise la morphologie du cône de croissance des CGRs et des neurones corticaux in vitro. De plus, la stimulation locale avec un agoniste de CB1R ou de CB2R modifie le comportement du cône de croissance entraînant sa répulsion. CB1R et CB2R modulent par le biais de la voie de signalisation AMPc/PKA, la mobilisation de DCC à la membrane plasmique. Par ailleurs, les résultats de cette recherche démontrent également l’implication de CB1R et CB2R dans la ségrégation des projections ipsi- et controlatérales et le développement de la voie rétinothalamique.

Organisation and tyrosine hydroxylase and calretinin immunoreactivity in the main olfactory bulb of paca (cuniculus paca): a large caviomorph rodent

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

A versatile toolkit to produce sensitive FRET biosensors to visualize signaling in time and space

Genetically encoded, ratiometric biosensors based on fluorescence resonance energy transfer (FRET) are powerful tools to study the spatiotemporal dynamics of cell signaling. However, many biosensors lack sensitivity. We present a biosensor library that contains circularly permutated mutants for both the donor and acceptor fluorophores, which alter the orientation of the dipoles and thus better accommodate structural constraints imposed by different signaling molecules while maintaining FRET efficiency. Our strategy improved the brightness and dynamic range of preexisting RhoA and extracellular signal-regulated protein kinase (ERK) biosensors. Using the improved RhoA biosensor, we found micrometer-sized zones of RhoA activity at the tip of F-actin bundles in growth cone filopodia during neurite extension, whereas RhoA was globally activated throughout collapsing growth cones. RhoA was also activated in filopodia and protruding membranes at the leading edge of motile fibroblasts. Using the improved ERK biosensor, we simultaneously measured ERK activation dynamics in multiple cells using low-magnification microscopy and performed in vivo FRET imaging in zebrafish. Thus, we provide a construction toolkit consisting of a vector set, which enables facile generation of sensitive biosensors.

Large scale multiphysics modeling of neurites

Situado en el límite entre Ingeniería, Informática y Biología, la mecánica computacional de las neuronas aparece como un nuevo campo interdisciplinar que potencialmente puede ser capaz de abordar problemas clínicos desde una perspectiva diferente. Este campo es multiescala por naturaleza, yendo desde la nanoescala (como, por ejemplo, los dímeros de tubulina) a la macroescala (como, por ejemplo, el tejido cerebral), y tiene como objetivo abordar problemas que son complejos, y algunas veces imposibles, de estudiar con medios experimentales. La modelización computacional ha sido ampliamente empleada en aplicaciones Neurocientíficas tan diversas como el crecimiento neuronal o la propagación de los potenciales de acción compuestos. Sin embargo, en la mayoría de los enfoques de modelización hechos hasta ahora, la interacción entre la célula y el medio/estímulo que la rodea ha sido muy poco explorada. A pesar de la tremenda importancia de esa relación en algunos desafíos médicos—como, por ejemplo, lesiones traumáticas en el cerebro, cáncer, la enfermedad del Alzheimer—un puente que relacione las propiedades electrofisiológicas-químicas y mecánicas desde la escala molecular al nivel celular todavía no existe. Con ese objetivo, esta investigación propone un marco computacional multiescala particularizado para dos escenarios respresentativos: el crecimiento del axón y el acomplamiento electrofisiológicomecánico de las neuritas. En el primer caso, se explora la relación entre los constituyentes moleculares del axón durante su crecimiento y sus propiedades mecánicas resultantes, mientras que en el último, un estímulo mecánico provoca deficiencias funcionales a nivel celular como consecuencia de sus alteraciones electrofisiológicas-químicas. La modelización computacional empleada en este trabajo es el método de las diferencias finitas, y es implementada en un nuevo programa llamado Neurite. Aunque el método de los elementos finitos es también explorado en parte de esta investigación, el método de las diferencias finitas tiene la flexibilidad y versatilidad necesaria para implementar mode los biológicos, así como la simplicidad matemática para extenderlos a simulaciones a gran escala con un coste computacional bajo. Centrándose primero en el efecto de las propiedades electrofisiológicas-químicas sobre las propiedades mecánicas, una versión adaptada de Neurite es desarrollada para simular la polimerización de los microtúbulos en el crecimiento del axón y proporcionar las propiedades mecánicas como función de la ocupación de los microtúbulos. Después de calibrar el modelo de crecimiento del axón frente a resultados experimentales disponibles en la literatura, las características mecánicas pueden ser evaluadas durante la simulación. Las propiedades mecánicas del axón muestran variaciones dramáticas en la punta de éste, donde el cono de crecimiento soporta las señales químicas y mecánicas. Bansándose en el conocimiento ganado con el modelo de diferencias finitas, y con el objetivo de ir de 1D a 3D, este esquema preliminar pero de una naturaleza innovadora allana el camino a futuros estudios con el método de los elementos finitos. Centrándose finalmente en el efecto de las propiedades mecánicas sobre las propiedades electrofisiológicas- químicas, Neurite es empleado para relacionar las cargas mecánicas macroscópicas con las deformaciones y velocidades de deformación a escala microscópica, y simular la propagación de la señal eléctrica en las neuritas bajo carga mecánica. Las simulaciones fueron calibradas con resultados experimentales publicados en la literatura, proporcionando, por tanto, un modelo capaz de predecir las alteraciones de las funciones electrofisiológicas neuronales bajo cargas externas dañinas, y uniendo lesiones mecánicas con las correspondientes deficiencias funcionales. Para abordar simulaciones a gran escala, aunque otras arquitecturas avanzadas basadas en muchos núcleos integrados (MICs) fueron consideradas, los solvers explícito e implícito se implementaron en unidades de procesamiento central (CPU) y unidades de procesamiento gráfico (GPUs). Estudios de escalabilidad fueron llevados acabo para ambas implementaciones mostrando resultados prometedores para casos de simulaciones extremadamente grandes con GPUs. Esta tesis abre la vía para futuros modelos mecánicos con el objetivo de unir las propiedades electrofisiológicas-químicas con las propiedades mecánicas. El objetivo general es mejorar el conocimiento de las comunidades médicas y de bioingeniería sobre la mecánica de las neuronas y las deficiencias funcionales que aparecen de los daños producidos por traumatismos mecánicos, como lesiones traumáticas en el cerebro, o enfermedades neurodegenerativas como la enfermedad del Alzheimer. ABSTRACT Sitting at the interface between Engineering, Computer Science and Biology, Computational Neuron Mechanics appears as a new interdisciplinary field potentially able to tackle clinical problems from a new perspective. This field is multiscale by nature, ranging from the nanoscale (e.g., tubulin dimers) to the macroscale (e.g., brain tissue), and aims at tackling problems that are complex, and sometime impossible, to study through experimental means. Computational modeling has been widely used in different Neuroscience applications as diverse as neuronal growth or compound action potential propagation. However, in the majority of the modeling approaches done in this field to date, the interactions between the cell and its surrounding media/stimulus have been rarely explored. Despite of the tremendous importance of such relationship in several medical challenges—e.g., traumatic brain injury (TBI), cancer, Alzheimer’s disease (AD)—a bridge between electrophysiological-chemical and mechanical properties of neurons from the molecular scale to the cell level is still lacking. To this end, this research proposes a multiscale computational framework particularized for two representative scenarios: axon growth and electrophysiological-mechanical coupling of neurites. In the former case, the relation between the molecular constituents of the axon during its growth and its resulting mechanical properties is explored, whereas in the latter, a mechanical stimulus provokes functional deficits at cell level as a consequence of its electrophysiological-chemical alterations. The computational modeling approach chosen in this work is the finite difference method (FDM), and was implemented in a new program called Neurite. Although the finite element method (FEM) is also explored as part of this research, the FDM provides the necessary flexibility and versatility to implement biological models, as well as the mathematical simplicity to extend them to large scale simulations with a low computational cost. Focusing first on the effect of electrophysiological-chemical properties on the mechanical proper ties, an adaptation of Neurite was developed to simulate microtubule polymerization in axonal growth and provide the axon mechanical properties as a function of microtubule occupancy. After calibrating the axon growth model against experimental results available in the literature, the mechanical characteristics can be tracked during the simulation. The axon mechanical properties show dramatic variations at the tip of the axon, where the growth cone supports the chemical and mechanical signaling. Based on the knowledge gained from the FDM scheme, and in order to go from 1D to 3D, this preliminary yet novel scheme paves the road for future studies with FEM. Focusing then on the effect of mechanical properties on the electrophysiological-chemical properties, Neurite was used to relate macroscopic mechanical loading to microscopic strains and strain rates, and simulate the electrical signal propagation along neurites under mechanical loading. The simulations were calibrated against experimental results published in the literature, thus providing a model able to predict the alteration of neuronal electrophysiological function under external damaging load, and linking mechanical injuries to subsequent acute functional deficits. To undertake large scale simulations, although other state-of-the-art architectures based on many integrated cores (MICs) were considered, the explicit and implicit solvers were implemented for central processing units (CPUs) and graphics processing units (GPUs). Scalability studies were done for both implementations showing promising results for extremely large scale simulations with GPUs. This thesis opens the avenue for future mechanical modeling approaches aimed at linking electrophysiological- chemical properties to mechanical properties. Its overarching goal is to enhance the bioengineering and medical communities knowledge on neuronal mechanics and functional deficits arising from damages produced by direct mechanical insults, such as TBI, or neurodegenerative evolving illness, such as AD.

A conventional myosin motor drives neurite outgrowth

Neuritic outgrowth is a striking example of directed motility, powered through the actions of molecular motors. Members of the myosin superfamily of actin-associated motors have been implicated in this complex process. Although conventional myosin II is known to be present in neurons, where it is localized at the leading edge of growth cones and in the cell cortex close to the plasma membrane, its functional involvement in growth cone motility has remained unproven. Here, we show that antisense oligodeoxyribonucleotides, complementary to a specific isoform of conventional myosin (myosin IIB), attenuate filopodial extension whereas sense and scrambled control oligodeoxyribonucleotides have no effect. Attenuation is shown to be reversible, neurite outgrowth being restored after cessation of the antisense regimen. Myosin IIB mRNA was present during active neurite extension, but levels were minimal in phenotypically rounded cells before neurite outgrowth and message levels decreased during antisense treatment. By contrast, the myosin IIA isoform is shown to be expressed constitutively both before and during neurite outgrowth and throughout exposure to myosin IIB antisense oligodeoxyribonucleotides. These results provide direct evidence that a conventional two-headed myosin is required for growth cone motility and is responsible, at least in part, for driving neuritic process outgrowth.

The cell surface metalloprotease/disintegrin Kuzbanian is required for axonal extension in Drosophila

It has long been suspected that proteolytic activity associated with advancing growth cones may be required for axon extension. We have isolated mutations in the kuzbanian (kuz) gene, which is expressed in the nervous system and encodes a putative zinc metalloprotease with a disintegrin domain. Drosophila embryos with loss-of-function mutations in kuz have dramatic defects in the development of central nervous system axon pathways, with many axons stalling and failing to extend through the nerve cord. This phenotype is rescued by panneural expression of kuz mRNA in the embryo. These results show that the Kuz metalloprotease is required for axon extension, suggesting a requirement for proteolytic activity at the growth cone surface.

Activation of RhoA by Lysophosphatidic Acid and Gα12/13 Subunits in Neuronal Cells: Induction of Neurite Retraction

Neuronal cells undergo rapid growth cone collapse, neurite retraction, and cell rounding in response to certain G protein–coupled receptor agonists such as lysophosphatidic acid (LPA). These shape changes are driven by Rho-mediated contraction of the actomyosin-based cytoskeleton. To date, however, detection of Rho activation has been hampered by the lack of a suitable assay. Furthermore, the nature of the G protein(s) mediating LPA-induced neurite retraction remains unknown. We have developed a Rho activation assay that is based on the specific binding of active RhoA to its downstream effector Rho-kinase (ROK). A fusion protein of GST and the Rho-binding domain of ROK pulls down activated but not inactive RhoA from cell lysates. Using GST-ROK, we show that in N1E-115 neuronal cells LPA activates endogenous RhoA within 30 s, concomitant with growth cone collapse. Maximal activation occurs after 3 min when neurite retraction is complete and the actin cytoskeleton is fully contracted. LPA-induced RhoA activation is completely inhibited by tyrosine kinase inhibitors (tyrphostin 47 and genistein). Activated Gα12 and Gα13 subunits mimic LPA both in activating RhoA and in inducing RhoA-mediated cytoskeletal contraction, thereby preventing neurite outgrowth. We conclude that in neuronal cells, LPA activates RhoA to induce growth cone collapse and neurite retraction through a G12/13-initiated pathway that involves protein-tyrosine kinase activity.

Real-Time Imaging of the Axonal Transport of Granules Containing a Tissue Plasminogen Activator/Green Fluorescent Protein Hybrid

A hybrid protein, tPA/GFP, consisting of rat tissue plasminogen activator (tPA) and green fluorescent protein (GFP) was expressed in PC12 cells and used to study the distribution, secretory behavior, and dynamics of secretory granules containing tPA in living cells with a neuronal phenotype. High-resolution images demonstrate that tPA/GFP has a growth cone-biased distribution in differentiated cells and that tPA/GFP is transported in granules of the regulated secretory pathway that colocalize with granules containing secretogranin II. Time-lapse images of secretion reveal that secretagogues induce substantial loss of cellular tPA/GFP fluorescence, most importantly from growth cones. Time-lapse images of the axonal transport of granules containing tPA/GFP reveal a surprising complexity to granule dynamics. Some granules undergo canonical fast axonal transport; others move somewhat more slowly, especially in highly fluorescent neurites. Most strikingly, granules traffic bidirectionally along neurites to an extent that depends on granule accumulation, and individual granules can reverse their direction of motion. The retrograde component of this bidirectional transport may help to maintain cellular homeostasis by transporting excess tPA/GFP back toward the cell body. The results presented here provide a novel view of the axonal transport of secretory granules. In addition, the results suggest that tPA is targeted for regulated secretion from growth cones of differentiated cells, strategically positioning tPA to degrade extracellular barriers or to activate other barrier-degrading proteases during axonal elongation.

Cloning, expression, and genetic mapping of Sema W, a member of the semaphorin family

The semaphorins comprise a large family of membrane-bound and secreted proteins, some of which have been shown to function in axon guidance. We have cloned a transmembrane semaphorin, Sema W, that belongs to the class IV subgroup of the semaphorin family. The mouse and rat forms of Sema W show 97% amino acid sequence identity with each other, and each shows about 91% identity with the human form. The gene for Sema W is divided into 15 exons, up to 4 of which are absent in the human cDNAs that we sequenced. Unlike many other semaphorins, Sema W is expressed at low levels in the developing embryo but was found to be expressed at high levels in the adult central nervous system and lung. Functional studies with purified membrane fractions from COS7 cells transfected with a Sema W expression plasmid showed that Sema W has growth-cone collapse activity against retinal ganglion-cell axons, indicating that vertebrate transmembrane semaphorins, like secreted semaphorins, can collapse growth cones. Genetic mapping of human SEMAW with human/hamster radiation hybrids localized the gene to chromosome 2p13. Genetic mapping of mouse Semaw with mouse/hamster radiation hybrids localized the gene to chromosome 6, and physical mapping placed the gene on bacteria artificial chromosomes carrying microsatellite markers D6Mit70 and D6Mit189. This localization places Semaw within the locus for motor neuron degeneration 2, making it an attractive candidate gene for this disease.

F-spondin is a contact-repellent molecule for embryonic motor neurons

The floor plate plays a key role in patterning axonal trajectory in the embryonic spinal cord by providing both long-range and local guidance cues that promote or inhibit axonal growth toward and across the ventral midline of the spinal cord, thus acting as an intermediate target for a number of crossing (commissural) and noncrossing (motor) axons. F-spondin, a secreted adhesion molecule expressed in the embryonic floor plate and the caudal somite of birds, plays a dual role in patterning the nervous system. It promotes adhesion and outgrowth of commissural axons and inhibits adhesion of neural crest cells. In the current study, we demonstrate that outgrowth of embryonic motor axons also is inhibited by F-spondin protein in a contact-repulsion fashion. Three independent lines of evidence support our hypothesis: substrate-attached F-spondin inhibits outgrowth of dissociated motor neurons in an outgrowth assay; F-spondin elicits acute growth cone collapse when applied to cultured motor neurons; and challenging ventral spinal cord explants with aggregates of HEK 293 cells expressing F-spondin, causes contact-repulsion of motor neurites. Structural–functional studies demonstrate that the processed carboxyl-half protein that contains the thrombospondin type 1 repeats is more prominent in inhibiting outgrowth, suggesting that the processing of F-spondin is important for enhancing its inhibitory activity.

Human semaphorins A(V) and IV reside in the 3p21.3 small cell lung cancer deletion region and demonstrate distinct expression patterns.

Semaphorins and collapsins make up a family of conserved genes that encode nerve growth cone guidance signals. We have identified two additional members of the human semaphorin family [human semaphorin A(V) and human semaphorin IV] in chromosome region 3p21.3, where several small cell lung cancer (SCLC) cell lines exhibit homozygous deletions indicative of a tumor suppressor gene. Human semaphorin A(V) has 86% amino acid homology with murine semaphorin A, whereas semaphorin IV is most closely related to murine semaphorin E, with 50% homology. These semaphorin genes are approximately 70 kb apart flanking two GTP-binding protein genes, GNAI-2 and GNAT-1. In contrast, other human semaphorin gene sequences (human semaphorin III and homologues of murine semaphorins B and C) are not located on chromosome 3. Human semaphorin A(V) is translated in vitro into a 90-kDa protein, which accumulates at the endoplasmic reticulum. The human semaphorin A(V) (3.4-kb mRNA) and IV (3.9- and 2.9-kb mRNAs) genes are expressed abundantly but differentially in a variety of human neural and nonneural tissues. Human semaphorin A(V) was expressed in only 1 out of 23 SCLCs and 7 out of 16 non-SCLCs, whereas semaphorin IV was expressed in 19 out of 23 SCLCs and 13 out of 16 non-SCLCs. Mutational analysis in semaphorin A(V) revealed mutations (germ line in one case) in 3 of 40 lung cancers. Our data suggest the need to determine the function of human semaphorins A(V) and IV in nonneural tissues and their role in the pathogenesis of lung cancer.

Antibodies against the T61 antigen inhibit neuronal migration in the chick optic tectum.

Cell migration in the central nervous system depends, in part, on receptors and extracellular matrix molecules that likewise support axonal outgrowth. We have investigated the influence of T61, a monoclonal antibody that has been shown to inhibit growth cone motility in vitro, on neuronal migration in the developing optic tectum. Intraventricular injections of antibody-producing hybridoma cells or ascites fluid were used to determine the action of this antibody in an in vivo environment. To document alterations in tectal layer formation, a combination of cell-nuclei staining and axonal immunolabeling methods was employed. In the presence of T61 antibody, cells normally destined for superficial layers accumulated in the ventricular zone instead, leading to a reduction of the cell-dense layer in the tectal plate. Experiments with 5-bromo-2'-deoxyuridine labeling followed by antibody staining confirmed that the nonmigrating cells remaining in the ventricular zone were postmitotic and had differentiated. The structure of radial glial cells, as judged by staining with a glia-specific antibody and the fluorescent tracer 1,1'-dioctadecyl-3,3,3',3'-tetramethylindocarbocyanine perchlorate (DiI), remained intact in these embryos. Our findings suggest that the T61 epitope is involved in a mechanism underlying axonal extension and neuronal migration, possibly by influencing the motility of the leading process.

Repelling class discrimination: ephrin-A5 binds to and activates EphB2 receptor signaling

The interactions between Eph receptor tyrosine kinases and their ephrin ligands regulate cell migration and axon pathfinding. The EphA receptors are generally thought to become activated by ephrin-A ligands, whereas the EphB receptors interact with ephrin-B ligands. Here we show that two of the most widely studied of these molecules, EphB2 and ephrin-A5, which have never been described to interact with each other, do in fact bind one another with high affinity. Exposure of EphB2-expressing cells to ephrin-A5 leads to receptor clustering, autophosphorylation and initiation of downstream signaling. Ephrin-A5 induces EphB2-mediated growth cone collapse and neurite retraction in a model system. We further show, using X-ray crystallography, that the ephrin-A5-EphB2 complex is a heterodimer and is architecturally distinct from the tetrameric EphB2-ephrin-B2 structure. The structural data reveal the molecular basis for EphB2-ephrin-A5 signaling and provide a framework for understanding the complexities of functional interactions and crosstalk between A- and B-subclass Eph receptors and ephrins.

Involvement of Islet-2 in the Slit signaling for axonal branching and defasciculation of the sensory neurons in embryonic zebrafish

In Drosophila melanogaster, Slit acts as a repulsive cue for the growth cones of the commissural axons which express a receptor for Slit, Roundabout (Robo), thus preventing the commissural axons from crossing the midline multiple times. Experiments using explant culture have shown that vertebrate Slit homologues also act repulsively for growth cone navigation and neural migration, and promote branching and elongation of sensory axons. Here, we demonstrate that overexpression of Slit2 in vivo in transgenic zebrafish embryos severely affected the behavior of the commissural reticulospinal neurons (Mauthner neurons), promoted branching of the peripheral axons of the trigeminal sensory ganglion neurons, and induced defasciculation of the medial longitudinal fascicles. In addition, Slit2 overexpression caused defasciculation and deflection of the central axons of the trigeminal sensory ganglion neurons from the hindbrain entry point. The central projection was restored by either functional repression or mutation of Robo2, supporting its role as a receptor mediating the Slit signaling in vertebrate neurons. Furthermore, we demonstrated that Islet-2, a LIM/homeodomain-type transcription factor, is essential for Slit2 to induce axonal branching of the trigeminal sensory ganglion neurons, suggesting that factors functioning downstream of Islet-2 are essential for mediating the Slit signaling for promotion of axonal branching. (C) 2004 Elsevier Ireland Ltd. All rights reserved.