Predicting subcellular localization of proteins based on their N-terminal amino acid sequence

Identification of anti-horn fly vaccine antigen candidates using a reverse vaccinology approach

Background

The horn fly, Haematobia irritans irritans, causes significant production losses to the cattle industry. Horn fly control relies on insecticides; however, alternative control methods such as vaccines are needed due to the fly's capacity to quickly develop resistance to insecticides, and the pressure for eco-friendly options.

Methods

We used a reverse vaccinology approach comprising three vaccine prediction and 11 annotation tools to evaluate and rank 79,542 translated open reading frames (ORFs) from the horn fly's transcriptome, and selected 10 transcript ORFs as vaccine candidates for expression in Pichia pastoris. The expression of the 10 selected transcripts and the proteins that they encoded were investigated in adult flies by reverse transcription polymerase chain reaction (RT-PCR) and mass spectrometry, respectively. Then, we evaluated the immunogenicity of a vaccine candidate in an immunization trial and the antigen’s effects on horn fly mortality and fecundity in an in vitro feeding assay.

Results

Six of the ten vaccine candidate antigens were successfully expressed in P. pastoris. RT-PCR confirmed the expression of all six ORFs in adult fly RNA. One of the vaccine candidate antigens, BI-HS009, was expressed in sufficient quantity for immunogenicity and efficacy trials. The IgG titers of animals vaccinated with BI-HS009 plus adjuvant were significantly higher than those of animals vaccinated with buffer plus adjuvant only from days 42 to 112, with a peak on day 56. Progeny of horn flies feeding upon blood from animals vaccinated with BI-HS009 plus adjuvant collected on day 56 had 63% lower pupariation rate and 57% lower adult emergence than the control group (ANOVA: F_{(1, 6)} = 8.221, P = 0.028 and F_{(1, 6)} = 8.299, P = 0.028, respectively).

Conclusions

The reverse vaccinology approach streamlined the discovery process by prioritizing possible vaccine antigen candidates. Through a thoughtful process of selection and in vivo and in vitro evaluations, we were able to identify a promising antigen for an anti-horn fly vaccine.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04938-5.

Identification of vaccine targets in pathogens and design of a vaccine using computational approaches

Antigen identification is an important step in the vaccine development process. Computational approaches including deep learning systems can play an important role in the identification of vaccine targets using genomic and proteomic information. Here, we present a new computational system to discover and analyse novel vaccine targets leading to the design of a multi-epitope subunit vaccine candidate. The system incorporates reverse vaccinology and immuno-informatics tools to screen genomic and proteomic datasets of several pathogens such as Trypanosoma cruzi, Plasmodium falciparum, and Vibrio cholerae to identify potential vaccine candidates (PVC). Further, as a case study, we performed a detailed analysis of the genomic and proteomic dataset of T. cruzi (CL Brenner and Y strain) to shortlist eight proteins as possible vaccine antigen candidates using properties such as secretory/surface-exposed nature, low transmembrane helix (< 2), essentiality, virulence, antigenic, and non-homology with host/gut flora proteins. Subsequently, highly antigenic and immunogenic MHC class I, MHC class II and B cell epitopes were extracted from top-ranking vaccine targets. The designed vaccine construct containing 24 epitopes, 3 adjuvants, and 4 linkers was analysed for its physicochemical properties using different tools, including docking analysis. Immunological simulation studies suggested significant levels of T-helper, T-cytotoxic cells, and IgG1 will be elicited upon administration of such a putative multi-epitope vaccine construct. The vaccine construct is predicted to be soluble, stable, non-allergenic, non-toxic, and to offer cross-protection against related Trypanosoma species and strains. Further, studies are required to validate safety and immunogenicity of the vaccine.

Mitochondrial Targeting of the Ammonia-Sensitive Uncoupler SLC4A11 by the Chaperone-Mediated Carrier Pathway in Corneal Endothelium

Purpose

SLC4A11, an electrogenic H⁺ transporter, is found in the plasma membrane and mitochondria of corneal endothelium. However, the underlying mechanism of SLC4A11 targeting to mitochondria is unknown.

Methods

The presence of mitochondrial targeting sequences was examined using in silico mitochondrial proteomic analyses. Thiol crosslinked peptide binding to SLC4A11 was screened by untargeted liquid chromatography/tandem mass spectrometry (LC-MS/MS) analysis. Direct protein interactions between SLC4A11 and chaperones were examined using coimmunoprecipitation analysis and proximity ligation assay. Knockdown or pharmacologic inhibition of chaperones in human corneal endothelial cells (HCECs) or mouse corneal endothelial cells (MCECs), ex vivo kidney, or HA-SLC4A11–transfected fibroblasts was performed to investigate the functional consequences of interfering with mitochondrial SLC4A11 trafficking.

Results

SLC4A11 does not contain canonical N-terminal mitochondrial targeting sequences. LC-MS/MS analysis showed that HSC70 and/or HSP90 are bound to HA-SLC4A11–transfected PS120 fibroblast whole-cell lysates or isolated mitochondria, suggesting trafficking through the chaperone-mediated carrier pathway. SLC4A11 and either HSP90 or HSC70 complexes are directly bound to the mitochondrial surface receptor, TOM70. Interference with this trafficking leads to dysfunctional mitochondrial glutamine catabolism and increased reactive oxygen species production. In addition, glutamine (Gln) use upregulated SLC4A11, HSP70, and HSP90 expression in whole-cell lysates or purified mitochondria of HCECs and HA-SLC4A11–transfected fibroblasts.

Conclusions

HSP90 and HSC70 are critical in mediating mitochondrial SLC4A11 translocation in corneal endothelial cells and kidney. Gln promotes SLC4A11 import to the mitochondria, and the continuous oxidative stress derived from Gln catabolism induced HSP70 and HSP90, protecting cells against oxidative stress.

Calmodulin7: recent insights into emerging roles in plant development and stress

Analyses of the function of Arabidopsis Calmodulin7 (CAM7) in concert with multiple regulatory proteins involved in various signal transduction processes. Calmodulin (CaM) plays various regulatory roles in multiple signaling pathways in eukaryotes. Arabidopsis CALMODULIN 7 (CAM7) is a unique member of the CAM family that works as a transcription factor in light signaling pathways. CAM7 works in concert with CONSTITUTIVE PHOTOMORPHOGENIC 1 and ELONGATED HYPOCOTYL 5, and plays an important role in seedling development. Further, it is involved in the regulation of the activity of various Ca²⁺-gated channels such as cyclic nucleotide gated channel 6 (CNGC6), CNGC14 and auto-inhibited Ca²⁺ ATPase 8. Recent studies further indicate that CAM7 is also an integral part of multiple signaling pathways including hormone, immunity and stress. Here, we review the recent advances in understanding the multifaceted role of CAM7. We highlight the open-ended questions, and also discuss the diverse aspects of CAM7 characterization that need to be addressed for comprehensive understanding of its cellular functions.

Auxin Metabolome Profiling in the Arabidopsis Endoplasmic Reticulum Using an Optimised Organelle Isolation Protocol

The endoplasmic reticulum (ER) is an extensive network of intracellular membranes. Its major functions include proteosynthesis, protein folding, post-transcriptional modification and sorting of proteins within the cell, and lipid anabolism. Moreover, several studies have suggested that it may be involved in regulating intracellular auxin homeostasis in plants by modulating its metabolism. Therefore, to study auxin metabolome in the ER, it is necessary to obtain a highly enriched (ideally, pure) ER fraction. Isolation of the ER is challenging because its biochemical properties are very similar to those of other cellular endomembranes. Most published protocols for ER isolation use density gradient ultracentrifugation, despite its suboptimal resolving power. Here we present an optimised protocol for ER isolation from Arabidopsis thaliana seedlings for the subsequent mass spectrometric determination of ER-specific auxin metabolite profiles. Auxin metabolite analysis revealed highly elevated levels of active auxin form (IAA) within the ER compared to whole plants. Moreover, samples prepared using our optimised isolation ER protocol are amenable to analysis using various “omics” technologies including analyses of both macromolecular and low molecular weight compounds from the same sample.

Molecular Characterization, Expression Analysis of Carotenoid, Xanthophyll, Apocarotenoid Pathway Genes, and Carotenoid and Xanthophyll Accumulation in Chelidonium majus L

Chelidonium majus L. is a perennial herbaceous plant that has various medicinal properties. However, the genomic information about its carotenoid biosynthesis pathway (CBP), xanthophyll biosynthesis pathway (XBP), and apocarotenoid biosynthesis pathway (ABP) genes were limited. Thus, the CBP, XBP, and ABP genes of C. majus were identified and analyzed. Among the 15 carotenoid pathway genes identified, 11 full and 4 partial open reading frames were determined. Phylogenetic analysis of these gene sequences showed higher similarity with higher plants. Through 3D structural analysis and multiple alignments, several distinct conserved motifs were identified, including dinucleotide binding motif, carotene binding motif, and aspartate or glutamate residues. Quantitative RT-PCR showed that CBP, XBP, and ABP genes were expressed in a tissue-specific manner; the highest expression levels were achieved in flowers, followed by those in leaves, roots, and stems. The HPLC analysis of the different organs showed the presence of eight different carotenoids. The highest total carotenoid content was found in leaves, followed by that in flowers, stems, and roots. This study provides information on the molecular mechanisms involved in CBP, XBP, and ABP genes, which might help optimize the carotenoid production in C. majus. The results could also be a basis of further studies on the molecular genetics and functional analysis of CBP, XBP, and ABP genes.

