A gene‐driven recovery mechanism: Drosophila larvae increase feeding activity for post‐stress weight recovery

Recovery from weight loss after stress is important for all organisms, although the recovery mechanisms are not fully understood. We are working to clarify these mechanisms. Here, we recorded enhanced feeding activity of Drosophila melanogaster larvae from 2 to 4 h after heat stress at 35°C for 1 h. During the post‐stress period, expression levels of sweet taste gustatory receptor genes (Grs), Gr5a, Gr43a, Gr64a, and Gr64f, were elevated, whereas bitter taste Grs, Gr66a, and Gr33a, were decreased in expression and expression of a non‐typical taste receptor Gr, Gr68a, was unchanged. Similar upregulation of Gr5a and downregulation of Gr66a was recorded after cold stress at 4°C. Expression levels of tropomyosin and ATP synthase ß subunit were significantly increased in larval mouth parts around 3 to 5 h after the heat stress. We infer that up‐regulation of post‐stress larval feeding activity, and weight recovery, is mediated by increasing capacity for mouth part muscular movements and changes in taste sensing physiology. We propose that Drosophila larvae, and likely insects generally, express an efficient mechanism to recover from weight loss during post‐stress periods.     

A<Emphasis Type="Italic"> Dickkopf</Emphasis>-<Emphasis Type="Italic">3</Emphasis>-related gene is expressed in differentiating nematocytes in the basal metazoan<Emphasis Type="Italic"> Hydra</Emphasis>

In vertebrate development the Dickkopf protein family carries out multiple functions and is represented by at least four different genes with distinct biological activities. In invertebrates such as Drosophila and Caenorhabditis, Dickkopf genes have so far not been identified. Here we describe the identification and characterization of a Dickkopf gene with a deduced amino acid sequence closely related to that of chicken Dkk-3 in the basal metazoan Hydra. HyDkk-3 appears to be the only Dickkopf gene in Hydra. The gene is expressed in the gastric region in nematocytes at a late differentiation stage. In silico searches of EST and genome databases indicated the absence of Dkk genes from the protostomes Drosophila and Caenorhabditis, whereas within the deuterostomes, a Dkk-3 gene could be identified in the genome of the urochordate Ciona intestinalis. The results indicate that at an early stage of evolution of multicellularity Dickkopf proteins have already played important roles as developmental signals. They also suggest that vertebrate Dkk-1, 2 and 4 may have originated from a common ancestor gene of Dkk-3.H. Fedders and R. Augustin contributed equally to this workEdited by D. Tautz     

Association Between Levels of Coding Sequence Divergence and Gene Misregulation in <Emphasis Type="Italic">Drosophila</Emphasis> Male Hybrids

Previous studies have shown widespread conservation of gene expression levels between species of the Drosophila melanogaster subgroup as well as a positive correlation between coding sequence divergence and expression level divergence between species. Meanwhile, large-scale misregulation of gene expression level has been described in interspecific sterile hybrids between D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Using data from gene expression analysis involving D. simulans, D. melanogaster, and their hybrids, we observed a significant positive correlation between protein sequence divergence and gene expression differences between hybrids and their parental species. Furthermore, we demonstrate that underexpressed misregulated genes in hybrids are evolving more rapidly at the protein sequence level than nonmisregulated genes or overexpressed misregulated genes, highlighting the possible effects of sexual and natural selection as male-biased genes and nonessential genes are the principal gene categories affected by interspecific hybrid misregulation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Carlo G. Artieri and Wilfried Haerty contributed equally to this publication.     

Comparative genomic analysis of teleost fish <Emphasis Type="Italic">bmal</Emphasis> genes

