In this study, we constructed a protein-protein interaction (PPI) network of B. licheniformis strain WX-02 with interolog method and domain-based method, which contained 15,864 edges and 2,448 nodes. Although computationally predicted networks have relatively low coverage and high false-positive rate, our prediction was confirmed from three perspectives: local structural features, functional similarities and transcriptional correlations. Further analysis of the COG heat map showed that protein interactions in B. licheniformis WX-02 mainly occurred in the same functional categories. By incorporating the transcriptome data, we found that the topological properties of the PPI network were robust under normal and high salt conditions. In addition, 267 different protein complexes were identified and 117 poorly characterized proteins were annotated with certain functions based on the PPI network. Furthermore, the sub-network showed that a hub protein CcpA jointed directly or indirectly many proteins related to γ-PGA synthesis and regulation, such as PgsB, GltA, GltB, ProB, ProJ, YcgM and two signal transduction systems ComP-ComA and DegS-DegU. Thus, CcpA might play an important role in the regulation of γ-PGA synthesis. This study therefore will facilitate the understanding of the complex cellular behaviors and mechanisms of γ-PGA synthesis in B. licheniformis WX-02.

Spiroplasmas are helical and motile members of a cell wall-less eubacterial group called Mollicutes. Although all spiroplasmas are associated with arthropods, they exhibit great diversity with respect to both their modes of transmission and their effects on their hosts; ranging from horizontally transmitted pathogens and commensals to endosymbionts that are transmitted transovarially (i.e., from mother to offspring). Here we provide the first genome sequence, along with proteomic validation, of an endosymbiotic inherited Spiroplasma bacterium, the Spiroplasma poulsonii MSRO strain harbored by Drosophila melanogaster. Comparison of the genome content of S. poulsonii with that of horizontally transmitted spiroplasmas indicates that S. poulsonii has lost many metabolic pathways and transporters, demonstrating a high level of interdependence with its insect host. Consistent with genome analysis, experimental studies showed that S. poulsonii metabolizes glucose but not trehalose. Notably, trehalose is more abundant than glucose in Drosophila hemolymph, and the inability to metabolize trehalose may prevent S. poulsonii from overproliferating. Our study identifies putative virulence genes, notably, those for a chitinase, the H2O2-producing glycerol-3-phosphate oxidase, and enzymes involved in the synthesis of the eukaryote-toxic lipid cardiolipin. S. poulsonii also expresses on the cell membrane one functional adhesion-related protein and two divergent spiralin proteins that have been implicated in insect cell invasion in other spiroplasmas. These lipoproteins may be involved in the colonization of the Drosophila germ line, ensuring S. poulsonii vertical transmission. The S. poulsonii genome is a valuable resource to explore the mechanisms of male killing and symbiont-mediated protection, two cardinal features of many facultative endosymbionts.

Pseudomonas fluorescens PCL1751 is a rod-shaped Gram-negative bacterium isolated from the rhizosphere of a greenhouse-grown tomato plant in Uzbekistan. It controls several plant root diseases caused by Fusarium fungi through the mechanism of competition for nutrients and niches (CNN). This mechanism does not rely on the production of antibiotics, so it avoids the concerns of resistance development and is environmentally safe. Additionally, this bacterium promotes plant growth by alleviating salt stress for its plant host. To investigate the genetic mechanisms that may explain these observations, we determined the complete genome sequence of this bacterium, examined its gene content, and performed comparative genomics analysis with other Pseudomonas strains. The genome of P. fluorescens PCL1751 consisted of one circular chromosome that is 6,143,950 base-pairs (bp) in size; no plasmid was found. The annotation included 19 rRNA, 70 tRNA, and 5,534 protein-coding genes. The gene content analysis identified a large number of genes involved in chemotaxis and motility, colonization of the rhizosphere, siderophore biosynthesis, and osmoprotectant production. In contrast, the pathways involved in the biosynthesis of phytohormones or antibiotics were not found. Comparison with other Pseudomonas genomes revealed extensive variations in their genome size and gene content. The presence and absence of secretion system genes were highly variable. As expected, the synteny conservation among strains decreased as a function of phylogenetic divergence. The integration of prophages appeared to be an important driver for genome rearrangements. The whole-genome gene content analysis of this plant growth-promoting rhizobacterium (PGPR) provided some genetic explanations to its phenotypic characteristics. The extensive and versatile substrate utilization pathways, together with the presence of many genes involved in competitive root colonization, provided further support for the finding that this strain achieves biological control of pathogens through effective competition for nutrients and niches.

