Genome-wide analysis of restriction-modification system in unicellular and filamentous cyanobacteria

Cyanobacteria are an ancient group of gram-negative bacteria with strong genome size variation ranging from 1.6 to 9.1 Mb. Here, we first retrieved all the putative restriction-modification (RM) genes in the draft genome of Spirulina and then performed a range of comparative and bioinformatic analyses on RM genes from unicellular and filamentous cyanobacterial genomes. We have identified 6 gene clusters containing putative Type I RMs and 11 putative Type II RMs or the solitary methyltransferases (MTases). RT-PCR analysis reveals that 6 of 18 MTases are not expressed in Spirulina, whereas one hsdM gene, with a mutated cognate hsdS, was detected to be expressed. Our results indicate that the number of RM genes in filamentous cyanobacteria is significantly higher than in unicellular species, and this expansion of RM systems in filamentous cyanobacteria may be related to their wide range of ecological tolerance. Furthermore, a coevolutionary pattern is found between hsdM and hsdR, with a large number of site pairs positively or negatively correlated, indicating the functional importance of these pairing interactions between their tertiary structures. No evidence for positive selection is found for the majority of RMs, e. g., hsdM, hsdS, hsdR, and Type II restriction endonuclease gene families, while a group of MTases exhibit a remarkable signature of adaptive evolution. Sites and genes identified here to have been under positive selection would provide targets for further research on their structural and functional evaluations.

The complete mitochondrial genome of the large yellow croaker, Larimichthys crocea (Perciformes, Sciaenidae): Unusual features of its control region and the phylogenetic position of the Sciaenidae

To understand the systematic status of Larimichthys crocea in the Percoidei, we determined the complete mitochondrial (mt) genome sequence using 454 sequencing-by-synthesis technology. The complete mt genome is 16,466 bp in length including the typical structure of 22 tRNAs, 2 rRNAs, 13 protein-coding genes and the noncoding control region (CR). Further sequencing for the complete CR was performed using the primers Cyt b-F and 12S-R on six L crocea individuals and two L polyactis individuals. Interestingly, all seven CR sequences from L crocea were identical while the three sequences from L polyactis were distinct (including one from GenBank). Although the conserved blocks such as TAS and CSB-1, -2, and -3 are readily identifiable in the control regions of the two species, the typical central conserved blocks CSB-D, -E, and -F could not be detected, while they are found in Cynoscion acoupa of Sciaenidae and other Percoidei species. Phylogenetic analysis shows that L crocea is a relatively recently emerged species in Sciaenidae and this family is closely related to family Pomacanthidae within the Percoidei. L crocea, as the first species of Sciaenidae with complete mitochondrial genome available, will provide important information on the molecular evolution of the group. Moreover, the genus-specific pair of primers designed in this study for amplifying the complete mt control region will be very useful in studies on the population genetics and conservation biology of Larimichthys. (c) 2008 Elsevier B.V. All rights reserved.

The complete mitochondrial genome sequence of the cutlassfish Trichiurus japonicus (Perciformes: Trichiuridae) : Genome characterization and phylogenetic considerations

Mitochondrial genome sequence and structure analysis has become a powerful tool for studying molecular evolution and phylogenetic relationships. To understand the systematic status of Trichiurus japonicus in suborder Scombroidei, we determined the complete mitochondrial genome (mitogenome) sequence using the long-polymerase chain reaction (long-PCR) and shotgun sequencing method. The entire mitogenome is 16,796 by in length and has three unusual features, including (1) the absence of tRNA(Pro) gene, (2) the possibly nonfunctional light-strand replication origin (O-L) showing a shorter loop in secondary structure and no conserved motif (5'-GCCGG-3'), (3) two sets of the tandem repeats at the 5' and 3' ends of the control region. The three features seem common for Trichiurus mitogenomes, as we have confirmed them in other three T. japonicus individuals and in T nanhaiensis. Phylogenetic analysis does not support the monophyly of Trichiuridae, which is against the morphological result. T. japonicus is most closely related to those species of family Scombridae; they in turn have a sister relationship with Perciformes members including suborders Acanthuroidei, Caproidei, Notothenioidei, Zoarcoidei, Trachinoidei, and some species of Labroidei, based on the current dataset of complete mitogenome. T japonicus together with T. brevis, T lepturus and Aphanopus carbo form a clade distinct from Lepidopus caudatus in terms of the complete Cyt b sequences. T. japonicus mitogenome, as the first discovered complete mitogenome of Trichiuridae, should provide important information on both genomics and phylogenetics of Trichiuridae. (C) 2009 Elsevier B.V. All rights reserved.