Investigation of cell wall proteins of C. sinensis leaves by combining cell wall proteomics and N-glycoproteomics

BACKGROUND

C. sinensis is an important economic crop with fluoride over-accumulation in its leaves, which poses a serious threat to human health due to its leaf consumption as tea. Recently, our study has indicated that cell wall proteins (CWPs) probably play a vital role in fluoride accumulation/detoxification in C. sinensis. However, there has been a lack in CWP identification and characterization up to now. This study is aimed to characterize cell wall proteome of C. sinensis leaves and to develop more CWPs related to stress response. A strategy of combined cell wall proteomics and N-glycoproteomics was employed to investigate CWPs. CWPs were extracted by sequential salt buffers, while N-glycoproteins were enriched by hydrophilic interaction chromatography method using C. sinensis leaves as a material. Afterwards all the proteins were subjected to UPLC-MS/MS analysis.

RESULTS

A total of 501 CWPs and 195 CWPs were identified respectively by cell wall proteomics and N-glycoproteomics profiling with 118 CWPs in common. Notably, N-glycoproteomics is a feasible method for CWP identification, and it can enhance CWP coverage. Among identified CWPs, proteins acting on cell wall polysaccharides constitute the largest functional class, most of which might be involved in cell wall structure remodeling. The second largest functional class mainly encompass various proteases related to CWP turnover and maturation. Oxidoreductases represent the third largest functional class, most of which (especially Class III peroxidases) participate in defense response. As expected, identified CWPs are mainly related to plant cell wall formation and defense response.

CONCLUSION

This was the first large-scale investigation of CWPs in C. sinensis through cell wall proteomics and N-glycoproteomics. Our results not only provide a database for further research on CWPs, but also an insight into cell wall formation and defense response in C. sinensis.

An essential, kinetoplastid-specific GDP-Fuc: β-D-Gal α-1,2-fucosyltransferase is located in the mitochondrion of Trypanosoma brucei

Fucose is a common component of eukaryotic cell-surface glycoconjugates, generally added by Golgi-resident fucosyltransferases. Whereas fucosylated glycoconjugates are rare in kinetoplastids, the biosynthesis of the nucleotide sugar GDP-Fuc has been shown to be essential in Trypanosoma brucei. Here we show that the single identifiable T. brucei fucosyltransferase (TbFUT1) is a GDP-Fuc: β-D-galactose α-1,2-fucosyltransferase with an apparent preference for a Galβ1,3GlcNAcβ1-O-R acceptor motif. Conditional null mutants of TbFUT1 demonstrated that it is essential for both the mammalian-infective bloodstream form and the insect vector-dwelling procyclic form. Unexpectedly, TbFUT1 was localized in the mitochondrion of T. brucei and found to be required for mitochondrial function in bloodstream form trypanosomes. Finally, the TbFUT1 gene was able to complement a Leishmania major mutant lacking the homologous fucosyltransferase gene (Guo et al., 2021). Together these results suggest that kinetoplastids possess an unusual, conserved and essential mitochondrial fucosyltransferase activity that may have therapeutic potential across trypanosomatids.

Comparative physiological, biochemical and transcriptomic analysis of hexaploid wheat (T. aestivum L.) roots and shoots identifies potential pathways and their molecular regulatory network during Fe and Zn starvation

The combined effect of iron (Fe) and zinc (Zn) starvation on their uptake and transportation and the molecular regulatory networks is poorly understood in wheat. To fill this gap, we performed a comprehensive physiological, biochemical and transcriptome analysis in two bread wheat genotypes, i.e. Narmada 195 and PBW 502, differing in inherent Fe and Zn content. Compared to PBW 502, Narmada 195 exhibited increased tolerance to Fe and Zn withdrawal by significantly modulating the critical physiological and biochemical parameters. We identified 25 core genes associated with four key pathways, i.e. methionine cycle, phytosiderophore biosynthesis, antioxidant and transport system, that exhibited significant up-regulation in both the genotypes with a maximum in Narmada 195. We also identified 26 microRNAs targeting 14 core genes across the four pathways. Together, core genes identified can serve as valuable resources for further functional research for genetic improvement of Fe and Zn content in wheat grain.

Channeling of Carbon Flux Towards Carotenogenesis in Botryococcus braunii: A Media Engineering Perspective

Microalgae, due to their unique properties, gained attention for producing promising feedstocks having high contents of proteins, antioxidants, carotenoids, and terpenoids for applications in nutraceutical and pharmaceutical industries. Optimizing production of the high-value renewables (HVRs) in microalgae requires an in-depth understanding of their functional relationship of the genes involved in these metabolic pathways. In the present study, bioinformatic tools were employed for characterization of the protein-encoding genes of methyl erythritol phosphate (MEP) pathway involved in carotenoid and squalene biosynthesis based upon their conserved motif/domain organization. Our analysis demonstrates nearly 200 putative genes showing a conservation pattern within divergent microalgal lineages. Furthermore, phylogenomic studies confirm the close evolutionary proximity among these microalgal strains in the carotenoid and squalene biosynthetic pathways. Further analysis employing STRING predicts interactions among two rate-limiting genes, i.e., phytoene synthase (PSY) and farnesyl diphosphate farnesyl synthase (FPPS), which are specifically involved in the synthesis of carotenoids and squalene. Experimentally, to understand the carbon flux of these rate-limiting genes involved in carotenogenesis, an industrial potential strain, namely, Botryococcus braunii, was selected in this study for improved biomass productivity (i.e., 100 mg L^-1 D^-1) along with enhanced carotenoid content [0.18% dry cell weight (DCW)] when subjected to carbon supplementation. In conclusion, our approach of media engineering demonstrates that the channeling of carbon flux favors carotenogenesis rather than squalene synthesis. Henceforth, employing omics perspectives will further provide us with new insights for engineering regulatory networks for enhanced production of high-value carbon biorenewables without compromising growth.

A Chromosome-Scale Assembly of the Wheat Leaf Rust Pathogen Puccinia triticina Provides Insights Into Structural Variations and Genetic Relationships With Haplotype Resolution

Despite the global economic importance of the wheat leaf rust pathogen Puccinia triticina (Pt), genomic resources for Pt are limited and chromosome-level assemblies of Pt are lacking. Here, we present a complete haplotype-resolved genome assembly at a chromosome-scale for Pt using the Australian pathotype 64-(6),(7),(10),11 (Pt64; North American race LBBQB) built upon the newly developed technologies of PacBio and Hi-C sequencing. PacBio reads with ∼200-fold coverage (29.8 Gb data) were assembled by Falcon and Falcon-unzip and subsequently scaffolded with Hi-C data using Falcon-phase and Proximo. This approach allowed us to construct 18 chromosome pseudomolecules ranging from 3.5 to 12.3 Mb in size for each haplotype of the dikaryotic genome of Pt64. Each haplotype had a total length of ∼147 Mb, scaffold N 50 of ∼9.4 Mb, and was ∼93% complete for BUSCOs. Each haplotype had ∼29,800 predicted genes, of which ∼2,000 were predicted as secreted proteins (SPs). The investigation of structural variants (SVs) between haplotypes A and B revealed that 10% of the total genome was spanned by SVs, highlighting variations previously undetected by short-read based assemblies. For the first time, the mating type (MAT) genes on each haplotype of Pt64 were identified, which showed that MAT loci a and b are located on two chromosomes (chromosomes 7 and 14), representing a tetrapolar type. Furthermore, the Pt64 assembly enabled haplotype-based evolutionary analyses for 21 Australian Pt isolates, which highlighted the importance of a haplotype resolved reference when inferring genetic relationships using whole genome SNPs. This Pt64 assembly at chromosome-scale with full phase information provides an invaluable resource for genomic and evolutionary research, which will accelerate the understanding of molecular mechanisms underlying Pt-wheat interactions and facilitate the development of durable resistance to leaf rust in wheat and sustainable control of rust disease.

iMLP, a predictor for internal matrix targeting-like sequences in mitochondrial proteins

Matrix targeting sequences (MTSs) direct proteins from the cytosol into mitochondria. Efficient targeting often relies on internal matrix targeting-like sequences (iMTS-Ls) which share structural features with MTSs. Predicting iMTS-Ls was tedious and required multiple tools and webservices. We present iMLP, a deep learning approach for the prediction of iMTS-Ls in protein sequences. A recurrent neural network has been trained to predict iMTS-L propensity profiles for protein sequences of interest. The iMLP predictor considerably exceeds the speed of existing approaches. Expanding on our previous work on iMTS-L prediction, we now serve an intuitive iMLP webservice available at http://iMLP.bio.uni-kl.de and a stand-alone command line tool for power user in addition.

