Prediction of whole-genome DNA G + C content within the genus Aeromonas based on housekeeping gene sequences

Different methods are available to determine the G + C content (e.g. thermal denaturation temperature or high performance liquid chromatography, HPLC), but obtained values may differ significantly between strains, as well as between laboratories. Recently, several authors have demonstrated that the genomic DNA G + C content of prokaryotes can be reliably estimated from one or several protein coding gene nucleotide sequences. Few G + C content values have been published for the Aeromonas species described and the data, when available, are often incomplete or provide only a range of values. Our aim in this current work was twofold. First, the genomic G + C content of the type or reference strains of all species and subspecies of the genus Aeromonas was determined with a traditional experimental method in the same laboratory. Second, we wanted to see if the sequence-based method to estimate the G + C content described by Fournier et al. [7] could be applied to determine the G + C content of the different species of Aeromonas from the sequences of the genes used in taxonomy or phylogeny for this genus.     

A new RNASeq-based reference transcriptome for sugar beet and its application in transcriptome-scale analysis of vernalization and gibberellin responses

Background

Actinoplanes sp. SE50/110 is known as the wild type producer of the alpha-glucosidase inhibitor acarbose, a potent drug used worldwide in the treatment of type-2 diabetes mellitus. As the incidence of diabetes is rapidly rising worldwide, an ever increasing demand for diabetes drugs, such as acarbose, needs to be anticipated. Consequently, derived Actinoplanes strains with increased acarbose yields are being used in large scale industrial batch fermentation since 1990 and were continuously optimized by conventional mutagenesis and screening experiments. This strategy reached its limits and is generally superseded by modern genetic engineering approaches. As a prerequisite for targeted genetic modifications, the complete genome sequence of the organism has to be known.

Results

Here, we present the complete genome sequence of Actinoplanes sp. SE50/110 [GenBank:CP003170], the first publicly available genome of the genus Actinoplanes, comprising various producers of pharmaceutically and economically important secondary metabolites. The genome features a high mean G + C content of 71.32% and consists of one circular chromosome with a size of 9,239,851 bp hosting 8,270 predicted protein coding sequences. Phylogenetic analysis of the core genome revealed a rather distant relation to other sequenced species of the family Micromonosporaceae whereas Actinoplanes utahensis was found to be the closest species based on 16S rRNA gene sequence comparison. Besides the already published acarbose biosynthetic gene cluster sequence, several new non-ribosomal peptide synthetase-, polyketide synthase- and hybrid-clusters were identified on the Actinoplanes genome. Another key feature of the genome represents the discovery of a functional actinomycete integrative and conjugative element.

Conclusions

The complete genome sequence of Actinoplanes sp. SE50/110 marks an important step towards the rational genetic optimization of the acarbose production. In this regard, the identified actinomycete integrative and conjugative element could play a central role by providing the basis for the development of a genetic transformation system for Actinoplanes sp. SE50/110 and other Actinoplanes spp. Furthermore, the identified non-ribosomal peptide synthetase- and polyketide synthase-clusters potentially encode new antibiotics and/or other bioactive compounds, which might be of pharmacologic interest.     

Genista tinctoria microsymbionts from Poland are new members of Bradyrhizobium japonicum bv. genistearum

The phylogeny and taxonomic position of slow-growing Genista tinctoria rhizobia from Poland, Ukraine and England were estimated by comparative 16S rDNA, atpD, and dnaK sequence analyses, PCR-RFLP of 16S rDNA, DNA G + C content, and DNA–DNA hybridization. Each core gene studied placed the G. tinctoria rhizobia in the genus Bradyrhizobium cluster with unequivocal bootstrap support. G. tinctoria symbionts and bradyrhizobial strains shared 96–99% similarity in 16S rDNA sequences. Their similarity for atpD and dnaK sequences was 93–99% and 89–99%, respectively. These data clearly showed that G. tinctoria rhizobia belonged to the genus Bradyrhizobium. 16S rDNA sequence analysis was in good agreement with the results of the PCR-RFLP of the 16S rRNA gene. Although the tested strains formed separate lineages to the reference bradyrhizobia their RFLP 16S rDNA patterns were quite similar. The genomic DNA G + C content of three G. tinctoria rhizobia was in the range from 60.64 to 62.83 mol%. Data for species identification were obtained from DNA–DNA hybridization experiments. G. tinctoria microsymbionts from Poland were classified within Bradyrhizobium japonicum genomospecies based on 56–82% DNA–DNA similarity.     

