Studies of OPA1 pathogenic mechanisms in Dominant Optic Atrophy and novel protein function in mitochondrial DNA stability

The mitochondrion is an essential cytoplasmic organelle that provides most of the energy necessary for eukaryotic cell physiology. Mitochondrial structure and functions are maintained by proteins of both mitochondrial and nuclear origin. These organelles are organized in an extended network that dynamically fuses and divides. Mitochondrial morphology results from the equilibrium between fusion and fission processes, controlled by a family of “mitochondria-shaping” proteins. It is becoming clear that defects in mitochondrial dynamics can impair mitochondrial respiration, morphology and motility, leading to apoptotic cell death in vitro and more or less severe neurodegenerative disorders in vivo in humans. Mutations in OPA1, a nuclear encoded mitochondrial protein, cause autosomal Dominant Optic Atrophy (DOA), a heterogeneous blinding disease characterized by retinal ganglion cell degeneration leading to optic neuropathy (Delettre et al., 2000; Alexander et al., 2000). OPA1 is a mitochondrial dynamin-related guanosine triphosphatase (GTPase) protein involved in mitochondrial network dynamics, cytochrome c storage and apoptosis. This protein is anchored or associated on the inner mitochondrial membrane facing the intermembrane space. Eight OPA1 isoforms resulting from alternative splicing combinations of exon 4, 4b and 5b have been described (Delettre et al., 2001). These variants greatly vary among diverse organs and the presence of specific isoforms has been associated with various mitochondrial functions. The different spliced exons encode domains included in the amino-terminal region and contribute to determine OPA1 functions (Olichon et al., 2006). It has been shown that exon 4, that is conserved throughout evolution, confers functions to OPA1 involved in maintenance of the mitochondrial membrane potential and in the fusion of the network. Conversely, exon 4b and exon 5b, which are vertebrate specific, are involved in regulation of cytochrome c release from mitochondria, and activation of apoptosis, a process restricted to vertebrates (Olichon et al., 2007). While Mgm1p has been identified thanks to its role in mtDNA maintenance, it is only recently that OPA1 has been linked to mtDNA stability. Missense mutations in OPA1 cause accumulation of multiple deletions in skeletal muscle. The syndrome associated to these mutations (DOA-1 plus) is complex, consisting of a combination of dominant optic atrophy, progressive external ophtalmoplegia, peripheral neuropathy, ataxia and deafness (Amati- Bonneau et al., 2008; Hudson et al., 2008). OPA1 is the fifth gene associated with mtDNA “breakage syndrome” together with ANT1, PolG1-2 and TYMP (Spinazzola et al., 2009). In this thesis we show for the first time that specific OPA1 isoforms associated to exon 4b are important for mtDNA stability, by anchoring the nucleoids to the inner mitochondrial membrane. Our results clearly demonstrate that OPA1 isoforms including exon 4b are intimately associated to the maintenance of the mitochondrial genome, as their silencing leads to mtDNA depletion. The mechanism leading to mtDNA loss is associated with replication inhibition in cells where exon 4b containing isoforms were down-regulated. Furthermore silencing of exon 4b associated isoforms is responsible for alteration in mtDNA-nucleoids distribution in the mitochondrial network. In this study it was evidenced that OPA1 exon 4b isoform is cleaved to provide a 10kd peptide embedded in the inner membrane by a second transmembrane domain, that seems to be crucial for mitochondrial genome maintenance and does correspond to the second transmembrane domain of the yeasts orthologue encoded by MGM1 or Msp1, which is also mandatory for this process (Diot et al., 2009; Herlan et al., 2003). Furthermore in this thesis we show that the NT-OPA1-exon 4b peptide co-immuno-precipitates with mtDNA and specifically interacts with two major components of the mitochondrial nucleoids: the polymerase gamma and Tfam. Thus, from these experiments the conclusion is that NT-OPA1- exon 4b peptide contributes to the nucleoid anchoring in the inner mitochondrial membrane, a process that is required for the initiation of mtDNA replication and for the distribution of nucleoids along the network. These data provide new crucial insights in understanding the mechanism involved in maintenance of mtDNA integrity, because they clearly demonstrate that, besides genes implicated in mtDNA replications (i.e. polymerase gamma, Tfam, twinkle and genes involved in the nucleotide pool metabolism), OPA1 and mitochondrial membrane dynamics play also an important role. Noticeably, the effect on mtDNA is different depending on the specific OPA1 isoforms down-regulated, suggesting the involvement of two different combined mechanisms. Over two hundred OPA1 mutations, spread throughout the coding region of the gene, have been described to date, including substitutions, deletions or insertions. Some mutations are predicted to generate a truncated protein inducing haploinsufficiency, whereas the missense nucleotide substitutions result in aminoacidic changes which affect conserved positions of the OPA1 protein. So far, the functional consequences of OPA1 mutations in cells from DOA patients are poorly understood. Phosphorus MR spectroscopy in patients with the c.2708delTTAG deletion revealed a defect in oxidative phosphorylation in muscles (Lodi et al., 2004). An energetic impairment has been also show in fibroblasts with the severe OPA1 R445H mutation (Amati-Bonneau et al., 2005). It has been previously reported by our group that OPA1 mutations leading to haploinsufficiency are associated in fibroblasts to an oxidative phosphorylation dysfunction, mainly involving the respiratory complex I (Zanna et al., 2008). In this study we have evaluated the energetic efficiency of a panel of skin fibroblasts derived from DOA patients, five fibroblast cell lines with OPA1 mutations causing haploinsufficiency (DOA-H) and two cell lines bearing mis-sense aminoacidic substitutions (DOA-AA), and compared with control fibroblasts. Although both types of DOA fibroblasts maintained a similar ATP content when incubated in a glucose-free medium, i.e. when forced to utilize the oxidative phosphorylation only to produce ATP, the mitochondrial ATP synthesis through complex I, measured in digitonin-permeabilized cells, was significantly reduced in cells with OPA1 haploinsufficiency only, whereas it was similar to controls in cells with the missense substitutions. Furthermore, evaluation of the mitochondrial membrane potential (DYm) in the two fibroblast lines DOA-AA and in two DOA-H fibroblasts, namely those bearing the c.2819-2A>C mutation and the c.2708delTTAG microdeletion, revealed an anomalous depolarizing response to oligomycin in DOA-H cell lines only. This finding clearly supports the hypothesis that these mutations cause a significant alteration in the respiratory chain function, which can be unmasked only when the operation of the ATP synthase is prevented. Noticeably, oligomycin-induced depolarization in these cells was almost completely prevented by preincubation with cyclosporin A, a well known inhibitor of the permeability transition pore (PTP). This results is very important because it suggests for the first time that the voltage threshold for PTP opening is altered in DOA-H fibroblasts. Although this issue has not yet been addressed in the present study, several are the mechanisms that have been proposed to lead to PTP deregulation, including in particular increased reactive oxygen species production and alteration of Ca2+ homeostasis, whose role in DOA fibroblasts PTP opening is currently under investigation. Identification of the mechanisms leading to altered threshold for PTP regulation will help our understanding of the pathophysiology of DOA, but also provide a strategy for therapeutic intervention.

