ORCAE: online resource for community annotation of eukaryotes

Phenylboronic Acid-Functionalized Micelles Dual-Targeting Boronic Acid Transporter and Polysaccharides for siRNA Delivery into Brown Algae

Brown algae play essential roles ecologically, practically, and evolutionarily because they maintain coastal areas, capture carbon dioxide, and produce valuable chemicals such as therapeutic drugs. To unlock their full potential, understanding the unique molecular biology of brown algae is imperative. Genetic engineering tools that regulate homeostasis in brown algae are essential for determining their biological mechanisms in detail. However, few methodologies have been developed to control gene expression due to the robust structural barriers of brown algae. To address this issue, we designed peptide-based, small interfering RNA (siRNA)-loaded micelles decorated with phenylboronic acid (PBA) ligands. The PBA ligands facilitated the cellular uptake of the micelles into a model brown alga, Ectocarpus siliculosus (E. Siliculosus), through chemical interaction with polysaccharides in the cell wall and biological recognition by boronic acid transporters on the plasma membrane. The micelles, featuring "kill two birds with one stone" ligands, effectively induced gene silencing related to auxin biosynthesis. As a result, the growth of E. siliculosus was temporarily inhibited without persistent genome editing. This study demonstrated the potential for exploring the characteristics of brown algae through a simple yet effective approach and presented a feasible system for delivering siRNA in brown algae.

DRANetSplicer: A Splice Site Prediction Model Based on Deep Residual Attention Networks

The precise identification of splice sites is essential for unraveling the structure and function of genes, constituting a pivotal step in the gene annotation process. In this study, we developed a novel deep learning model, DRANetSplicer, that integrates residual learning and attention mechanisms for enhanced accuracy in capturing the intricate features of splice sites. We constructed multiple datasets using the most recent versions of genomic data from three different organisms, Oryza sativa japonica, Arabidopsis thaliana and Homo sapiens. This approach allows us to train models with a richer set of high-quality data. DRANetSplicer outperformed benchmark methods on donor and acceptor splice site datasets, achieving an average accuracy of (96.57%, 95.82%) across the three organisms. Comparative analyses with benchmark methods, including SpliceFinder, Splice2Deep, Deep Splicer, EnsembleSplice, and DNABERT, revealed DRANetSplicer's superior predictive performance, resulting in at least a (4.2%, 11.6%) relative reduction in average error rate. We utilized the DRANetSplicer model trained on O. sativa japonica data to predict splice sites in A. thaliana, achieving accuracies for donor and acceptor sites of (94.89%, 94.25%). These results indicate that DRANetSplicer possesses excellent cross-organism predictive capabilities, with its performance in cross-organism predictions even surpassing that of benchmark methods in non-cross-organism predictions. Cross-organism validation showcased DRANetSplicer's excellence in predicting splice sites across similar organisms, supporting its applicability in gene annotation for understudied organisms. We employed multiple methods to visualize the decision-making process of the model. The visualization results indicate that DRANetSplicer can learn and interpret well-known biological features, further validating its overall performance. Our study systematically examined and confirmed the predictive ability of DRANetSplicer from various levels and perspectives, indicating that its practical application in gene annotation is justified.

Repeated co-option of HMG-box genes for sex determination in brown algae and animals

In many eukaryotes, genetic sex determination is not governed by XX/XY or ZW/ZZ systems but by a specialized region on the poorly studied U (female) or V (male) sex chromosomes. Previous studies have hinted at the existence of a dominant male-sex factor on the V chromosome in brown algae, a group of multicellular eukaryotes distantly related to animals and plants. The nature of this factor has remained elusive. Here, we demonstrate that an HMG-box gene acts as the male-determining factor in brown algae, mirroring the role HMG-box genes play in sex determination in animals. Over a billion-year evolutionary timeline, these lineages have independently co-opted the HMG box for male determination, representing a paradigm for evolution's ability to recurrently use the same genetic "toolkit" to accomplish similar tasks.

Contrasting roles of cytochrome P450s in amitraz and chlorfenapyr resistance in the crop pest Tetranychus urticae

The molecular mechanisms of amitraz and chlorfenapyr resistance remain only poorly understood for major agricultural pests and vectors of human diseases. This study focusses on a multi-resistant field strain of the crop pest Tetranychus urticae, which could be readily selected in the laboratory to high levels of amitraz and chlorfenapyr resistance. Toxicity experiments using tralopyril, the active toxophore of chlorfenapyr, suggested decreased activation as a likely mechanism underlying resistance. Starting from the same parental strain, transcriptome profiling revealed that a cluster of detoxifying genes was upregulated after amitraz selection, but unexpectedly downregulated after chlorfenapyr selection. Further functional validation associated the upregulation of CYP392A16 with amitraz metabolism and the downregulation of CYP392D8 with reduced activation of chlorfenapyr to tralopyril. Genetic mapping (QTL analysis by BSA) was conducted in an attempt to unravel the genetic mechanisms of expression variation and resistance. This revealed that chlorfenapyr resistance was associated with a single QTL, while 3 QTLs were uncovered for amitraz resistance. Together with the observed contrasting gene expression patterns, we argue that transcriptional regulators most likely underly the distinct expression profiles associated with resistance, but these await further functional validation.

A glutamate-gated chloride channel as the mite-specific target-site of dicofol and other diphenylcarbinol acaricides

Dicofol has been widely used to control phytophagous mites. Although dicofol is chemically related to DDT, its mode of action has remained elusive. Here, we mapped dicofol resistance in the spider mite Tetranychus urticae to two genomic regions. Each region harbored a glutamate-gated chloride channel (GluCl) gene that contained a mutation-G314D or G326E-known to confer resistance against the unrelated acaricide abamectin. Using electrophysiology assays we showed that dicofol and other diphenylcarbinol acaricides-bromopropylate and chlorobenzilate-induce persistent currents in Xenopus oocytes expressing wild-type T. urticae GluCl3 receptors and potentiate glutamate responses. In contrast, the G326E substitution abolished the agonistic activity of all three compounds. Assays with the wild-type Drosophila GluClα revealed that this receptor was unresponsive to dicofol. Homology modeling combined with ligand-docking confirmed the specificity of electrophysiology assays. Altogether, this work elucidates the mode of action of diphenylcarbinols as mite-specific agonists of GluCl.

A link between energy metabolism and plant host adaptation states in the two-spotted spider mite, Tetranychus urticae (Koch)

Energy metabolism is a highly conserved process that balances generation of cellular energy and maintenance of redox homeostasis. It consists of five interconnected pathways: glycolysis, tricarboxylic acid cycle, pentose phosphate, trans-sulfuration, and NAD+ biosynthesis pathways. Environmental stress rewires cellular energy metabolism. Type-2 diabetes is a well-studied energy metabolism rewiring state in human pancreatic β-cells where glucose metabolism is uncoupled from insulin secretion. The two-spotted spider mite, Tetranychus urticae (Koch), exhibits a remarkable ability to adapt to environmental stress. Upon transfer to unfavourable plant hosts, mites experience extreme xenobiotic stress that dramatically affects their survivorship and fecundity. However, within 25 generations, mites adapt to the xenobiotic stress and restore their fitness. Mites' ability to withstand long-term xenobiotic stress raises a question of their energy metabolism states during host adaptation. Here, we compared the transcriptional responses of five energy metabolism pathways between host-adapted and non-adapted mites while using responses in human pancreatic islet donors to model these pathways under stress. We found that non-adapted mites and human pancreatic β-cells responded in a similar manner to host plant transfer and diabetogenic stress respectively, where redox homeostasis maintenance was favoured over energy generation. Remarkably, we found that upon host-adaptation, mite energy metabolic states were restored to normal. These findings suggest that genes involved in energy metabolism can serve as molecular markers for mite host-adaptation.

The host plant strongly modulates acaricide resistance levels to mitochondrial complex II inhibitors in a multi-resistant field population of Tetranychus urticae

