Homology-based Searches Research Articles

Functional characterization of four antenna-enriched odorant binding proteins in Rhaphuma horsfieldi reveals the importance of RhorOBP1 in odorant recognition and insecticide resistance

The cerambycid beetles are key players for the sustenance of biodiversity in the forest ecosystem, but in most cases are well known due to their harmfulness to agricultural and forest plants. Here, we characterized the odorant binding protein (OBP) gene family in Rhaphuma horsfieldi, emphasizing the roles of RhorOBP1 in odorant reception and insecticide sequestering. A homology-based search led to the identification of 35 RhorOBP genes with a major distribution in the Minus-C OBPs clade (21/35 genes). Expression profiles showed that RhorOBP1–RhorOBP4 had the abundant expression in antennae. Binding assays revealed that the four RhorOBPs exhibited diverse odorant response profiles tuned differentially to various classes of plant odorants, comprising walnut-derived host volatiles and ordinary floral scents. Two broadly tuned RhorOBP1 and RhorOBP2 exhibited different chain length-dependent binding properties to 10C12C alcohols, aldehydes or acetates. Compared with other three proteins, RhorOBP1 reduced the binding to ligands with high affinities at pH 5.0 (1.27–6.72-fold differences relative to pH 7.4). Molecular docking and point-mutation experiments confirmed that Ser107, Tyr118, Tyr119 and Phe120 situated in the binding pocket of RhorOBP1 were critical determinants for the recognition of 14, 15, 10 and 10 compounds, respectively. On the other hand, RhorOBP1 could strongly bind six insecticides, particularly chlorpyrifos (dissociation constant, Ki = 3.69 ± 0.74 μM). This study has provided insights into different binding properties of four antenna-enriched RhorOBPs in R. horsfieldi and identifies a dual role of RhorOBP1 in the binding of odorants and insecticides.

Bioinformatic Prohormone Discovery in Basal Metazoans: Insights from Trichoplax.

Experimental discovery of neuropeptides and peptide hormones is a long and tedious task. Mining the genomic and transcriptomic sequence data with robust secretory peptide prediction tools can significantly facilitate subsequent experiments. We describe the application of various in silico neuropeptide discovery methods for the placozoan Trichopax adhaerens as an illustrated example and a powerful experimental paradigm for cellular and evolutionary biology. In total, 33 placozoan (neuro)peptide-like hormone precursors were found using homology-based BLAST search and repeat-based and comparative evolutionary methods. Some of the discovered precursors are homologous to insulins and RFamide precursors from Cnidaria and other animal phyla.

The interaction networks of small rubber particle proteins in the latex of Taraxacum koksaghyz reveal diverse functions in stress responses and secondary metabolism.

The Russian dandelion (Taraxacum koksaghyz) is a promising source of natural rubber (NR). The synthesis of NR takes place on the surface of organelles known as rubber particles, which are found in latex - the cytoplasm of specialized cells known as laticifers. As well as the enzymes directly responsible for NR synthesis, the rubber particles also contain small rubber particle proteins (SRPPs), the most abundant of which are SRPP3, 4 and 5. These three proteins support NR synthesis by maintaining rubber particle stability. We used homology-based searches to identify the whole TkSRPP gene family and qPCR to create their spatial expression profiles. Affinity enrichment-mass spectrometry was applied to identify TkSRPP3/4/5 protein interaction partners in T. koksaghyz latex and selected interaction partners were analyzed using qPCR, confocal laser scanning microscopy and heterologous expression in yeast. We identified 17 SRPP-like sequences in the T. koksaghyz genome, including three apparent pseudogenes, 10 paralogs arranged as an inverted repeat in a cluster with TkSRPP3/4/5, and one separate gene (TkSRPP6). Their sequence diversity and different expression profiles indicated distinct functions and the latex interactomes obtained for TkSRPP3/4/5 suggested that TkSRPP4 is a promiscuous hub protein that binds many partners from different compartments, whereas TkSRPP3 and 5 have more focused interactomes. Two interactors shared by TkSRPP3/4/5 (TkSRPP6 and TkUGT80B1) were chosen for independent validation and detailed characterization. TkUGT80B1 triterpenoid glycosylating activity provided first evidence for triterpenoid saponin synthesis in T. koksaghyz latex. Based on its identified interaction partners, TkSRPP4 appears to play a special role in the endoplasmic reticulum, interacting with lipidmodifying enzymes that may facilitate rubber particle formation. TkSRPP5 appears to be involved in GTPase-dependent signaling and TkSRPP3 may act as part of a kinase signaling cascade, with roles in stress tolerance. TkSRPP interaction with TkUGT80B1 draws a new connection between TkSRPPs and triterpenoid saponin synthesis in T. koksaghyz latex. Our data contribute to the functional differentiation between TkSRPP paralogs and demonstrate unexpected interactions that will help to further elucidate the network of proteins linking TkSRPPs, stress responses and NR biosynthesis within the cellular complexity of latex.