Bmal1 (Brain and muscle ARNT like 1) gene is a key circadian clock gene. Tetrapods also have the second Bmal gene, Bmal2. Fruit fly has only one bmal1/cycle gene. Interrogation of the five teleost fish genome sequences coupled with phylogenetic and splice site analyses found that zebrafish have two bmal1 genes, bmal1a and bmal1b, and bmal2a; Japanese pufferfish (fugu), green spotted pufferfish (tetraodon) and Japanese medaka fish each have two bmal2 genes, bmal2a and bmal2b, and bmal1a; and three-spine stickleback have bmal1a and bmal2b. Syntenic analysis further indicated that zebrafish bmal1a/bmal1b, and fugu, tetraodon and medaka bmal2a/bmal2b are ancient duplicates. Although the dN/dS ratios of these four fish bmal duplicates are all <1, implicating they have been under purifying selection, the Tajima relative rate test showed that fugu, tetraodon and medaka bmal2a/bmal2b have asymmetric evolutionary rates, suggesting that one of these duplicates have been subject to positive selection or relaxed functional constraint. These results support the notion that teleost fish bmal genes were derived from the fish-specific genome duplication (FSGD), divergent resolution following the duplication led to retaining different ancient bmal duplicates in different fishes, which could have shaped the evolution of the complex teleost fish timekeeping mechanisms. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evolutionary History and Mode of the <Emphasis Type="Italic">amylase</Emphasis> Multigene Family in <Emphasis Type="Italic">Drosophila</Emphasis>

Previous studies indicate that the tandemly repeated members of the amylase (Amy) gene family evolved in a concerted manner in the melanogaster subgroup and in some other species. In this paper, we analyzed all of the 49 active and complete Amy gene sequences in Drosophila, mostly from subgenus Sophophora. Phylogenetic analysis indicated that the two types of diverged Amy genes in the Drosophila montium subgroup and Drosophila ananassae, which are located in distant chromosomal regions from each other, originated independently in different evolutionary lineages of the melanogaster group after the split of the obscura and melanogaster groups. One of the two clusters was lost after duplication in the melanogaster subgroup. Given the time, 24.9 mya, of divergence between the obscura and the melanogaster groups (Russo et al. 1995), the two duplication events were estimated to occur at about 13.96 ± 1.93 and 12.38 ± 1.76 mya in the montium subgroup and D. ananassae, respectively. An accelerated rate of amino acid changes was not observed in either lineage after these gene duplications. However, the G+C contents at the third codon positions (GC3) decreased significantly along one of the two Amy clusters both in the montium subgroup and in D. ananassae right after gene duplication. Furthermore, one of the two types of the Amy genes with a lower GC3 content has lost a specific regulatory element within the montium subgroup species and D. ananassae. While the tandemly repeated members evolved in a concerted manner, the two types of diverged Amy genes in Drosophila experienced frequent gene duplication, gene loss, and divergent evolution following the model of a birth-and-death process.     

<Emphasis Type="Italic">Drosophila</Emphasis><Emphasis Type="Italic">P</Emphasis> transposons of the urochordata <Emphasis Type="Italic">Ciona intestinalis</Emphasis>

P transposons belong to the eukaryotic DNA transposons, which are transposed by a cut and paste mechanism using a P-element-coded transposase. They have been detected in Drosophila, and reside as single copies and stable homologous sequences in many vertebrate species. We present the P elements Pcin1, Pcin2 and Pcin3 from Ciona intestinalis, a species of the most primitive chordates, and compare them with those from Ciona savignyi. They showed typical DNA transposon structures, namely terminal inverted repeats and target site duplications. The coding region of Pcin1 consisted of 13 small exons that could be translated into a P-transposon-homologous protein. C. intestinalis and C. savignyi displayed nearly the same phenotype. However, their P elements were highly divergent and the assumed P transposase from C. intestinalis was more closely related to the transposase from Drosophila melanogaster than to the transposase of C. savignyi. The present study showed that P elements with typical features of transposable DNA elements may be found already at the base of the chordate lineage. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Molecular evolution of <Emphasis Type="Italic">PKD2</Emphasis> gene family in mammals

PKD2 gene encodes a critical cation channel protein that plays important roles in various developmental processes and is usually evolutionarily conserved. In the present study, we analyzed the evolutionary patterns of PKD2 and its homologous genes (PKD2L1, PKD2L2) from nine mammalian species. In this study, we demonstrated the orthologs of PKD2 gene family evolved under a dominant purifying selection force. Our results in combination with the reported evidences from functional researches suggested the entire PKD2 gene family are conserved and perform essential biological roles during mammalian evolution. In rodents, PKD2 gene family members appeared to have evolved more rapidly than other mammalian lineages, probably resulting from relaxation of purifying selection. However, positive selection imposed on synonymous sites also potentially contributed to this case. For the paralogs, our results implied that PKD2L2 genes evolved under a weaker purifying selection constraint than PKD2 and PKD2L1 genes. Interestingly, some loop regions of transmembrane domain of PKD2L2 exhibited higher P _N/P _S ratios than expected, suggesting these regions are more functional divergent in organisms and worthy of special attention. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Chun Ye, Huan Sun have contributed equally to this work.     

