HAMAP: a database of completely sequenced microbial proteome sets and manually curated microbial protein families in UniProtKB/Swiss-Prot

CANDy: Automated analysis of domain architectures in carbohydrate-active enzymes

Carbohydrate-active enzymes (CAZymes) can be found in all domains of life and play a crucial role in metabolic and physiological processes. CAZymes often possess a modular structure, comprising not only catalytic domains but also associated domains such as carbohydrate-binding modules (CBMs) and linker domains. By exploring the modular diversity of CAZy families, catalysts with novel properties can be discovered and further insight in their biological functions and evolutionary relationships can be obtained. Here we present the carbohydrate-active enzyme domain analysis tool (CANDy), an assembly of several novel scripts, tools and databases that allows users to analyze the domain architecture of all protein sequences in a given CAZy family. CANDy's usability is shown on glycoside hydrolase family 48, a small yet underexplored family containing multi-domain enzymes. Our analysis reveals the existence of 35 distinct domain assemblies, including eight known architectures, with the remaining assemblies awaiting characterization. Moreover, we substantiate the occurrence of horizontal gene transfer from prokaryotes to insect orthologs and provide evidence for the subsequent removal of auxiliary domains, likely through a gene fission event. CANDy is available at https://github.com/PyEED/CANDy.

Genetic insights: mapping sex-specific loci in Siamese cobra (Naja kaouthia) sheds light on the putative sex determining region

The location of female-specific/linked loci identified in Siamese cobra (Naja kaouthia) previously has been determined through in silico chromosome mapping of the Indian cobra genome (N. naja) as a reference genome. In the present study, we used in silico chromosome mapping to identify sex-specific and linked loci in Siamese cobra. Many sex-specific and sex-linked loci were successfully mapped on the Z sex chromosome, with 227 of the 475 specific loci frequently mapped in a region covering 57 Mb and positioned at 38,992,675-95,561,177 bp of the Indian cobra genome (N. naja). This suggested the existence of a putative sex-determining region (SDR), with one specific locus (PA100000600) homologous to the TOPBP1 gene. The involvement of TOPBP1 gene may lead to abnormal synaptonemal complexes and meiotic chromosomal defects, resulting in male infertility. These findings offer valuable insights into the genetic basis and functional aspects of sex-specific traits in the Siamese cobra, which will contribute to our understanding of snake genetics and evolutionary biology.

Genome-wide SNP analysis provides insights into the XX/XY sex-determination system in silver barb (Barbonymus gonionotus)

Silver barb (Barbonymus gonionotus) is among the most economically important freshwater fish species in Thailand. It ranks fourth in economic value and third in production weight for fisheries and culture in Thailand. An XX/XY sex-determination system based on gynogenesis was previously reported for this fish. In this study, the molecular basis underlying the sex-determination system was further investigated. Genome-wide single-nucleotide polymorphism data were generated for 32 captive-bred silver barb individuals, previously scored by phenotypic sex, to identify sex-linked regions associated with sex determination. Sixty-three male-linked loci, indicating putative XY chromosomes, were identified. Male-specific loci were not observed, which indicates that the putative Y chromosome is young and the sex determination region is cryptic. A homology search revealed that most male-linked loci were homologous to the Mariner/Tc1 and Gypsy transposable elements and are probably the remnants of an initial accumulation of repeats on the Y chromosome from the early stages of sex chromosome differentiation. This research provides convincing insights into the mechanism of sex determination and reveals the potential sex determination regions in silver barb. The study provides the basic data necessary for increasing the commercial value of silver barbs through genetic improvements.

Identification and cultivation of anaerobic bacterial scavengers of dead cells

The cycle of life and death and Earth's carbon cycle(s) are intimately linked, yet how bacterial cells, one of the largest pools of biomass on Earth, are recycled back into the carbon cycle remains enigmatic. In particular, no bacteria capable of scavenging dead cells in oxygen-depleted environments have been reported thus far. In this study, we discover the first anaerobes that scavenge dead cells and the two isolated strains use distinct strategies. Based on live-cell imaging, transmission electron microscopy, and hydrolytic enzyme assays, one strain (designated CYCD) relied on cell-to-cell contact and cell invagination for degrading dead food bacteria where as the other strain (MGCD) degraded dead food bacteria via excretion of lytic extracellular enzymes. Both strains could degrade dead cells of differing taxonomy (bacteria and archaea) and differing extents of cell damage, including those without artificially inflicted physical damage. In addition, both depended on symbiotic metabolic interactions for maximizing cell degradation, representing the first cultured syntrophic Bacteroidota. We collectively revealed multiple symbiotic bacterial decomposition routes of dead prokaryotic cells, providing novel insight into the last step of the carbon cycle.

Preserving Pure Siamese Crocodile Populations: A Comprehensive Approach Using Multi-Genetic Tools

Hybrids between the critically endangered Siamese crocodile (Crocodylus siamensis) and least-concern saltwater crocodile (C. porosus) in captive populations represent a serious challenge for conservation and reintroduction programs due to the impact of anthropogenic activities. A previous study used microsatellite and mitochondrial DNA data to establish the criteria for identifying species and their hybrids; however, the results may have been influenced by biased allelic frequencies and genetic drift within the examined population. To overcome these limitations and identify the true signals of selection, alternative DNA markers and a diverse set of populations should be employed. Therefore, this study used DArT sequencing to identify genome-wide single nucleotide polymorphisms (SNPs) in both species and confirm the genetic scenario of the parental species and their hybrids. A population of saltwater crocodiles from Australia was used to compare the distribution of species-diagnostic SNPs. Different analytical approaches were compared to diagnose the level of hybridization when an admixture was present, wherein three individuals had potential backcrossing. Approximately 17.00-26.00% of loci were conserved between the Siamese and saltwater crocodile genomes. Species-diagnostic SNP loci for Siamese and saltwater crocodiles were identified as 8051 loci and 1288 loci, respectively. To validate the species-diagnostic SNP loci, a PCR-based approach was used by selecting 20 SNP loci for PCR primer design, among which 3 loci were successfully able to differentiate the actual species and different hybridization levels. Mitochondrial and nuclear genetic information, including microsatellite genotyping and species-diagnostic DNA markers, were combined as a novel method that can compensate for the limitations of each method. This method enables conservation prioritization before release into the wild, thereby ensuring sustainable genetic integrity for long-term species survival through reintroduction and management programs.

Something Fishy about Siamese Fighting Fish (Betta splendens) Sex: Polygenic Sex Determination or a Newly Emerged Sex-Determining Region?

