Suppression subtractive hybridization profiles of radial growth phase and metastatic melanoma cell lines reveal novel potential targets

Background: Melanoma progression occurs through three major stages: radial growth phase (RGP), confined to the epidermis; vertical growth phase (VGP), when the tumor has invaded into the dermis; and metastasis. In this work, we used suppression subtractive hybridization (SSH) to investigate the molecular signature of melanoma progression, by comparing a group of metastatic cell lines with an RGP-like cell line showing characteristics of early neoplastic lesions including expression of the metastasis suppressor KISS1, lack of alpha v beta 3-integrin and low levels of RHOC. Methods: Two subtracted cDNA collections were obtained, one (RGP library) by subtracting the RGP cell line (WM1552C) cDNA from a cDNA pool from four metastatic cell lines (WM9, WM852, 1205Lu and WM1617), and the other (Met library) by the reverse subtraction. Clones were sequenced and annotated, and expression validation was done by Northern blot and RT-PCR. Gene Ontology annotation and searches in large-scale melanoma expression studies were done for the genes identified. Results: We identified 367 clones from the RGP library and 386 from the Met library, of which 351 and 368, respectively, match human mRNA sequences, representing 288 and 217 annotated genes. We confirmed the differential expression of all genes selected for validation. In the Met library, we found an enrichment of genes in the growth factors/receptor, adhesion and motility categories whereas in the RGP library, enriched categories were nucleotide biosynthesis, DNA packing/repair, and macromolecular/vesicular trafficking. Interestingly, 19% of the genes from the RGP library map to chromosome 1 against 4% of the ones from Met library. Conclusion: This study identifies two populations of genes differentially expressed between melanoma cell lines from two tumor stages and suggests that these sets of genes represent profiles of less aggressive versus metastatic melanomas. A search for expression profiles of melanoma in available expression study databases allowed us to point to a great potential of involvement in tumor progression for several of the genes identified here. A few sequences obtained here may also contribute to extend annotated mRNAs or to the identification of novel transcripts.

Functional characterization of the sciarid BhC4-1 core promoter in transgenic Drosophila

Background: Core promoters are cis-regulatory modules to which bind the basal transcriptional machinery and which participate in the regulation of transcription initiation. Although core promoters have not been extensively investigated through functional assays in a chromosomal context, the available data suggested that the response of a given core promoter might vary depending on the promoter context. Previous studies suggest that a (-57/+40) fragment constitutes the core promoter of the BhC4-1 gene which is located in DNA puff C4 of the sciarid fly Bradysia hygida. Here we tested this (-57/+40) fragment in distinct regulatory contexts in order to verify if promoter context affects its core promoter activity. Results: Consistent with the activity of a core promoter, we showed that in the absence of upstream regulatory sequences the (-57/+40) fragment drives low levels of reporter gene mRNA expression throughout development in transgenic Drosophila. By assaying the (-57/+40) fragment in two distinct regulatory contexts, either downstream of the previously characterized Fbp1 enhancer or downstream of the UAS element, we showed that the BhC4-1 core promoter drives regulated transcription in both the germline and in various tissues throughout development. Furthermore, the use of the BhC4-1 core promoter in a UAS construct significantly reduced salivary gland ectopic expression in third instar larvae, which was previously described to occur in the context of the GAL4/UAS system. Conclusions: Our results from functional analysis in transgenic Drosophila show that the BhC4-1 core promoter drives gene expression regardless of the promoter context that was assayed. New insights into the functioning of the GAL4/UAS system in Drosophila were obtained, indicating that the presence of the SV40 sequence in the 3' UTR of a UAS construct does not preclude expression in the germline. Furthermore, our analysis indicated that ectopic salivary gland expression in the GAL4/UAS system does not depend only on sequences present in the GAL4 construct, but can also be affected by the core promoter sequences in the UAS construct. In this context, we propose that the sciarid BhC4-1 core promoter constitutes a valuable core promoter which can be employed in functional assays in insects.

Genome Survey Sequence Analysis and Identification of Homologs of Major Surface Protease (gp63) Genes in Trypanosoma rangeli

