Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs

Comparative genomics of the major parasitic worms

Parasitic nematodes (roundworms) and platyhelminths (flatworms) cause debilitating chronic infections of humans and animals, decimate crop production and are a major impediment to socioeconomic development. Here we report a broad comparative study of 81 genomes of parasitic and non-parasitic worms. We have identified gene family births and hundreds of expanded gene families at key nodes in the phylogeny that are relevant to parasitism. Examples include gene families that modulate host immune responses, enable parasite migration though host tissues or allow the parasite to feed. We reveal extensive lineage-specific differences in core metabolism and protein families historically targeted for drug development. From an in silico screen, we have identified and prioritized new potential drug targets and compounds for testing. This comparative genomics resource provides a much-needed boost for the research community to understand and combat parasitic worms.

Complete Genome Sequence of Lactococcus lactis subsp. cremoris 3107, Host for the Model Lactococcal P335 Bacteriophage TP901-1

The complete genome sequence of Lactococcus lactis subsp. cremoris 3107, a dairy starter strain and a host for the model lactococcal P335 bacteriophage TP901-1, is reported here. The circular chromosome of L. lactis subsp. cremoris 3107 is among the smallest genomes of currently sequenced lactococcal strains. L. lactis subsp. cremoris 3107 harbors a complement of six plasmids, which appears to be a reflection of its adaptation to the nutrient-rich dairy environment.

Recent polyploidization events in three Saccharum founding species

The complexity of polyploid Saccharum genomes hindered progress of genome research and crop improvement in sugarcane. To understand their genome structure, transcriptomes of 59 F₁ individuals derived from S. officinarumLA Purple and S. robustum Molokai 5829 (2n = 80, x = 10 for both) were sequenced, yielding 11 157 and 8998 SNPs and 83 and 105 linkage groups, respectively. Most markers in each linkage group aligned to single sorghum chromosome. However, 71 interchromosomal rearrangements were detected between sorghum and S. officinarum or S. robustum, and 24 (33.8%) of them were shared between S. officinarum and S. robustum, indicating their occurrence before the speciation event that separated these two species. More than 2000 gene pairs from S. spontaneum, S. officinarum and S. robustum were analysed to estimate their divergence time. Saccharum officinarum and S. robustum diverged about 385 thousand years ago, and the whole-genome duplication events occurred after the speciation event because of shared interchromosomal rearrangements. The ancestor of these two species diverged from S. spontaneum about 769 thousand years ago, and the reduction in basic chromosome number from 10 to 8 in S. spontaneum occurred after the speciation event but before the two rounds of whole-genome duplication. Our results proved that S. officinarum is a legitimate species in its own right and not a selection from S. robustum during the domestication process in the past 10 000 years. Our findings rejected a long-standing hypothesis and clarified the timing of speciation and whole-genome duplication events in Saccharum.

Evolutionary Analysis of Plastid Genomes of Seven Lonicera L. Species: Implications for Sequence Divergence and Phylogenetic Relationships

Plant plastomes play crucial roles in species evolution and phylogenetic reconstruction studies due to being maternally inherited and due to the moderate evolutionary rate of genomes. However, patterns of sequence divergence and molecular evolution of the plastid genomes in the horticulturally- and economically-important Lonicera L. species are poorly understood. In this study, we collected the complete plastomes of seven Lonicera species and determined the various repeat sequence variations and protein sequence evolution by comparative genomic analysis. A total of 498 repeats were identified in plastid genomes, which included tandem (130), dispersed (277), and palindromic (91) types of repeat variations. Simple sequence repeat (SSR) elements analysis indicated the enriched SSRs in seven genomes to be mononucleotides, followed by tetra-nucleotides, dinucleotides, tri-nucleotides, hex-nucleotides, and penta-nucleotides. We identified 18 divergence hotspot regions (rps15, rps16, rps18, rpl23, psaJ, infA, ycf1, trnN-GUU-ndhF, rpoC2-rpoC1, rbcL-psaI, trnI-CAU-ycf2, psbZ-trnG-UCC, trnK-UUU-rps16, infA-rps8, rpl14-rpl16, trnV-GAC-rrn16, trnL-UAA intron, and rps12-clpP) that could be used as the potential molecular genetic markers for the further study of population genetics and phylogenetic evolution of Lonicera species. We found that a large number of repeat sequences were distributed in the divergence hotspots of plastid genomes. Interestingly, 16 genes were determined under positive selection, which included four genes for the subunits of ribosome proteins (rps7, rpl2, rpl16, and rpl22), three genes for the subunits of photosystem proteins (psaJ, psbC, and ycf4), three NADH oxidoreductase genes (ndhB, ndhH, and ndhK), two subunits of ATP genes (atpA and atpB), and four other genes (infA, rbcL, ycf1, and ycf2). Phylogenetic analysis based on the whole plastome demonstrated that the seven Lonicera species form a highly-supported monophyletic clade. The availability of these plastid genomes provides important genetic information for further species identification and biological research on Lonicera.

Advances in genotyping microsatellite markers through sequencing and consequences of scoring methods for Ceratonia siliqua (Leguminosae)

PREMISE OF THE STUDY

Simple sequence repeat (SSR) or microsatellite markers have been used in a broad range of studies mostly scoring alleles on the basis of amplicon size as a proxy for the number of repeat units of an SSR motif. However, additional sources of variation within the SSR or in the flanking regions have largely remained undetected.

METHODS

In this study, we implemented a next-generation sequencing-based genotyping approach in a newly characterized set of 18 nuclear SSR markers for the carob tree, Ceratonia siliqua. Our aim was to evaluate the effect of three different methods of scoring molecular variation present within microsatellite markers on the genetic diversity and structure results.

RESULTS

The analysis of the sequences of 77 multilocus genotypes from four populations revealed SSR variation and additional sources of polymorphism in 87% of the loci analyzed (42 single-nucleotide polymorphisms and five insertion/deletion polymorphisms), as well as divergent paralog copies in two loci. Ignoring sequence variation under standard amplicon size genotyping resulted in incorrect identification of 69% of the alleles, with important effects on the genetic diversity and structure estimates.

DISCUSSION

Next-generation sequencing allows the detection and scoring of SSRs, single-nucleotide polymorphisms, and insertion/deletion polymorphisms to increase the resolution of population genetic studies.

Oropouche orthobunyavirus: Genetic characterization of full-length genomes and development of molecular methods to discriminate natural reassortments

Oropouche orthobunyavirus (OROV) has significant impact in public health in Amazon region. This arbovirus is one of the most common causes of febrile illness in Brazil, and is responsible for several epidemics since 1960's. In this study, we sequenced and characterized the complete coding sequences (S-, M- and L-RNA) of 35 OROV isolates from Brazil. Here, we classified 20 strains in genotype I from Pará and Maranhão states, nine as genotype II from Pará and Rondônia states confirmed, four classified into genotype III from Acre, Maranhão, Minas Gerais and Rondônia states and two genotype IV from Amazonas State. Also, we did not observe reassortment events involving the OROV isolates. In addition, we developed novel RT-PCR tools to identify reassortment events among OROV strains. These data will be useful to better understand the molecular epidemiology and diagnostic of OROV infections.