The two-spotted spider mite Tetranychus urticae is a polyphagous pest with an extraordinary ability to develop acaricide resistance. Here, we characterize the resistance mechanisms in a T. urticae population (VR-BE) collected from a Belgian tomato greenhouse, where the grower was unsuccessful in chemically controlling the mite population resulting in crop loss. Upon arrival in the laboratory, the VR-BE population was established both on bean and tomato plants as hosts. Toxicity bioassays on both populations confirmed that the population was highly multi-resistant, recording resistance to 12 out of 13 compounds tested from various mode of action groups. DNA sequencing revealed the presence of multiple target-site resistance mutations, but these could not explain resistance to all compounds. In addition, striking differences in toxicity for six acaricides were observed between the populations on bean and tomato. The highest difference was recorded for the complex II inhibitors cyenopyrafen and cyflumetofen, which were 4.4 and 3.3-fold less toxic for VR-BE mites on tomato versus bean. PBO synergism bioassays suggested increased P450 based detoxification contribute to the host-dependent toxicity. Given the involvement of increased detoxification, we subsequently determined genome-wide gene expression levels of VR-BE on both hosts, in comparison to a reference susceptible population, revealing overexpression of a large set of detoxification genes in VR-BE on both hosts compared to the reference. In addition, a number of mainly detoxification genes with higher expression in VR-BE on tomato compared to bean was identified, including several cytochrome P450s. Together, our work suggests that multi-resistant field populations can accumulate a striking number of target-site resistance mutations. We also show that the host plant can have a profound effect on the P450-associated resistance levels to cyenopyrafen and cyflumetofen.

A nuclear receptor HR96-related gene underlies large trans-driven differences in detoxification gene expression in a generalist herbivore

The role, magnitude, and molecular nature of trans-driven expression variation underlying the upregulation of detoxification genes in pesticide resistant arthropod populations has remained enigmatic. In this study, we performed expression quantitative trait locus (eQTL) mapping (n = 458) between a pesticide resistant and a susceptible strain of the generalist herbivore and crop pest Tetranychus urticae. We found that a single trans eQTL hotspot controlled large differences in the expression of a subset of genes in different detoxification gene families, as well as other genes associated with host plant use. As established by additional genetic approaches including RNAi gene knockdown, a duplicated gene with a nuclear hormone receptor HR96-related ligand-binding domain was identified as causal for the expression differences between strains. The presence of a large family of HR96-related genes in T. urticae may enable modular control of detoxification and host plant use genes, facilitating this species' known and rapid evolution to diverse pesticides and host plants.

Identification and characterization of CYP71 subclade cytochrome P450 enzymes involved in the biosynthesis of bitterness compounds in Cichorium intybus

Industrial chicory (Cichorium intybus var. sativum) and witloof (C. intybus var. foliosum) are crops with an important economic value, mainly cultivated for inulin production and as a leafy vegetable, respectively. Both crops are rich in nutritionally relevant specialized metabolites with beneficial effects for human health. However, their bitter taste, caused by the sesquiterpene lactones (SLs) produced in leaves and taproot, limits wider applications in the food industry. Changing the bitterness would thus create new opportunities with a great economic impact. Known genes encoding enzymes involved in the SL biosynthetic pathway are GERMACRENE A SYNTHASE (GAS), GERMACRENE A OXIDASE (GAO), COSTUNOLIDE SYNTHASE (COS) and KAUNIOLIDE SYNTHASE (KLS). In this study, we integrated genome and transcriptome mining to further unravel SL biosynthesis. We found that C. intybus SL biosynthesis is controlled by the phytohormone methyl jasmonate (MeJA). Gene family annotation and MeJA inducibility enabled the pinpointing of candidate genes related with the SL biosynthetic pathway. We specifically focused on members of subclade CYP71 of the cytochrome P450 family. We verified the biochemical activity of 14 C. intybus CYP71 enzymes transiently produced in Nicotiana benthamiana and identified several functional paralogs for each of the GAO, COS and KLS genes, pointing to redundancy in and robustness of the SL biosynthetic pathway. Gene functionality was further analyzed using CRISPR/Cas9 genome editing in C. intybus. Metabolite profiling of mutant C. intybus lines demonstrated a successful reduction in SL metabolite production. Together, this study increases our insights into the C. intybus SL biosynthetic pathway and paves the way for the engineering of C. intybus bitterness.

Industrial chicory genome gives insights into the molecular timetable of anther development and male sterility

Industrial chicory (Cichorium intybus var. sativum) is a biannual crop mostly cultivated for extraction of inulin, a fructose polymer used as a dietary fiber. F1 hybrid breeding is a promising breeding strategy in chicory but relies on stable male sterile lines to prevent self-pollination. Here, we report the assembly and annotation of a new industrial chicory reference genome. Additionally, we performed RNA-Seq on subsequent stages of flower bud development of a fertile line and two cytoplasmic male sterile (CMS) clones. Comparison of fertile and CMS flower bud transcriptomes combined with morphological microscopic analysis of anthers, provided a molecular understanding of anther development and identified key genes in a range of underlying processes, including tapetum development, sink establishment, pollen wall development and anther dehiscence. We also described the role of phytohormones in the regulation of these processes under normal fertile flower bud development. In parallel, we evaluated which processes are disturbed in CMS clones and could contribute to the male sterile phenotype. Taken together, this study provides a state-of-the-art industrial chicory reference genome, an annotated and curated candidate gene set related to anther development and male sterility as well as a detailed molecular timetable of flower bud development in fertile and CMS lines.

Low Spontaneous Mutation Rate in Complex Multicellular Eukaryotes with a Haploid-Diploid Life Cycle

The spontaneous mutation rate µ is a crucial parameter to understand evolution and biodiversity. Mutation rates are highly variable across species, suggesting that µ is susceptible to selection and drift and that species life cycle and life history may impact its evolution. In particular, asexual reproduction and haploid selection are expected to affect the mutation rate, but very little empirical data are available to test this expectation. Here, we sequence 30 genomes of a parent-offspring pedigree in the model brown alga Ectocarpus sp.7, and 137 genomes of an interspecific cross of the closely related brown alga Scytosiphon to have access to the spontaneous mutation rate of representative organisms of a complex multicellular eukaryotic lineage outside animals and plants, and to evaluate the potential impact of life cycle on the mutation rate. Brown algae alternate between a haploid and a diploid stage, both multicellular and free living, and utilize both sexual and asexual reproduction. They are, therefore, excellent models to empirically test expectations of the effect of asexual reproduction and haploid selection on mutation rate evolution. We estimate that Ectocarpus has a base substitution rate of µbs = 4.07 × 10-10 per site per generation, whereas the Scytosiphon interspecific cross had µbs = 1.22 × 10-9. Overall, our estimations suggest that these brown algae, despite being multicellular complex eukaryotes, have unusually low mutation rates. In Ectocarpus, effective population size (Ne) could not entirely explain the low µbs. We propose that the haploid-diploid life cycle, combined with extensive asexual reproduction, may be additional key drivers of the mutation rate in these organisms.

Increased metabolism in combination with the novel cytochrome b target-site mutation L258F confers cross-resistance between the Qo inhibitors acequinocyl and bifenazate in Tetranychus urticae

Acequinocyl and bifenazate are potent acaricides acting at the Q_o site of complex III of the electron transport chain, but frequent applications of these acaricides have led to the development of resistance in spider mites. Target-site resistance caused by mutations in the conserved cd1- and ef-helices of the Q_o pocket of cytochrome b has been elucidated as the main resistance mechanism. We therefore monitored Q_o pocket mutations in European field populations of Tetranychus urticae and uncovered a new mutation, L258F. The role of this mutation was validated by revealing patterns of maternal inheritance and by the independently replicated introgression in an unrelated susceptible genetic background. However, the parental strain exhibited higher resistance levels than conferred by the mutation alone in isogenic lines, especially for acequinocyl, implying the involvement of strong additional resistance mechanisms. This was confirmed by revealing a polygenic inheritance pattern with classical genetic crosses and via synergism experiments. Therefore, a genome-wide expression analysis was conducted that identified a number of highly overexpressed detoxification genes, including many P450s. Functional expression revealed that the P450 CYP392A11 can metabolize bifenazate by hydroxylation of the ring structure. In conclusion, the novel cytochrome b target-site mutation L258F was uncovered in a recently collected field strain and its role in acequinocyl and bifenazate resistance was validated. However, the high level of resistance in this strain is most likely caused by a combination of target-site resistance and P450-based increased detoxification, potentially acting in synergism.

Intraspecific diversity in the mechanisms underlying abamectin resistance in a cosmopolitan pest

Pesticide resistance relies on a myriad of mechanisms, ranging from single mutations to a complex and polygenic architecture, and it involves mechanisms such as target-site insensitivity, metabolic detoxification, or a combination of these, with either additive or synergistic effects. Several resistance mechanisms against abamectin, a macrocyclic lactone widely used in crop protection, have been reported in the cosmopolitan pest Tetranychus urticae. However, it has been shown that a single mechanism cannot account for the high levels of abamectin resistance found across different mite populations. Here, we used experimental evolution combined with bulked segregant analyses to map quantitative trait loci (QTL) associated with abamectin resistance in two genetically unrelated populations of T. urticae. In these two independent QTL mapping experiments, three and four QTLs were identified, of which three were shared between experiments. Shared QTLs contained genes encoding subunits of the glutamate-gated chloride channel (GluCl) and harboured previously reported mutations, including G314D in GluCl1 and G326E in GluCl3, but also novel resistance candidate loci, including DNA helicases and chemosensory receptors. Surprisingly, the fourth QTL, present only in only one of the experiments and thus unique for one resistant parental line, revealed a non-functional variant of GluCl2, suggesting gene knock-out as resistance mechanism. Our study uncovers the complex basis of abamectin resistance, and it highlights the intraspecific diversity of genetic mechanisms underlying resistance in a cosmopolitan pest.