We reconstructed the first genome-scale metabolic network of the plant pathogen Pectobacterium carotovorum subsp. carotovorum PC1 based on its genomic sequence, annotation, and physiological data. Metabolic characteristics were analyzed using flux balance analysis (FBA), and the results were afterwards validated by phenotype microarray (PM) experiments. The reconstructed genome-scale metabolic model, iPC1209, contains 2235 reactions, 1113 metabolites and 1209 genes. We identified 19 potential bactericide targets through a comprehensive in silico gene-deletion study. Next, we performed virtual screening to identify candidate inhibitors for an important potential drug target, alkaline phosphatase, and experimentally verified that three lead compounds were able to inhibit both bacterial cell viability and the activity of alkaline phosphatase in vitro. This study illustrates a new strategy for the discovery of agricultural bactericides.

We describe the discovery of sno-lncRNAs, a class of nuclear-enriched intron-derived long noncoding RNAs (lncRNAs) that are processed on both ends by the snoRNA machinery. During exonucleolytic trimming, the sequences between the snoRNAs are not degraded, leading to the accumulation of lncRNAs flanked by snoRNA sequences but lacking 5' caps and 3' poly(A) tails. Such RNAs are widely expressed in cells and tissues and can be produced by either box C/D or box H/ACA snoRNAs. Importantly, the genomic region encoding one abundant class of sno-lncRNAs (15q11-q13) is specifically deleted in Prader-Willi Syndrome (PWS). The PWS region sno-lncRNAs do not colocalize with nucleoli or Cajal bodies, but rather accumulate near their sites of synthesis. These sno-lncRNAs associate strongly with Fox family splicing regulators and alter patterns of splicing. These results thus implicate a previously unannotated class of lncRNAs in the molecular pathogenesis of PWS.

Citrus is one of the most important and widely grown fruit crop with global production ranking firstly among all the fruit crops in the world. Sweet orange accounts for more than half of the Citrus production both in fresh fruit and processed juice. We have sequenced the draft genome of a double-haploid sweet orange (C. sinensis cv. Valencia), and constructed the Citrus sinensis annotation project (CAP) to store and visualize the sequenced genomic and transcriptome data. CAP provides GBrowse-based organization of sweet orange genomic data, which integrates ab initio gene prediction, EST, RNA-seq and RNA-paired end tag (RNA-PET) evidence-based gene annotation. Furthermore, we provide a user-friendly web interface to show the predicted protein-protein interactions (PPIs) and metabolic pathways in sweet orange. CAP provides comprehensive information beneficial to the researchers of sweet orange and other woody plants, which is freely available at http://citrus.hzau.edu.cn/.

Natural antisense transcripts (NATs), as one type of regulatory RNAs, occur prevalently in plant genomes and play significant roles in physiological and pathological processes. Although their important biological functions have been reported widely, a comprehensive database is lacking up to now. Consequently, we constructed a plant NAT database (PlantNATsDB) involving approximately 2 million NAT pairs in 69 plant species. GO annotation and high-throughput small RNA sequencing data currently available were integrated to investigate the biological function of NATs. PlantNATsDB provides various user-friendly web interfaces to facilitate the presentation of NATs and an integrated, graphical network browser to display the complex networks formed by different NATs. Moreover, a 'Gene Set Analysis' module based on GO annotation was designed to dig out the statistical significantly overrepresented GO categories from the specific NAT network. PlantNATsDB is currently the most comprehensive resource of NATs in the plant kingdom, which can serve as a reference database to investigate the regulatory function of NATs. The PlantNATsDB is freely available at http://bis.zju.edu.cn/pnatdb/.