TransFlow: a modular framework for assembling and assessing accurate de novo transcriptomes in non-model organisms

BACKGROUND

The advances in high-throughput sequencing technologies are allowing more and more de novo assembling of transcriptomes from many new organisms. Some degree of automation and evaluation is required to warrant reproducibility, repetitivity and the selection of the best possible transcriptome. Workflows and pipelines are becoming an absolute requirement for such a purpose, but the issue of assembling evaluation for de novo transcriptomes in organisms lacking a sequenced genome remains unsolved. An automated, reproducible and flexible framework called TransFlow to accomplish this task is described.

RESULTS

TransFlow with its five independent modules was designed to build different workflows depending on the nature of the original reads. This architecture enables different combinations of Illumina and Roche/454 sequencing data, and can be extended to other sequencing platforms. Its capabilities are illustrated with the selection of reliable plant reference transcriptomes and the assembling six transcriptomes (three case studies for grapevine leaves, olive tree pollen, and chestnut stem, and other three for haustorium, epiphytic structures and their combination for the phytopathogenic fungus Podosphaera xanthii). Arabidopsis and poplar transcriptomes revealed to be the best references. A common result regarding de novo assemblies is that Illumina paired-end reads of 100 nt in length assembled with OASES can provide reliable transcriptomes, while the contribution of longer reads is noticeable only when they complement a set of short, single-reads.

CONCLUSIONS

TransFlow can handle up to 181 different assembling strategies. Evaluation based on principal component analyses allows its self-adaptation to different sets of reads to provide a suitable transcriptome for each combination of reads and assemblers. As a result, each case study has its own behaviour, prioritises evaluation parameters, and gives an objective and automated way for detecting the best transcriptome within a pool of them. Sequencing data type and quantity (preferably several hundred millions of 2×100 nt or longer), assemblers (OASES for Illumina, MIRA4 and EULER-SR reconciled with CAP3 for Roche/454) and strategy (preferably scaffolding with OASES, and probably merging with Roche/454 when available) arise as the most impacting factors.

De novo assembly of wheat root transcriptomes and transcriptional signature of longitudinal differentiation

Hidden underground, root systems constitute an important part of the plant for its development, nourishment and sensing the soil environment around it, but we know very little about its genetic regulation in crop plants like wheat. In the present study, we de novo assembled the root transcriptomes in reference cultivar Chinese Spring from RNA-seq reads generated by the 454-GS-FLX and HiSeq platforms. The FLX reads were assembled into 24,986 transcripts with completeness of 54.84%, and the HiSeq reads were assembled into 91,543 high-confidence protein-coding transcripts, 2,404 low-confidence protein-coding transcripts, and 13,181 non-coding transcripts with the completeness of >90%. Combining the FLX and HiSeq assemblies, we assembled a root transcriptome of 92,335 ORF-containing transcripts. Approximately 7% of the coding transcripts and ~2% non-coding transcripts are not present in the current wheat genome assembly. Functional annotation of both assemblies showed similar gene ontology patterns and that ~7% coding and >5% non-coding transcripts are root-specific. Transcription quantification identified 1,728 differentially expressed transcripts between root tips and maturation zone, and functional annotation of these transcripts captured a transcriptional signature of longitudinal development of wheat root. With the transcriptomic resources developed, this study provided the first view of wheat root transcriptome under different developmental zones and laid a foundation for molecular studies of wheat root development and growth using a reverse genetic approach.

Rapid Divergence of Genome Architectures Following the Origin of an Ectomycorrhizal Symbiosis in the Genus Amanita

Fungi are evolutionary shape shifters and adapt quickly to new environments. Ectomycorrhizal (EM) symbioses are mutualistic associations between fungi and plants and have evolved repeatedly and independently across the fungal tree of life, suggesting lineages frequently reconfigure genome content to take advantage of open ecological niches. To date analyses of genomic mechanisms facilitating EM symbioses have involved comparisons of distantly related species, but here, we use the genomes of three EM and two asymbiotic (AS) fungi from the genus Amanita as well as an AS outgroup to study genome evolution following a single origin of symbiosis. Our aim was to identify the defining features of EM genomes, but our analyses suggest no clear differentiation of genome size, gene repertoire size, or transposable element content between EM and AS species. Phylogenetic inference of gene gains and losses suggests the transition to symbiosis was dominated by the loss of plant cell wall decomposition genes, a confirmation of previous findings. However, the same dynamic defines the AS species A. inopinata, suggesting loss is not strictly associated with origin of symbiosis. Gene expansions in the common ancestor of EM Amanita were modest, but lineage specific and large gene family expansions are found in two of the three EM extant species. Even closely related EM genomes appear to share few common features. The genetic toolkit required for symbiosis appears already encoded in the genomes of saprotrophic species, and this dynamic may explain the pervasive, recurrent evolution of ectomycorrhizal associations.

De novo transcriptome based on next-generation sequencing reveals candidate genes with sex-specific expression in Arapaima gigas (Schinz, 1822), an ancient Amazonian freshwater fish

Background

The Arapaima (Arapaima gigas) is one of the world's largest freshwater bony fish, and is found in the rivers of the Amazon basin. This species is a potential aquaculture resource, although reproductive management in captivity is limited in particular due to the lack of external sexual dimorphism. In this study, using the 454 Roche platform (pyrosequencing) techniques, we evaluated a major portion of the transcriptome of this important Amazonian species.

Results

Four libraries obtained from the liver and skin tissue of juvenile specimens (representing males and females separately) were sequenced, yielding 5,453,919 high-quality reads. The de novo transcriptome assembly resulted in 175,792 contigs, with 51,057 significant blast hits. A total of 38,586 transcripts were mapped by Gene Ontology using Blast2GO. We identified 20,219 genes in the total transcriptome (9,551 in the liver and 16,818 in the skin). The gene expression analyses indicated 105 genes in the liver and 204 in the skin with differentiated expression profiles, with 95 being over-expressed in the females and 214 in the males. The log2 Fold Change and heatmap based on Reads Per Kilobase per Million mapped reads (RPKM) revealed that the gene expression in the skin is highly differentiated between male and female arapaima, while the levels of expression in the liver are similar between the sexes.

Conclusion

Transcriptome analysis based on pyrosequencing proved to be a reliable tool for the identification of genes with differentiated expression profiles between male and female arapaima. These results provide useful insights into the molecular pathways of sexual dimorphism in this important Amazonian species, and for comparative analyses with other teleosts.

Distribution and Evolution of Peroxisomes in Alveolates (Apicomplexa, Dinoflagellates, Ciliates)

The peroxisome was the last organelle to be discovered and five decades later it is still the Cinderella of eukaryotic compartments. Peroxisomes have a crucial role in the detoxification of reactive oxygen species, the beta-oxidation of fatty acids, and the biosynthesis of etherphospholipids, and they are assumed to be present in virtually all aerobic eukaryotes. Apicomplexan parasites including the malaria and toxoplasmosis agents were described as the first group of mitochondriate protists devoid of peroxisomes. This study was initiated to reassess the distribution and evolution of peroxisomes in the superensemble Alveolata (apicomplexans, dinoflagellates, ciliates). We established transcriptome data from two chromerid algae (Chromera velia, Vitrella brassicaformis), and two dinoflagellates (Prorocentrum minimum, Perkinsus olseni) and identified the complete set of essential peroxins in all four reference species. Our comparative genome analysis provides unequivocal evidence for the presence of peroxisomes in Toxoplasma gondii and related genera. Our working hypothesis of a common peroxisomal origin of all alveolates is supported by phylogenetic analyses of essential markers such as the import receptor Pex5. Vitrella harbors the most comprehensive set of peroxisomal proteins including the catalase and the glyoxylate cycle and it is thus a promising model organism to investigate the functional role of this organelle in Apicomplexa.