6-10× pyrosequencing is a practical approach for whole prokaryote genome studies

Next generation 454 pyrosequencing technology for whole bacterial genome sequencing involves a deep sequencing strategy with at least 15-20 × in depth proposed by official protocols but usually done with over 20 × in practices. In this study, we carried out a comprehensive evaluation of quality of the de novo assemblies based on realistic pyrosequencing simulated data from 1480 prokaryote genomes and 7 runs of machine-generated data. Our results demonstrated that for most of the prokaryote genomes, 6-10 × sequencing in qualified runs with 400 bp reads could produce high quality draft assembly (> 98% genome coverage, < 100 contigs with N50 size > 100 kb, single base accuracy > 99.99, indel error rate < 0.01%, false gene loss/duplication rate < 0.5%). Our study proves the power of low depth pyrosequencing strategy, which provides a cost-effective way for sequencing whole prokaryote genomes in a short time and enables further studies in microbial population diversity and comparative genomics.     

Crystal structures of the bacterial solute receptor AcbH displaying an exclusive substrate preference for β-D-galactopyranose

Solute receptors (binding proteins) are indispensable components of canonical ATP-binding cassette importers in prokaryotes. Here, we report on the characterization and crystal structures in the closed and open conformations of AcbH, the solute receptor of the putative carbohydrate transporter AcbFG which is encoded in the acarbose (acarviosyl-1,4-maltose) biosynthetic gene cluster from Actinoplanes sp. SE50/110. Binding assays identified AcbH as a high-affinity monosaccharide-binding protein with a dissociation constant (K_d) for β-d-galactopyranose of 9.8 ± 1.0 nM. Neither galactose-containing di- and trisaccharides, such as lactose and raffinose, nor monosaccharides including d-galacturonic acid, l-arabinose, d-xylose and l-rhamnose competed with [¹⁴C]galactose for binding to AcbH. Moreover, AcbH does not bind d-glucose, which is a common property of all but one d-galactose-binding proteins characterized to date. Strikingly, determination of the X-ray structure revealed that AcbH is structurally homologous to maltose-binding proteins rather than to glucose-binding proteins. Two helices are inserted in the substrate-binding pocket, which reduces the cavity size and allows the exclusive binding of monosaccharides, specifically β-d-galactopyranose, in the ⁴C₁ conformation. Site-directed mutagenesis of three residues from the binding pocket (Arg82, Asp361 and Arg362) that interact with the axially oriented O4-H hydroxyl of the bound galactopyranose and subsequent functional analysis indicated that these residues are crucial for galactose binding. To our knowledge, this is the first report of the tertiary structure of a solute receptor with exclusive affinity for β-d-galactopyranose. The putative role of a galactose import system in the context of acarbose metabolism in Actinoplanes sp. is discussed.     

Prediction of pathogenicity islands in Enterohemorrhagic Escherichia coli O157:H7 using genomic barcodes

The genome of lethal animal pathogenic bacterium Enterohemorrhagic Escherichia coli (EHEC) O157:H7 is characterized by the presence of multiple pathogenicity islands (PAIs). Computational methods have been developed to identify PAIs based on the distinguishing G + C levels in some PAI versus non-PAI regions. We observed that PAIs can have a very similar G + C level to that of the host chromosome, which may have led to false negative predictions using these methods. We have applied a novel method of genomic barcodes to identify PAIs. Using this technique, we have successfully identified both known and novel PAIs in the genomes of three strains of EHEC O157:H7.     

Species-genomic relationships among the tribasic diploid and polyploid Carthamus taxa based on physical mapping of active and inactive 18S–5.8S–26S and 5S ribosomal RNA gene families,and the two tandemly repeated DNA sequences