RNA viruses alter house dust mite physiology and allergen production with no detected consequences for allergenicity

RNA viruses have recently been detected in association with house dust mites, including laboratory cultures, dust samples, and mite-derived pharmaceuticals used for allergy diagnosis. This study aimed to assess the incidence of viral infection on Dermatophagoides pteronyssinus physiology and on the allergenic performance of extracts derived from its culture. Transcriptional changes between genetically identical control and virus-infected mite colonies were analysed by RNAseq with the support of a new D. pteronyssinus high-quality annotated genome (56.8 Mb, 108 scaffolds, N50 = 2.73 Mb, 96.7% BUSCO-completeness). Extracts of cultures and bodies from both colonies were compared by inspecting major allergen accumulation by enzyme-linked immunosorbent assay (ELISA), allergen-related enzymatic activities by specific assays, airway inflammation in a mouse model of allergic asthma, and binding to allergic patient's sera IgE by ImmunoCAP. Viral infection induced a significant transcriptional response, including several immunity and stress-response genes, and affected the expression of seven allergens, putative isoallergens and allergen orthologs. Major allergens were unaffected except for Der p 23 that was upregulated, increasing ELISA titers up to 29% in infected-mite extracts. By contrast, serine protease allergens Der p 3, 6 and 9 were downregulated, being trypsin and chymotrypsin enzymatic activities reduced up to 21% in extracts. None of the parameters analysed in our mouse model, nor binding to human IgE were significantly different when comparing control and infected-mite extracts. Despite the described physiological impact of viral infection on the mites, no significant consequences for the allergenicity of derived extracts or their practical use in allergy diagnosis have been detected.

Differences by origin in methylome suggest eco-phenotypes in the kelp Saccharina latissima

Most kelp species are of high ecological and economic importance worldwide, but are highly susceptible to rising ocean temperatures due to their sessile lifestyle. Due to interference with reproduction, development and growth, natural kelp forests have vanished in multiple regions after extreme summer heat waves. Furthermore, increasing temperatures are likely to decrease biomass production and, thus, reduce production security of farmed kelp. Epigenetic variation, and cytosine methylation as a heritable epigenetic trait, is a rapid means of acclimation and adaptation to environmental conditions, including temperature. While the first methylome of brown macroalgae has been recently described in the kelp Saccharina japonica, its functional relevance and contribution to environmental acclimation is currently unknown. The main objective of our study was to identify the importance of the methylome in the congener kelp species Saccharina latissima for temperature acclimation. Our study is the first to compare DNA methylation in kelp between wild populations of different latitudinal origin, and the first to investigate the effect of cultivation and rearing temperature on genome-wide cytosine methylation. Origin appears to determine many traits in kelp, but it is unknown to what extent the effects of thermal acclimation may be overruled by lab-related acclimation. Our results suggest that seaweed hatchery conditions have strong effects on the methylome and, thus, putatively on the epigenetically controlled characteristics of young kelp sporophytes. However, culture origin could best explain epigenetic differences in our samples suggesting that epigenetic mechanisms contribute to local adaptation of eco-phenotypes. Our study is a first step to understand whether DNA methylation marks (via their effect on gene regulation) may be used as biological regulators to enhance production security and kelp restoration success under rising temperatures, and highlights the importance to match hatchery conditions to origin.

Trans-driven variation in expression is common among detoxification genes in the extreme generalist herbivore Tetranychus urticae

The extreme adaptation potential of the generalist herbivore Tetranychus urticae (the two-spotted spider mite) to pesticides as well as diverse host plants has been associated with clade-specific gene expansions in known detoxifying enzyme families, and with extensive and rapid transcriptional responses. However, how this broad transcriptional potential is regulated remains largely unknown. Using a parental/F1 design in which four inbred strains were crossed to a common inbred strain, we assessed the genetic basis and inheritance of gene expression variation in T. urticae. Mirroring known phenotypic variation in the progenitor strains of the inbreds, we confirmed that the inbred strains we created were genetically distinct, varied markedly in pesticide resistance, and also captured variation in host plant fitness as is commonly observed in this species. By examining differences in gene expression between parents and allele-specific expression in F1s, we found that variation in RNA abundance was more often explained in trans as compared to cis, with the former associated with dominance in inheritance. Strikingly, in a gene ontology analysis, detoxification genes of the cytochrome P450 monooxygenase (CYP) family, as well as dioxygenases (DOGs) acquired from horizontal gene transfer from fungi, were specifically enriched at the extremes of trans-driven up- and downregulation. In particular, multiple CYPs and DOGs with broad substrate-specificities for pesticides or plant specialized compounds were exceptionally highly upregulated as a result of trans-regulatory variation, or in some cases synergism of cis and trans, in the most multi-pesticide resistant strains. Collectively, our findings highlight the potential importance of trans-driven expression variation in genes associated with xenobiotic metabolism and host plant use for rapid adaptation in T. urticae, and also suggests modular control of these genes, a regulatory architecture that might ameliorate negative pleiotropic effects.

High-Resolution Genetic Mapping Combined with Transcriptome Profiling Reveals That Both Target-Site Resistance and Increased Detoxification Confer Resistance to the Pyrethroid Bifenthrin in the Spider Mite Tetranychus urticae

Pyrethroids are widely applied insecticides in agriculture, but their frequent use has provoked many cases of resistance, in which mutations in the voltage-gated sodium channel (VGSC), the pyrethroid target-site, were shown to play a major role. However, for the spider mite Tetranychus urticae, it has also been shown that increased detoxification contributes to resistance against the pyrethroid bifenthrin. Here, we performed QTL-mapping to identify the genomic loci underlying bifenthrin resistance in T. urticae. Two loci on chromosome 1 were identified, with the VGSC gene being located near the second QTL and harboring the well-known L1024V mutation. In addition, the presence of an L925M mutation in the VGSC of a highly bifenthrin-resistant strain and its loss in its derived, susceptible, inbred line indicated the importance of target-site mutations in bifenthrin resistance. Further, RNAseq experiments revealed that genes encoding detoxification enzymes, including carboxyl/choline esterases (CCEs), cytochrome P450 monooxygenases and UDP-glycosyl transferases (UGTs), were overexpressed in resistant strains. Toxicity bioassays with bifenthrin (ester pyrethroid) and etofenprox (non-ester pyrethroid) also indicated a possible role for CCEs in bifenthrin resistance. A selection of CCEs and UGTs were therefore functionally expressed, and CCEinc18 was shown to metabolize bifenthrin, while teturUGT10 could glycosylate bifenthrin-alcohol. To conclude, our findings suggest that both target-site and metabolic mechanisms underlie bifenthrin resistance in T. urticae, and these might synergize high levels of resistance.

Improving helminth genome resources in the post-genomic era

Rapid advancement in high-throughput sequencing and analytical approaches has seen a steady increase in the generation of genomic resources for helminth parasites. Now, helminth genomes and their annotations are a cornerstone of numerous efforts to compare genetic and transcriptomic variation, from single cells to populations of globally distributed parasites, to genome modifications to understand gene function. Our understanding of helminths is increasingly reliant on these genomic resources, which are primarily static once published and vary widely in quality and completeness between species. This article seeks to highlight the cause and effect of this variation and argues for the continued improvement of these genomic resources - even after their publication - which is necessary to provide a more accurate and complete understanding of the biology of these important pathogens.

DeepGenGrep: a general deep learning-based predictor for multiple genomic signals and regions

MOTIVATION

Accurate annotation of different genomic signals and regions (GSRs) from DNA sequences is fundamentally important for understanding gene structure, regulation and function. Numerous efforts have been made to develop machine learning-based predictors for in silico identification of GSRs. However, it remains a great challenge to identify GSRs as the performance of most existing approaches is unsatisfactory. As such, it is highly desirable to develop more accurate computational methods for GSRs prediction.

RESULTS

