The novel mitochondrial gene arrangement of the cattle tick, Boophilus microplus: Fivefold tandem repetition of a coding region

We sequenced across all of the gene boundaries in the mitochondrial genome of the cattle tick, Boophilus microplus, to determine the arrangement of its genes. The mtDNA of B. microplus has a coding region, composed of tRNA(Glu) and 60 bp of the 3' end of ND1, that is repeated five times. Boophilus microplus is the first coelomate animal known to have more than two copies of a coding sequence. The mitochondrial genome of B, microplus has other unusual features, including (1) reduced T arms in tRNAs, (2) an AT bias in codon use, (3) two control regions that have evolved in concert, (4) three gene rearrangements, and (5) a stem-loop between tRNA(Gln) and tRNA(Phe). The short T arms and small control regions (CRs) of B. microplus and other ticks suggest strong selection for small genomes. Imprecise termination of replication beyond its origin, which can account for the evolution of tandem repeats of coding regions in other mitochondrial genomes, cannot explain the evolution of the fivefold repeated sequence in the mitochondrial genome of B. microplus. Instead, slipped-strand mispairing or recombination are the most plausible explanations for the evolution of these tandem repeats.

Mitochondrial genome data alone are not enough to unambiguously resolve the relationships of Entognatha, Insecta and Crustacea sensu lato (Arthropoda)

An analysis of the relationships of the major arthropod groups Was undertaken using mitochondrial genome data to examine the hypotheses that Hexapoda is polyphyletic and that Collembola is more closely related to branchiopod crustaceans than insects. We sought to examine the sensitivity of this relationship to outgroup choice, data treatment. gene choice and optimality criteria used in the phylogenetic analysis of mitochondrial genome data. Additionally we sequenced the mitochondrial genome of ail archaeognathan, Nesomachilis australica. to improve taxon selection in the apterygote insects, a group poorly represented in previous mitochondrial phylogenies. The sister group of the Collembola was rarely resolved in our analyses with a significant level of support. The use of different outgroups (myriapods, nematodes, or annelids + mollusks) resulted in many different placements of Collembola. The way in which the dataset was coded for analysis (DNA, DNA with the exclusion of third codon position and as amino acids) also had marked affects on tree topology. We found that nodal Support was spread evenly throughout the 13 mitochondrial genes and the exclusion of genes resulted in significantly less resolution in the inferred trees. Optimality criteria had a much lesser effect on topology than the preceding factors; parsimony and Bayesian trees for a given data set and treatment were quite similar. We therefore conclude that the relationships of the extant arthropod groups as inferred by mitochondrial genomes are highly vulnerable to outgroup choice, data treatment and gene choice, and no consistent alternative hypothesis of Collembola's relationships is supported. Pending the resolution of these identified problems with the application of mitogenomic data to basal arthropod relationships, it is difficult to justify the rejection of hexapod monophyly, which is well supported on morphological grounds. (c) The Willi Hennig Society 2004.

Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera)

The complete arrangement of genes in the mitochondrial (mt) genome is known for 12 species of insects, and part of the gene arrangement in the mt genome is known for over 300 other species of insects. The arrangement of genes in the mt genome is very conserved in insects studied, since all of the protein-coding and rRNA genes and most of the tRNA genes are arranged in the same way. We sequenced the entire mt genome of the wallaby louse, Heterodoxus macropus, which is 14,670 bp long and has the 37 genes typical of animals and some noncoding regions. The largest noncoding region is 73 bp long (93% A+T), and the second largest is 47 bp long (92% AST). Both of these noncoding regions seem to be able to form stem-loop structures. The arrangement of genes in the mt genome of this louse is unlike that of any other animal studied. All tRNA genes have moved and/or inverted relative to the ancestral gene arrangement of insects, which is present in the fruit fly Drosophila yakuba. At least nine protein-coding genes (atp6, atp8, cox2, cob, nad1-nad3, nad5, and nad6) have moved; moreover, four of these genes (atp6, atp8, nad1, and nad3) have inverted. The large number of gene rearrangements in the mt genome of H. macropus is unprecedented for an arthropod.