Mosquitoes are hosts of several Spiroplasma species that belong to different serogroups. To investigate the genetic mechanisms that may be involved in the utilization of similar hosts in these phylogenetically distinct bacteria, we determined the complete genome sequences of Spiroplasma diminutum and S. taiwanense for comparative analysis. The genome alignment indicates that their chromosomal organization is highly conserved, which is in sharp contrast to the elevated genome instabilities observed in other Spiroplasma lineages. Examination of the substrate utilization strategies revealed that S. diminutum can use a wide range of carbohydrates, suggesting that it is well suited to living in the gut (and possibly the circulatory system) of its mosquito hosts. In comparison, S. taiwanense has lost several carbohydrate utilization genes and acquired additional sets of oligopeptide transporter genes through tandem duplications, suggesting that proteins from digested blood meal or lysed host cells may be an important nutrient source. Moreover, one glycerol-3-phosphate oxidase gene (glpO) was found in S. taiwanense but not S. diminutum. This gene is linked to the production of reactive oxygen species and has been shown to be a major virulence factor in Mycoplasma mycoides. This finding may explain the pathogenicity of S. taiwanense observed in previous artificial infection experiments, while no apparent effect was found for S. diminutum. To infer the gene content evolution at deeper divergence levels, we incorporated other Mollicutes genomes for comparative analyses. The results suggest that the losses of biosynthetic pathways are a recurrent theme in these host-associated bacteria.

The genus Spiroplasma contains a group of helical, motile, and wall-less bacteria in the class Mollicutes. Similar to other members of this class, such as the animal-pathogenic Mycoplasma and the plant-pathogenic 'Candidatus Phytoplasma', all characterized Spiroplasma species were found to be associated with eukaryotic hosts. While most of the Spiroplasma species appeared to be harmless commensals of insects, a small number of species have evolved pathogenicity toward various arthropods and plants. In this study, we isolated a novel strain of honeybee-associated S. melliferum and investigated its genetic composition and evolutionary history by whole-genome shotgun sequencing and comparative analysis with other Mollicutes genomes.

Adenosine-to-inosine (A-to-I) RNA editing is recognized as a cellular mechanism for generating both RNA and protein diversity. Inosine base pairs with cytidine during reverse transcription and therefore appears as guanosine during sequencing of cDNA. Current approaches of RNA editing identification largely depend on the comparison between transcriptomes and genomic DNA (gDNA) sequencing datasets from the same individuals, and it has been challenging to identify editing candidates from transcriptomes in the absence of gDNA information.

Phytoplasmas and mycoplasmas are two groups of important pathogens in the bacterial class Mollicutes. Because of their economical and clinical importance, these obligate pathogens have attracted much research attention. However, difficulties involved in the empirical study of these bacteria, particularly the fact that phytoplasmas have not yet been successfully cultivated outside of their hosts despite decades of attempts, have greatly hampered research progress. With the rapid advancements in genome sequencing, comparative genome analysis provides a new approach to facilitate our understanding of these bacteria. In this study, our main focus is to investigate the evolution of gene content in phytoplasmas, mycoplasmas, and their common ancestor. By using a phylogenetic framework for comparative analysis of 12 complete genome sequences, we characterized the putative gains and losses of genes in these obligate parasites. Our results demonstrated that the degradation of metabolic capacities in these bacteria has occurred predominantly in the common ancestor of Mollicutes, prior to the evolutionary split of phytoplasmas and mycoplasmas. Furthermore, we identified a list of genes that are acquired by the common ancestor of phytoplasmas and are conserved across all strains with complete genome sequences available. These genes include several putative effectors for the interactions with hosts and may be good candidates for future functional characterization.

