Sequencing and phylogenetic analysis of the Pyrgilauda ruficollis(Aves,Passeridae) complete mitochondrial genome

In this study,both long PCR and conserved primers walking sequencing methods were used to determine the complete sequence of the of Pyrgilauda ruficollis mitochondrial genome(KC836121).The results showed that the complete mitochondrial genome of P.ruficollis is 16909 bp in length with 55.0%A+T content,harboring the typical 37 genes.The mitogenome had the same gene order with that of Podoces hendersoni.All protein coding genes started with ATG codon,except ND3 with GTG.For the stop codon usage,most genes terminate with codons TAA or TAG,but ND5 terminated with AGA,while ND1 and COI genes with AGG,and both the genes COIII and ND4 have an incomplete termination codon(T).The secondary structures of 22 tRNA genes were also predicted,showing that all tRNAs can form typical clover-leaf secondary structures,except for the tRNASer(AGN)which loses the DHU arm,while tRNAPhe harbor an extra nucleotide inserted in the TψC arm.The predicted secondary structures of 12S rRNA and16S rRNA exhibit 47 helices in 4 domains and 60 helices in 6 domains respectively.The control region of P.ruficollis with the length of 1 305 bp was located between tRNAGlu and tRNAPhe,and typical domains of which could be found as other bird groups.Using the data from 13 mitochondrial protein-coding genes,results of a final phylogenetic analysis strongly supports the traditional view that P.ruficollis is closely related with Passeridae and Fringillidae.     

Complete Mitochondrial Genome of the Steppe Ribbon Racer(Psammophis lineolatus): The First Representative from the Snake Family Psammophiidae and its Phylogenetic Implications

Comparisons of vertebrate mitochondrial genomes(mitogenomes) may yield significant insights into the evolution of organisms and genomes. However,no complete mitogenome from the snake family Psammophiidae has been reported. In this study, we sequenced and annotated the complete mitogenome of Psammophis lineolatus, representing the first mitogenome of Psammophiidae. The total length is 17 166 bp,consisting of 13 protein-coding genes(PCGs), 22 transfer RNAs(tRNAs), two ribosomal RNAs(12 S rRNA and16 S rRNA), and duplicate control regions(CRs). This gene arrangement belongs to the Type Ⅲ pattern,which is a widely shared gene order in Alethinophidian snakes. All tRNAs exhibit cloverleaf structures with the exception of tRNA-SerAGY and tRNA-Pro, which lack a dihydrouridine(DHU) arm/stem and TΨC loop,respectively. The 13 PCGs include five start codons(ATG,GTG, ATA, ATT, and ATC), two complete stop codons(TAA and AGG), and two incomplete stop codon(T--and TA-). In addition, the Ka/Ks ratios indicate that all PCGs had undergone a strong purifying selection.Four types of CR domains rearrangement occurred among eight species of Elapoidea. The phylogenetic reconstructions with both Bayesian inference and maximum likelihood methods support the placement of Psammophiidae in the Elapoidea superfamily,with Homalopsidae being the sister taxon to Elapoidea and Colubroidea. However, the sister taxon of Psammophiidae is unclear due to the availability of Elapoidea mitogenomes being limited to the family Elapidae. More mitogenomes from different taxonomic groups in Elapoidea are needed to better understand the phylogenetic relationships within Elapoidea.     

New selenoproteins identified in silico from the genome of Anopheles gambiae

Selenoprotein is biosynthesized by the incorporation of selenocysteine into proteins,where the TGA codon in the open reading frame does not act as a stop signal but is translated into selenocysteine.The dual functions of TGA result in mis-annotation or lack of selenoproteins in the sequenced genomes of many species.Available computational tools fail to correctly predict selenoproteins.Thus,we devel-oped a new method to identify selenoproteins from the genome of Anopheles gambiae computationally.Based on released genomic information,several programs were edited with PERL language to identify selenocysteine insertion sequence(SECIS)element,the coding potential of TGA codons,and cys-teine-containing homologs of selenoprotein genes.Our results showed that 11365 genes were termi-nated with TGA codons,918 of which contained SECIS elements.Similarity search revealed that 58 genes contained Sec/Cys pairs and similar flanking regions around in-frame TGA codons.Finally,7 genes were found to fully meet requirements for selenoproteins,although they have not been anno-tated as selenoproteins in NCBI databases.Deduced from their basic properties,the newly found se-lenoproteins in the genome of Anopheles gambiae are possibly related to in vivo oxidation tolerance and protein regulation in order to interfere with anopheles' vectorial capacity of Plasmodium.This study may also provide theoretical bases for the prevention of malaria from anopheles transmission.     