In this study, we propose a general deep learning framework termed DeepGenGrep, a general predictor for the systematic identification of multiple different GSRs from genomic DNA sequences. DeepGenGrep leverages the power of hybrid neural networks comprising a three-layer convolutional neural network and a two-layer long short-term memory to effectively learn useful feature representations from sequences. Benchmarking experiments demonstrate that DeepGenGrep outperforms several state-of-the-art approaches on identifying polyadenylation signals, translation initiation sites and splice sites across four eukaryotic species including Homo sapiens, Mus musculus, Bos taurus and Drosophila melanogaster. Overall, DeepGenGrep represents a useful tool for the high-throughput and cost-effective identification of potential GSRs in eukaryotic genomes.

AVAILABILITY AND IMPLEMENTATION

The webserver and source code are freely available at http://bigdata.biocie.cn/deepgengrep/home and Github (https://github.com/wx-cie/DeepGenGrep/).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

β-Cyanoalanine synthase protects mites against Arabidopsis defenses

Glucosinolates are antiherbivory chemical defense compounds in Arabidopsis (Arabidopsis thaliana). Specialist herbivores that feed on brassicaceous plants have evolved various mechanisms aimed at preventing the formation of toxic isothiocyanates. In contrast, generalist herbivores typically detoxify isothiocyanates through glutathione conjugation upon exposure. Here, we examined the response of an extreme generalist herbivore, the two-spotted spider mite Tetranychus urticae (Koch), to indole glucosinolates. Tetranychus urticae is a composite generalist whose individual populations have a restricted host range but have an ability to rapidly adapt to initially unfavorable plant hosts. Through comparative transcriptomic analysis of mite populations that have differential susceptibilities to Arabidopsis defenses, we identified β-cyanoalanine synthase of T. urticae (TuCAS), which encodes an enzyme with dual cysteine and β-cyanoalanine synthase activities. We combined Arabidopsis genetics, chemical complementation and mite reverse genetics to show that TuCAS is required for mite adaptation to Arabidopsis through its β-cyanoalanine synthase activity. Consistent with the β-cyanoalanine synthase role in detoxification of hydrogen cyanide (HCN), we discovered that upon mite herbivory, Arabidopsis plants release HCN. We further demonstrated that indole glucosinolates are sufficient for cyanide formation. Overall, our study uncovered Arabidopsis defenses that rely on indole glucosinolate-dependent cyanide for protection against mite herbivory. In response, Arabidopsis-adapted mites utilize the β-cyanoalanine synthase activity of TuCAS to counter cyanide toxicity, highlighting the mite's ability to activate resistant traits that enable this extreme polyphagous herbivore to exploit cyanogenic host plants.

Hydroxynitrile lyase defends Arabidopsis against Tetranychus urticae

Plant-pest interactions involve multifaceted processes encompassing a complex crosstalk of pathways, molecules, and regulators aimed at overcoming defenses developed by each interacting organism. Among plant defensive compounds against phytophagous arthropods, cyanide-derived products are toxic molecules that directly target pest physiology. Here, we identified the Arabidopsis (Arabidopsis thaliana) gene encoding hydroxynitrile lyase (AtHNL, At5g10300) as one gene induced in response to spider mite (Tetranychus urticae) infestation. AtHNL catalyzes the reversible interconversion between cyanohydrins and derived carbonyl compounds with free cyanide. AtHNL loss- and gain-of-function Arabidopsis plants showed that specific activity of AtHNL using mandelonitrile as substrate was higher in the overexpressing lines than in wild-type (WT) and mutant lines. Concomitantly, mandelonitrile accumulated at higher levels in mutant lines than in WT plants and was significantly reduced in the AtHNL overexpressing lines. After mite infestation, mandelonitrile content increased in WT and overexpressing plants but not in mutant lines, while hydrogen cyanide (HCN) accumulated in the three infested Arabidopsis genotypes. Feeding bioassays demonstrated that the AtHNL gene participated in Arabidopsis defense against T. urticae. The reduced leaf damage detected in the AtHNL overexpressing lines reflected the mite's reduced ability to feed on leaves, which consequently restricted mite fecundity. In turn, mites upregulated TuCAS1 encoding β-cyanoalanine synthase to avoid the respiratory damage produced by HCN. This detoxification effect was functionally demonstrated by reduced mite fecundity observed when dsRNA-TuCAS-treated mites fed on WT plants and hnl1 mutant lines. These findings add more players in the Arabidopsis-T. urticae interplay to overcome mutual defenses.

A functionally conserved STORR gene fusion in Papaver species that diverged 16.8 million years ago

The STORR gene fusion event is considered essential for the evolution of the promorphinan/morphinan subclass of benzylisoquinoline alkaloids (BIAs) in opium poppy as the resulting bi-modular protein performs the isomerization of (S)- to (R)-reticuline essential for their biosynthesis. Here, we show that of the 12 Papaver species analysed those containing the STORR gene fusion also contain promorphinans/morphinans with one important exception. P. californicum encodes a functionally conserved STORR but does not produce promorphinans/morphinans. We also show that the gene fusion event occurred only once, between 16.8-24.1 million years ago before the separation of P. californicum from other Clade 2 Papaver species. The most abundant BIA in P. californicum is (R)-glaucine, a member of the aporphine subclass of BIAs, raising the possibility that STORR, once evolved, contributes to the biosynthesis of more than just the promorphinan/morphinan subclass of BIAs in the Papaveraceae.

Intradiol ring cleavage dioxygenases from herbivorous spider mites as a new detoxification enzyme family in animals

Background

Generalist herbivores such as the two-spotted spider mite Tetranychus urticae thrive on a wide variety of plants and can rapidly adapt to novel hosts. What traits enable polyphagous herbivores to cope with the diversity of secondary metabolites in their variable plant diet is unclear. Genome sequencing of T. urticae revealed the presence of 17 genes that code for secreted proteins with strong homology to “intradiol ring cleavage dioxygenases (DOGs)” from bacteria and fungi, and phylogenetic analyses show that they have been acquired by horizontal gene transfer from fungi. In bacteria and fungi, DOGs have been well characterized and cleave aromatic rings in catecholic compounds between adjacent hydroxyl groups. Such compounds are found in high amounts in solanaceous plants like tomato, where they protect against herbivory. To better understand the role of this gene family in spider mites, we used a multi-disciplinary approach to functionally characterize the various T. urticae DOG genes.

Results

We confirmed that DOG genes were present in the T. urticae genome and performed a phylogenetic reconstruction using transcriptomic and genomic data to advance our understanding of the evolutionary history of spider mite DOG genes. We found that DOG expression differed between mites from different plant hosts and was induced in response to jasmonic acid defense signaling. In consonance with a presumed role in detoxification, expression was localized in the mite’s gut region. Silencing selected DOGs expression by dsRNA injection reduced the mites’ survival rate on tomato, further supporting a role in mitigating the plant defense response. Recombinant purified DOGs displayed a broad substrate promiscuity, cleaving a surprisingly wide array of aromatic plant metabolites, greatly exceeding the metabolic capacity of previously characterized microbial DOGs.

Conclusion

Our findings suggest that the laterally acquired spider mite DOGs function as detoxification enzymes in the gut, disarming plant metabolites before they reach toxic levels. We provide experimental evidence to support the hypothesis that this proliferated gene family in T. urticae is causally linked to its ability to feed on an extremely wide range of host plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01323-1.

QTL mapping suggests that both cytochrome P450-mediated detoxification and target-site resistance are involved in fenbutatin oxide resistance in Tetranychus urticae

The organotin acaricide fenbutatin oxide (FBO) - an inhibitor of mitochondrial ATP-synthase - has been one of the most extensively used acaricides for the control of spider mites, and is still in use today. Resistance against FBO has evolved in many regions around the world but only few studies have investigated the molecular and genetic mechanisms of resistance to organotin acaricides. Here, we found that FBO resistance is polygenic in two genetically distant, highly resistant strains of the spider mite Tetranychus urticae, MAR-AB and MR-VL. To identify the loci underlying FBO resistance, two independent bulked segregant analysis (BSA) based QTL mapping experiments, BSA MAR-AB and BSA MR-VL, were performed. Two QTLs on chromosome 1 were associated with FBO resistance in each mapping experiment. At the second QTL of BSA MAR-AB, several cytochrome P450 monooxygenase (CYP) genes were located, including CYP392E4, CYP392E6 and CYP392E11, the latter being overexpressed in MAR-AB. Synergism tests further implied a role for CYPs in FBO resistance. Subunit c of mitochondrial ATP-synthase was located near the first QTL of both mapping experiments and harbored a unique V89A mutation enriched in the resistant parents and selected BSA populations. Marker-assisted introgression into a susceptible strain demonstrated a moderate but significant effect of the V89A mutation on toxicity of organotin acaricides. The impact of the mutation on organotin inhibition of ATP synthase was also functionally confirmed by ATPase assays on mitochondrial preparations. To conclude, our findings suggest that FBO resistance in the spider mite T. urticae is a complex interplay between CYP-mediated detoxification and target-site resistance.