Natural selection in the <Emphasis Type="Italic">TLR-</Emphasis>related genes in the course of primate evolution

The innate immune system constitutes the front line of host defense against pathogens. Toll-like receptors (TLRs) recognize molecules derived from pathogens and play crucial roles in the innate immune system. Here, we provide evidence that the TLR-related genes have come under natural selection pressure in the course of primate evolution. We compared the nucleotide sequences of 16 TLR-related genes, including TLRs (TLR1–10), MYD88, TILAP, TICAM1, TICAM2, MD2, and CD14, among seven primate species. Analysis of the non-synonymous/synonymous substitution ratio revealed the presence of both strictly conserved and rapidly evolving regions in the TLR-related genes. The genomic segments encoding the intracellular Toll/interleukin 1 receptor domains, which exhibited lower rates of non-synonymous substitution, have undergone purifying selection. In contrast, TLR4, which carried a high proportion of non-synonymous substitutions in the part of extracellular domain spanning 200 amino acids, was found to have been the suggestive target of positive Darwinian selection in primate evolution. However, sequence analyses from 25 primate species, including eight hominoids, six Old World monkeys, eight New World monkeys, and three prosimians, showed no evidence that the pressure of positive Darwinian selection has shaped the pattern of sequence variations in TLR4 among New World monkeys and prosimians. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

TGFβ signaling in Tribolium: vertebrate-like components in a beetle

The cytokines of the TGFβ superfamily are highly conserved in evolution and elicit a diverse range of cellular responses in all metazoa. In Drosophila, the signaling pathways of the two TGFβ subfamilies, Activins and Bone Morphogenetic Proteins (BMPs), have been well studied. To address the question of whether the findings from Drosophila are representative of insects in general, we analyzed the components of TGFβ-signaling present in the genome of the beetle Tribolium castaneum. We were able to identify orthologs of the BMPs Decapentaplegic and Glass bottom boat, of the Activins Activinβ and Dawdle, as well as orthologs of the less well-known ligands Myoglianin and Maverick, together with orthologs of all TGFβ receptors and cytoplasmic signal transducers present in Drosophila. This indicates that the diversity of TGFβ signaling components is generally well conserved between Drosophila and Tribolium. However, the genome of the beetle—and of the bee Apis mellifera—lacks an ortholog of the Drosophila BMP Screw but does contain a vertebrate-like BMP10 homolog which is not found in Drosophila. Concerning BMP inhibitors, Tribolium displays an even more vertebrate-like ensemble of components. We found two orthologs of the vertebrate DAN family, Dan and Gremlin, and show embryonic expression of a vertebrate-like BAMBI ortholog, all of which are absent in Drosophila. This suggests that Tribolium might have retained a more ancestral composition of TGFβ signaling components and that TGFβ signaling underwent considerable change in the Drosophila lineage. Tribolium is an excellent model to study the function of these ancestral signaling components in insects. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

<Emphasis Type="Italic">odd-skipped</Emphasis> homologs function during gut development in<Emphasis Type="Italic"> C. elegans</Emphasis>

Genes in the odd-skipped (odd) family encode a discrete subset of C2H2 zinc finger proteins that are widely distributed among metazoan phyla. Although the initial member (odd) was identified as a Drosophila pair-rule gene, various homologs are expressed within each of the three germ layers in complex patterns that suggest roles in many pathways beyond segmentation. To further investigate the evolutionary history and extant functions of genes in this family, we have initiated a characterization of two homologs, odd-1 and odd-2, identified in the genome of the nematode, Caenorhabditis elegans. Sequence comparisons with homologs from insects (Drosophila and Anopheles) and mammals suggest that two paralogs were present within an ancestral metazoan; additional insect paralogs and both extant mammalian genes likely resulted from gene duplications that occurred after the split between the arthropods and chordates. Analyses of gene function using RNAi indicate that odd-1 and odd-2 play essential and distinct roles during gut development. Specific expression of both genes in the developing intestine and other cells in the vicinity of the gut was shown using GFP-reporters. These results indicate primary functions for both genes that are most like those of the Drosophila paralogs bowel and drumstick, and support a model in which gut specification represents the ancestral role for genes in this family.Edited by C. Desplan     

