Genetic Variability and Geographical Diversity of the Main Chagas' Disease Vector Panstrongylus megistus (Hemiptera: Triatominae) in Brazil Based on Ribosomal DNA Intergenic Sequences

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

Genetic Variability of Beauveria bassiana and a DNA Marker for Environmental Monitoring of a Highly Virulent Isolate Against Cosmopolites sordidus

The banana weevil Cosmopolites sordidus (Germar) is one of a number of pests that attack banana crops. The use of the entomopathogenic fungus Beauveria bassiana as a biological control agent for this pest may contribute towards reducing the application of chemical insecticides on banana crops. In this study, the genetic variability of a collection of Brazilian isolates of B. bassiana was evaluated. Samples were obtained from various geographic regions of Brazil, and from different hosts of the Curculionidae family. Based on the DNA fingerprints generated by RAPD and AFLP, we found that 92 and 88 % of the loci were polymorphic, respectively. The B. bassiana isolates were attributed to two genotypic clusters based on the RAPD data, and to three genotypic clusters, when analyzed with AFLP. The nucleotide sequences of nuclear ribosomal DNA intergenic spacers confirmed that all isolates are in fact B. bassiana. Analysis of molecular variance showed that variability among the isolates was not correlated with geographic origin or hosts. A RAPD-specific marker for isolate CG 1024, which is highly virulent to C. sordidus, was cloned and sequenced. Based on the sequences obtained, specific PCR primers BbasCG1024F (5'-TGC GGC TGA GGA GGA CT-3') and BbasCG1024R (5'-TGC GGC TGA GTG TAG AAC-3') were designed for detecting and monitoring this isolate in the field.

Unprecedented loss of ammonia assimilation capability in a urease-encoding bacterial mutualist.

BACKGROUND: Blochmannia are obligately intracellular bacterial mutualists of ants of the tribe Camponotini. Blochmannia perform key nutritional functions for the host, including synthesis of several essential amino acids. We used Illumina technology to sequence the genome of Blochmannia associated with Camponotus vafer. RESULTS: Although Blochmannia vafer retains many nutritional functions, it is missing glutamine synthetase (glnA), a component of the nitrogen recycling pathway encoded by the previously sequenced B. floridanus and B. pennsylvanicus. With the exception of Ureaplasma, B. vafer is the only sequenced bacterium to date that encodes urease but lacks the ability to assimilate ammonia into glutamine or glutamate. Loss of glnA occurred in a deletion hotspot near the putative replication origin. Overall, compared to the likely gene set of their common ancestor, 31 genes are missing or eroded in B. vafer, compared to 28 in B. floridanus and four in B. pennsylvanicus. Three genes (queA, visC and yggS) show convergent loss or erosion, suggesting relaxed selection for their functions. Eight B. vafer genes contain frameshifts in homopolymeric tracts that may be corrected by transcriptional slippage. Two of these encode DNA replication proteins: dnaX, which we infer is also frameshifted in B. floridanus, and dnaG. CONCLUSIONS: Comparing the B. vafer genome with B. pennsylvanicus and B. floridanus refines the core genes shared within the mutualist group, thereby clarifying functions required across ant host species. This third genome also allows us to track gene loss and erosion in a phylogenetic context to more fully understand processes of genome reduction.

Discordances between phylogenetic and morphological patterns in alpine leaf beetles attest to an intricate biogeographic history of lineages in postglacial Europe.

Pleistocene glacial and interglacial periods have moulded the evolutionary history of European cold-adapted organisms. The role of the different mountain massifs has, however, not been accurately investigated in the case of high-altitude insect species. Here, we focus on three closely related species of non-flying leaf beetles of the genus Oreina (Coleoptera, Chrysomelidae), which are often found in sympatry within the mountain ranges of Europe. After showing that the species concept as currently applied does not match barcoding results, we show, based on more than 700 sequences from one nuclear and three mitochondrial genes, the role of biogeography in shaping the phylogenetic hypothesis. Dating the phylogeny using an insect molecular clock, we show that the earliest lineages diverged more than 1 Mya and that the main shift in diversification rate occurred between 0.36 and 0.18 Mya. By using a probabilistic approach on the parsimony-based dispersal/vicariance framework (MP-DIVA) as well as a direct likelihood method of state change optimization, we show that the Alps acted as a cross-roads with multiple events of dispersal to and reinvasion from neighbouring mountains. However, the relative importance of vicariance vs. dispersal events on the process of rapid diversification remains difficult to evaluate because of a bias towards overestimation of vicariance in the DIVA algorithm. Parallels are drawn with recent studies of cold-adapted species, although our study reveals novel patterns in diversity and genetic links between European mountains, and highlights the importance of neglected regions, such as the Jura and the Balkanic range.

