Predicted Coding Sequences Research Articles

Inferring bona fide transfrags in RNA-Seq derived-transcriptome assemblies of non-model organisms.

BackgroundDe novo transcriptome assembly of short transcribed fragments (transfrags) produced from sequencing-by-synthesis technologies often results in redundant datasets with differing levels of unassembled, partially assembled or mis-assembled transcripts. Post-assembly processing intended to reduce redundancy typically involves reassembly or clustering of assembled sequences. However, these approaches are mostly based on common word heuristics and often create clusters of biologically unrelated sequences, resulting in loss of unique transfrags annotations and propagation of mis-assemblies.ResultsHere, we propose a structured framework that consists of a few steps in pipeline architecture for Inferring Functionally Relevant Assembly-derived Transcripts (IFRAT). IFRAT combines 1) removal of identical subsequences, 2) error tolerant CDS prediction, 3) identification of coding potential, and 4) complements BLAST with a multiple domain architecture annotation that reduces non-specific domain annotation. We demonstrate that independent of the assembler, IFRAT selects bona fide transfrags (with CDS and coding potential) from the transcriptome assembly of a model organism without relying on post-assembly clustering or reassembly. The robustness of IFRAT is inferred on RNA-Seq data of Neurospora crassa assembled using de Bruijn graph-based assemblers, in single (Trinity and Oases-25) and multiple (Oases-Merge and additive or pooled) k-mer modes. Single k-mer assemblies contained fewer transfrags compared to the multiple k-mer assemblies. However, Trinity identified a comparable number of predicted coding sequence and gene loci to Oases pooled assembly. IFRAT selects bona fide transfrags representing over 94% of cumulative BLAST-derived functional annotations of the unfiltered assemblies. Between 4-6% are lost when orphan transfrags are excluded and this represents only a tiny fraction of annotation derived from functional transference by sequence similarity. The median length of bona fide transfrags ranged from 1.5kb (Trinity) to 2kb (Oases), which is consistent with the average coding sequence length in fungi. The fraction of transfrags that could be associated with gene ontology terms ranged from 33-50%, which is also high for domain based annotation. We showed that unselected transfrags were mostly truncated and represent sequences from intronic, untranslated (5′ and 3′) regions and non-coding gene loci.ConclusionsIFRAT simplifies post-assembly processing providing a reference transcriptome enriched with functionally relevant assembly-derived transcripts for non-model organism.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-015-0492-5) contains supplementary material, which is available to authorized users.

Comparative genomic analysis and phenazine production of Pseudomonas chlororaphis, a plant growth-promoting rhizobacterium.

Pseudomonas chlororaphis HT66, a plant growth-promoting rhizobacterium that produces phenazine-1-carboxamide with high yield, was compared with three genomic sequenced P. chlororaphis strains, GP72, 30–84 and O6. The genome sizes of four strains vary from 6.66 to 7.30 Mb. Comparisons of predicted coding sequences indicated 4833 conserved genes in 5869–6455 protein-encoding genes. Phylogenetic analysis showed that the four strains are closely related to each other. Its competitive colonization indicates that P. chlororaphis can adapt well to its environment. No virulence or virulence-related factor was found in P. chlororaphis. All of the four strains could synthesize antimicrobial metabolites including different phenazines and insecticidal protein FitD. Some genes related to the regulation of phenazine biosynthesis were detected among the four strains. It was shown that P. chlororaphis is a safe PGPR in agricultural application and could also be used to produce some phenazine antibiotics with high-yield.

Genome sequence and description of the mosquitocidal and heavy metal tolerant strain Lysinibacillus sphaericus CBAM5.

Lysinibacillus sphaericus CBAM5, was isolated from subsurface soil of oil well explorations in the Easter Planes of Colombia. This strain has potential in bioremediation of heavy-metal polluted environments and biological control of Culex quinquefasciatus. According to the phylogenetic analysis of 16S rRNA gene sequences, the strain CBAM5 was assigned to the Lysinibacillus sphaericus taxonomic group 1 that comprises mosquito pathogenic strains. After a combination assembly-integration, alignment and gap-filling steps, we propose a 4,610,292 bp chromosomal scaffold. The whole genome (consisting of 5,146,656 bp long, 60 contigs and 5,209 predicted-coding sequences) revealed strong functional and syntenial similarities to the L. sphaericus C3-41 genome. Mosquitocidal (Mtx), binary (Bin) toxins, cereolysin O, and heavy metal resistance clusters from nik, ars, czc, mnt, ter, cop, cad, and znu operons were identified.