In this study, 222 genome survey sequences were generated for Trypanosoma rangeli strain P07 isolated from an opossum (Didelphis albiventris) in Minas Gerais State, Brazil. T. rangeli sequences were compared by BLASTX (Basic Local Alignment Search Tool X) analysis with the assembled contigs of Leishmania braziliensis, Leishmania infantum, Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi. Results revealed that 82% (182/222) of the sequences were associated with predicted proteins described, whereas 18% (40/222) of the sequences did not show significant identity with sequences deposited in databases, suggesting that they may represent T. rangeli-specific sequences. Among the 182 predicted sequences, 179 (80.6%) had the highest similarity with T. cruzi, 2 (0.9%) with T. brucei, and 1 (0.5%) with L. braziliensis. Computer analysis permitted the identification of members of various gene families described for trypanosomatids in the genome of T. rangeli, such as trans-sialidases, mucin-associated surface proteins, and major surface proteases (MSP or gp63). This is the first report identifying sequences of the MSP family in T. rangeli. Multiple sequence alignments showed that the predicted MSP of T. rangeli presented the typical characteristics of metalloproteases, such as the presence of the HEXXH motif, which corresponds to a region previously associated with the catalytic site of the enzyme, and various cysteine and proline residues, which are conserved among MSPs of different trypanosomatid species. Reverse transcriptase-polymerase chain reaction analysis revealed the presence of MSP transcripts in epimastigote forms of T. rangeli.

Suppression subtractive hybridization analysis reveals expression of conserved and novel genes in male accessory glands of the ant Leptothorax gredleri

Background: During mating, insect males eject accessory gland proteins (Acps) into the female genital tract. These substances are known to affect female post-mating behavior and physiology. In addition, they may harm the female, e. g., in reducing its lifespan. This is interpreted as a consequence of sexual antagonistic co-evolution. Whereas sexual conflict abounds in non-social species, the peculiar life history of social insects (ants, bees, wasps) with lifelong pair-bonding and no re-mating aligns the reproductive interests of the sexes. Harming the female during mating would negatively affect male fitness and sexual antagonism is therefore not expected. Indeed, mating appears to increase female longevity in at least one ant species. Acps are presumed to play a role in this phenomenon, but the underlying mechanisms are unknown. In this study, we investigated genes, which are preferentially expressed in male accessory glands of the ant Leptothorax gredleri, to determine which proteins might be transferred in the seminal fluid. Results: By a suppression subtractive hybridization protocol we obtained 20 unique sequences (USs). Twelve had mutual best matches with genes predicted for Apis mellifera and Nasonia vitripennis. Functional information (Gene Ontology) was available only for seven of these, including intracellular signaling, energy-dependent transport and metabolic enzyme activities. The remaining eight USs did not match sequences from other species. Six genes were further analyzed by quantitative RT-PCR in three life cycle stages of male ants. A gene with carboxy-lyase activity and one of unpredicted function were significantly overexpressed in accessory glands of sexually mature males. Conclusions: Our study is the first one to investigate differential gene expression in ants in a context related to mating. Our findings indicate that male accessory glands of L. gredleri express a series of genes that are unique to this species, possibly representing novel genes, in addition to conserved ones for which functions can be predicted. Identifying differentially expressed genes might help to better understand molecular mechanisms involved in reproductive processes in eusocial Hymenoptera. While the novel genes could account for rapidly evolving ones driven by intra-sexual conflict between males, conserved genes imply that rather beneficial traits might get fixed by a process described as inter-sexual cooperation between males and females.

Novel Primate-Specific Genes, RMEL 1, 2 and 3, with Highly Restricted Expression in Melanoma, Assessed by New Data Mining Tool

Melanoma is a highly aggressive and therapy resistant tumor for which the identification of specific markers and therapeutic targets is highly desirable. We describe here the development and use of a bioinformatic pipeline tool, made publicly available under the name of EST2TSE, for the in silico detection of candidate genes with tissue-specific expression. Using this tool we mined the human EST (Expressed Sequence Tag) database for sequences derived exclusively from melanoma. We found 29 UniGene clusters of multiple ESTs with the potential to predict novel genes with melanoma-specific expression. Using a diverse panel of human tissues and cell lines, we validated the expression of a subset of three previously uncharacterized genes (clusters Hs.295012, Hs.518391, and Hs.559350) to be highly restricted to melanoma/melanocytes and named them RMEL1, 2 and 3, respectively. Expression analysis in nevi, primary melanomas, and metastatic melanomas revealed RMEL1 as a novel melanocytic lineage-specific gene up-regulated during melanoma development. RMEL2 expression was restricted to melanoma tissues and glioblastoma. RMEL3 showed strong up-regulation in nevi and was lost in metastatic tumors. Interestingly, we found correlations of RMEL2 and RMEL3 expression with improved patient outcome, suggesting tumor and/or metastasis suppressor functions for these genes. The three genes are composed of multiple exons and map to 2q12.2, 1q25.3, and 5q11.2, respectively. They are well conserved throughout primates, but not other genomes, and were predicted as having no coding potential, although primate-conserved and human-specific short ORFs could be found. Hairpin RNA secondary structures were also predicted. Concluding, this work offers new melanoma-specific genes for future validation as prognostic markers or as targets for the development of therapeutic strategies to treat melanoma.