Complete mitochondrial genome of the sea cucumber Apostichopus japonicus (Echinodermata: Holothuroidea): The first representative from the subclass Aspidochirotacea with the echinoderm ground pattern

Complete mitochondrial genome plays an important role in the accurate revelation of phylogenetic relationships among metazoans. Here we present the complete mitochondrial genome sequence from a sea cucumber Apostichopus japonicus (Echinodermata: Holothuroidea), which is the first representative from the subclass Aspidochirotacea. The mitochondrial genome of A. japonicus is 16,096 bp in length. The heavy strand consists of 31.8% A, 20.2% C, 17.9% G, and 30.1% T bases (AT skew = 0.027: GC skew = 0.062). It contains thirteen protein-coding genes (PCGs), twenty-two transfer RNA genes, and two ribosomal RNA genes. There are a total of 3793 codons in all thirteen mitochondrial PCGs, excluding incomplete termination codons. The most frequently used amino acid is Leu (15.77%), followed by Set (9.73%), Met (8.62%), Phe (7.94%), and Ala (7.28%). Intergenetic regions in the mitochondrial genome of A. japonicus are 839 bp in total, with three relatively large regions of Unassigned Sequences (UAS) greater than 100 bp. The gene order of A. japonicus is identical to that observed in the five studied sea urchins, which confirms that the gene order shared by the two classes (Holothuroidea and Echinoidea) is a ground pattern of echinoderm mitochondrial genomes. Bayesian tree based on the cob gene supports the following relationship: (outgroup, (Crinoids, (Asteroids, Ophiuroids, (Echinoids, Holothuroids)))). (C) 2009 Elsevier B.V. All rights reserved.

Analysis of tandem repeats in the genome of Chinese shrimp Fenneropenaeus chinensis

Through random sequencing, we found a total of 884000 base-pairs (bp) of random genomic sequences in the genome of Chinese shrimp (Fenneropenaeus chinensis). Using bio-soft Tandem Repeat Finder (TRF) software, 2159 tandem repeats were found, in which there were 1714 microsatellites and 445 minisatellites, accounting for 79.4% and 20.6% of repeat sequences, respectively. The cumulative length of repeat sequences was found to be 116685 bp, accounting for 13.2% of the total DNA sequence; the cumulative length of microsatellites occupied 9.78% of the total DNA sequence, and that of minisatellites occupied 3.42%. In decreasing order, the 20 most abundant repeat sequence classes were as follows: AT (557), AC (471), AG (274), AAT (92), A (56), AAG (28), ATC (27), ATAG (27), AGG (18), ACT (15), C (11), AAC (11), ACAT (11), CAGA (10), AGAA (9), AGGG (7), CAAA (7), CGCA (6), ATAA (6), AGAGAA (6). Dinucleotide repeats, not only in the aspect of the number, but also in cumulative length, were the preponderant repeat type. There were few classes and low copy numbers of repeat units of the pentanucleotide repeat type, which included only three classes: AGAGA, GAGGC and AAAGA. The classes and copy numbers of heptanucleotide, eleven-nucleotide and thirteen-nucleotide primer-number-composed repeats were distinctly less than that of repeat types beside them.

Complete mitochondrial genome of the Asian paddle crab Charybdis japonica (Crustacea: Decapoda: Portunidae): gene rearrangement of the marine brachyurans and phylogenetic considerations of the decapods

Given the commercial and ecological importance of the Asian paddle crab, Charybdis japonica, there is a clearly need for genetic and molecular research on this species. Here, we present the complete mitochondrial genome sequence of C. japonica, determined by the long-polymerase chain reaction and primer walking sequencing method. The entire genome is 15,738 bp in length, encoding a standard set of 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes, plus the putative control region, which is typical for metazoans. The total A+T content of the genome is 69.2%, lower than the other brachyuran crabs except for Callinectes sapidus. The gene order is identical to the published marine brachyurans and differs from the ancestral pancrustacean order by only the position of the tRNA (His) gene. Phylogenetic analyses using the concatenated nucleotide and amino acid sequences of 13 protein-coding genes strongly support the monophyly of Dendrobranchiata and Pleocyemata, which is consistent with the previous taxonomic classification. However, the systematic status of Charybdis within subfamily Thalamitinae of family Portunidae is not supported. C. japonica, as the first species of Charybdis with complete mitochondrial genome available, will provide important information on both genomics and molecular ecology of the group.