Genome mining for macrolactam-encoding gene clusters allowed for the network-guided isolation of β-amino acid-containing cyclic derivatives and heterologous production of ciromicin A

β-Amino acid-containing macrolactams represent a structurally diverse group of bioactive natural products derived from polyketides; however we are currently lacking a comprehensive overview about their abundance across bacterial families and the underlying biosynthetic diversity. In this study, we employed a targeted β-amino acid-specific homology-based multi-query search to identify potential bacterial macrolactam producers. Here we demonstrate that approximately 10% of each of the identified actinobacterial genera harbor a biosynthetic gene cluster (BGC) encoding macrolactam production. Based on our comparative study, we propose that mutations occurring in specific regions of polyketide synthases (PKS) are the primary drivers behind the variation in macrolactam ring sizes. We successfully validated two producers of ciromicin A from the genus Amycolatopsis, revised the composition of the biosynthetic gene cluster region mte of macrotermycins, and confirmed the ciromicin biosynthetic pathway through heterologous expression. Additionally, network-based metabolomic analysis uncovered three previously unreported macrotermycin congeners from Amycolatopsis sp. M39. The combination of targeted mining and network-based analysis serves as a powerful tool for identifying macrolactam producers and our studies will catalyze the future discovery of yet unreported macrolactams.

Comparative transcriptome analysis of two contrasting genotypes provides new insights into the drought response mechanism in pigeon pea (Cajanus cajan L. Millsp.).

Despite plant's ability to adapt and withstand challenging environments, drought poses a severe threat to their growth and development. Although pigeon peaisalready quite resistant to drought, the prolonged dehydration induced by the aberrant climate poses a serious threat to their survival and productivity. Comparative physiological and transcriptome analyses of drought-tolerant (CO5) and drought-sensitive (CO1) pigeon pea genotypes subjected to drought stress were carried out in order to understand the molecular basis of drought tolerance in pigeon pea. The transcriptomic analysis allowed us to examine how drought affects the gene expression of C. cajan. Using bioinformatics tools, the unigenes were de novo assembled, annotated, and functionally evaluated. Additionally, a homology-based sequence search against the droughtDB database was performedto identify the orthologs of the DEGs. 1102 potential drought-responsive genes werefound to be differentially expressed genes (DEGs) between drought-tolerant and drought-sensitive genotypes. These included Abscisic acid insensitive 5 (ABI5), Nuclear transcription factor Y subunit A-7 (NF-YA7), WD40 repeat-containing protein 55 (WDR55), Anthocyanidin reductase (ANR) and Zinc-finger homeodomain protein 6 (ZF-HD6) andwere highly expressed in the tolerant genotype. Further, GO analysis revealed that the most enriched classes belonged to biosynthetic and metabolic processes in the biological process category, binding and catalytic activity in the molecular function category and nucleus and protein-containing complex in the cellular component category. Results of KEGG pathway analysis revealed that the DEGs were significantly abundant in signalling pathways such as plant hormone signal transduction and MAPK signalling pathways. Consequently, in our investigation, we have identified and validated by qPCR a group of genes involved in signal reception and propagation, stress-specific TFs, and basal regulatory genes associated with drought response. In conclusion, our comprehensive transcriptome dataset enabled the discovery of candidate genes connected to pathways involved in pigeon pea drought response. Our research uncovered a number of unidentified genes and transcription factors that could be used to understand andimprove susceptibility to drought.