Genome-wide identification of alcohol dehydrogenase (ADH) gene family under waterlogging stress in wheat (Triticum aestivum)

Background

Alcohol dehydrogenase (ADH) plays an important role in plant survival under anaerobic conditions. Although some research about ADH in many plants have been carried out, the bioinformatics analysis of the ADH gene family from Triticum aestivum and their response to abiotic stress is unclear.

Methods

A total of 22 ADH genes were identified from the wheat genome, and these genes could be divided into two subfamilies (subfamily I and subfamily II). All TaADH genes belonged to the Medium-chain ADH subfamily. Sequence alignment analysis showed that all TaADH proteins contained a conservative GroES-like domain and Zinc-binding domain. A total of 64 duplicated gene pairs were found, and the Ka/Ks value of these gene pairs was less than 1, which indicated that these genes were relatively conservative and did not change greatly in the process of duplication.

Results

The organizational analysis showed that nine TaADH genes were highly expressed in all organs, and the rest of TaADH genes had tissue specificity. Cis-acting element analysis showed that almost all of the TaADH genes contained an anaerobic response element. The expression levels of ADH gene in waterlogging tolerant and waterlogging sensitive wheat seeds were analyzed by quantitative real-time PCR (qRT-PCR). This showed that some key ADH genes were significantly responsive to waterlogging stress at the seed germination stage, and the response of waterlogging tolerant and waterlogging sensitive wheat seeds to waterlogging stress was regulated by different ADH genes. The results may be helpful to further study the function of TaADH genes and to determine the candidate gene for wheat stress resistance breeding.

The complete mitochondrial genome of Cycas debaoensis revealed unexpected static evolution in gymnosperm species

Mitochondrial genomes of vascular plants are well known for their liability in architecture evolution. However, the evolutionary features of mitogenomes at intra-generic level are seldom studied in vascular plants, especially among gymnosperms. Here we present the complete mitogenome of Cycas debaoensis, an endemic cycad species to the Guangxi region in southern China. In addition to assemblage of draft mitochondrial genome, we test the conservation of gene content and mitogenomic stability by comparing it to the previously published mitogenome of Cycas taitungensis. Furthermore, we explored the factors such as structural rearrangements and nuclear surveillance of double-strand break repair (DSBR) proteins in Cycas in comparison to other vascular plant groups. The C. debaoensis mitogenome is 413,715 bp in size and encodes 69 unique genes, including 40 protein coding genes, 26 tRNAs, and 3 rRNA genes, similar to that of C. taitungensis. Cycas mitogenomes maintained the ancestral intron content of seed plants (26 introns), which is reduced in other lineages of gymnosperms, such as Ginkgo biloba, Taxus cuspidata and Welwitschia mirabilis due to selective pressure or retroprocessing events. C. debaoensis mitogenome holds 1,569 repeated sequences (> 50 bp), which partially account for fairly large intron size (1200 bp in average) of Cycas mitogenome. The comparison of RNA-editing sites revealed 267 shared non-silent editing site among predicted vs. empirically observed editing events. Another 33 silent editing sites from empirical data increase the total number of editing sites in Cycas debaoensis mitochondrial protein coding genes to 300. Our study revealed unexpected conserved evolution between the two Cycas species. Furthermore, we found strict collinearity of the gene order along with the identical set of genomic content in Cycas mt genomes. The stability of Cycas mt genomes is surprising despite the existence of large number of repeats. This structural stability may be related to the relative expansion of three DSBR protein families (i.e., RecA, OSB, and RecG) in Cycas nuclear genome, which inhibit the homologous recombinations, by monitoring the accuracy of mitochondrial chromosome repair.

Genomics and transcriptomics yields a system-level view of the biology of the pathogen Naegleria fowleri

BACKGROUND

The opportunistic pathogen Naegleria fowleri establishes infection in the human brain, killing almost invariably within 2 weeks. The amoeba performs piece-meal ingestion, or trogocytosis, of brain material causing direct tissue damage and massive inflammation. The cellular basis distinguishing N. fowleri from other Naegleria species, which are all non-pathogenic, is not known. Yet, with the geographic range of N. fowleri advancing, potentially due to climate change, understanding how this pathogen invades and kills is both important and timely.

RESULTS

Here, we report an -omics approach to understanding N. fowleri biology and infection at the system level. We sequenced two new strains of N. fowleri and performed a transcriptomic analysis of low- versus high-pathogenicity N. fowleri cultured in a mouse infection model. Comparative analysis provides an in-depth assessment of encoded protein complement between strains, finding high conservation. Molecular evolutionary analyses of multiple diverse cellular systems demonstrate that the N. fowleri genome encodes a similarly complete cellular repertoire to that found in free-living N. gruberi. From transcriptomics, neither stress responses nor traits conferred from lateral gene transfer are suggested as critical for pathogenicity. By contrast, cellular systems such as proteases, lysosomal machinery, and motility, together with metabolic reprogramming and novel N. fowleri proteins, are all implicated in facilitating pathogenicity within the host. Upregulation in mouse-passaged N. fowleri of genes associated with glutamate metabolism and ammonia transport suggests adaptation to available carbon sources in the central nervous system.

CONCLUSIONS

In-depth analysis of Naegleria genomes and transcriptomes provides a model of cellular systems involved in opportunistic pathogenicity, uncovering new angles to understanding the biology of a rare but highly fatal pathogen.

Evolutionary Implications of a Peroxidase with High Affinity for Cinnamyl Alcohols from Physcomitrium patens, a Non-Vascular Plant

Physcomitrium (Physcomitrella) patens is a bryophyte highly tolerant to different stresses, allowing survival when water supply is a limiting factor. This moss lacks a true vascular system, but it has evolved a primitive water-conducting system that contains lignin-like polyphenols. By means of a three-step protocol, including ammonium sulfate precipitation, adsorption chromatography on phenyl Sepharose and cationic exchange chromatography on SP Sepharose, we were able to purify and further characterize a novel class III peroxidase, PpaPrx19, upregulated upon salt and H₂O₂ treatments. This peroxidase, of a strongly basic nature, shows surprising homology to angiosperm peroxidases related to lignification, despite the lack of true lignins in P. patens cell walls. Moreover, PpaPrx19 shows catalytic and kinetic properties typical of angiosperm peroxidases involved in oxidation of monolignols, being able to efficiently use hydroxycinnamyl alcohols as substrates. Our results pinpoint the presence in P. patens of peroxidases that fulfill the requirements to be involved in the last step of lignin biosynthesis, predating the appearance of true lignin.

EDGAR3.0: comparative genomics and phylogenomics on a scalable infrastructure

The EDGAR platform, a web server providing databases of precomputed orthology data for thousands of microbial genomes, is one of the most established tools in the field of comparative genomics and phylogenomics. Based on precomputed gene alignments, EDGAR allows quick identification of the differential gene content, i.e. the pan genome, the core genome, or singleton genes. Furthermore, EDGAR features a wide range of analyses and visualizations like Venn diagrams, synteny plots, phylogenetic trees, as well as Amino Acid Identity (AAI) and Average Nucleotide Identity (ANI) matrices. During the last few years, the average number of genomes analyzed in an EDGAR project increased by two orders of magnitude. To handle this massive increase, a completely new technical backend infrastructure for the EDGAR platform was designed and launched as EDGAR3.0. For the calculation of new EDGAR3.0 projects, we are now using a scalable Kubernetes cluster running in a cloud environment. A new storage infrastructure was developed using a file-based high-performance storage backend which ensures timely data handling and efficient access. The new data backend guarantees a memory efficient calculation of orthologs, and parallelization has led to drastically reduced processing times. Based on the advanced technical infrastructure new analysis features could be implemented including POCP and FastANI genomes similarity indices, UpSet intersecting set visualization, and circular genome plots. Also the public database section of EDGAR was largely updated and now offers access to 24,317 genomes in 749 free-to-use projects. In summary, EDGAR 3.0 provides a new, scalable infrastructure for comprehensive microbial comparative gene content analysis. The web server is accessible at http://edgar3.computational.bio.

The branched-chain amino acid aminotransferase TaBCAT1 modulates amino acid metabolism and positively regulates wheat rust susceptibility

Plant pathogens suppress defense responses to evade recognition and promote successful colonization. Although identifying the genes essential for pathogen ingress has traditionally relied on screening mutant populations, the post-genomic era provides an opportunity to develop novel approaches that accelerate identification. Here, RNA-seq analysis of 68 pathogen-infected bread wheat (Triticum aestivum) varieties, including three (Oakley, Solstice and Santiago) with variable levels of susceptibility, uncovered a branched-chain amino acid aminotransferase (termed TaBCAT1) as a positive regulator of wheat rust susceptibility. We show that TaBCAT1 is required for yellow and stem rust infection and likely functions in branched-chain amino acid (BCAA) metabolism, as TaBCAT1 disruption mutants had elevated BCAA levels. TaBCAT1 mutants also exhibited increased levels of salicylic acid (SA) and enhanced expression of associated defense genes, indicating that BCAA regulation, via TaBCAT1, has a key role in SA-dependent defense activation. We also identified an association between the levels of BCAAs and resistance to yellow rust infection in wheat. These findings provide insight into SA-mediated defense responses in wheat and highlight the role of BCAA metabolism in the defense response. Furthermore, TaBCAT1 could be manipulated to potentially provide resistance to two of the most economically damaging diseases of wheat worldwide.