Assessing the mitochondrial DNA diversity of the Chagas disease vector Triatoma sordida (Hemiptera: Reduviidae)

Triatoma sordida is a species that transmits Trypanosoma cruzi to humans. In Brazil, T. sordida currently deserves special attention because of its wide distribution, tendency to invade domestic environments and vectorial competence. For the planning and execution of control protocols to be effective against Triatominae, they must consider its population structure. In this context, this study aimed to characterise the genetic variability of T. sordida populations collected in areas with persistent infestations from Minas Gerais, Brazil. Levels of genetic variation and population structure were determined in peridomestic T. sordida by sequencing a polymorphic region of the mitochondrial cytochrome b gene. Low nucleotide and haplotype diversity were observed for all 14 sampled areas; π values ranged from 0.002-0.006. Most obtained haplotypes occurred at low frequencies, and some were exclusive to only one of the studied populations. Interpopulation genetic diversity analysis revealed strong genetic structuring. Furthermore, the genetic variability of Brazilian populations is small compared to that of Argentinean and Bolivian specimens. The possible factors related to the reduced genetic variability and strong genetic structuring obtained for studied populations are discussed in this paper.

Geographic variation of sex pheromone and mitochondrial DNA in Diatraea saccharalis (Fab., 1794) (Lepidoptera: Crambidae)