Genome-wide identification and characterization of parthenocarpic fruit set-related gene homologs in cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.), a major horticultural crop, in the family Cucurbitaceae is grown and consumed globally. Parthenocarpy is an ideal trait for many fruit and vegetables which produces seedless fruit desired by consumers. The seedlessness occurs when fruit develops without fertilization which can be either natural or induced. So far, a limited number of genes regulating parthenocarpic fruit set have been reported in several fruit or vegetable crops, most of which are involved in hormone biosynthesis or signalling. Although parthenocarpic cucumber has been widely used in commercial production for a long time; its genetic basis is not well understood. In this study, we retrieved thirty five parthenocarpy fruit-set related genes (PRGs) from bibliomic data in various plants. Thirty-five PRG homologs were identified in the cucumber genome via homology-based search. An in silico analysis was performed on phylogenetic tree, exon–intron structure, cis-regulatory elements in the promoter region, and conserved domains of their deduced proteins, which provided insights into the genetic make-up of parthenocarpy-related genes in cucumber. Simple sequence repeat (SSR) sequences were mined in these PRGs, and 31 SSR markers were designed. SSR genotyping identified three SSRs in two polymorphic genes. Quantitative real-time PCR of selected genes was conducted in five cucumber lines with varying degrees of parthenocarpic fruit set capacities, which revealed possible association of their expression with parthenocarpy. The results revealed that homologs CsWD40 and CsPIN-4 could be considered potential genes for determination of parthenocarpy as these genes showed parental polymorphism and differential gene expression in case of parthenocarpic and non-parthenocarpic parents.

Characterization of Latex-Clearing Protein and Aldehyde Dehydrogenases Involved in the Utilization of poly(cis-1,4-isoprene) by Nocardia farcinica NBRC 15532

Microbial degradation of natural rubber and synthetic poly(cis-1,4-isoprene) is expected to become an alternative treatment system for waste from poly(cis-1,4-isoprene) products including scrap tires. Nocardia farcinica NBRC 15,532, a gram-positive rubber-degrading bacterium, can utilize poly(cis-1,4-isoprene) as the sole source of carbon and energy to produce oligo-isoprene metabolites containing aldehyde and keto end groups. A homology-based search of the genome revealed a gene encoding a latex-clearing protein (Lcp). Gene disruption analysis indicated that this gene is essential for the utilization of poly(cis-1,4-isoprene) in this strain. Further analysis of the genome sequence identified aldehyde dehydrogenase (ALDH) genes as potential candidates for oxidative degradation of oligo-isoprene aldehydes. Based on the enzymatic activity of the ALDH candidates, NF2_RS14000 and NF2_RS14385 may be involved in the degradation of oligo-isoprene aldehydes. Analysis of the reaction products revealed that these ALDHs oxidized tri- to penta-isoprene aldehydes, which were generated by the reaction of Lcp. Based on the inability of ALDH gene deletion mutants, we concluded that NF2_RS14000 is mainly involved in the utilization of poly(cis-1,4-isoprene) and the oxidative degradation of oligo-isoprene aldehydes in Nocardia farcinica NBRC 15,532.

Carotenoid analysis and functional characterization of lycopene cyclases in Zinnia elegans L.

