Viral impacts on honey bee populations: A review

Honey bee is vital for pollination and ecological services, boosting crops productivity in terms of quality and quantity and production of colony products: wax, royal jelly, bee venom, honey, pollen and propolis. Honey bees are most important plant pollinators and almost one third of diet depends on bee’s pollination, worth billions of dollars. Hence the role that honey bees have in environment and their economic importance in food production, their health is of dominant significance. Honey bees can be infected by various pathogens like: viruses, bacteria, fungi, or infested by parasitic mites. At least more than 20 viruses have been identified to infect honey bees worldwide, generally from Dicistroviridae as well as Iflaviridae families, like ABPV (Acute Bee Paralysis Virus), BQCV (Black Queen Cell Virus), KBV (Kashmir Bee Virus), SBV (Sacbrood Virus), CBPV (Chronic bee paralysis virus), SBPV (Slow Bee Paralysis Virus) along with IAPV (Israeli acute paralysis virus), and DWV (Deformed Wing Virus) are prominent and cause infections harmful for honey bee colonies health. This issue about honey bee viruses demonstrates remarkably how diverse this field is, and considerable work has to be done to get a comprehensive interpretation of the bee virology.     

FragGeneScan: predicting genes in short and error-prone reads

The advances of next-generation sequencing technology have facilitated metagenomics research that attempts to determine directly the whole collection of genetic material within an environmental sample (i.e. the metagenome). Identification of genes directly from short reads has become an important yet challenging problem in annotating metagenomes, since the assembly of metagenomes is often not available. Gene predictors developed for whole genomes (e.g. Glimmer) and recently developed for metagenomic sequences (e.g. MetaGene) show a significant decrease in performance as the sequencing error rates increase, or as reads get shorter. We have developed a novel gene prediction method FragGeneScan, which combines sequencing error models and codon usages in a hidden Markov model to improve the prediction of protein-coding region in short reads. The performance of FragGeneScan was comparable to Glimmer and MetaGene for complete genomes. But for short reads, FragGeneScan consistently outperformed MetaGene (accuracy improved ∼62% for reads of 400 bases with 1% sequencing errors, and ∼18% for short reads of 100 bases that are error free). When applied to metagenomes, FragGeneScan recovered substantially more genes than MetaGene predicted (>90% of the genes identified by homology search), and many novel genes with no homologs in current protein sequence database.     

A curated data resource of 214K metagenomes for characterization of the global antimicrobial resistome

The growing threat of antimicrobial resistance (AMR) calls for new epidemiological surveillance methods, as well as a deeper understanding of how antimicrobial resistance genes (ARGs) have been transmitted around the world. The large pool of sequencing data available in public repositories provides an excellent resource for monitoring the temporal and spatial dissemination of AMR in different ecological settings. However, only a limited number of research groups globally have the computational resources to analyze such data. We retrieved 442 Tbp of sequencing reads from 214,095 metagenomic samples from the European Nucleotide Archive (ENA) and aligned them using a uniform approach against ARGs and 16S/18S rRNA genes. Here, we present the results of this extensive computational analysis and share the counts of reads aligned. Over 6.76∙10⁸ read fragments were assigned to ARGs and 3.21∙10⁹ to rRNA genes, where we observed distinct differences in both the abundance of ARGs and the link between microbiome and resistome compositions across various sampling types. This collection is another step towards establishing global surveillance of AMR and can serve as a resource for further research into the environmental spread and dynamic changes of ARGs.

The growing threat of antimicrobial resistance (AMR) calls for new epidemiological surveillance methods and a deeper understanding of how resistance genes are transmitted around the world. This study presents a large-scale remapping of sequencing reads of publicly available metagenomic datasets that can be used to monitor the global prevalence of AMR genes.     

A molecular phylogeny of the genus Apis suggests that the Giant Honey Bee of the Philippines,A. breviligula Maa,and the Plains Honey Bee of southern India,A. indica Fabricius,are valid species

Two new taxa have been added to the existing molecular phylogenies of the genus Apis. The new phylogeny supports the accepted phylogenetic relationships of {dwarf honey bees [giant honey bees (cavity‐nesting honey bees)]}. Based on Bayesian and maximum parsimony trees, our analysis supports recognition of Apis indica, the Plains Honey Bee of south India, as a separate species from A. cerana. Our analysis also supports recognition of the Giant Philippines Honey Bee, A. breviligula, as a separate species from A. dorsata.     

In-depth Proteome of the Hypopharyngeal Glands of Honeybee Workers Reveals Highly Activated Protein and Energy Metabolism in Priming the Secretion of Royal Jelly

1. Download : Download high-res image (120KB)
2. Download : Download full-size image

Highlights

• •In-depth proteome underpins the gland ontogeny and age-specific activity of the HGs.
• •The well-developed acini in the HGs of NBs promote the RJ secretary activities.
• •The enhanced protein and energy metabolism in the HGs boost the stronger RJ secretion of RJBs.

     