Metagenomics of Coral Reefs Under Phase Shift and High Hydrodynamics

Local and global stressors have affected coral reef ecosystems worldwide. Switches from coral to algal dominance states and microbialization are the major processes underlying the global decline of coral reefs. However, most of the knowledge concerning microbialization has not considered physical disturbances (e.g., typhoons, waves, and currents). Southern Japan reef systems have developed under extreme physical disturbances. Here, we present analyses of a three-year investigation on the coral reefs of Ishigaki Island that comprised benthic and fish surveys, water quality analyses, metagenomics and microbial abundance data. At the four studied sites, inorganic nutrient concentrations were high and exceeded eutrophication thresholds. The dissolved organic carbon (DOC) concentration (up to 233.3 μM) and microbial abundance (up to 2.5 × 105 cell/mL) values were relatively high. The highest vibrio counts coincided with the highest turf cover (∼55-85%) and the lowest coral cover (∼4.4-10.2%) and fish biomass (0.06 individuals/m2). Microbiome compositions were similar among all sites and were dominated by heterotrophs. Our data suggest that a synergic effect among several regional stressors are driving coral decline. In a high hydrodynamics reef environment, high algal/turf cover, stimulated by eutrophication and low fish abundance due to overfishing, promote microbialization. Together with crown-of-thorns starfish (COTS) outbreaks and possible of climate changes impacts, theses coral reefs are likely to collapse.

Comparative Analyses of Chloroplast Genomes of Cucurbitaceae Species: Lights into Selective Pressures and Phylogenetic Relationships

Cucurbitaceae is the fourth most important economic plant family with creeping herbaceous species mainly distributed in tropical and subtropical regions. Here, we described and compared the complete chloroplast genome sequences of ten representative species from Cucurbitaceae. The lengths of the ten complete chloroplast genomes ranged from 155,293 bp (C. sativus) to 158,844 bp (M. charantia), and they shared the most common genomic features. 618 repeats of three categories and 813 microsatellites were found. Sequence divergence analysis showed that the coding and IR regions were highly conserved. Three protein-coding genes (accD, clpP, and matK) were under selection and their coding proteins often have functions in chloroplast protein synthesis, gene transcription, energy transformation, and plant development. An unconventional translation initiation codon of psbL gene was found and provided evidence for RNA editing. Applying BI and ML methods, phylogenetic analysis strongly supported the position of Gomphogyne, Hemsleya, and Gynostemma as the relatively original lineage in Cucurbitaceae. This study suggested that the complete chloroplast genome sequences were useful for phylogenetic studies. It would also determine potential molecular markers and candidate DNA barcodes for coming studies and enrich the valuable complete chloroplast genome resources of Cucurbitaceae.

PdumBase: a transcriptome database and research tool for Platynereis dumerilii and early development of other metazoans

BACKGROUND

The marine polychaete annelid Platynereis dumerilii has recently emerged as a prominent organism for the study of development, evolution, stem cells, regeneration, marine ecology, chronobiology and neurobiology within metazoans. Its phylogenetic position within the spiralian/ lophotrochozoan clade, the comparatively high conservation of ancestral features in the Platynereis genome, and experimental access to any stage within its life cycle, make Platynereis an important model for elucidating the complex regulatory and functional molecular mechanisms governing early development, later organogenesis, and various features of its larval and adult life. High resolution RNA-seq gene expression data obtained from specific developmental stages can be used to dissect early developmental mechanisms. However, the potential for discovery of these mechanisms relies on tools to search, retrieve, and compare genome-wide information within Platynereis, and across other metazoan taxa.

RESULTS

To facilitate exploration and discovery by the broader scientific community, we have developed a web-based, searchable online research tool, PdumBase, featuring the first comprehensive transcriptome database for Platynereis dumerilii during early stages of development (2 h ~ 14 h). Our database also includes additional stages over the P. dumerilii life cycle and provides access to the expression data of 17,213 genes (31,806 transcripts) along with annotation information sourced from Swiss-Prot, Gene Ontology, KEGG pathways, Pfam domains, TmHMM, SingleP, and EggNOG orthology. Expression data for each gene includes the stage, the normalized FPKM, the raw read counts, and information that can be leveraged for statistical analyses of differential gene expression and the construction of genome-wide co-expression networks. In addition, PdumBase offers early stage transcriptome expression data from five further species as a valuable resource for investigators interested in comparing early development in different organisms. To understand conservation of Platynereis gene models and to validate gene annotation, most Platynereis gene models include a comprehensive phylogenetic analysis across 18 species representing diverse metazoan taxa.

CONCLUSIONS

PdumBase represents the first online resource for the early developmental transcriptome of Platynereis dumerilii. It serves as a research platform for discovery and exploration of gene expression during early stages, throughout the Platynereis life cycle, and enables comparison to other model organisms. PdumBase is freely available at http://pdumbase.gdcb.iastate.edu .

Gomesin inhibits melanoma growth by manipulating key signaling cascades that control cell death and proliferation

Consistent with their diverse pharmacology, peptides derived from venomous animals have been developed as drugs to treat disorders as diverse as hypertension, diabetes and chronic pain. Melanoma has a poor prognosis due in part to its metastatic capacity, warranting further development of novel targeted therapies. This prompted us to examine the anti-melanoma activity of the spider peptides gomesin (AgGom) and a gomesin-like homolog (HiGom). AgGom and HiGom dose-dependently reduced the viability and proliferation of melanoma cells whereas it had no deleterious effects on non-transformed neonatal foreskin fibroblasts. Concordantly, gomesin-treated melanoma cells showed a reduced G0/G1 cell population. AgGom and HiGom compromised proliferation of melanoma cells via activation of the p53/p21 cell cycle check-point axis and the Hippo signaling cascade, together with attenuation of the MAP kinase pathway. We show that both gomesin peptides exhibit antitumoral activity in melanoma AVATAR-zebrafish xenograft tumors and that HiGom also reduces tumour progression in a melanoma xenograft mouse model. Taken together, our data highlight the potential of gomesin for development as a novel melanoma-targeted therapy.

The complete nucleotide sequence of chloroplast genome of Euphorbia kansui (Euphorbiaceae), an endemic herb in China

Euphorbia kansui T.N. Liou ex S.B. Ho (Euphorbiaceae) is a perennial herb plant endemic to China. This species has important economic and medicinal values. In this study, we first characterized the complete nucleotide sequence of chloroplast (cp) genome of E. kansui using the Illumina Hiseq platform. The cp genome was 161,061 bp in length, comprising of a large single copy (LSC) region of 91,288 bp, a small single copy (SSC) region of 17,086 bp, and two inverted repeat regions of 26,343 bp each. The cp genome contains 130 genes, including 86 protein-coding genes, 8 ribosomal RNAs (rRNAs), and 36 transfer RNAs (tRNAs). The phylogenetic analysis indicated that E. kansui was placed as a sister to the congeneric Euphorbia esula.