The complete mitochondrial genome of the mantid shrimp Oratosquilla oratoria (Crustacea: Malacostraca: Stomatopoda): Novel non-coding regions features and phylogenetic implications of the Stomatopoda

The complete mitochondrial (mt) genome sequence of Oratosquilla oratoria (Crustacea: Malacostraca: Stomatopoda) was determined; a circular molecule of 15,783 bp in length. The gene content and arrangement are consistent with the pancrustacean ground pattern. The mt control region of O. oratoria is characterized by no GA-block near the 3' end and different position of [TA(A)]n-blocks compared with other reported Stomatopoda species. The sequence of the second hairpin structure is relative conserved which suggests this region may be a synapomorphic character for the Stomatopoda. In addition, a relative large intergenic spacer (101 bp) with higher A + T content than that in control region was identified between the tRNA(Glu) and tRNA(Phe) genes. Phylogenetic analyses based on the current dataset of complete mt genomes strongly support the Stomatopoda is closely related to Euphausiacea. They in turn cluster with Penaeoidea and Caridea clades while other decapods form a separate group, which rejects the monophyly of Decapoda. This challenges the suitability of Stomatopoda as an outgroup of Decapoda in phylogenetic analyses. The basal position of Stomatopoda within Eumalacostraca according to the morphological characters is also questioned. (C) 2010 Elsevier Inc. All rights reserved.

A close phylogenetic relationship between Sipuncula and Annelida evidenced from the complete mitochondrial genome sequence of Phascolosoma esculenta

Background: There are many advantages to the application of complete mitochondrial (mt) genomes in the accurate reconstruction of phylogenetic relationships in Metazoa. Although over one thousand metazoan genomes have been sequenced, the taxonomic sampling is highly biased, left with many phyla without a single representative of complete mitochondrial genome. Sipuncula (peanut worms or star worms) is a small taxon of worm-like marine organisms with an uncertain phylogenetic position. In this report, we present the mitochondrial genome sequence of Phascolosoma esculenta, the first complete mitochondrial genome of the phylum. Results: The mitochondrial genome of P. esculenta is 15,494 bp in length. The coding strand consists of 32.1% A, 21.5% C, 13.0% G, and 33.4% T bases (AT = 65.5%; AT skew = -0.019; GC skew = -0.248). It contains thirteen protein-coding genes (PCGs) with 3,709 codons in total, twenty-two transfer RNA genes, two ribosomal RNA genes and a non-coding AT-rich region (AT = 74.2%). All of the 37 identified genes are transcribed from the same DNA strand. Compared with the typical set of metazoan mt genomes, sipunculid lacks trnR but has an additional trnM. Maximum Likelihood and Bayesian analyses of the protein sequences show that Myzostomida, Sipuncula and Annelida (including echiurans and pogonophorans) form a monophyletic group, which supports a closer relationship between Sipuncula and Annelida than with Mollusca, Brachiopoda, and some other lophotrochozoan groups. Conclusion: This is the first report of a complete mitochondrial genome as a representative within the phylum Sipuncula. It shares many more similar features with the four known annelid and one echiuran mtDNAs. Firstly, sipunculans and annelids share quite similar gene order in the mitochondrial genome, with all 37 genes located on the same strand; secondly, phylogenetic analyses based on the concatenated protein sequences also strongly support the sipunculan + annelid clade (including echiurans and pogonophorans). Hence annelid "key-characters" including segmentation may be more labile than previously assumed.

The complete mitochondrial genome of the ridgetail white prawn Exopalaemon carinicauda Holthuis, 1950 (Crustacean: Decapoda: Palaemonidae) revealed a novel rearrangement of tRNA genes