Hantavirus Pulmonary Syndrome, Central Plateau, Southeastern, and Southern Brazil

Hantavirus pulmonary syndrome (HPS) is an increasing health problem in Brazil because of encroachment of sprawling urban, agricultural, and cattle-raising areas into habitats of subfamily Sigmodontinae rodents, which serve as hantavirus reservoirs. From 1993 through June 2007, a total of 884 cases of HIPS were reported in Brazil (case-fatality rate 39%). To better understand this emerging disease, we collected 89 human serum samples and 68 rodent lung samples containing antibodies to hantavirus from a 2,500-km-wide area in Brazil. RNA was isolated from human samples and rodent lung tissues and subjected to reverse transcription-PCR. Partial sequences of nucleocapsid protein and glycoprotein genes from 22 human and 16 rodent sources indicated only Araraquara virus and Juquitiba virus lineages. The case-fatality rate of HPS was higher in the area with Araraquara virus. This virus, which may be the most virulent hantavirus in Brazil, was associated with areas that have had greater anthropogenic changes.

Genotypic Characteristics of HIV Type 1 Based on gp120 Hypervariable Region 3 of Isolates from Southern Brazil

The aim of this study was to investigate HIV-1 molecular diversity and the epidemiological profile of HIV-1-infected patients from Ribeirao Preto, Brazil. A nested PCR followed by sequencing of a 302-base pair fragment of the env gene (C2-V3 region) was performed in samples from HIV-1-positive patients. A total of 45 sequences were aligned with final manual adjustments. The phylogenetic analyses showed a higher prevalence of HIV-1 subtype B in the studied population (97.8%) with only one sample yielding an F1 subtype. The viral genotyping prediction showed that CCR5 tropism was the most prevalent in the studied cohort. Geno2pheno analysis showed that R5 and CXCR4 prediction were 69% and 31%, respectively. There was no statistical significance, either in viral load or in CD4(+) T cell count when R5 and X4 prediction groups were compared. Moreover, the GPGR tetramer was the most common V3 loop core motif identified in the HIV-1 strains studied (34.1%) followed by GWGR, identified in 18.1% of the samples. The high level of B subtype in this Brazilian population reinforces the nature of the HIV epidemic in Brazil, and corroborates previous data obtained in the Brazilian HIV-infected population.

Chagas Cardiomyopathy Manifestations and Trypanosoma cruzi Genotypes Circulating in Chronic Chagasic Patients

Chagas disease caused by Trypanosoma cruzi is a complex disease that is endemic and an important problem in public health in Latin America. The T. cruzi parasite is classified into six discrete taxonomic units (DTUs) based on the recently proposed nomenclature (TcI, TcII, TcIII, TcIV, TcV and TcVI). The discovery of genetic variability within TcI showed the presence of five genotypes (Ia, Ib, Ic, Id and Ie) related to the transmission cycle of Chagas disease. In Colombia, TcI is more prevalent but TcII has also been reported, as has mixed infection by both TcI and TcII in the same Chagasic patient. The objectives of this study were to determine the T. cruzi DTUs that are circulating in Colombian chronic Chagasic patients and to obtain more information about the molecular epidemiology of Chagas disease in Colombia. We also assessed the presence of electrocardiographic, radiologic and echocardiographic abnormalities with the purpose of correlating T. cruzi genetic variability and cardiac disease. Molecular characterization was performed in Colombian adult chronic Chagasic patients based on the intergenic region of the mini-exon gene, the 24S alpha and 18S regions of rDNA and the variable region of satellite DNA, whereby the presence of T. cruzi I, II, III and IV was detected. In our population, mixed infections also occurred, with TcI-TcII, TcI-TcIII and TcI-TcIV, as well as the existence of the TcI genotypes showing the presence of genotypes Ia and Id. Patients infected with TcI demonstrated a higher prevalence of cardiac alterations than those infected with TcII. These results corroborate the predominance of TcI in Colombia and show the first report of TcIII and TcIV in Colombian Chagasic patients. Findings also indicate that Chagas cardiomyopathy manifestations are more correlated with TcI than with TcII in Colombia.