A H258Y mutation in subunit B of the succinate dehydrogenase complex of the spider mite Tetranychus urticae confers resistance to cyenopyrafen and pyflubumide, but likely reinforces cyflumetofen binding and toxicity

Succinate dehydrogenase (SDH) inhibitors such as cyflumetofen, cyenopyrafen and pyflubumide, are selective acaricides that control plant-feeding spider mite pests. Resistance development to SDH inhibitors has been investigated in a limited number of populations of the spider mite Tetranychus urticae and is associated with cytochrome P450 based detoxification and target-site mutations such as I260 T/V in subunit B and S56L in subunit C of SDH. Here, we report the discovery of a H258Y substitution in subunit B of SDH in a highly pyflubumide resistant population of T. urticae. As this highly conserved residue corresponds to one of the ubiquinone binding residues in fungi and bacteria, we hypothesized that H258Y could have a strong impact on SDH inhibitors toxicity. Marker assisted introgression and toxicity bioassays revealed that H258Y caused high cross resistance between cyenopyrafen and pyflubumide, but increased cyflumetofen toxicity. Resistance associated with H258Y was determined as dominant for cyenopyrafen, but recessive for pyflubumide. In vitro SDH assays with extracted H258 mitochondria showed that cyenopyrafen and the active metabolites of pyflubumide and cyflumetofen, interacted strongly with complex II. However, a clear shift in IC50s was observed for cyenopyrafen and the metabolite of pyflubumide when Y258 mitochondria were investigated. In contrast, the mutation slightly increased affinity of the cyflumetofen metabolite, likely explaining its increased toxicity for the mite lines carrying the substitution. Homology modeling and ligand docking further revealed that, although the three acaricides share a common binding motif in the Q-site of SDH, H258Y eliminated an important hydrogen bond required for cyenopyrafen and pyflubumide binding. In addition, the hydrogen bond between cyenopyrafen and Y117 in subunit D was also lost upon mutation. In contrast, cyflumetofen affinity was enhanced due to an additional hydrogen bond to W215 and hydrophobic interactions with the introduced Y258 in subunit B. Altogether, our findings not only highlight the importance of the highly conserved histidine residue in the binding of SDH inhibitors, but also reveal that a resistance mutation can provide both positive and negative cross-resistance within the same acaricide mode of action group.

The APAF1_C/WD40 repeat domain-encoding gene from the sea lettuce Ulva mutabilis sheds light on the evolution of NB-ARC domain-containing proteins in green plants

We advance Ulva's genetic tractability and highlight its value as a model organism by characterizing its APAF1_C/WD40 domain-encoding gene, which belongs to a family that bears homology to R genes. The multicellular chlorophyte alga Ulva mutabilis (Ulvophyceae, Ulvales) is native to coastal ecosystems worldwide and attracts both high socio-economic and scientific interest. To further understand the genetic mechanisms that guide its biology, we present a protocol, based on adapter ligation-mediated PCR, for retrieving flanking sequences in U. mutabilis vector-insertion mutants. In the created insertional library, we identified a null mutant with an insertion in an apoptotic protease activating factor 1 helical domain (APAF1_C)/WD40 repeat domain-encoding gene. Protein domain architecture analysis combined with phylogenetic analysis revealed that this gene is a member of a subfamily that arose early in the evolution of green plants (Viridiplantae) through the acquisition of a gene that also encoded N-terminal nucleotide-binding adaptor shared by APAF-1, certain R-gene products and CED-4 (NB-ARC) and winged helix-like (WH-like) DNA-binding domains. Although phenotypic analysis revealed no mutant phenotype, gene expression levels in control plants correlated to the presence of bacterial symbionts, which U. mutabilis requires for proper morphogenesis. In addition, our analysis led to the discovery of a putative Ulva nucleotide-binding site and leucine-rich repeat (NBS-LRR) Resistance protein (R-protein), and we discuss how the emergence of these R proteins in green plants may be linked to the evolution of the APAF1_C/WD40 protein subfamily.

Environmental and Circadian Regulation Combine to Shape the Rhythmic Selenoproteome

The circadian clock orchestrates an organism’s endogenous processes with environmental 24 h cycles. Redox homeostasis and the circadian clock regulate one another to negate the potential effects of our planet’s light/dark cycle on the generation of reactive oxygen species (ROS) and attain homeostasis. Selenoproteins are an important class of redox-related enzymes that have a selenocysteine residue in the active site. This study reports functional understanding of how environmental and endogenous circadian rhythms integrate to shape the selenoproteome in a model eukaryotic cell. We mined quantitative proteomic data for the 24 selenoproteins of the picoeukaryote Ostreococcus tauri across time series, under environmentally rhythmic entrained conditions of light/dark (LD) cycles, compared to constant circadian conditions of constant light (LL). We found an overrepresentation of selenoproteins among rhythmic proteins under LL, but an underrepresentation under LD conditions. Rhythmic selenoproteins under LL that reach peak abundance later in the day showed a greater relative amplitude of oscillations than those that peak early in the day. Under LD, amplitude did not correlate with peak phase; however, we identified high-amplitude selenium uptake rhythms under LD but not LL conditions. Selenium deprivation induced strong qualitative defects in clock gene expression under LD but not LL conditions. Overall, the clear conclusion is that the circadian and environmental cycles exert differential effects on the selenoproteome, and that the combination of the two enables homeostasis. Selenoproteins may therefore play an important role in the cellular response to reactive oxygen species that form as a consequence of the transitions between light and dark.

De novo transcriptome sequencing of the northern fowl mite, Ornithonyssus sylviarum, shed light on parasitiform poultry mites evolution and its chemoreceptor repertoires

The northern fowl mite (NFM), Ornithonyssus sylviarum, and the poultry red mite (PRM), Dermanyssus gallinae, are the most serious pests of poultry, both of which have an expanding global prevalence. Research on NFM has been constrained by a lack of genomic and transcriptomic data. Here, we report and analyze the first global transcriptome data across all mite live stages and sexes. A total of 28,999 unigenes were assembled, of which 19,750 (68.10%) were annotated using seven functional databases. The biological function of these unigenes was classified using the GO, KOG, and KEGG databases. To gain insight into the chemosensory receptor-based system of parasitiform mites, we furthermore assessed the gene repertoire of gustatory receptors (GRs) and ionotropic receptors (IRs), both of which encode putative ligand-gated ion channel proteins. While these receptors are well characterized in insect model species, our understanding of chemosensory detection in mites and ticks is in its infancy. To address this paucity of data, we identified 9 IR/iGluRs and 2 GRs genes by analyzing transcriptome data in the NFM, while 9 GRs and 41 IR/iGluRs genes were annotated in the PRM genome. Taken together, the transcriptomic and genomic annotation of these two species provide a valuable reference for studies of parasitiform mites and also help to understand how chemosensory gene family expansion/contraction events may have been reshaped by an obligate parasitic lifestyle compared with their free-living closest relatives. Future studies should include additional species to validate this observation and functional characterization of the identified proteins as a step forward in identifying tools for controlling these poultry pests.

OUP accepted manuscript

Recent studies have suggested that species of the Kazachstania genus may be interesting models of yeast domestication. Among these, Kazachstania barnettii has been isolated from various microbially transformed foodstuffs such as sourdough bread and kefir. In the present work, we sequence, assemble, and annotate the complete genomes of two K. barnettii strains: CLIB 433, being one of the two reference strains for K. barnettii that was isolated as a spoilage organism in soft drink, and CLIB 1767, recently isolated from artisan bread-making sourdough. Both assemblies are of high quality with N50 statistics greater than 1.3 Mb and BUSCO score greater than 99%. An extensive comparison of the two obtained genomes revealed very few differences between the two K. barnettii strains, considering both genome structure and gene content. The proposed genome assemblies will constitute valuable references for future comparative genomic, population genomic, or transcriptomic studies of the K. barnettii species.

Adaptation of Temperate Seagrass to Arctic Light Relies on Seasonal Acclimatization of Carbon Capture and Metabolism

