Évaluation neurobiologique des souris spontanément hypertendues : Du vieillissement à la génomique comparative

Le but de cette thèse est premièrement d’évaluer l’effet du vieillissement sur les fonctions psychomotrices des souches de souris sélectionnées génétiquement en fonction de leur tension artérielle (TA); deuxièmement, de localiser les déterminants génétiques des phénotypes psychophysiologiques à partir de souches recombinantes congéniques (RCS). Ces travaux ont mené à la publication de 4 articles. Le premier article décrit l’évaluation des fonctions psychomotrices des souches avec une tension artérielle élevée (HBP), basse (LBP) et normale (NBP). La performance aux épreuves d’exploration, d’habiletés motrices et d’apprentissage spatial, a été mesurée sur deux cohortes âgées respectivement de 12 mois et de trois mois. Indépendamment de l’âge, les HBPs sont hyperactives dans l’open-field (OF), mais pas dans le test d’exploration de trous. Inversement, les LBP explorent moins d’espaces que les NBP et, à trois mois seulement, sont hypoactives dans l’OF. Par ailleurs, les HBPs et les LBP présentent des déficits précoces de coordination motrice et des fonctions visuo-motrices. Le second article concerne l’évaluation longitudinale de la coordination motrice, de l’anxiété et de l’apprentissage spatial des souches HBP, LBP et NBP, à l’âge de deux mois et de 12 mois. Le vieillissement accentue l’hyperactivité des HBPs dans l’OF. Par contre, l’hypoactivité des souris LBP est détectable seulement à l’âge de deux mois. Indépendamment de l’âge, les souris HBP et LBP montrent une perception réduite du danger dans l’épreuve d’anxiété et des dysfonctions visuo-motrices au labyrinthe aquatique. Enfin, des déficits précoces de coordination motrice se manifestent seulement chez les HBPs. Il reste à déterminer si les déficits observés sont liés à des déterminants génétiques indépendants ou secondaires aux altérations de la tension artérielle. Le troisième article présente la comparaison entre les souches consanguines A/J et C57Bl/6J (B6) aux épreuves de l’OF, de la planche à trous, du labyrinthe aquatique et du cintre (coordination motrice). Les B6 explore d’avantage l’OF et la planche à trous. Les B6 sont moins rapides sur le cintre, mais supérieurs aux A/J dans le labyrinthe aquatique, avec une plate-forme invisible ou visible. Ces résultats démontrent l’implication de déterminants génétiques. Cette thèse se termine par un quatrième article sur la localisation des déterminants génétiques de la susceptibilité au stress dans les RCS, dérivées de A/J et B6, et présentant un agencement spécifique de 12.5% du génome. La réactivité émotionnelle est évaluée dans l’OF et le plus-maze; la réponse de stress est mesurée par radio télémétrie de la température interne pendant le stress d’immobilisation (SI) sous diète régulière et riche en sel; l’excrétion des électrolytes urinaires est dosée après 24 heures de diète salée. Les loci les plus significatifs sont situés dans les régions suivantes: de l’émotionalité dans l’OF (Emo1) sur le chr. 1 (LOD=4.6) correspondant à la région homologue impliquée dans la cohorte d’hypertension familiale du Saguenay; de la dopa décarboxylase (ddc) sur le chr. 11 pour l’émergence du plus-maze (LOD=4.7); de la protéine liant l’endotoxine (lbp) sur le chr. 2 pour l’hypothermie initiale en réponse au SI (LOD=4); et de HSP90 sur le chr. 12 pour l’excrétion de Ca++ (LOD=4.6). Des banques de données sont ensuite interrogées pour recenser les polymorphismes des régions régulatrices ou codantes des gènes candidats chez les souches ancestrales A/J et B6, dont les séquences sont disponibles pour le génome entier. Des utilitaires web permettent de dévoiler les changements dans la structure secondaire de l’ARNm, l’interférence avec des microARN ou avec d’autres motifs de liaison. Plusieurs SNPs fonctionnels ont été identifiés pour le QTL du chr. 1, particulièrement dans les éléments de régulation; ceux-ci impliquant des gènes reliés avec les réponses inflammatoire/immunitaire ou avec le système cardiovasculaire. La quantification par la PCR confirme une régulation à la baisse d’atp1a2 dans le cœur et le cerveau des souches susceptibles à l’anxiété. Ces résultats confirment l’intrication des altérations de la susceptibilité au stress et de la régulation de la TA.