<Emphasis Type="Italic">hedgehog</Emphasis> is a segment polarity gene in a crustacean and a chelicerate

The evolution of arthropod segmentation has been studied by comparing expression patterns of pair-rule and segment polarity genes in various species. In Drosophila, the formation and maintenance of the parasegmental boundaries depend on the interactions between the wingless (wg), engrailed (en) and hedgehog (hh) genes. Until now, the expression pattern of hh has not been analysed to such a great extent as en or wg. We report the cloning and expression analysis of hh genes from Euscorpius flavicaudis, a chelicerate, and Artemia franciscana, a branchiopod crustacean. Our data provide evidence that hh, being expressed in the posterior part of every segment, is a segment polarity gene in both organisms. Additional hh expression sites were observed in the rostrum and appendages of Euscorpius and in the gut of Artemia. From the available data on hh expression in various bilaterians, we review the various hypotheses on the evolution of hh function and we suggest an ancestral role of hh in proctodeum specification and gut formation.Edited by D. Tautz     

Duplicate Gene Evolution Toward Multiple Fates at the <Emphasis Type="Italic">Drosophila melanogaster HIP</Emphasis>/<Emphasis Type="Italic">HIP-Replacement</Emphasis> Locus

Hsc/Hsp70-interacting protein (HIP) is a rapidly evolving Hsp70 cofactor. Analyses of multiple Drosophila species indicate that the HIP gene is duplicated only in D. melanogaster. The HIP region, in fact, contains seven distinctly evolving duplicated genes. The regional duplication occurred in two steps, fixed rapidly, and illustrates multiple modes of duplicate gene evolution. HIP and its duplicate HIP-R are adaptively evolving in a manner unique to the region: they exhibit elevated divergence from other drosophilids and low polymorphism within D. melanogaster. HIP and HIP-R are virtually identical, share polymorphisms, and are subject to gene conversion. In contrast, two other duplicate genes in the region, CG33221 and GP-CG32779, are pseudogenes, and the chimeric gene Crg1 is subject to balancing selection. HIP and HIP-R are evolving rapidly and adaptively; however, positive selection is not sufficient to explain the molecular evolution of the region as a whole. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Gene Conversion and Positive Selection Driving the Evolution of the <Emphasis Type="Italic">Caenorhabditis</Emphasis> ssp. ZIM/HIM-8 Protein Family

In C. elegans, four C2H2 zinc-finger proteins (ZIM-1, ZIM-2, ZIM-3, and HIM-8), which are arranged in tandem, mediate chromosome-specific pairing and synapsis during meiosis. The zim/him-8 genes from three Caenorhabditis species were contrasted in an effort to investigate the mechanisms driving their evolution. Here it is shown that the preservation of higher degree of sequence similarity in the N-terminal portion, particularly in several regions within the second exon between paralogous zim genes (especially between zim-1 and zim-3), is due to independent interparalogue gene conversions. However, the evolutionary force is not uniformly strong across species. The present data reveal that more frequent gene conversion events have occurred in C. elegans, whereas only gene conversions between zim-1 and zim-3 are detected in C. remanei. Although gene conversions are predicted to be present among zim-1, zim-2, and zim-3 in C. briggsae, the conversion tracts between zim-1/zim-2 and zim-2/zim-3 are very short. Moreover, positive selection analysis was performed on the basis of the significantly discordant phylogenies reconstructed using the N- and C-terminal sequences, respectively. Several codon sites located in the regions that are supposed not to have experienced gene conversions are predicted to be under the influence of positive selection. In comparison, stronger positive selection has acted on the C-terminal region relative to the N-terminal region. Thus, the zim/him-8 genes that evolve concertedly have also been shown to undergo adaptive diversifying selection. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Correlated Evolution of Synonymous and Nonsynonymous Sites in <Emphasis Type="Italic">Drosophila</Emphasis>