Due to rising global surface temperatures, Arctic habitats are becoming thermally suitable for temperate species. Whether a temperate species can immigrate into an ice-free Arctic depends on its ability to tolerate extreme seasonal fluctuations in daylength. Thus, understanding adaptations to polar light conditions can improve the realism of models predicting poleward range expansions in response to climate change. Plant adaptations to polar light have rarely been studied and remain unknown in seagrasses. If these ecosystem engineers can migrate polewards, seagrasses will enrich biodiversity, and carbon capture potential in shallow coastal regions of the Arctic. Eelgrass (Zostera marina) is the most widely distributed seagrass in the northern hemisphere. As the only seagrass species growing as far north as 70°N, it is the most likely candidate to first immigrate into an ice-free Arctic. Here, we describe seasonal (and diurnal) changes in photosynthetic characteristics, and in genome-wide gene expression patterns under strong annual fluctuations of daylength. We compared PAM measurements and RNA-seq data between two populations at the longest and shortest day of the year: (1) a Mediterranean population exposed to moderate annual fluctuations of 10-14 h daylength and (2) an Arctic population exposed to high annual fluctuations of 0-24 h daylength. Most of the gene expression specificities of the Arctic population were found in functions of the organelles (chloroplast and mitochondrion). In winter, Arctic eelgrass conserves energy by repressing respiration and reducing photosynthetic energy fluxes. Although light-reactions, and genes involved in carbon capture and carbon storage were upregulated in summer, enzymes involved in CO2 fixation and chlorophyll-synthesis were upregulated in winter, suggesting that winter metabolism relies not only on stored energy resources but also on active use of dim light conditions. Eelgrass is unable to use excessive amounts of light during summer and demonstrates a significant reduction in photosynthetic performance under long daylengths, possibly to prevent photoinhibition constrains. Our study identified key mechanisms that allow eelgrass to survive under Arctic light conditions and paves the way for experimental research to predict whether and up to which latitude eelgrass can potentially migrate polewards in response to climate change.

Deep-coverage spatiotemporal proteome of the picoeukaryote Ostreococcus tauri reveals differential effects of environmental and endogenous 24-hour rhythms

The cellular landscape changes dramatically over the course of a 24 h day. The proteome responds directly to daily environmental cycles and is additionally regulated by the circadian clock. To quantify the relative contribution of diurnal versus circadian regulation, we mapped proteome dynamics under light:dark cycles compared with constant light. Using Ostreococcus tauri, a prototypical eukaryotic cell, we achieved 85% coverage, which allowed an unprecedented insight into the identity of proteins that facilitate rhythmic cellular functions. The overlap between diurnally- and circadian-regulated proteins was modest and these proteins exhibited different phases of oscillation between the two conditions. Transcript oscillations were generally poorly predictive of protein oscillations, in which a far lower relative amplitude was observed. We observed coordination between the rhythmic regulation of organelle-encoded proteins with the nuclear-encoded proteins that are targeted to organelles. Rhythmic transmembrane proteins showed a different phase distribution compared with rhythmic soluble proteins, indicating the existence of a circadian regulatory process specific to the biogenesis and/or degradation of membrane proteins. Our observations argue that the cellular spatiotemporal proteome is shaped by a complex interaction between intrinsic and extrinsic regulatory factors through rhythmic regulation at the transcriptional as well as post-transcriptional, translational, and post-translational levels.

Holly Kay, Ellen Grünewald, et al. provide an in-depth examination of the proteome in the eukaryotic green alga, Ostreococcus tauri, under circadian constant light or cycling diurnal light-dark conditions. They observe that there is little overlap between mRNA and protein expression rhythms, or the diurnal and circadian proteome, suggesting that the cellular spatiotemporal proteome is shaped through rhythmic regulation at multiple stages of transcription and translation.

Targeted CRISPR-Cas9-based gene knockouts in the model brown alga Ectocarpus

Brown algae are an important group of multicellular eukaryotes, phylogenetically distinct from both the animal and land plant lineages. Ectocarpus has emerged as a model organism to study diverse aspects of brown algal biology, but this system currently lacks an effective reverse genetics methodology to analyse the functions of selected target genes. Here, we report that mutations at specific target sites are generated following the introduction of CRISPR-Cas9 ribonucleoproteins into Ectocarpus cells, using either biolistics or microinjection as the delivery method. Individuals with mutations affecting the ADENINE PHOSPHORIBOSYL TRANSFERASE (APT) gene were isolated following treatment with 2-fluoroadenine, and this selection system was used to isolate individuals in which mutations had been introduced simultaneously at APT and at a second gene. This double mutation approach could potentially be used to isolate mutants affecting any Ectocarpus gene, providing an effective reverse genetics tool for this model organism. The availability of this tool will significantly enhance the utility of Ectocarpus as a model organism for this ecologically and economically important group of marine organisms. Moreover, the methodology described here should be readily transferable to other brown algal species.

The genome of the extremophile Artemia provides insight into strategies to cope with extreme environments

BACKGROUND

Brine shrimp Artemia have an unequalled ability to endure extreme salinity and complete anoxia. This study aims to elucidate its strategies to cope with these stressors.

RESULTS AND DISCUSSION

Here, we present the genome of an inbred A. franciscana Kellogg, 1906. We identified 21,828 genes of which, under high salinity, 674 genes and under anoxia, 900 genes were differentially expressed (42%, respectively 30% were annotated). Under high salinity, relevant stress genes and pathways included several Heat Shock Protein and Leaf Embryogenesis Abundant genes, as well as the trehalose metabolism. In addition, based on differential gene expression analysis, it can be hypothesized that a high oxidative stress response and endocytosis/exocytosis are potential salt management strategies, in addition to the expression of major facilitator superfamily genes responsible for transmembrane ion transport. Under anoxia, genes involved in mitochondrial function, mTOR signalling and autophagy were differentially expressed. Both high salt and anoxia enhanced degradation of erroneous proteins and protein chaperoning. Compared with other branchiopod genomes, Artemia had 0.03% contracted and 6% expanded orthogroups, in which 14% of the genes were differentially expressed under high salinity or anoxia. One phospholipase D gene family, shown to be important in plant stress response, was uniquely present in both extremophiles Artemia and the tardigrade Hypsibius dujardini, yet not differentially expressed under the described experimental conditions.

CONCLUSIONS

A relatively complete genome of Artemia was assembled, annotated and analysed, facilitating research on its extremophile features, and providing a reference sequence for crustacean research.

The Mastigamoeba balamuthi Genome and the Nature of the Free-Living Ancestor of Entamoeba

The transition of free-living organisms to parasitic organisms is a mysterious process that occurs in all major eukaryotic lineages. Parasites display seemingly unique features associated with their pathogenicity; however, it is important to distinguish ancestral preconditions to parasitism from truly new parasite-specific functions. Here, we sequenced the genome and transcriptome of anaerobic free-living Mastigamoeba balamuthi and performed phylogenomic analysis of four related members of the Archamoebae, including Entamoeba histolytica, an important intestinal pathogen of humans. We aimed to trace gene histories throughout the adaptation of the aerobic ancestor of Archamoebae to anaerobiosis and throughout the transition from a free-living to a parasitic lifestyle. These events were associated with massive gene losses that, in parasitic lineages, resulted in a reduction in structural features, complete losses of some metabolic pathways, and a reduction in metabolic complexity. By reconstructing the features of the common ancestor of Archamoebae, we estimated preconditions for the evolution of parasitism in this lineage. The ancestor could apparently form chitinous cysts, possessed proteolytic enzyme machinery, compartmentalized the sulfate activation pathway in mitochondrion-related organelles, and possessed the components for anaerobic energy metabolism. After the split of Entamoebidae, this lineage gained genes encoding surface membrane proteins that are involved in host–parasite interactions. In contrast, gene gains identified in the M. balamuthi lineage were predominantly associated with polysaccharide catabolic processes. A phylogenetic analysis of acquired genes suggested an essential role of lateral gene transfer in parasite evolution (Entamoeba) and in adaptation to anaerobic aquatic sediments (Mastigamoeba).

Signatures of Transcription Factor Evolution and the Secondary Gain of Red Algae Complexity

Red algae (Rhodophyta) belong to the superphylum Archaeplastida, and are a species-rich group exhibiting diverse morphologies. Theory has it that the unicellular red algal ancestor went through a phase of genome contraction caused by adaptation to extreme environments. More recently, the classes Porphyridiophyceae, Bangiophyceae, and Florideophyceae experienced genome expansions, coinciding with an increase in morphological complexity. Transcription-associated proteins (TAPs) regulate transcription, show lineage-specific patterns, and are related to organismal complexity. To better understand red algal TAP complexity and evolution, we investigated the TAP family complement of uni- and multi-cellular red algae. We found that the TAP family complement correlates with gain of morphological complexity in the multicellular Bangiophyceae and Florideophyceae, and that abundance of the C2H2 zinc finger transcription factor family may be associated with the acquisition of morphological complexity. An expansion of heat shock transcription factors (HSF) occurred within the unicellular Cyanidiales, potentially as an adaption to extreme environmental conditions.

Genome Wide MeDIP-Seq Profiling of Wild and Cultivated Olives Trees Suggests DNA Methylation Fingerprint on the Sensory Quality of Olive Oil