The identification of novel and differentially expressed apple-tree genes under low-temperature stress using high-throughput Illumina sequencing.

Low temperature is a major environmental constraint on the production of apples worldwide. An apple rootstock with high cold tolerance was selected to identify genes related to stress tolerance. The transcriptional profiles of the genes in the leaves were examined after various intervals of exposure to cold stress. We obtained three libraries of 14,219, 11,176 and 16,116 tag-mapped predicted coding sequences at three time points (0, 1 and 6 h) during cold stress. In the two time periods, which were from 0 to 1 h and from 1 to 6 h, 139 and 1,085 genes were upregulated, and 1,499 and 381 genes were downregulated, respectively. These groups included a large number of unknown genes. The distribution of genes indicated cold adaptation in the plant. Most of the differential expression genes (DEGs) had temporal specificity and significantly different expression changes at different time points. The classification of DEGs by GO category and KEGG pathway analysis revealed that the DEGs are involved in numerous biological pathways, including metabolism, plant-pathogen interaction and signal transduction. Eleven randomly selected tag-mapped genes were examined by qRT-PCR. The results of the qRT-PCR were in accordance with the transcriptional profiles. The most upregulated gene (MDP0000198054) from 0 to 1 h encodes a dehydration-responsive element-binding protein/C-repeat factor (DREB/CBF). In this study, MDP0000198054 and related genes involved in the cold stress response were discussed. These results could provide new insights into the stress tolerance mechanisms of apple rootstocks.

Hepatic patatin-like phospholipase domain-containing protein 3 sequence, single nucleotide polymorphism presence, protein confirmation, and responsiveness to energy balance in dairy cows

Patatin-like phospholipase domain-containing protein 3 (PNPLA3), commonly known as adiponutrin, is part of a novel subfamily of triglyceride lipase enzymes with potential effects on triglyceride metabolism in adipose and hepatic tissues. The predicted bovine PNPLA3 sequence has been identified, but expression of the gene had not been examined. The objectives of this study were to confirm the predicted bovine PNPLA3 gene sequence, determine expression of the bovine PNPLA3 gene in response to whole-animal energy balance, identify single nucleotide polymorphisms present in dairy cows, and verify the presence of the protein in the liver. Using liver biopsy samples collected from cows at +28d relative to calving (DRTC), RNA was isolated and used to generate a cDNA template for amplification of the entire predicted coding sequence of PNPLA3 via PCR. To determine if energy balance alters the expression of PNPLA3, RNA was isolated and mRNA expression quantified in liver samples from mid-lactation cows after a 5-d ad libitum period (n=5) and after a subsequent 5-d 50% feed restriction period (n=5), and in samples collected from cows at −14, +1, +14, and +28 DRTC (n=16). The presence of PNPLA3 protein was detected by Western blot in liver protein samples collected at +28 DRTC. Expression of hepatic PNPLA3 was decreased after a period of feed restriction (8.14 vs. 1.08±2.17 arbitrary units, ad libitum vs. fasted). Expression of PNPLA3 mRNA was decreased at +1 and +14 DRTC compared with −14 DRTC (23.35, 7.28, 10.17, and 14.5±4.9 arbitrary units, −14, +1, +14, and +28 DRTC, respectively). The presence of PNPLA3 protein was detected as a 55-kDa band in hepatic protein isolations from liver tissue collected at +28 DRTC. These data confirm the presence and sequence of the bovine hepatic PNPLA3 gene and single nucleotide polymorphisms. Furthermore, these data indicate responsiveness of bovine hepatic PNPLA3 to energy balance.

Patterns of simple sequence repeats in cultivated blueberries (Vaccinium section Cyanococcus spp.) and their use in revealing genetic diversity and population structure