The Mitochondrial Genome of Raphanus sativus and Gene Evolution of Cruciferous Mitochondrial Types

To explore the mitochondrial genes of the Cruciferae family, the mitochondrial genome of Raphanus sativus (sat) was sequenced and annotated. The circular mitochondrial genome of sat is 239,723 bp and includes 33 protein-coding genes, three rRNA genes and 17 tRNA genes. The mitochondrial genome also contains a pair of large repeat sequences 5.9 kb in length, which may mediate genome reorga-nization into two sub-genomic circles, with predicted sizes of 124.8 kb and 115.0 kb, respectively. Furthermore, gene evolution of mitochondrial genomes within the Cruciferae family was analyzed using sat mitochondrial type (mitotype), together with six other re-ported mitotypes. The cruciferous mitochondrial genomes have maintained almost the same set of functional genes. Compared with Cycas taitungensis (a representative gymnosperm), the mitochondrial genomes of the Cruciferae have lost nine protein-coding genes and seven mitochondrial-like tRNA genes, but acquired six chloroplast-like tRNAs. Among the Cruciferae, to maintain the same set of genes that are necessary for mitochondrial function, the exons of the genes have changed at the lowest rates, as indicated by the numbers of single nucleotide polymorphisms. The open reading frames (ORFs) of unknown function in the cruciferous genomes are not conserved. Evolutionary events, such as mutations, genome reorganizations and sequence insertions or deletions (indels), have resulted in the non- conserved ORFs in the cruciferous mitochondrial genomes, which is becoming significantly different among mitotypes. This work represents the first phylogenic explanation of the evolution of genes of known function in the Cruciferae family. It revealed significant variation in ORFs and the causes of such variation.     

Phylogenetic placement of Cynomorium in Rosales inferred from sequences of the inverted repeat region of the chloroplast genome

Cynomorium is a herbaceous holoparasite that has been placed in Santalales, Saxifragales, Myrtales, or Sapindales. The inverted repeat （IR） region of the chloroplast genome region is slow evolving and, unlike mitochondrial genes, the chloroplast genome experiences few horizontal gene transfers between the host and parasite. Thus, in the present study, we used sequences of the IR region to test the phylogenetic placements of Cynomorium. Phylogenetic analyses of the chloroplast IR sequences generated largely congruent ordinal relationships with those from previous studies of angiosperm phylogeny based on single or multiple genes. Santalales was closely related to Caryophyllales and asterids. Saxifragales formed a clade where Peridiscus was sister to the remainder of the order, whereas Paeonia was sister to the woody clade of Saxifragales. Cynomorium is not closely related to Santalales, Saxifragales, Myrtales, or Sapindales; instead, it is included in Rosales and sister to Rosaceae. The various placements of the holoparasite on the basis of different regions of the mitochondrial genome may indicate the heterogeneous nature of the genome in the parasite. However, it is unlikely that the placement of Cynomorium in Rosales is the result of chloroplast gene transfer because Cynomorium does not parasitize on rosaceous plants and there is no chloroplast gene transfer between Cynomorium and Nitraria, a confirmed host of Cynomorium and a member of Sapindales.     

The complete nucleotide sequence of the mitochondrial genome of the cabbage butterfly,Artogeia melete （Lepidoptera： Pieridae）

The complete mitochondrial genome （mitogenome） of Artogeia melete was determined as being composed of 15,140 bp, including 13 protein-coding genes （PCGs）, 2 rRNA genes, 22 tRNA genes, and one control region. The gene order of A. melete mitogenome is typical of Lepidoptera and differs from the insect ancestral type in the location of trnM. The A. melete mitogenome has a total of 119 bp of intergenic spacer sequences spread over 10 regions, ranging in sizes between 1 and 48 bp. The nucleotide composition of the A. melete mitogenome is also biased toward A ＋ T nucleotides （79.77%）, which is higher than that of Ochrogaster lunifer （77.84%）, but lower than nine other lepidopterans sequenced. The PCGs have typical mitochondrial start codons, except for coxl, which contains the unusual CGA. The coxl, cox2, nad2, and nad5 genes of the A. melete mitogenome have incomplete stop codons （T）. The A. melete A ＋ T-rich region contains some conserved structures that are similar to those found in other lepidopteran mitogenomes, including a structure combining the motif ‘ATAGA＇, a 19-bp poly（T） stretch, a microsatellite （AT）n element, and a 9-bp poly（A） upstream trnM. The A. melete mitogenome contains a duplicated 36-bp repeat element, which consists of a 26- bp core sequence flanked by 10-bp perfectly inverted repeats.     