A fault-tolerant method for HLA typing with PacBio data

Background

Human leukocyte antigen (HLA) genes are critical genes involved in important biomedical aspects, including organ transplantation, autoimmune diseases and infectious diseases. The gene family contains the most polymorphic genes in humans and the difference between two alleles is only a single base pair substitution in many cases. The next generation sequencing (NGS) technologies could be used for high throughput HLA typing but in silico methods are still needed to correctly assign the alleles of a sample. Computer scientists have developed such methods for various NGS platforms, such as Illumina, Roche 454 and Ion Torrent, based on the characteristics of the reads they generate. However, the method for PacBio reads was less addressed, probably owing to its high error rates. The PacBio system has the longest read length among available NGS platforms, and therefore is the only platform capable of having exon 2 and exon 3 of HLA genes on the same read to unequivocally solve the ambiguity problem caused by the “phasing” issue.

Results

We proposed a new method BayesTyping1 to assign HLA alleles for PacBio circular consensus sequencing reads using Bayes’ theorem. The method was applied to simulated data of the three loci HLA-A, HLA-B and HLA-DRB1. The experimental results showed its capability to tolerate the disturbance of sequencing errors and external noise reads.

Conclusions

The BayesTyping1 method could overcome the problems of HLA typing using PacBio reads, which mostly arise from sequencing errors of PacBio reads and the divergence of HLA genes, to some extent.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-296) contains supplementary material, which is available to authorized users.     

SAP-A Sequence Mapping and Analyzing Program for Long Sequence Reads Alignment and Accurate Variants Discovery

The third-generation of sequencing technologies produces sequence reads of 1000 bp or more that may contain high polymorphism information. However, most currently available sequence analysis tools are developed specifically for analyzing short sequence reads. While the traditional Smith-Waterman (SW) algorithm can be used to map long sequence reads, its naive implementation is computationally infeasible. We have developed a new Sequence mapping and Analyzing Program (SAP) that implements a modified version of SW to speed up the alignment process. In benchmarks with simulated and real exon sequencing data and a real E. coli genome sequence data generated by the third-generation sequencing technologies, SAP outperforms currently available tools for mapping short and long sequence reads in both speed and proportion of captured reads. In addition, it achieves high accuracy in detecting SNPs and InDels in the simulated data. SAP is available at https://github.com/davidsun/SAP.     

Enterocin AS-48RJ: a variant of enterocin AS-48 chromosomally encoded by Enterococcus faecium RJ16 isolated from food

The bacteriocinogenic strain RJ16 isolated from goat cheese has been identified as Enterococcusfaecium by species-specific PCR, DNA-rRNA hybridization and rDNA sequencing. Purified bacteriocin from strain RJ16 is a carboxypeptidase A-resistant peptide with a molecular mass (7125 Da) very close to the cyclic peptide enterocin AS-48. Bacteriocin from strain RJ16 and AS-48 show identical antibacterial spectra, although the former is slightly less active on strains of Listeria monocytogenes and Bacillus cereus. Producer strains show cross-immunity. PCR amplification of total DNA from strain RJ16 with primers for the AS-48 structural gene and sequencing of the amplified fragment revealed an almost identical sequence (99.5%), except for a single mutation that predicts the change of Glu residue at position 20 of AS-48 to Val. Therefore, bacteriocin produced by E. faecium RJ16 should be considered a variant of AS-48, which we call AS-48RJ. PCR amplification revealed that strain RJ16 contains the complete as-48. gene cluster. Hybridization with probes for as-48 gene cluster revealed a chromosomal location of as-48 genes in strain RJ16, being the first example of a chromosomal location of this bacteriocin trait. Strain RJ16 produced enzymes of interest in food processing (esterase, esterase lipase and phytase activities), and did not decarboxylate amino acids precursors for biogenic amines. Strain RJ16 did not exhibit haemolytic or gelatinase activities, and PCR amplification revealed the lack of genes encoding for known virulence determinants (aggregation substance, collagen adhesin, enterococcal surface protein, endocarditis antigens, as well as haemolysin and gelatinase production). Strain RJ16 was resistant to ciprofloxacin (MIC > 2 mgl(-1)) and levofloxacin (MIC > 4 mgl(-1)) and showed intermediate resistance to nitrofurantoin and erythromycin, but was sensitive to ampicillin, penicillin, streptomycin, gentamicin, rifampicin, chloramphenicol, tetracycline, quinupristin/dalfopristin, vancomycin and teicoplanin. Altogether, results from this study suggest that this broad-spectrum bacteriocin-producing strain may have a potential use in food preservation.     

Decoding the oak genome: public release of sequence data,assembly, annotation and publication strategies