Secondary metabolites are particularly important to humans due to their pharmaceutical properties. Moreover, secondary metabolites are key compounds in climate change adaptation in long-living trees. Recently, it has been described that the domestication of Olea subspecies had no major selection signature on coding variants and was mainly related to changes in gene expression. In addition, the phenotypic plasticity in Olea subspecies was linked to the activation of transposable elements in the genes neighboring. Here, we investigated the imprint of DNA methylation in the unassigned fraction of the phenotypic plasticity of the Olea subspecies, using methylated DNA immuno-precipitation sequencing (MeDIP-seq) for a high-resolution genome-wide DNA methylation profiling of leaves and fruits during fruit development in wild and cultivated olives from Turkey. Notably, the methylation profiling showed a differential DNA methylation in secondary metabolism responsible for the sensory quality of olive oil. Here, we highlight for the first time the imprint of DNA methylation in modulating the activity of the Linoleate 9S lipoxygenase in the biosynthesis of volatile aromatic compounds. Unprecedently, the current study reveals the methylation status of the olive genome during fruit ripening.

Cytogenomic analysis unveils mixed molecular evolution and recurrent chromosomal rearrangements shaping the multigene families on Schistocerca grasshopper genomes

Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat-enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long-read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth-and-death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra- and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus.

A chromosome-level Amaranthus cruentus genome assembly highlights gene family evolution and biosynthetic gene clusters that may underpin the nutritional value of this traditional crop

Traditional crops have historically provided accessible and affordable nutrition to millions of rural dwellers but have been neglected, with most modern agricultural systems over-reliant on a small number of internationally traded crops. Traditional crops are typically well-adapted to local agro-ecological conditions and many are nutrient-dense. They can play a vital role in local food systems through enhanced nutrition (particularly where diets are dominated by starch crops), food security and livelihoods for smallholder farmers, and a climate-resilient and biodiverse agriculture. Using short-read, long-read and phased sequencing technologies, we generated a high-quality chromosome-level genome assembly for Amaranthus cruentus, an under-researched crop with micronutrient- and protein-rich leaves and gluten-free seed, but lacking improved varieties, with respect to productivity and quality traits. The 370.9 Mb genome demonstrates a shared whole genome duplication with a related species, Amaranthus hypochondriacus. Comparative genome analysis indicates chromosomal loss and fusion events following genome duplication that are common to both species, as well as fission of chromosome 2 in A. cruentus alone, giving rise to a haploid chromosome number of 17 (versus 16 in A. hypochondriacus). Genomic features potentially underlying the nutritional value of this crop include two A. cruentus-specific genes with a likely role in phytic acid synthesis (an anti-nutrient), expansion of ion transporter gene families, and identification of biosynthetic gene clusters conserved within the amaranth lineage. The A. cruentus genome assembly will underpin much-needed research and global breeding efforts to develop improved varieties for economically viable cultivation and realization of the benefits to global nutrition security and agrobiodiversity.

Transcriptomic analysis of the poultry red mite, Dermanyssus gallinae, across all stages of the lifecycle

BACKGROUND

The blood feeding poultry red mite (PRM), Dermanyssus gallinae, causes substantial economic damage to the egg laying industry worldwide, and is a serious welfare concern for laying hens and poultry house workers. In this study we have investigated the temporal gene expression across the 6 stages/sexes (egg, larvae, protonymph and deutonymph, adult male and adult female) of this neglected parasite in order to understand the temporal expression associated with development, parasitic lifestyle, reproduction and allergen expression.

RESULTS

RNA-seq transcript data for the 6 stages were mapped to the PRM genome creating a publicly available gene expression atlas (on the OrcAE platform in conjunction with the PRM genome). Network analysis and clustering of stage-enriched gene expression in PRM resulted in 17 superclusters with stage-specific or multi-stage expression profiles. The 6 stage specific superclusters were clearly demarked from each other and the adult female supercluster contained the most stage specific transcripts (2725), whilst the protonymph supercluster the fewest (165). Fifteen pairwise comparisons performed between the different stages resulted in a total of 6025 Differentially Expressed Genes (DEGs) (P > 0.99). These data were evaluated alongside a Venn/Euler analysis of the top 100 most abundant genes in each stage. An expanded set of cuticle proteins and enzymes (chitinase and metallocarboxypeptidases) were identified in larvae and underpin cuticle formation and ecdysis to the protonymph stage. Two mucin/peritrophic-A salivary proteins (DEGAL6771g00070, DEGAL6824g00220) were highly expressed in the blood-feeding stages, indicating peritrophic membrane formation during feeding. Reproduction-associated vitellogenins were the most abundant transcripts in adult females whilst, in adult males, an expanded set of serine and cysteine proteinases and an epididymal protein (DEGAL6668g00010) were highly abundant. Assessment of the expression patterns of putative homologues of 32 allergen groups from house dust mites indicated a bias in their expression towards the non-feeding larval stage of PRM.

CONCLUSIONS

This study is the first evaluation of temporal gene expression across all stages of PRM and has provided insight into developmental, feeding, reproduction and survival strategies employed by this mite. The publicly available PRM resource on OrcAE offers a valuable tool for researchers investigating the biology and novel interventions of this parasite.

Untangling a Gordian knot: the role of a GluCl3 I321T mutation in abamectin resistance in Tetranychus urticae

BACKGROUND

The cys-loop ligand-gated ion channels, including the glutamate-gated chloride channel (GluCl) and GABA-gated chloride channel (Rdl) are important targets for drugs and pesticides. The macrocyclic lactone abamectin primarily targets GluCl and is commonly used to control the spider mite Tetranychus urticae, an economically important crop pest. However, abamectin resistance has been reported for multiple T. urticae populations worldwide, and in several cases was associated with the mutations G314D in GluCl1 and G326E in GluCl3. Recently, an additional I321T mutation in GluCl3 was identified in several abamectin resistant T. urticae field populations. Here, we aim to functionally validate this mutation and determine its phenotypic strength.

RESULTS

The GluCl3 I321T mutation was introgressed into a T. urticae susceptible background by marker-assisted backcrossing, revealing contrasting results in phenotypic strength, ranging from almost none to 50-fold. Next, we used CRISPR-Cas9 to introduce I321T, G314D and G326E in the orthologous Drosophila GluCl. Genome modified flies expressing GluCl I321T were threefold less susceptible to abamectin, while CRISPRed GluCl G314D and G326E flies were lethal. Last, functional analysis in Xenopus oocytes revealed that the I321T mutation might reduce GluCl3 sensitivity to abamectin, but also suggested that all three T. urticae Rdls are affected by abamectin.

CONCLUSION

Three different techniques were used to characterize the role of I321T in GluCl3 in abamectin resistance and, combining all results, our analysis suggests that the I321T mutation has a complex role in abamectin resistance. Given the reported subtle effect, additional synergistic factors in resistance warrant more investigation. © 2020 Society of Chemical Industry.

Genome sequence and genetic diversity analysis of an under-domesticated orphan crop, white fonio (Digitaria exilis)

Background

Digitaria exilis, white fonio, is a minor but vital crop of West Africa that is valued for its resilience in hot, dry, and low-fertility environments and for the exceptional quality of its grain for human nutrition. Its success is hindered, however, by a low degree of plant breeding and improvement.

Findings

We sequenced the fonio genome with long-read SMRT-cell technology, yielding a ∼761 Mb assembly in 3,329 contigs (N50, 1.73 Mb; L50, 126). The assembly approaches a high level of completion, with a BUSCO score of >99%. The fonio genome was found to be a tetraploid, with most of the genome retained as homoeologous duplications that differ overall by ∼4.3%, neglecting indels. The 2 genomes within fonio were found to have begun their independent divergence ∼3.1 million years ago. The repeat content (>49%) is fairly standard for a grass genome of this size, but the ratio of Gypsy to Copia long terminal repeat retrotransposons (∼6.7) was found to be exceptionally high. Several genes related to future improvement of the crop were identified including shattering, plant height, and grain size. Analysis of fonio population genetics, primarily in Mali, indicated that the crop has extensive genetic diversity that is largely partitioned across a north-south gradient coinciding with the Sahel and Sudan grassland domains.

Conclusions

We provide a high-quality assembly, annotation, and diversity analysis for a vital African crop. The availability of this information should empower future research into further domestication and improvement of fonio.

RNAi gene knockdown in the poultry red mite, Dermanyssus gallinae (De Geer 1778), a tool for functional genomics

Background