Nelumbonaceae:Systematic position and species diversification revealed by the complete chloroplast genome

<正>Nelumbonaceae is a morphologically unique family of angiosperms and was traditionally placed in Nymphaeales;more recently,it was placed in Proteales based on molecular data,or in an order of its own,Nelumbonales. To determine the systematic position of the family and to date the divergence time of the family and the divergence time of its two intercontinentally disjunct species,we sequenced the entire chloroplast genome of Nelumbo lutea and most of the chloroplast genes of,N.nucifera.We carried out phylogenetic and molecular dating analyses of the two species and representatives of 47 other plant families,representing the major lineages of angiosperms, using 83 plastid genes.The N.lutea genome was 163 510 bp long,with a total of 130 coding genes and an overall GC content of 38%.No significant structural differences among the genomes of N.lutea,Nymphaea alba, and Platanus occidentalis were observed.The phylogenetic relationships based on the 83 plastid genes revealed a close relationship between Nelumbonaceae and Platanaceae.The divergence times were estimated to be 109 Ma between the two families and 1.5 Ma between the two Nelumbo species.The estimated time was only slightly longer than the age of known Nelumbo fossils,suggesting morphological stasis within Nelumbonaceae.We conclude that Nelumbonaceae holds a position in or close to Proteales.We further conclude that the two species of Nelumbo diverged recently from a common ancestor and do not represent ancient relicts on different continents.     

Comparative Genomic Study of the Thioredoxin Family in Photosynthetic Organisms with Emphasis on Populus trichocarpa

The recent genome sequencing of Populus trichocarpa and Vitis vinifera, two models of woody plants, of Sorghum bicolor, a model of monocot using C4 metabolism, and of the moss Physcomitrella patens, together with the availability of photosynthetic organism genomes allows performance of a comparative genomic study with organisms having different ways of life, reproduction modes, biological traits, and physiologies. Thioredoxins （Trxs） are small ubiq- uitous proteins involved in the reduction of disulfide bridges in a variety of target enzymes present in all sub-cellular compartments and involved in many biochemical reactions. The genes coding for these enzymes have been identified in these newly sequenced genomes and annotated. The gene content, organization and distribution were compared to other photosynthetic organisms, leading to a refined classification. This analysis revealed that higher plants and bryo- phytes have a more complex family compared to algae and cyanobacteria and to non-photosynthetic organisms, since poplar exhibits 49 genes coding for typical and atypical thioredoxins and thioredoxin reductases, namely one-third more than monocots such as Oryza sativa and S. bicolor. The higher number of Trxs in poplar is partially explained by gene duplication in the Trx m, h, and nucleoredoxin classes. Particular attention was paid to poplar genes with emphasis on Trx-like classes called Clot, thioredoxin-like, thioredoxins of the lilium type and nucleoredoxins, which were not described in depth in previous genomic studies.     

Complete mitochondrial genome of Oxya chinensis （Orthoptera, Acridoidea）

The complete sequence of Oxya chinensis （0. chinensis） mitochondrial genome is reported here. It is 15,443 bp in length and contains 75.9% A＋T. The protein-coding genes have a similar A＋T content （75.2%）. The initiation codon of the cytochrome oxidase subunit I gene in the mitochondrial genome of O. chinensis appears to be ATC, instead of the tetranucleotides that have been reported in Locusta migratoria （L migratoria） mitochondrial genome. The sizes of the large and small ribosomal RNA genes are 1319 and 850 bp, respectively. The transfer RNA genes have been modeled and showed strong resemblance to the dipteran transfer RNAs, and all anticodons are identical to those of dipteran. The A＋T-rich region is 562 bp, shorter than that of other known Orthoptera insects. The six conserved domains were identified within the A＋T-rich region by comparing its sequence with those of other grasshoppers. The result of phylogenetic analysis based on the dataset containing 12 concatenated protein sequences confirms the close relation-ship of O. chinensis with L migratoria.     

The complete mitochondrial genome sequence of Geisha distinctissima （Hemiptera： Flatidae） and comparison with other hemipteran insects