In the genus Carthamus (2n = 20, 22, 24, 44, 64; x = 10, 11, 12), most of the homologues within and between the chromosome complements are difficult to be identified. In the present work, we used fluorescent in situ hybridisation (FISH) to determine the chromosome distribution of the two rRNA gene families, and the two isolated repeated DNA sequences in the 14 Carthamus taxa. The distinctive variability in the distribution, number and signal intensity of hybridisation sites for 18S–26S and 5S rDNA loci could generally distinguish the 14 Carthamus taxa. Active 18S–26S rDNA sites were generally associated with NOR loci on the nucleolar chromosomes. The two A genome taxa, C. glaucus ssp. anatolicus and C. boissieri with 2n = 20, and the two botanical varieties of B genome C. tinctorius (2n = 24) had diagnostic FISH patterns. The present results support the origin of C. tinctorius from C. palaestinus. FISH patterns of C. arborescens vis-à-vis the other taxa indicate a clear division of Carthamus taxa into two distinct lineages. Comparative distribution and intensity pattern of 18S–26S rDNA sites could distinguish each of the tetraploid and hexaploid taxa. The present results indicate that C. boissieri (2n = 20) is one of the genome donors for C. lanatus and C. lanatus ssp. lanatus (2n = 44), and C. lanatus is one of the progenitors for the hexaploid (2n = 64) taxa. The association of pCtKpnI-2 repeated sequence with rRNA gene cluster (orphon) in 2–10 nucleolar and non-nucleolar chromosomes and the consistent occurrence of pCtKpnI-1 repeated sequence at the subtelomeric region in all the taxa analysed indicate some functional role of these sequences.     

Next‐generation DNA barcoding: using next‐generation sequencing to enhance and accelerate DNA barcode capture from single specimens

DNA barcoding is an efficient method to identify specimens and to detect undescribed/cryptic species. Sanger sequencing of individual specimens is the standard approach in generating large‐scale DNA barcode libraries and identifying unknowns. However, the Sanger sequencing technology is, in some respects, inferior to next‐generation sequencers, which are capable of producing millions of sequence reads simultaneously. Additionally, direct Sanger sequencing of DNA barcode amplicons, as practiced in most DNA barcoding procedures, is hampered by the need for relatively high‐target amplicon yield, coamplification of nuclear mitochondrial pseudogenes, confusion with sequences from intracellular endosymbiotic bacteria (e.g. Wolbachia) and instances of intraindividual variability (i.e. heteroplasmy). Any of these situations can lead to failed Sanger sequencing attempts or ambiguity of the generated DNA barcodes. Here, we demonstrate the potential application of next‐generation sequencing platforms for parallel acquisition of DNA barcode sequences from hundreds of specimens simultaneously. To facilitate retrieval of sequences obtained from individual specimens, we tag individual specimens during PCR amplification using unique 10‐mer oligonucleotides attached to DNA barcoding PCR primers. We employ 454 pyrosequencing to recover full‐length DNA barcodes of 190 specimens using 12.5% capacity of a 454 sequencing run (i.e. two lanes of a 16 lane run). We obtained an average of 143 sequence reads for each individual specimen. The sequences produced are full‐length DNA barcodes for all but one of the included specimens. In a subset of samples, we also detected Wolbachia, nontarget species, and heteroplasmic sequences. Next‐generation sequencing is of great value because of its protocol simplicity, greatly reduced cost per barcode read, faster throughout and added information content.     

Potential use of DNA barcoding for the identification of Salvia based on cpDNA and nrDNA sequences

An effective DNA marker for authenticating the genus Salvia was screened using seven DNA regions (rbcL, matK, trnL–F, and psbA–trnH from the chloroplast genome, and ITS, ITS1, and ITS2 from the nuclear genome) and three combinations (rbcL + matK, psbA–trnH + ITS1, and trnL–F + ITS1). The present study collected 232 sequences from 27 Salvia species through DNA sequencing and 77 sequences within the same taxa from the GenBank. The discriminatory capabilities of these regions were evaluated in terms of PCR amplification success, intraspecific and interspecific divergence, DNA barcoding gaps, and identification efficiency via a tree-based method. ITS1 was superior to the other marker for discriminating between species, with an accuracy of 81.48%. The three combinations did not increase species discrimination. Finally, we found that ITS1 is a powerful barcode for identifying Salvia species, especially Salvia miltiorrhiza.     

Use of an adenosine triphosphate assay, and simultaneous staining with fluorescein diacetate and propidium iodide, to evaluate the effects of cryoprotectants on hard coral (Echinopora spp.) oocytes