Blueberry (Vaccinium section Cyanococcus) is an important small fruit crop native to North America. Relationships among the primary species known as ‘blueberry’ has been a source of speculation due to the out-breeding nature of the crop, the use of intra-specific hybridization and open-pollinated selection in breeding programs, and the lack of genomic resources to adequately address the issue. The objectives of this study were to characterize simple sequence repeats (SSRs) from an emerging genomic draft sequence, develop useful molecular markers and provide an in-depth analysis of genetic diversity and population structure in blueberry using a broad range of cultivated blueberry accessions representing multiple species, ploidy levels and sources of origin. Genomic scaffolding was assembled from the whole-genome sequencing of a diploid V. corymbosum accession ‘W8520.’ From the assembled 358 Mb sequence, a total number of 43,594 SSRs were identified (122 per Mb). Among genomic regions, SSRs were longest and occurred most frequently in predicted 5′ untranslated regions (5′ UTR), while SSRs were shortest and least common in the predicted coding sequences. AG/CT and AAG/CTT were the most frequent motifs while CG/CG and CCG/CGG motifs were the rarest di- and trinucleotide motifs, respectively. For analysis of genetic diversity and population structure, 42 genomic SSR and EST–SSR markers with an average of 14.2 alleles and 56.0 allele phenotypes per locus were used to genotype a diverse blueberry population of 150 accessions. Cluster analysis grouped the 150 accessions in a manner consistent with known information regarding species, ploidy levels and pedigree. The analysis of population structure among blueberry accessions revealed inter- and intra-specific levels of stratification. Rabbiteye blueberry (V. virgatum) represents a genetically distinct subgroup within Cyanococcus. Three additional subpopulations were detected among highbush varieties that are largely attributable to distinctions between northern and southern highbush and founder effects of a single cultivar (‘Weymouth’). The identification of substructure that correlates with known pedigree information, and the availability of new genomic molecular markers will facilitate future evolutionary and genetic studies in blueberry.

Gene Content and Diversity of the Loci Encoding Biosynthesis of Capsular Polysaccharides of the 15 Serovar Reference Strains of Haemophilus parasuis

Haemophilus parasuis is the causative agent of Glässer's disease, a systemic disease of pigs, and is also associated with pneumonia. H. parasuis can be classified into 15 different serovars. Here we report, from the 15 serotyping reference strains, the DNA sequences of the loci containing genes for the biosynthesis of the group 1 capsular polysaccharides, which are potential virulence factors of this bacterium. We contend that these loci contain genes for polysaccharide capsule structures, and not a lipopolysaccharide O antigen, supported by the fact that they contain genes such as wza, wzb, and wzc, which are associated with the export of polysaccharide capsules in the current capsule classification system. A conserved region at the 3' end of the locus, containing the wza, ptp, wzs, and iscR genes, is consistent with the characteristic export region 1 of the model group 1 capsule locus. A potential serovar-specific region (region 2) has been found by comparing the predicted coding sequences (CDSs) in all 15 loci for synteny and homology. The region is unique to each reference strain with the exception of those in serovars 5 and 12, which are identical in terms of gene content. The identification and characterization of this locus among the 15 serovars is the first step in understanding the genetic, molecular, and structural bases of serovar specificity in this poorly studied but important pathogen and opens up the possibility of developing an improved molecular serotyping system, which would greatly assist diagnosis and control of Glässer's disease.

Comparative genomic analysis of four representative plant growth-promoting rhizobacteria in Pseudomonas