Fishes provide a unique and intriguing model system for studying the genomic origin and evolutionary mechanisms underlying sex determination and high sex-chromosome turnover. In this study, the mode of sex determination was investigated in Siamese fighting fish, a species of commercial importance. Genome-wide SNP analyses were performed on 75 individuals (40 males and 35 females) across commercial populations to determine candidate sex-specific/sex-linked loci. In total, 73 male-specific loci were identified and mapped to a 5.6 kb region on chromosome 9, suggesting a putative male-determining region (pMDR) containing localized dmrt1 and znrf3 functional sex developmental genes. Repeat annotations of the pMDR revealed an abundance of transposable elements, particularly Ty3/Gypsy and novel repeats. Remarkably, two out of the 73 male-specific loci were located on chromosomes 7 and 19, implying the existence of polygenic sex determination. Besides male-specific loci, five female-specific loci on chromosome 9 were also observed in certain populations, indicating the possibility of a female-determining region and the polygenic nature of sex determination. An alternative explanation is that male-specific loci derived from other chromosomes or female-specific loci in Siamese fighting fish recently emerged as new sex-determining loci during domestication and repeated hybridization.

Chromosome-scale genome assembly of an important medicinal plant honeysuckle

Lonicera japonica (honeysuckle) is one of the most important medicinal plants and widely utilized in traditional Chinese medicine. At present, there are many varieties of honeysuckle used in cultivation, among which Sijihua variety are widely cultivated due to its wide adaptability, stress resistance, early flowering and high yield. In this study, we assembled the genome of Sijihua, which was approximately 886.04 Mb in size with a scaffold N50 of 79.5 Mb. 93.28% of the total assembled sequences were anchored to 9 pseudo-chromosomes by using PacBio long reads and Hi-C sequencing data. We predicted 39,320 protein-coding genes and 92.87% of them could be annotated in NR, GO, KOG, KEGG and other databases. In addition, we identified 644 tRNAs, 2,156 rRNAs, 109 miRNAs and 5,502 pseudogenes from the genome. The chromosome-scale genome of Sijihua will be a significant resource for understanding the genetic basis of high stress-resistance, which will facilitate further study of the genetic diversity and accelerate the genetic improvement and breeding of L. japonica.

Measurement(s)	Lonicera japonica • RNA sequencing • genome assembly • sequence annotation
Technology Type(s)	SMRT Sequencing • RNA sequencing • Hi-C • biomolecular annotation design
Factor Type(s)	Genotype
Sample Characteristic - Organism	Lonicera japonica
Sample Characteristic - Environment	occurrence
Sample Characteristic - Location	Shandong Province

Genome-Wide Survey and Development of the First Microsatellite Markers Database (AnCorDB) in Anemone coronaria L

Anemone coronaria L. (2n = 2x = 16) is a perennial, allogamous, highly heterozygous plant marketed as a cut flower or in gardens. Due to its large genome size, limited efforts have been made in order to develop species-specific molecular markers. We obtained the first draft genome of the species by Illumina sequencing an androgenetic haploid plant of the commercial line “MISTRAL^® Magenta”. The genome assembly was obtained by applying the MEGAHIT pipeline and consisted of 2 × 10⁶ scaffolds. The SciRoKo SSR (Simple Sequence Repeats)-search module identified 401.822 perfect and 188.987 imperfect microsatellites motifs. Following, we developed a user-friendly “Anemone coronaria Microsatellite DataBase” (AnCorDB), which incorporates the Primer3 script, making it possible to design couples of primers for downstream application of the identified SSR markers. Eight genotypes belonging to eight cultivars were used to validate 62 SSRs and a subset of markers was applied for fingerprinting each cultivar, as well as to assess their intra-cultivar variability. The newly developed microsatellite markers will find application in Breeding Rights disputes, developing genetic maps, marker assisted breeding (MAS) strategies, as well as phylogenetic studies.

Integration of Aspergillus niger transcriptomic profile with metabolic model identifies potential targets to optimise citric acid production from lignocellulosic hydrolysate

BACKGROUND

Citric acid is typically produced industrially by Aspergillus niger-mediated fermentation of a sucrose-based feedstock, such as molasses. The fungus Aspergillus niger has the potential to utilise lignocellulosic biomass, such as bagasse, for industrial-scale citric acid production, but realising this potential requires strain optimisation. Systems biology can accelerate strain engineering by systematic target identification, facilitated by methods for the integration of omics data into a high-quality metabolic model. In this work, we perform transcriptomic analysis to determine the temporal expression changes during fermentation of bagasse hydrolysate and develop an evolutionary algorithm to integrate the transcriptomic data with the available metabolic model to identify potential targets for strain engineering.

RESULTS

The novel integrated procedure matures our understanding of suboptimal citric acid production and reveals potential targets for strain engineering, including targets consistent with the literature such as the up-regulation of citrate export and pyruvate carboxylase as well as novel targets such as the down-regulation of inorganic diphosphatase.

CONCLUSIONS

In this study, we demonstrate the production of citric acid from lignocellulosic hydrolysate and show how transcriptomic data across multiple timepoints can be coupled with evolutionary and metabolic modelling to identify potential targets for further engineering to maximise productivity from a chosen feedstock. The in silico strategies employed in this study can be applied to other biotechnological goals, assisting efforts to harness the potential of microorganisms for bio-based production of valuable chemicals.

Chromosome-level genome assembly of the shuttles hoppfish, Periophthalmus modestus

Background

The shuttles hoppfish (mudskipper), Periophthalmus modestus, is one of the mudskippers, which are the largest group of amphibious teleost fishes, which are uniquely adapted to live on mudflats. Because mudskippers can survive on land for extended periods by breathing through their skin and through the lining of the mouth and throat, they were evaluated as a model for the evolutionary sea-land transition of Devonian protoamphibians, ancestors of all present tetrapods.

Results

A total of 39.6, 80.2, 52.9, and 33.3 Gb of Illumina, Pacific Biosciences, 10X linked, and Hi-C data, respectively, was assembled into 1,419 scaffolds with an N50 length of 33 Mb and BUSCO score of 96.6%. The assembly covered 117% of the estimated genome size (729 Mb) and included 23 pseudo-chromosomes anchored by a Hi-C contact map, which corresponded to the top 23 longest scaffolds above 20 Mb and close to the estimated one. Of the genome, 43.8% were various repetitive elements such as DNAs, tandem repeats, long interspersed nuclear elements, and simple repeats. Ab initio and homology-based gene prediction identified 30,505 genes, of which 94% had homology to the 14 Actinopterygii transcriptomes and 89% and 85% to Pfam familes and InterPro domains, respectively. Comparative genomics with 15 Actinopterygii species identified 59,448 gene families of which 12% were only in P. modestus.

Conclusions

We present the high quality of the first genome assembly and gene annotation of the shuttles hoppfish. It will provide a valuable resource for further studies on sea-land transition, bimodal respiration, nitrogen excretion, osmoregulation, thermoregulation, vision, and mechanoreception.

Chromosome-level genome assembly of Bactrocera dorsalis reveals its adaptation and invasion mechanisms