Ribosomal DNA ITS-1 intergenic spacer polymorphism in triatomines (Triatominae, Heteroptera)

The length polymorphism of ribosomal DNA ITS-1 intergenic spacer was analyzed in eight species of triatomines belonging to Triatoma, Rhodnius, and Panstrongylus genera. The analyzed species were Rhodnius domesticus, R. neivai, R. robustus, Triatoma brasiliensis, T. infestans, T. vitticeps, Panstrongylus megistus, and P. herreri. These insects are vectors of Chagas' disease, one of the most prominent public health problems among South American countries. This work allowed the differentiation between species of the Triatomini and Rhodniini tribes through the analysis of ITS-1 length polymorphism by PCR and RFLP techniques. The species of the Triatoma and Panstrongylus genera presented an amplified ITS-1 fragment between 600 and 1000 bp, whereas Rhodnius presented a less variable ITS-1 length fragment, around 300 bp, which could reflect the monophyletic origin of the Rhodniini tribe. Species belonging to this genus were further differentiated by RFLP with HaeIII and AluI endonucleases. Our results corroborate the hypothesis of polyphyletic origin in this group of insects and contribute to knowledge about evolutionary relationships in triatomines.

Identification of long intergenic region sequences involved in maize streak virus replication

The main cis-acting control regions for replication of the single-stranded DNA genome of maize streak virus (MSV) are believed to reside within an approximately 310 nt long intergenic region (LIR). However, neither the minimum LIR sequence required nor the sequence determinants of replication specificity have been determined experimentally. There are iterated sequences, or iterons, both within the conserved inverted-repeat sequences with the potential to form a stem-loop structure at the origin of virion-strand replication, and upstream of the rep gene TATA box (the rep-proximal iteron or RPI). Based on experimental analyses of similar iterons in viruses from other geminivirus genera and their proximity to known Rep-binding sites in the distantly related mastrevirus wheat dwarf virus, it has been hypothesized that the iterons may be Rep-binding and/or -recognition sequences. Here, a series of LIR deletion mutants was used to define the upper bounds of the LIR sequence required for replication. After identifying MSV strains and distinct mastreviruses with incompatible replication-specificity determinants (RSDs), LIR chimaeras were used to map the primary MSV RSD to a 67 nt sequence containing the RPI. Although the results generally support the prevailing hypothesis that MSV iterons are functional analogues of those found in other geminivirus genera, it is demonstrated that neither the inverted-repeat nor RPI sequences are absolute determinants of replication specificity. Moreover, widely divergent mastreviruses can trans-replicate one another. These results also suggest that sequences in the 67 nt region surrounding the RPI interact in a sequence-specific manner with those of the inverted repeat.

Specificity of the ModA11, ModA12 and ModD1 epigenetic regulator N-6-adenine DNA methyltransferases of Neisseria meningitidis

Phase variation (random ON/OFF switching) of gene expression is a common feature of host-adapted pathogenic bacteria. Phase variably expressed N-6-adenine DNA methyltransferases (Mod) alter global methylation patterns resulting in changes in gene expression. These systems constitute phase variable regulons called phasevarions. Neisseria meningitidis phasevarions regulate genes including virulence factors and vaccine candidates, and alter phenotypes including antibiotic resistance. The target site recognized by these Type III N-6-adenine DNA methyltransferases is not known. Single molecule, real-time (SMRT) methylome analysis was used to identify the recognition site for three key N. meningitidis methyltransferases: ModA11 (exemplified by M.NmeMC58I) (5'-CGY(m6)AG-3'), ModA12 (exemplified by M.Nme77I, M.Nme18I and M.Nme579II) (5'-AC(m6)ACC-3') and ModD1 (exemplified by M.Nme579I) (5'-CC(m6)AGC-3'). Restriction inhibition assays and mutagenesis confirmed the SMRT methylome analysis. The ModA11 site is complex and atypical and is dependent on the type of pyrimidine at the central position, in combination with the bases flanking the core recognition sequence 5'-CGY(m6)AG-3'. The observed efficiency of methylation in the modA11 strain (MC58) genome ranged from 4.6% at 5'-GCGC(m6)AGG-3' sites, to 100% at 5'-ACGT(m6)AGG-3' sites. Analysis of the distribution of modified sites in the respective genomes shows many cases of association with intergenic regions of genes with altered expression due to phasevarion switching.

Genome-wide high-resolution mapping and functional analysis of DNA methylation in arabidopsis.