The complete nucleotide sequence of the mitochondrial genome （mitogenome） of Geisha distinctissima （Hemiptera： Flatidae） has been determined in this study. The genome is a circular molecule of 15,971 bp with a total A＋T content of 75.1%. The gene content, order, and structure are consistent with the Drosophila yakuba genome structure and the hypothesized ancestral arthro- pod genome arrangement. All 13 protein-coding genes are observed to have a putative, inframe ATR methionine or ATT isoleucine codons as start signals. Canonical TAA and TAG termination codons are found in nine protein-coding genes, and the remaining four （cox1, atp6, cox3, and nad4） have incomplete termination codons. The anticodons of all transfer RNA （tRNAs） are identical to those observed in D. yakuba and Philaenus spumarius, and can be folded in the form of a typical clover-leaf structure except for tRNASer（AGN）. The major non-coding region （the A ＋ T-rich region or putative control region） between the small ribosomal subunit and the tRNAne gene includes two sets of repeat regions. The first repeat region consists of a direct 152-bp repetitive unit located near the srRNA gene end, and the second repeat region is composed of a direct repeat unit of 19 bp located toward tRNAIle gene. Comparisons of gene variability across the order suggest that the gene content and arrangement of G. distinctissima mitogenome are similar to other hemipteran insects.     

The complete mitochondrial genome and phylogenetic analysis of Nyctereutes procyonoides

The complete mitochondrial genome sequence of the raccoon dog (Nyctereutes procyonoides) was determined by using the long and accurate polymerase chain reaction. The entire mitochondrial genome sequence is 16,713 bp in length contains two ribosomal RNA genes, 13 protein-coding genes, 22 transfer RNA genes and 1 control region. Most mitochondrial genes are encoded on the H strand, except for the ND6 gene and 8 tRNA genes. The base compositions of mitochondrial genomes present clearly A–T skew. All the transfer RNA genes can be folded into the typical cloverleaf-shaped structure except tRNA-Ser (AGY), which lacks the dihydrouridine arm. Protein-coding genes mainly initiate with ATG and terminate with TAA. Some reading frame intervals and overlaps are found in the mitochondrial genome. The control region can be divided into three domains: the extended termination associated sequences (ETASs) domain, the central conserved domain and the conserved sequence blocks (CSBs) domain. Three conserved sequence blocks (CSBs) and one extended termination associated sequences (ETAS-1) is found in the control region. The phylogenetic analysis based on the concatenated data set of 14 genes in the mitochondrial genome of Canidae shows that the raccoon dog has close phylogenetic position with the red fox (Vulpes vulpes) and they constitute a clade which has an equil evolutionary position with the clade formed by the genera Canis and Cuon.     

SCGPred： A Score-based Method for Gene Structure Prediction by Combining Multiple Sources of Evidence

Predicting protein-coding genes still remains a significant challenge. Although a variety of computational programs that use commonly machine learning methods have emerged, the accuracy of predictions remains a low level when implementing in large genomic sequences. Moreover, computational gene finding in newly se- quenced genomes is especially a difficult task due to the absence of a training set of abundant validated genes. Here we present a new gene-finding program, SCGPred, to improve the accuracy of prediction by combining multiple sources of evidence. SCGPred can perform both supervised method in previously well-studied genomes and unsupervised one in novel genomes. By testing with datasets composed of large DNA sequences from human and a novel genome of Ustilago maydi, SCGPred gains a significant improvement in comparison to the popular ab initio gene predictors. We also demonstrate that SCGPred can significantly improve prediction in novel genomes by combining several foreign gene finders with similarity alignments, which is superior to other unsupervised methods. Therefore, SCGPred can serve as an alternative gene-finding tool for newly sequenced eukaryotic genomes. The program is freely available at http：//bio.scu.edu.cn/SCGPred/.     

Complete chloroplast genome sequence of Magnolia grandiflora and comparative analysis with related species

Magnolia grandiflora is an important medicinal,ornamental and horticultural plant species.The chloroplast(cp) genome of M.grandiflora was sequenced using a 454 sequencing platform and the genome structure was compared with other related species.The complete cp genome of M.grandiflora was 159623 bp in length and contained a pair of inverted repeats(IR) of 26563 bp separated by large and small single copy(LSC,SSC) regions of 87757 and 18740 bp,respectively.A total of 129 genes were successfully annotated,18 of which included introns.The identity,number and GC content of M.grandiflora cp genes were similar to those of other Magnoliaceae species genomes.Analysis revealed 218 simple sequence repeat(SSR) loci,most composed of A or T,contributing to a bias in base composition.The types and abundances of repeat units in Magnoliaceae species were relatively conserved and these loci will be useful for developing M.grandiflora cp genome vectors.In addition,results indicated that the cp genome size in Magnoliaceae species and the position of the IR border were closely related to the length of the ycf1 gene.Phylogenetic analyses based on 66 shared genes from 30 species using maximum parsimony(MP) and maximum likelihood(ML) methods provided strong support for the phylogenetic position of Magnolia.The availability of the complete cp genome sequence of M.grandiflora provides valuable information for breeding of desirable varieties,cp genetic engineering,developing useful molecular markers and phylogenetic analyses in Magnoliaceae.     