Complete DNA sequence and gene organization of the mitochondrial genome of the liverfluke, Fasciola hepatica L. (Platyhelminthes; Trematoda)

The complete nucleotide sequence of the mitochondrial (mt) DNA molecule of the liverfluke, Fasciola hepatica (phylum Platyhelminthes, class Trematoda, family Fasciolidae), was determined, It comprises 14462 bp, contains 12 protein-encoding, 2 ribosomal and 22 transfer RNA genes, and is the second complete flatworm (and the first trematode) mitochondrial sequence to be described in detail. All of the genes are transcribed from the same strand. Of the genes typically found in mitochondrial genomes of eumetazoans, only atp8 is absent. The nad4L and nad4 genes overlap by 40 nt. Most intergenic sequences are very short. Two larger non-coding regions are present. The longer one (817 nt) is located between trnG and cox3 and consists of 8 identical tandem repeats of 85 nt, rich in G and C, followed by 1 imperfect repeat. The shorter non-coding region (187 nt) exhibits no special features and is separated from the longer region by trnG. The gene arrangement resembles that of some other trematodes including the eastern Asian Schistosoma species (and cyclophyllidean cestode species) but it is strikingly different from that of the African schistosomes, represented by Schistosoma mansoni. The genetic code is as inferred previously for flatworms. Transfer RNA genes range in length from 58 to 70 nt, their products producing characteristic 'clover leaf' structures, except for tRNA(S-VON) and tRNA(S-AGN) lacking the DHU arm.

The value of idiosyncratic markers and changes to conserved tRNA sequences from the mitochondrial genome of hard ticks (Acari: Ixodida: Ixodidae) for phylogenetic inference

Idiosyncratic markers are features of genes and genomes that are so unusual that it is unlikely that they evolved more than once in a lineage of organisms. Here we explore further the potential of idiosyncratic markers and changes to typically conserved tRNA sequences for phylogenetic inference. Hard ticks were chosen as the model group because their phylogeny has been studied extensively. Fifty-eight candidate markers from hard ticks ( family Ixodidae) and 22 markers from the subfamily Rhipicephalinae sensu lato were mapped onto phylogenies of these groups. Two of the most interesting markers, features of the secondary structure of two different tRNAs, gave strong support to the hypothesis that species of the Prostriata ( Ixodes spp.) are monophyletic. Previous analyses of genes and morphology did not strongly support this relationship, instead suggesting that the Prostriata is paraphyletic with respect to the Metastriata ( the rest of the hard ticks). Parallel or convergent evolution was not found in the arrangements of mitochondrial genes in ticks nor were there any reversals to the ancestral arthropod character state. Many of the markers identified were phylogenetically informative, whereas others should be informative with study of additional taxa. Idiosyncratic markers and changes to typically conserved nucleotides in tRNAs that are phylogenetically informative were common in this data set, and thus these types of markers might be found in other organisms.

Novel organization of the mitochondrial genome in the deep-sea coral, Madrepora oculata (Hexacorallia, Scleractinia, Oculinidae) and its taxonomic implications

Madrepora is one of the most ecologically important genera of reef-building scleractinians in the deep sea, occurring from tropical to high-latitude regions. Despite this, the taxonomic affinities and relationships within the genus Madrepora remain unclear. To clarify these issues, we sequenced the mitochondrial (mt) genome of the most widespread Madrepora species, M. oculata, and compared this with data for other scleractinians. The architecture of the M. oculara mt genome was very similar to that of other scleractinians, except for a novel gene rearrangement affecting only cox2 and cox3. This pattern of gene organization was common to four geographically distinct M. oculata individuals as well as the congeneric species M. minutiseptum, but was not shared by other genera that are closely related on the basis of cox1 sequence analysis nor other oculinids, suggesting that it might be unique to Madrepora. (C) 2012 Elsevier Inc. All rights reserved.

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver.

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.