Carotenoids, a group of essential pigments for petal coloration in many species, could be generally categorized into linear and cyclic groups. The cyclization of the linear substrate lycopene at both end groups by β- and ε-cyclases (LCYBs and LCYEs, respectively) is a key branching point of the biosynthesis of cyclic carotenoids in flowers. Common zinnia (Zinnia elegans L.), an Asteraceae family plant that has been widely used in landscape greening, is best known for its brilliant flower colors. However, the mechanism of its petal coloration is not well studied, especially with respect to carotenoid pigmentation. In this study, total carotenoid content and carotenoid compositions of petals were analyzed in two common zinnia cultivars, 'Dreamland Red’ (DRE) and 'Dreamland Yellow’ (DY). Despite the lighter color of DY petals, the total carotenoid content in DY flowers was significantly higher than that in DRE. Metabolomic analysis revealed that the major carotenoids of both cultivars are cyclic carotenoids, with β-carotene being the most abundant. In addition, one LCYB, LCYE, and capsanthin/capsorubin synthase (CCS) gene were identified from Z. elegans by homology-based search in transcriptome. After cloning and sequencing, full-length open reading frame of each gene is exactly the same between two cultivars. Bacterial pigment complementation experiments revealed that ZeLCYB and ZeCCS can cyclize lycopene at both ends to produce β-carotene. Unlike most LCYEs characterized from higher plants, ZeLCYE produced predominantly ε-carotene in E. coli system. However, ZeLCYE is likely to have a lower affinity for lycopene or lower activity compared with ZeLCYB. Gene expression analysis showed that DY has higher expression levels of ZeLCYB and ZeCCS than DRE at specific developmental stages, consistent with its higher total carotenoid content. This work would lay a solid foundation for future studies on carotenoid accumulation in Z. elegans.

Production of retinoic acid by engineered Saccharomyces cerevisiae using an endogenous aldehyde dehydrogenase

Retinoic acid (RA), a vitamin A (retinol)-derived lipophilic compound, is involved in various physiological functions. The demand for RA is growing in the pharmaceutical industry, but RA biosynthesis is still in its infancy compared to other forms of retinoids such as retinol and retinal, largely due to the lack of efficient retinal dehydrogenases. To achieve effective biosynthesis of RA, the catalytic activities of exogenous retinal dehydrogenases were comparatively analyzed in a previously constructed retinoids-producing Saccharomyces cerevisiae strain, followed by mining of endogenous enzymes with higher retinal dehydrogenase activities using homology-based search. After confirming the retinal oxidation activity of the endogenous aldehyde dehydrogenase Hfd1 using in vivo and in vitro experiments, it was overexpressed in multiple copies, and the resulting strain produced 99.71 mg/L of RA in shake-flask cultures. Finally, 545.28 mg/L of RA was produced in fed-batch fermentation. This study suggests the yeast endogenous Hfd1 as a potent catalyst for RA biosynthesis, and demonstrates the potential of yeast as a platform for RA production.

An Integrated Bioinformatics and Functional Approach for miRNA Validation.

MicroRNAs (miRNAs) are small (20-24 nucleotides) non-coding ribo-regulatory molecules with significant roles in regulating target mRNA and long non-coding RNAs at transcriptional and post-transcriptional levels. Rapid advancement in the small RNA sequencing methods with integration of degradome sequencing has accelerated the understanding of miRNA-mediated regulatory hubs in plants and yielded extensive annotation of miRNAs and corresponding targets. However, it is becoming clear that large numbers of such annotations are questionable. Therefore, it is imperative to adopt reliable and strict bioinformatics pipelines for miRNA identification. Furthermore, sensitive methods are needed for validation and functional characterization of miRNA and its target(s). In this chapter, we have provided a comprehensive and streamlined methodology for miRNA identification and its functional validation in plants. This includes a combination of various in silico and experimental methodologies. To identify miRNA compendium from large-scale Next-Generation Sequencing (NGS) small RNA datasets, the miR-PREFeR (miRNA PREdiction From small RNA-Seq data) bioinformatics tool has been described. Also, a homology-based search protocol for finding members of a specific miRNA family has been discussed. The chapter also includes techniques to ascertain miRNA:target pair specificity using in silico target prediction from degradome NGS libraries using CleaveLand pipeline, miRNA:target validation by in planta transient assays, 5' RLM-RACE and expression analysis as well as functional techniques like miRNA overexpression, short tandem target mimic and resistant target approaches. The proposed strategy offers a reliable and sensitive way for miRNA:target identification and validation. Additionally, we strongly promulgate the use of multiple methodologies to validate a miRNA as well as its target.