Recent work has shown that Drosophila melanogaster genes with fast-evolving nonsynonymous sites have lower codon usage bias. This pattern has been attributed to interference between positive selection at nonsynonymous sites and weak selection on codon usage. Here we have looked for this correlation in a much larger and less biased dataset, comprising 630 gene pairs from D. melanogaster and D. yakuba. We confirmed that there is a negative correlation between the rate of nonsynonymous substitutions (d_N) and codon bias in D. melanogaster. We then tested the interference hypothesis and other alternative explanations, including one involving gene expression. We found that d_N indeed correlates with the level of gene expression. Given that gene expression is a strong determinant of codon bias, the relationship between d_N and codon bias might be a by-product of gene expression. However, our tests show that none of the hypotheses we consider seem to explain the data fully.This article contains online supplementary material.Reviewing Editor: Dr. John Huelsenbeck     

Rapid Intraspecific Evolution of miRNA and siRNA Genes in the Mosquito Aedes aegypti

RNA silencing, or RNA interference (RNAi) in metazoans mediates development, reduces viral infection and limits transposon mobility. RNA silencing involves 21–30 nucleotide RNAs classified into microRNA (miRNA), exogenous and endogenous small interfering RNAs (siRNA), and Piwi-interacting RNA (piRNA). Knock-out, silencing and mutagenesis of genes in the exogenous siRNA (exo-siRNA) regulatory network demonstrate the importance of this RNAi pathway in antiviral immunity in Drosophila and mosquitoes. In Drosophila, genes encoding components for processing exo-siRNAs are among the fastest evolving 3% of all genes, suggesting that infection with pathogenic RNA viruses may drive diversifying selection in their host. In contrast, paralogous miRNA pathway genes do not evolve more rapidly than the genome average. Silencing of exo-siRNA pathway genes in mosquitoes orally infected with arboviruses leads to increased viral replication, but little is known about the comparative patterns of molecular evolution among the exo-siRNA and miRNA pathways genes in mosquitoes. We generated nearly complete sequences of all exons of major miRNA and siRNA pathway genes dicer-1 and dicer-2, argonaute-1 and argonaute-2, and r3d1 and r2d2 in 104 Aedes aegypti mosquitoes collected from six distinct geographic populations and analyzed their genetic diversity. The ratio of replacement to silent amino acid substitutions was 1.4 fold higher in dicer-2 than in dicer-1, 27.4 fold higher in argonaute-2 than in argonaute-1 and similar in r2d2 and r3d1. Positive selection was supported in 32% of non-synonymous sites in dicer-1, in 47% of sites in dicer-2, in 30% of sites in argonaute-1, in all sites in argonaute-2, in 22% of sites in r3d1 and in 55% of sites in r2d2. Unlike Drosophila, in Ae. aegypti, both exo-siRNA and miRNA pathway genes appear to be undergoing rapid, positive, diversifying selection. Furthermore, refractoriness of mosquitoes to infection with dengue virus was significantly positively correlated for nucleotide diversity indices in dicer-2.     

The involvement of <Emphasis Type="Italic">engrailed</Emphasis> and <Emphasis Type="Italic">wingless</Emphasis> during segmentation in the onychophoran <Emphasis Type="Italic">Euperipatoides kanangrensis</Emphasis> (Peripatopsidae: Onychophora) (Reid 1996)

As the putative sister group to the arthropods, onychophorans can provide insight into ancestral developmental mechanisms in the panarthropod clade. Here, we examine the expression during segmentation of orthologues of wingless (Wnt1) and engrailed, two genes that play a key role in defining segment boundaries in Drosophila and that appear to play a role in segmentation in many other arthropods. Both are expressed in segmentally reiterated stripes in all forming segments except the first (brain) segment, which only shows an engrailed stripe. Engrailed is expressed before segments are morphologically visible and is expressed in both mesoderm and ectoderm. Segmental wingless expression is not detectable until after mesodermal somites are clearly distinct. Early engrailed expression lies in and extends to both sides of the furrow that first demarcates segments in the ectoderm, but is largely restricted to the posterior part of somites. Wingless expression lies immediately anterior to engrailed expression, as it does in many arthropods, but there is no precise cellular boundary between the two expression domains analogous to the overt parasegment boundary seen in Drosophila. Engrailed stripes extend along the posterior part of each limb bud, including the antenna, while wingless is restricted to the distal tip of the limbs and the neurectoderm basal to the limbs. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Tempo and Mode in the Endocannaboinoid System