The 1.5 Gbp/2C genome of pedunculate oak (Quercus robur) has been sequenced. A strategy was established for dealing with the challenges imposed by the sequencing of such a large, complex and highly heterozygous genome by a whole‐genome shotgun (WGS) approach, without the use of costly and time‐consuming methods, such as fosmid or BAC clone‐based hierarchical sequencing methods. The sequencing strategy combined short and long reads. Over 49 million reads provided by Roche 454 GS‐FLX technology were assembled into contigs and combined with shorter Illumina sequence reads from paired‐end and mate‐pair libraries of different insert sizes, to build scaffolds. Errors were corrected and gaps filled with Illumina paired‐end reads and contaminants detected, resulting in a total of 17 910 scaffolds (>2 kb) corresponding to 1.34 Gb. Fifty per cent of the assembly was accounted for by 1468 scaffolds (N50 of 260 kb). Initial comparison with the phylogenetically related Prunus persica gene model indicated that genes for 84.6% of the proteins present in peach (mean protein coverage of 90.5%) were present in our assembly. The second and third steps in this project are genome annotation and the assignment of scaffolds to the oak genetic linkage map. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement, the oak genome data have been released into public sequence repositories in advance of publication. In this presubmission paper, the oak genome consortium describes its principal lines of work and future directions for analyses of the nature, function and evolution of the oak genome.     

Proteome comparison of hypopharyngeal gland development between Italian and royal jelly producing worker honeybees (Apis mellifera L.)

The hypopharyngeal gland (HG) of the honeybee (Apis mellifera L.) produces royal jelly (RJ) that is essential to feed and raise broods and queens. A strain of bees (high royal jelly producing bee, RJb) has been selected for its high RJ production, but the mechanisms of its higher yield are not understood. In this study, we compared HG acini size, RJ production, and protein differential expressions between the RJb and nonselected honeybee (Italian bee, ITb) using proteomics in combination with an electron microscopy, Western blot, and quantitative real-time PCR (qRT-PCR). Generally, the HG of both bees showed age-dependent changes in acini sizes and protein expression as worker behaviors changed from brood nursing to nectar ripening, foraging, and storage activities. The electron microscopic analysis revealed that the HG acini diameter of the RJb strain was large and produced 5 times more RJ than the ITb, demonstrating a positive correlation between the yield and HG acini size. In addition, the proteomic analysis showed that RJb significantly upregulated a large group of proteins involved in carbohydrate metabolism and energy production, those involved in protein biosynthesis, development, amino acid metabolism, nucleotide and fatty acid, transporter, protein folding, cytoskeleton, and antioxidation, which coincides with the fact that the HGs of the RJb strain produce more RJ than the ITb strain that is owing to selection pressure. We also observed age-dependent major royal jelly proteins (MRJPs) changing both in form and expressional intensity concurrent with task-switching. In addition to MRJPs, the RJb overexpressed proteins such as enolase and transitional endoplasmic reticulum ATPase, protein biosynthesis, and development proteins compared to the ITb strain to support its large HG growth and RJ secretion. Because of selection pressure, RJb pursued a different strategy of increased RJ production by involving additional proteins compared to its original counterpart ITb. To our knowledge, this morphological and proteomic comparison study on the HG of the two strains of worker honeybees associated with their age-dependent division of labor is the first of its kind. The study provided not only the quantity and quality differences in the HG from the RJb and the ITb, but also addressed the cellular and behavioral biology development question of how the RJb strain can produce RJ more efficiently than its wild type strain (ITb).     

A Next Generation Semiconductor Based Sequencing Approach for the Identification of Meat Species in DNA Mixtures

The identification of the species of origin of meat and meat products is an important issue to prevent and detect frauds that might have economic, ethical and health implications. In this paper we evaluated the potential of the next generation semiconductor based sequencing technology (Ion Torrent Personal Genome Machine) for the identification of DNA from meat species (pig, horse, cattle, sheep, rabbit, chicken, turkey, pheasant, duck, goose and pigeon) as well as from human and rat in DNA mixtures through the sequencing of PCR products obtained from different couples of universal primers that amplify 12S and 16S rRNA mitochondrial DNA genes. Six libraries were produced including PCR products obtained separately from 13 species or from DNA mixtures containing DNA from all species or only avian or only mammalian species at equimolar concentration or at 1:10 or 1:50 ratios for pig and horse DNA. Sequencing obtained a total of 33,294,511 called nucleotides of which 29,109,688 with Q20 (87.43%) in a total of 215,944 reads. Different alignment algorithms were used to assign the species based on sequence data. Error rate calculated after confirmation of the obtained sequences by Sanger sequencing ranged from 0.0003 to 0.02 for the different species. Correlation about the number of reads per species between different libraries was high for mammalian species (0.97) and lower for avian species (0.70). PCR competition limited the efficiency of amplification and sequencing for avian species for some primer pairs. Detection of low level of pig and horse DNA was possible with reads obtained from different primer pairs. The sequencing of the products obtained from different universal PCR primers could be a useful strategy to overcome potential problems of amplification. Based on these results, the Ion Torrent technology can be applied for the identification of meat species in DNA mixtures.