ReGSP: a visualized application for homology-based gene searching and plotting using multiple reference sequences.

The massively parallel nature of next-generation sequencing technologies has contributed to the generation of massive sequence data in the last two decades. Deciphering the meaning of each generated sequence requires multiple analysis tools, at all stages of analysis, from the reads stage all the way up to the whole-genome level. Homology-based approaches based on related reference sequences are usually the preferred option for gene and transcript prediction in newly sequenced genomes, resulting in the popularity of a variety of BLAST and BLAST-based tools. For organelle genomes, a single-reference–based gene finding tool that uses grouping parameters for BLAST results has been implemented in the Genome Search Plotter (GSP). However, this tool does not accept multiple and user-customized reference sequences required for a broad homology search. Here, we present multiple Reference–based Gene Search and Plot (ReGSP), a simple and convenient web tool that accepts multiple reference sequences for homology-based gene search. The tool incorporates cPlot, a novel dot plot tool, for illustrating nucleotide sequence similarity between the query and the reference sequences. ReGSP has an easy-to-use web interface and is freely accessible at https://ds.mju.ac.kr/regsp.

Identification of putative miRNAs from expressed sequence tags of Gnetum gnemon L. and their cross-kingdom targets.

Wild edible plants are often found to be rich sources of nutrients and medicinally beneficial compounds with pharmacological activities. Gnetum gnemon is a nutritionally important plant and a popular food source in parts of Assam and North-East India. Various microRNAs (miRNAs) have been recently identified in many plants; however, there are no records of identification of miRNAs in any species of Gnetum. The prediction of miRNA-target associations in G. gnemon is an important step to facilitate functional genomics studies in this species. In the present study, all known miRNAs from plants available in public domain were used to search for the conserved G. gnemon miRNA homologues in publicly accessible expressed sequence tags (ESTs) in NCBI database. An aggregate of 20 new potential miRNAs belonging to two diverse miRNA families (miR399 and miR5021) were identified through a homology-based search by following stringent filtering criteria. To investigate the potential cross-kingdom effects of the identified miRNAs, we further identified the putative target genes of G. gnemon miRNAs in human transcriptome and analyzed them against the NCBI non-redundant protein database. The KEGG analysis of the target genes indicated that these genes were involved in different metabolic pathways such as caffeine metabolism, drug metabolism, and nitrotoluene degradation. The target genes of G. gnemon miRNAs in humans were found to be associated with various disorders of both hereditary and non-hereditary origin. These results could help to shed new light on understanding of miRNA-mRNAs functional networks in this species and its potential use as a small RNA-based therapy against some human diseases.

Tracing the active genetic diversity of Microcystis and Microcystis phage through a temporal survey of Taihu.