Assembly and Analysis of Unmapped Genome Sequence Reads Reveal Novel Sequence and Variation in Dogs

Dogs are excellent animal models for human disease. They have extensive veterinary histories, pedigrees, and a unique genetic system due to breeding practices. Despite these advantages, one factor limiting their usefulness is the canine genome reference (CGR) which was assembled using a single purebred Boxer. Although a common practice, this results in many high-quality reads remaining unmapped. To address this whole-genome sequence data from three breeds, Border Collie (n = 26), Bearded Collie (n = 7), and Entlebucher Sennenhund (n = 8), were analyzed to identify novel, non-CGR genomic contigs using the previously validated pseudo-de novo assembly pipeline. We identified 256,957 novel contigs and paired-end relationships together with BLAT scores provided 126,555 (49%) high-quality contigs with genomic coordinates containing 4.6 Mb of novel sequence absent from the CGR. These contigs close 12,503 known gaps, including 2.4 Mb containing partially missing sequences for 11.5% of Ensembl, 16.4% of RefSeq and 12.2% of canFam3.1+ CGR annotated genes and 1,748 unmapped contigs containing 2,366 novel gene variants. Examples for six disease-associated genes (SCARF2, RD3, COL9A3, FAM161A, RASGRP1 and DLX6) containing gaps or alternate splice variants missing from the CGR are also presented. These findings from non-reference breeds support the need for improvement of the current Boxer-only CGR to avoid missing important biological information. The inclusion of the missing gene sequences into the CGR will facilitate identification of putative disease mutations across diverse breeds and phenotypes.

Metagenomics reveals functional synergy and novel polysaccharide utilization loci in the Castor canadensis fecal microbiome

The North American beaver (Castor canadensis) has long been considered an engineering marvel, transforming landscapes and shaping biological diversity through its dam building behavior. While the beaver possesses conspicuous morphological features uniquely adapted for the use of woody plants as construction materials and dietary staples, relatively little is known about the specialized microorganisms inhabiting the beaver gastrointestinal tract and their functional roles in determining host nutrition. Here we use a combination of shotgun metagenomics, functional screening and carbohydrate biochemistry to chart the community structure and metabolic power of the beaver fecal microbiome. We relate this information to the metabolic capacity of other wood feeding and hindgut fermenting organisms and profile the functional repertoire of glycoside hydrolase (GH) families distributed among and between population genome bins. Metagenomic screening revealed novel mechanisms of xylan oligomer degradation involving GH43 enzymes from uncharacterized subfamilies and divergent polysaccharide utilization loci, indicating the potential for synergistic biomass deconstruction. Together, these results open a functional metagenomic window on less conspicuous adaptations enabling the beaver microbiome to efficiently convert woody plants into host nutrition and point toward rational design of enhanced enzyme mixtures for biorefining process streams.

The complete mitochondrial genome of Xizicus (Haploxizicus) maculatus revealed by Next-Generation Sequencing and phylogenetic implication (Orthoptera, Meconematinae)

Xizicus Gorochov, 1993, the quiet-calling katydid, is a diverse genus with 68 species in world, which includes more than 45 species in China, has undergone numerous taxonomic revisions with contradicting conclusions. In this study the complete mitochondrial genome of Xizicus (Haploxizicus) maculatus collected from Hainan for the first time was sequenced using the Next-Generation Sequencing (NGS) technology. The length of whole mitogenome is 16,358 bp and contains the typical gene arrangement, base composition, and codon usage found in other related species. The overall base composition of the mitochondrial genome is 37.0 % A, 32.2 % T, 20.2 % C, and 10.6 % G. All 13 protein-coding genes (PCGs) began with typical ATN initiation codon. Nine of the 13 PCGs have a complete termination codon, but the remaining four genes (COI, COIII, ND5, and ND4) terminate with an incomplete T. Phylogenetic analyses are carried out based on the concatenated dataset of 13 PCGs and two rRNAs of Tettigoniidae species available in GenBank. Both Bayesian inference and Maximum Likelihood analyses recovered each subfamily as a monophyletic group. Regardless of the position of Lipotactinae, the relationships among the subfamilies of Tettigoniidae were as follows: ((((Tettigoniinae, Bradyporinae) Meconematinae) Conocephalinae) Hexacentrinae). The topological structure of the phylogeny trees showed that the Xizicus (Haploxizicus) maculatus is closer to Xizicus (Xizicus) fascipes than Xizicus (Eoxizicus) howardi.

Complete mitochondrial genome of Xiphidiopsis (Xiphidiopsis) gurneyi (Orthoptera, Tettigoniidae, Meconematinae)

Xiphidiopsis (Xiphidiopsis) gurneyi belongs to Meconematinae. The complete mitochondrial genome of X. (X.) gurneyi was sequenced by the next-generation sequencing (NGS) technology. The total length of the mitogenome was 16,225 bp and contains the typical gene arrangement, base composition, codon usage found in Meconematinae species. Phylogenetic tree was constructed based on concatenated datasets of PCGs and rRNAs of X. (X.) gurneyi and 19 Tettigoniidae species to assess their phylogenic relationship. Phylogenetic analysis showed that X. (X.) gurneyi was more closely related to the genus of Xizicus.

Characterization of the whole plastid genome sequence of Abies chensiensis (Pinaceae), an endangered endemic conifer in China

Abies chensiensis Van Tiegh. (Pinaceae) is a rare and endangered endemic conifer in China. In this study, using the Illumina sequencing platform, we firstly characterized its whole plastid genome sequence. Our study revealed that A. chensiensis have a typical plastid genome of 121,498 bp in length, comprised of a large single copy region of 76,484 bp, a small single copy region of 42,654 bp and two inverted repeat regions of 1180 bp. A total of 113 genes, 74 protein-coding genes, 35 tRNA, and 4 rRNA genes were identified. The phylogenetic analysis indicated that A. chensiensis was placed as a sister to the congeneric A. sibirica.

Identification and manipulation of a novel locus to improve cell tolerance to short-chain alcohols in Escherichia coli

Escherichia coli KO11 is a popular ethanologenic strain, but is more sensitive to ethanol than other producers. Here, an ethanol-tolerant mutant EM was isolated from ultraviolet mutagenesis library of KO11. Comparative genomic analysis added by piecewise knockout strategy and complementation assay revealed EKO11_3023 (espA) within the 36.6-kb deletion from KO11 was the only locus responsible for ethanol sensitivity. Interestingly, when espA was deleted in strain W (the parent strain of KO11), ethanol tolerance was dramatically elevated to the level of espA-free hosts [e.g., MG1655 and BL21(DE3)]. And overexpression of espA in strains MG1655 and BL21(DE3) led to significantly enhanced ethanol sensitivity. In addition to ethanol, deletion of espA also improved cell tolerance to other short-chain (C2–C4) alcohols, including methanol, isopropanol, n-butanol, isobutanol and 2-butanol. Therefore, espA was responsible for short-chain alcohol sensitivity of W-strains compared to other cells, which provides a potential engineering target for alcohols production.