BackgroundSome Pseudomonas strains function as predominant plant growth-promoting rhizobacteria (PGPR). Within this group, Pseudomonas chlororaphis and Pseudomonas fluorescens are non-pathogenic biocontrol agents, and some Pseudomonas aeruginosa and Pseudomonas stutzeri strains are PGPR. P. chlororaphis GP72 is a plant growth-promoting rhizobacterium with a fully sequenced genome. We conducted a genomic analysis comparing GP72 with three other pseudomonad PGPR: P. fluorescens Pf-5, P. aeruginosa M18, and the nitrogen-fixing strain P. stutzeri A1501. Our aim was to identify the similarities and differences among these strains using a comparative genomic approach to clarify the mechanisms of plant growth-promoting activity.ResultsThe genome sizes of GP72, Pf-5, M18, and A1501 ranged from 4.6 to 7.1 M, and the number of protein-coding genes varied among the four species. Clusters of Orthologous Groups (COGs) analysis assigned functions to predicted proteins. The COGs distributions were similar among the four species. However, the percentage of genes encoding transposases and their inactivated derivatives (COG L) was 1.33% of the total genes with COGs classifications in A1501, 0.21% in GP72, 0.02% in Pf-5, and 0.11% in M18. A phylogenetic analysis indicated that GP72 and Pf-5 were the most closely related strains, consistent with the genome alignment results. Comparisons of predicted coding sequences (CDSs) between GP72 and Pf-5 revealed 3544 conserved genes. There were fewer conserved genes when GP72 CDSs were compared with those of A1501 and M18. Comparisons among the four Pseudomonas species revealed 603 conserved genes in GP72, illustrating common plant growth-promoting traits shared among these PGPR. Conserved genes were related to catabolism, transport of plant-derived compounds, stress resistance, and rhizosphere colonization. Some strain-specific CDSs were related to different kinds of biocontrol activities or plant growth promotion. The GP72 genome contained the cus operon (related to heavy metal resistance) and a gene cluster involved in type IV pilus biosynthesis, which confers adhesion ability.ConclusionsComparative genomic analysis of four representative PGPR revealed some conserved regions, indicating common characteristics (metabolism of plant-derived compounds, heavy metal resistance, and rhizosphere colonization) among these pseudomonad PGPR. Genomic regions specific to each strain provide clues to its lifestyle, ecological adaptation, and physiological role in the rhizosphere.

MetaSAMS—A novel software platform for taxonomic classification, functional annotation and comparative analysis of metagenome datasets

Metagenomics aims at exploring microbial communities concerning their composition and functioning. Application of high-throughput sequencing technologies for the analysis of environmental DNA-preparations can generate large sets of metagenome sequence data which have to be analyzed by means of bioinformatics tools to unveil the taxonomic composition of the analyzed community as well as the repertoire of genes and gene functions. A bioinformatics software platform is required that allows the automated taxonomic and functional analysis and interpretation of metagenome datasets without manual effort. To address current demands in metagenome data analyses, the novel platform MetaSAMS was developed. MetaSAMS automatically accomplishes the tasks necessary for analyzing the composition and functional repertoire of a given microbial community from metagenome sequence data by implementing two software pipelines: (i) the first pipeline consists of three different classifiers performing the taxonomic profiling of metagenome sequences and (ii) the second functional pipeline accomplishes region predictions on assembled contigs and assigns functional information to predicted coding sequences. Moreover, MetaSAMS provides tools for statistical and comparative analyses based on the taxonomic and functional annotations. The capabilities of MetaSAMS are demonstrated for two metagenome datasets obtained from a biogas-producing microbial community of a production-scale biogas plant. The MetaSAMS web interface is available at https://metasams.cebitec.uni-bielefeld.de.

Characterization and differential expression analysis of complete coding sequences of Vitis vinifera L. sirtuin genes

The sirtuin/Sir2 (Silent information regulator 2) family of NAD+-dependent deacetylases and mono-ADP-ribosyltransferases plays an important role in several cellular processes including gene silencing, cell cycle regulation and life span extension in yeast and animals. Compared to other eukaryotes, plants have relatively fewer SIR2 related genes encoding only two putative SIR2 family proteins. Recently, two putative sirtuin genes were identified also in the grapevine genome. Starting from the predicted coding sequences present in the database, we have been able to obtain two truly expressed coding sequences from the start to the stop codon for both sirtuin genes that were named VvSRT1 and VvSRT2. The search for the expressed coding sequences was performed by comparing the predicted sequences with the recently available grape RNA seq database with the aim to develop the primers to be used in reverse transcriptase PCR reactions to amplify the genes of interest. Finally, in order to better understand the physiological role of both sirtuins, we investigated the expression of these genes in young leaves, mature leaves, and berries sampled at different growing stages. In leaves, usually it has been observed that VvSRT1 is less expresses than VvSRT2, moreover in young leaves VvSRT2 showed the higher expression during setting while in mature leaves during the flowering time. No particular variations have been observed concerning VvSRT1. In berries the two genes showed more similar expression level, and they showed the highest expression during the flowering time. Finally, the expression of VvSRT2 in berries is smaller than in leaves.