Harmful algal blooms are commonly thought to be dominated by a single genus, but they are not homogenous communities. Current approaches, both molecular and culture-based, often overlook fine-scale variations in community composition that can influence bloom dynamics. We combined homology-based searches (BLASTX) and phylogenetics to distinguish and quantify Microcystis host and phage members across a summer season during a 2014 Microcystis- dominated bloom that occurred in Lake Tai (Taihu), China. We found 47 different genotypes of the Microcystis-specific DNA-dependent RNA polymerase (rpoB), which included several morphospecies. Microcystis flos-aquae and Microcystis wesenbergii accounted for ~86% of total Microcystis transcripts, while the more commonly studied Microcystis aeruginosa only accounted for ~7%. Microcystis genotypes were classified into three temporal groups according to their expression patterns across the course of the bloom: early, constant and late. All Microcystis morphospecies were present in each group, indicating that expression patterns were likely dictated by competition driven by environmental factors, not phylogeny. We identified three primary Microcystis-infecting phages based on the viral terminase, including a novel Siphoviridae phage that may be capable of lysogeny. Within our dataset, Myoviridae phages consistent with those infecting Microcystis in a lytic manner were positively correlated to the early host genotypes, while the Siphoviridae phages were positively correlated to the late host genotypes, when the Myoviridae phages express putative genetic markers for lysogeny. The expression of genes in the microcystin-encoding mcy cassette was estimated using mcyA, which revealed 24 Microcystis-specific genotypes that were negatively correlated to the early host genotypes. Of all environmental factors measured, pH best described the temporal shift in the Microcystis community genotypic composition, promoting hypotheses regarding carbon concentration mechanisms and oxidative stress. Our work expounds on the complexity of HAB events, using a well-studied dataset to highlight the need for increased resolution of community dynamics.

Identifying Genes Involved in alkaloid Biosynthesis in Vinca minor Through Transcriptomics and Gene Co-Expression Analysis.

The lesser periwinkle Vinca minor accumulates numerous monoterpene indole alkaloids (MIAs) including the vasodilator vincamine. While the biosynthetic pathway of MIAs has been largely elucidated in other Apocynaceae such as Catharanthus roseus, the counterpart in V. minor remains mostly unknown, especially for reactions leading to MIAs specific to this plant. As a consequence, we generated a comprehensive V. minor transcriptome elaborated from eight distinct samples including roots, old and young leaves exposed to low or high light exposure conditions. This optimized resource exhibits an improved completeness compared to already published ones. Through homology-based searches using C. roseus genes as bait, we predicted candidate genes for all common steps of the MIA pathway as illustrated by the cloning of a tabersonine/vincadifformine 16-O-methyltransferase (Vm16OMT) isoform. The functional validation of this enzyme revealed its capacity of methylating 16-hydroxylated derivatives of tabersonine, vincadifformine and lochnericine with a Km 0.94 ± 0.06 µM for 16-hydroxytabersonine. Furthermore, by combining expression of fusions with yellow fluorescent proteins and interaction assays, we established that Vm16OMT is located in the cytosol and forms homodimers. Finally, a gene co-expression network was performed to identify candidate genes of the missing V. minor biosynthetic steps to guide MIA pathway elucidation.

Synthetic Peptide Libraries Designed From a Minimal Alpha-Helical Domain of AS-48-Bacteriocin Homologs Exhibit Potent Antibacterial Activity.

The circularized bacteriocin enterocin AS-48 produced by Enterococcus sp. exhibits antibacterial activity through membrane disruption. The membrane-penetrating activity of enterocin AS-48 has been attributed to a specific alpha-helical region on the circular peptide. Truncated, linearized forms containing these domains have been shown to preserve limited bactericidal activity. We utilized the amino acid sequence of the active helical domain of enterocin AS-48 to perform a homology-based search of similar sequences in other bacterial genomes. We identified similar domains in three previously uncharacterized AS-48-like bacteriocin genes in Clostridium sordellii, Paenibacillus larvae, and Bacillus xiamenensis. Enterocin AS-48 and homologs from these bacterial species were used as scaffolds for the design of a minimal peptide library based on the active helical domain of each bacteriocin sequence. 95 synthetic peptide variants of each scaffold peptide, designated Syn-enterocin, Syn-sordellicin, Syn-larvacin, and Syn-xiamensin, were designed and synthesized from each scaffold sequence based on defined biophysical parameters. A total of 384 total peptides were assessed for antibacterial activity against Gram-negative and Gram-positive bacteria. Minimal Inhibitory Concentrations (MICs) as low as 15.6 nM could be observed for the most potent peptide candidate tested, with no significant cytotoxicity to eukaryotic cells. Our work demonstrates for the first time a general workflow of using minimal domains of natural bacteriocin sequences as scaffolds to design and rapidly synthesize a library of bacteriocin-based antimicrobial peptide variants for evaluation.