The complete mitochondrial (mt) DNA sequence was determined for a ridgetail white prawn, Exopalaemon carinicauda Holthuis, 1950 (Crustacea: Decopoda: Palaemonidae). The mt genome is 15,730 bp in length, encoding a standard set of 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes, which is typical for metazoans. The majority-strand consists of 33.6% A, 23.0% C, 13.4% G, and 30.0% T bases (AT skew = 0.057: GC skew = -0.264). A total of 1045 bp of non-coding nucleotides were observed in 16 intergenic regions,,including a major A+ T rich (79.7%) noncoding region (886 bp). A novel translocation of tRNA(Pro) and tRNA(Thr) was found when comparing this genome with the pancrustacean ground pattern indicating that gene order is not conserved among caridean mitochondria. Furthermore, the rate of Ka/Ks in 13 protein-coding genes between three caridean species is Much less than 1, which indicates a strong Purifying selection within this group. To investigate the phylogenetic relationship within Malacostraca, phylogenetic trees based oil Currently available malacostracan complete mitochondrial sequences were built with the maximum likelihood and Bayesian models. All analyses based oil nucleotide and amino acid data strongly support the monophyly of Decapoda. The Penaeidae, Reptantia, Caridea, and Meiura clades were also recovered as monophyletic groups with Strong Statistical Support. However, the phylogenetic relationships within Pleocyemata are unstable, as represented by the inclusion or exclusion of Caridea. (C) 2009 Elsevier B.V. All rights reserved.

Isolation and characterization of a novel variant of HMW glutenin subunit gene from the St genome of Pseudoroegneria stipifolia

The x- and y-type high molecular weight (HMW) glutenin subunits are conserved seed storage proteins in wheat and related species. Here we describe investigations on the HMW glutenin subunits from several Pseudoroegneria accessions. The electrophoretic mobilities of the HMW glutenin subunits from Pd. stipifolia, Pd tauri and Pd strigosa were much faster than those of orthologous wheat subunits, indicating that their protein size may be smaller than that of wheat subunits. The coding sequence of the Glu-1St1 subunit (encoded by the Pseudoroegneria stipifolia accession PI325181) was isolated, and found to represent the native open reading frame (ORF) by in vitro expression. The deduced amino acid sequence of Glu-1St1 matched with that determined from the native subunit by mass spectrometric analysis. The domain organization in Glu-1St1 showed high similarity with that of typical HMW glutenin subunits. However, Glu-1St1 exhibited several distinct characteristics. First, the length of its repetitive domain was substantially smaller than that of conventional subunits, which explains its much faster electrophoretic mobility in SDS-PAGE. Second, although the N-terminal domain of Glu-1St1 resembled that of y-type subunit, its C-terminal domain was more similar to that of x-type subunit. Third, the N- and C-terminat domains of Glu-1St1 shared conserved features with those of barley D-hordein, but the repeat motifs and the organization of its repetitive domain were more similar to those of HMW glutenin subunits than to D-hordein. We conclude that Glu-1St1 is a novel variant of HMW glutenin subunits. The analysis of Glu-1St1 may provide new insight into the evolution of HMW glutenin subunits in Triticeae species. (C) 2007 Elsevier Ltd. All rights reserved.

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.

Mitochondrial and chloroplast phylogeography of Picea crassifolia Kom. (Pinaceae) in the Qinghai-Tibetan Plateau and adjacent highlands