The objective was to examine the effects of cryoprotectants on oocytes of hard corals (Echinopora spp.) to obtain basic knowledge for cryopreservation procedures. Oocytes were exposed to various concentrations of cryoprotectants (0.25 to 5.0 M) for 20 min at room temperature (25 °C). Two tests were used to assess ovarian follicle viability: fluorescein diacetate (FDA) + propidium iodide (PI) staining, and adenosine triphosphate (ATP) assay. Both FDA + PI staining and ATP assay indicated that cryoprotectant toxicity to oocytes increased in the order methanol, dimethyl sulfoxide (DMSO), propylene glycol (PG), and ethylene glycol (EG). The no observed effect concentrations for Echinopora spp. oocytes were 1.0, 0.5, 0.25, and 0.25 M for methanol, DMSO, PG, and EG, respectively, when assessed with FDA + PI. The ATP assay was more sensitive than FDA + PI staining (P < 0.05). Oocyte viability after 1.0 M methanol, DMSO, EG, or PG treatment for 20 min at room temperature assessed with FDA + PI tests and ATP assay were 88.9 ± 3.1% and 72.2 ± 4.4%, 66.2 ± 5.0% and 23.2 ± 4.9%, 58.9 ± 5.4% and 1.1 ± 0.7%, and 49.1 ± 5.1% and 0.9 ± 0.5%, respectively. We inferred that the ATP assay was a valuable measure of cellular injury after cryoprotectant incubation. The results of this study provided a basis for development of protocols to cryopreserve coral oocytes.     

Influence of interleukin-6 gene − 174G>C polymorphism on development of atherosclerosis: A meta-analysis of 50 studies involving 33,514 subjects

Increasing epidemiological studies have focused on the associations between interleukin-6 (IL-6) gene − 174G>C polymorphism and atherosclerotic diseases, but the results are still controversial. This meta-analysis was designed to identify whether this association exists. PubMed, Embase, Web of Science, Cochrane database, Clinicaltrials.gov and Current Controlled Trials, Chinese Clinical Trial Registry, CBMdisc, CNKI and Google Scholar were searched to get the genetic association studies. The crude odds ratios (ORs) and their corresponding 95% confidence intervals (CIs) were used to estimate the association between the IL-6 gene − 174G>C polymorphism and atherosclerosis ( AS ) risk. The subgroup analyses were made on the following: ethnicity, atherosclerotic diseases and source of controls. Finally, 50 studies (15,029 cases and 18,485 controls) were included in this meta-analysis. Overall, no significant association was found between the IL-6 gene − 174G>C polymorphism and AS risk (for C allele vs. G allele: OR = 1.02, 95% CI = 0.94–1.11, p = 0.64; for C/C vs. G/G: OR = 1.01, 95% CI = 0.85–1.21, p = 0.88; for C/C vs. C/G + G/G: OR = 0.97, 95% CI = 0.84–1.12, p = 0.68; for C/C + C/G vs. G/G: OR = 1.07, 95% CI = 0.97–1.17, p = 0.18). In the subgroup analyses, significant associations were found between the IL-6 gene − 174G>C polymorphism and AS in non-Caucasian group (for CC + CG vs. GG: OR = 1.22, 95% CI = 1.06–1.41, p = 0.005), other atherosclerotic diseases group (for C allele vs. G allele: OR = 0.75, 95% CI = 0.61–0.93, p = 0.008; for C/C vs. G/G: OR = 0.56, 95% CI = 0.38–0.81, p = 0.002; for C/C vs. C/G + G/G: OR = 0.60, 95% CI = 0.45–0.79, p = 0.0004) and population-based group (for C allele vs. G allele: OR = 1.09, 95% CI = 1.00–1.18, p = 0.04; for CC + CG vs. GG: OR = 1.15, 95% CI = 1.04–1.27, p = 0.005). In summary, the present meta-analysis suggests that the IL-6 gene − 174G C polymorphism is associated with the susceptibility to AS. However, due to the high heterogeneity in the meta-analysis, the results should be interpreted with caution.     

Hydrogen atom abstraction from HOOOCl by chlorine atom and OH radical