Characterization of the cork formation and production transcriptome in Quercus cerris × suber hybrids

Cork oak is the main cork-producing species worldwide, and plays a significant economic, ecological and social role in the Mediterranean countries, in particular in Portugal and Spain. The ability to produce cork is limited to a few species, hence it must involve specific regulation mechanisms that are unique to these species. However, to date, these mechanisms remain largely understudied, especially with approaches involving the use of high-throughput sequencing technology. In this study, the transcriptome of cork-producing and non-cork-producing Quercus cerris × suber hybrids was analyzed in order to elucidate the differences between the two groups of trees displaying contrasting phenotypes for cork production. The results revealed the presence of a significant number of genes exclusively associated with cork production, in the trees that developed cork. Moreover, several gene ontology subcategories, such as cell wall biogenesis, lipid metabolic processes, metal ion binding and apoplast/cell wall, were only detected in the trees with cork production. These results indicate the existence, at the transcriptome level, of mechanisms that seem to be unique and necessary for cork production, which is an advancement in our knowledge regarding the genetic regulation behind cork formation and production.

Recovery of Five Complete Influenza A(H1N1)pdm09 Genome Sequences from the 2015 Influenza Outbreak in India by Metagenomic Sequencing

Five complete (H1N1)pdm09 viral sequences were recovered from hospitalized individuals during the 2015 influenza outbreak by metagenomic sequencing. Four of the genomes are from oropharyngeal swabs, and one is from an isolate.

Five complete (H1N1)pdm09 viral sequences were recovered from hospitalized individuals during the 2015 influenza outbreak by metagenomic sequencing. Four of the genomes are from oropharyngeal swabs, and one is from an isolate. All five sequences belong to an emerging 6B clade. Studying them further is critical for outbreak preparedness.

Metagenomic insights into the microbial diversity in manganese-contaminated mine tailings and their role in biogeochemical cycling of manganese

To extend the knowledge on the microbial diversity of manganese rich environments, we performed a clone library based study using metagenomic approach. Pyrosequencing based analysis of 16S rRNA genes were carried out on an Illumina platform to gain insights into the bacterial community inhabiting in a manganese mining site and the taxonomic profiles were correlated with the inherent capacities of these strains to solubilise manganese. The application of shot gun sequencing in this study yielded results which revealed the highest prevalence of Proteobacteria (42.47%), followed by Actinobacteria (23.99%) in the area of study. Cluster of orthologous group (COG) functional category has 85,066 predicted functions. Out of which 11% are involved in metabolism of amino acid, 9% are involved in production and conversion of energy while Keto Encyclopedia of Gene and Genomes (KEGG) functional category has 107,388 predicted functions, out of which 55% are involved in cellular metabolism, 15% are environmental and information processing and 12% are genetic information processing in nature. The isolated microbial consortia demonstrated visible growth in presence of high concentrations of Mn. Solubilisation studies resulted in 86% of manganese recovery after 20 days. The result presented in this study has important implications in understanding the microbial diversity in manganese contaminated mine tailings and their role in natural geochemical cycling of Mn.

Comparison of Compression-Based Measures with Application to the Evolution of Primate Genomes

An efficient DNA compressor furnishes an approximation to measure and compare information quantities present in, between and across DNA sequences, regardless of the characteristics of the sources. In this paper, we compare directly two information measures, the Normalized Compression Distance (NCD) and the Normalized Relative Compression (NRC). These measures answer different questions; the NCD measures how similar both strings are (in terms of information content) and the NRC (which, in general, is nonsymmetric) indicates the fraction of one of them that cannot be constructed using information from the other one. This leads to the problem of finding out which measure (or question) is more suitable for the answer we need. For computing both, we use a state of the art DNA sequence compressor that we benchmark with some top compressors in different compression modes. Then, we apply the compressor on DNA sequences with different scales and natures, first using synthetic sequences and then on real DNA sequences. The last include mitochondrial DNA (mtDNA), messenger RNA (mRNA) and genomic DNA (gDNA) of seven primates. We provide several insights into evolutionary acceleration rates at different scales, namely, the observation and confirmation across the whole genomes of a higher variation rate of the mtDNA relative to the gDNA. We also show the importance of relative compression for localizing similar information regions using mtDNA.

The genomes of three Bradyrhizobium sp. isolated from root nodules of Lupinus albescens grown in extremely poor soils display important genes for resistance to environmental stress

Lupinus albescens is a resistant cover plant that establishes symbiotic relationships with bacteria belonging to the Bradyrhizobium genus. This symbiosis helps the development of these plants in adverse environmental conditions, such as the ones found in arenized areas of Southern Brazil. This work studied three Bradyrhizobium sp. (AS23, NAS80 and NAS96) isolated from L. albescens plants that grow in extremely poor soils (arenized areas and adjacent grasslands). The genomes of these three strains were sequenced in the Ion Torrent platform using the IonXpress library preparation kit, and presented a total number of bases of 1,230,460,823 for AS23, 1,320,104,022 for NAS80, and 1,236,105,093 for NAS96. The genome comparison with closest strains Bradyrhizobium japonicum USDA6 and Bradyrhizobium diazoefficiens USDA110 showed important variable regions (with less than 80% of similarity). Genes encoding for factors for resistance/tolerance to heavy metal, flagellar motility, response to osmotic and oxidative stresses, heat shock proteins (present only in the three sequenced genomes) could be responsible for the ability of these microorganisms to survive in inhospitable environments. Knowledge about these genomes will provide a foundation for future development of an inoculant bioproduct that should optimize the recovery of degraded soils using cover crops.

Draft Genome Sequence of a Fungus (Fusarium tricinctum) Cultured from a Monoisolate Native to the Himalayas

Here, we report the draft de novo genome sequence assembly of Fusarium tricinctum (strain T6), using IonTorrent sequencing chemistry and an Ion 530 chip ExT kit for sequencing. The genome assembly resulted in 42,732,204 bp from a total 6.62 Gb, with a median read length of 386 bp.

Microbial communities in natural rubber coagula during maturation: impacts on technological properties of dry natural rubber

AIM

To characterize microbial communities present in natural rubber (NR) coagula from Hevea brasiliensis latex during maturation and identify microbial taxa (bacteria and fungi) having an impact on dry NR properties.

METHODS AND RESULTS

Microbial community dynamics in NR coagula maturated under controlled conditions were compared and related with the evolution of dry NR properties. The pyrosequencing of 16S (119 837 effective reads) and 18S (131 879 effective reads) rRNA gene regions was performed on 21 samples covering different maturation times and two aeration conditions. Results showed a relatively high bacterial richness (Chao1 estimates of 200-1000) associated with significant bacterial dynamics. Lactic acid bacteria (LAB) were dominant in the first days of maturation. Then, in aerobic conditions, development of Actinobacteria represented by the family Microbacteriaceae was associated with alkalinization of the samples and a higher sensitivity of NR to thermo-oxidation as evaluated by its plasticity retention index (PRI). In anaerobiosis, the reduced development of bacteria, mostly LAB present, was associated with improved NR properties (higher initial plasticity P₀ and PRI).

CONCLUSIONS