(9Z,11E)-hexadecadienal and (Z11)-hexadecenal, the main sex pheromone components of the sugarcane borer, Diatraea saccharalis, were identified and quantified from four Brazilian and one Colombian populations using GC-EAD, GC-MS and GC analyses. Three different ratios were observed, 9:1,6:1, and 3:1. The pheromone concentration for the major component, (9Z,11E)-hexadecadienal, varied from 6.8 ng/gland to 21.9 ng/gland and from 1.7 ng/gland to 6.5 to the minor component, (Z11)-hexadecenal. The 25 D. saccharalis cytochrome oxidase II sequences that were analyzed showed low intra-specific variation and represented only 11 haplotypes, with the most frequent being the one represented by specimens from Sao Paulo, Parana, and Pernambuco states. Specimens from Colombia showed the highest genetic divergence from the others haplotypes studied. Data on the genetic variability among specimens, more than their geographic proximity, were in agreement with data obtained from analyses of the pheromone extracts. Our data demonstrate a variation in pheromone composition and a covariation in haplotypes of the D. saccharalis populations studied. (C) 2010 Elsevier Ltd. All rights reserved.

Evidence from mitochondrial DNA that head lice and body lice of humans (Phthiraptera : Pediculidae) are conspecific

The specific status of the head and body lice of humans has been debated for more than 200 yr. To clarify the specific status of head and body lice, we sequenced 524 base pairs (bp) of the cytochrome oxidase I (COI) gene of 28 head and 28 body lice from nine countries. Ten haplotypes that differed by 1-5 bp at II nucleotide positions were identified. A phylogeny of these sequences indicates that these head and body lice are not from reciprocally monophyletic lineages. Indeed, head and body lice share three of the 10 haplotypes we found. F-ST values and exact tests of haplotype frequencies showed significant differences between head and body lice. However, the same tests also showed significant differences among lice from different countries. Indeed, more of the variation in haplotype frequencies was explained by differences among lice from different countries than by differences between head and body lice. Our results indicate the following: (1) bead and body lice do not represent reciprocally monophyletic lineages and are conspecific; (2) gene flow among populations of lice from different countries is limited; and (3) frequencies of COI haplotypes can be used to study maternal gene flow among populations of head and body lice and thus transmission of lice among their human hosts.

Regulation of the DNA-binding and transcriptional activities of Xenopus laevis NFI-X by a novel C-terminal domain.

The nuclear factor I (NFI) family consists of sequence-specific DNA-binding proteins that activate both transcription and adenovirus DNA replication. We have characterized three new members of the NFI family that belong to the Xenopus laevis NFI-X subtype and differ in their C-termini. We show that these polypeptides can activate transcription in HeLa and Drosophila Schneider line 2 cells, using an activation domain that is subdivided into adjacent variable and subtype-specific domains each having independent activation properties in chimeric proteins. Together, these two domains constitute the full NFI-X transactivation potential. In addition, we find that the X. laevis NFI-X proteins are capable of activating adenovirus DNA replication through their conserved N-terminal DNA-binding domains. Surprisingly, their in vitro DNA-binding activities are specifically inhibited by a novel repressor domain contained within the C-terminal part, while the dimerization and replication functions per se are not affected. However, inhibition of DNA-binding activity in vitro is relieved within the cell, as transcriptional activation occurs irrespective of the presence of the repressor domain. Moreover, the region comprising the repressor domain participates in transactivation. Mechanisms that may allow the relief of DNA-binding inhibition in vivo and trigger transcriptional activation are discussed.

False phylogenies on wood mice due to cryptic cytochrome-b pseudogene.

The phylogeny and phylogeography of the Old World wood mice (subgenus Sylvaemus, genus Apodemus, Muridae) are well-documented. Nevertheless, the distributions of species, such as A. fulvipectus and A. ponticus remain dubious, as well as their phylogenetic relationships with A. sylvaticus. We analysed samples of Apodemus spp. across Europe using the mitochondrial cytochrome-b gene (cyt-b) and compared the DNA and amino-acid compositions of previously published sequences. The main result stemming from this study is the presence of a well-differentiated lineage of Sylvaemus including samples of various species (A. sylvaticus, A. fulvipectus, A. ponticus) from distant locations, which were revealed to be nuclear copies of the mitochondrial cyt-b. The presence of this cryptic pseudogene in published sequences is supported by different pathways. This has led to important errors in previous molecular trees and hence to partial misinterpretations in the phylogeny of Apodemus.

Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews.