Report of a chimeric origin of transposable elements in a bovine-coding gene

Despite the wide distribution of transposable elements (TEs) in mammalian genomes, part of their evolutionary significance remains to be discovered. Today there is a substantial amount of evidence showing that TEs are involved in the generation of new exons in different species. In the present study, we searched 22,805 genes and reported the occurrence of TE-cassettes in coding sequences of 542 cow genes using the RepeatMasker program. Despite the significant number (542) of genes with TE insertions in exons only 14 (2.6%) of them were translated into protein, which we characterized as chimeric genes. From these chimeric genes, only the FAST kinase domains 3 (FASTKD3) gene, present on chromosome BTA 20, is a functional gene and showed evidence of the exaptation event. The genome sequence analysis showed that the last exon coding sequence of bovine FASTKD3 is similar to 85% similar to the ART2A retrotransposon sequence. In addition, comparison among FASTKD3 proteins shows that the last exon is very divergent from those of Homo sapiens, Pan troglodytes and Canis familiares. We suggest that the gene structure of bovine FASTKD3 gene could have originated by several ectopic recombinations between TE copies. Additionally, the absence of TE sequences in all other species analyzed suggests that the TE insertion is clade-specific, mainly in the ruminant lineage.

Tick-Borne Bacteria in Free-Living Jaguars (Panthera onca) in Pantanal, Brazil

Tick-borne bacteria were investigated in 10 free-living jaguars and their ticks in the Pantanal biome, Brazil. Jaguar sera were tested by indirect fluorescent antibody assays using Rickettsia rickettsii, Rickettsia parkeri, Rickettsia amblyommii, Rickettsia rhipicephali, Rickettsia felis, Rickettsia bellii, Ehrlichia canis, and Coxiella burnetii as crude antigens. All 10 jaguar sera reacted (titer >= 64) to at least one Rickettsia species; 4 and 3 sera reacted with E. canis and C. burnetii, respectively. One jaguar presented antibody titer to R. parkeri at least fourfold higher than those to any of the other five Rickettsia antigens, suggesting that this animal was infected by R. parkeri. Ticks collected from jaguars included the species Amblyomma cajennense, Amblyomma triste, and Rhipicephalus (Boophilus) microplus. No Rickettsia DNA was detected in jaguar blood samples, but an A. triste specimen collected on a jaguar was shown by PCR to be infected by R. parkeri. The blood of two jaguars and samples of A. triste, A. cajennense, and Amblyomma sp. yielded Ehrlichia DNA by PCR targeting the ehrlichial genes 16S rRNA and dsb. Partial DNA sequences obtained from PCR products resulted in a new ehrlichial strain, here designated as Ehrlichia sp. strain Jaguar. A partial DNA sequence of the 16S rRNA gene of this novel strain showed to be closest (99.0%) to uncultured strains of Ehrlichia sp. from Japan and Russia and 98.7% identical to different strains of Ehrlichia ruminantium. The ehrlichial dsb partial sequence of strain jaguar showed to be at most 80.7% identical to any Ehrlichia species or genotype available in GenBank. Through phylogenetic analysis, Ehrlichia sp. strain jaguar grouped in a cluster, albeit distantly, with different genotypes of E. ruminantium. Results highlight risks for human and animal health, considering that cattle ranching and ecotourism are major economic activities in the Pantanal region of Brazil.

New Epidemiological Data on Brazilian Spotted Fever in an Endemic Area of the State of Sao Paulo, Brazil

