Sequenceserver: A Modern Graphical User Interface for Custom BLAST Databases

TeaAS: a comprehensive database for alternative splicing in tea plants (Camellia sinensis)

Alternative splicing (AS) increases the diversity of transcripts and proteins through the selection of different splice sites and plays an important role in the growth, development and stress tolerance of plants. With the release of the reference genome of the tea plant (Camellia sinensis) and the development of transcriptome sequencing, researchers have reported the existence of AS in tea plants. However, there is a lack of a platform, centered on different RNA-seq datasets, that provides comprehensive information on AS.To facilitate access to information on AS and reveal the molecular function of AS in tea plants, we established the first comprehensive AS database for tea plants (TeaAS, http://www.teaas.cn/index.php ). In this study, 3.96 Tb reads from 66 different RNA-seq datasets were collected to identify AS events. TeaAS supports four methods of retrieval of AS information based on gene ID, gene name, annotation (non-redundant/Kyoto encyclopedia of genes and genomes/gene ontology annotation or chromosomal location) and RNA-seq data. It integrates data pertaining to genome annotation, type of AS event, transcript sequence, and isoforms expression levels from 66 RNA-seq datasets. The AS events resulting from different environmental conditions and that occurring in varied tissue types, and the expression levels of specific transcripts can be clearly identified through this online database. Moreover, it also provides two useful tools, Basic Local Alignment Search Tool and Generic Genome Browser, for sequence alignment and visualization of gene structure.The features of the TeaAS database make it a comprehensive AS bioinformatics platform for researchers, as well as a reference for studying AS events in woody crops. It could also be helpful for revealing the novel biological functions of AS in gene regulation in tea plants.

AniProtDB: A Collection of Consistently Generated Metazoan Proteomes for Comparative Genomics Studies

To address the void in the availability of high-quality proteomic data traversing the animal tree, we have implemented a pipeline for generating de novo assemblies based on publicly available data from the NCBI Sequence Read Archive, yielding a comprehensive collection of proteomes from 100 species spanning 21 animal phyla. We have also created the Animal Proteome Database (AniProtDB), a resource providing open access to this collection of high-quality metazoan proteomes, along with information on predicted proteins and protein domains for each taxonomic classification and the ability to perform sequence similarity searches against all proteomes generated using this pipeline. This solution vastly increases the utility of these data by removing the barrier to access for research groups who do not have the expertise or resources to generate these data themselves and enables the use of data from non-traditional research organisms that have the potential to address key questions in biomedicine.

Database Resources for Functional Circular RNAs

Circular RNA (or circRNA) is a type of single-stranded covalently closed circular RNA molecule and play important roles in diverse biological pathways. A comprehensive functionally annotated circRNA database will help to understand the circRNAs and their functions. CircFunBase is such a web-accessible database that aims to provide a high-quality functional circRNA resource including experimentally validated and computationally predicted functions. CircFunBase provides visualized circRNA-miRNA interaction networks. In addition, a genome browser is provided to visualize the genome context of circRNA. In this chapter, we illustrate examples of searching for circRNA and getting detailed information of circRNA. Moreover, other circRNA related databases are outlined.

MolluscDB: a genome and transcriptome database for molluscs

As sequencing becomes more accessible and affordable, the analysis of genomic and transcriptomic data has become a cornerstone of many research initiatives. Communities with a focus on particular taxa or ecosystems need solutions capable of aggregating genomic resources and serving them in a standardized and analysis-friendly manner. Taxon-focussed resources can be more flexible in addressing the needs of a research community than can universal or general databases. Here, we present MolluscDB, a genome and transcriptome database for molluscs. MolluscDB offers a rich ecosystem of tools, including an Ensembl browser, a BLAST server for homology searches and an HTTP server from which any dataset present in the database can be downloaded. To demonstrate the utility of the database and verify the quality of its data, we imported data from assembled genomes and transcriptomes of 22 species, estimated the phylogeny of Mollusca using single-copy orthologues, explored patterns of gene family size change and interrogated the data for biomineralization-associated enzymes and shell matrix proteins. MolluscDB provides an easy-to-use and openly accessible data resource for the research community. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Phylogenomics of Porites from the Arabian Peninsula

The advent of high throughput sequencing technologies provides an opportunity to resolve phylogenetic relationships among closely related species. By incorporating hundreds to thousands of unlinked loci and single nucleotide polymorphisms (SNPs), phylogenomic analyses have a far greater potential to resolve species boundaries than approaches that rely on only a few markers. Scleractinian taxa have proved challenging to identify using traditional morphological approaches and many groups lack an adequate set of molecular markers to investigate their phylogenies. Here, we examine the potential of Restriction-site Associated DNA sequencing (RADseq) to investigate phylogenetic relationships and species limits within the scleractinian coral genus Porites. A total of 126 colonies were collected from 16 localities in the seas surrounding the Arabian Peninsula and ascribed to 12 nominal and two unknown species based on their morphology. Reference mapping was used to retrieve and compare nearly complete mitochondrial genomes, ribosomal DNA, and histone loci. De novo assembly and reference mapping to the P. lobata coral transcriptome were compared and used to obtain thousands of genome-wide loci and SNPs. A suite of species discovery methods (phylogenetic, ordination, and clustering analyses) and species delimitation approaches (coalescent-based, species tree, and Bayesian Factor delimitation) suggested the presence of eight molecular lineages, one of which included six morphospecies. Our phylogenomic approach provided a fully supported phylogeny of Porites from the Arabian Peninsula, suggesting the power of RADseq data to solve the species delineation problem in this speciose coral genus.

Rust expression browser: an open source database for simultaneous analysis of host and pathogen gene expression profiles with expVIP

Background

Transcriptomics is being increasingly applied to generate new insight into the interactions between plants and their pathogens. For the wheat yellow (stripe) rust pathogen (Puccinia striiformis f. sp. tritici, Pst) RNA-based sequencing (RNA-Seq) has proved particularly valuable, overcoming the barriers associated with its obligate biotrophic nature. This includes the application of RNA-Seq approaches to study Pst and wheat gene expression dynamics over time and the Pst population composition through the use of a novel RNA-Seq based surveillance approach called “field pathogenomics”. As a dual RNA-Seq approach, the field pathogenomics technique also provides gene expression data from the host, giving new insight into host responses. However, this has created a wealth of data for interrogation.

Results

Here, we used the field pathogenomics approach to generate 538 new RNA-Seq datasets from Pst-infected field wheat samples, doubling the amount of transcriptomics data available for this important pathosystem. We then analysed these datasets alongside 66 RNA-Seq datasets from four Pst infection time-courses and 420 Pst-infected plant field and laboratory samples that were publicly available. A database of gene expression values for Pst and wheat was generated for each of these 1024 RNA-Seq datasets and incorporated into the development of the rust expression browser (http://www.rust-expression.com). This enables for the first time simultaneous ‘point-and-click’ access to gene expression profiles for Pst and its wheat host and represents the largest database of processed RNA-Seq datasets available for any of the three Puccinia wheat rust pathogens. We also demonstrated the utility of the browser through investigation of expression of putative Pst virulence genes over time and examined the host plants response to Pst infection.

Conclusions

The rust expression browser offers immense value to the wider community, facilitating data sharing and transparency and the underlying database can be continually expanded as more datasets become publicly available.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07488-3.

Mating Type Idiomorphs, Heterothallism, and High Genetic Diversity in Venturia carpophila, Cause of Peach Scab

Scab (caused by Venturia carpophila) is a major disease affecting peach in the eastern United States. The aims of the study were to characterize the mating-type loci in V. carpophila, determine whether they are in equilibrium, and assess the population genetic diversity and structure of the pathogen. The mating-type gene MAT1-1-1 was identified in isolate JP3-5 in an available genome sequence, and the MAT1-2-1 gene was PCR amplified from isolate PS1-1, thus indicating a heterothallic structure. Mating-type loci structures were consistent with those of other Venturia spp. (V. effusa and V. inaequalis): the mating-type gene is positioned between APN2 encoding a DNA lyase and a gene encoding a Pleckstrin homology domain. Primers designed to each of the mating-type genes and a reference gene TUB2 were used as a multiplex PCR to screen a population (n = 81) of V. carpophila from various locations in the eastern United States. Mating types in five of the nine populations studied were in equilibrium. Among the 81 isolates, there were 69 multilocus genotypes. A population genetic analysis of the populations with >10 individuals (four populations) showed them to be genetically diverse. Linkage disequilibrium was found in five of nine populations with ≥4 isolates. A discriminant analysis of principal components indicated three genetic clusters, although extensive admixture was observed. Mating-type identification in V. carpophila provides a basis for understanding reproductive methods of the pathogen and can be a basis for further studies of the genetics of the peach scab pathogen.

Cyanorak v2.1: a scalable information system dedicated to the visualization and expert curation of marine and brackish picocyanobacteria genomes

Cyanorak v2.1 (http://www.sb-roscoff.fr/cyanorak) is an information system dedicated to visualizing, comparing and curating the genomes of Prochlorococcus, Synechococcus and Cyanobium, the most abundant photosynthetic microorganisms on Earth. The database encompasses sequences from 97 genomes, covering most of the wide genetic diversity known so far within these groups, and which were split into 25,834 clusters of likely orthologous groups (CLOGs). The user interface gives access to genomic characteristics, accession numbers as well as an interactive map showing strain isolation sites. The main entry to the database is through search for a term (gene name, product, etc.), resulting in a list of CLOGs and individual genes. Each CLOG benefits from a rich functional annotation including EggNOG, EC/K numbers, GO terms, TIGR Roles, custom-designed Cyanorak Roles as well as several protein motif predictions. Cyanorak also displays a phyletic profile, indicating the genotype and pigment type for each CLOG, and a genome viewer (Jbrowse) to visualize additional data on each genome such as predicted operons, genomic islands or transcriptomic data, when available. This information system also includes a BLAST search tool, comparative genomic context as well as various data export options. Altogether, Cyanorak v2.1 constitutes an invaluable, scalable tool for comparative genomics of ecologically relevant marine microorganisms.

A novel terpene synthase controls differences in anti-aphrodisiac pheromone production between closely related Heliconius butterflies

Plants and insects often use the same compounds for chemical communication, but not much is known about the genetics of convergent evolution of chemical signals. The terpene (E)-β-ocimene is a common component of floral scent and is also used by the butterfly Heliconius melpomene as an anti-aphrodisiac pheromone. While the biosynthesis of terpenes has been described in plants and microorganisms, few terpene synthases (TPSs) have been identified in insects. Here, we study the recent divergence of 2 species, H. melpomene and Heliconius cydno, which differ in the presence of (E)-β-ocimene; combining linkage mapping, gene expression, and functional analyses, we identify 2 novel TPSs. Furthermore, we demonstrate that one, HmelOS, is able to synthesise (E)-β-ocimene in vitro. We find no evidence for TPS activity in HcydOS (HmelOS ortholog of H. cydno), suggesting that the loss of (E)-β-ocimene in this species is the result of coding, not regulatory, differences. The TPS enzymes we discovered are unrelated to previously described plant and insect TPSs, demonstrating that chemical convergence has independent evolutionary origins.

Plants and insects often use the same compounds for chemical communication, but little is known about the convergent evolution of such chemical signals. This study identifies a novel terpene synthase involved in production of an anti-aphrodisiac pheromone by the butterfly Heliconius melpomene. This enzyme is unrelated to other insect terpene synthases, providing evidence that the ability to synthesise terpenes has arisen multiple times independently within the insects.

PSDX: A Comprehensive Multi-Omics Association Database of Populus trichocarpa With a Focus on the Secondary Growth in Response to Stresses

Populus trichocarpa (P. trichocarpa) is a model tree for the investigation of wood formation. In recent years, researchers have generated a large number of high-throughput sequencing data in P. trichocarpa. However, no comprehensive database that provides multi-omics associations for the investigation of secondary growth in response to diverse stresses has been reported. Therefore, we developed a public repository that presents comprehensive measurements of gene expression and post-transcriptional regulation by integrating 144 RNA-Seq, 33 ChIP-seq, and six single-molecule real-time (SMRT) isoform sequencing (Iso-seq) libraries prepared from tissues subjected to different stresses. All the samples from different studies were analyzed to obtain gene expression, co-expression network, and differentially expressed genes (DEG) using unified parameters, which allowed comparison of results from different studies and treatments. In addition to gene expression, we also identified and deposited pre-processed data about alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI). The post-transcriptional regulation, differential expression, and co-expression network datasets were integrated into a new P. trichocarpa Stem Differentiating Xylem (PSDX) database (http://forestry.fafu.edu.cn/db/SDX), which further highlights gene families of RNA-binding proteins and stress-related genes. The PSDX also provides tools for data query, visualization, a genome browser, and the BLAST option for sequence-based query. Much of the data is also available for bulk download. The availability of PSDX contributes to the research related to the secondary growth in response to stresses in P. trichocarpa, which will provide new insights that can be useful for the improvement of stress tolerance in woody plants.

Comparative Genome Analysis Reveals Differences in Biocontrol Efficacy According to Each Individual Isolate Belonging to Rhizospheric Fluorescent Pseudomonads

Biocontrol fluorescent pseudomonads produce a number of antibiotic organic compounds, including 2,4-diacetylphloroglucinol, pyoluteorin, pyrrolnitrin, and phenazine. We previously classified rhizospheric fluorescent pseudomonads harboring antibiotic biosynthetic gene clusters into 10 operational taxonomic units (OTUs). In the present study, we report the complete genome sequences of selected strains from these OTUs. The genetic diversity of antibiotic biosynthetic gene clusters and their surrounding sequences correlated with the OTU classification. In comparisons of the biocontrol activity and distribution of antibiotic biosynthetic gene clusters, we found that the pyrrolnitrin biosynthetic gene cluster more effectively controlled the growth of Rhizoctonia solani.

NBIGV-DB: A dedicated database of non-B cell derived immunoglobulin variable region

Immunoglobulins (Ig) are important immune molecules that possess highly diverse variable region sequences enabling antigen recognition. According to classical immune theory, B lymphocytes have been considered the only source of Ig production (B-Igs). However, accumulating evidence have suggested that Igs are also produced by many non-B cells (non-B Igs), including epithelial cells, neurons, germ cells, as well as myeloid cells of hemopoietic system. Besides acting as bona fide antibodies, Non-B Igs have alternative cellular functions, such as promotion of cell survival, adhesion and migration. More importantly, Unlike the unlimited sequence diversity of B-Igs, the non-B Igs exhibit conserved V(D)J patterns across the same lineages. To support the analysis and comparison of variable region sequences from Igs, produced by B and non-B cells, we established a database (NBIGV) constituted by a non-B Ig variable region repertoire, which includes 727,989 V_HDJ_H and VκJκ recombination sequences of non-B Igs sequenced from mouse samples. Upon database search, users can view, browse and investigate the variable region sequences of non-B Igs according to respective mice strains and tissues as well as Ig classes. Moreover, users can easily download selected sequences and/or compare sequences of interest with known non-B Ig sequences present in the database using NCBI-BLAST algorithms. Additionally, our database integrates a submission page and supplementary sample information. The NBIGV database may serve as a valuable resource for sequence analyses of Non-B Igs. NBIGV database is freely available at http://nbigv.org.

Using online tools at the Bovine Genome Database to manually annotate genes in the new reference genome

With the availability of a new highly contiguous Bos taurus reference genome assembly (ARS-UCD1.2), it is the opportune time to upgrade the bovine gene set by seeking input from researchers. Furthermore, advances in graphical genome annotation tools now make it possible for researchers to leverage sequence data generated with the latest technologies to collaboratively curate genes. For many years the Bovine Genome Database (BGD) has provided tools such as the Apollo genome annotation editor to support manual bovine gene curation. The goal of this paper is to explain the reasoning behind the decisions made in the manual gene curation process while providing examples using the existing BGD tools. We will describe the sources of gene annotation evidence provided at the BGD, including RNA-seq and Iso-Seq data. We will also explain how to interpret various data visualizations when curating gene models, and will demonstrate the value of manual gene annotation. The process described here can be applied to manual gene curation for other species with similar tools. With a better understanding of manual gene annotation, researchers will be encouraged to edit gene models and contribute to the enhancement of livestock gene sets.

Mechanistic concepts of iron-sulfur protein biogenesis in Biology

Iron-sulfur (Fe/S) proteins are present in virtually all living organisms and are involved in numerous cellular processes such as respiration, photosynthesis, metabolic reactions, nitrogen fixation, radical biochemistry, protein synthesis, antiviral defense, and genome maintenance. Their versatile functions may go back to the proposed role of their Fe/S cofactors in the origin of life as efficient catalysts and electron carriers. More than two decades ago, it was discovered that the in vivo synthesis of cellular Fe/S clusters and their integration into polypeptide chains requires assistance by complex proteinaceous machineries, despite the fact that Fe/S proteins can be assembled chemically in vitro. In prokaryotes, three Fe/S protein biogenesis systems are known; ISC, SUF, and the more specialized NIF. The former two systems have been transferred by endosymbiosis from bacteria to mitochondria and plastids, respectively, of eukaryotes. In their cytosol, eukaryotes use the CIA machinery for the biogenesis of cytosolic and nuclear Fe/S proteins. Despite the structural diversity of the protein constituents of these four machineries, general mechanistic concepts underlie the complex process of Fe/S protein biogenesis. This review provides a comprehensive and comparative overview of the various known biogenesis systems in Biology, and summarizes their common or diverging molecular mechanisms, thereby illustrating both the conservation and diverse adaptions of these four machineries during evolution and under different lifestyles. Knowledge of these fundamental biochemical pathways is not only of basic scientific interest, but is important for the understanding of human 'Fe/S diseases' and can be used in biotechnology.

A gene prioritization method based on a swine multi-omics knowledgebase and a deep learning model

The analyses of multi-omics data have revealed candidate genes for objective traits. However, they are integrated poorly, especially in non-model organisms, and they pose a great challenge for prioritizing candidate genes for follow-up experimental verification. Here, we present a general convolutional neural network model that integrates multi-omics information to prioritize the candidate genes of objective traits. By applying this model to Sus scrofa, which is a non-model organism, but one of the most important livestock animals, the model precision was 72.9%, recall 73.5%, and F1-Measure 73.4%, demonstrating a good prediction performance compared with previous studies in Arabidopsis thaliana and Oryza sativa. Additionally, to facilitate the use of the model, we present ISwine (http://iswine.iomics.pro/), which is an online comprehensive knowledgebase in which we incorporated almost all the published swine multi-omics data. Overall, the results suggest that the deep learning strategy will greatly facilitate analyses of multi-omics integration in the future.

Yuhua Fu et al. develop a CNN model that integrates multi-omics information to prioritize candidate genes of objective traits. Their model performs well when applied to important livestock non-model animals like Sus scrofa. Finally, the authors present ISwine, an online comprehensive knowledgebase which includes all published swine omics data to facilitate the integration of heterogeneous data.

Molecular and functional characterization of somatostatin-type signalling in a deuterostome invertebrate

Somatostatin (SS) and allatostatin-C (ASTC) are structurally and evolutionarily related neuropeptides that act as inhibitory regulators of physiological processes in mammals and insects, respectively. Here, we report the first molecular and functional characterization of SS/ASTC-type signalling in a deuterostome invertebrate—the starfish Asterias rubens (phylum Echinodermata). Two SS/ASTC-type precursors were identified in A. rubens (ArSSP1 and ArSSP2) and the structures of neuropeptides derived from these proteins (ArSS1 and ArSS2) were analysed using mass spectrometry. Pharmacological characterization of three cloned A. rubens SS/ASTC-type receptors (ArSSR1–3) revealed that ArSS2, but not ArSS1, acts as a ligand for all three receptors. Analysis of ArSS2 expression in A. rubens using mRNA in situ hybridization and immunohistochemistry revealed stained cells/fibres in the central nervous system, the digestive system (e.g. cardiac stomach) and the body wall and its appendages (e.g. tube feet). Furthermore, in vitro pharmacological tests revealed that ArSS2 causes dose-dependent relaxation of tube foot and cardiac stomach preparations, while injection of ArSS2 in vivo causes partial eversion of the cardiac stomach. Our findings provide new insights into the molecular evolution of SS/ASTC-type signalling in the animal kingdom and reveal an ancient role of SS-type neuropeptides as inhibitory regulators of muscle contractility.

A Natural Mutational Event Uncovers a Life History Trade-Off via Hormonal Pleiotropy

Environmental signals often control central life history decisions, including the choice between reproduction and somatic maintenance. Such adaptive developmental plasticity occurs in the nematode Caenorhabditis elegans, where environmental cues govern whether larvae will develop directly into reproducing adults or arrest their development to become stress-resistant dauer larvae. Here, we identified a natural variant underlying enhanced sensitivity to dauer-inducing cues in C. elegans: a 92-bp deletion in the cis-regulatory region of the gene eak-3. This deletion reduces synthesis or activity of the steroid hormone dafachronic acid (DA), thereby increasing environmental sensitivity for dauer induction. Consistent with known pleiotropic roles of DA, this eak-3 variant significantly slows down reproductive growth. We experimentally show that, although the eak-3 deletion can provide a fitness advantage through facilitated dauer production in stressful environments, this allele becomes rapidly outcompeted in favorable environments. The identified eak-3 variant therefore reveals a trade-off in how hormonal responses influence both the pace of developmental timing and the way in which environmental sensitivity controls adaptive plasticity. Together, our results show how a single mutational event altering hormonal signaling can lead to the emergence of a complex life history trade-off.

The immunogenetics of sexual parasitism

Sexual parasitism has evolved as a distinctive mode of reproduction among deep-sea anglerfishes. The permanent attachment of males to host females observed in these species represents a form of anatomical joining, which is otherwise unknown in nature. Pronounced modifications to immune facilities are associated with this reproductive trait. The genomes of species with temporarily attaching males lack functional aicda genes that underpin affinity maturation of antibodies. Permanent attachment is associated with additional alterations, culminating in the loss of functional rag genes in some species, abolishing somatic diversification of antigen receptor genes, the hallmark of canonical adaptive immunity. In anglerfishes, coevolution of innate and adaptive immunity has been disentangled, implying that an alternative form of immunity supported the emergence of this evolutionarily successful group of vertebrates.

Laboratory culture of the California Sea Firefly Vargula tsujii (Ostracoda: Cypridinidae): Developing a model system for the evolution of marine bioluminescence

Bioluminescence, or the production of light by living organisms via chemical reaction, is widespread across Metazoa. Laboratory culture of bioluminescent organisms from diverse taxonomic groups is important for determining the biosynthetic pathways of bioluminescent substrates, which may lead to new tools for biotechnology and biomedicine. Some bioluminescent groups may be cultured, including some cnidarians, ctenophores, and brittle stars, but those use luminescent substrates (luciferins) obtained from their diets, and therefore are not informative for determination of the biosynthetic pathways of the luciferins. Other groups, including terrestrial fireflies, do synthesize their own luciferin, but culturing them is difficult and the biosynthetic pathway for firefly luciferin remains unclear. An additional independent origin of endogenous bioluminescence is found within ostracods from the family Cypridinidae, which use their luminescence for defense and, in Caribbean species, for courtship displays. Here, we report the first complete life cycle of a luminous ostracod (Vargula tsujii Kornicker & Baker, 1977, the California Sea Firefly) in the laboratory. We also describe the late-stage embryogenesis of Vargula tsujii and discuss the size classes of instar development. We find embryogenesis in V. tsujii ranges from 25–38 days, and this species appears to have five instar stages, consistent with ontogeny in other cypridinid lineages. We estimate a complete life cycle at 3–4 months. We also present the first complete mitochondrial genome for Vargula tsujii. Bringing a luminous ostracod into laboratory culture sets the stage for many potential avenues of study, including learning the biosynthetic pathway of cypridinid luciferin and genomic manipulation of an autogenic bioluminescent system.

Dynamic transcriptional and chromatin accessibility landscape of medaka embryogenesis

Medaka (Oryzias latipes) has become an important vertebrate model widely used in genetics, developmental biology, environmental sciences, and many other fields. A high-quality genome sequence and a variety of genetic tools are available for this model organism. However, existing genome annotation is still rudimentary, as it was mainly based on computational prediction and short-read RNA-seq data. Here we report a dynamic transcriptome landscape of medaka embryogenesis profiled by long-read RNA-seq, short-read RNA-seq, and ATAC-seq. By integrating these data sets, we constructed a much-improved gene model set including about 17,000 novel isoforms and identified 1600 transcription factors, 1100 long noncoding RNAs, and 150,000 potential cis-regulatory elements as well. Time-series data sets provided another dimension of information. With the expression dynamics of genes and accessibility dynamics of cis-regulatory elements, we investigated isoform switching, as well as regulatory logic between accessible elements and genes, during embryogenesis. We built a user-friendly medaka omics data portal to present these data sets. This resource provides the first comprehensive omics data sets of medaka embryogenesis. Ultimately, we term these three assays as the minimum ENCODE toolbox and propose the use of it as the initial and essential profiling genomic assays for model organisms that have limited data available. This work will be of great value for the research community using medaka as the model organism and many others as well.

Echinoderms provide missing link in the evolution of PrRP/sNPF-type neuropeptide signalling

Neuropeptide signalling systems comprising peptide ligands and cognate receptors are evolutionarily ancient regulators of physiology and behaviour. However, there are challenges associated with determination of orthology between neuropeptides in different taxa. Orthologs of vertebrate neuropeptide-Y (NPY) known as neuropeptide-F (NPF) have been identified in protostome invertebrates, whilst prolactin-releasing peptide (PrRP) and short neuropeptide-F (sNPF) have been identified as paralogs of NPY/NPF in vertebrates and protostomes, respectively. Here we investigated the occurrence of NPY/NPF/PrRP/sNPF-related signalling systems in a deuterostome invertebrate phylum - the Echinodermata. Analysis of transcriptome/genome sequence data revealed loss of NPY/NPF-type signalling, but orthologs of PrRP-type neuropeptides and sNPF/PrRP-type receptors were identified in echinoderms. Furthermore, experimental studies revealed that the PrRP-type neuropeptide pQDRSKAMQAERTGQLRRLNPRF-NH₂ is a potent ligand for a sNPF/PrRP-type receptor in the starfish Asterias rubens. Our findings indicate that PrRP-type and sNPF-type signalling systems are orthologous and originated as a paralog of NPY/NPF-type signalling in Urbilateria.

ZEAMAP, a Comprehensive Database Adapted to the Maize Multi-Omics Era

As one of the most extensively cultivated crops, maize (Zea mays L.) has been extensively studied by researchers and breeders for over a century. With advances in high-throughput detection of various omics data, a wealth of multi-dimensional and multi-omics information has been accumulated for maize and its wild relative, teosinte. Integration of this information has the potential to accelerate genetic research and generate improvements in maize agronomic traits. To this end, we constructed ZEAMAP, a comprehensive database incorporating multiple reference genomes, annotations, comparative genomics, transcriptomes, open chromatin regions, chromatin interactions, high-quality genetic variants, phenotypes, metabolomics, genetic maps, genetic mapping loci, population structures, and populational DNA methylation signals within maize inbred lines. ZEAMAP is user friendly, with the ability to interactively integrate, visualize, and cross-reference multiple different omics datasets.

•

Functional annotations and comparative genomics of maize and teosinte genomes

•

Multi-omics data generated from the same maize inbred lines panel

•

Interactive tools to query, cross-refer, and visualize the omics data

Bovine Genome Database: new annotation tools for a new reference genome

The Bovine Genome Database (BGD) (http://bovinegenome.org) has been the key community bovine genomics database for more than a decade. To accommodate the increasing amount and complexity of bovine genomics data, BGD continues to advance its practices in data acquisition, curation, integration and efficient data retrieval. BGD provides tools for genome browsing (JBrowse), genome annotation (Apollo), data mining (BovineMine) and sequence database searching (BLAST). To augment the BGD genome annotation capabilities, we have developed a new Apollo plug-in, called the Locus-Specific Alternate Assembly (LSAA) tool, which enables users to identify and report potential genome assembly errors and structural variants. BGD now hosts both the newest bovine reference genome assembly, ARS-UCD1.2, as well as the previous reference genome, UMD3.1.1, with cross-genome navigation and queries supported in JBrowse and BovineMine, respectively. Other notable enhancements to BovineMine include the incorporation of genomes and gene annotation datasets for non-bovine ruminant species (goat and sheep), support for multiple assemblies per organism in the Regions Search tool, integration of additional ontologies and development of many new template queries. To better serve the research community, we continue to focus on improving existing tools, developing new tools, adding new datasets and encouraging researchers to use these resources.

MaGenDB: a functional genomics hub for Malvaceae plants

Malvaceae is a family of flowering plants containing many economically important plant species including cotton, cacao and durian. Recently, the genomes of several Malvaceae species have been decoded, and many omics data were generated for individual species. However, no integrative database of multiple species, enabling users to jointly compare and analyse relevant data, is available for Malvaceae. Thus, we developed a user-friendly database named MaGenDB (http://magen.whu.edu.cn) as a functional genomics hub for the plant community. We collected the genomes of 13 Malvaceae species, and comprehensively annotated genes from different perspectives including functional RNA/protein element, gene ontology, KEGG orthology, and gene family. We processed 374 sets of diverse omics data with the ENCODE pipelines and integrated them into a customised genome browser, and designed multiple dynamic charts to present gene/RNA/protein-level knowledge such as dynamic expression profiles and functional elements. We also implemented a smart search system for efficiently mining genes. In addition, we constructed a functional comparison system to help comparative analysis between genes on multiple features in one species or across closely related species. This database and associated tools will allow users to quickly retrieve large-scale functional information for biological discovery.

Coexpression Analysis Reveals Dynamic Modules Regulating the Growth and Development of Cirri in the Rattans (Calamus simplicifolius and Daemonorops jenkinsiana)

Rattan is regarded as one of the major non-timber forest products, second only to wood and bamboo, worldwide. Although the published genomes of Calamus simplicifolius and Daemonorops jenkinsiana have facilitated genome-wide gene functional analyses, coexpression networks (CENs) provide more comprehensive and complete annotations of gene function at the transcriptome level. Thus, we analyzed the CENs of the two rattans, C. simplicifolius and D. jenkinsiana, by integrating the genome sequences and analyzing in-house transcriptome data from different development stages of their cirri using a well-developed strategy. A total of 3,504 and 3,027 functional modules were identified in C. simplicifolius and D. jenkinsiana, respectively, based on a combination of CENs, gene family classification, and function enrichment tools. These modules covered the major developmental processes, including photosynthesis, lignin biosynthesis, flavonoid biosynthesis, and phenylpropanoid biosynthesis. Reference annotations were refined using CENs and functional modules. Moreover, we obtained novel insights into the regulation of cirrus growth and development in rattans. Furthermore, Rattan-NET (http://rattan.bamboogdb.org/), an online database with analysis tools for gene set enrichment analysis, module enrichment, network comparison analysis, and cis-element analysis, was constructed for the easy analysis of gene function and regulation modules involved in the growth and development of cirri in rattans.

The Mnemiopsis Genome Project Portal: integrating new gene expression resources and improving data visualization

Following the completion of the genome sequencing and gene prediction of Mnemiopsis leidyi, a lobate ctenophore that is native to the coastal waters of the western Atlantic Ocean, we developed and implemented the Mnemiopsis Genome Project Portal (MGP Portal), a comprehensive Web-based data portal for navigating the genome sequence and gene annotations. In the years following the first release of the MGP Portal, it has become evident that the inclusion of data from significant published studies on Mnemiopsis has been critical to its adoption as the centralized resource for this emerging model organism. With this most recent update, the Portal has significantly expanded to include in situ images, temporal developmental expression profiles and single-cell expression data. Recent enhancements also include implementations of an updated BLAST interface, new graphical visualization tools and updates to gene pages that integrate all new data types. Database URL: https://research.nhgri.nih.gov/mnemiopsis/.

The chloroalkaloid (-)-acutumine is biosynthesized via a Fe(II)- and 2-oxoglutarate-dependent halogenase in Menispermaceae plants

Plant halogenated natural products are rare and harbor various interesting bioactivities, yet the biochemical basis for the involved halogenation chemistry is unknown. While a handful of Fe(II)- and 2-oxoglutarate-dependent halogenases (2ODHs) have been found to catalyze regioselective halogenation of unactivated C–H bonds in bacteria, they remain uncharacterized in the plant kingdom. Here, we report the discovery of dechloroacutumine halogenase (DAH) from Menispermaceae plants known to produce the tetracyclic chloroalkaloid (−)-acutumine. DAH is a 2ODH of plant origin and catalyzes the terminal chlorination step in the biosynthesis of (−)-acutumine. Phylogenetic analyses reveal that DAH evolved independently in Menispermaceae plants and in bacteria, illustrating an exemplary case of parallel evolution in specialized metabolism across domains of life. We show that at the presence of azide anion, DAH also exhibits promiscuous azidation activity against dechloroacutumine. This study opens avenues for expanding plant chemodiversity through halogenation and azidation biochemistry.

Halogenated plant natural products are rare and plant halogenation enzymes are thus far unknown. Here Kim et al. identify a dechloroacutumine halogenase from Common Moonseed that catalyzes the final chlorination step in the biosynthesis of acutumine, a chloroalkaloid with selective cytotoxicity to cultured T cells.

Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced-complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over-represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini - the largest radiation of mimetic butterflies.

Genomic, metabolic and phenotypic variability shapes ecological differentiation and intraspecies interactions of Alteromonas macleodii

Ecological differentiation between strains of bacterial species is shaped by genomic and metabolic variability. However, connecting genotypes to ecological niches remains a major challenge. Here, we linked bacterial geno- and phenotypes by contextualizing pangenomic, exometabolomic and physiological evidence in twelve strains of the marine bacterium Alteromonas macleodii, illuminating adaptive strategies of carbon metabolism, microbial interactions, cellular communication and iron acquisition. In A. macleodii strain MIT1002, secretion of amino acids and the unique capacity for phenol degradation may promote associations with Prochlorococcus cyanobacteria. Strain 83-1 and three novel Pacific isolates, featuring clonal genomes despite originating from distant locations, have profound abilities for algal polysaccharide utilization but without detrimental implications for Ecklonia macroalgae. Degradation of toluene and xylene, mediated via a plasmid syntenic to terrestrial Pseudomonas, was unique to strain EZ55. Benzoate degradation by strain EC673 related to a chromosomal gene cluster shared with the plasmid of A. mediterranea EC615, underlining that mobile genetic elements drive adaptations. Furthermore, we revealed strain-specific production of siderophores and homoserine lactones, with implications for nutrient acquisition and cellular communication. Phenotypic variability corresponded to different competitiveness in co-culture and geographic distribution, indicating linkages between intraspecific diversity, microbial interactions and biogeography. The finding of "ecological microdiversity" helps understanding the widespread occurrence of A. macleodii and contributes to the interpretation of bacterial niche specialization, population ecology and biogeochemical roles.

SAGER: a database of Symbiodiniaceae and Algal Genomic Resource

Symbiodiniaceae dinoflagellates are essential endosymbionts of reef building corals and some other invertebrates. Information of their genome structure and function is critical for understanding coral symbiosis and bleaching. With the rapid development of sequencing technology, genome draft assemblies of several Symbiodiniaceae species and diverse marine algal genomes have become publicly available but spread in multiple separate locations. Here, we present a Symbiodiniaceae and Algal Genomic Resource Database (SAGER), a user-friendly online repository for integrating existing genomic data of Symbiodiniaceae species and diverse marine algal gene sets from MMETSP and PhyloDB databases. Relevant algal data are included to facilitate comparative analyses. The database is freely accessible at http://sampgr.org.cn. It provides comprehensive tools for studying gene function, expression and comparative genomics, including search tools to identify gene information from Symbiodiniaceae species, and BLAST tool to find orthologs from marine algae and protists. Moreover, SAGER integrates transcriptome datasets derived from diverse culture conditions of corresponding Symbiodiniaceae species. SAGER was developed with the capacity to incorporate future Symbiodiniaceae and algal genome and transcriptome data, and will serve as an open-access and sustained platform providing genomic and molecular tools that can be conveniently used to study Symbiodiniaceae and other marine algae. Database URL: http://sampgr.org.cn.

PncStress: a manually curated database of experimentally validated stress-responsive non-coding RNAs in plants

Non-coding RNAs (ncRNAs) are recognized as key regulatory molecules in many biological processes. Accumulating evidence indicates that ncRNA-related mechanisms play important roles in plant stress responses. Although abundant plant stress-responsive ncRNAs have been identified, these experimentally validated results have not been gathered into a single public domain archive. Therefore, we established PncStress by curating experimentally validated stress-responsive ncRNAs in plants, including microRNAs, long non-coding RNAs and circular RNAs. The current version of PncStress contains 4227 entries from 114 plants covering 48 biotic and 91 abiotic stresses. For each entry, PncStress has biological information and network visualization. Serving as a manually curated database, PncStress will become a valuable resource in support of plant stress response research.

RNA-seq, de novo transcriptome assembly and flavonoid gene analysis in 13 wild and cultivated berry fruit species with high content of phenolics

BACKGROUND

Flavonoids are produced in all flowering plants in a wide range of tissues including in berry fruits. These compounds are of considerable interest for their biological activities, health benefits and potential pharmacological applications. However, transcriptomic and genomic resources for wild and cultivated berry fruit species are often limited, despite their value in underpinning the in-depth study of metabolic pathways, fruit ripening as well as in the identification of genotypes rich in bioactive compounds.

RESULTS

To access the genetic diversity of wild and cultivated berry fruit species that accumulate high levels of phenolic compounds in their fleshy berry(-like) fruits, we selected 13 species from Europe, South America and Asia representing eight genera, seven families and seven orders within three clades of the kingdom Plantae. RNA from either ripe fruits (ten species) or three ripening stages (two species) as well as leaf RNA (one species) were used to construct, assemble and analyse de novo transcriptomes. The transcriptome sequences are deposited in the BacHBerryGEN database (http://jicbio.nbi.ac.uk/berries) and were used, as a proof of concept, via its BLAST portal (http://jicbio.nbi.ac.uk/berries/blast.html) to identify candidate genes involved in the biosynthesis of phenylpropanoid compounds. Genes encoding regulatory proteins of the anthocyanin biosynthetic pathway (MYB and basic helix-loop-helix (bHLH) transcription factors and WD40 repeat proteins) were isolated using the transcriptomic resources of wild blackberry (Rubus genevieri) and cultivated red raspberry (Rubus idaeus cv. Prestige) and were shown to activate anthocyanin synthesis in Nicotiana benthamiana. Expression patterns of candidate flavonoid gene transcripts were also studied across three fruit developmental stages via the BacHBerryEXP gene expression browser (http://www.bachberryexp.com) in R. genevieri and R. idaeus cv. Prestige.

CONCLUSIONS

We report a transcriptome resource that includes data for a wide range of berry(-like) fruit species that has been developed for gene identification and functional analysis to assist in berry fruit improvement. These resources will enable investigations of metabolic processes in berries beyond the phenylpropanoid biosynthetic pathway analysed in this study. The RNA-seq data will be useful for studies of berry fruit development and to select wild plant species useful for plant breeding purposes.

High-density linkage map construction and QTL analyses for fiber quality, yield and morphological traits using CottonSNP63K array in upland cotton (Gossypium hirsutum L.)

Background

Improving fiber quality and yield are the primary research objectives in cotton breeding for enhancing the economic viability and sustainability of Upland cotton production. Identifying the quantitative trait loci (QTL) for fiber quality and yield traits using the high-density SNP-based genetic maps allows for bridging genomics with cotton breeding through marker assisted and genomic selection. In this study, a recombinant inbred line (RIL) population, derived from cross between two parental accessions, which represent broad allele diversity in Upland cotton, was used to construct high-density SNP-based linkage maps and to map the QTLs controlling important cotton traits.

Results

Molecular genetic mapping using RIL population produced a genetic map of 3129 SNPs, mapped at a density of 1.41 cM. Genetic maps of the individual chromosomes showed good collinearity with the sequence based physical map. A total of 106 QTLs were identified which included 59 QTLs for six fiber quality traits, 38 QTLs for four yield traits and 9 QTLs for two morphological traits. Sub-genome wide, 57 QTLs were mapped in A sub-genome and 49 were mapped in D sub-genome. More than 75% of the QTLs with favorable alleles were contributed by the parental accession NC05AZ06. Forty-six mapped QTLs each explained more than 10% of the phenotypic variation. Further, we identified 21 QTL clusters where 12 QTL clusters were mapped in the A sub-genome and 9 were mapped in the D sub-genome. Candidate gene analyses of the 11 stable QTL harboring genomic regions identified 19 putative genes which had functional role in cotton fiber development.

Conclusion

We constructed a high-density genetic map of SNPs in Upland cotton. Collinearity between genetic and physical maps indicated no major structural changes in the genetic mapping populations. Most traits showed high broad-sense heritability. One hundred and six QTLs were identified for the fiber quality, yield and morphological traits. Majority of the QTLs with favorable alleles were contributed by improved parental accession. More than 70% of the mapped QTLs shared the similar map position with previously reported QTLs which suggest the genetic relatedness of Upland cotton germplasm. Identification of QTL clusters could explain the correlation among some fiber quality traits in cotton. Stable and major QTLs and QTL clusters of traits identified in the current study could be the targets for map-based cloning and marker assisted selection (MAS) in cotton breeding. The genomic region on D12 containing the major stable QTLs for micronaire, fiber strength and lint percentage could be potential targets for MAS and gene cloning of fiber quality traits in cotton.

First mitochondrial genome of a periwinkle from the genus Littoraria: Littoraria sinensis

The nearly complete sequence of the mitochondrial genome of the periwinkle Littoraria sinensis is presented. This is the first mitochondrial genome of the genus Littoraria and the fifth sequenced Littorinidae. The mitogenome of L. sinensis is 16,320 bp in size, with 37 genes and one partial control region. The genes order is identical to those observed in other Littorinidae mitogenomes. The overall base composition shows that AT content (64.76%) is higher than GC content (35.24%). This study should provide new information for the further phylogenetic studies of this species.

Independent evolution of rosmarinic acid biosynthesis in two sister families under the Lamiids clade of flowering plants

As a means to maintain their sessile lifestyle amid challenging environments, plants produce an enormous diversity of compounds as chemical defenses against biotic and abiotic insults. The underpinning metabolic pathways that support the biosynthesis of these specialized chemicals in divergent plant species provide a rich arena for understanding the molecular evolution of complex metabolic traits. Rosmarinic acid (RA) is a phenolic natural product first discovered in plants of the mint family (Lamiaceae) and is recognized for its wide range of medicinal properties and potential applications in human dietary and medical interventions. Interestingly, the RA chemotype is present sporadically in multiple taxa of flowering plants as well as some hornworts and ferns, prompting the question whether its biosynthesis arose independently across different lineages. Here we report the elucidation of the RA biosynthetic pathway in Phacelia campanularia (desert bells). This species represents the borage family (Boraginaceae), an RA-producing family closely related to the Lamiaceae within the Lamiids clade. Using a multi-omics approach in combination with functional characterization of candidate genes both in vitro and in vivo, we found that RA biosynthesis in P. campanularia involves specific activities of a BAHD acyltransferase and two cytochrome P450 hydroxylases. Further phylogenetic and comparative structure-function analyses of the P. campanularia RA biosynthetic enzymes clearly indicate that RA biosynthesis has evolved independently at least twice in the Lamiids, an exemplary case of chemotypic convergence through disparate evolutionary trajectories.

Greenlip Abalone (Haliotis laevigata) Genome and Protein Analysis Provides Insights into Maturation and Spawning

Wild abalone (Family Haliotidae) populations have been severely affected by commercial fishing, poaching, anthropogenic pollution, environment and climate changes. These issues have stimulated an increase in aquaculture production; however production growth has been slow due to a lack of genetic knowledge and resources. We have sequenced a draft genome for the commercially important temperate Australian ‘greenlip’ abalone (Haliotis laevigata, Donovan 1808) and generated 11 tissue transcriptomes from a female adult abalone. Phylogenetic analysis of the greenlip abalone with reference to the Pacific abalone (Haliotis discus hannai) indicates that these abalone species diverged approximately 71 million years ago. This study presents an in-depth analysis into the features of reproductive dysfunction, where we provide the putative biochemical messenger components (neuropeptides) that may regulate reproduction including gonad maturation and spawning. Indeed, we isolate the egg-laying hormone neuropeptide and under trial conditions induce spawning at 80% efficiency. Altogether, we provide a solid platform for further studies aimed at stimulating advances in abalone aquaculture production. The H. laevigata genome and resources are made available to the public on the abalone ‘omics website, http://abalonedb.org.

Peptidoglycan Production by an Insect-Bacterial Mosaic

Peptidoglycan (PG) is a defining feature of bacteria, involved in cell division, shape, and integrity. We previously reported that several genes related to PG biosynthesis were horizontally transferred from bacteria to the nuclear genome of mealybugs. Mealybugs are notable for containing a nested bacteria-within-bacterium endosymbiotic structure in specialized insect cells, where one bacterium, Moranella, lives in the cytoplasm of another bacterium, Tremblaya. Here we show that horizontally transferred genes on the mealybug genome work together with genes retained on the Moranella genome to produce a PG layer exclusively at the Moranella cell periphery. Furthermore, we show that an insect protein encoded by a horizontally transferred gene of bacterial origin is transported into the Moranella cytoplasm. These results provide a striking parallel to the genetic and biochemical mosaicism found in organelles, and prove that multiple horizontally transferred genes can become integrated into a functional pathway distributed between animal and bacterial endosymbiont genomes.

•

Mealybugs have two bacterial endosymbionts; one symbiont lives inside the other

•

The mealybug genome has acquired some bacterial peptidoglycan (PG)-related genes

•

This insect-symbiont mosaic pathway produces a PG layer at the innermost symbiont

•

Endosymbionts and organelles have evolved similar levels of biochemical integration