The disjunct distribution of forests in the Qinghai-Tibetan Plateau (QTP) and adjacent Helan Shan and Daqing Shan highlands provides an excellent model to examine vegetation shifts, glacial refugia and gene flow of key species in this complex landscape region in response to past climatic oscillations and human disturbance. In this study, we examined maternally inherited mitochondrial DNA (nad1 intron b/c and nad5 intron 1) and paternally inherited chloroplast DNA (trnC-trnD) sequence variation within a dominant forest species, Picea crassifolia Kom. We recovered nine mitotypes and two chlorotypes in a survey of 442 individuals from 32 populations sampled throughout the species' range. Significant mitochondrial DNA population subdivision was detected (G(ST) = 0.512; N-ST = 0.679), suggesting low levels of recurrent gene flow through seeds among populations and significant phylogeographical structure (N-ST > GST, P < 0.05). Plateau haplotypes differed in sequence from those in the adjacent highlands, suggesting a long period of allopatric fragmentation between the species in the two regions and the presence of independent refugia in each region during Quaternary glaciations. On the QTP platform, all but one of the disjunct populations surveyed were fixed for the same mitotype, while most populations at the plateau edge contained more than one haplotype with the mitotype that was fixed in plateau platform populations always present at high frequency. This distribution pattern suggests that present-day disjunct populations on the QTP platform experienced a common recolonization history. The same phylogeographical pattern, however, was not detected for paternally inherited chloroplast DNA haplotypes. Two chlorotypes were distributed throughout the range of the species with little geographical population differentiation (G(ST) = N-ST = 0.093). This provides evidence for highly efficient pollen-mediated gene flow among isolated forest patches, both within and between the QTP and adjacent highland populations. A lack of isolation to pollen-mediated gene flow between forests on the QTP and adjacent highlands is surprising given that the Tengger Desert has been a geographical barrier between these two regions for approximately the last 1.8 million years.

Genome constitutions of Hystrix patula, H-duthiei ssp duthiei and H-duthiei ssp longearistata (Poaceae : Triticeae) revealed by meiotic pairing behavior and genomic in-situ hybridization

Genomic constitutions of three taxa of Hystrix Moench, H. patula, H. duthiei ssp. duthiei and H. duthiei ssp. longearistata, were examined by meiotic pairing behavior and genomic in-situ hybridization. Meiotic pairing in hybrids of H. patula x Pseudoroegneria spicata (St), H. patula x Elymus wawawaiensis (StH), H. patula x H. duthiei ssp. longearistata, H. patula x Psathyrostachys huashanica (Ns(h)), H. duthiei ssp. duthiei x Psa. huashanica, H. duthiei ssp. longearistata x Psa. huashanica, Leymus multicaulis (NsXm) x H. duthiei ssp. longearistata averaged 6.53, 12.83, 1.32, 0.29, 5.18, 5.11 and 10.47 bivalents per cell, respectively. The results indicate that H. patula has the StH genome and H. duthiei ssp. duthiei and H. duthiei ssp. longearistata have the NsXm genome. Results of genomic in-situ hybridization analysis strongly supported the chromosome pairing data; therefore it is concluded that the type species of Hystrix, H. patula, should be included in Elymus, and that H. duthiei ssp. duthiei and H. duthiei ssp. longearistata should be transferred to Leymus.

Phylogeography of the Qinghai-Tibetan Plateau endemic Juniperus przewalskii (Cupressaceae) inferred from chloroplast DNA sequence variation

The vegetation of the northeast Qinghai-Tibetan Plateau is dominated by alpine meadow and desert-steppe with sparse forests scattered within it. To obtain a better understanding of the phylogeography of one constituent species of the forests in this region, we examined chloroplast trnT-trnF and trnS-trnG sequence variation within Juniperus przewalskii, a key endemic tree species. Sequence data were obtained from 392 trees in 20 populations covering the entire distribution range of the species. Six cpDNA haplotypes were identified. Significant population subdivision was detected (G(ST) = 0.772, N-ST = 0.834), suggesting low levels of recurrent gene flow among populations and significant phylogeographic structure (N-ST > G(ST), P < 0.05). Eight of the nine disjunct populations surveyed on the high-elevation northeast plateau were fixed for a single haplotype (A), while the remaining, more westerly population, contained the same haplotype at high frequency together with two low frequency haplotypes (C and F). In contrast, most populations that occurred at lower altitudes at the plateau edge were fixed or nearly fixed for one of two haplotypes, A or E. However, two plateau edge populations had haplotype compositions different from the rest. In one, four haplotypes (A, B, D and E) were present at approximately equivalent frequencies, which might reflect a larger refugium in the area of this population during the last glacial period. Phylogenetic analysis indicated that the most widely distributed haplotype A is not ancestral to other haplotypes. The contrasting phylogeographic structures of the haplotype-rich plateau edge area and the almost haplotype-uniform plateau platform region indicate that the plateau platform was recolonized by J. przewalskii during the most recent postglacial period. This is supported by the findings of a nested clade analysis, which inferred that postglacial range expansion from the plateau edge followed by recent fragmentation is largely responsible for the present-day spatial distribution of cpDNA haplotypes within the species.