Host Gamma-Aminobutyric Acid Metabolic Pathway Is Involved in Resistance Against Rhizoctonia solani

Rhizoctonia solani is a highly destructive necrotrophic fungal pathogen having a diverse host range, including rice and tomato. Previously R. solani infection has been found to cause large-scale readjustment in host primary metabolism and accumulation of various stress-associated metabolites such as gamma-aminobutyric acid (GABA) in rice. In this study, we report upregulation of GABA pathway genes during pathogenesis of R. solani in rice and tomato. The exogenous application of GABA provided partial resistance against R. solani infection in both the hosts. Furthermore, by using the virus-induced gene silencing approach, we knocked down the expression of some of the tomato genes involved in GABA biosynthesis (glutamate decarboxylase) and GABA catabolism (GABA-transaminase and succinic semialdehyde dehydrogenase) to study their role in host defense against R. solani infection. The silencing of each of these genes increased disease susceptibility in tomato. Overall the results from gene expression analysis, exogenous chemical application, and gene silencing studies suggest that the GABA pathway plays a positive role in plant defense against necrotrophic pathogen R. solani.

Genome analysis of Candida subhashii reveals its hybrid nature and dual mitochondrial genome conformations

Candida subhashii belongs to the CUG-Ser clade, a group of phylogenetically closely related yeast species that includes some human opportunistic pathogens, such as Candida albicans. Despite being present in the environment, C. subhashii was initially described as the causative agent of a case of peritonitis. Considering the relevance of whole-genome sequencing and analysis for our understanding of genome evolution and pathogenicity, we sequenced, assembled and annotated the genome of C. subhashii type strain. Our results show that C. subhashii presents a highly heterozygous genome and other signatures that point to a hybrid ancestry. The presence of functional pathways for assimilation of hydroxyaromatic compounds goes in line with the affiliation of this yeast with soil microbial communities involved in lignin decomposition. Furthermore, we observed that different clones of this strain may present circular or linear mitochondrial DNA. Re-sequencing and comparison of strains with differential mitochondrial genome topology revealed five candidate genes potentially associated with this conformational change: MSK1, SSZ1, ALG5, MRPL9 and OYE32.

YbeY is required for ribosome small subunit assembly and tRNA processing in human mitochondria

Mitochondria contain their own translation apparatus which enables them to produce the polypeptides encoded in their genome. The mitochondrially-encoded RNA components of the mitochondrial ribosome require various post-transcriptional processing steps. Additional protein factors are required to facilitate the biogenesis of the functional mitoribosome. We have characterized a mitochondrially-localized protein, YbeY, which interacts with the assembling mitoribosome through the small subunit. Loss of YbeY leads to a severe reduction in mitochondrial translation and a loss of cell viability, associated with less accurate mitochondrial tRNA^Ser(AGY) processing from the primary transcript and a defect in the maturation of the mitoribosomal small subunit. Our results suggest that YbeY performs a dual, likely independent, function in mitochondria being involved in precursor RNA processing and mitoribosome biogenesis. Issue Section: Nucleic Acid Enzymes.

Duplication and Functional Divergence of Branched-Chain Amino Acid Biosynthesis Genes in Aspergillus nidulans

Fungi, bacteria, and plants, but not animals, synthesize the branched-chain amino acids: leucine, isoleucine, and valine. While branched-chain amino acid (BCAA) biosynthesis has been well characterized in the yeast Saccharomyces cerevisiae, it is incompletely understood in filamentous fungi. The three BCAAs share several early biosynthesis steps before divergence into specific pathways. In Aspergillus nidulans, the genes for the first two dedicated steps in leucine biosynthesis have been characterized, but the final two have not. We used sequence searches of the A. nidulans genome to identify two genes encoding β-isopropylmalate dehydrogenase, which catalyzes the penultimate step of leucine biosynthesis, and six genes encoding BCAA aminotransferase, which catalyzes the final step in biosynthesis of all three BCAA. We have used combinations of gene knockouts to determine the relative contribution of each of these genes to BCAA biosynthesis. While both β-isopropylmalate dehydrogenase genes act in leucine biosynthesis, the two most highly expressed BCAA aminotransferases are responsible for BCAA biosynthesis. We have also characterized the expression of leucine biosynthesis genes using reverse transcriptase-quantitative PCR and found regulation in response to leucine availability is mediated through the Zn(II)₂Cys₆ transcription factor LeuB. IMPORTANCE Branched-chain amino acid (BCAA) biosynthesis is important for pathogenic fungi to successfully cause disease in human and plant hosts. The enzymes for their production are absent from humans and, therefore, provide potential antifungal targets. While BCAA biosynthesis is well characterized in yeasts, it is poorly understood in filamentous fungal pathogens. Developing a thorough understanding of both the genes encoding the metabolic enzymes for BCAA biosynthesis and how their expression is regulated will inform target selection for antifungal drug development.

Heterotrophic euglenid Rhabdomonas costata resembles its phototrophic relatives in many aspects of molecular and cell biology

Euglenids represent a group of protists with diverse modes of feeding. To date, only a partial genomic sequence of Euglena gracilis and transcriptomes of several phototrophic and secondarily osmotrophic species are available, while primarily heterotrophic euglenids are seriously undersampled. In this work, we begin to fill this gap by presenting genomic and transcriptomic drafts of a primary osmotroph, Rhabdomonas costata. The current genomic assembly length of 100 Mbp is 14× smaller than that of E. gracilis. Despite being too fragmented for comprehensive gene prediction it provided fragments of the mitochondrial genome and comparison of the transcriptomic and genomic data revealed features of its introns, including several candidates for nonconventional types. A set of 39,456 putative R. costata proteins was predicted from the transcriptome. Annotation of the mitochondrial core metabolism provides the first data on the facultatively anaerobic mitochondrion of R. costata, which in most respects resembles the mitochondrion of E. gracilis with a certain level of streamlining. R. costata can synthetise thiamine by enzymes of heterogenous provenances and haem by a mitochondrial-cytoplasmic C4 pathway with enzymes orthologous to those found in E. gracilis. The low percentage of green algae-affiliated genes supports the ancestrally osmotrophic status of this species.

Genomic Signatures of Coevolution between Nonmodel Mammals and Parasitic Roundworms

Antagonistic coevolution between host and parasite drives species evolution. However, most of the studies only focus on parasitism adaptation and do not explore the coevolution mechanisms from the perspective of both host and parasite. Here, through the de novo sequencing and assembly of the genomes of giant panda roundworm, red panda roundworm, and lion roundworm parasitic on tiger, we investigated the genomic mechanisms of coevolution between nonmodel mammals and their parasitic roundworms and those of roundworm parasitism in general. The genome-wide phylogeny revealed that these parasitic roundworms have not phylogenetically coevolved with their hosts. The CTSZ and prolyl 4-hydroxylase subunit beta (P4HB) immunoregulatory proteins played a central role in protein interaction between mammals and parasitic roundworms. The gene tree comparison identified that seven pairs of interactive proteins had consistent phylogenetic topology, suggesting their coevolution during host–parasite interaction. These coevolutionary proteins were particularly relevant to immune response. In addition, we found that the roundworms of both pandas exhibited higher proportions of metallopeptidase genes, and some positively selected genes were highly related to their larvae’s fast development. Our findings provide novel insights into the genetic mechanisms of coevolution between nonmodel mammals and parasites and offer the valuable genomic resources for scientific ascariasis prevention in both pandas.

Cloning and Expression of the Neuropeptide F and Neuropeptide F Receptor Genes and Their Regulation of Food Intake in the Chinese White Pine Beetle Dendroctonus armandi

Neuropeptide F (NPF) is an important signaling molecule that acts as a neuromodulator to regulate a diversity of physiological and behavioral processes from vertebrates to invertebrates by interaction with NPF receptors, which are G protein-coupled receptors (GPCR). However, nothing is known about NPF in Chinese white pine beetle, Dendroctonus armandi, a destructive pest of natural and coniferous forests in the middle Qinling Mountains of China. We have cloned and characterized cDNAs encoding one NPF precursor and two NPF receptors in D. armandi and made bioinformatics predictions according to the deduced amino acid sequences. They were highly similar to that of Dendroctonus ponderosa. The transcription levels of these genes were different between larvae and adults of sexes, and there were significant differences among the different developmental stages and tissues and between beetles under starvation and following re-feeding states. Additionally, downregulation of NPF and NPFR by injecting dsRNA into beetles reduced their food intake, caused increases of mortality and decreases of body weight, and also resulted in a decrease of glycogen and free fatty acid and an increase of trehalose. These results indicate that the NPF signaling pathway plays a significant positive role in the regulation of food intake and provides a potential target for the sustainable management of this pest.

Genome-Wide Characterization and Expression Analysis of CAMTA Gene Family Under Salt Stress in Cucurbita moschata and Cucurbita maxima