Régulation de la proprotéine convertase subtilisine / kexine type 9 (PCSK9) dans les cellules intestinales Caco-2/15

La proprotéine convertase subtilisine/kexine type 9 (PCSK9) favorise la dégradation post-transcriptionnelle du récepteur des lipoprotéines de faible densité (LDLr) dans les hépatocytes et augmente le LDL-cholestérol dans le plasma. Cependant, il n’est pas clair si la PCSK9 joue un rôle dans l’intestin. Dans cette étude, nous caractérisons les variations de la PCSK9 et du LDLr dans les cellules Caco-2/15 différentiées en fonction d’une variété d’effecteurs potentiels. Le cholestérol (100 µM) lié à l’albumine ou présenté en micelles a réduit de façon significative l’expression génique (30%, p<0,05) et l’expression protéique (50%, p<0,05) de la PCSK9. Étonnamment, une diminution similaire dans le LDLr protéique a été enregistrée (45%, p<0,05). Les cellules traitées avec le 25-hydroxycholestérol (50 µM) présentent également des réductions significatives dans l’ARNm (37%, p<0,01) et la protéine (75%, p<0,001) de la PCSK9. Une baisse des expressions génique (30%, p<0,05) et protéique (57%, p<0,01) a également été constatée dans le LDLr. Des diminutions ont aussi été observées pour la HMG CoA réductase et la protéine liant l’élément de réponse aux stérols SREBP-2. Il a été démontré que le SREBP-2 peut activer transcriptionnellement la PCSK9 par le biais de la liaison de SREBP-2 à son élément de réponse aux stérols situé dans la région proximale du promoteur de la PCSK9. Inversement, la déplétion du contenu cellulaire en cholestérol par l’hydroxypropyl-β-cyclodextrine a augmenté l’expression génique de la PCSK9 (20%, p<0,05) et son contenu protéique (540%, p<0,001), en parallèle avec les niveaux protéiques de SREBP-2. L’ajout des acides biliaires taurocholate et déoxycholate dans le milieu apical des cellules intestinales Caco-2/15 a provoqué une baisse d’expression génique (30%, p<0,01) et une hausse d’expression protéique (43%, p<0,01) de la PCSK9 respectivement, probablement via la modulation du FXR (farnesoid X receptor). Ces données combinées semblent donc indiquer que la PCSK9 fonctionne comme un senseur de stérols dans le petit intestin.

Caractérisation et étude d'un élément régulateur du gène codant pour le récepteur à la vasopressine de type 2

Thèse réalisée en cotutelle avec l'Université Pierre et Marie Curie, Paris VI, France

La dynamique chromatinienne induite par le pic de LH dans les cellules de granulosa chez la souris

La régulation transcriptionnelle des gènes est un processus indispensable sans lequel la diversité phénotypique des cellules ainsi que l’adaptation à leur environnement serait inexistant. L’identification des éléments de régulation dans le génome est d’une importance capitale afin de comprendre les mécanismes gouvernant l’expression des gènes spécifiques à un type cellulaire donné. Ainsi, suite au pic de LH, le follicule ovarien entre dans un programme intensif de différentiation cellulaire, orchestré par des modifications majeures du profile transcriptionnel des cellules de granulosa, déclenchant ultimement l’ovulation et la lutéinisation, processus indispensables à la fertilité femelle. L’hypothèse supportée par cette étude stipule qu’une réorganisation de la structure chromatinienne survient aux régions régulatrices d’une panoplie de gènes dans les heures suivant le pic de LH et qu’en isolant et identifiant ces régions, il serait possible de retrouver des éléments essentiels aux processus d’ovulation et de lutéinisation. Ainsi, en utilisant un protocole standard de superovulation chez la souris, les éléments de régulation se modifiant 4h suivant l’administration de hCG ont été isolés et identifiés dans les cellules de granulosa en utilisant la méthode FAIRE (Formaldehyde-Assisted Isolation of Regulatory Elements) combinée à un séquençage haut débit. Cette étude a démontré que suite au stimulus ovulatoire, les cellules de granulosa subissent une reprogrammation majeure des éléments de régulation, qui est corrélée avec une modification drastique de leurs fonctions biologiques. De plus, cette étude a mis en évidence une association majoritaire des éléments de régulation à des régions intergéniques distales et à des introns, indiquant que ces régions ont une importance capitale dans la régulation transcriptionnelle dans les cellules de granulosa. Cette étude a également permis d’identifier une panoplie de régulateurs transcriptionnels reconnus pour être essentiels à la fonction ovarienne, ainsi que leur sites de liaison dans le génome, démontrant que la méthode FAIRE est une méthode assez puissante pour permettre la prédiction d’événements moléculaires précis ayant un sens physiologique réel.