Gas-phase reactions of HOOOCl with both Cl atom and OH radical are investigated using ab initio methods. The structures of all reactants, products, intermediates, and transition states have been optimized and characterized with the quadratic configuration interaction (QCISD) method. The overall mechanism for the Cl + HOOOCl and OH + HOOOCl reaction is the formation of HCl + O₂ + ClO and H₂O + O₂ + ClO, respectively. The rate-limiting step in each reaction is the abstraction of hydrogen from HOOOCl by either Cl or OH radicals and the barrier height is predicted to be 1.9 kcal mol⁻¹ and 8.1 kcal mol⁻¹ for abstraction by Cl atom and OH radical, respectively. Since both barriers for hydrogen abstraction are high, the reaction is suggested to be slow. These results also suggest that an atmospheric removal mechanism for HOOOCl may result from reaction with Cl atoms rather than with OH radicals, and that photolysis of HOOOCl may be the major removal mechanism for the intermediate.     

Determining the repertoire of immunodominant proteins via whole-genome amplification of intracellular pathogens

Culturing many obligate intracellular bacteria is difficult or impossible. However, these organisms have numerous adaptations allowing for infection persistence and immune system evasion, making them some of the most interesting to study. Recent advancements in genome sequencing, pyrosequencing and Phi29 amplification, have allowed for examination of whole-genome sequences of intracellular bacteria without culture. We have applied both techniques to the model obligate intracellular pathogen Anaplasma marginale and the human pathogen Anaplasma phagocytophilum, in order to examine the ability of phi29 amplification to determine the sequence of genes allowing for immune system evasion and long-term persistence in the host. When compared to traditional pyrosequencing, phi29-mediated genome amplification had similar genome coverage, with no additional gaps in coverage. Additionally, all msp2 functional pseudogenes from two strains of A. marginale were detected and extracted from the phi29-amplified genomes, highlighting its utility in determining the full complement of genes involved in immune evasion.     

The analysis of microsatellites and compound microsatellites in 56 complete genomes of Herpesvirales

Simple sequence repeats (SSRs), or microsatellites, are special DNA/RNA sequences with repeated unit of 1–6 bp. The genomes of Herpesvirales have many repeating structures, which is an excellent system to study the evolution and roles of microsatellites and compound microsatellites in viruses. Therefore, 56 genomes of Herpesvirales were selected and the occurrence, composition and complexity of different repeats were investigated in the genomes. A total of 63,939 microsatellites and 5825 compound microsatellites were extracted from 56 genomes. It found that GC content has a significant strong correlation with both the counts of microsatellites (CM) and the counts of compound microsatellites (CCM). However, genome size has a moderate correlation only with CM and almost no correlation with CCM. The compound microsatellites occurring in genic regions are obviously more than that in intergenic regions. In general, the number of compound microsatellite decreases with the increase of complexity (C) (the count of individual microsatellites being part of a compound microsatellite) and the complexity hardly exceeds C = 4. The vast majority of compound microsatellites exist in intergenic regions, when C ≥ 10. The distributions of SSRs tend to be organism-specific rather than host-specific in herpesvirus genomes. The diversity of microsatellites and compound microsatellites may be helpful for a better understanding of the viral genetic diversity, genotyping, and evolutionary biology in herpesviruses genomes.     

Simplified process for ethanol production from sugarcane bagasse using hydrolysate-resistant Escherichia coli strain MM160

Hexose and pentose sugars from phosphoric acid pretreated sugarcane bagasse were co-fermented to ethanol in a single vessel (SScF), eliminating process steps for solid-liquid separation and sugar cleanup. An initial liquefaction step (L) with cellulase was included to improve mixing and saccharification (L + SScF), analogous to a corn ethanol process. Fermentation was enabled by the development of a hydrolysate-resistant mutant of Escherichia coli LY180, designated MM160. Strain MM160 was more resistant than the parent to inhibitors (furfural, 5-hydroxymethylfurfural, and acetate) formed during pretreatment. Bagasse slurries containing 10% and 14% dry weight (fiber plus solubles) were tested using pretreatment temperatures of 160-190 °C (1% phosphoric acid, 10 min). Enzymatic saccharification and inhibitor production both increased with pretreatment temperature. The highest titer (30 g/L ethanol) and yield (0.21 g ethanol/g bagasse dry weight) were obtained after incubation for 122 h using 14% dry weight slurries of pretreated bagasse (180 °C).     

BIPES,a cost-effective high-throughput method for assessing microbial diversity