Cucurbita Linn. vegetables have a long history of cultivation and have been cultivated all over the world. With the increasing area of saline–alkali soil, Cucurbita Linn. is affected by salt stress, and calmodulin-binding transcription activator (CAMTA) is known for its important biological functions. Although the CAMTA gene family has been identified in several species, there is no comprehensive analysis on Cucurbita species. In this study, we analyzed the genome of Cucurbita maxima and Cucurbita moschata. Five C. moschata calmodulin-binding transcription activators (CmoCAMTAs) and six C. maxima calmodulin-binding transcription activators (CmaCAMTAs) were identified, and they were divided into three subfamilies (Subfamilies I, II, and III) based on the sequence identity of amino acids. CAMTAs from the same subfamily usually have similar exon–intron distribution and conserved domains (CG-1, TIG, IQ, and Ank_2). Chromosome localization analysis showed that CmoCAMTAs and CmaCAMTAs were unevenly distributed across four and five out of 21 chromosomes, respectively. There were a total of three duplicate gene pairs, and all of which had experienced segmental duplication events. The transcriptional profiles of CmoCAMTAs and CmaCAMTAs in roots, stems, leaves, and fruits showed that these CAMTAs have tissue specificity. Cis-acting elements analysis showed that most of CmoCAMTAs and CmaCAMTAs responded to salt stress. By analyzing the transcriptional profiles of CmoCAMTAs and CmaCAMTAs under salt stress, it was shown that both C. moschata and C. maxima shared similarities against salt tolerance and that it is likely to contribute to the development of these species. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) further demonstrated the key role of CmoCAMTAs and CmaCAMTAs under salt stress. This study provided a theoretical basis for studying the function and mechanism of CAMTAs in Cucurbita Linn.

Using Diatom and Apicomplexan Models to Study the Heme Pathway of Chromera velia

Heme biosynthesis is essential for almost all living organisms. Despite its conserved function, the pathway’s enzymes can be located in a remarkable diversity of cellular compartments in different organisms. This location does not always reflect their evolutionary origins, as might be expected from the history of their acquisition through endosymbiosis. Instead, the final subcellular localization of the enzyme reflects multiple factors, including evolutionary origin, demand for the product, availability of the substrate, and mechanism of pathway regulation. The biosynthesis of heme in the apicomonad Chromera velia follows a chimeric pathway combining heme elements from the ancient algal symbiont and the host. Computational analyses using different algorithms predict complex targeting patterns, placing enzymes in the mitochondrion, plastid, endoplasmic reticulum, or the cytoplasm. We employed heterologous reporter gene expression in the apicomplexan parasite Toxoplasma gondii and the diatom Phaeodactylum tricornutum to experimentally test these predictions. 5-aminolevulinate synthase was located in the mitochondria in both transfection systems. In T. gondii, the two 5-aminolevulinate dehydratases were located in the cytosol, uroporphyrinogen synthase in the mitochondrion, and the two ferrochelatases in the plastid. In P. tricornutum, all remaining enzymes, from ALA-dehydratase to ferrochelatase, were placed either in the endoplasmic reticulum or in the periplastidial space.

Comparative Genomics of Eight Fusarium graminearum Strains with Contrasting Aggressiveness Reveals an Expanded Open Pangenome and Extended Effector Content Signatures

Fusarium graminearum, the primary cause of Fusarium head blight (FHB) in small-grain cereals, demonstrates remarkably variable levels of aggressiveness in its host, producing different infection dynamics and contrasted symptom severity. While the secreted proteins, including effectors, are thought to be one of the essential components of aggressiveness, our knowledge of the intra-species genomic diversity of F. graminearum is still limited. In this work, we sequenced eight European F. graminearum strains of contrasting aggressiveness to characterize their respective genome structure, their gene content and to delineate their specificities. By combining the available sequences of 12 other F. graminearum strains, we outlined a reference pangenome that expands the repertoire of the known genes in the reference PH-1 genome by 32%, including nearly 21,000 non-redundant sequences and gathering a common base of 9250 conserved core-genes. More than 1000 genes with high non-synonymous mutation rates may be under diverse selection, especially regarding the trichothecene biosynthesis gene cluster. About 900 secreted protein clusters (SPCs) have been described. Mostly localized in the fast sub-genome of F. graminearum supposed to evolve rapidly to promote adaptation and rapid responses to the host's infection, these SPCs gather a range of putative proteinaceous effectors systematically found in the core secretome, with the chloroplast and the plant nucleus as the main predicted targets in the host cell. This work describes new knowledge on the intra-species diversity in F. graminearum and emphasizes putative determinants of aggressiveness, providing a wealth of new candidate genes potentially involved in the Fusarium head blight disease.

Historical genomics reveals the evolutionary mechanisms behind multiple outbreaks of the host-specific coffee wilt pathogen Fusarium xylarioides

BACKGROUND

Nearly 50% of crop yields are lost to pests and disease, with plants and pathogens locked in an amplified co-evolutionary process of disease outbreaks. Coffee wilt disease, caused by Fusarium xylarioides, decimated coffee production in west and central Africa following its initial outbreak in the 1920s. After successful management, it later re-emerged and by the 2000s comprised two separate epidemics on arabica coffee in Ethiopia and robusta coffee in east and central Africa.

RESULTS

Here, we use genome sequencing of six historical culture collection strains spanning 52 years to identify the evolutionary processes behind these repeated outbreaks. Phylogenomic reconstruction using 13,782 single copy orthologs shows that the robusta population arose from the initial outbreak, whilst the arabica population is a divergent sister clade to the other strains. A screen for putative effector genes involved in pathogenesis shows that the populations have diverged in gene content and sequence mainly by vertical processes within lineages. However, 15 putative effector genes show evidence of horizontal acquisition, with close homology to genes from F. oxysporum. Most occupy small regions of homology within wider scaffolds, whereas a cluster of four genes occupy a 20Kb scaffold with strong homology to a region on a mobile pathogenicity chromosome in F. oxysporum that houses known effector genes. Lacking a match to the whole mobile chromosome, we nonetheless found close associations with DNA transposons, especially the miniature impala type previously proposed to facilitate horizontal transfer of pathogenicity genes in F. oxysporum. These findings support a working hypothesis that the arabica and robusta populations partly acquired distinct effector genes via transposition-mediated horizontal transfer from F. oxysporum, which shares coffee as a host and lives on other plants intercropped with coffee.

CONCLUSION

Our results show how historical genomics can help reveal mechanisms that allow fungal pathogens to keep pace with our efforts to resist them. Our list of putative effector genes identifies possible future targets for fungal control. In turn, knowledge of horizontal transfer mechanisms and putative donor taxa might help to design future intercropping strategies that minimize the risk of transfer of effector genes between closely-related Fusarium taxa.

Red light-regulated interaction of Per-Arnt-Sim histidine kinases with partner histidine-containing phosphotransfer proteins in Physcomitrium patens

Multi-step phosphorelay (MSP) is a broadly distributed signaling system in organisms. In MSP, histidine kinases (HKs) receive various environmental signals and transmit them by autophosphorylation followed by phosphotransfer to partner histidine-containing phosphotransfer proteins (HPts). Previously, we reported that Per-Arnt-Sim (PAS) domain-containing HK1 (PHK1) and PHK2 of the moss Physcomitrium (Physcomitrella) patens repressed red light-induced protonema branching, a critical step in the moss life cycle. In plants, PHK homolog-encoding genes are conserved only in early-diverging lineages such as bryophytes and lycophytes. PHKs-mediated signaling machineries attract attention especially from an evolutionary viewpoint, but they remain uninvestigated. Here, we studied the P. patens PHKs focusing on their subcellular patterns of localization and interaction with HPts. Yeast two-hybrid analysis, a localization assay with a green fluorescent protein, and a bimolecular fluorescence complementation analysis together showed that PHKs are localized and interact with partner HPts mostly in the nucleus, as unprecedented features for plant HKs. Additionally, red light triggered the interactions between PHKs and HPts in the cytoplasm, and light co-repressed the expression of PHK1 and PHK2 as well as genes encoding their partner HPts. Our results emphasize the uniqueness of PHKs-mediated signaling machineries, and functional implications of this uniqueness are discussed.

Zone of Interaction Between the Parasite and the Host: Protein Profile of the Body Cavity Fluid of Gasterosteus aculeatus L. Infected with the Cestode Schistocephalus solidus (Muller, 1776)

PURPOSE

During infection, the host and the parasite "communicate" with each other through various molecules, including proteins. The aim of this study was to describe the excretory-secretory proteins from the helminth Schistocephalus solidus and its intermediate host, the three-spined stickleback Gasterosteus aculeatus L., which are likely to be involved in interactions between them.

METHODS

Combined samples of washes from the G. aculeatus sticklebacks cavity infected with the S. solidus, and washes from the parasite surface were used as experimental samples, while washes from the uninfected fish body cavity were used as control. The obtained samples were analyzed using mass-spectrometry nLC-MS/MS.

RESULTS