Bactrocera dorsalis is an invasive polyphagous pest causing considerable ecological and economic damage worldwide. We report a high-quality chromosome-level genome assembly and combine various transcriptome data to explore the molecular mechanisms of its rapid adaptation to new environments. The expansions of the DDE transposase superfamily and key gene families related to environmental adaptation and enrichment of the expanded and unique gene families in metabolism and defence response pathways explain its environmental adaptability. The relatively high but not significantly different expression of heat-shock proteins, regardless of the environmental conditions, suggests an intrinsic mechanism underlying its adaptation to high temperatures. The mitogen-activated protein kinase pathway plays a key role in adaptation to new environments. The prevalence of duplicated genes in its genome explains the diversity in the B. dorsalis complex. These findings provide insights into the genetic basis of the invasiveness and diversity of B. dorsalis, explaining its rapid adaptation and expansion.

Jiang et al. sequence the genome of Bactrocera dorsalis, a destructive and invasive agricultural pest. Insights from this chromosome-level assembly shed light on molecular adaptations that allow for the global invasion and expansion of this pest.

Whole genome sequencing and bulked segregant analysis suggest a new mechanism of amitraz resistance in the citrus red mite, Panonychus citri (Acari: Tetranychidae)

BACKGROUND

Amitraz is a broad-spectrum insecticide/acaricide for the control of aphids, psyllids, ticks and mites. Current evidence suggests that ticks and phytophagous mites have developed strong resistance to amitraz. Previous studies have shown that multiple mechanisms are associated with amitraz resistance in ticks, but very few reports have involved Panonychus citri. We therefore used whole genome sequencing and bulked segregant analysis (BSA) to identify the mechanism underlying P. citri's resistance to amitraz.

RESULTS

High-quality assembly of the whole P. citri genome was completed, resulting in a genome of approximately 83.97 Mb and a contig N50 of approximately 1.81 Mb. Gene structure predictions revealed 11 577 genes, of which 10 940 genes were annotated. Trait-associated regions in the genome were mapped with bulked segregant analysis and 38 candidate SNPs were obtained, of which T752C had the strongest correlation with the resistant trait, located at the 5' untranslated region (UTR) of the β-2R adrenergic-like octopamine receptor gene. The mutation resulted in the formation of a short hairpin loop structure in mRNA and gene expression was down-regulated by more than 50% in the amitraz-resistant strain. Validation of the T752C mutation in field populations of P. citri found that the correlation between the resistance ratio and the base mutation was 94.40%.

CONCLUSION

Our results suggest that this 5' UTR mutation of the β-2R octopamine receptor gene, confers amitraz resistance in P. citri. This discovery provides a new explanation for the mechanism of pest resistance: base mutations in the 5' untranslated region of target gene may regulate the susceptibility of pests to pesticides.

A high-quality reference genome sequence of Salvia miltiorrhiza provides insights into tanshinone synthesis in its red rhizomes

Salvia miltiorrhiza Bunge, also known as red sage or Danshen, is an important traditional Chinese medicine (TCM) that has been used for thousands of years to treat cardiovascular and other diseases. It is also considered an important model TCM plant. Here, a high-quality reference genome of S. miltiorrhiza was generated by combining PacBio long-read sequencing and chromatin interaction mapping (Hi-C) technologies, resulting in the chromosome-scale assembly of a 594.75-Mb genome sequence with a contig N50 of 2.70 Mb. This assembly shows the highest level of continuity for a Danshen genome generated thus far. The S. miltiorrhiza genome contained 32,483 protein-coding genes, with a repetitive DNA content of approximately 64.84%. The high percentage of young LTRs suggests that multiple TE transposition bursts occurred recently in S. miltiorrhiza. Genes unique to secondary metabolism pathways were expanded in the S. miltiorrhiza genome. A new CYP450 gene cluster was identified in the phloem of red roots where active components were synthesized. This reference genome sequence will facilitate future studies aimed at the elucidation of the secondary metabolism synthesis pathway and the genetic improvement of S. miltiorrhiza.

"Mind the Gap": Hi-C Technology Boosts Contiguity of the Globe Artichoke Genome in Low-Recombination Regions