The complete mitochondrial DNA (mtDNA) control region was amplified and directly sequenced in two species of shrew, Crocidura russula and Sorex araneus (Insectivora, Mammalia). The general organization is similar to that found in other mammals: a central conserved region surrounded by two more variable domains. However, we have found in shrews the simultaneous presence of arrays of tandem repeats in potential locations where repeats tend to occur separately in other mammalian species. These locations correspond to regions which are associated with a possible interruption of the replication processes, either at the end of the three-stranded D-loop structure or toward the end of the heavy-strand replication. In the left domain the repeated sequences (R1 repeats) are 78 bp long, whereas in the right domain the repeats are 12 bp long in C. russula and 14 bp long in S. araneus (R2 repeats). Variation in the copy number of these repeated sequences results in mtDNA control region length differences. Southern blot analysis indicates that level of heteroplasmy (more than one mtDNA form within an individual) differs between species. A comparative study of the R2 repeats in 12 additional species representing three shrew subfamilies provides useful indications for the understanding of the origin and the evolution of these homologous tandemly repeated sequences. An asymmetry in the distribution of variants within the arrays, as well as the constant occurrence of shorter repeated sequences flanking only one side of the R2 arrays, could be related to asymmetry in the replication of each strand of the mtDNA molecule. The pattern of sequence and length variation within and between species, together with the capability of the arrays to form stable secondary structures, suggests that the dominant mechanism involved in the evolution of these arrays in unidirectional replication slippage.

Effects of exercise on mitochondrial DNA content in skeletal muscle of patients with COPD