As a result of mass-spectrometry analysis 215 proteins were identified. Comparative quantitative analysis revealed significant differences in LFQ intensity between experimental and control samples for 20 stickleback proteins. In the experimental samples, we found an increase in the content of serpins, plasminogen, angiotensin 1-10, complement component C9, and a decrease in the content of triosephosphate isomerase, creatine kinase, fructose-biphosphate aldolase, superoxide dismutase, peroxidoxin-1, homocysteine-binding and fatty acid-binding proteins, compared to uninfected fish samples. In the experimental group washes, 30 S. solidus proteins were found, including malate dehydrogenase, annexin family proteins, serpins, peptidyl-prolyl cis-trans isomerase and fatty acid-binding protein.

CONCLUSIONS

Thus, the protein composition of washes from the helminth S. solidus surface and the body cavity of infected and uninfected stickleback G. aculeatus were studied. As a result, it was shown that various components of the immune defense system predominated in the washes of infected fish and helminths.

Expression and analysis of the SAM-dependent RNA methyltransferase Rsm22 from Saccharomyces cerevisiae

Rsm22-family proteins are conserved putative SAM-dependent methyltransferases with important functions in mitochondrial translation. Here, the results of a comparative bioinformatics analysis of Rsm22-type proteins are presented, the expression, biophysical characterization and crystallization of Saccharomyces cerevisiae Rsm22 are reported, a low-resolution SAXS structure of the protein is revealed, and SAM-dependent RNA methyl transferase activity of the protein is demonstrated.

The Saccharomyces cerevisiae Rsm22 protein (Sc-Rsm22), encoded by the nuclear RSM22 (systematic name YKL155c) gene, is a distant homologue of Rsm22 from Trypanosoma brucei (Tb-Rsm22) and METTL17 from mouse (Mm-METTL17). All three proteins have been shown to be associated with mitochondrial gene expression, and Sc-Rsm22 has been documented to be essential for mitochondrial respiration. The Sc-Rsm22 protein comprises a polypeptide of molecular weight 72.2 kDa that is predicted to harbor an N-terminal mitochondrial targeting sequence. The precise physiological function of Rsm22-family proteins is unknown, and no structural information has been available for Sc-Rsm22 to date. In this study, Sc-Rsm22 was expressed and purified in monomeric and dimeric forms, their folding was confirmed by circular-dichroism analyses and their low-resolution structures were determined using a small-angle X-ray scattering (SAXS) approach. The solution structure of the monomeric form of Sc-Rsm22 revealed an elongated three-domain arrangement, which differs from the shape of Tb-Rsm22 in its complex with the mitochondrial small ribosomal subunit in T. brucei (PDB entry 6sg9). A bioinformatic analysis revealed that the core domain in the middle (Leu117–Asp462 in Sc-Rsm22) resembles the corresponding region in Tb-Rsm22, including a Rossmann-like methyltransferase fold followed by a zinc-finger-like structure. The latter structure is not present in this position in other methyltransferases and is therefore a unique structural motif for this family. The first half of the C-terminal domain is likely to form an OB-fold, which is typically found in RNA-binding proteins and is also seen in the Tb-Rsm22 structure. In contrast, the N-terminal domain of Sc-Rsm22 is predicted to be fully α-helical and shares no sequence similarity with other family members. Functional studies demonstrated that the monomeric variant of Sc-Rsm22 methylates mitochondrial tRNAs in vitro. These data suggest that Sc-Rsm22 is a new and unique member of the RNA methyltransferases that is important for mitochondrial protein synthesis.

Phylogeny, tissue-specific expression, and activities of root-secreted purple acid phosphatases for P uptake from ATP in P starved poplar

Plants secrete purple acid phosphatases (PAPs) under phosphorus (P) shortage but the contribution of plant PAPs to P acquisition is not well understood. The goals of this study were to investigate comprehensively the transcription patterns of PAPs under P shortage in poplar (Populus × canescens), to identify secreted PAPs and to characterize their contribution to mobilize organic P. Phylogenetic analyses of the PAP family revealed 33 putative members. In this study, distinct, tissue-specific P responsive expression patterns could be shown for 23 PAPs in roots and leaves. Root-associated PAP activities were localized on the root surface by in-vivo staining. The activities of root-surface PAPs increased significantly under low P availability, but were suppressed by a PAP inhibitor and corresponded to elevated P uptake from ATP as an organic P source. By proteomic analyses of the root apoplast, we identified three newly secreted proteins under P shortage: PtPAP1 (Potri.005G233400) and two proteins with unknown functions (Potri.013G100800 and Potri.001G209300). Our results, based on the combination of transcriptome and proteome analyses with phosphatase activity assays, support that PtPAP1 plays a central role in enhanced P acquisition from organic sources, when the phosphate concentrations in soil are limited.

In Silico Analysis of Functionalized Hydrocarbon Production Using Ehrlich Pathway and Fatty Acid Derivatives in an Endophytic Fungus

Functionalized hydrocarbons have various ecological and industrial uses, from signaling molecules and antifungal/antibacterial agents to fuels and specialty chemicals. The potential to produce functionalized hydrocarbons using the cellulolytic, endophytic fungus, Ascocoryne sarcoides, was quantified using genome-enabled, stoichiometric modeling. In silico analysis identified available routes to produce these hydrocarbons, including both anabolic- and catabolic-associated strategies, and determined correlations between the type and size of the hydrocarbons and culturing conditions. The analysis quantified the limits of the wild-type metabolic network to produce functionalized hydrocarbons from cellulose-based substrates and identified metabolic engineering targets, including cellobiose phosphorylase (CP) and cytosolic pyruvate dehydrogenase complex (PDHcyt). CP and PDHcyt activity increased the theoretical production limits under anoxic conditions where less energy was extracted from the substrate. The incorporation of both engineering targets resulted in near-complete conservation of substrate electrons in functionalized hydrocarbons. The in silico framework was integrated with in vitro fungal batch growth experiments to support O₂ limitation and functionalized hydrocarbon production predictions. The metabolic reconstruction of this endophytic filamentous fungus describes pathways for both specific and general production strategies of 161 functionalized hydrocarbons applicable to many eukaryotic hosts.

StPIP1, a PAMP-induced peptide in potato, elicits plant defenses and is associated with disease symptom severity in a compatible interaction with Potato virus Y

The role of small secreted peptides in plant defense responses to viruses has seldom been investigated. Here, we report a role for potato (Solanum tuberosum) PIP1, a gene predicted to encode a member of the pathogen-associated molecular pattern (PAMP)-induced peptide (PIP) family, in the response of potato to Potato virus Y (PVY) infection. We show that exogenous application of synthetic StPIP1 to potato leaves and nodes increased the production of reactive oxygen species and the expression of plant defense-related genes, revealing that StPIP1 triggers early defense responses. In support of this hypothesis, transgenic potato plants that constitutively overexpress StPIP1 had higher levels of leaf callose deposition and, based on measurements of viral RNA titers, were less susceptible to infection by a compatible PVY strain. Interestingly, systemic infection of StPIP1-overexpressing lines with PVY resulted in clear rugose mosaic symptoms that were absent or very mild in infected non-transgenic plants. A transcriptomics analysis revealed that marker genes associated with both pattern-triggered immunity and effector-triggered immunity were induced in infected StPIP1 overexpressors but not in non-transgenic plants. Together, our results reveal a role for StPIP1 in eliciting plant defense responses and in regulating plant antiviral immunity.

RUS6, a DUF647-containing protein, is essential for early embryonic development in Arabidopsis thaliana

BACKGROUND

The Arabidopsis RUS (ROOT UV-B SENSITIVE) gene family contains six members, each of which encodes a protein containing a DUF647 (domain of unknown function 647) that is commonly found in eukaryotes. Previous studies have demonstrated that RUS1 and RUS2 play critical roles in early seedling development. All six RUS genes are expressed throughout the plant, but little is known about the functional roles of RUS3, RUS4, RUS5 and RUS6.

RESULTS

We used a reverse-genetic approach to identify knockout mutants for RUS3, RUS4, RUS5 and RUS6. Each mutant was confirmed by direct DNA sequencing and genetic segregation analysis. No visible phenotypic differences were observed in rus3, rus4, or rus5 knockout mutants under standard growth conditions, but rus6 knockout mutants displayed a strong embryo-lethal phenotype. Two independent knockout lines for RUS6 were characterized. The rus6 mutations could only be maintained through a heterozygote, because rus6 homozygous mutants did not survive. Closer examinations of homozygous rus6 embryos from rus6/ + parent plants revealed that RUS6 is required for early embryo development. Loss of RUS6 resulted in embryo lethality, specifically at the mid-globular stage. The embryo-lethality phenotype was complemented by a RUS6::RUS6-GFP transgene, and GFP signal was detected throughout the embryo. Histological analyses with the β-glucuronidase reporter gene driven by the RUS6 promoter showed tissue- and development-specific expression of RUS6, which was highest in floral tissues.

CONCLUSION

Our data revealed that RUS6 is essential for early embryo development in Arabidopsis, and that the RUS gene family functions in multiple stages of plant development.

In silico identification of excretory/secretory proteins and drug targets in monogenean parasites