The avian haematophagous ectoparasite Dermanyssus gallinae, commonly known as the poultry red mite, causes significant economic losses to the egg-laying industry worldwide and also represents a significant welfare threat. Current acaricide-based controls are unsustainable due to the mite’s ability to rapidly develop resistance, thus developing a novel sustainable means of control for D. gallinae is a priority. RNA interference (RNAi)-mediated gene silencing is a valuable tool for studying gene function in non-model organisms, but is also emerging as a novel tool for parasite control.

Methods

Here we use an in silico approach to identify core RNAi pathway genes in the recently sequenced D. gallinae genome. In addition we utilise an in vitro feeding device to deliver double-stranded (ds) RNA to D. gallinae targeting the D. gallinae vATPase subunit A (Dg vATPase A) gene and monitor gene knockdown using quantitative PCR (qPCR).

Results

Core components of the small interfering RNA (siRNA) and microRNA (miRNA) pathways were identified in D. gallinae, which indicates that these gene silencing pathways are likely functional. Strikingly, the P-element-induced wimpy testis (PIWI)-interacting RNA (piRNA) pathway was absent in D. gallinae. In addition, feeding Dg vATPase A dsRNA to adult female D. gallinae resulted in silencing of the targeted gene compared to control mites fed non-specific lacZ dsRNA. In D. gallinae, dsRNA-mediated gene knockdown was rapid, being detectable 24 h after oral delivery of the dsRNA, and persisted for at least 120 h.

Conclusions

This study shows the presence of core RNAi machinery components in the D. gallinae genome. In addition, we have developed a robust RNAi methodology for targeting genes in D. gallinae that will be of value for studying genes of unknown function and validating potential control targets in D. gallinae.

Histone modifications during the life cycle of the brown alga Ectocarpus

BACKGROUND

Brown algae evolved complex multicellularity independently of the animal and land plant lineages and are the third most developmentally complex phylogenetic group on the planet. An understanding of developmental processes in this group is expected to provide important insights into the evolutionary events necessary for the emergence of complex multicellularity. Here, we focus on mechanisms of epigenetic regulation involving post-translational modifications of histone proteins.

RESULTS

A total of 47 histone post-translational modifications are identified, including a novel mark H2AZR38me1, but Ectocarpus lacks both H3K27me3 and the major polycomb complexes. ChIP-seq identifies modifications associated with transcription start sites and gene bodies of active genes and with transposons. H3K79me2 exhibits an unusual pattern, often marking large genomic regions spanning several genes. Transcription start sites of closely spaced, divergently transcribed gene pairs share a common nucleosome-depleted region and exhibit shared histone modification peaks. Overall, patterns of histone modifications are stable through the life cycle. Analysis of histone modifications at generation-biased genes identifies a correlation between the presence of specific chromatin marks and the level of gene expression.

CONCLUSIONS

The overview of histone post-translational modifications in the brown alga presented here will provide a foundation for future studies aimed at understanding the role of chromatin modifications in the regulation of brown algal genomes.

The medium-size noncoding RNA transcriptome of Ostreococcus tauri, the smallest living eukaryote, reveals a large family of small nucleolar RNAs displaying multiple genomic expression strategies

The small nucleolar RNAs (snoRNAs), essential for ribosome biogenesis, constitute a major family of medium-size noncoding RNAs (mncRNAs) in all eukaryotes. We present here, for the first time in a marine unicellular alga, the characterization of the snoRNAs family in Ostreococcus tauri, the smallest photosynthetic eukaryote. Using a transcriptomic approach, we identified 131 O. tauri snoRNAs (Ot–snoRNA) distributed in three classes: the C/D snoRNAs, the H/ACA snoRNAs and the MRP RNA. Their genomic organization revealed a unique combination of both the intronic organization of animals and the polycistronic organization of plants. Remarkably, clustered genes produced Ot–snoRNAs with unusual structures never previously described in plants. Their abundances, based on quantification of reads and northern blots, showed extreme differences in Ot–snoRNA accumulation, mainly determined by their differential stability. Most of these Ot–snoRNAs were predicted to target rRNAs or snRNAs. Seventeen others were orphan Ot–snoRNAs that would not target rRNA. These were specific to O. tauri or Mamiellophyceae and could have functions unrelated to ribosome biogenesis. Overall, these data reveal an ‘evolutionary response’ adapted to the extreme compactness of the O. tauri genome that accommodates the essential Ot–snoRNAs, developing multiple strategies to optimize their coordinated expression with a minimal cost on regulatory circuits.

Resistance to pyridaben in Canadian greenhouse populations of two-spotted spider mites, Tetranychus urticae (Koch)

Two-spotted spider mite (TSSM) Tetranychus urticae (Koch) is an important agricultural pest that causes considerable yield losses to over 150 field and greenhouse crops. Mitochondrial electron transport inhibitors (METI) acaricides are commonly used to control mite species in commercial Canadian greenhouses. Development of resistance to METIs in TSSM populations have been reported worldwide, but not until recently in Canada. The objectives of this study were to: 1) monitor the acaricide-susceptibility in greenhouse TSSM populations, and 2) investigate the resistance to pyridaben, a METI acaricide, in greenhouse resistant and pyridaben-selected (SRS) mite strains. The increased mortality to the pyridaben sub-lethal concentration (LC₃₀) when SRS mites were exposed to piperonyl butoxide (PBO), a general cytochrome P450 monooxygenase inhibitor, and higher P450 activity compared to the greenhouse strain (RS) mites, indicated that P450s may be at least partially responsible for the resistance. The molecular mechanisms of target site insensitivity-mediated resistance in the pyridaben resistant strain of TSSM were investigated by comparing the DNA sequence of NADH dehydrogenase subunits TYKY and PSST, NADH-ubiquinone oxidoreductase chain 1 and 5 (ND1, ND5) and the NADH-ubiquinone oxidoreductase subunit 49 kDa from SRS to the reference strain (SS) and RS. Despite a number of nucleotide substitutions, none correlated with the pyridaben resistance. Understanding the underlying mechanisms of TSSM adaptation to acaricides is an essential part of resistance management strategy in any IPM program. The findings of this study will encourage growers to apply acaricides with different modes of action to reduce the rate at which acaricide resistance will occur in greenhouse TSSM populations.

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.

A Genomic Perspective on the Evolutionary Diversity of the Plant Cell Wall

The plant cell wall is a complex and dynamic structure composed of numerous different molecules that play multiple roles in all aspects of plant life. Currently, a new frontier in biotechnology is opening up, which is providing new insights into the structural and functional diversity of cell walls, and is thus serving to re-emphasize the significance of cell wall divergence in the evolutionary history of plant species. The ever-increasing availability of plant genome datasets will thus provide an invaluable basis for enhancing our knowledge regarding the diversity of cell walls among different plant species. In this review, as an example of a comparative genomics approach, I examine the diverse patterns of cell wall gene families among 100 species of green plants, and illustrate the evident benefits of using genome databases for studying cell wall divergence. Given that the growth and development of all types of plant cells are intimately associated with cell wall dynamics, gaining a further understanding of the functional diversity of cell walls in relation to diverse biological events will make significant contributions to a broad range of plant sciences.

The Seagrass Methylome Is Associated With Variation in Photosynthetic Performance Among Clonal Shoots

Evolutionary theory predicts that clonal organisms are more susceptible to extinction than sexually reproducing organisms, due to low genetic variation and slow rates of evolution. In agreement, conservation management considers genetic variation as the ultimate measure of a population's ability to survive over time. However, clonal plants are among the oldest living organisms on our planet. Here, we test the hypothesis that clonal seagrass meadows display epigenetic variation that complements genetic variation as a source of phenotypic variation. In a clonal meadow of the seagrass Zostera marina, we characterized DNA methylation among 42 shoots. We also sequenced the whole genome of 10 shoots to correlate methylation patterns with photosynthetic performance under exposure to and recovery from 27°C, while controlling for somatic mutations. Here, we show for the first time that clonal seagrass shoots display DNA methylation variation that is independent from underlying genetic variation, and associated with variation in photosynthetic performance under experimental conditions. It remains unknown to what degree this association could be influenced by epigenetic responses to transplantation-related stress, given that the methylomes showed a strong shift under acclimation to laboratory conditions. The lack of untreated control samples in the heat stress experiment did not allow us to distinguish methylome shifts induced by acclimation from such induced by heat stress. Notwithstanding, the co-variation in DNA methylation and photosynthetic performance may be linked via gene expression because methylation patterns varied in functionally relevant genes involved in photosynthesis, and in the repair and prevention of heat-induced protein damage. While genotypic diversity has been shown to enhance stress resilience in seagrass meadows, we suggest that epigenetic variation plays a similar role in meadows dominated by a single genotype. Consequently, conservation management of clonal plants should consider epigenetic variation as indicator of resilience and stability.