Background Exhausting exercise reduces the mitochondrial DNA (mtDNA) content in the skeletal muscle of healthy subjects due to oxidative damage. Since patients with chronic obstructive pulmonary disease (COPD) suffer enhanced oxidative stress during exercise, it was hypothesised that the mtDNA content will be further reduced. Objective To investigate the effects of exercise above and below the lactate threshold (LT) on the mtDNA content of skeletal muscle of patients with COPD. Methods Eleven patients with COPD (676 8 years; forced expiratory volume in 1s (FEV1)456 8%ref) and 10 healthy controls (666 4 years; FEV1 906 7% ref) cycled 45 min above LT (65% peak oxygen uptake (V9O2 peak)and another 7 patients (656 6 years; FEV1 506 4%ref)and 7 controls (566 9 years;FEV1 926 6%ref) cycled 45 min below their LT (50% V9O2 peak). Biopsies from the vastus lateralis muscle were obtained before exercise, immediately after and 1 h, 1 day and 1 week later to determine by PCR the mtDNA/nuclear DNA (nDNA) ratio (a marker of mtDNA content) and the expression of the peroxisome proliferator-activated receptor- g coactivator-1 a (PGC-1a)mRNA and the amount of reactive oxygen species produced during exercise was estimated from total V9O2. Results Skeletal muscle mtDNA/nDNA fell significantly after exercise above the LT both in controls and in patients with COPD, but the changes were greater in those with COPD. These changes correlated with production of reactive oxygen species, increases in manganese superoxide dismutase and PGC-1 a mRNA and returned to baseline values 1 week later. This pattern of response wa was also observed, albeit minimised, in patients exercising below the LT. Conclusions In patients with COPD, exercise enhances the decrease in mtDNA content of skeletal muscle and the expression of PGC-1 a mRNA seen in healthy subjects probably due to oxidative stress.

Geographic patterns of genetic variation in a broadly distributed marine vertebrate: new insights into loggerhead turtle stock structure from expanded mitochondrial DNA sequences

Previous genetic studies have demonstrated that natal homing shapes the stock structure of marine turtle nesting populations. However, widespread sharing of common haplotypes based on short segments of the mitochondrial control region often limits resolution of the demographic connectivity of populations. Recent studies employing longer control region sequences to resolve haplotype sharing have focused on regional assessments of genetic structure and phylogeography. Here we synthesize available control region sequences for loggerhead turtles from the Mediterranean Sea, Atlantic, and western Indian Ocean basins. These data represent six of the nine globally significant regional management units (RMUs) for the species and include novel sequence data from Brazil, Cape Verde, South Africa and Oman. Genetic tests of differentiation among 42 rookeries represented by short sequences (380 bp haplotypes from 3,486 samples) and 40 rookeries represented by long sequences (~800 bp haplotypes from 3,434 samples) supported the distinction of the six RMUs analyzed as well as recognition of at least 18 demographically independent management units (MUs) with respect to female natal homing. A total of 59 haplotypes were resolved. These haplotypes belonged to two highly divergent global lineages, with haplogroup I represented primarily by CC-A1, CC-A4, and CC-A11 variants and haplogroup II represented by CC-A2 and derived variants. Geographic distribution patterns of haplogroup II haplotypes and the nested position of CC-A11.6 from Oman among the Atlantic haplotypes invoke recent colonization of the Indian Ocean from the Atlantic for both global lineages. The haplotypes we confirmed for western Indian Ocean RMUs allow reinterpretation of previous mixed stock analysis and further suggest that contemporary migratory connectivity between the Indian and Atlantic Oceans occurs on a broader scale than previously hypothesized. This study represents a valuable model for conducting comprehensive international cooperative data management and research in marine ecology.

Restriction site heteroplasmy in the mitochondrial DNA of Brycon opalinus (Cuvier, 1819) (Characiformes, Characidae, Bryconiae)

Homoplasmy is a feature usually found in the mtDNA of higher animal taxa. On the other hand, the presence of two classes of mtDNA in the same cell or organism is rare and may appear in length or site variation. Data from mtDNA RFLP analysis of Brycon opalinus populations (Cuvier, 1819; Characiformes, Characidae, Bryconinae) revealed site heteroplasmy from endonuclease NheI digestion. Southern blotting hybridization was used to survey a total of 257 specimens with 24 restriction enzymes. Three different restriction fragment patterns of mtDNA were obtained from NheI digestion. Two individuals from hatchery broodstock were found to have two of them. NheI digests of heteroplasmic individuals yielded two fragments of approximately 1180 and 1260 bp. Despite the low frequency of this type of heteroplasmy in the whole B. opalinus population, the presence of site heteroplasmy in this species supports the evidence of this phenomenon in lower vertebrate groups.

Detection of transgenic DNA and protein, in rumen fluid, duodenal digesta, milk, blood and faeces of lactating dairy cows

The objective was to determine the presence or absence of transgenic and endogenous plant DNA in ruminal fluid, duodenal digesta, milk, blood, and feces, and if found, to determine fragment size. Six multiparous lactating Holstein cows fitted with ruminal and duodenal cannulas received a total mixed ration. There were two treatments (T). In T1, the concentrate contained genetically modified (GM) soybean meal (cp4epsps gene) and GM corn grain (cry1a[b] gene), whereas T2 contained the near isogenic non-GM counterparts. Polymerase chain reaction analysis was used to determine the presence or absence of DNA sequences. Primers were selected to amplify small fragments from single-copy genes (soy lectin and corn high-mobility protein and cp4epsps and cry1a[b] genes from the GM crops) and multicopy genes (bovine mitochondrial cytochrome b and rubisco). Single-copy genes were only detected in the solid phase of rumen and duodenal digesta. In contrast, fragments of the rubisco gene were detected in the majority of samples analyzed in both the liquid and solid phases of ruminal and duodenal digesta, milk, and feces, but rarely in blood. The size of the rubisco gene fragments detected decreased from 1176 bp in ruminal and duodenal digesta to 351 bp in fecal samples.

Wolbachia infection and dramatic intraspecific mitochondrial dna divergence in a fig wasp

Mitochondria and Wolbachia are maternally inherited genomes that exhibit strong linkage disequilibrium in many organisms. We surveyed Wolbachia infections in 187 specimens of the fig wasp species, Ceratosolen solmsi, and found an infection prevalence of 89.3%. DNA sequencing of 20 individuals each from Wolbachia-infected and -uninfected subpopulations revealed extreme mtDNA divergence (up to 9.2% and 15.3% in CO1 and cytochrome b, respectively) between infected and uninfected wasps. Further, mtDNA diversity was significantly reduced within the infected group. Our sequencing of a large part of the mitochondrial genome from both Wolbachia-infected and -uninfected individuals revealed that high sequence divergence is common throughout the mitochondrial genome. These patterns suggest a partial selective sweep of mitochondria subsequent to the introduction of Wolbachia into C. solsmi, by hybrid introgression from a related species.

Clonal expansion of early to mid-life mitochondrial DNA point mutations drives mitochondrial dysfunction during human ageing

Age-related decline in the integrity of mitochondria is an important contributor to the human ageing process. In a number of ageing stem cell populations, this decline in mitochondrial function is due to clonal expansion of individual mitochondrial DNA (mtDNA) point mutations within single cells. However the dynamics of this process and when these mtDNA mutations occur initially are poorly understood. Using human colorectal epithelium as an exemplar tissue with a well-defined stem cell population, we analysed samples from 207 healthy participants aged 17-78 years using a combination of techniques (Random Mutation Capture, Next Generation Sequencing and mitochondrial enzyme histochemistry), and show that: 1) non-pathogenic mtDNA mutations are present from early embryogenesis or may be transmitted through the germline, whereas pathogenic mtDNA mutations are detected in the somatic cells, providing evidence for purifying selection in humans, 2) pathogenic mtDNA mutations are present from early adulthood (<20 years of age), at both low levels and as clonal expansions, 3) low level mtDNA mutation frequency does not change significantly with age, suggesting that mtDNA mutation rate does not increase significantly with age, and 4) clonally expanded mtDNA mutations increase dramatically with age. These data confirm that clonal expansion of mtDNA mutations, some of which are generated very early in life, is the major driving force behind the mitochondrial dysfunction associated with ageing of the human colorectal epithelium.