Cytosine methylation is important for transposon silencing and epigenetic regulation of endogenous genes, although the extent to which this DNA modification functions to regulate the genome is still unknown. Here we report the first comprehensive DNA methylation map of an entire genome, at 35 base pair resolution, using the flowering plant Arabidopsis thaliana as a model. We find that pericentromeric heterochromatin, repetitive sequences, and regions producing small interfering RNAs are heavily methylated. Unexpectedly, over one-third of expressed genes contain methylation within transcribed regions, whereas only approximately 5% of genes show methylation within promoter regions. Interestingly, genes methylated in transcribed regions are highly expressed and constitutively active, whereas promoter-methylated genes show a greater degree of tissue-specific expression. Whole-genome tiling-array transcriptional profiling of DNA methyltransferase null mutants identified hundreds of genes and intergenic noncoding RNAs with altered expression levels, many of which may be epigenetically controlled by DNA methylation.

Epigenetic remodeling of meiotic crossover frequency in Arabidopsis thaliana DNA methyltransferase mutants.

Meiosis is a specialized eukaryotic cell division that generates haploid gametes required for sexual reproduction. During meiosis, homologous chromosomes pair and undergo reciprocal genetic exchange, termed crossover (CO). Meiotic CO frequency varies along the physical length of chromosomes and is determined by hierarchical mechanisms, including epigenetic organization, for example methylation of the DNA and histones. Here we investigate the role of DNA methylation in determining patterns of CO frequency along Arabidopsis thaliana chromosomes. In A. thaliana the pericentromeric regions are repetitive, densely DNA methylated, and suppressed for both RNA polymerase-II transcription and CO frequency. DNA hypomethylated methyltransferase1 (met1) mutants show transcriptional reactivation of repetitive sequences in the pericentromeres, which we demonstrate is coupled to extensive remodeling of CO frequency. We observe elevated centromere-proximal COs in met1, coincident with pericentromeric decreases and distal increases. Importantly, total numbers of CO events are similar between wild type and met1, suggesting a role for interference and homeostasis in CO remodeling. To understand recombination distributions at a finer scale we generated CO frequency maps close to the telomere of chromosome 3 in wild type and demonstrate an elevated recombination topology in met1. Using a pollen-typing strategy we have identified an intergenic nucleosome-free CO hotspot 3a, and we demonstrate that it undergoes increased recombination activity in met1. We hypothesize that modulation of 3a activity is caused by CO remodeling driven by elevated centromeric COs. These data demonstrate how regional epigenetic organization can pattern recombination frequency along eukaryotic chromosomes.

An improved PCR method for direct identification of Porphyra (Bangiales, Rhodophyta) using conchocelis based on a RUBISCO intergenic spacer

An improved method of PCR in which the small segment of conchocelis is amplified directly without DNA extraction was used to amplify a RUBISCO intergenic spacer DNA fragment from nine species of red algal genus Porphyra (Bangiales, Rhodophyta), including Porphyra yezoensis (Jiangsu, China), P. haitanensis (Fujian, China), P. oligospermatangia (Qingdao, China), P. katadai (Qingdao, China), P. tenera (Qingdao, China), P. suborboculata (Fujian, China), P. pseudolinearis (Kogendo, Korea), P. linearis (Devon, England), and P. fallax (Seattle, USA). Standard PCR and the method developed here were both conducted using primers specific for the RUBISCO spacer region, after which the two PCR products were sequenced. The sequencing data of the amplicons obtained using both methods were identical, suggesting that the improved PCR method was functional. These findings indicate that the method developed here may be useful for the rapid identification of species of Porphyra in a germplasm bank. In addition, a phylogenetic tree was constructed using the RUBISCO spacer and partial rbcS sequence, and the results were in concordant with possible alternative phylogenies based on traditional morphological taxonomic characteristics, indicating that the RUBISCO spacer is a useful region for phylogenetic studies.

Potential refugium on the Qinghai-Tibet Plateau revealed by the chloroplast DNA phylogeography of the alpine species Metagentiana striata (Gentianaceae)

Metagentiana striata is an alpine annual herbaceous plant endemic to the east of the Qinghai-Tibet (Q-T) Plateau and adjacent areas. The phylogeography of M. striata was studied by sequencing the chloroplast DNA (cpDNA) trnS-trnG intergenic spacer. Ten haplotypes were identified from an investigation of 232 individuals of M. striata from 14 populations covering the entire geographical range of this species. The level of differentiation amongst populations was very high (G(ST) = 0.746; N-ST = 0.774) and a significant phylogeographical structure was observed (P < 0.05). An analysis of molecular variance found a high variation amongst populations (76%), with F-ST = 0.762 (highly significant, P < 0.001), indicating that little gene flow occurred amongst the different regions; this was explained by the isolation of populations by high mountains along the Q-T Plateau and adjacent areas (N-m = 0.156). Only one ancestral haplotype (A) was common and widespread throughout the distributional range of M. striata. The populations of the Hengduan Mountains region of the south-eastern Q-T Plateau showed high diversity and uniqueness of haplotypes. It is suggested that this region was the potential refugium of M. striata during the Quaternary glaciation, and that interglacial and postglacial range expansion occurred from this refugium. This scenario was in good agreement with the results of nested clade analysis, which inferred that the current spatial distribution of cpDNA haplotypes and populations resulted from range expansion, together with past allopatric fragmentation events. (c) 2008 The Linnean Society of London.