Genome-Wide Definition of the SigF Regulon in Mycobacterium tuberculosis

In Mycobacterium tuberculosis the alternative sigma factor SigF controls the expression of a particular subset of genes by altering RNA polymerase specificity. Here, we utilize two genome-wide approaches to identify SigF-binding sites: chromatin immunoprecipitation (ChIP-on-chip) and microarray analysis of SigF-mediated transcripts. Since SigF is not an abundant protein in the logarithmic phase of growth, a pristinamyin IA-inducible system was used to control its expression. We identified 67 high-affinity SigF-binding sites and 16 loci where a SigF promoter directs the expression of a transcript. These loci include sigF itself, genes involved in lipid and intermediary metabolism and virulence, and at least one transcriptional regulator (Rv2884), possibly acting downstream of SigF. In addition, SigF was also found to direct the transcription of the gene for small RNA F6. Many loci were also found where SigF may be involved in antisense transcription, and in two cases (Rv1358 and Rv1870c) the SigF-dependent promoter was located within the predicted coding sequence. Quantitative PCR confirmed the microarray findings and 5'-rapid amplification of cDNA ends was used to map the SigF-specific transcriptional start points. A canonical SigF consensus promoter sequence GGTTT-N((15-17))-GGGTA was found prior to 11 genes. Together, these data help to define the SigF regulon and show that SigF not only governs expression of proteins such as the virulence factor, HbhA, but also impacts novel functions, such as noncoding RNAs and antisense transcripts.

Comparative analysis of a plant pseudoautosomal region (PAR) in Silene latifolia with the corresponding S. vulgaris autosome

BackgroundThe sex chromosomes of Silene latifolia are heteromorphic as in mammals, with females being homogametic (XX) and males heterogametic (XY). While recombination occurs along the entire X chromosome in females, recombination between the X and Y chromosomes in males is restricted to the pseudoautosomal region (PAR). In the few mammals so far studied, PARs are often characterized by elevated recombination and mutation rates and high GC content compared with the rest of the genome. However, PARs have not been studied in plants until now. In this paper we report the construction of a BAC library for S. latifolia and the first analysis of a > 100 kb fragment of a S. latifolia PAR that we compare to the homologous autosomal region in the closely related gynodioecious species S. vulgaris.ResultsSix new sex-linked genes were identified in the S. latifolia PAR, together with numerous transposable elements. The same genes were found on the S. vulgaris autosomal segment, with no enlargement of the predicted coding sequences in S. latifolia. Intergenic regions were on average 1.6 times longer in S. latifolia than in S. vulgaris, mainly as a consequence of the insertion of transposable elements. The GC content did not differ significantly between the PAR region in S. latifolia and the corresponding autosomal region in S. vulgaris.ConclusionsOur results demonstrate the usefulness of the BAC library developed here for the analysis of plant sex chromosomes and indicate that the PAR in the evolutionarily young S. latifolia sex chromosomes has diverged from the corresponding autosomal region in the gynodioecious S. vulgaris mainly with respect to the insertion of transposable elements. Gene order between the PAR and autosomal region investigated is conserved, and the PAR does not have the high GC content observed in evolutionarily much older mammalian sex chromosomes.

Genome-wide Association Study Identifies Four Genetic Loci Associated with Thyroid Volume and Goiter Risk