Globe artichoke (Cynara cardunculus var. scolymus; 2n2x=34) is cropped largely in the Mediterranean region, being Italy the leading world producer; however, over time, its cultivation has spread to the Americas and China. In 2016, we released the first (v1.0) globe artichoke genome sequence (http://www.artichokegenome.unito.it/). Its assembly was generated using ∼133-fold Illumina sequencing data, covering 725 of the 1,084 Mb genome, of which 526 Mb (73%) were anchored to 17 chromosomal pseudomolecules. Based on v1.0 sequencing data, we generated a new genome assembly (v2.0), obtained from a Hi-C (Dovetail) genomic library, and which improves the scaffold N₅₀ from 126 kb to 44.8 Mb (∼356-fold increase) and N₉₀ from 29 kb to 17.8 Mb (∼685-fold increase). While the L₉₀ of the v1.0 sequence included 6,123 scaffolds, the new v2.0 just 15 super-scaffolds, a number close to the haploid chromosome number of the species. The newly generated super-scaffolds were assigned to pseudomolecules using reciprocal blast procedures. The cumulative size of unplaced scaffolds in v2.0 was reduced of 165 Mb, increasing to 94% the anchored genome sequence. The marked improvement is mainly attributable to the ability of the proximity ligation-based approach to deal with both heterochromatic (e.g.: peri-centromeric) and euchromatic regions during the assembly procedure, which allowed to physically locate low recombination regions. The new high-quality reference genome enhances the taxonomic breadth of the data available for comparative plant genomics and led to a new accurate gene prediction (28,632 genes), thus promoting the map-based cloning of economically important genes.

The Isolation and Genome Sequencing of Five Novel Bacteriophages From the Rumen Active Against Butyrivibrio fibrisolvens

Although the prokaryotic communities of the rumen microbiome are being uncovered through genome sequencing, little is known about the resident viral populations. Whilst temperate phages can be predicted as integrated prophages when analyzing bacterial and archaeal genomes, the genetics underpinning lytic phages remain poorly characterized. To the five genomes of bacteriophages isolated from rumen-associated samples sequenced and analyzed previously, this study adds a further five novel genomes and predictions gleaned from them to further the understanding of the rumen phage population. Lytic bacteriophages isolated from fresh ovine and bovine fecal and rumen fluid samples were active against the predominant fibrolytic ruminal bacterium Butyrivibrio fibrisolvens. The double stranded DNA genomes were sequenced and reconstructed into single circular complete contigs. Based on sequence similarity and genome distances, the five phages represent four species from three separate genera, consisting of: (1) Butyrivibrio phages Arian and Bo-Finn; (2) Butyrivibrio phages Idris and Arawn; and (3) Butyrivibrio phage Ceridwen. They were predicted to all belong to the Siphoviridae family, based on evidence in the genomes such as size, the presence of the tail morphogenesis module, genes that share similarity to those in other siphovirus isolates and phylogenetic analysis using phage proteomes. Yet, phylogenomic analysis and sequence similarity of the entire phage genomes revealed that these five phages are unique and novel. These phages have only been observed undergoing the lytic lifecycle, but there is evidence in the genomes of phages Arawn and Idris for the potential to be temperate. However, there is no evidence in the genome of the bacterial host Butyrivibrio fibrisolvens of prophage genes or genes that share similarity with the phage genomes.

Survival of the Halophilic Archaeon Halovarius luteus after Desiccation, Simulated Martian UV Radiation and Vacuum in Comparison to Bacillus atrophaeus

Extraterrestrial environments influence the biochemistry of organisms through a variety of factors, including high levels of radiation and vacuum, temperature extremes and a lack of water and nutrients. A wide variety of terrestrial microorganisms, including those counted amongst the most ancient inhabitants of Earth, can cope with high levels of salinity, extreme temperatures, desiccation and high levels of radiation. Key among these are the haloarchaea, considered particularly relevant for astrobiological studies due to their ability to thrive in hypersaline environments. In this study, a novel haloarchaea isolated from Urmia Salt Lake, Iran, Halovarius luteus strain DA50^T, was exposed to varying levels of simulated extraterrestrial conditions and compared to that of the bacteria Bacillus atrophaeus. Bacillus atrophaeus was selected for comparison due to its well-described resistance to extreme conditions and its ability to produce strong spore structures. Thin films were produced to investigate viability without the protective influence of cell multi-layers. Late exponential phase cultures of Hvr. luteus and B. atrophaeus were placed in brine and phosphate buffered saline media, respectively. The solutions were allowed to evaporate and cells were encapsulated and exposed to radiation, desiccation and vacuum conditions, and their post-exposure viability was studied by the Most Probable Number method. The protein profile using High Performance Liquid Chromatography and Matrix Assisted Laser Desorption/Ionization bench top reflector time-of-flight are explored after vacuum and UV-radiation exposure. Results showed that the change in viability of the spore-forming bacteria B. atrophaeus was only minor whereas Hvr. luteus demonstrated a range of viability under different conditions. At the peak radiation flux of 10⁵ J/m² under nitrogen flow and after two weeks of desiccation, Hvr. luteus demonstrated the greatest decrease in viability. This study further expands our understanding of the boundary conditions of astrobiologically relevant organisms in the harsh space environment.

Whole genome resequencing of four Italian sweet pepper landraces provides insights on sequence variation in genes of agronomic value

Sweet pepper (Capsicum annuum L.) is a high value crop and one of the most widely grown vegetables belonging to the Solanaceae family. In addition to commercial varieties and F₁ hybrids, a multitude of landraces are grown, whose genetic combination is the result of hundreds of years of random, environmental, and farmer selection. High genetic diversity exists in the landrace gene pool which however has scarcely been studied, thus bounding their cultivation. We re-sequenced four pepper inbred lines, within as many Italian landraces, which representative of as many fruit types: big sized blocky with sunken apex ('Quadrato') and protruding apex or heart shaped ('Cuneo'), elongated ('Corno') and smaller sized sub-spherical ('Tumaticot'). Each genomic sequence was obtained through Illumina platform at coverage ranging from 39 to 44×, and reconstructed at a chromosome scale. About 35.5k genes were predicted in each inbred line, of which 22,017 were shared among them and the reference genome (accession 'CM334'). Distinctive variations in miRNAs, resistance gene analogues (RGAs) and susceptibility genes (S-genes) were detected. A detailed survey of the SNP/Indels occurring in genes affecting fruit size, shape and quality identified the highest frequencies of variation in regulatory regions. Many structural variations were identified as presence/absence variations (PAVs), notably in resistance gene analogues (RGAs) and in the capsanthin/capsorubin synthase (CCS) gene. The large allelic diversity observed in the four inbred lines suggests their potential use as a pre-breeding resource and represents a one-stop resource for C. annuum genomics and a key tool for dissecting the path from sequence variation to phenotype.

Large-scale network analysis captures biological features of bacterial plasmids

Many bacteria can exchange genetic material through horizontal gene transfer (HGT) mediated by plasmids and plasmid-borne transposable elements. Here, we study the population structure and dynamics of over 10,000 bacterial plasmids, by quantifying their genetic similarities and reconstructing a network based on their shared k-mer content. We use a community detection algorithm to assign plasmids into cliques, which correlate with plasmid gene content, bacterial host range, GC content, and existing classifications based on replicon and mobility (MOB) types. Further analysis of plasmid population structure allows us to uncover candidates for yet undescribed replicon genes, and to identify transposable elements as the main drivers of HGT at broad phylogenetic scales. Our work illustrates the potential of network-based analyses of the bacterial 'mobilome' and opens up the prospect of a natural, exhaustive classification framework for bacterial plasmids.

A High-Quality Genome Sequence of Model Legume Lotus japonicus (MG-20) Provides Insights into the Evolution of Root Nodule Symbiosis

Lotus japonicus is an important model legume for studying symbiotic nitrogen fixation as well as plant development. A genomic sequence of L. japonicus (MG20) has been available for more than ten years. However, the low quality of the genome limits its application in functional genomic studies. Therefore, it is necessary to assemble high-quality chromosome sequences of L. japonicus using new sequencing technology to facilitate the study of functional genomics. In this report, we used the third-generation sequencing combined with the Illumina HiSeq platform to sequence the genome of L. japonicus (MG20). We obtained 544 Mb of genomic sequence using third-generation assembly. Based on sequence analysis, 357 Mb of repeats, 28,251 genes, 626 tRNAs, 1409 rRNAs, and 1233 pseudogenes were predicted in the genome. A total of 27,991 genes were annotated into databases. Compared to the previously published data, the new genome database contains complete L. japonicus sequences in the proper order and orientation with a contig N50 2.81Mb and an excellent genome coverage, which provides more accurate genome information and more precise assembly for functional genomic study.

First draft genome for the sand-hopper Trinorchestia longiramus

Crustacean amphipods are important trophic links between primary producers and higher consumers. Although most amphipods occur in or around aquatic environments, the family Talitridae is the only family found in terrestrial and semi-terrestrial habitats. The sand-hopper Trinorchestia longiramus is a talitrid species often found in the sandy beaches of South Korea. In this study, we present the first draft genome assembly and annotation of this species. We generated ~380.3 Gb of sequencing data assembled in a 0.89 Gb draft genome. Annotation analysis estimated 26,080 protein-coding genes, with 89.9% genome completeness. Comparison with other amphipods showed that T. longiramus has 327 unique orthologous gene clusters, many of which are expanded gene families responsible for cellular transport of toxic substances, homeostatic processes, and ionic and osmotic stress tolerance. This first talitrid genome will be useful for further understanding the mechanisms of adaptation in terrestrial environments, the effects of heavy metal toxicity, as well as for studies of comparative genomic variation across amphipods.

Draft Genome of the Asian Buffalo Leech Hirudinaria manillensis

The Asian Buffalo leech, Hirudinaria manillensis, is an aquatic sanguivorous species distributed widely in Southeast Asia. H. manillensis has long been used clinically for bloodletting and other medical purposes. Recent studies have focused on artificial culturing, strain optimization, and the identification and development new drugs based on the anticoagulant effects of H. manillensis bites; however, data regarding its genome remain unclear. This study aimed to determine the genome sequence of an adult Asian Buffalo leech. We generated a draft assembly of 151.8 Mb and a N50 scaffold of 2.28 Mb. Predictions indicated that the assembled genome contained 21,005 protein-coding genes. Up to 17,865 genes were annotated in multiple databases including Gene Ontology. Sixteen anticoagulant proteins with a Hirudin or Antistasin domain were identified. This study is the first to report the whole-genome sequence of the Asian Buffalo leech, an important sanguivorous leech of clinical significance. The quality of the assembly is comparable to those of other annelids. These data will help further the current understanding of the biological mechanisms and genetic characteristics of leeches and serve as a valuable resource for future studies.

Isolation of an archaeon at the prokaryote-eukaryote interface

The origin of eukaryotes remains unclear1-4. Current data suggest that eukaryotes may have emerged from an archaeal lineage known as 'Asgard' archaea5,6. Despite the eukaryote-like genomic features that are found in these archaea, the evolutionary transition from archaea to eukaryotes remains unclear, owing to the lack of cultured representatives and corresponding physiological insights. Here we report the decade-long isolation of an Asgard archaeon related to Lokiarchaeota from deep marine sediment. The archaeon-'Candidatus Prometheoarchaeum syntrophicum' strain MK-D1-is an anaerobic, extremely slow-growing, small coccus (around 550 nm in diameter) that degrades amino acids through syntrophy. Although eukaryote-like intracellular complexes have been proposed for Asgard archaea6, the isolate has no visible organelle-like structure. Instead, Ca. P. syntrophicum is morphologically complex and has unique protrusions that are long and often branching. On the basis of the available data obtained from cultivation and genomics, and reasoned interpretations of the existing literature, we propose a hypothetical model for eukaryogenesis, termed the entangle-engulf-endogenize (also known as E3) model.

The persimmon (Diospyros oleifera Cheng) genome provides new insights into the inheritance of astringency and ancestral evolution

Persimmon (Diospyros kaki) is an oriental perennial woody fruit tree whose popular fruit is produced and consumed worldwide. The persimmon fruit is unique because of the hyperaccumulation of proanthocyanidins during fruit development, causing the mature fruit of most cultivars to have an astringent taste. In this study, we obtained a chromosome-scale genome assembly for 'Youshi' (Diospyros oleifera, 2n = 2x = 30), the diploid species of persimmon, by integrating Illumina sequencing, single-molecule real-time sequencing, and high-throughput chromosome conformation capture techniques. The assembled D. oleifera genome consisted of 849.53 Mb, 94.14% (799.71 Mb) of which was assigned to 15 pseudochromosomes, and is the first assembled genome for any member of the Ebenaceae. Comparative genomic analysis revealed that the D. oleifera genome underwent an ancient γ whole-genome duplication event. We studied the potential genetic basis for astringency development (proanthocyanidin biosynthesis) and removal (proanthocyanidin insolublization). Proanthocyanidin biosynthesis genes were mainly distributed on chromosome 1, and the clustering of these genes is responsible for the genetic stability of astringency heredity. Genome-based RNA-seq identified deastringency genes, and promoter analysis showed that most of their promoters contained large numbers of low oxygen-responsive motifs, which is consistent with the efficient industrial application of high CO2 treatment to remove astringency. Using the D. oleifera genome as the reference, SLAF-seq indicated that 'Youshi' is one of the ancestors of the cultivated persimmon (2n = 6x = 90). Our study provides significant insights into the genetic basis of persimmon evolution and the development and removal astringency, and it will facilitate the improvement of the breeding of persimmon fruit.

Draft Genome Sequences of 12 Monophasic Salmonella enterica subsp. enterica Serotype Typhimurium 1,4,[5],12:i:- Strains Isolated from Wild Griffon Vultures in Eastern Spain

Monophasic Salmonella enterica subsp. enterica serovar Typhimurium is one of the most common zoonotic pathogens. Salmonella species reside in a wide variety of hosts, including wild animals. Thus, we report here the genome sequences of 12 monophasic S. Typhimurium strains isolated from healthy wild vultures to gain better insight into their epidemiology and host-pathogen interactions.

Genomic Analysis of Pseudomonas sp. Strain SCT, an Iodate-Reducing Bacterium Isolated from Marine Sediment, Reveals a Possible Use for Bioremediation

Strain SCT is an iodate-reducing bacterium isolated from marine sediment in Kanagawa Prefecture, Japan. In this study, we determined the draft genome sequence of strain SCT and compared it to complete genome sequences of other closely related bacteria, including Pseudomonas stutzeri. A phylogeny inferred from concatenation of core genes revealed that strain SCT was closely related to marine isolates of P. stutzeri. Genes present in the SCT genome but absent from the other analyzed P. stutzeri genomes comprised clusters corresponding to putative prophage regions and possible operons. They included pil genes, which encode type IV pili for natural transformation; the mer operon, which encodes resistance systems for mercury; and the pst operon, which encodes a Pi-specific transport system for phosphate uptake. We found that strain SCT had more prophage-like genes than the other P. stutzeri strains and that the majority (70%) of them were SCT strain-specific. These genes, encoded on distinct prophage regions, may have been acquired after branching from a common ancestor following independent phage transfer events. Thus, the genome sequence of Pseudomonas sp. strain SCT can provide detailed insights into its metabolic potential and the evolution of genetic elements associated with its unique phenotype.

Long-term investigation of microbial community composition and transcription patterns in a biogas plant undergoing ammonia crisis

Ammonia caused disturbance of biogas production is one of the most frequent incidents in regular operation of biogas reactors. This study provides a detailed insight into the microbial community of a mesophilic, full-scale biogas reactor (477 kWh h-1 ) fed with maize silage, dried poultry manure and cow manure undergoing initial process disturbance by increased ammonia concentration. Over a time period of 587 days, the microbial community of the reactor was regularly monitored on a monthly basis by high-throughput amplicon sequencing of the archaeal and bacterial 16S rRNA genes. During this sampling period, the total ammonia concentrations varied between 2.7 and 5.8 g l-1 [NH4 + -N]. To gain further inside into the active metabolic pathways, for selected time points metatranscriptomic shotgun analysis was performed allowing the quantification of marker genes for methanogenesis, hydrolysis and syntrophic interactions. The results obtained demonstrated a microbial community typical for a mesophilic biogas plant. However in response to the observed changing process conditions (e.g. increasing NH4 + levels, changing feedstock composition), the microbial community reacted highly flexible by changing and adapting the community composition. The Methanosarcina-dominated archaeal community was shifted to a Methanomicrobiales-dominated archaeal community in the presence of increased ammonia conditions. A similar trend as in the phylogenetic composition was observed in the transcription activity of genes coding for enzymes involved in acetoclastic methanogenesis and syntrophic acetate oxidations (Codh/Acs and Fthfs). In accordance, Clostridia simultaneously increased under elevated ammonia concentrations in abundance and were identified as the primary syntrophic interaction partner with the now Methanomicrobiales-dominated archaeal community. In conclusion, overall stable process performance was maintained during increased ammonia concentration in the studied reactor based on the microbial communities' ability to flexibly respond by reorganizing the community composition while remaining functionally stable.

GROOLS: reactive graph reasoning for genome annotation through biological processes

BACKGROUND

High quality functional annotation is essential for understanding the phenotypic consequences encoded in a genome. Despite improvements in bioinformatics methods, millions of sequences in databanks are not assigned reliable functions. The curation of protein functions in the context of biological processes is a way to evaluate and improve their annotation.

RESULTS

We developed an expert system using paraconsistent logic, named GROOLS (Genomic Rule Object-Oriented Logic System), that evaluates the completeness and the consistency of predicted functions through biological processes like metabolic pathways. Using a generic and hierarchical representation of knowledge, biological processes are modeled in a graph from which observations (i.e. predictions and expectations) are propagated by rules. At the end of the reasoning, conclusions are assigned to biological process components and highlight uncertainties and inconsistencies. Results on 14 microbial organisms are presented.

CONCLUSIONS

GROOLS software is designed to evaluate the overall accuracy of functional unit and pathway predictions according to organism experimental data like growth phenotypes. It assists biocurators in the functional annotation of proteins by focusing on missing or contradictory observations.

A Practical Guide for Comparative Genomics of Mobile Genetic Elements in Prokaryotic Genomes

Mobile genetic elements (MGEs) are an important feature of prokaryote genomes but are seldom well annotated and, consequently, are often underestimated. MGEs include transposons (Tn), insertion sequences (ISs), prophages, genomic islands (GEIs), integrons, and integrative and conjugative elements (ICEs). They are intimately involved in genome evolution and promote phenomena such as genomic expansion and rearrangement, emergence of virulence and pathogenicity, and symbiosis. In spite of the annotation bottleneck, there are so far at least 75 different programs and databases dedicated to prokaryotic MGE analysis and annotation, and this number is rapidly growing. Here, we present a practical guide to explore, compare, and visualize prokaryote MGEs using a combination of available software and databases tailored to small scale genome analyses. This protocol can be coupled with expert MGE annotation and exploited for evolutionary and comparative genomic analyses.

Proteomic characterization and bio-informatic analysis of differentially expressed E. coli Nissle 1917 proteins with response to cocoti wine stress

The present study emphases the comparative proteomic analysis of Escherichia coli Nissle 1917 under cocoti palm wine stress and identified differentially expressed proteins. Protein samples were analyzed by 2-D, MALDI-TOF combined with MS access. In 2-D electrophoresis, eight differentially expressed proteins were identified: five up-regulated, two down-regulated and one newly expressed protein. Protein spots were digested with trypsin for MALDI-TOF-MS analysis; protein sequences were obtained from MASCOT search. Sequences were aligned with template using Swiss Model server. Phyre-2 was used to predict homology modeling, RasMol was used to analyze the modeling structures, PSVS server was utilized to validate the protein structure by Ramachandran's plot analysis, physical and chemical properties were analyzed using ProtParam server, Phylogenetic tree was constructed by Mega4. UniProt search helps to find protein functional information of differentially expressed proteins, involved in catalytic activities, regulation mechanisms, DNA damage stimulus, anti-termination and termination process, protein binding, electron transport mechanism, and cell signaling process functions. A detailed exploration of the proteins under cocoti palm wine stress have provided the composition, structure and functions of the expressed proteins for further investigation.

Draft Genome Sequence of Roseovarius sp. A-2, an Iodide-Oxidizing Bacterium Isolated from Natural Gas Brine Water, Chiba, Japan

Roseovarius sp. A-2 is a heterotrophic iodide (I^-)-oxidizing bacterium isolated from iodide-rich natural gas brine water in Chiba, Japan. This strain oxidizes iodide to molecular iodine (I₂) by means of an extracellular multicopper oxidase. Here we report the draft genome sequence of strain A-2. The draft genome contained 46 tRNA genes, 1 copy of a 16S-23S-5S rRNA operon, and 4,514 protein coding DNA sequences, of which 1,207 (27%) were hypothetical proteins. The genome contained a gene encoding IoxA, a multicopper oxidase previously found to catalyze the oxidation of iodide in Iodidimonas sp. Q-1. This draft genome provides detailed insights into the metabolism and potential application of Roseovarius sp. A-2.

The landscape of microbial phenotypic traits and associated genes

Bacteria and Archaea display a variety of phenotypic traits and can adapt to diverse ecological niches. However, systematic annotation of prokaryotic phenotypes is lacking. We have therefore developed ProTraits, a resource containing ∼545 000 novel phenotype inferences, spanning 424 traits assigned to 3046 bacterial and archaeal species. These annotations were assigned by a computational pipeline that associates microbes with phenotypes by text-mining the scientific literature and the broader World Wide Web, while also being able to define novel concepts from unstructured text. Moreover, the ProTraits pipeline assigns phenotypes by drawing extensively on comparative genomics, capturing patterns in gene repertoires, codon usage biases, proteome composition and co-occurrence in metagenomes. Notably, we find that gene synteny is highly predictive of many phenotypes, and highlight examples of gene neighborhoods associated with spore-forming ability. A global analysis of trait interrelatedness outlined clusters in the microbial phenotype network, suggesting common genetic underpinnings. Our extended set of phenotype annotations allows detection of 57 088 high confidence gene-trait links, which recover many known associations involving sporulation, flagella, catalase activity, aerobicity, photosynthesis and other traits. Over 99% of the commonly occurring gene families are involved in genetic interactions conditional on at least one phenotype, suggesting that epistasis has a major role in shaping microbial gene content.

Draft Genome Sequence of Arenibacter sp. Strain C-21, an Iodine-Accumulating Bacterium Isolated from Surface Marine Sediment

Arenibacter sp. strain C-21, isolated from surface marine sediment of Japan, accumulates iodine in the presence of glucose and iodide (I^-). We report here the draft genome sequence of this strain to provide insight into the molecular mechanism underlying its iodine-accumulating ability.

Global transcriptional analysis suggests Lasiodiplodia theobromae pathogenicity factors involved in modulation of grapevine defensive response

BACKGROUND

Lasiodiplodia theobromae is a fungus of the Botryosphaeriaceae that causes grapevine vascular disease, especially in regions with hot climates. Fungi in this group often remain latent within their host and become virulent under abiotic stress. Transcriptional regulation analysis of L. theobromae exposed to heat stress (HS) was first carried out in vitro in the presence of grapevine wood (GW) to identify potential pathogenicity genes that were later evaluated for in planta expression.

RESULTS

A total of 19,860 de novo assembled transcripts were obtained, forty-nine per cent of which showed homology to the Botryosphaeriaceae fungi, Neofusicoccum parvum or Macrophomina phaseolina. Three hundred ninety-nine have homology with genes involved in pathogenic processes and several belonged to expanded gene families in others fungal grapevine vascular pathogens. Gene expression analysis showed changes in fungal metabolism of phenolic compounds; where genes encoding for enzymes, with the ability to degrade salicylic acid (SA) and plant phenylpropanoid precursors, were up-regulated during in vitro HS response, in the presence of GW. These results suggest that the fungal L-tyrosine catabolism pathway could help the fungus to remove phenylpropanoid precursors thereby evading the host defense response. The in planta up-regulation of salicylate hydroxylase, intradiol ring cleavage dioxygenase and fumarylacetoacetase encoding genes, further supported this hypothesis. Those genes were even more up-regulated in HS-stressed plants, suggesting that fungus takes advantage of the increased phenylpropanoid precursors produced under stress. Pectate lyase was up-regulated while a putative amylase was down-regulated in planta, this could be associated with an intercellular growth strategy during the first stages of colonization.

CONCLUSIONS

L. theobromae transcriptome was established and validated. Its usefulness was demonstrated through the identification of genes expressed during the infection process. Our results support the hypothesis that heat stress facilitates fungal colonization, because of the fungus ability to use the phenylpropanoid precursors and SA, both compounds known to control host defense.

The use of matrix coating assisted by an electric field (MCAEF) to enhance mass spectrometric imaging of human prostate cancer biomarkers

In this work, we combined a newly developed matrix coating technique - matrix coating assisted by an electric field (MCAEF) and matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) to enhance the imaging of peptides and proteins in tissue specimens of human prostate cancer. MCAEF increased the signal-to-noise ratios of the detected proteins by a factor of 2 to 5, and 232 signals were detected within the m/z 3500-37500 mass range on a time-of-flight mass spectrometer and with the sinapinic acid MALDI matrix. Among these species, three proteins (S100-A9, S100-A10, and S100-A12) were only observed in the cancerous cell region and 14 proteins, including a fragment of mitogen-activated protein kinase/extracellular signal-regulated kinase kinase kinase 2, a fragment of cAMP-regulated phosphoprotein 19, 3 apolipoproteins (C-I, A-I, and A-II), 2 S100 proteins (A6 and A8), β-microseminoprotein, tumor protein D52, α-1-acid glycoprotein 1, heat shock protein β-1, prostate-specific antigen, and 2 unidentified large peptides at m/z 5002.2 and 6704.2, showed significantly differential distributions at the p < 0.05 (t-test) level between the cancerous and the noncancerous regions of the tissue. Among these 17 species, the distributions of apolipoprotein C-I, S100-A6, and S100-A8 were verified by immunohistological staining. In summary, this study resulted in the imaging of the largest group of proteins in prostate cancer tissues by MALDI-MS reported thus far, and is the first to show a correlation between S100 proteins and prostate cancer in a MS imaging study. The successful imaging of the three proteins only found in the cancerous tissues, as well as those showing differential expressions demonstrated the potential of MCAEF-MALDI/MS for the in situ detection of potential cancer biomarkers. Copyright © 2015 John Wiley & Sons, Ltd.

