The multicopy appearance of a large inverted duplication and the sequence at the inversion joint suggest a new model for gene amplification.

The amplified DNA of HC50474, a Chinese hamster fibroblast cell line selected in three steps for high resistance to coformycin, consists chiefly of 150 copies of a large inverted duplication including the adenylate deaminase gene. Most if not all of these units are more than 2 x 120 kb long. The inverted duplication was first detected in the cells recovered from the second selection step, at the same chromosomal location as the first step amplified units. Its formation and amplification appear to be coupled since the second step cell line already contained 40 copies of this novel structure. Reamplification of the inverted duplication occurred at the third step of selection concomitant with the loss of amplified DNA acquired during the first step. The head-to-head junction has been formed by recombination within a recombinational hotspot described previously [Hyrien, O., Debatisse, M., Buttin, G. and Robert de Saint Vincent, B. (1987) EMBO J., 6, 2401-2408]. Sequences at the joint and in the corresponding wild-type region reveal that the crossover sites, one of which occurs in the putative promoter region of B2 repeat, are located at the top of significant stem-loop structures and that patchy homologies between the parental molecules on one side of the breakpoints allow alignment of these crossover sites. We present a model which explains the formation and amplification of this and other large inverted duplications by errors in DNA replication.     

Patterns of gene duplication and intron loss in the ENCODE regions suggest a confounding factor

The exon-intron structure of eukaryotic genes allows for phenomena such as alternative splicing, nonsense-mediated decay, and regulation through untranslated regions. However, the evolution of the exon structure of genes is not well elucidated because of limited and phylogenetically sparse data sets. In this study, we use the phylogenetically diverse sequencing of the ENCODE regions to study gene structure evolution in mammalian genomes. This first phylogenetically diverse study of gene structure changes offers insights into the mode and tempo of mammalian gene structure evolution. The genes undergoing structure changes appear to be moderately to highly expressed in germline cells and show levels of selection similar to those of other ENCODE genes. Patterns of gene duplication of the affected genes are more complex than expected. The number of sampled genomes is sufficiently dense to infer that certain gene duplications happened after intron loss. Thus, although gene duplication is highly correlated with intron loss, we conclude that structural changes in genes are not necessarily due to a loss of constraint following gene duplication as previously suggested.     

Multiple sequence rearrangements accompanying the duplication of a tRNAPro gene in wheat mitochondrial DNA

In the course of isolating tRNA genes from wheat mtDNA, we have found the same tRNA^Pro gene in two different Hind III restriction fragments, H-P1 (0.7 kbp) and H-P2 (1.7 kbp). Sequences immediately flanking these duplicate genes are closely related, although not identical; sequence comparisons suggest that multiple rearrangements have occurred in the vicinity of the H-P2 tRNA^Pro gene, relative to the H-P1 version. The chimeric nature of H-P2 is emphasized by the presence of sequences that are also found upstream of the wheat mitochondrial 26S rRNA gene, as well as sequences derived from chloroplast DNA. Comparison of H-P2 with H-P1 plus upstream sequences provides some insight into possible molecular events that might have generated H-P2. In particular, such comparisons suggest a model in which the homologous sequences in H-P2 are seen to be derived from H-P1 plus upstream sequences as a result of an intragenomic, site-specific rearrangement event, followed by amplification of the product, its fixation in the mitochondrial genome, and subsequent sequence divergence (single base changes as well as insertions/deletions of up to 50 nucleotides). The results reported here implicate particular primary sequence motifs in certain of the rearrangements that characterize H-P2.     

Intron sequences of arginine kinase in an intertidal snail suggest an ecotype-specific selective sweep and a gene duplication

Many species with restricted gene flow repeatedly respond similarly to local selection pressures. To fully understand the genetic mechanisms behind this process, the phylogeographic history of the species (inferred from neutral markers) as well as the loci under selection need to be known. Here we sequenced an intron in the arginine kinase gene (Ark), which shows strong clinal variation between two locally adapted ecotypes of the flat periwinkle, Littorina fabalis. The ‘small-sheltered'' ecotype was almost fixed for one haplotype, H1, in populations on both sides of the North Sea, unlike the ‘large-moderately exposed ecotype'', which segregated for ten different haplotypes. This contrasts with neutral markers, where the two ecotypes are equally variable. H1 could have been driven to high frequency in an ancestral population and then repeatedly spread to sheltered habitats due to local selection pressures with the colonization of both sides of the North Sea, after the last glacial maximum (∼18 000 years ago). An alternative explanation is that a positively selected mutation, in or linked to Ark, arose after the range expansion and secondarily spread through sheltered populations throughout the distribution range, causing this ecotype to evolve in a concerted fashion. Also, we were able to sequence up to four haplotypes consistently from some individuals, suggesting a gene duplication in Ark.     