The involvement of micro-organisms in the evolution of dry NR properties during the maturation of NR coagula was confirmed. The importance of the structure and dynamics of microbial communities is specifically highlighted.

SIGNIFICANCE AND IMPACT OF THE STUDY

Natural rubber is a key elastomer for the tyre industry and for a variety of other applications. The majority of raw NR is obtained by natural coagulation of H. brasiliensis latex under the activity of micro-organisms. An improved understanding of the microbial communities involved in the maturation of NR coagula may lead to an improvement in the production process of raw NR to provide a better consistency in NR quality.

Complete assembly of a dengue virus type 3 genome from a recent genotype III clade by metagenomic sequencing of serum

Background: Mosquito-borne flaviviruses, such as dengue and Japanese encephalitis virus (JEV), cause life-threatening diseases, particularly in the tropics. Methods: Here we performed unbiased metagenomic sequencing of RNA extracted from the serum of four patients and the plasma of one patient, all hospitalized at a tertiary care centre in South India with severe or prolonged febrile illness, together with the serum from one healthy control, in 2014. Results: We identified and assembled a complete dengue virus type 3 sequence from a case of severe dengue fever. We also identified a small number of JEV sequences in the serum of two adults with febrile illness, including one with severe dengue. Phylogenetic analysis revealed that the dengue sequence belonged to genotype III. It has an estimated divergence time of 13.86 years from the most highly related Indian strains. In total, 11 amino acid substitutions were predicted for this strain in the antigenic envelope protein, when compared to the parent strain used for development of the first commercial dengue vaccine.  Conclusions: We demonstrate that both genome assembly and detection of a low number of viral sequences are possible through the unbiased sequencing of clinical material. These methods may help ascertain causal agents for febrile illnesses with no known cause.

The complete chloroplast genome sequence of Echinocodon Lobophyllus (Campanulaceae) revealed by next-generation sequencing and phylogenetic implication

The complete chloroplast genome sequence of Echinocodon lobophyllus was determined by Illumina pair-end sequencing. The results showed that the complete plastid genome was 169,419 bp in length, containing a large single copy (LSC) of 85,599 bp and a small single copy (SSC) of 8054 bp, which were separated by a pair of 37,883 bp inverted repeats (IRs). A total of 109 unique genes were annotated, including 75 protein coding genes, 30 tRNA genes, and 4 rRNA genes. Among these genes, 16 genes contained one or two introns. The overall GC contents of the plastid genome were 38.4%. A maximum likelihood phylogenetic analysis showed that E. lobophyllus and Codonopsis minima formed one clade at the base of the phylogenetic tree of Campanulaceae with a high support value.

Comparative Plastid Genomes of Primula Species: Sequence Divergence and Phylogenetic Relationships

Compared to traditional DNA markers, genome-scale datasets can provide mass information to effectively address historically difficult phylogenies. Primula is the largest genus in the family Primulaceae, with members distributed mainly throughout temperate and arctic areas of the Northern Hemisphere. The phylogenetic relationships among Primula taxa still maintain unresolved, mainly due to intra- and interspecific morphological variation, which was caused by frequent hybridization and introgression. In this study, we sequenced and assembled four complete plastid genomes (Primula handeliana, Primula woodwardii, Primula knuthiana, and Androsace laxa) by Illumina paired-end sequencing. A total of 10 Primula species (including 7 published plastid genomes) were analyzed to investigate the plastid genome sequence divergence and their inferences for the phylogeny of Primula. The 10 Primula plastid genomes were similar in terms of their gene content and order, GC content, and codon usage, but slightly different in the number of the repeat. Moderate sequence divergence was observed among Primula plastid genomes. Phylogenetic analysis strongly supported that Primula was monophyletic and more closely related to Androsace in the Primulaceae family. The phylogenetic relationships among the 10 Primula species showed that the placement of P. knuthiana–P. veris clade was uncertain in the phylogenetic tree. This study indicated that plastid genome data were highly effective to investigate the phylogeny.

Different Selected Mechanisms Attenuated the Inhibitory Interaction of KIR2DL1 with C2+ HLA-C in Two Indigenous Human Populations in Southern Africa

The functions of human NK cells in defense against pathogens and placental development during reproduction are modulated by interactions of killer cell Ig-like receptors (KIRs) with HLA-A, -B and -C class I ligands. Both receptors and ligands are highly polymorphic and exhibit extensive differences between human populations. Indigenous to southern Africa are the KhoeSan, the most ancient group of modern human populations, who have highest genomic diversity worldwide. We studied two KhoeSan populations, the Nama pastoralists and the ≠Khomani San hunter-gatherers. Comprehensive next-generation sequence analysis of HLA-A, -B, and -C and all KIR genes identified 248 different KIR and 137 HLA class I, which assort into ∼200 haplotypes for each gene family. All 74 Nama and 78 ≠Khomani San studied have different genotypes. Numerous novel KIR alleles were identified, including three arising by intergenic recombination. On average, KhoeSan individuals have seven to eight pairs of interacting KIR and HLA class I ligands, the highest diversity and divergence of polymorphic NK cell receptors and ligands observed to date. In this context of high genetic diversity, both the Nama and the ≠Khomani San have an unusually conserved, centromeric KIR haplotype that has arisen to high frequency and is different in the two KhoeSan populations. Distinguishing these haplotypes are independent mutations in KIR2DL1, which both prevent KIR2DL1 from functioning as an inhibitory receptor for C2⁺ HLA-C. The relatively high frequency of C2⁺ HLA-C in the Nama and the ≠Khomani San appears to have led to natural selection against strong inhibitory C2-specific KIR.

Genomic characterization and phylogenetic analysis of the novel Pseudomonas phage PPSC2

We isolated a Pseudomonas phage infecting Pseudomonas fluorescens SA1 separated from a soil sample collected in Sichuan Province, China. This phage, which we named PPSC2, has a genome that is composed of a 97,330-bp-long linear double-stranded DNA with 47.51% G+C content and 168 putative protein-coding genes. We identified 20 tRNA genes in the genome of PPSC2, and the tRNA GC content ranged from 44.2% to 58.4%. Phylogenetic and BLASTn analysis revealed that the Pseudomonas phage PPSC2 should be considered a new member of the family Myoviridae.

Genomic characterization of the novel Ralstonia phage RPSC1

In this paper, I describe the genomic characteristics of a Ralstonia phage infecting Ralstonia solanacearum. The Ralstonia phage RPSC1 was isolated from a soil sample collected in Sichuan Province, in southwestern China. The complete genome of RPSC1 is composed of a linear double-stranded DNA 39,628 bp in length, with G+C content of 61.55%, and 43 putative protein-coding genes. All the putative protein-coding genes were on the same strand. No tRNA-encoding genes were identified. Phylogenetic and comparative genomics analyses indicate that Ralstonia phage RPSC1 should be considered a new member of the family Podoviridae. The wide host range contributes to the potential of Ralstonia phage RPSC1 as a biocontrol agent.

Molecular Evolution of Chloroplast Genomes of Orchid Species: Insights into Phylogenetic Relationship and Adaptive Evolution