Thyroid disorders such as goiters represent important diseases, especially in iodine-deficient areas. Sibling studies have demonstrated that genetic factors substantially contribute to the interindividual variation of thyroid volume. We performed a genome-wide association study of this phenotype by analyzing a discovery cohort consisting of 3620 participants of the Study of Health in Pomerania (SHIP). Four genetic loci were associated with thyroid volume on a genome-wide level of significance. Of these, two independent loci are located upstream of and within CAPZB, which encodes the β subunit of the barbed-end F-actin binding protein that modulates actin polymerization, a process crucial in the colloid engulfment during thyroglobulin mobilization in the thyroid. The third locus marks FGF7, which encodes fibroblast growth factor 7. Members of this protein family have been discussed as putative signal molecules involved in the regulation of thyroid development. The fourth locus represents a "gene desert" on chromosome 16q23, located directly downstream of the predicted coding sequence LOC440389, which, however, had already been removed from the NCBI database as a result of the standard genome annotation processing at the time that this study was initiated. Experimental proof of the formerly predicted mature mRNA, however, demonstrates that LOC440389 indeed represents a real gene. All four associations were replicated in an independent sample of 1290 participants of the KORA study. These results increase the knowledge about genetic factors and physiological mechanisms influencing thyroid volume.

Purification to homogeneity and characterization of nonproteolyzed potato (Solanum tuberosum) tuber hexokinase 1

We have developed an extraction procedure that improves the stability of potato ( Solanum tuberosum L.) tuber hexokinase (HK) after extraction. Using this protocol, we showed that at least four HK isoforms are present in this tissue, and they can be separated by hydrophobic-interaction chromatography on a butyl-Sepharose™ 4 Fast Flow column. One of the main HK isoforms was purified to homogeneity using further chromatographic separations on red dye, DEAE Fractogel, hydroxyapatite, cibacron blue, and MonoQ matrices. HK-specific activity of this fraction (10.2 U·mg protein–1) corresponds to an enrichment of more than 5500-fold, with a yield of 0.9%. This is the highest reported HK-specific activity from a plant source. The purified enzyme consisted of a monomer with a subunit apparent Mr of 51 kDa when analyzed by SDS–PAGE. This polypeptide was recognized by affinity-purified anti- Solanum chacoense Bitt. recombinant HK IgGs. The protein was digested with trypsin and its digestion products were subjected to MS – MS sequencing after HPLC separation. The sequences of these tryptic peptides matched the predicted coding sequence of the S. tuberosum HK1 gene with a coverage of 57%. Examination of the kinetic properties of the purified protein HK1 indicates that it may be regulated by the internal O2 concentration of the tuber because of its sensitivity to acidic pHs and inhibition by ADP.

Complete genome sequence of Crohn's disease-associated adherent-invasive E. coli strain LF82.

BackgroundIleal lesions of Crohn's disease (CD) patients are abnormally colonized by pathogenic adherent-invasive Escherichia coli (AIEC) able to invade and to replicate within intestinal epithelial cells and macrophages.Principal FindingsWe report here the complete genome sequence of E. coli LF82, the reference strain of adherent-invasive E. coli associated with ileal Crohn's disease. The LF82 genome of 4,881,487 bp total size contains a circular chromosome with a size of 4,773,108 bp and a plasmid of 108,379 bp. The analysis of predicted coding sequences (CDSs) within the LF82 flexible genome indicated that this genome is close to the avian pathogenic strain APEC_01, meningitis-associated strain S88 and urinary-isolated strain UTI89 with regards to flexible genome and single nucleotide polymorphisms in various virulence factors. Interestingly, we observed that strains LF82 and UTI89 adhered at a similar level to Intestine-407 cells and that like LF82, APEC_01 and UTI89 were highly invasive. However, A1EC strain LF82 had an intermediate killer phenotype compared to APEC-01 and UTI89 and the LF82 genome does not harbour most of specific virulence genes from ExPEC. LF82 genome has evolved from those of ExPEC B2 strains by the acquisition of Salmonella and Yersinia isolated or clustered genes or CDSs located on pLF82 plasmid and at various loci on the chromosome.ConclusionLF82 genome analysis indicated that a number of genes, gene clusters and pathoadaptative mutations which have been acquired may play a role in virulence of AIEC strain LF82.

Genome Analysis of Moraxella catarrhalis Strain RH4, a Human Respiratory Tract Pathogen