Une nouvelle approche computationnelle pour la découverte des sites de fixation de facteurs de transcription à l’ADN, adaptée aux données de ChIP-chip et de ChIP-séquençage

Les facteurs de transcription sont des protéines spécialisées qui jouent un rôle important dans différents processus biologiques tel que la différenciation, le cycle cellulaire et la tumorigenèse. Ils régulent la transcription des gènes en se fixant sur des séquences d’ADN spécifiques (éléments cis-régulateurs). L’identification de ces éléments est une étape cruciale dans la compréhension des réseaux de régulation des gènes. Avec l’avènement des technologies de séquençage à haut débit, l’identification de tout les éléments fonctionnels dans les génomes, incluant gènes et éléments cis-régulateurs a connu une avancée considérable. Alors qu’on est arrivé à estimer le nombre de gènes chez différentes espèces, l’information sur les éléments qui contrôlent et orchestrent la régulation de ces gènes est encore mal définie. Grace aux techniques de ChIP-chip et de ChIP-séquençage il est possible d’identifier toutes les régions du génome qui sont liées par un facteur de transcription d’intérêt. Plusieurs approches computationnelles ont été développées pour prédire les sites fixés par les facteurs de transcription. Ces approches sont classées en deux catégories principales: les algorithmes énumératifs et probabilistes. Toutefois, plusieurs études ont montré que ces approches génèrent des taux élevés de faux négatifs et de faux positifs ce qui rend difficile l’interprétation des résultats et par conséquent leur validation expérimentale. Dans cette thèse, nous avons ciblé deux objectifs. Le premier objectif a été de développer une nouvelle approche pour la découverte des sites de fixation des facteurs de transcription à l’ADN (SAMD-ChIP) adaptée aux données de ChIP-chip et de ChIP-séquençage. Notre approche implémente un algorithme hybride qui combine les deux stratégies énumérative et probabiliste, afin d’exploiter les performances de chacune d’entre elles. Notre approche a montré ses performances, comparée aux outils de découvertes de motifs existants sur des jeux de données simulées et des jeux de données de ChIP-chip et de ChIP-séquençage. SAMD-ChIP présente aussi l’avantage d’exploiter les propriétés de distributions des sites liés par les facteurs de transcription autour du centre des régions liées afin de limiter la prédiction aux motifs qui sont enrichis dans une fenêtre de longueur fixe autour du centre de ces régions. Les facteurs de transcription agissent rarement seuls. Ils forment souvent des complexes pour interagir avec l’ADN pour réguler leurs gènes cibles. Ces interactions impliquent des facteurs de transcription dont les sites de fixation à l’ADN sont localisés proches les uns des autres ou bien médier par des boucles de chromatine. Notre deuxième objectif a été d’exploiter la proximité spatiale des sites liés par les facteurs de transcription dans les régions de ChIP-chip et de ChIP-séquençage pour développer une approche pour la prédiction des motifs composites (motifs composés par deux sites et séparés par un espacement de taille fixe). Nous avons testé ce module pour prédire la co-localisation entre les deux demi-sites ERE qui forment le site ERE, lié par le récepteur des œstrogènes ERα. Ce module a été incorporé à notre outil de découverte de motifs SAMD-ChIP.

Análisis de la incidencia de la producción de biocombustibles a través de cultivos transgénicos en la seguridad alimentaria. Estudio de caso: Cultivos de soja transgénica en Argentina