DNA barcoding of a large genus, Aspalathus L. (Fabaceae)

The DNA barcode potential of three regions (the nuclear ribosomal ITS and the plastid psbA-trnH and trnT-trnL intergenic spacers) was investigated for the plant genus Aspalathus L. (Fabaceac: Crotalarieae). Aspalathus is a large genus (278 species) that revealed low levels of DNA variation in phylogenetic studies. In a 51-species dataset for the psbA-trnH and ITS regions, 45%, and 16% of sequences respectively were identical to the sequence of at least one other species, with two species undiscriminated even when the two regions were combined. In contrast, trnT-trnL, discriminated between all species in this dataset. In a larger ITS and trnT-trnL dataset. including a further 82 species. 7 species in five pairwise comparisons remained Undiscriminated when the two regions were combined. Four of the five pairs of species not discriminated by sequence data were readily distinguished using a combination of qualitative and quantitative morphological data. The difficulty of barcoding in this group is increased by the presence of intraspecific variation in all three regions studied. In the case of psbA-trnH, three intraspecific samples had a sequence identical to at least one other species. Overall, psbA-trnH. currently a candidate for plant barcoding, was the least discriminatory region in our study.

Variation in 16S-23S rRNA intergenic spacer regions in Photobacterium damselae: a mosaic-like structure

Phenotypically, Photobacterium damselae subsp. piscicida and P. damselae subsp. damselae are easily distinguished. However, their 16S rRNA gene sequences are identical, and attempts to discriminate these two subspecies by molecular tools are hampered by their high level of DNA-DNA similarity. The 16S-23S rRNA internal transcribed spacers (ITS) were sequenced in two strains of Photobacterium damselae subsp. piscicida and two strains of P. damselae subsp. damselae to determine the level of molecular diversity in this DNA region. A total of 17 different ITS variants, ranging from 803 to 296 bp were found, some of which were subspecies or strain specific. The largest ITS contained four tRNA genes (tDNAs) coding for tRNA(Glu(UUC)), tRNA(LyS(UUU)), tRNA(Val(UAC)), and tRNA(Ala(GGC)). Five amplicons contained tRNA(Glu(UUC)) combined with two additional tRNA genes, including tRNA(Lys(UUU)), tRNA(Val(UAC)), or tRNA(Ala(UGC)). Five amplicons contained tRNA(Ile(GAU)) and tRNA(Ala(UGC)). Two amplicons contained tRNA(Glu(UUC)) and tRNA(Val(UGC)). Two different isoacceptor tRNA(Ala) genes (GGC and UGC anticodons) were found. The five smallest amplicons contained no tRNA genes. The tRNA-gene combinations tRNA(Glu(UUC)) -tRNA(Val(UAC)) -tRNA(Ala(UGC)) and tRNA(Glu(UUC)) -tRNA(Ala(UGC)) have not been previously reported in bacterial ITS regions. The number of copies of the ribosomal operon (rrn) in the P. damselae chromosome ranged from at least 9 to 12. For ITS variants coexisting in two strains of different subspecies or in strains of the same subspecies, nucleotide substitution percentages ranged from 0 to 2%. The main source of variation between ITS variants was due to different combinations of DNA sequence blocks, constituting a mosaic-like structure.

Complete mitochondrial DNA sequences of the decapod crustaceans pseudocarcinus gigas (Menippidae) and macrobrachium rosenbergii (Palaemonidae)

The complete mitochondrial DNA sequence was determined for the Australian giant crab Pseudocarcinns gigas (Crustacea: Decapoda: Menippidae) and the giant freshwater shrimp Macrobrachium rosenbergii (Crustacea: Decapoda: Palaemonidae). The Pse gigas and Mrosenbergii mitochondrial genomes are circular molecules, 15,515 and 15,772 bp in length, respectively, and have the same gene composition as found in other metazoans. The gene arrangement of M. rosenbergii corresponds with that of the presumed ancestral arthropod gene order, represented by Limulus polyphemus, except for the position of the tRNA^Leu(UUR) gene. The Pse. gigas gene arrangement corresponds exactly with that reported for another brachyuran, Portunus trituberculatus, and differs from the M. rosenbergii gene order by only the position of the tRNA^His gene. Given the relative positions of intergenic nonoding nucleotides, the “duplication/random loss” model appears to be the most plausible mechanism for the translocation of this gene. These data represent the first caridean and only the second brachyuran complete mtDNA sequences, and a source of information that will facilitate surveys of intraspecific variation within these commercially important decapod species.