Complete nucleotide sequences and gene organization of mitochondrial genome of Bufo gargarizans

The complete mitochondrial genome of Bufo gargarizans was sequenced using overlapping polymerase chain reaction (PCR) amplicons (GenBank Accession No. DQ275350). The genome is 17,277 base pairs in length, containing 13 protein-coding genes (ATP6, ATP8, COI-III, ND1-6, ND4L, Cyt b), 2 ribosomal RNAs (12S rRNA and 16S rRNA), 22 transfer RNAs and a putative control region. We analyzed the sequence using bioinformatics methods comparing the obtained mtDNA sequence with other toads and frogs. Based on the concatenated nucleotide sequences of protein-coding genes, we constructed a phylogenetic tree with maximum likelihood (ML) and maximum parsimony (MP) methods and discussed the phylogenetic relationships among 11 species of Anura.     

Complete mitochondrial genomes and novel gene rearrangements in two dicroglossid frogs, Hoplobatrachus tigerinus and Euphlyctis hexadactylus, from Bangladesh

We determined the complete nucleotide sequences of mitochondrial (mt) genomes from two dicroglossid frogs, Hoplobatrachus tigerinus (Indian Bullfrog) and Euphlyctis hexadactylus (Indian Green frog). The genome sizes are 20462 bp in H. tigerinus and 20280 bp in E. hexadactylus. Although both genomes encode the typical 37 mt genes, the following unique features are observed: 1) the ND5 genes are duplicated in H. tigerinus that have completely identical sequences, whereas duplicated ND5 genes in E. hexadactylus possessed dissimilar substitutions; 2) duplicated control region (CR) in H. tigerinus has almost identical sequences whereas single control region (CR) was found in E. hexadactylus; 3) the tRNA-Leu (CUN) gene is translocated from the LTPF tRNA cluster to downstream of ND5-1 in H. tigerinus, and the tRNA-Pro gene is translocated from the LTPF tRNA cluster to downstream of CR in E. hexadactylus; 4) pseudo tRNA-Leu (CUN) and tRNA-Pro genes are observed in E. hexadactylus; and 5) two tRNA-Met genes are encoded in both species, as observed in the previously reported dicroglossid mt genomes. Almost all observed gene rearrangements in H. tigerinus and E. hexadactylus can be explained by the tandem duplication and random loss model, except translocation of tRNA-Pro in E. hexadactylus. The novel mt genomic features found in this study may be useful for future phylogenetic studies in the dicroglossid taxa. However, the mt genome with interesting features found in the present study reveal a high level of variation of gene order and gene content, inspiring more research to understand the mechanisms behind gene and genome evolution in the dicroglossid and as well as in the amphibian taxa in future studies.     

Gene rearrangements in snake mitochondrial genomes: highly concerted evolution of control-region-like sequences duplicated and inserted into a tRNA gene cluster

Mitochondrial DNA (mtDNA) regions corresponding to two major tRNA gene clusters were amplified and sequenced for the Japanese pit viper, himehabu. In one of these clusters, which in most vertebrates characterized to date contains three tightly connected genes for tRNA(Ile), and tRNA(Gln), and tRNA(Met), a sequence of approximately 1.3 kb was found to be inserted between the genes for tRNA(Ile) and tRNA(Gln). The insert consists of a control-region-like sequence possessing some conserved sequence blocks, and short flanking sequences which may be folded into tRNA(Pro), tRNA(Phe), and tRNA(Leu) genes. Several other snakes belonging to different families were also found to possess a control-region-like sequence and tRNA(Leu) gene between the tRNA(Ile)and tRNA(Gln) genes. We also sequenced a region surrounded by genes for cytochrome b and 12S rRNA, where the control region and genes for tRNA(Pro) and tRNA(Phe) are normally located in the mtDNAs of most vertebrates. In this region of three examined snakes, a control-region- like sequence exists that is almost completely identical to the one found between the tRNA(Ile) and tRNA(Gln) genes. The mtDNAs of these snakes thus possess two nearly identical control-region-like sequences which are otherwise divergent to a large extent between the species. These results suggest that the duplicate state of the control-region- like sequences has long persisted in snake mtDNAs, possibly since the original insertion of the control-region-like sequence and tRNA(Leu) gene into the tRNA gene cluster, which occurred in the early stage of the divergence of snakes. It is also suggested that the duplicated control-region-like sequences at two distant locations of mtDNA have evolved concertedly by a mechanism such as frequent gene conversion. The secondary structures of the determined tRNA genes point to the operation of simplification pressure on the T psi C arm of snake mitochondrial tRNAs.      