La búsqueda de los Estados por mitigar su dependencia energética de las fuentes fósiles, ha traído consigo la búsqueda de energías alternativas, desencadenando en el uso y producción de biocombustibles. A su vez, la producción de estos últimos a través de cultivos transgénicos ha ido cobrando importancia en el escenario internacional. Esta opción se ha considerado como una salida al dilema de utilización de tierras "Biocombustible vs. Alimentos". En este contexto, el caso de Argentina, como uno de los mayores productores de cultivos transgénicos del mundo, entre los cuales se destaca la soja, se analiza en esta investigación por ser importante para determinar cuál es el impacto de los biocombustibles producidos a través de cultivos transgénicos en la seguridad alimentaria de la población.

Genomic architecture of adaptive color pattern divergence and convergence in Heliconius butterflies

Identifying the genetic changes driving adaptive variation in natural populations is key to understanding the origins of biodiversity. The mosaic of mimetic wing patterns in Heliconius butterflies makes an excellent system for exploring adaptive variation using next-generation sequencing. In this study, we use a combination of techniques to annotate the genomic interval modulating red color pattern variation, identify a narrow region responsible for adaptive divergence and convergence in Heliconius wing color patterns, and explore the evolutionary history of these adaptive alleles. We use whole genome resequencing from four hybrid zones between divergent color pattern races of Heliconius erato and two hybrid zones of the co-mimic Heliconius melpomene to examine genetic variation across 2.2 Mb of a partial reference sequence. In the intergenic region near optix, the gene previously shown to be responsible for the complex red pattern variation in Heliconius, population genetic analyses identify a shared 65-kb region of divergence that includes several sites perfectly associated with phenotype within each species. This region likely contains multiple cis-regulatory elements that control discrete expression domains of optix. The parallel signatures of genetic differentiation in H. erato and H. melpomene support a shared genetic architecture between the two distantly related co-mimics; however, phylogenetic analysis suggests mimetic patterns in each species evolved independently. Using a combination of next-generation sequencing analyses, we have refined our understanding of the genetic architecture of wing pattern variation in Heliconius and gained important insights into the evolution of novel adaptive phenotypes in natural populations.

Physiological changes in Campylobacter jejuni on entry into stationary phase

Campylobacter jejuni NCTC 11168 does not exhibit the general increase in cellular stress resistance on entry into stationary phase that is seen in most other bacteria. This is consistent with the lack of global stationary phase regulatory elements in this organism. deduced from an analysis of its genome sequence. We now show that C. jejuni NCTC 11168 does undergo certain changes in stationary phase, of a pattern not previously described. As cells entered stationary phase there was a change in membrane fatty acid composition, principally a decrease in the proportion of unsaturated fatty acids and an increase in the content of cyclopropane and short-chain fatty acids. These changes in membrane composition were accompanied by an increase in the resilience of the cell membrane towards loss of integrity caused by pressure and an increase in cellular pressure resistance. By contrast. there were no major changes in resistance to acid or heat treatment. A similar pattern of changes in stress resistance on entry, into stationary phase was seen in C. jejuni NCTC 11351, the type strain. These changes appear to represent a restricted Physiological response to the conditions existing in stationary phase cultures, in an organism having limited capacity for genetic regulation and adaptation to environment. &COPY; 2004 Elsevier B.V. All rights reserved.

Elevated levels of corticotrophin-releasing factor binding protein in the blood of patients suffering from arthritis and septicaemia and the presence of novel ligands in synovial fluid.

In view of the reported inflammatory effects of corticotrophin-releasing factor (CRF) and the associated regulatory elements in the gene of its binding protein (BP), we postulate that both BP as well as novel BP-ligands other than CRF may be involved in inflammatory disease. We have investigated BP in the blood of patients with arthritis and septicaemia and have attempted to identify CRF and other BP-ligands in synovial fluid. The BP was found to be significantly elevated in the blood of patients with rheumatoid arthritis and septicaemia. There was less BP-ligand and CRF in synovial fluid from patients with rheumatoid arthritis that from those with osteo- or psoriatic arthritis. There was at least 10-fold more BP-ligand than CRF in the fluid of all three groups of patients. A small amount of immunoreactive human (h)CRF, eluting in the expected position of CRF-41, was detected after high-pressure liquid chromatography of arthritic synovial fluid; however, the bulk of material with BP-ligand binding activity eluted earlier, suggesting that synovial fluid contained novel peptides that interacted with the BP. These results would suggest that the BP and its ligands could play an endocrine immunomodulatory role in inflammatory disease.

Changes in gene expression in Arabidopsis shoots during phosphate starvation and the potential for developing smart plants