The E protein is a multifunctional membrane protein of SARS-CoV

The E (envelope) protein is the smallest structural protein in all coronaviruses and is the only viral structural protein in which no variation has been detected. We conducted genome sequencing and phylogenetic analyses of SARS-CoV. Based on genome sequencing, we predicted the E protein is a transmembrane (TM) protein characterized by a TM region with strong hydrophobicity and α-helix conformation. We identified a segment (NH2-_L-Cys-A-Y-Cys-Cys-N_-COOH) in the carboxyl-terminal region of the E protein that appears to form three disulfide bonds with another segment of corresponding cysteines in the carboxyl-terminus of the S (spike) protein. These bonds point to a possible structural association between the E and S proteins. Our phylogenetic analyses of the E protein sequences in all published coronaviruses place SARS-CoV in an independent group in Coronaviridae and suggest a non-human animal origin.     

The transfer RNA genes in Oryza sativa L. ssp. indica

The availability of the draft genome sequence of Oryza sativa L ssp. indica has made it possible to study the rice tRNA genes. A total of 596 tRNA genes, including 3 selenocysteine tRNA genes and one suppressor tRNA gene are identified in 127551 rice contigs. There are 45 species of tRNA genes and the revised wobble hypothesis proposed by Guthrie and Abelson is perfectly obeyed. The relationship between codon usage and the number of corresponding tRNA genes is discussed. Redundancy may exist in the present list of tRNA genes and novel ones may be found in the future. A set of 33 tRNA genes is discovered in the complete chloroplast genome of Oryza sativa L. ssp. indica. These tRNA genes are identical to those in ssp. japonica identified by us independently from the origional annotation.     

Molecular phylogenetic systematics of twelve species of Acipenseriformes based on mtDNA ND4L -ND4 gene sequence analysis

Acipenseriformes is an endangered primitive fish group, which occupies a special place in the history of ideas concerning fish evolution, even in vertebrate evolution. However, the classification and evolution of the fishes have been debated. The mitochondrial DMA (mtDNA) ND4L and partial A7D4 genes were first sequenced in twelve species of the order Acipenseriformes, including endemic Chinese species. The following points were drawn from DNA sequences analysis: (i) the two species of Huso can be ascribed to Acipenser; (ii) A. dabryanus is the mostly closely related to A. sinensis, and most likely the landlocked form of A. sinensis; (iii) genus Acipenser in trans-Pacific region might have a common origin; (iv) mtDNA ND4L and ND4 genes are the ideal genetic markers for phylogenetic analysis of the order Acipenseriformes.     

The revolution of the biology of the genome

Sequence data of entire eukaryotic genomes and their detailed comparison have provided new evidence on genome evolution. The major mechanisms involved in the increase of genome sizes are polyploidization and gene duplication.Subsequent gene silencing or mutations, preferentially in regulatory sequences of genes, modify the genome and permit the development of genes with new properties. Mechanisms such as lateral gene transfer, exon shuffling or the creation of new genes by transposition contribute to the evolution of a genome, but remain of relatively restricted relevance.Mechanisms to decrease genome sizes and, in particular, to remove specific DNA sequences, such as blocks of satellite DNAs, appear to involve the action of RNA interference (RNAi). RNAi mechanisms have been proven to be involved in chromatin packaging related with gene inactivation as well as in DNA excision during the macronucleus development in ciliates.     

Study of Completed Archaeal Genomes and Proteomes： Hypothesis of Strong Mutational AT Pressure Existed in Their Common Predecessor

The number of completely sequenced archaeal genomes has been sufficient for a large-scale bioinformatic study.We have conducted analyses for each coding region from 36 archaeal genomes using the original CGS algorithm by calculating the total GC content(G+C),GC content in first,second and third codon positions as well as in fourfold and twofold degenerated sites from third codon positions,levels of arginine codon usage(Arg2:AGA/G;Arg4:CGX),levels of amino acid usage and the entropy of amino acid content distribution.In archaeal genomes with strong GC pressure,arginine is coded preferably by GC-rich Arg4 codons,whereas in most of archaeal genomes with G+C0.6,arginine is coded preferably by AT-rich Arg2 codons.In the genome of Haloquadratum walsbyi,which is closely related to GC-rich archaea,GC content has decreased mostly in third codon positions,while Arg4Arg2 bias still persists.Proteomes of archaeal species carry characteristic amino acid biases:levels of isoleucine and lysine are elevated,while levels of alanine,histidine,glutamine and cytosine are relatively decreased.Numerous genomic and proteomic biases observed can be explained by the hypothesis of previously existed strong mutational AT pressure in the common predecessor of all archaea.