Organization of lin genes and IS6100 among different strains of hexachlorocyclohexane-degrading Sphingomonas paucimobilis: evidence for horizontal gene transfer.

The organization of lin genes and IS6100 was studied in three strains of Sphingomonas paucimobilis (B90A, Sp+, and UT26) which degraded hexachlorocyclohexane (HCH) isomers but which had been isolated at different geographical locations. DNA-DNA hybridization data revealed that most of the lin genes in these strains were associated with IS6100, an insertion sequence classified in the IS6 family and initially found in Mycobacterium fortuitum. Eleven, six, and five copies of IS6100 were detected in B90A, Sp+, and UT26, respectively. IS6100 elements in B90A were sequenced from five, one, and one regions of the genomes of B90A, Sp+, and UT26, respectively, and were found to be identical. DNA-DNA hybridization and DNA sequencing of cosmid clones also revealed that S. paucimobilis B90A contains three and two copies of linX and linA, respectively, compared to only one copy of these genes in strains Sp+ and UT26. Although the copy number and the sequence of the remaining genes of the HCH degradative pathway (linB, linC, linD, and linE) were nearly the same in all strains, there were striking differences in the organization of the linA genes as a result of replacement of portions of DNA sequences by IS6100, which gave them a strange mosaic configuration. Spontaneous deletion of linD and linE from B90A and of linA from Sp+ occurred and was associated either with deletion of a copy of IS6100 or changes in IS6100 profiles. The evidence gathered in this study, coupled with the observation that the G+C contents of the linA genes are lower than that of the remaining DNA sequence of S. paucimobilis, strongly suggests that all these strains acquired the linA gene through horizontal gene transfer mediated by IS6100. The association of IS6100 with the rest of the lin genes further suggests that IS6100 played a role in shaping the current lin gene organization.

Nuclear Domains Enriched in RNA 3′-processing Factors Associate with Coiled Bodies and Histone Genes in a Cell Cycle–dependent Manner

Nuclear domains, called cleavage bodies, are enriched in the RNA 3′-processing factors CstF 64 kDa and and CPSF 100 kDa. Cleavage bodies have been found either overlapping with or adjacent to coiled bodies. To determine whether the spatial relationship between cleavage bodies and coiled bodies was influenced by the cell cycle, we performed cell synchronization studies. We found that in G1 phase cleavage bodies and coiled bodies were predominantly coincident, whereas in S phase they were mostly adjacent to each other. In G2 cleavage bodies were often less defined or absent, suggesting that they disassemble at this point in the cell cycle. A small number of genetic loci have been reported to be juxtaposed to coiled bodies, including the genes for U1 and U2 small nuclear RNA as well as the two major histone gene clusters. Here we show that cleavage bodies do not overlap with small nuclear RNA genes but do colocalize with the histone genes next to coiled bodies. These findings demonstrate that the association of cleavage bodies and coiled bodies is both dynamic and tightly regulated and suggest that the interaction between these nuclear neighbors is related to the cell cycle–dependent expression of histone genes.

Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA

Two small RNAs regulate the timing of Caenorhabditis elegans development(1,2). Transition from the first to the second larval stage fates requires the 22-nucleotide lin-4 RNA(1,3,4), and transition from late larval to adult cell fates requires the 21-nucleotide let-7 RNA 2. The lin-4 and let-7 RNA genes are not homologous to each other, but are each complementary to sequences in the 3' untranslated regions of a set of protein-coding target genes that are normally negatively regulated by the RNAs1,2,5,6. Here we have detected let-7 RNAs of similar to 21 nucleotides in samples from a wide range of animal species, including vertebrate, ascidian, hemichordate, mollusc, annelid and arthropod, but not in RNAs from several cnidarian and poriferan species, Saccharomyces cerevisiae, Escherichia coli or Arabidopsis. We did not detect lin-4 RNA in these species. We found that let-7 temporal regulation is also conserved: let-7 RNA expression is first detected at late larval stages in C. elegans and Drosophila, at 48 hours after fertilization in zebrafish, and in adult stages of annelids and molluscs. The let-7 regulatory RNA may control late temporal transitions during development across animal phylogeny.

