Repertoire-wide gene structure analyses: a case study comparing automatically predicted and manually annotated gene models

Manual correction of genome annotation improved alternative splicing identification of Artemisia annua

Gene annotation is essential for genome-based studies. However, algorithm-based genome annotation is difficult to fully and correctly reveal genomic information, especially for species with complex genomes. Artemisia annua L. is the only commercial resource of artemisinin production though the content of artemisinin is still to be improved. Genome-based genetic modification and breeding are useful strategies to boost artemisinin content and therefore, ensure the supply of artemisinin and reduce costs, but better gene annotation is urgently needed. In this study, we manually corrected the newly released genome annotation of A. annua using second- and third-generation transcriptome data. We found that incorrect gene information may lead to differences in structural, functional, and expression levels compared to the original expectations. We also identified alternative splicing events and found that genome annotation information impacted identifying alternative splicing genes. We further demonstrated that genome annotation information and alternative splicing could affect gene expression estimation and gene function prediction. Finally, we provided a valuable version of A. annua genome annotation and demonstrated the importance of gene annotation in future research.

Systematic Characterization of GATA Transcription Factors in Liriodendron chinense and Functional Validation in Abiotic Stresses

The Liriodendron chinense in the Magnoliaceae family is an endangered tree species useful for its socio-economic and ecological benefits. Abiotic stresses (cold, heat, and drought stress), among other factors, affect its growth, development, and distribution. However, GATA transcription factors (TFs) respond to various abiotic stresses and play a significant role in plant acclimatization to abiotic stresses. To determine the function of GATA TFs in L. chinense, we investigated the GATA genes in the genome of L. chinense. In this study, a total of 18 GATA genes were identified, which were randomly distributed on 12 of the total 17 chromosomes. These GATA genes clustered together in four separate groups based on their phylogenetic relationships, gene structures, and domain conservation arrangements. Detailed interspecies phylogenetic analyses of the GATA gene family demonstrated a conservation of the GATAs and a probable diversification that prompted gene diversification in plant species. In addition, the LcGATA gene family was shown to be evolutionarily closer to that of O. sativa, giving an insight into the possible LcGATA gene functions. Investigations of LcGATA gene duplication showed four gene duplicate pairs by the segmental duplication event, and these genes were a result of strong purified selection. Analysis of the cis-regulatory elements demonstrated a significant representation of the abiotic stress elements in the promoter regions of the LcGATA genes. Additional gene expressions through transcriptome and qPCR analyses revealed a significant upregulation of LcGATA17, and LcGATA18 in various stresses, including heat, cold, and drought stress in all time points analyzed. We concluded that the LcGATA genes play a pivotal role in regulating abiotic stress in L. chinense. In summary, our results provide new insights into understanding of the LcGATA gene family and their regulatory functions during abiotic stresses.

Persistent Parental RNAi in the Beetle Tribolium castaneum Involves Maternal Transmission of Long Double-Stranded RNA

Parental RNA interference (pRNAi) is a powerful and widely used method for gene-specific knockdown. Yet in insects its efficacy varies between species, and how the systemic response is transmitted from mother to offspring remains elusive. Using the beetle Tribolium castaneum, an RT-qPCR strategy to distinguish the presence of double-stranded RNA (dsRNA) from endogenous mRNA is reported. It is found that injected dsRNA is directly transmitted into the egg and persists throughout embryogenesis. Despite this depletion of dsRNA from the mother, it is shown that strong pRNAi can persist for months before waning at strain-specific rates. In seeking the receptor proteins for cellular uptake of long dsRNA into the egg, a phylogenomics profiling approach of candidate proteins is also presented. A visualization strategy based on taxonomically hierarchical assessment of orthology clustering data to rapidly assess gene age and copy number changes, refined by sequence-based evidence, is demonstrated. Repeated losses of SID-1-like channel proteins in the arthropods, including wholesale loss in the Heteroptera (true bugs), which are nonetheless highly sensitive to pRNAi, are thereby documented. Overall, practical considerations for insect pRNAi against a backdrop of outstanding questions on the molecular mechanism of dsRNA transmission for long-term, systemic knockdown are elucidated.

Molecular evidence of intertidal habitats selecting for repeated ice-binding protein evolution in invertebrates

Ice-binding proteins (IBPs) have evolved independently in multiple taxonomic groups to improve their survival at sub-zero temperatures. Intertidal invertebrates in temperate and polar regions frequently encounter sub-zero temperatures, yet there is little information on IBPs in these organisms. We hypothesized that there are far more IBPs than are currently known and that the occurrence of freezing in the intertidal zone selects for these proteins. We compiled a list of genome-sequenced invertebrates across multiple habitats and a list of known IBP sequences and used BLAST to identify a wide array of putative IBPs in those invertebrates. We found that the probability of an invertebrate species having an IBP was significantly greater in intertidal species than in those primarily found in open ocean or freshwater habitats. These intertidal IBPs had high sequence similarity to fish and tick antifreeze glycoproteins and fish type II antifreeze proteins. Previously established classifiers based on machine learning techniques further predicted ice-binding activity in the majority of our newly identified putative IBPs. We investigated the potential evolutionary origin of one putative IBP from the hard-shelled mussel Mytilus coruscus and suggest that it arose through gene duplication and neofunctionalization. We show that IBPs likely readily evolve in response to freezing risk and that there is an array of uncharacterized IBPs, and highlight the need for broader laboratory-based surveys of the diversity of ice-binding activity across diverse taxonomic and ecological groups.