Moraxella catarrhalis is an emerging human-restricted respiratory tract pathogen that is a common cause of childhood otitis media and exacerbations of chronic obstructive pulmonary disease in adults. Here, we report the first completely assembled and annotated genome sequence of an isolate of M. catarrhalis, strain RH4, which originally was isolated from blood of an infected patient. The RH4 genome consists of 1,863,286 nucleotides that form 1,886 protein-encoding genes. Comparison of the RH4 genome to the ATCC 43617 contigs demonstrated that the gene content of both strains is highly conserved. In silico phylogenetic analyses based on both 16S rRNA and multilocus sequence typing revealed that RH4 belongs to the seroresistant lineage. We were able to identify almost the entire repertoire of known M. catarrhalis virulence factors and mapped the members of the biosynthetic pathways for lipooligosaccharide, peptidoglycan, and type IV pili. Reconstruction of the central metabolic pathways suggested that RH4 relies on fatty acid and acetate metabolism, as the genes encoding the enzymes required for the glyoxylate pathway, the tricarboxylic acid cycle, the gluconeogenic pathway, the nonoxidative branch of the pentose phosphate pathway, the beta-oxidation pathway of fatty acids, and acetate metabolism were present. Moreover, pathways important for survival under challenging in vivo conditions, such as the iron-acquisition pathways, nitrogen metabolism, and oxidative stress responses, were identified. Finally, we showed by microarray expression profiling that approximately 88% of the predicted coding sequences are transcribed under in vitro conditions. Overall, these results provide a foundation for future research into the mechanisms of M. catarrhalis pathogenesis and vaccine development.

Comparative genomic analysis of 1047 completely sequenced cDNAs from an Arabidopsis-related model halophyte, Thellungiella halophila

BackgroundThellungiella halophila (also known as T. salsuginea) is a model halophyte with a small size, short life cycle, and small genome. Thellungiella genes exhibit a high degree of sequence identity with Arabidopsis genes (90% at the cDNA level). We previously generated a full-length enriched cDNA library of T. halophila from various tissues and from whole plants treated with salinity, chilling, freezing stress, or ABA. We determined the DNA sequences of 20 000 cDNAs at both the 5'- and 3' ends, and identified 9569 distinct genes.ResultsHere, we completely sequenced 1047 Thellungiella full-length cDNAs representing abiotic-stress-related genes, transcription factor genes, and protein phosphatase 2C genes. The predicted coding sequences, 5'-UTRs, and 3'-UTRs were compared with those of orthologous genes from Arabidopsis for length, sequence similarity, and structure. The 5'-UTR sequences of Thellungiella and Arabidopsis orthologs shared a significant level of similarity, although the motifs were rearranged. While examining the stress-related Thellungiella coding sequences, we found a short splicing variant of T. halophila salt overly sensitive 1 (ThSOS1), designated ThSOS1S. ThSOS1S contains the transmembrane domain of ThSOS1 but lacks the C-terminal hydrophilic region. The expression level of ThSOS1S under normal growth conditions was higher than that of ThSOS1. We also compared the expression levels of Na+-transport-system genes between Thellungiella and Arabidopsis by using full-length cDNAs from each species as probes. Several genes that play essential roles in Na+ excretion, compartmentation, and diffusion (SOS1, SOS2, NHX1, and HKT1) were expressed at higher levels in Thellungiella than in Arabidopsis.ConclusionsThe full-length cDNA sequences obtained in this study will be essential for the ongoing annotation of the Thellungiella genome, especially for further improvement of gene prediction. Moreover, they will enable us to find splicing variants such as ThSOS1S (AB562331).

Genomic Analysis of an Attenuated Chlamydia abortus Live Vaccine Strain Reveals Defects in Central Metabolism and Surface Proteins