The ionotropic receptor gene family in Lepidoptera and Trichoptera: Annotation, evolutionary and functional perspectives

Lepidoptera (moths and butterflies) and Trichoptera (caddisflies), belonging to the superorder Amphiesmenoptera, are the most diverse insect orders as representatives of the terrestrial and aquatic insects, respectively. The insects of the two orders possess different biological and behavioral characteristics, especially their larvae, presumably resulting in the differences of the ionotropic receptor (IR) genes in numbers, sequence characteristics or gene structure. Here, we employed genomics, transcriptomics, bioinformatics, phylogenetics and molecular biology strategies to characterize the IR gene repertoire in Lepidoptera and Trichoptera. Genome and transcriptome analyses with exhaustive homology-based searches and manual efforts, in 32 lepidopterans and five trichopterans, led to the identification of 1449 genes encoding IRs with 1170 full-length sequences, representing the most comprehensive set of chemoreceptor superfamilies across the Amphiesmenoptera. Analysis of gene gains and losses in orthologous groups implied that some IRs were lost in related species, and multiple gene copies occurred mainly in divergent IRs (D-IRs) by gene duplications. Phylogenetic analysis of 2442 IR proteins from 67 species revealed that Lepidoptera and Trichoptera IRs could be classified into three subfamilies, i.e., 14 antennal IRs (A-IRs), five Lepidoptera-specific IRs (LS-IRs) and four D-IRs. Of the three subfamilies, A-IRs and LS-IRs members within orthologous groups exhibited high conservation of gene structure, but D-IRs shared extremely low amino acid identities (below 30%). Expression profiles revealed functional diversities of IRs from Bombyx mori and Papilio xuthus involving smell, taste or reproduction, in which some genes displayed sex-biased expression in antennae associated with specific chemosensory behaviors of female or male adults. Our current study has provided insights into the evolution, conservation and divergence of IRs between/within Lepidoptera and Trichoptera, and allows for further experiments to investigate IR functions.

TERL: classification of transposable elements by convolutional neural networks.

Transposable elements (TEs) are the most represented sequences occurring in eukaryotic genomes. Few methods provide the classification of these sequences into deeper levels, such as superfamily level, which could provide useful and detailed information about these sequences. Most methods that classify TE sequences use handcrafted features such as k-mers and homology-based search, which could be inefficient for classifying non-homologous sequences. Here we propose an approach, called transposable elements pepresentation learner (TERL), that preprocesses and transforms one-dimensional sequences into two-dimensional space data (i.e., image-like data of the sequences) and apply it to deep convolutional neural networks. This classification method tries to learn the best representation of the input data to classify it correctly. We have conducted six experiments to test the performance of TERL against other methods. Our approach obtained macro mean accuracies and F1-score of 96.4% and 85.8% for superfamilies and 95.7% and 91.5% for the order sequences from RepBase, respectively. We have also obtained macro mean accuracies and F1-score of 95.0% and 70.6% for sequences from seven databases into superfamily level and 89.3% and 73.9% for the order level, respectively. We surpassed accuracy, recall and specificity obtained by other methods on the experiment with the classification of order level sequences from seven databases and surpassed by far the time elapsed of any other method for all experiments. Therefore, TERL can learn how to predict any hierarchical level of the TEs classification system and is about 20 times and three orders of magnitude faster than TEclass and PASTEC, respectively https://github.com/muriloHoracio/TERL. Contact:murilocruz@alunos.utfpr.edu.br.

SERK genes identification and expression analysis during somatic embryogenesis and sporogenesis of sexual and apomictic Brachiaria brizantha (Syn. Urochloa brizantha).