iLoci: robust evaluation of genome content and organization for provisional and mature genome assemblies

We introduce a new framework for genome analyses based on parsing an annotated genome assembly into distinct interval loci (iLoci), available as open-source software as part of the AEGeAn Toolkit (https://github.com/BrendelGroup/AEGeAn). We demonstrate that iLoci provide an alternative coordinate system that is robust to changes in assembly and annotation versions and facilitates granular quality control of genome data. We discuss how statistics computed on iLoci reflect various characteristics of genome content and organization and illustrate how these statistics can be used to establish a baseline for assessment of the completeness and accuracy of the data. We also introduce a well-defined measure of relative genome compactness and compute other iLocus statistics that reveal genome-wide characteristics of gene arrangements in the whole genome context. Given the fast pace of assembly/annotation updates, our AEGeAn Toolkit fills a niche in computational genomics based on deriving persistent and species-specific genome statistics. Gene structure model-centric iLoci provide a precisely defined coordinate system that can be used to store assembly/annotation updates that reflect either stable or changed assessments. Large-scale application of the approach revealed species- and clade-specific genome organization in precisely defined computational terms, promising intriguing forays into the forces of shaping genome structure as more and more genome assemblies are being deposited.

Twelve quick steps for genome assembly and annotation in the classroom

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.

Genome-enabled insights into the biology of thrips as crop pests

BACKGROUND

The western flower thrips, Frankliniella occidentalis (Pergande), is a globally invasive pest and plant virus vector on a wide array of food, fiber, and ornamental crops. The underlying genetic mechanisms of the processes governing thrips pest and vector biology, feeding behaviors, ecology, and insecticide resistance are largely unknown. To address this gap, we present the F. occidentalis draft genome assembly and official gene set.

RESULTS

We report on the first genome sequence for any member of the insect order Thysanoptera. Benchmarking Universal Single-Copy Ortholog (BUSCO) assessments of the genome assembly (size = 415.8 Mb, scaffold N50 = 948.9 kb) revealed a relatively complete and well-annotated assembly in comparison to other insect genomes. The genome is unusually GC-rich (50%) compared to other insect genomes to date. The official gene set (OGS v1.0) contains 16,859 genes, of which ~ 10% were manually verified and corrected by our consortium. We focused on manual annotation, phylogenetic, and expression evidence analyses for gene sets centered on primary themes in the life histories and activities of plant-colonizing insects. Highlights include the following: (1) divergent clades and large expansions in genes associated with environmental sensing (chemosensory receptors) and detoxification (CYP4, CYP6, and CCE enzymes) of substances encountered in agricultural environments; (2) a comprehensive set of salivary gland genes supported by enriched expression; (3) apparent absence of members of the IMD innate immune defense pathway; and (4) developmental- and sex-specific expression analyses of genes associated with progression from larvae to adulthood through neometaboly, a distinct form of maturation differing from either incomplete or complete metamorphosis in the Insecta.

CONCLUSIONS

Analysis of the F. occidentalis genome offers insights into the polyphagous behavior of this insect pest that finds, colonizes, and survives on a widely diverse array of plants. The genomic resources presented here enable a more complete analysis of insect evolution and biology, providing a missing taxon for contemporary insect genomics-based analyses. Our study also offers a genomic benchmark for molecular and evolutionary investigations of other Thysanoptera species.

A benchmark study of ab initio gene prediction methods in diverse eukaryotic organisms

Background

The draft genome assemblies produced by new sequencing technologies present important challenges for automatic gene prediction pipelines, leading to less accurate gene models. New benchmark methods are needed to evaluate the accuracy of gene prediction methods in the face of incomplete genome assemblies, low genome coverage and quality, complex gene structures, or a lack of suitable sequences for evidence-based annotations.

Results

We describe the construction of a new benchmark, called G3PO (benchmark for Gene and Protein Prediction PrOgrams), designed to represent many of the typical challenges faced by current genome annotation projects. The benchmark is based on a carefully validated and curated set of real eukaryotic genes from 147 phylogenetically disperse organisms, and a number of test sets are defined to evaluate the effects of different features, including genome sequence quality, gene structure complexity, protein length, etc. We used the benchmark to perform an independent comparative analysis of the most widely used ab initio gene prediction programs and identified the main strengths and weaknesses of the programs. More importantly, we highlight a number of features that could be exploited in order to improve the accuracy of current prediction tools.

Conclusions

The experiments showed that ab initio gene structure prediction is a very challenging task, which should be further investigated. We believe that the baseline results associated with the complex gene test sets in G3PO provide useful guidelines for future studies.