Bacterial ferrous iron transport: the Feo system

To maintain iron homeostasis within the cell, bacteria have evolved various types of iron acquisition systems. Ferric iron (Fe(3+)) is the dominant species in an oxygenated environment, while ferrous iron (Fe(2+)) is more abundant under anaerobic conditions or at low pH. For organisms that must combat oxygen limitation for their everyday survival, pathways for the uptake of ferrous iron are essential. Several bacterial ferrous iron transport systems have been described; however, only the Feo system appears to be widely distributed and is exclusively dedicated to the transport of iron. In recent years, many studies have explored the role of the FeoB and FeoA proteins in ferrous iron transport and their contribution toward bacterial virulence. The three-dimensional structures for the Feo proteins have recently been determined and provide insight into the molecular details of the transport system. A highly select group of bacteria also express the FeoC protein from the same operon. This review will provide a comprehensive look at the structural and functional aspects of the Feo system. In addition, bioinformatics analyses of the feo operon and the Feo proteins have been performed to complement our understanding of this ubiquitous bacterial uptake system, providing a new outlook for future studies.

Complete genome sequencing of Dehalococcoides sp. strain UCH007 using a differential reads picking method

A novel Dehalococcoides sp. strain UCH007 was isolated from the groundwater polluted with chlorinated ethenes in Japan. This strain is capable of dechlorinating trichloroethene, cis-1,2-dichloroethene and vinyl chloride to ethene. Dehalococcoides bacteria are hardly cultivable, so genome sequencing has presented a challenge. In this study, we developed a differential reads picking method for mixed genomic DNA obtained from a co-culture, and applied it to the sequencing of strain UCH007. The genome of strain UCH007 consists of a 1,473,548-bp chromosome that encodes 1509 coding sequences including 29 putative reductive dehalogenase genes. Strain UCH007 is the first strain in the Victoria subgroup found to possess the pceA, tceA and vcrA genes.

