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.

Purification of High-quality Plastid DNA from the Brown Alga Laminaria japonica Sporophytes

Extracting DNA from a variety of algae is rather difficult because of high levels of polysaccharides, tannins, and phenolics as these interfere with DNA isolation and downstream applications. High-quality plastid DNA (ptDNA) purification is particularly difficult because of its small proportion in total genomic DNA. This report describes an improved protocol for ptDNA purification that efficiently produces high-quality ptDNA from sporophytes of Laminaria japonica and several other algae. This improved protocol simplifies procedures for ptDNA purification and improves yield to 150-200 mu g of ptDNA per 100 g of frozen algal tissue. Polymerase chain reaction (PCR) amplification of conserved sequences has been used to verify purity of the ptDNA product.

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.

Ulva and Enteromorpha (Ulvaceae, Chlorophyta) from two sides of the Yellow Sea: analysis of nuclear rDNA ITS and plastid rbcL sequence data

Ulvacean green seaweeds are common worldwide; they formed massive green tides in the Yellow Sea in recent years, which caused marine ecological problems as well as a social issue. We investigated two major genera of the Ulvaceae, Ulva and Enteromorpha, and collected the plastid rbcL and nuclear ITS sequences of specimens of the genera in two sides of the Yellow Sea and analyzed them. Phylogenetic trees of rbcL data show the occurrence of five species of Enteromorpha (E. compressa, E. flexuosa, E. intestinalis, E. linza and E. prolifera) and three species of Ulva (U. pertusa, U. rigida and U. ohnoi). However, we found U. ohnoi, which is known as a subtropical to tropical species, at two sites on Jeju Island, Korea. Four ribotypes in partial sequences of 5.8S rDNA and ITS2 from E. compressa were also found. Ribotype network analysis revealed that the common ribotype, occurring in China, Korea and Europe, is connected with ribotypes from Europe and China/Japan. Although samples of the same species were collected from both sides of the Yellow Sea, intraspecific genetic polymorphism of each species was low among samples collected worldwide.

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.

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.

Molecular systematics and biogeography of Crawfurdia, Metagentiana and Tripterospermum (Gentianaceae) based on nuclear ribosomal and plastid DNA sequences

Background and Aims The systematic position of the genus Metagentiana and its phylogenetic relationships with Crawfurdia, Gentiana and Tripterospermum have not been explicitly addressed. These four genera belong to one of two subtribes (Gentianinae) of Gentianeae. The aim of this paper is to examine the systematic position of Crawfurdia, Metagentiana and Tripterospermum and to clarify their phylogenetic affinities more clearly using ITS and trnL intron sequences.Methods Nucleotide sequences from the internal transcribed spacers (ITS) of nuclear ribosomal DNA and the plastid DNA trnL (UAA) intron were analysed phylogenetically. Ten of fourteen Metagentiana species were sampled, together with 40 species of other genera in the subtribe Gentianinae.Key Results The data support several previously published conclusions relating to the separation of Metagentiana from Gentiana and its closer relationships to Crawfurdia and Tripterospermum based on studies of gross morphology, floral anatomy, chromosomes, palynology, embryology and previous molecular data. The molecular clock hypothesis for the tested sequences in subtribe Gentianinae was not supported by the data (P < 0.05), so the clock-independent non-parametric rate smoothing method was used to estimate divergence time. This indicates that the separation of Crawfurdia, Metagentiana and Tripterospermum from Gentiana occurred about 11.4-21.4 Mya (million years ago), and the current species of these three genera diverged at times ranging from 0.4 to 6.2 Mya.Conclusions The molecular analyses revealed that Crawfurdia, Metagentiana and Tripterospermum do not merit status as three separate genera, because sampled species of Crawfurdia and Tripterospermum are embedded within Metagentiana. The speciation and rapid radiation of these three genera is likely to have occurred in western China as a result of upthrust of the Himalayas during the late Miocene and the Pleistocene.