The best-known endocannabinoid ligands, anandamide and 2-AG, signal at least seven receptors and involve ten metabolic enzymes. Genes for the receptors and enzymes were examined for heterogeneities in tempo (relative rate of evolution, RRE) and mode (selection pressure, Ka/Ks) in six organisms with sequenced genomes. BLAST identified orthologs as reciprocal best hits, and nucleotide alignments were performed with ClustalX and MacClade. Two bioinformatics platforms, LiKaKs (a distance-based LWL85 model) and SNAP (a parsimony-based NG86 model) made pairwise comparisons of orthologs in murids (rat and mouse) and primates (human and macaque). Mean RRE of the 18 endocannabinoid genes was significantly greater in murids than primates, whereas mean Ka/Ks did not differ significantly. Next we used FUGE (tree-based maximum-likelihood model) to compute human lineage-specific Ka/Ks calculations for 18 genes, which ranged from 1.11 to 0.00, in rank order from highest to lowest: PTPN22, NAAA, TRPV1, TRPA1, NAPE-PLD, MAGL, PPARγ, FAAH1, COX2, FAAH2, ABDH4, CB2, GPR55, DAGLβ, PPARα, TRPV4, CB1, DAGLα; differences were significant (p < 0.0001). Rat and mouse presented different rank orders (e.g., GPR55 generated the greatest Ka/Ks ratio). The 18 genes were then tested for recent positive selection (within 10,000 yr) using an extended haplotype homozygosity analysis of SNP data from the HapMap database. Significant evidence (p < 0.05) of a positive “selective sweep” was exhibited by PTPN22, TRPV1, NAPE-PLD, and DAGLα. In conclusion, the endocannabinoid system is collectively under strong purifying selection, although some genes show evidence of adaptive evolution. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users. Reviewing Editor: Dr. Bryan Fry     

Comparative genomic analysis of the major histocompatibility complex class I region in the teleost genus <Emphasis Type="Italic">Oryzias</Emphasis>

The major histocompatibility complex (MHC) class I region of teleosts harbors a tight cluster of the class IA genes and several other genes directly involved in class I antigen presentation. Moreover, the dichotomous haplotypic lineages (termed d- and N- lineages) of the proteasome subunit beta genes, PSMB8 and PSMB10, are present in this region of the medaka, Oryzias latipes. To understand the evolution of the Oryzias MHC class I region at the nucleotide sequence level, we analyzed bacterial artificial chromosome clones covering the MHC class I region containing the d- lineage of Oryzias luzonensis and the d- and N- lineages of Oryzias dancena. Comparison among these three elucidated sequences and the published sequences of the d- and N- lineages of O. latipes indicated that the order and orientation of the encoded genes were completely conserved among these five genomic regions, except for the class IA genes, which showed species-specific variation in copy number. The PSMB8 and PSMB10 genes showed trans-species dimorphism. The remaining regions flanking the PSMB10, PSMB8, and class IA genes showed high degrees of sequence conservation at interspecies as well as intraspecies levels. Thus, the three independent evolutionary patterns under apparently distinctive selective pressures are recognized in the Oryzias MHC class I region. Electronic Supplementary Material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The chloroplast genome from a lycophyte (microphyllophyte), <Emphasis Type="Italic">Selaginella uncinata</Emphasis>, has a unique inversion,transpositions and many gene losses

We determined the complete nucleotide sequence of the chloroplast genome of Selaginella uncinata, a lycophyte belonging to the basal lineage of the vascular plants. The circular double-stranded DNA is 144,170 bp, with an inverted repeat of 25,578 bp separated by a large single copy region (LSC) of 77,706 bp and a small single copy region (SSC) of 40,886 bp. We assigned 81 protein-coding genes including four pseudogenes, four rRNA genes and only 12 tRNA genes. Four genes, rps15, rps16, rpl32 and ycf10, found in most chloroplast genomes in land plants were not present in S. uncinata. While gene order and arrangement of the chloroplast genome of another lycophyte, Hupertzia lucidula, are almost the same as those of bryophytes, those of S. uncinata differ considerably from the typical structure of bryophytes with respect to the presence of a unique 20 kb inversion within the LSC, transposition of two segments from the LSC to the SSC and many gene losses. Thus, the organization of the S. uncinata chloroplast genome provides a new insight into the evolution of lycophytes, which were separated from euphyllophytes approximately 400 million years ago. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.