Comparative genomic analysis of a wild-type strain of the ovine pathogen Chlamydia abortus and its nitrosoguanidine-induced, temperature-sensitive, virulence-attenuated live vaccine derivative identified 22 single nucleotide polymorphisms unique to the mutant, including nine nonsynonymous mutations, one leading to a truncation of pmpG, which encodes a polymorphic membrane protein, and two intergenic mutations potentially affecting promoter sequences. Other nonsynonymous mutations mapped to a pmpG pseudogene and to predicted coding sequences encoding a putative lipoprotein, a sigma-54-dependent response regulator, a PhoH-like protein, a putative export protein, two tRNA synthetases, and a putative serine hydroxymethyltransferase. One of the intergenic mutations putatively affects transcription of two divergent genes encoding pyruvate kinase and a putative SOS response nuclease, respectively. These observations suggest that the temperature-sensitive phenotype and associated virulence attenuation of the vaccine strain result from disrupted metabolic activity due to altered pyruvate kinase expression and/or alteration in the function of one or more membrane proteins, most notably PmpG and a putative lipoprotein.

Co-evolution of genomes and plasmids within Chlamydia trachomatis and the emergence in Sweden of a new variant strain.

BackgroundChlamydia trachomatis is the most common cause of sexually transmitted infections globally and the leading cause of preventable blindness in the developing world. There are two biovariants of C. trachomatis: 'trachoma', causing ocular and genital tract infections, and the invasive 'lymphogranuloma venereum' strains. Recently, a new variant of the genital tract C. trachomatis emerged in Sweden. This variant escaped routine diagnostic tests because it carries a plasmid with a deletion. Failure to detect this strain has meant it has spread rapidly across the country provoking a worldwide alert. In addition to being a key diagnostic target, the plasmid has been linked to chlamydial virulence. Analysis of chlamydial plasmids and their cognate chromosomes was undertaken to provide insights into the evolutionary relationship between chromosome and plasmid. This is essential knowledge if the plasmid is to be continued to be relied on as a key diagnostic marker, and for an understanding of the evolution of Chlamydia trachomatis.ResultsThe genomes of two new C. trachomatis strains were sequenced, together with plasmids from six C. trachomatis isolates, including the new variant strain from Sweden. The plasmid from the new Swedish variant has a 377 bp deletion in the first predicted coding sequence, abolishing the site used for PCR detection, resulting in negative diagnosis. In addition, the variant plasmid has a 44 bp duplication downstream of the deletion. The region containing the second predicted coding sequence is the most highly conserved region of the plasmids investigated. Phylogenetic analysis of the plasmids and chromosomes are fully congruent. Moreover this analysis also shows that ocular and genital strains diverged from a common C. trachomatis progenitor.ConclusionThe evolutionary pathways of the chlamydial genome and plasmid imply that inheritance of the plasmid is tightly linked with its cognate chromosome. These data suggest that the plasmid is not a highly mobile genetic element and does not transfer readily between isolates. Comparative analysis of the plasmid sequences has revealed the most conserved regions that should be used to design future plasmid based nucleic acid amplification tests, to avoid diagnostic failures.

Identification of SFR6, a key component in cold acclimation acting post‐translationally on CBF function

The sfr6-1 mutant of Arabidopsis thaliana was identified previously on the basis of its failure to undergo acclimation to freezing temperatures following exposure to low positive temperatures. This failure is attributed to a defect in the pathway leading to cold on-regulated (COR) gene expression via CBF (C-box binding factor) transcription factors. We identified a region of chromosome 4 containing SFR6 by positional mapping. Fine mapping of the sfr6-1 mutation proved impossible as the locus resides very close to the centromere. Therefore, we screened 380 T-DNA lines with insertions in genes within the large region to which sfr6-1 mapped. This resulted in the identification of two further mutant alleles of SFR6 (sfr6-2 and sfr6-3); like the original sfr6-1 mutation, these disrupt freezing tolerance and COR gene expression. To determine the protein sequence, we cloned an SFR6 cDNA based on the predicted coding sequence, but this offered no indication as to the mechanism by which SFR6 acts. The SFR6 gene itself is not strongly regulated by cold, thus discounting regulation of SFR6 activity at the transcriptional level. We show that over-expression of CBF1 or CBF2 transcription factors, which constitutively activate COR genes in the wild-type, cannot do so in sfr6-1. We demonstrate that CBF protein accumulates to wild-type levels in response to cold in sfr6-1. These results indicate a role for the SFR6 protein in the CBF pathway -downstream of CBF translation. The fact that the SFR6 protein is targeted to the nucleus may suggest a direct role in modulating gene expression.