Taxonomic and Phylogenetic Relationships Between Species of the Genus Culex (Diptera: Culicidae) From Brazil Inferred From the Cytochrome c Oxidase I Mitochondrial Gene

Species of the genus Culex Linnaeus have been incriminated as the main vectors of lymphatic filariases and are important vectors of arboviruses, including West Nile virus. Sequences corresponding to a fragment of 478 bp of the cytochrome c oxidase subunit I gene, which includes part of the barcode region, of 37 individuals of 17 species of genus Culex were generated to establish relationships among five subgenera, Culex, Phenacomyia, Melanoconion, Microculex, and Carrollia, and one species of the genus Lutzia that occurs in Brazil. Bayesian methods were employed for the phylogenetic analyses. Results of sequence comparisons showed that individuals identified as Culex dolosus, Culex mollis, and Culex imitator possess high intraspecific divergence (3.1, 2.3, and 3.5%, respectively) when using the Kimura two parameters model. These differences were associated either with distinct morphological characteristics of the male genitalia or larval and pupal stages, suggesting that these may represent species complexes. The Bayesian topology suggested that the genus and subgenus Culex are paraphyletic relative to Lutzia and Phenacomyia, respectively. The cytochrome c oxidase subunit I sequences may be a useful tool to both estimate phylogenetic relationships and identify morphologically similar species of the genus Culex.

Genetic structure, phylogeny, and biogeography of Brazilian eyelid-less lizards of genera Calyptommatus and Nothobachia (Squamata, Gymnophthalmidae) as inferred from mitochondrial DNA sequences

Calyptommatus and Nothobachia genera of gymnophthalmid lizards are restricted to sandy open habitats on Sao Francisco River margins, northeastern Brazil. Phylogenetic relationships and geographic distribution of the four recognized species of Calyptommatus were analyzed from partial mitochondrial cyt b, 12S, and 16S rRNA genes sequencing, taking allopatric populations of the monotypic Nothobachia ablephara as the outgroup. In Calyptommatus a basal split separated C. sinebrachiatus, a species restricted to the eastern bank of the river, from the three other species. In this clade, C. confusionibus, found on western margin, was recovered as the sister group of the two other species, C. leiolepis and C. nicterus, from opposite margins. According to approximate date estimations, C. sinebrachiatus would have separated from the other congeneric species by 4.4-6.5 my, and C. nicterus, also from eastern bank, would be diverging by 1.8-2.6 my from C. leiolepis, the sister species on the opposite margin. C. confusionibus and C. leiolepis, both from western sandy areas, would be differentiating by 2.8-5.0 my. Divergence times of about 3.0-4.0 my were estimated for allopatric populations of Nothobachia restricted to western margin. Significant differences in 16S rRNA secondary structure relatively to other vertebrates are reported. Distinct evolutionary patterns are proposed for different taxa in those sandy areas, probably related to historical changes in the course of Sao Francisco River. (C) 2010 Elsevier Inc. All rights reserved.