Pyrosequencing of 16S rRNA (16S) variable tags has become the most popular method for assessing microbial diversity, but the method remains costly for the evaluation of large numbers of environmental samples with high sequencing depths. We developed a barcoded Illumina paired-end (PE) sequencing (BIPES) method that sequences each 16S V6 tag from both ends on the Illumina HiSeq 2000, and the PE reads are then overlapped to obtain the V6 tag. The average accuracy of Illumina single-end (SE) reads was only 97.9%, which decreased from ∼99.9% at the start of the read to less than 85% at the end of the read; nevertheless, overlapping of the PE reads significantly increased the sequencing accuracy to 99.65% by verifying the 3′ end of each SE in which the sequencing quality was degraded. After the removal of tags with two or more mismatches within the medial 40–70 bases of the reads and of tags with any primer errors, the overall base sequencing accuracy of the BIPES reads was further increased to 99.93%. The BIPES reads reflected the amounts of the various tags in the initial template, but long tags and high GC tags were underestimated. The BIPES method yields 20–50 times more 16S V6 tags than does pyrosequencing in a single-flow cell run, and each of the BIPES reads costs less than 1/40 of a pyrosequencing read. As a laborsaving and cost-effective method, BIPES can be routinely used to analyze the microbial ecology of both environmental and human microbiomes.     

PyroClean: Denoising Pyrosequences from Protein-Coding Amplicons for the Recovery of Interspecific and Intraspecific Genetic Variation

High-throughput parallel sequencing is a powerful tool for the quantification of microbial diversity through the amplification of nuclear ribosomal gene regions. Recent work has extended this approach to the quantification of diversity within otherwise difficult-to-study metazoan groups. However, nuclear ribosomal genes present both analytical challenges and practical limitations that are a consequence of the mutational properties of nuclear ribosomal genes. Here we exploit useful properties of protein-coding genes for cross-species amplification and denoising of 454 flowgrams. We first use experimental mixtures of species from the class Collembola to amplify and pyrosequence the 5′ region of the COI barcode, and we implement a new algorithm called PyroClean for the denoising of Roche GS FLX pyrosequences. Using parameter values from the analysis of experimental mixtures, we then analyse two communities sampled from field sites on the island of Tenerife. Cross-species amplification success of target mitochondrial sequences in experimental species mixtures is high; however, there is little relationship between template DNA concentrations and pyrosequencing read abundance. Homopolymer error correction and filtering against a consensus reference sequence reduced the volume of unique sequences to approximately 5% of the original unique raw reads. Filtering of remaining non-target sequences attributed to PCR error, sequencing error, or numts further reduced unique sequence volume to 0.8% of the original raw reads. PyroClean reduces or eliminates the need for an additional, time-consuming step to cluster reads into Operational Taxonomic Units, which facilitates the detection of intraspecific DNA sequence variation. PyroCleaned sequence data from field sites in Tenerife demonstrate the utility of our approach for quantifying evolutionary diversity and its spatial structure. Comparison of our sequence data to public databases reveals that we are able to successfully recover both interspecific and intraspecific sequence diversity.     

Jumping and gliding rodents: Mitogenomic affinities of Pedetidae and Anomaluridae deduced from an RNA-Seq approach

An RNA-Seq strategy was used to obtain the complete set of protein-coding mitochondrial genes from two rodent taxa. Thanks to the next generation sequencing (NGS) 454 approach, we determined the complete mitochondrial DNA genome from Graphiurus kelleni (Mammalia: Rodentia: Gliridae) and partial mitogenome from Pedetes capensis (Pedetidae), and compared them with published rodent and outgroup mitogenomes. We finished the mitogenome sequencing by a series of amplicons using conserved PCR primers to fill the gaps corresponding to tRNA, rRNA and control regions. Phylogenetic analyses of the mitogenomes suggest a well-supported rodent phylogeny in agreement with nuclear gene trees. Pedetes groups with Anomalurus into the clade Anomaluromorpha, while Graphiurus branches within the squirrel-related clade. Moreover, Pedetes + Anomalurus branch with Castor into the mouse-related clade. Our study demonstrates the utility of NGS for obtaining new mitochondrial genomes as well as the importance of choosing adequate models of sequence evolution to infer the phylogeny of rodents.