Mitochondrial gene content, arrangement and composition compared in African and Asian schistosomes

Complete sequences were obtained for the coding portions of the mitochondrial (mt) genomes of Schistosoma mansoni (NMRI strain, Puerto Rico; 14415 bp), S. japonicum (Anhui strain, China; 14085 bp) and S. mekongi (Khong Island, Laos; 14072 bp). Each comprises 36 genes: 12 protein-encoding genes (cox1-3, nad1-6, nad4L, atp6 and cob); two ribosomal RNAs, rrnL (large subunit rRNA or 16S) and rrnS (small subunit rRNA or 12S); as well as 22 transfer RNA (tRNA) genes. The atp8 gene is absent. A large segment (9.6 kb) of the coding region (comprising 14 tRNAs, eight complete and two incomplete protein-encoding genes) for S. malayensis (Baling, Malaysian Peninsula) was also obtained. Each genome also possesses a long non-coding region that is divided into two parts (a small and a large non-coding region, the latter not fully sequenced in any species) by one or more tRNAs. The protein-encoding genes are similar in size, composition and codon usage in all species except for cox1 in S. mansoni (609 aa) and cox2 in S. mekongi (219 an), both of which are longer than homologues in other species. An unexpected finding in all the Schistosoma species was the presence of a leucine zipper motif in the nad4L gene. The gene order in S. mansoni is strikingly different from that seen in the S. japonicum group and other flatworms. There is a high level of identity (87-94% at both the nucleotide and amino acid levels) for all protein-encoding genes of S. mekongi and S. malayensis. The identity between genes of these two species and those of S. japonicum is less (56-83% for amino acids and 73-79 for nucleotides). The identity between the genes of S. mansoni and the Asian schistosomes is far less (33-66% for amino acids and 54-68% for nucleotides), an observation consistent with the known phylogenetic distance between S. mansoni and the other species. (C) 2001 Elsevier Science B.V. All rights reserved.

The highly rearranged mitochondrial genome of the plague thrips, Thrips imaginis (Insecta : thysanoptera): Convergence of two novel gene boundaries and an extraordinary arrangement of rRNA genes

To help understand the mechanisms of gene rearrangement in the mitochondrial (mt) genomes of hemipteroid insects, we sequenced the mt genome of the plague thrips, Thrips imaginis (Thysanoptera). This genome is circular, 15,407 by long, and has many unusual features, including (1) rRNA genes inverted and distant from one another, (2) an extra gene for tRNA-Ser, (3) a tRNA-Val lacking a D-arm, (4) two pseudo-tRNA genes, (5) duplicate control regions, and (6) translocations and/or inversions of 24 of the 37 genes. The mechanism of rRNA gene transcription in T. imaginis may be different from that of other arthropods since the two rRNA genes have inverted and are distant from one another. Further, the rRNA genes are not adjacent or even close to either of the two control regions. Tandem duplication and deletion is a plausible model for the evolution of duplicate control regions and for the gene translocations, but intramitochondrial recombination may account for the gene inversions in T. imaginis. All the 18 genes between control regions #1 and #2 have translocated and/or inverted, whereas only six of the 20 genes outside this region have translocated and/or inverted. Moreover, the extra tRNA gene and the two pseudo-tRNA genes are either in this region or immediately adjacent to one of the control regions. These observations suggest that tandem duplication and deletion may be facilitated by the duplicate control regions and may have occurred a number of times in the lineage leading to T. imaginis. T. imaginis shares two novel gene boundaries with a lepidopsocid species from another order of hemipteroid insects, the Psocoptera. The evidence available suggests that these shared gene boundaries evolved by convergence and thus are not informative for the interordinal phylogeny of hemipteroid insects. We discuss the potential of hemipteroid insects as a model system for studies of the evolution of animal rut genomes and outline some fundamental questions that may be addressed with this system.