We devised software tools to systematically investigate the contents and contexts of bacterial tRNA and tmRNA genes, which are known insertion hotspots for genomic islands (GIs). The strategy, based on MAUVE-facilitated multigenome comparisons, was used to examine 87 Escherichia coli MG1655 tRNA and tmRNA genes and their orthologues in E.coli EDL933, E.coli CFT073 and Shigella flexneri Sf301. Our approach identified 49 GIs occupying approximately 1.7 Mb that mapped to 18 tRNA genes, missing 2 but identifying a further 30 GIs as compared with Islander [Y. Mantri and K. P. Williams (2004), Nucleic Acids Res., 32, D55-D58]. All these GIs had many strain-specific CDS, anomalous GC contents and/or significant dinucleotide biases, consistent with foreign origins. Our analysis demonstrated marked conservation of sequences flanking both empty tRNA sites and tRNA-associated GIs across all four genomes. Remarkably, there were only 2 upstream and 5 downstream deletions adjacent to the 328 loci investigated. In silico PCR analysis based on conserved flanking regions was also used to interrogate hotspots in another eight completely or partially sequenced E.coli and Shigella genomes. The tools developed are ideal for the analysis of other bacterial species and will lead to in silico and experimental discovery of new genomic islands.

The protistan phylum Apicomplexa contains many important pathogens and is the subject of intense genome sequencing efforts. Based upon the genome sequences from seven apicomplexan species and a ciliate outgroup, we identified 268 single-copy genes suitable for phylogenetic inference. Both concatenation and consensus approaches inferred the same species tree topology. This topology is consistent with most prior conceptions of apicomplexan evolution based upon ultrastructural and developmental characters, that is, the piroplasm genera Theileria and Babesia form the sister group to the Plasmodium species, the coccidian genera Eimeria and Toxoplasma are monophyletic and are the sister group to the Plasmodium species and piroplasm genera, and Cryptosporidium forms the sister group to the above mentioned with the ciliate Tetrahymena as the outgroup. The level of incongruence among gene trees appears to be high at first glance; only 19% of the genes support the species tree, and a total of 48 different gene-tree topologies are observed. Detailed investigations suggest that the low signal-to-noise ratio in many genes may be the main source of incongruence. The probability of being consistent with the species tree increases as a function of the minimum bootstrap support observed at tree nodes for a given gene tree. Moreover, gene sequences that generate high bootstrap support are robust to the changes in alignment parameters or phylogenetic method used. However, caution should be taken in that some genes can infer a "wrong" tree with strong support because of paralogy, model violations, or other causes. The importance of examining multiple, unlinked genes that possess a strong phylogenetic signal cannot be overstated.

Protein translation regulation has essential roles in inflammatory responses, cancer initiation and the pathogenesis of several neurodegenerative disorders. However, the role of the regulation of protein translation in mammalian skeleton development has been rarely elaborated. Here we report that the lack of the RNA-binding protein sterile alpha motif domain containing protein 4 (SAMD4) resulted in multiple developmental defects in mice, including delayed bone development and decreased osteogenesis. Samd4-deficient mesenchymal progenitors exhibit impaired osteoblast differentiation and function. Mechanism study demonstrates that SAMD4 binds the Mig6 mRNA and inhibits MIG6 protein synthesis. Consistent with this, Samd4-deficient cells have increased MIG6 protein level and knockdown of Mig6 rescues the impaired osteogenesis in Samd4-deficient cells. Furthermore, Samd4-deficient mice also display chondrocyte defects, which is consistent with the regulation of MIG6 protein level by SAMD4. These findings define SAMD4 as a previously unreported key regulator of osteoblastogenesis and bone development, implying that regulation of protein translation is an important mechanism governing skeletogenesis and that control of protein translation could have therapeutic potential in metabolic bone diseases, such as osteoporosis.