The present work evaluated rickettsial infection in dogs and their ticks in an area endemic for Brazilian spotted fever (BSF) in the metropolitan area of Sao Paulo, Brazil, where the tick Amblyomma aureolatum was presumed to be the vector of the disease. Ticks were collected on dogs from 185 houses, encompassing single infestations by Rhipicephalus sanguineus, Amblyomma aureolatum, Amblyomma longirostre, or Amblyomma sp. in dogs from 60 (32.4%), 77 (41.6%), 2 (1.1%), and 25 (13.5%) houses, respectively; 19 (10.3%) houses had dogs with mixed infestations by R. sanguineus and A. aureolatum; 1 (0.5%) house had dogs with infestations by A. aureolatum and A. longirostre; and 1 (0.5%) house had dogs with infestations by R. sanguineus and Amblyomma sp. Overall, A. aureolatum was present in dogs from 97 (52.4%) houses, and R. sanguineus in dogs from 80 (43.2%) houses. A total of 287 ticks (130 A. aureolatum and 157 R. sanguineus) infesting dogs from 98 houses were selected for testing by polymerase chain reaction (PCR) targeting rickettsial genes. Overall, 3.1% of the A. aureolatum ticks were infected by Rickettsia bellii, and 1.3% of the R. sanguineus were infected by Ricketttsii rickettsii. For serology, we selected 23 dogs living in and in the vicinity of the house where the R. rickettsii-infected ticks were collected. The indirect fluorescent antibody (IFA) test detected antibodies reactive with R. rickettsii in sera from 16 (69.6%) dogs, with titers ranging from 256 to 32,768. It is established that Amblyomma aureolatum is a vector of R. rickettsii in the Sao Paulo metropolitan area, but our results highlight for the first time in Brazil, a possible role of R. sanguineus in the epidemiology of R. rickettsii, corroborating previous findings in Mexico and the United States, where R. sanguineus has been implicated in the transmission of R. rickettsii to humans.

ISOLATION OF Rickettsia bellii FROM Amblyomma ovale AND Amblyomma incisum TICKS FROM SOUTHERN BRAZIL

Objective. To isolate and characterize rickettsiae from the ticks Amblyomma ovale and Amblyomma incisum collected in the state of Sao Paulo. Materials and methods. Adult, free-living A. ovale and A. incisum were collected in an Atlantic rainforest area in the state of Sao Paulo, Brazil. Each tick was tested using the hemolymph assay; samples from positive ticks were placed in shell vials in order to isolate rickettsiae and subsequently grown in Vero cells. Amplification of three rickettsial genes ( gltA, htrA and ompA) was attempted using polymerase chain reaction (PCR) for each isolate obtained. Amplicons were subsequently sequenced. Results. A total of 388 A. incisum and 50 A. ovale were collected. Only one A. incisum and one A. ovale were hemolymph-test positive. Rickettsiae were successfully isolated from these ticks; however establishment in Vero cell culture was successful only for the isolate from A. ovale. Bacterial contamination in the first cell passage of the A. incisum isolate precluded successful isolation of the organism. PCR products were obtained with the gltA and htrA primers for the two isolates, but no product was obtained with the ompA primers. By BLAST analysis, partial gltA and htrA sequences of isolates from A. ovale and A. incisum were similar to the corresponding sequences of R. bellii. Conclusions. This is the first report of R. bellii infecting A. incisum and the first successful isolation from A. ovale.

Determination and Molecular Analysis of the Complete Genome Sequence of Two Wild-Type Rabies Viruses Isolated from a Haematophagous Bat and a Frugivorous Bat in Brazil

The complete genome sequences of two Brazilian wild-type rabies viruses (RABV), a BR-DR1 isolate from a haematophagous bat (Desmodus rotundus) and a BR-AL1 isolate from a frugivorous bat (Artibeus lituratus), were determined. The genomes of the BR-DR1 and RR-AL1 had 11,923 and 11,922 nt, respectively, and both encoded the five standard genes of rhabdoviruses. The complete nucleotide sequence identity between the BR-DR1 and BR-AL1 isolates was 97%. The BR-DR1 and BR-AL1 isolates had some conserved functional sites revealed by the fixed isolates, whereas both isolates had unique amino acid substitutions in the antigenic region IV of the nucleocapsid gene. Therefore, it is speculated that both isolates were nearly identical in virologic character. According to our phylogenetic analysis based on the complete genomes, both isolates belonged to genotype 1, and to the previously defined ""vampire bat-related RABV lineage"" which consisted of mainly D. rotundus- and A. lituratus- isolates; however, a branch pattern with high bootstrap values suggested that BR-DR1 was more closely related to the 9001FRA isolate, which was collected from a dog bitten by a bat in French Guiana, than to BR-AL1. This result suggests that the vampire bat-related RABV lineage includes Brazilian vampire bat and Brazilian frugivorous bat RABV and is further divided into Brazilian vampire bat and Brazilian frugivorous bat RABV sub-lineages. The phylogenetic analysis based on the complete genomes was valuable in discriminating among very closely related isolates.