The Excretory/Secretory (ES) proteins of parasites are involved in invasion and colonization of their hosts. In addition, since ES proteins circulate in the extracellular space, they can be more accessible to drugs than other proteins, which makes ES proteins optimal targets for the development of new and better pharmacological strategies. Monogeneans are a group of parasitic Platyhelminthes that includes some pathogenic species problematic for finfish aquaculture. In the present study, 8297 putative ES proteins from four monogenean species which genomic resources are publicly available were identified and functionally annotated by bioinformatic tools. Additionally, for comparative purposes, ES proteins in other parasitic and free-living platyhelminths were identified. Based on data from the monogenean Gyrodactylus salaris, 15 ES proteins are considered potential drug targets. One of them showed homology to 10 cathepsins with known 3D structure. A docking molecular analysis uncovered that the anthelmintic emodepside shows good affinity to these cathepsins suggesting that emodepside can be experimentally tested as a monogenean's cathepsin inhibitor.

Analysis of diverse eukaryotes suggests the existence of an ancestral mitochondrial apparatus derived from the bacterial type II secretion system

The type 2 secretion system (T2SS) is present in some Gram-negative eubacteria and used to secrete proteins across the outer membrane. Here we report that certain representative heteroloboseans, jakobids, malawimonads and hemimastigotes unexpectedly possess homologues of core T2SS components. We show that at least some of them are present in mitochondria, and their behaviour in biochemical assays is consistent with the presence of a mitochondrial T2SS-derived system (miT2SS). We additionally identified 23 protein families co-occurring with miT2SS in eukaryotes. Seven of these proteins could be directly linked to the core miT2SS by functional data and/or sequence features, whereas others may represent different parts of a broader functional pathway, possibly also involving the peroxisome. Its distribution in eukaryotes and phylogenetic evidence together indicate that the miT2SS-centred pathway is an ancestral eukaryotic trait. Our findings thus have direct implications for the functional properties of the early mitochondrion.

Bacteria use the type 2 secretion system to secrete enzymes and toxins across the outer membrane to the environment. Here the authors analyse the T2SS pathway in three protist lineages and suggest that the early mitochondrion may have been capable of secreting proteins into the cytosol.

Allelopathy in rice: a story of momilactones, kin recognition, and weed management

In the struggle to secure nutrient access and to outperform competitors, some plant species have evolved a biochemical arsenal with which they inhibit the growth or development of neighbouring plants. This process, known as allelopathy, exists in many of today's major crops, including rice. Rice synthesizes momilactones, diterpenoids that are released into the rhizosphere and inhibit the growth of numerous plant species. While the allelopathic potential of rice was recognized decades ago, many questions remain unresolved regarding the biosynthesis, exudation, and biological activity of momilactones. Here, we review current knowledge on momilactones, their role in allelopathy, and their potential to serve as a basis for sustainable weed management. We emphasize the gaps in our current understanding of when and how momilactones are produced and of how they act in plant cells, and outline what we consider the next steps in momilactone and rice allelopathy research.

Structural and mechanistic insights into the bifunctional HISN2 enzyme catalyzing the second and third steps of histidine biosynthesis in plants

The second and third steps of the histidine biosynthetic pathway (HBP) in plants are catalyzed by a bifunctional enzyme–HISN2. The enzyme consists of two distinct domains, active respectively as a phosphoribosyl-AMP cyclohydrolase (PRA-CH) and phosphoribosyl-ATP pyrophosphatase (PRA-PH). The domains are analogous to single-domain enzymes encoded by bacterial hisI and hisE genes, respectively. The calculated sequence similarity networks between HISN2 analogs from prokaryotes and eukaryotes suggest that the plant enzymes are closest relatives of those in the class of Deltaproteobacteria. In this work, we obtained crystal structures of HISN2 enzyme from Medicago truncatula (MtHISN2) and described its architecture and interactions with AMP. The AMP molecule bound to the PRA-PH domain shows positioning of the N1-phosphoribosyl relevant to catalysis. AMP bound to the PRA-CH domain mimics a part of the substrate, giving insights into the reaction mechanism. The latter interaction also arises as a possible second-tier regulatory mechanism of the HBP flux, as indicated by inhibition assays and isothermal titration calorimetry.

Protein expression plasticity contributes to heat and drought tolerance of date palm

Climate change is increasing the frequency and intensity of warming and drought periods around the globe, currently representing a threat to many plant species. Understanding the resistance and resilience of plants to climate change is, therefore, urgently needed. As date palm (Phoenix dactylifera) evolved adaptation mechanisms to a xeric environment and can tolerate large diurnal and seasonal temperature fluctuations, we studied the protein expression changes in leaves, volatile organic compound emissions, and photosynthesis in response to variable growth temperatures and soil water deprivation. Plants were grown under controlled environmental conditions of simulated Saudi Arabian summer and winter climates challenged with drought stress. We show that date palm is able to counteract the harsh conditions of the Arabian Peninsula by adjusting the abundances of proteins related to the photosynthetic machinery, abiotic stress and secondary metabolism. Under summer climate and water deprivation, these adjustments included efficient protein expression response mediated by heat shock proteins and the antioxidant system to counteract reactive oxygen species formation. Proteins related to secondary metabolism were downregulated, except for the P. dactylifera isoprene synthase (PdIspS), which was strongly upregulated in response to summer climate and drought. This study reports, for the first time, the identification and functional characterization of the gene encoding for PdIspS, allowing future analysis of isoprene functions in date palm under extreme environments. Overall, the current study shows that reprogramming of the leaf protein profiles confers the date palm heat- and drought tolerance. We conclude that the protein plasticity of date palm is an important mechanism of molecular adaptation to environmental fluctuations.

Genome-Wide Transcriptomic Identification and Functional Insight of Lily WRKY Genes Responding to Botrytis Fungal Disease

Genome-wide transcriptome analysis using RNA-Seq of Lilium longiflorum revealed valuable genes responding to biotic stresses. WRKY transcription factors are regulatory proteins playing essential roles in defense processes under environmental stresses, causing considerable losses in flower quality and production. Thirty-eight WRKY genes were identified from the transcriptomic profile from lily genotypes, exhibiting leaf blight caused by Botrytis elliptica. Lily WRKYs have a highly conserved motif, WRKYGQK, with a common variant, WRKYGKK. Phylogeny of LlWRKYs with homologous genes from other representative plant species classified them into three groups- I, II, and III consisting of seven, 22, and nine genes, respectively. Base on functional annotation, 22 LlWRKY genes were associated with biotic stress, nine with abiotic stress, and seven with others. Sixteen unique LlWRKY were studied to investigate responses to stress conditions using gene expression under biotic and abiotic stress treatments. Five genes-LlWRKY3, LlWRKY4, LlWRKY5, LlWRKY10, and LlWRKY12-were substantially upregulated, proving to be biotic stress-responsive genes in vivo and in vitro conditions. Moreover, the expression patterns of LlWRKY genes varied in response to drought, heat, cold, and different developmental stages or tissues. Overall, our study provides structural and molecular insights into LlWRKY genes for use in the genetic engineering in Lilium against Botrytis disease.

ZmOrphan94 Transcription Factor Downregulates ZmPEPC1 Gene Expression in Maize Bundle Sheath Cells

Spatial separation of the photosynthetic reactions is a key feature of C4 metabolism. In most C4 plants, this separation requires compartmentation of photosynthetic enzymes between mesophyll (M) and bundle sheath (BS) cells. The upstream region of the gene encoding the maize PHOSPHOENOLPYRUVATE CARBOXYLASE 1 (ZmPEPC1) has been shown sufficient to drive M-specific ZmPEPC1 gene expression. Although this region has been well characterized, to date, only few trans-factors involved in the ZmPEPC1 gene regulation were identified. Here, using a yeast one-hybrid approach, we have identified three novel maize transcription factors ZmHB87, ZmCPP8, and ZmOrphan94 as binding to the ZmPEPC1 upstream region. Bimolecular fluorescence complementation assays in maize M protoplasts unveiled that ZmOrphan94 forms homodimers and interacts with ZmCPP8 and with two other ZmPEPC1 regulators previously reported, ZmbHLH80 and ZmbHLH90. Trans-activation assays in maize M protoplasts unveiled that ZmHB87 does not have a clear transcriptional activity, whereas ZmCPP8 and ZmOrphan94 act as activator and repressor, respectively. Moreover, we observed that ZmOrphan94 reduces the trans-activation activity of both activators ZmCPP8 and ZmbHLH90. Using the electromobility shift assay, we showed that ZmOrphan94 binds to several cis-elements present in the ZmPEPC1 upstream region and one of these cis-elements overlaps with the ZmbHLH90 binding site. Gene expression analysis revealed that ZmOrphan94 is preferentially expressed in the BS cells, suggesting that ZmOrphan94 is part of a transcriptional regulatory network downregulating ZmPEPC1 transcript level in the BS cells. Based on both this and our previous work, we propose a model underpinning the importance of a regulatory mechanism within BS cells that contributes to the M-specific ZmPEPC1 gene expression.

Carbonic Anhydrases: An Ancient Tool in Calcareous Sponge Biomineralization