Complete mitochondrial genome of Chinese big-headed turtle, Platysternon megacephalum, with a novel gene organization in vertebrate mtDNA

The mitochondrial genome of the Chinese big-headed turtle, Platysternon megacephalum, was obtained using polymerase chain reaction (PCR). The entire mtDNA sequence, the longest mitochondrial genome in turtles reported so far, is 19161 bp. This mitochondrial genome exhibits a novel gene order, which greatly differs from that of any other vertebrates. It is characterized by four distinctive features: 1) the translocation of a gene cluster including three tRNA genes (tRNAHis, tRNASer, tRNALeu(CUN)) and ND5 gene, 2) two tRNAThr pseudogenes, 3) a duplication of pseudo tRNAThr/tRNAPro/D-loop region and 4) 3 non-coding spacers. These unique identities represent a new mitogenomic gene order in vertebrates. The TDRL model was proposed to account for the generation of the gene order in P. megacephalum.     

Complete nucleotide sequences of mitochondrial genomes of two solitary entoprocts, Loxocorone allax and Loxosomella aloxiata: implications for lophotrochozoan phylogeny

The complete nucleotide sequences of the mitochondrial (mt) genomes of the entoprocts Loxocorone allax and Loxosomella aloxiata were determined. Both species carry the typical gene set of metazoan mt genomes and have similar organizations of their mt genes. However, they show differences in the positions of two tRNA(Leu) genes. Additionally, the tRNA(Val) gene, and half of the long non-coding region, is duplicated and inverted in the Loxos. aloxiata mt genome. The initiation codon of the Loxos. aloxiata cytochrome oxidase subunit I gene is expected to be ACG rather than AUG. The mt gene organizations in these two entoproct species most closely resemble those of mollusks such as Katharina tunicata and Octopus vulgaris, which have the most evolutionarily conserved mt gene organization reported to date in mollusks. Analyses of the mt gene organization in the lophotrochozoan phyla (Annelida, Brachiopoda, Echiura, Entoprocta, Mollusca, Nemertea, and Phoronida) suggested a close phylogenetic relationship between Brachiopoda, Annelida, and Echiura. However, Phoronida was excluded from this grouping. Molecular phylogenetic analyses based on the sequences of mt protein-coding genes suggested a possible close relationship between Entoprocta and Phoronida, and a close relationship among Brachiopoda, Annelida, and Echiura.     

Maximum-likelihood divergence date estimates based on rRNA gene sequences suggest two scenarios of Trypanosoma cruzi intraspecific evolution.

The phylogenetic relationships of Trypanosoma cruzi strains were inferred using maximum-likelihood from complete 18S rDNA sequences and D7-24Salpha rDNA regions from 20 representative strains of T. cruzi. For this we sequenced the 18S rDNA of 14 strains and the D7-24Salpha rDNA of four strains and aligned them to previously published sequences. Phylogenies inferred from these data sets identified four groups, named Riboclades 1, 2, 3, and 4, and a basal dichotomy that separated Riboclade 1 from Riboclades 2, 3, and 4. Substitution models and other parameters were optimized by hierarchical likelihood tests, and our analysis of the 18S rDNA molecular clock by the likelihood ratio test suggests that a taxa subset encompassing all 2,150 positions in the alignment supports rate constancy among lineages. The present analysis supports the notion that divergence dates of T. cruzi Riboclades can be estimated from 18S rDNA sequences and therefore, we present alternative evolutionary scenarios based on two different views of T. cruzi intraspecific divergence. The first assumes a faster evolutionary rate, which suggests that the divergence between T. cruzi I and II and the extant strains occurred in the Tertiary period (37-18 MYA). The other, which supports the hypothesis that the divergence between T. cruzi I and II occurred in the Cretaceous period (144-65 MYA) and the divergence of the extant strains occurred in the Tertiary period of the Cenozoic era (65-1.8 MYA), is consistent with our previously proposed hypothesis of divergence by geographical isolation and mammalian host coevolution.     

Complete mitochondrial genome of Dacus conopsoides (Insecta: Tephritidae) with tRNA gene duplication and molecular phylogeny of Dacini tribe