Our aim was to generate and prove the concept of "smart" plants to monitor plant phosphorus (P) status in Arabidopsis. Smart plants can be genetically engineered by transformation with a construct containing the promoter of a gene up-regulated specifically by P starvation in an accessible tissue upstream of a marker gene such as beta-glucuronidase (GUS). First, using microarrays, we identified genes whose expression changed more than 2.5-fold in shoots of plants growing hydroponically when P, but not N or K, was withheld from the nutrient solution. The transient changes in gene expression occurring immediately (4 h) after P withdrawal were highly variable, and many nonspecific, shock-induced genes were up-regulated during this period. However, two common putative cis-regulatory elements (a PHO-like element and a TATA box-like element) were present significantly more often in the promoters of genes whose expression increased 4 h after the withdrawal of P compared with their general occurrence in the promoters of all genes represented on the microarray. Surprisingly, the expression of only four genes differed between shoots of P-starved and -replete plants 28 h after P was withdrawn. This lull in differential gene expression preceded the differential expression of a new group of 61 genes 100 h after withdrawing P. A literature survey indicated that the expression of many of these "late" genes responded specifically to P starvation. Shoots had reduced P after 100 h, but growth was unaffected. The expression of SQD1, a gene involved in the synthesis of sulfolipids, responded specifically to P starvation and was increased 100 h after withdrawing P. Leaves of Arabidopsis bearing a SQD1::GUS construct showed increased GUS activity after P withdrawal, which was detectable before P starvation limited growth. Hence, smart plants can monitor plant P status. Transferring this technology to crops would allow precision management of P fertilization, thereby maintaining yields while reducing costs, conserving natural resources, and preventing pollution.

Genetic responses to phosphorus deficiency

Background: Phosphorus (P) is an essential macronutrient for plants. Plants take up P as phosphate (Pi) from the soil solution. Since little Pi is available in most soils, P fertilizers are applied to crops. However, the use of P fertilizers is unsustainable and may cause pollution. Consequently, there is a need to develop more P-use-efficient (PUE) crops and precise methods to monitor crop P-status. Scope: Manipulating the expression of genes to improve the PUE of crops could reduce their P fertilizer requirement. This has stimulated research towards the identification of genes and signalling cascades involved in plant responses to P deficiency. Genes that respond to P deficiency can be grouped into 'early' genes that respond rapidly and often non-specifically to P deficiency, or 'late' genes that impact on the morphology, physiology or metabolism of plants upon Prolonged P deficiency. Summary: The use of micro-array technology has allowed researchers to catalogue the genetic responses of plants to P deficiency. Genes whose expression is altered by P deficiency include various transcription factors, which are thought to coordinate plant responses to P deficiency, and other genes involved in P acquisition and tissue P economy. Several common cis-regulatory elements have been identified in the promoters of these genes, suggesting that their expression might be coordinated. It is suggested that knowledge of the genes whose expression changes in response to P deficiency might allow the development of crops with improved PUE, and could be used in diagnostic techniques to monitor P deficiency in crops either directly using 'smart' indicator plants or indirectly through transcript profiling. The development of crops with improved PUE and the adoption of diagnostic technology could reduce production costs, minimize the use of a non-renewable resource, reduce pollution and enhance biodiversity.

An intronic G run within HIV-1 intron 2 is critical for splicing regulation of vif mRNA

Within target T lymphocytes, human immunodeficiency virus type I (HIV-1) encounters the retroviral restriction factor APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G; A3G), which is counteracted by the HIV-1 accessory protein Vif. Vif is encoded by intron-containing viral RNAs that are generated by splicing at 3' splice site (3'ss) A1 but lack splicing at 5'ss D2, which results in the retention of a large downstream intron. Hence, the extents of activation of 3'ss A1 and repression of D2, respectively, determine the levels of vif mRNA and thus the ability to evade A3G-mediated antiviral effects. The use of 3'ss A1 can be enhanced or repressed by splicing regulatory elements that control the recognition of downstream 5'ss D2. Here we show that an intronic G run (G(I2)-1) represses the use of a second 5'ss, termed D2b, that is embedded within intron 2 and, as determined by RNA deep-sequencing analysis, is normally inefficiently used. Mutations of G(I2)-1 and activation of D2b led to the generation of transcripts coding for Gp41 and Rev protein isoforms but primarily led to considerable upregulation of vif mRNA expression. We further demonstrate, however, that higher levels of Vif protein are actually detrimental to viral replication in A3G-expressing T cell lines but not in A3G-deficient cells. These observations suggest that an appropriate ratio of Vif-to-A3G protein levels is required for optimal virus replication and that part of Vif level regulation is effected by the novel G run identified here.