Orchidaceae is the 3rd largest family of angiosperms, an evolved young branch of monocotyledons. This family contains a number of economically-important horticulture and flowering plants. However, the limited availability of genomic information largely hindered the study of molecular evolution and phylogeny of Orchidaceae. In this study, we determined the evolutionary characteristics of whole chloroplast (cp) genomes and the phylogenetic relationships of the family Orchidaceae. We firstly characterized the cp genomes of four orchid species: Cremastra appendiculata, Calanthe davidii, Epipactis mairei, and Platanthera japonica. The size of the chloroplast genome ranged from 153,629 bp (C. davidi) to 160,427 bp (E. mairei). The gene order, GC content, and gene compositions are similar to those of other previously-reported angiosperms. We identified that the genes of ndhC, ndhI, and ndhK were lost in C. appendiculata, in that the ndh I gene was lost in P. japonica and E. mairei. In addition, the four types of repeats (forward, palindromic, reverse, and complement repeats) were examined in orchid species. E. mairei had the highest number of repeats (81), while C. davidii had the lowest number (57). The total number of Simple Sequence Repeats is at least 50 in C. davidii, and, at most, 78 in P. japonica. Interestingly, we identified 16 genes with positive selection sites (the psbH, petD, petL, rpl22, rpl32, rpoC1, rpoC2, rps12, rps15, rps16, accD, ccsA, rbcL, ycf1, ycf2, and ycf4 genes), which might play an important role in the orchid species’ adaptation to diverse environments. Additionally, 11 mutational hotspot regions were determined, including five non-coding regions (ndhB intron, ccsA-ndhD, rpl33-rps18, ndhE-ndhG, and ndhF-rpl32) and six coding regions (rps16, ndhC, rpl32, ndhI, ndhK, and ndhF). The phylogenetic analysis based on whole cp genomes showed that C. appendiculata was closely related to C. striata var. vreelandii, while C. davidii and C. triplicate formed a small monophyletic evolutionary clade with a high bootstrap support. In addition, five subfamilies of Orchidaceae, Apostasioideae, Cypripedioideae, Epidendroideae, Orchidoideae, and Vanilloideae, formed a nested evolutionary relationship in the phylogenetic tree. These results provide important insights into the adaptive evolution and phylogeny of Orchidaceae.

Mating-Induced Differential Peptidomics of Neuropeptides and Protein Hormones in Agrotis ipsilon Moths

In many insects, mating induces drastic changes in male and female responses to sex pheromones or host-plant odors. In the male moth Agrotis ipsilon, mating induces a transient inhibition of behavioral and neuronal responses to the female sex pheromone. As neuropeptides and peptide hormones regulate most behavioral processes, we hypothesize that they could be involved in this mating-dependent olfactory plasticity. Here we used next-generation RNA sequencing and a combination of liquid chromatography, matrix assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry, and direct tissue profiling to analyze the transcriptome and peptidome of different brain compartments in virgin and mated males and females of A. ipsilon. We identified 37 transcripts encoding putative neuropeptide precursors and 54 putative bioactive neuropeptides from 23 neuropeptide precursors (70 sequences in total, 25 neuropeptide precursors) in different areas of the central nervous system including the antennal lobes, the gnathal ganglion, and the corpora cardiaca-corpora allata complex. Comparisons between virgin and mated males and females revealed tissue-specific differences in peptide composition between sexes and according to physiological state. Mated males showed postmating differences in neuropeptide occurrence, which could participate in the mating-induced olfactory plasticity.

Characterization of the Gamboa Virus Serogroup (Orthobunyavirus Genus, Peribunyaviridae Family)

Comprehensive comparative phylogenetic analyses were performed on 17 Gamboa serogroup viruses (GAMSVs) from distinct geographic regions in the Americas and other representative members of the genus Orthobunyavirus (Peribunyaviridae), based on small (S), medium (M), and large (L) open reading frame full-length and partial sequences. Genome characterization showed that the GAMSVs divide into four clades or genotypes. The GAMSVs have a genetic organization similar to other orthobunyaviruses, except that they have a larger NSm protein than other orthobunyaviruses. A serosurvey for Gamboa virus antibodies was performed in plasma from birds, other wild animals, and humans living around the Tucuruí hydroelectric dam in Pará state, northern Brazil, a known focus of GAMSV activity. Newborn chicks (Gallus gallus domesticus) were experimentally infected with a GAMSV, and the pathogenesis is described. Histopathological changes were primarily in the lungs and liver. Also, a review of the ecology of the GAMSVs in the Americas is included. In sum, this study presents the genomic and evolutionary characterization of the Gamboa group and the potential model of pathogenesis, which would be helpful for diagnostic purposes, epidemiology, and immunopathogenesis studies.

Exploring the Genomic Traits of Non-toxigenic Vibrio parahaemolyticus Strains Isolated in Southern Chile

Vibrio parahaemolyticus is the leading cause of seafood-borne gastroenteritis worldwide. As reported in other countries, after the rise and fall of the pandemic strain in Chile, other post-pandemic strains have been associated with clinical cases, including strains lacking the major toxins TDH and TRH. Since the presence or absence of tdh and trh genes has been used for diagnostic purposes and as a proxy of the virulence of V. parahaemolyticus isolates, the understanding of virulence in V. parahaemolyticus strains lacking toxins is essential to detect these strains present in water and marine products to avoid possible food-borne infection. In this study, we characterized the genome of four environmental and two clinical non-toxigenic strains (tdh-, trh-, and T3SS2-). Using whole-genome sequencing, phylogenetic, and comparative genome analysis, we identified the core and pan-genome of V. parahaemolyticus of strains of southern Chile. The phylogenetic tree based on the core genome showed low genetic diversity but the analysis of the pan-genome revealed that all strains harbored genomic islands carrying diverse virulence and fitness factors or prophage-like elements that encode toxins like Zot and RTX. Interestingly, the three strains carrying Zot-like toxin have a different sequence, although the alignment showed some conserved areas with the zot sequence found in V. cholerae. In addition, we identified an unexpected diversity in the genetic architecture of the T3SS1 gene cluster and the presence of the T3SS2 gene cluster in a non-pandemic environmental strain. Our study sheds light on the diversity of V. parahaemolyticus strains from the southern Pacific which increases our current knowledge regarding the global diversity of this organism.

Plastid Genome Comparative and Phylogenetic Analyses of the Key Genera in Fagaceae: Highlighting the Effect of Codon Composition Bias in Phylogenetic Inference

Fagaceae is one of the largest and economically important taxa within Fagales. Considering the incongruence among inferences from plastid and nuclear genes in the previous Fagaceae phylogeny studies, we assess the performance of plastid phylogenomics in this complex family. We sequenced and assembled four complete plastid genomes (Fagus engleriana, Quercus spinosa, Quercus aquifolioides, and Quercus glauca) using reference-guided assembly approach. All of the other 12 published plastid genomes in Fagaceae were retrieved for genomic analyses (including repeats, sequence divergence and codon usage) and phylogenetic inference. The genomic analyses reveal that plastid genomes in Fagaceae are conserved. Comparing the phylogenetic relationships of the key genera in Fagaceae inferred from different codon positions and gene function datasets, we found that the first two codon sites dataset recovered nearly all relationships and received high support. Thus, the result suggested that codon composition bias had great influence on Fagaceae phylogenetic inference. Our study not only provides basic understanding of Fagaceae plastid genomes, but also illuminates the effectiveness of plastid phylogenomics in resolving relationships of this intractable family.