Functional classification of CATH superfamilies: a domain-based approach for protein function annotation

Motivation: Computational approaches that can predict protein functions are essential to bridge the widening function annotation gap especially since <1.0% of all proteins in UniProtKB have been experimentally characterized. We present a domain-based method for protein function classification and prediction of functional sites that exploits functional sub-classification of CATH superfamilies. The superfamilies are sub-classified into functional families (FunFams) using a hierarchical clustering algorithm supervised by a new classification method, FunFHMMer.

Results: FunFHMMer generates more functionally coherent groupings of protein sequences than other domain-based protein classifications. This has been validated using known functional information. The conserved positions predicted by the FunFams are also found to be enriched in known functional residues. Moreover, the functional annotations provided by the FunFams are found to be more precise than other domain-based resources. FunFHMMer currently identifies 110 439 FunFams in 2735 superfamilies which can be used to functionally annotate > 16 million domain sequences.

Availability and implementation: All FunFam annotation data are made available through the CATH webpages (http://www.cathdb.info). The FunFHMMer webserver (http://www.cathdb.info/search/by_funfhmmer) allows users to submit query sequences for assignment to a CATH FunFam.

Contact:sayoni.das.12@ucl.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

The functional landscape bound to the transcription factors of Escherichia coli K-12

Motivated by the experimental evidences accumulated in the last ten years and based on information deposited in RegulonDB, literature look up, and sequence analysis, we analyze the repertoire of 304 DNA-binding Transcription factors (TFs) in Escherichia coli K-12. These regulators were grouped in 78 evolutionary families and are regulating almost half of the total genes in this bacterium. In structural terms, 60% of TFs are composed by two-domains, 30% are monodomain, and 10% three- and four-structural domains. As previously noticed, the most abundant DNA-binding domain corresponds to the winged helix-turn-helix, with few alternative DNA-binding structures, resembling the hypothesis of successful protein structures with the emergence of new ones at low scales. In summary, we identified and described the characteristics associated to the DNA-binding TF in E. coli K-12. We also identified twelve functional modules based on a co-regulated gene matrix. Finally, diverse regulons were predicted based on direct associations between the TFs and potential regulated genes. This analysis should increase our knowledge about the gene regulation in the bacterium E. coli K-12, and provide more additional clues for comprehensive modelling of transcriptional regulatory networks in other bacteria.

Complete genome sequence of Selenomonas ruminantium subsp. lactilytica will accelerate further understanding of the nature of the class Negativicutes

Selenomonas ruminantium subsp. lactilytica, a strictly anaerobic ruminal bacterium, possesses typical Gram-negative cell surface structure comprising cytoplasmic membrane, peptidoglycan layer and outer membrane, whereas its 16S rRNA-based taxonomy shows that the bacteria belongs to Gram-positive Firmicutes. Complete genome analysis showed that genes or gene clusters involved in Gram-negative cell structure were scattered in the S. ruminantium genome, and might provide the new insight of phylogenetic relationship between the bacterium and other bacterial species.

An integrative computational model for large-scale identification of metalloproteins in microbial genomes: a focus on iron-sulfur cluster proteins

Metalloproteins represent a ubiquitous group of molecules which are crucial to the survival of all living organisms. While several metal-binding motifs have been defined, it remains challenging to confidently identify metalloproteins from primary protein sequences using computational approaches alone. Here, we describe a comprehensive strategy based on a machine learning approach to design and assess a penalized generalized linear model. We used this strategy to detect members of the iron-sulfur cluster protein family. A new category of descriptors, whose profile is based on profile hidden Markov models, encoding structural information was combined with public descriptors into a linear model. The model was trained and tested on distinct datasets composed of well-characterized iron-sulfur protein sequences, and the resulting model provided higher sensitivity compared to a motif-based approach, while maintaining a good level of specificity. Analysis of this linear model allows us to detect and quantify the contribution of each descriptor, providing us with a better understanding of this complex protein family along with valuable indications for further experimental characterization. Two newly-identified proteins, YhcC and YdiJ, were functionally validated as genuine iron-sulfur proteins, confirming the prediction. The computational model was then applied to over 550 prokaryotic genomes to screen for iron-sulfur proteomes; the results are publicly available at: . This study represents a proof-of-concept for the application of a penalized linear model to identify metalloprotein superfamilies on a large-scale. The application employed here, screening for iron-sulfur proteomes, provides new candidates for further biochemical and structural analysis as well as new resources for an extensive exploration of iron-sulfuromes in the microbial world.

Gene loss rather than gene gain is associated with a host jump from monocots to dicots in the Smut Fungus Melanopsichium pennsylvanicum

Smut fungi are well-suited to investigate the ecology and evolution of plant pathogens, as they are strictly biotrophic, yet cultivable on media. Here we report the genome sequence of Melanopsichium pennsylvanicum, closely related to Ustilago maydis and other Poaceae-infecting smuts, but parasitic to a dicot plant. To explore the evolutionary patterns resulting from host adaptation after this huge host jump, the genome of Me. pennsylvanicum was sequenced and compared with the genomes of U. maydis, Sporisorium reilianum, and U. hordei. Although all four genomes had a similar completeness in CEGMA (Core Eukaryotic Genes Mapping Approach) analysis, gene absence was highest in Me. pennsylvanicum, and most pronounced in putative secreted proteins, which are often considered as effector candidates. In contrast, the amount of private genes was similar among the species, highlighting that gene loss rather than gene gain is the hallmark of adaptation after the host jump to the dicot host. Our analyses revealed a trend of putative effectors to be next to another putative effector, but the majority of these are not in clusters and thus the focus on pathogenicity clusters might not be appropriate for all smut genomes. Positive selection studies revealed that Me. pennsylvanicum has the highest number and proportion of genes under positive selection. In general, putative effectors showed a higher proportion of positively selected genes than noneffector candidates. The 248 putative secreted effectors found in all four smut genomes might constitute a core set needed for pathogenicity, whereas those 92 that are found in all grass-parasitic smuts but have no ortholog in Me. pennsylvanicum might constitute a set of effectors important for successful colonization of grass hosts.

In search of genome annotation consistency: solid gene clusters and how to use them

Maintaining consistency in genome annotations is important for supporting many computational tasks, particularly metabolic modeling. The SEED project has implemented a process that improves annotation consistencies across microbial genomes for proteins with conserved sequences and genomic context. In this research report, we describe this process and show how this effort has resulted in improvements to microbial genome annotations in the SEED. We also compare SEED annotation consistencies with other commonly used resources such as IMG (the Joint Genome Institute’s Integrated Microbial Genomes system), RefSeq (the National Center for Biotechnology Information’s Reference Sequence Database), Swiss-Prot (the annotated protein sequence database of the Swiss Institute of Bioinformatics, European Molecular Biology Laboratory and the European Bioinformatics Institute) and TrEMBL (Translated European Molecular Biology Laboratory nucleotide sequence data Library). Our analysis indicates that manual and computational efforts are paying off for the databases where consistency is a major goal.

ENZYMAP: exploiting protein annotation for modeling and predicting EC number changes in UniProt/Swiss-Prot

The volume and diversity of biological data are increasing at very high rates. Vast amounts of protein sequences and structures, protein and genetic interactions and phenotype studies have been produced. The majority of data generated by high-throughput devices is automatically annotated because manually annotating them is not possible. Thus, efficient and precise automatic annotation methods are required to ensure the quality and reliability of both the biological data and associated annotations. We proposed ENZYMatic Annotation Predictor (ENZYMAP), a technique to characterize and predict EC number changes based on annotations from UniProt/Swiss-Prot using a supervised learning approach. We evaluated ENZYMAP experimentally, using test data sets from both UniProt/Swiss-Prot and UniProt/TrEMBL, and showed that predicting EC changes using selected types of annotation is possible. Finally, we compared ENZYMAP and DETECT with respect to their predictions and checked both against the UniProt/Swiss-Prot annotations. ENZYMAP was shown to be more accurate than DETECT, coming closer to the actual changes in UniProt/Swiss-Prot. Our proposal is intended to be an automatic complementary method (that can be used together with other techniques like the ones based on protein sequence and structure) that helps to improve the quality and reliability of enzyme annotations over time, suggesting possible corrections, anticipating annotation changes and propagating the implicit knowledge for the whole dataset.