Enzymes of the α-carbonic anhydrase gene family (CAs) are essential for the deposition of calcium carbonate biominerals. In calcareous sponges (phylum Porifera, class Calcarea), specific CAs are involved in the formation of calcite spicules, a unique trait and synapomorphy of this class. However, detailed studies on the CA repertoire of calcareous sponges exist for only two species of one of the two Calcarea subclasses, the Calcaronea. The CA repertoire of the second subclass, the Calcinea, has not been investigated so far, leaving a considerable gap in our knowledge about this gene family in Calcarea. Here, using transcriptomic analysis, phylogenetics, and in situ hybridization, we study the CA repertoire of four additional species of calcareous sponges, including three from the previously unsampled subclass Calcinea. Our data indicate that the last common ancestor of Calcarea had four ancestral CAs with defined subcellular localizations and functions (mitochondrial/cytosolic, membrane-bound, and secreted non-catalytic). The evolution of membrane-bound and secreted CAs involved gene duplications and losses, whereas mitochondrial/cytosolic and non-catalytic CAs are evidently orthologous genes. Mitochondrial/cytosolic CAs are biomineralization-specific genes recruited for biomineralization in the last common ancestor of calcareous sponges. The spatial-temporal expression of these CAs differs between species, which may reflect differences between subclasses or be related to the secondary thickening of spicules during biomineralization that does not occur in all species. With this study, we extend the understanding of the role and the evolution of a key biomineralization gene in calcareous sponges.

The evolutionary fate of rpl32 and rps16 losses in the Euphorbia schimperi (Euphorbiaceae) plastome

Gene transfers from mitochondria and plastids to the nucleus are an important process in the evolution of the eukaryotic cell. Plastid (pt) gene losses have been documented in multiple angiosperm lineages and are often associated with functional transfers to the nucleus or substitutions by duplicated nuclear genes targeted to both the plastid and mitochondrion. The plastid genome sequence of Euphorbia schimperi was assembled and three major genomic changes were detected, the complete loss of rpl32 and pseudogenization of rps16 and infA. The nuclear transcriptome of E. schimperi was sequenced to investigate the transfer/substitution of the rpl32 and rps16 genes to the nucleus. Transfer of plastid-encoded rpl32 to the nucleus was identified previously in three families of Malpighiales, Rhizophoraceae, Salicaceae and Passifloraceae. An E. schimperi transcript of pt SOD-1-RPL32 confirmed that the transfer in Euphorbiaceae is similar to other Malpighiales indicating that it occurred early in the divergence of the order. Ribosomal protein S16 (rps16) is encoded in the plastome in most angiosperms but not in Salicaceae and Passifloraceae. Substitution of the E. schimperi pt rps16 was likely due to a duplication of nuclear-encoded mitochondrial-targeted rps16 resulting in copies dually targeted to the mitochondrion and plastid. Sequences of RPS16-1 and RPS16-2 in the three families of Malpighiales (Salicaceae, Passifloraceae and Euphorbiaceae) have high sequence identity suggesting that the substitution event dates to the early divergence within Malpighiales.

ENLARGED STARCH GRAIN1 affects amyloplast development and starch biosynthesis in rice endosperm

Cereal crops accumulate large amounts of starch which is synthesized and stored in amyloplasts in the form of starch grains (SGs). Despite significant progress in deciphering starch biosynthesis, our understanding of amyloplast development in rice (Oryza sativa) endosperm remains largely unknown. Here, we report a novel rice floury mutant named enlarged starch grain1 (esg1). The mutant has decreased starch content, altered starch physicochemical properties, slower grain-filling rate and reduced 1000-grain weight. A distinctive feature in esg1 endosperm is that SGs are much larger, mainly due to an increased number of starch granules per SG. Spherical and loosely assembled granules, together with those weakly stained SGs may account for decreased starch content in esg1. Map-based cloning revealed that ESG1 encodes a putative permease subunit of a bacterial-type ABC (ATP-binding cassette) lipid transporter. ESG1 is constitutively expressed in various tissues. It encodes a protein localized to the chloroplast and amyloplast membranes. Mutation of ESG1 causes defective galactolipid synthesis. The overall study indicates that ESG1 is a newly identified protein affecting SG development and subsequent starch biosynthesis, which provides novel insights into amyloplast development in rice.

First Draft Genome Sequence Resource of Colletotrichum liriopes Causing Leaf Anthracnose on Ophiopogon japonicus

Leaf anthracnose of Ophiopogon japonicus is an important disease that can significantly reduce the quality and economic value of this traditional Chinese medicinal plant. The disease is caused by Colletotrichum liriopes, a necrotrophic fungus that belongs to the Glomerellaceae family of the Sordariomycetes class. Here, we present the draft whole-genome sequence of the C. liriopes strain A2 that caused leaf anthracnose on O. japonicus. The assembly consists of 407 contigs with an estimated genome size of 53.1 Mb. Furthermore, we identified 670 carbohydrate-active enzymes, 1,377 secreted proteins, and 60 secondary metabolite gene clusters, which may be associated with the pathogenicity of this pathogen. This genome resource will provide a valuable resource for future research on the pathogenesis of C. liriopes and comparative genome analyses within genus Colletotrichum.

Ensemble of Multiple Classifiers for Multilabel Classification of Plant Protein Subcellular Localization

The accurate prediction of protein localization is a critical step in any functional genome annotation process. This paper proposes an improved strategy for protein subcellular localization prediction in plants based on multiple classifiers, to improve prediction results in terms of both accuracy and reliability. The prediction of plant protein subcellular localization is challenging because the underlying problem is not only a multiclass, but also a multilabel problem. Generally, plant proteins can be found in 10–14 locations/compartments. The number of proteins in some compartments (nucleus, cytoplasm, and mitochondria) is generally much greater than that in other compartments (vacuole, peroxisome, Golgi, and cell wall). Therefore, the problem of imbalanced data usually arises. Therefore, we propose an ensemble machine learning method based on average voting among heterogeneous classifiers. We first extracted various types of features suitable for each type of protein localization to form a total of 479 feature spaces. Then, feature selection methods were used to reduce the dimensions of the features into smaller informative feature subsets. This reduced feature subset was then used to train/build three different individual models. In the process of combining the three distinct classifier models, we used an average voting approach to combine the results of these three different classifiers that we constructed to return the final probability prediction. The method could predict subcellular localizations in both single- and multilabel locations, based on the voting probability. Experimental results indicated that the proposed ensemble method could achieve correct classification with an overall accuracy of 84.58% for 11 compartments, on the basis of the testing dataset.

Identification and Expression Characterization of ATP-Binding Cassette (ABC) Transporter Genes in Melon Fly

Simple Summary

The melon fly, Zeugodacus cucurbitae, is an important agricultural pest. At present, chemical pesticide treatment is the main method for field control, but this promotes pesticide resistance by Z. cucurbitae, because of its frequent use. ABC transporters are involved in detoxification metabolism, but few studies have yet considered their expression in melon fly. In this study, we identified the ABC transporters genes at a genome-wide level in melon fly, and analysed their spatiotemporal expression patterns, as well as changes in expression after insecticides treatments. A total of 49 ABC transporters were identified, and their expression levels varied at different developmental stages and between tissues. After three insecticides treatment, ZcABCB7 and ZcABCC2 were up-regulated. After β-cypermethrin induction, tissues were dissected at 12, 24 and 48 h, and the expression levels of a number of ABC genes were highly expressed within the fat body. From these results, we conclude that ZcABCB7 and ZcABCC2 may be involved in detoxification metabolism, and that the fat body is the main tissue that plays this role.

Abstract

The ATP-binding cassette (ABC) transporter is a protein superfamily that transports specific substrate molecules across lipid membranes in all living species. In insects, ABC transporter is one of the major transmembrane protein families involved in the development of xenobiotic resistance. Here, we report 49 ABC transporter genes divided into eight subfamilies (ABCA-ABCH), including seven ABCAs, seven ABCBs, 10 ABCCs, two ABCDs, one ABCE, three ABCFs, 16 ABCGs, and three ABCHs according to phylogenetic analysis in Zeugodacus cucurbitae, a highly destructive insect pest of cucurbitaceous and other related crops. The expressions level of 49 ABC transporters throughout various developmental stages and within different tissues were evaluated by quantitative transcriptomic analysis, and their expressions in response to three different insecticides were evaluated by quantitative real-time polymerase chain reaction (qRT-PCR). These ABC transporter genes were widely expressed at developmental stages but most highly expressed in tissues of the midgut, fat body and Malpighian tube. When challenged by exposure to three insecticides, abamectin, β-cypermethrin, and dinotefuran, the expressions of ZcABCB7 and ZcABCC2 were significantly up-regulated. ZcABCB1, ZcABCB6, ZcABCB7, ZcABCC2, ZcABCC3, ZcABCC4, ZcABCC5, and ZcABCC7 were significantly up-regulated in the fat body at 24 h after β-cypermethrin exposure. These data suggest that ZcABCB7 and ZcABCC2 might play key roles in xenobiotic metabolism in Z. cucurbitae. Collectively, these data provide a foundation for further analysis of ABCs in Z. cucurbitae.