Dominant negative MyD88 proteins inhibit interleukin-1beta /interferon-gamma -mediated induction of nuclear factor kappa B-dependent nitrite production and apoptosis in beta cells.

Insulin-dependent diabetes mellitus is an autoimmune disease in which pancreatic islet beta cells are destroyed by a combination of immunological and inflammatory mechanisms. In particular, cytokine-induced production of nitric oxide has been shown to correlate with beta cell apoptosis and/or inhibition of insulin secretion. In the present study, we investigated whether the interleukin (IL)-1beta intracellular signal transduction pathway could be blocked by overexpression of dominant negative forms of the IL-1 receptor interacting protein MyD88. We show that overexpression of the Toll domain or the lpr mutant of MyD88 in betaTc-Tet cells decreased nuclear factor kappaB (NF-kappaB) activation upon IL-1beta and IL-1beta/interferon (IFN)-gamma stimulation. Inducible nitric oxide synthase mRNA accumulation and nitrite production, which required the simultaneous presence of IL-1beta and IFN-gamma, were also suppressed by approximately 70%, and these cells were more resistant to cytokine-induced apoptosis as compared with parental cells. The decrease in glucose-stimulated insulin secretion induced by IL-1beta and IFN-gamma was however not prevented. This was because these dysfunctions were induced by IFN-gamma alone, which decreased cellular insulin content and stimulated insulin exocytosis. These results demonstrate that IL-1beta is involved in inducible nitric oxide synthase gene expression and induction of apoptosis in mouse beta cells but does not contribute to impaired glucose-stimulated insulin secretion. Furthermore, our data show that IL-1beta cellular actions can be blocked by expression of MyD88 dominant negative proteins and, finally, that cytokine-induced beta cell secretory dysfunctions are due to the action of IFN-gamma.

Molecular characterisation of intermediate snail hosts and the search for resistance genes

The relationship between schistosomes and their intermediate hosts is an extremely intricate one with strains and species of the parasite depending on particular species of snail, which in turn may vary in their susceptibility to the parasites. In order to gain a better understanding of the epidemiology of the disease we have been investigating the use of molecular markers for snail identification and for studying host-parasite relationships. In this paper we will draw on examples concerning schistosomiasis in West and East Africa to illustrate how a molecular analysis can be used as part of a "total evidence" approach to characterisation of Bulinus species and provide insights into parasite transmission. Particular emphasis is given to ribosomal RNA genes (rRNA), random amplified polymorphic DNA (RAPDs) and the mitochondrial gene cytochrome oxidase I (COI). Snails resistant to infection occur naturally and there is a genetic basis for this resistance. In Biomphalaria glabrata resistance to Schistosoma mansoni is known to be a polygenic trait and we have initiated a preliminary search for snail genomic regions linked to, or involved in, resistance by using a RAPD based approach in conjunction with progeny pooling methods. We are currently characterising a variety of STSs (sequence tagged sites) associated with resistance. These can be used for local linkage and interval mapping to define genomic regions associated with the resistance trait. The development of such markers into simple dot-blot or specific PCR-based assays may have a direct and practical application for the identification of resistant snails in natural populations.

Identification of Biomphalaria havanensis and Biomphalaria obstructa populations from Cuba using polymerase chain reaction and restriction fragment length polymorphism of the ribosomal RNA intergenic spacer

In Cuba, several Biomphalaria species have been reported such as B. orbignyi, B. schrammi, B. helophila, B. havanensis and B. peregrina; only the latter three are considered as potential hosts of Schistosoma mansoni. The specific identification of Biomphalaria species is based on anatomical and morphological characters of genital organs and shells. The correct identification of these snails is complicated by the high variation in these characters, similarity among species and in some cases by the small size of the snails. In this paper, we reported the classical morphological identification, the use of PCR and RFLP analysis of the internal transcribed spacer region of the ribosomal RNA genes for molecular identification of seven snail populations from different localities in Cuba. Using morphological and molecular analysis, we showed that among the studied Cuban Biomphalaria populations only B. havanensis and B. obstructa species were found.