DESCRIPTION OF ADULTS AND NYMPH, AND REDESCRIPTION OF THE LARVA, OF ORNITHODOROS MARINKELLEI (ACARI: ARGASIDAE), WITH DATA ON ITS PHYLOGENETIC POSITION

The argasid tick Ornithodoros marinkellei Kohls, Clifford, and Jones, 1969 was described 4 decades ago based on larval specimens collected from bats (Pteronotus spp.) in Colombia and Panama. Thereafter, larval O. marinkellei parasitizing bats were reported from Venezuela, Guyana, and Brazil. Herein, we describe the adults and nymph, and redescribe the larva of O. marinkellei based on specimens recently collected in the western Brazilian Amazon region. In contrast to all other known adult argasids, the idiosoma of both males and females of O. marinkellei is covered with sclerotized plaques. The idiosoma of the nymph of O. marinkellei is entirely micromamillated, and differs from the adults by the absence of plaques. The larva of O. marinkellei is morphologically similar to the larvae of the 2 other species belonging to the subgenus Subparmatus, i.e., Ornithodoros viguerasi Cooley and Kohls, 1941 and Ornithodoros mormoops Kohls, Clifford, and Jones, 1969. Because of the long and narrow dorsal plate, the larva of O. marinkellei is readily distinguished from O. viguerasi and O. mormoops. Comparison of our larvae from Brazil with O. marinkellei paratype specimens from Colombia confirmed their taxonomic identification. However, a few morphological differences, particularly in the size of the gnathosoma, were observed. Further studies are necessary to clarify whether O. marinkellei is a complex of different species, or a single species represented by morphologically polymorphic, and geographically distinct populations. Partial mitochondrial 16S rDNA gene sequences were generated for O. marinkellei specimens from Brazil, and compared with available homologous sequences in GenBank. Phylogenetic analyses revealed O. marinkellei to be distinct from the remaining argasid species available in GenBank, including other bat-associated tick species that are found in sympatry with O. marinkellei in the Neotropical region.

Lineage divergence detected in the malaria vector Anopheles marajoara (Diptera: Culicidae) in Amazonian Brazil

Background: Cryptic species complexes are common among anophelines. Previous phylogenetic analysis based on the complete mtDNA COI gene sequences detected paraphyly in the Neotropical malaria vector Anopheles marajoara. The ""Folmer region"" detects a single taxon using a 3% divergence threshold. Methods: To test the paraphyletic hypothesis and examine the utility of the Folmer region, genealogical trees based on a concatenated (white + 3' COI sequences) dataset and pairwise differentiation of COI fragments were examined. The population structure and demographic history were based on partial COI sequences for 294 individuals from 14 localities in Amazonian Brazil. 109 individuals from 12 localities were sequenced for the nDNA white gene, and 57 individuals from 11 localities were sequenced for the ribosomal DNA (rDNA) internal transcribed spacer 2 (ITS2). Results: Distinct A. marajoara lineages were detected by combined genealogical analysis and were also supported among COI haplotypes using a median joining network and AMOVA, with time since divergence during the Pleistocene (< 100,000 ya). COI sequences at the 3' end were more variable, demonstrating significant pairwise differentiation (3.82%) compared to the more moderate 2.92% detected by the Folmer region. Lineage 1 was present in all localities, whereas lineage 2 was restricted mainly to the west. Mismatch distributions for both lineages were bimodal, likely due to multiple colonization events and spatial expansion (similar to 798 - 81,045 ya). There appears to be gene flow within, not between lineages, and a partial barrier was detected near Rio Jari in Amapa state, separating western and eastern populations. In contrast, both nDNA data sets (white gene sequences with or without the retention of the 4th intron, and ITS2 sequences and length) detected a single A. marajoara lineage. Conclusions: Strong support for combined data with significant differentiation detected in the COI and absent in the nDNA suggest that the divergence is recent, and detectable only by the faster evolving mtDNA. A within subgenus threshold of >2% may be more appropriate among sister taxa in cryptic anopheline complexes than the standard 3%. Differences in demographic history and climatic changes may have contributed to mtDNA lineage divergence in A. marajoara.