Current Bioinformatics resources in combating infectious diseases

Bioinformatics tools and techniques analyzing next-generation sequencing (NGS) data are increasingly used for the diagnosis and monitoring of infectious diseases. It is of interest to review the application of bioinformatics tools, commonly used databases and NGS data in clinical microbiology, focusing on molecular identification, genotypic, microbiome research, antimicrobial resistance analysis and detection of unknown disease-associated pathogens in clinical specimens. This review documents available bioinformatics resources and databases that are used by medical microbiology scientists and physicians to control emerging infectious pathogens.

Phylogenetic classification of six novel species belonging to the genus Bifidobacterium comprising Bifidobacterium anseris sp. nov., Bifidobacterium criceti sp. nov., Bifidobacterium imperatoris sp. nov., Bifidobacterium italicum sp. nov., Bifidobacterium margollesii sp. nov. and Bifidobacterium parmae sp. nov

Six Bifidobacterium strains, i.e., Goo31D, Ham19E, Rab10A, Tam1G, Uis4E and Uis1B, were isolated from domestic goose (Anser domesticus), European hamster (Cricetus cricetus), European rabbit (Oryctolagus cuniculus), emperor tamarin (Saguinus imperator) and pygmy marmoset (Callithrix pygmaea). Cells are Gram-positive, non-motile, non-sporulating, facultative anaerobic and fructose 6-phosphate phosphoketolase-positive. Phylogenetic analyses based on 16S rRNA, ITS-, multilocus- sequences and the core genome revealed that bifidobacterial strains Goo31D, Ham19E, Rab10A, Tam1G, Uis4E and Uis1B exhibit close phylogenetic relatedness with Bifidobacterium choerinum LMG 10510, Bifidobacterium hapali DSM 100202, Bifidobacterium saguini DSM 23967 and Bifidobacterium stellenboschense DSM 23968. Genotyping based on the genome sequence of the isolated strains combined with phenotypic analyses, clearly show that these strains are distinct from each of the type strains of the so far recognized Bifidobacterium species. Thus, Bifidobacterium anseris sp. nov. (Goo31D=LMG 30189^T=CCUG 70960^T), Bifidobacterium criceti sp. nov. (Ham19E=LMG 30188^T=CCUG 70962^T), Bifidobacterium imperatoris sp. nov. (Tam1G=LMG 30297^T=CCUG 70961^T), Bifidobacterium italicum sp. nov. (Rab10A=LMG 30187^T=CCUG 70963^T), Bifidobacterium margollesii sp. nov. (Uis1B=LMG 30296^T=CCUG 70959^T) and Bifidobacterium parmae sp. nov. (Uis4E=LMG 30295^T=CCUG 70964^T) are proposed as novel Bifidobacterium species.

Tracking the Taxonomy of the Genus Bifidobacterium Based on a Phylogenomic Approach

For decades, bacterial taxonomy has been based on in vitro molecular biology techniques and comparison of molecular marker sequences to measure the degree of genetic similarity and deduce phylogenetic relatedness of novel bacterial species to reference microbial taxa. Due to the advent of the genomic era, access to complete bacterial genome contents has become easier, thereby presenting the opportunity to precisely investigate the overall genetic diversity of microorganisms. Here, we describe a high-accuracy phylogenomic approach to assess the taxonomy of members of the genus Bifidobacterium and identify apparent misclassifications in current bifidobacterial taxonomy. The developed method was validated by the classification of seven novel taxa belonging to the genus Bifidobacterium by employing their overall genetic content. The results of this study demonstrate the potential of this whole-genome approach to become the gold standard for phylogenomics-based taxonomic classification of bacteria.IMPORTANCE Nowadays, next-generation sequencing has given access to genome sequences of the currently known bacterial taxa. The public databases constructed by means of these new technologies allowed comparison of genome sequences between microorganisms, providing information to perform genomic, phylogenomic, and evolutionary analyses. In order to avoid misclassifications in the taxonomy of novel bacterial isolates, new (bifido)bacterial taxons should be validated with a phylogenomic assessment like the approach presented here.

Glutamate synthases from conifers: gene structure and phylogenetic studies

BACKGROUND

Plants synthesize glutamate from ammonium by the combined activity of the enzymes glutamine synthetase (GS) and glutamate synthase (GOGAT) through the glutamate synthase cycle. In plants, there are two forms of glutamate synthases that differ in their electron donors, NADH-GOGAT (EC 1.4.1.14) and Fd-GOGAT (EC 1.4.7.1), which have differential roles either in primary ammonia assimilation or in the reassimilation of ammonium from different catabolic processes. Glutamate synthases are complex iron-sulfur flavoproteins containing functional domains involved in the control and coordination of their catalytic activities in annual plants. In conifers, partial cDNA sequences for GOGATs have been isolated and used for gene expression studies. However, knowledge of the gene structure and of phylogenetic relationships with other plant enzymes is quite scant.

RESULTS

Technological advances in conifer megagenomes sequencing have made it possible to obtain full-length cDNA sequences encoding Fd- and NADH-GOGAT from maritime pine, as well as BAC clones containing sequences for NADH-GOGAT and Fd-GOGAT genes. In the current study, we studied the genomic organization of pine GOGAT genes, the size of their exons/introns, copy numbers in the pine genome and relationships with other plant genes. Phylogenetic analysis was performed, and the degree of preservation and dissimilarity of key domains for the catalytic activities of these enzymes in different taxa were determined.

CONCLUSIONS

Fd- and NADH-GOGAT are encoded by single-copy genes in the maritime pine genome. The Fd-GOGAT gene is extremely large spanning more than 330 kb and the presence of very long introns highlights the important contribution of LTR retrotransposons to the gene size in conifers. In contrast, the structure of the NADH-GOGAT gene is similar to the orthologous genes in angiosperms. Our phylogenetic analysis indicates that these two genes had different origins during plant evolution. The results provide new insights into the structure and molecular evolution of these essential genes.

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains

BACKGROUND

Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation.

RESULTS

We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains.

CONCLUSIONS

In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.

Discovering viral genomes in human metagenomic data by predicting unknown protein families

Massive amounts of metagenomics data are currently being produced, and in all such projects a sizeable fraction of the resulting data shows no or little homology to known sequences. It is likely that this fraction contains novel viruses, but identification is challenging since they frequently lack homology to known viruses. To overcome this problem, we developed a strategy to detect ORFan protein families in shotgun metagenomics data, using similarity-based clustering and a set of filters to extract bona fide protein families. We applied this method to 17 virus-enriched libraries originating from human nasopharyngeal aspirates, serum, feces, and cerebrospinal fluid samples. This resulted in 32 predicted putative novel gene families. Some families showed detectable homology to sequences in metagenomics datasets and protein databases after reannotation. Notably, one predicted family matches an ORF from the highly variable Torque Teno virus (TTV). Furthermore, follow-up from a predicted ORFan resulted in the complete reconstruction of a novel circular genome. Its organisation suggests that it most likely corresponds to a novel bacteriophage in the microviridae family, hence it was named bacteriophage HFM.