To date there is only a single report on the complete mitochondrial genome of the Dacus fruit flies. We report here the whole mitogenome of Dacus conopsoides with first report of tRNA gene duplication in tephritid fruit flies determined using next-generation sequencing and discuss the molecular phylogeny of Dacini tribe. It had a total length of 15,852 bp, comprising 13 protein coding genes, 2 rRNA genes, 23 tRNA genes, and a non-coding region (A + T-rich control region). The 65-bp trnF gene was duplicated, and the 68-bp trnE gene was partially duplicated resulting in a 31-bp pseudogene. The cloverleaf structure for trnN, trnH, and trnF lacked the TΨC-loop, while trnS lacked the D-stem. The start codons for the protein coding genes included 6 ATG, 3 ATC, 2 ATA, and 1 each of ATT and TCG. Seven PCGs had TAA stop codon, two had TAG and four had incomplete T stop codon. Molecular phylogeny based on 15 mt-genes (13 PCGs +2 rRNA genes) and 30 taxa of Tephritidae indicated D. conopsoides forming a monophyletic sister group with D. longicornis supported by high bootstrap value. The lineage containing also the monophyletic genus Zeugodacus. The Dacini and Ceratitidini tribes of the subfamily Dacinae were monophyletic but the subfamilies Dacinae and Trypetinae were paraphyletic. A broader taxa sampling of the Tephritidae is needed to better elucidate the phylogenetics and systematics of the tribes and subfamilies of tephritid fruit flies.     

A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.     

Fructosyltransferase and invertase genes evolved by gene duplication and rearrangements: rice, perennial ryegrass, and wheat gene families.

The invertase enzyme family is responsible for carbohydrate metabolism in rice, perennial ryegrass, and wheat. Fructan molecules accumulate in cell vacuoles of perennial ryegrass and wheat and are associated with abiotic stress tolerance. High levels of amino acid similarity between the fructosyltransferases responsible for fructan accumulation indicates that they may have evolved from invertase-like ancestral genes. In this study, we have applied comparative genomics to determine the mechanisms that lead to the evolution of fructosytransferase and invertase genes in rice, perennial ryegrass, and wheat. Duplications and rearrangements have been inferred to generate variant forms of the rice invertases since divergence from a common grass progenitor. The occurrence of multiple copies of fructosyltransferase genes indicated that duplication events continued during evolution of the wheat and perennial ryegrass lineages. Further gene rearrangements were evident in perennial ryegrass genes, albeit at a reduced level compared with the rice invertases. Gene orthologs were largely static after duplication during evolution of the wheat lineage. This study details evolutionary events that contribute to fructosyltransferase and invertase gene variation in grasses.     

Complete mitochondrial genomes of Ceratobaeus sp. and Idris sp. (Hymenoptera: Scelionidae): shared gene rearrangements as potential phylogenetic markers at the tribal level

We sequenced the complete mitochondrial genomes of two sceliond taxa (Ceratobaeus sp. and Idris sp.). An atypical tRNA-Arg which lacks a D-stem was identified in both taxa, and represents a potentially derived character of sceliond wasps. A number of tRNA genes have rearranged in the two mitochondrial genomes compared with the ancestral organization. Some of these derived genome organizations are shared, and thus have much potential as phylogenetic markers at the tribal level in the subfamily Scelioninae. We test the influence of third codon inclusion/exclusion, alignment methods and partition schemes on the reconstruction of phylogenetic relationships. The results show that inclusion of third codon positions does not appear to be problematic when investigating the phylogeny of closely related taxa. Muscle and PartitionFinder schemes significantly improve the likelihood scores.     

Topology and maturation of filamentous haemagglutinin suggest a new model for two-partner secretion

Two-partner secretion (TPS) is the most widely distributed secretion pathway known. These systems export large exoproteins through highly conserved channel-forming beta-barrel proteins. Filamentous haemagglutinin (FHA), expressed by Bordetella species, is the prototypical TPS family member. Here we show that the C-terminus of mature FHA, as opposed to the N-terminus as previously proposed, is exposed on the cell surface and is required for mediating adherence to cultured epithelial cells. We show that the C-terminus of the FHA pro-protein (FhaB) is required for FHA function in vitro and in vivo and we show that cleavage of FhaB to form FHA is not the mechanism by which FHA is released from the cell. Based on these data, we propose a new model for TPS. This model provides an explanation for the energetics of export of globular protein domains across membranes in the absence of ATP and it suggests a new mechanism for the control of protein folding.