The Cochaperone BAG2 Sweeps Paired Helical Filament-Insoluble Tau from the Microtubule

Tau inclusions are a prominent feature of many neurodegenerative diseases including Alzheimer`s disease. Their accumulation in neurons as ubiquitinated filaments suggests a failure in the degradation limb of the Tau pathway. The components of a Tau protein triage system consisting of CHIP/Hsp70 and other chaperones have begun to emerge. However, the site of triage and the master regulatory elements are unknown. Here, we report an elegant mechanism of Tau degradation involving the cochaperone BAG2. The BAG2/Hsp70 complex is tethered to the microtubule and this complex can capture and deliver Tau to the proteasome for ubiquitin-independent degradation. This complex preferentially degrades Sarkosyl insoluble Tau and phosphorylated Tau. BAG2 levels in cells are under the physiological control of the microRNA miR-128a, which can tune paired helical filament Tau levels in neurons. Thus, we propose that ubiquitinated Tau inclusions arise due to shunting of Tau degradation toward a less efficient ubiquitin-dependent pathway.

Genetic and Functional Evidence Implicating DLL1 as the Gene That Influences Susceptibility to Visceral Leishmaniasis at Chromosome 6q27

Background. Visceral leishmaniasis (VL) is caused by Leishmania donovani and Leishmania infantum chagasi. Genome-wide linkage studies from Sudan and Brazil identified a putative susceptibility locus on chromosome 6q27. Methods. Twenty-two single-nucleotide polymorphisms (SNPs) at genes PHF10, C6orf70, DLL1, FAM120B, PSMB1, and TBP were genotyped in 193 VL cases from 85 Sudanese families, and 8 SNPs at genes PHF10, C6orf70, DLL1, PSMB1, and TBP were genotyped in 194 VL cases from 80 Brazilian families. Family-based association, haplotype, and linkage disequilibrium analyses were performed. Multispecies comparative sequence analysis was used to identify conserved noncoding sequences carrying putative regulatory elements. Quantitative reverse-transcription polymerase chain reaction measured expression of candidate genes in splenic aspirates from Indian patients with VL compared with that in the control spleen sample. Results. Positive associations were observed at PHF10, C6orf70, DLL1, PSMB1, and TBP in Sudan, but only at DLL1 in Brazil (combined P = 3 x 10(-4) at DLL1 across Sudan and Brazil). No functional coding region variants were observed in resequencing of 22 Sudanese VL cases. DLL1 expression was significantly (P = 2 x 10(-7)) reduced (mean fold change, 3.5 [SEM, 0.7]) in splenic aspirates from patients with VL, whereas other 6q27 genes showed higher levels (1.27 x 10(-6) < P < .01) than did the control spleen sample. A cluster of conserved noncoding sequences with putative regulatory variants was identified in the distal promoter of DLL1. Conclusions. DLL1, which encodes Delta-like 1, the ligand for Notch3, is strongly implicated as the chromosome 6q27 VL susceptibility gene.

Identificação de genes MUTM em BACs da cana-de-açucar e caracterização preliminar destes genes e seus promotores

Sugarcane has an importance in Brazil due to sugar and biofuel production. Considering this aspect, there is basic research being done in order to understand its physiology to improve production. The aim of this research is the Base Excision Repair pathway, in special the enzyme MUTM DNA-glycosylase (formamidopyrimidine) which recognizes oxidized guanine in DNA. The sugarcane scMUTM genes were analyzed using four BACs (Bacterial Artificial Chromosome) from a sugarcane genomic library from R570 cultivar. The resulted showed the presence in the region that had homology to scMUTM the presence of transposable elements. Comparing the similarity, it was observed a highest similarity to Sorghum bicolor sequence, both nucleotide and peptide sequences. Furthermore, promoter regions from MUTM genes in some grass showed different cis-regulatory elements, among which, most were related to oxidative stress, suggesting a gene regulation by oxidative stress