In Brachiaria brizantha BbrizSERK1, BbrizSERK2 and BbrizSERK3 were identified. SERK expression marks somatic embryogenesis, sexual MMC, and sexual and apomictic PMC. BbrizSERK3 might have a regulatory role in reproductive development. Somatic embryogenesis receptor-like kinase (SERK) consists of plasma membrane receptor genes that have been characterized in various species, associated with several aspects of plant development, including reproduction. SERK genes are involved in anther development and in early embryo development in sexual and asexual seed formation. To comprehend the complexity of the SERK genes and their function in Brachiaria reproduction, we performed a homology-based search in a genomic database of a sexual B. brizantha and identified sequences of three SERK genes, BbrizSERK1, BbrizSERK2, and BbrizSERK3. RNASeq data showed equivalent abundance of BbrizSERK1 and BbrizSERK2 transcripts in ovaries at early megasporogenesis of sexuals and apomicts, while BbrizSERK3 transcripts were more abundant in ovaries of sexuals than in apomicts. BbrizSERK3 results in three coding sequences due to alternative splicing, among them Variant 1 results in a protein with all the predicted domains of a SERK. BbrizSERK transcripts were detected in male reproductive tissues of both sexual and apomictic plants, suggesting a role in controlling anther development. BbrizSERK transcripts were detected early in ovule development, in the integuments, and in the megaspore mother cell of the sexual plant, butnot in the cells that give rise to apomictic embryo sacs, suggesting a role in female reproductive development of sexuals. This paper provides evidences that SERK genes plays a role in the onset and establishment of somatic embryogenesis and in the reproductive development of B. brizantha and suggests a distinct role of BbrizSERK in apomixis initiation.

The First Draft Genome Assembly of Snow Sheep (Ovis nivicola).

The snow sheep, Ovis nivicola, which is endemic to the mountain ranges of northeastern Siberia, are well adapted to the harsh cold climatic conditions of their habitat. In this study, using long reads of Nanopore sequencing technology, whole-genome sequencing, assembly, and gene annotation of a snow sheep were carried out. Additionally, RNA-seq reads from several tissues were also generated to supplement the gene prediction in snow sheep genome. The assembled genome was ∼2.62 Gb in length and was represented by 7,157 scaffolds with N50 of about 2 Mb. The repetitive sequences comprised of 41% of the total genome. BUSCO analysis revealed that the snow sheep assembly contained full-length or partial fragments of 97% of mammalian universal single-copy orthologs (n = 4,104), illustrating the completeness of the assembly. In addition, a total of 20,045 protein-coding sequences were identified using comprehensive gene prediction pipeline. Of which 19,240 (∼96%) sequences were annotated using protein databases. Moreover, homology-based searches and de novo identification detected 1,484 tRNAs; 243 rRNAs; 1,931 snRNAs; and 782 miRNAs in the snow sheep genome. To conclude, we generated the first de novo genome of the snow sheep using long reads; these data are expected to contribute significantly to our understanding related to evolution and adaptation within the Ovis genus.

MusiteDeep: a deep-learning based webserver for protein post-translational modification site prediction and visualization.

MusiteDeep is an online resource providing a deep-learning framework for protein post-translational modification (PTM) site prediction and visualization. The predictor only uses protein sequences as input and no complex features are needed, which results in a real-time prediction for a large number of proteins. It takes less than three minutes to predict for 1000 sequences per PTM type. The output is presented at the amino acid level for the user-selected PTM types. The framework has been benchmarked and has demonstrated competitive performance in PTM site predictions by other researchers. In this webserver, we updated the previous framework by utilizing more advanced ensemble techniques, and providing prediction and visualization for multiple PTMs simultaneously for users to analyze potential PTM cross-talks directly. Besides prediction, users can interactively review the predicted PTM sites in the context of known PTM annotations and protein 3D structures through homology-based search. In addition, the server maintains a local database providing pre-processed PTM annotations from Uniport/Swiss-Prot for users to download. This database will be updated every three months. The MusiteDeep server is available at https://www.musite.net. The stand-alone tools for locally using MusiteDeep are available at https://github.com/duolinwang/MusiteDeep_web.