Seeds of the desert shrub, jojoba (Simmondsia chinensis), are an abundant, renewable source of liquid wax esters, which are valued additives in cosmetic products and industrial lubricants. Jojoba is relegated to its own taxonomic family, and there is little genetic information available to elucidate its phylogeny. Here, we report the high-quality, 887-Mb genome of jojoba assembled into 26 chromosomes with 23,490 protein-coding genes. The jojoba genome has only the whole-genome triplication (γ) shared among eudicots and no recent duplications. These genomic resources coupled with extensive transcriptome, proteome, and lipidome data helped to define heterogeneous pathways and machinery for lipid synthesis and storage, provided missing evolutionary history information for this taxonomically segregated dioecious plant species, and will support efforts to improve the agronomic properties of jojoba.

Visualizing the location and dynamics of RNAs in live cells is key to understanding their function. Here, we identify two endonuclease-deficient, single-component programmable RNA-guided and RNA-targeting Cas13 RNases (dCas13s) that allow robust real-time imaging and tracking of RNAs in live cells, even when using single 20- to 27-nt-long guide RNAs. Compared to the aptamer-based MS2-MCP strategy, an optimized dCas13 system is user friendly, does not require genetic manipulation, and achieves comparable RNA-labeling efficiency. We demonstrate that the dCas13 system is capable of labeling NEAT1, SatIII, MUC4, and GCN4 RNAs and allows the study of paraspeckle-associated NEAT1 dynamics. Applying orthogonal dCas13 proteins or combining dCas13 and MS2-MCP allows dual-color imaging of RNAs in single cells. Further combination of dCas13 and dCas9 systems allows simultaneous visualization of genomic DNA and RNA transcripts in living cells.

Phytoplasmas are uncultivated plant-pathogenic bacteria with agricultural importance. Those belonging to the 16SrII group, represented by 'Candidatus P. aurantifolia', have a wide range of plant hosts and cause significant yield losses in valuable crops, such as pear, sweet potato, peanut, and soybean. In this study, a method that combines immunoprecipitation-based enrichment and MinION long-read DNA sequencing was developed to solve the challenge of phytoplasma genome studies. This approach produced long reads with high mapping rates and high genomic coverage that can be combined with Illumina reads to produce complete genome assemblies with high accuracy. We applied this method to strain NCHU2014 and determined its complete genome sequence, which consists of one circular chromosome with 635,584 bp and one plasmid with 4,224 bp. Although 'Ca. P. aurantifolia' NCHU2014 has a small chromosome with only 471 protein-coding genes, it contains 33 transporter genes and 27 putative effector genes, which may contribute to obtaining nutrients from hosts and manipulating host developments for their survival and multiplication. Two effectors, the homologs of SAP11 and SAP54/PHYL1 identified in 'Ca. P. aurantifolia' NCHU2014, have the biochemical activities in destabilizing host transcription factors, which can explain the disease symptoms observed in infected plants. Taken together, this study provides the first complete genome available for the 16SrII phytoplasmas and contributes to the understanding of phytoplasma pathogenicity.

Most circular RNAs are produced from the back-splicing of exons of precursor mRNAs. Recent technological advances have in part overcome problems with their circular conformation and sequence overlap with linear cognate mRNAs, allowing a better understanding of their cellular roles. Depending on their localization and specific interactions with DNA, RNA, and proteins, circular RNAs can modulate transcription and splicing, regulate stability and translation of cytoplasmic mRNAs, interfere with signaling pathways, and serve as templates for translation in different biological and pathophysiological contexts. Emerging applications of RNA circles to interfere with cellular processes, modulate immune responses, and direct translation into proteins shed new light on biomedical research. In this review, we discuss approaches used in circular RNA studies and the current understanding of their regulatory roles and potential applications.

Understanding gene transcription and mRNA-protein (mRNP) dynamics in single cells in a multicellular organism has been challenging. The catalytically dead CRISPR-Cas13 (dCas13) system has been used to visualize RNAs in live cells without genetic manipulation. We optimize this system to track developmentally expressed mRNAs in zebrafish embryos and to understand features of endogenous transcription kinetics and mRNP export.

Long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) can play important roles in many biological processes. However, no study of the influence of epigenetics factors or the 3D structure of the genome in their regulation is available in plants.