Locating proteins in the cell using TargetP, SignalP and related tools

Study of lipase producing gene in wheat - an in silico approach

BACKGROUND

Lipases (EC 3.1.1.3) catalyze the hydrolysis of oil into free fatty acids and glycerol forming the 3rd largest group of commercialized enzymes. Plant lipases grab attention recently because of their specificity, less production and purified cost, and easy availability. In silico approach is the first step to identify different genes coding for lipase in a most common indigenous plant, wheat, to explore the possibility of this plant as an alternative source for commercial lipase production. As the hierarchy organization of genes reflects an ancient process of gene duplication and divergence, many of the theoretical and analytical tools of the phylogenetic systematics can be utilized for comparative genomic studies. Also, in addition to experimental identification and characterization of genes, for computational genomic analysis, Arabidopsis has become a popular strategy to identify crop genes which are economically important, as Arabidopsis genes had been well identified and characterized for lipase. A number of articles had been reported in which genes of wheat have shown strong homology with Arabidopsis. The complete genome sequences of rice and Arabidopsis constitute a valuable resource for comparative genome analysis as they are representatives of the two major evolutionary lineages within the angiosperms. Here, in this in silico approach, Arabidopsis and Oryza sativa serve as models for dicotyledonous and monocotyledonous species, respectively, and the genomic sequence data available was used to identify the lipase genes in wheat.

RESULTS

In this present study, Ensembl Plants database was explored for lipase producing gene present in wheat genome and 21 genes were screened down as they contain specific domain and motif for lipase (GXSXG). According to the evolutionary analysis, it was found that the gene TraesCS5B02G157100, located in 5B chromosome, has 58.35% sequence similarity with the reported lipase gene of Arabidopsis thaliana and gene TraesCS3A02G463500 located in the 3A chromosome has 51.74% sequence similarity with the reported lipase gene of Oryza sativa. Homology modeling was performed using protein sequences coded by aforementioned genes and optimized by molecular dynamic simulations. Further with the help of molecular docking of modeled structures with tributyrin, binding efficiency was checked, and the difference in energies (DE) was -9.83 kcal/mol and -6.67 kcal/mol, respectively.

CONCLUSIONS

The present work provides a basic understanding of the gene-encoding lipase in wheat, which could be easily accessible and used as a potent industrial enzyme. The study enlightens another direction which can be used further to explore plant lipases.

Potential role of chimeric genes in pathway-related gene co-expression modules

Background

Gene fusion has epigenetic modification functions. The novel proteins encoded by gene fusion products play a role in cancer development. Therefore, a better understanding of the novel protein products may provide insights into the pathogenesis of tumors. However, the characteristics of chimeric genes are rarely studied. Here, we used weighted co-expression network analysis to investigate the biological roles and underlying mechanisms of chimeric genes.

Methods

Download the pig transcriptome data, we screened chimeric genes and parental genes from 688 sequences and 153 samples, predict their domains, and analyze their associations. We constructed a co-expression network of chimeric genes in pigs and conducted Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analysis on the generated modules using DAVID to identify key networks and modules related to chimeric genes.

Results

Our findings showed that most of the protein domains of chimeric genes were derived from fused pre-genes. Chimeric genes were enriched in modules involved in the negative regulation of cell proliferation and protein localization to centrosomes. In addition, the chimeric genes were related to the growth factor-β superfamily, which regulates cell growth and differentiation. Furthermore, in helper T cells, chimeric genes regulate the specific recognition of T cell receptors, implying that chimeric genes play a key role in the regulation pathway of T cells. Chimeric genes can produce new domains, and some chimeric genes are a key role involved in pathway-related function.

Conclusions

Most chimeric genes show binding activity. Domains of chimeric genes are derived from several combinations of parent genes. Chimeric genes play a key role in the regulation of several cellular pathways. Our findings may provide new directions to explore the roles of chimeric genes in tumors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12957-021-02248-9.

HBD1 protein with a tandem repeat of two HMG-box domains is a DNA clip to organize chloroplast nucleoids in Chlamydomonas reinhardtii

Compaction of bulky DNA is a universal issue for all DNA-based life forms. Chloroplast nucleoids (chloroplast DNA-protein complexes) are critical for chloroplast DNA maintenance and transcription, thereby supporting photosynthesis, but their detailed structure remains enigmatic. Our proteomic analysis of chloroplast nucleoids of the green alga Chlamydomonas reinhardtii identified a protein (HBD1) with a tandem repeat of two DNA-binding high mobility group box (HMG-box) domains, which is structurally similar to major mitochondrial nucleoid proteins transcription factor A, mitochondrial (TFAM), and ARS binding factor 2 protein (Abf2p). Disruption of the HBD1 gene by CRISPR-Cas9-mediated genome editing resulted in the scattering of chloroplast nucleoids. This phenotype was complemented when intact HBD1 was reintroduced, whereas a truncated HBD1 with a single HMG-box domain failed to complement the phenotype. Furthermore, ectopic expression of HBD1 in the mitochondria of yeast Δabf2 mutant successfully complemented the defects, suggesting functional similarity between HBD1 and Abf2p. Furthermore, in vitro assays of HBD1, including the electrophoretic mobility shift assay and DNA origami/atomic force microscopy, showed that HBD1 is capable of introducing U-turns and cross-strand bridges, indicating that proteins with two HMG-box domains would function as DNA clips to compact DNA in both chloroplast and mitochondrial nucleoids.

OsGERLP: A novel aluminum tolerance rice gene isolated from a local cultivar in Indonesia

There is a decrease in the land available for rice cultivation due to the rapid conversion to urban uses. Subsequently, acid soil could be an alternative land cultivating rice, but will require the use of aluminum (Al)-tolerant rice varieties. This Al tolerance trait is genetically controlled, and there is a need to discover more genes needed to develop Al-tolerant rice. Therefore, the objective of this study was to clone and characterize a novel Al tolerance gene isolated from a local cultivar of Indonesian rice. The gene cloning was conducted based on the rye/rice microsynteny relationship. In addition, the root growth and gene expression analyses were performed to verify the role of the gene on Al tolerance in gene-silenced rice and in overexpressed transgenic tobacco. The results showed an Al tolerance candidate gene, OsGERLP, was successfully cloned from rice cv. Hawara Bunar, with its gene encoding a protein similar to a bacterial ribosomal L32 protein. Additionally, the analysis showed that low gene expression caused the gene-silenced rice to be sensitive to Al, while high expression induced the Al tolerance in transgenic tobacco. Furthermore, it was discovered that the gene expression level in both plants was in line with the lower expression of the OsFRDL4 gene in the silenced rice and the high expression of the MATE gene in transgenic tobacco also with the higher citrate secretion from transgenic tobacco roots. In conclusion, the OsGERLP gene could act as a regulator for other Al tolerance genes, with the potential to develop Al-tolerant rice varieties.

Uncovering genetic mechanisms of hypertension through multi-omic analysis of the kidney

The kidney is an organ of key relevance to blood pressure (BP) regulation, hypertension and antihypertensive treatment. However, genetically mediated renal mechanisms underlying susceptibility to hypertension remain poorly understood. We integrated genotype, gene expression, alternative splicing and DNA methylation profiles of up to 430 human kidneys to characterize the effects of BP index variants from genome-wide association studies (GWASs) on renal transcriptome and epigenome. We uncovered kidney targets for 479 (58.3%) BP-GWAS variants and paired 49 BP-GWAS kidney genes with 210 licensed drugs. Our colocalization and Mendelian randomization analyses identified 179 unique kidney genes with evidence of putatively causal effects on BP. Through Mendelian randomization, we also uncovered effects of BP on renal outcomes commonly affecting patients with hypertension. Collectively, our studies identified genetic variants, kidney genes, molecular mechanisms and biological pathways of key relevance to the genetic regulation of BP and inherited susceptibility to hypertension.

Plasma membrane phylloquinone biosynthesis in nonphotosynthetic parasitic plants

Nonphotosynthetic holoparasites exploit flexible targeting of phylloquinone biosynthesis to facilitate plasma membrane redox signaling. Phylloquinone is a lipophilic naphthoquinone found predominantly in chloroplasts and best known for its function in photosystem I electron transport and disulfide bridge formation of photosystem II subunits. Phylloquinone has also been detected in plasma membrane (PM) preparations of heterotrophic tissues with potential transmembrane redox function, but the molecular basis for this noncanonical pathway is unknown. Here, we provide evidence of PM phylloquinone biosynthesis in a nonphotosynthetic holoparasite Phelipanche aegyptiaca. A nonphotosynthetic and nonplastidial role for phylloquinone is supported by transcription of phylloquinone biosynthetic genes during seed germination and haustorium development, by PM-localization of alternative terminal enzymes, and by detection of phylloquinone in germinated seeds. Comparative gene network analysis with photosynthetically competent parasites revealed a bias of P. aegyptiaca phylloquinone genes toward coexpression with oxidoreductases involved in PM electron transport. Genes encoding the PM phylloquinone pathway are also present in several photoautotrophic taxa of Asterids, suggesting an ancient origin of multifunctionality. Our findings suggest that nonphotosynthetic holoparasites exploit alternative targeting of phylloquinone for transmembrane redox signaling associated with parasitism.

Genome-wide identification, characterization, and expression analysis of the monovalent cation-proton antiporter superfamily in maize, and functional analysis of its role in salt tolerance

Na⁺, K⁺ and pH homeostasis are important for plant life and they are controlled by the monovalent cation proton antiporter (CPA) superfamily. The roles of ZmCPAs in salt tolerance are not fully elucidated. In this study, we identified 35 ZmCPAs comprising 13 Na⁺/H⁺ exchangers (ZmNHXs), 16 cation/H⁺ exchanger (ZmCHXs), and 6 K⁺ efflux antiporters (ZmKEAs). All ZmCPAs have transmembrane domains and most of them were localized to plasma membrane or tonoplast. ZmCHXs were specifically highly expressed in anthers, while ZmNHXs and ZmKEAs showed high expression in various tissues. ZmNHX5 and ZmKEA2 were up-regulated in maize seedlings under both NaCl and KCl stresses. Yeast complementation experiments revealed the roles of ZmNHX5, ZmKEA2 in NaCl tolerance. Analysis of the maize mutants further validated the salt tolerance functions of ZmNHX5 and ZmKEA2. Our study highlights comprehensive information of ZmCPAs and provides new gene targets for salt tolerance maize breeding.

The Chemosensory Repertoire of the Eastern Diamondback Rattlesnake (Crotalus adamanteus) Reveals Complementary Genetics of Olfactory and Vomeronasal-Type Receptors

Pitviper sensory perception incorporates diverse stimuli through the integration of trichromatic color vision, bifocal heat-sensing, and dual-system chemoperception. Chemoperception, or olfaction, is mediated by chemoreceptors in the olfactory bulb and the vomeronasal organ, but the true genomic complexity of the gene families and their relative contributions is unknown. A full genomic accounting of pitviper chemoperception directly complements our current understanding of their venoms by generating a more complete polyphenic representation of their predatory arsenal. To characterize the genetic repertoire of pitviper chemoperception, we analyzed a full-genome assembly for Crotalus adamanteus, the eastern diamondback rattlesnake. We identified hundreds of genes encoding both olfactory receptors (ORs; 362 full-length genes) and type-2 vomeronasal receptors (V2Rs; 430 full-length genes). Many chemoreceptor genes are organized into large tandem repeat arrays. Comparative analysis of V2R orthologs across squamates demonstrates how gene array expansion and contraction underlies the evolution of the chemoreceptor repertoire, which likely reflects shifts in life history traits. Chromosomal assignments of chemosensory genes identified sex chromosome specific chemoreceptor genes, providing gene candidates underlying observed sex-specific chemosensory-based behaviors. We detected widespread episodic evolution in the extracellular, ligand-binding domains of both ORs and V2Rs, suggesting the diversification of chemoreceptors is driven by transient periods of positive selection. We provide a robust genetic framework for studying pitviper chemosensory ecology and evolution.

CwlQ Is Required for Swarming Motility but Not Flagellar Assembly in Bacillus subtilis

Lytic enzymes play an essential role in the remodeling of bacterial peptidoglycan (PG), an extracellular mesh-like structure that retains the membrane in the context of high internal osmotic pressure. Peptidoglycan must be unfailingly stable to preserve cell integrity, but must also be dynamically remodeled for the cell to grow, divide, and insert macromolecular machines. The flagellum is one such macromolecular machine that transits the PG, and flagellar insertion is aided by localized activity of a dedicated PG lyase in Gram-negative bacteria. To date, there is no known dedicated lyase in Gram-positive bacteria for the insertion of flagella. Here, we take a reverse-genetic candidate-gene approach and find that cells mutated for the lytic transglycosylase CwlQ exhibit a severe defect in flagellum-dependent swarming motility. We further show that CwlQ is expressed by the motility sigma factor SigD and is secreted by the type III secretion system housed inside the flagellum. Nonetheless, cells with mutations of CwlQ remain proficient for flagellar biosynthesis even when mutated in combination with four other lyases related to motility (LytC, LytD, LytF, and CwlO). The PG lyase (or lyases) essential for flagellar synthesis in B. subtilis, if any, remains unknown.IMPORTANCE Bacteria are surrounded by a wall of peptidoglycan and early work in Bacillus subtilis was the first to suggest that bacteria needed to enzymatically remodel the wall to permit insertion of the flagellum. No PG remodeling enzyme alone or in combination, however, has been found to be essential for flagellar assembly in B. subtilis Here, we take a reverse-genetic candidate-gene approach and find that the PG lytic transglycosylase CwlQ is required for swarming motility. Subsequent characterization determined that while CwlQ was coexpressed with motility genes and is secreted by the flagellar secretion apparatus, it was not required for flagellar synthesis. The PG lyase needed for flagellar assembly in B. subtilis remains unknown.

Anti-inflammatory effects of FS48, the first potassium channel inhibitor from the salivary glands of the flea Xenopsylla cheopis

The voltage-gated potassium (Kv) 1.3 channel plays a crucial role in the immune responsiveness of T-lymphocytes and macrophages, presenting a potential target for treatment of immune- and inflammation related-diseases. FS48, a protein from the rodent flea Xenopsylla cheopis, shares the three disulfide bond feature of scorpion toxins. However, its three-dimensional structure and biological function are still unclear. In the present study, the structure of FS48 was evaluated by circular dichroism and homology modeling. We also described its in vitro ion channel activity using patch clamp recording and investigated its anti-inflammatory activity in LPS-induced Raw 264.7 macrophage cells and carrageenan-induced paw edema in mice. FS48 was found to adopt a common αββ structure and contain an atypical dyad motif. It dose-dependently exhibited the Kv1.3 channel in Raw 264.7 and HEK 293T cells, and its ability to block the channel pore was demonstrated by the kinetics of activation and competition binding with tetraethylammonium. FS48 also downregulated the secretion of proinflammatory molecules NO, IL-1β, TNF-α, and IL-6 by Raw 264.7 cells in a manner dependent on Kv1.3 channel blockage and the subsequent inactivation of the MAPK/NF-κB pathways. Finally, we observed that FS48 inhibited the paw edema formation, tissue myeloperoxidase activity, and inflammatory cell infiltrations in carrageenan-treated mice. We therefore conclude that FS48 identified from the flea saliva is a novel potassium channel inhibitor displaying anti-inflammatory activity. This discovery will promote understanding of the bloodsucking mechanism of the flea and provide a new template molecule for the design of Kv1.3 channel blockers.

Selective axonal translation of the mRNA isoform encoding prenylated Cdc42 supports axon growth

The small Rho-family GTPase Cdc42 has long been known to have a role in cell motility and axon growth. The eukaryotic Ccd42 gene is alternatively spliced to generate mRNAs with two different 3' untranslated regions (UTRs) that encode proteins with distinct C-termini. The C-termini of these Cdc42 proteins include CaaX and CCaX motifs for post-translational prenylation and palmitoylation, respectively. Palmitoyl-Cdc42 protein was previously shown to contribute to dendrite maturation, while the prenyl-Cdc42 protein contributes to axon specification and its mRNA was detected in neurites. Here, we show that the mRNA encoding prenyl-Cdc42 isoform preferentially localizes into PNS axons and this localization selectively increases in vivo during peripheral nervous system (PNS) axon regeneration. Functional studies indicate that prenyl-Cdc42 increases axon length in a manner that requires axonal targeting of its mRNA, which, in turn, needs an intact C-terminal CaaX motif that can drive prenylation of the encoded protein. In contrast, palmitoyl-Cdc42 has no effect on axon growth but selectively increases dendrite length. Together, these data show that alternative splicing of the Cdc42 gene product generates an axon growth promoting, locally synthesized prenyl-Cdc42 protein. This article has an associated First Person interview with one of the co-first authors of the paper.

Deciphering the Monilinia fructicola Genome to Discover Effector Genes Possibly Involved in Virulence

Brown rot is the most economically important fungal disease of stone fruits and is primarily caused by Monilinia laxa and Monlinia fructicola. Both species co-occur in European orchards although M. fructicola is considered to cause the most severe yield losses in stone fruit. This study aimed to generate a high-quality genome of M. fructicola and to exploit it to identify genes that may contribute to pathogen virulence. PacBio sequencing technology was used to assemble the genome of M. fructicola. Manual structural curation of gene models, supported by RNA-Seq, and functional annotation of the proteome yielded 10,086 trustworthy gene models. The genome was examined for the presence of genes that encode secreted proteins and more specifically effector proteins. A set of 134 putative effectors was defined. Several effector genes were cloned into Agrobacterium tumefaciens for transient expression in Nicotiana benthamiana plants, and some of them triggered necrotic lesions. Studying effectors and their biological properties will help to better understand the interaction between M. fructicola and its stone fruit host plants.

Vestiges of the Bacterial Signal Recognition Particle-Based Protein Targeting in Mitochondria

The main bacterial pathway for inserting proteins into the plasma membrane relies on the signal recognition particle (SRP), composed of the Ffh protein and an associated RNA component, and the SRP-docking protein FtsY. Eukaryotes use an equivalent system of archaeal origin to deliver proteins into the endoplasmic reticulum, whereas a bacteria-derived SRP and FtsY function in the plastid. Here we report on the presence of homologs of the bacterial Ffh and FtsY proteins in various unrelated plastid-lacking unicellular eukaryotes, namely Heterolobosea, Alveida, Goniomonas, and Hemimastigophora. The monophyly of novel eukaryotic Ffh and FtsY groups, predicted mitochondrial localization experimentally confirmed for Naegleria gruberi, and a strong alphaproteobacterial affinity of the Ffh group, collectively suggest that they constitute parts of an ancestral mitochondrial signal peptide-based protein-targeting system inherited from the last eukaryotic common ancestor, but lost from the majority of extant eukaryotes. The ability of putative signal peptides, predicted in a subset of mitochondrial-encoded N. gruberi proteins, to target a reporter fluorescent protein into the endoplasmic reticulum of Trypanosoma brucei, likely through their interaction with the cytosolic SRP, provided further support for this notion. We also illustrate that known mitochondrial ribosome-interacting proteins implicated in membrane protein targeting in opisthokonts (Mba1, Mdm38, and Mrx15) are broadly conserved in eukaryotes and nonredundant with the mitochondrial SRP system. Finally, we identified a novel mitochondrial protein (MAP67) present in diverse eukaryotes and related to the signal peptide-binding domain of Ffh, which may well be a hitherto unrecognized component of the mitochondrial membrane protein-targeting machinery.

Coding and non-coding roles of MOCCI (C15ORF48) coordinate to regulate host inflammation and immunity

Mito-SEPs are small open reading frame-encoded peptides that localize to the mitochondria to regulate metabolism. Motivated by an intriguing negative association between mito-SEPs and inflammation, here we screen for mito-SEPs that modify inflammatory outcomes and report a mito-SEP named "Modulator of cytochrome C oxidase during Inflammation" (MOCCI) that is upregulated during inflammation and infection to promote host-protective resolution. MOCCI, a paralog of the NDUFA4 subunit of cytochrome C oxidase (Complex IV), replaces NDUFA4 in Complex IV during inflammation to lower mitochondrial membrane potential and reduce ROS production, leading to cyto-protection and dampened immune response. The MOCCI transcript also generates miR-147b, which targets the NDUFA4 mRNA with similar immune dampening effects as MOCCI, but simultaneously enhances RIG-I/MDA-5-mediated viral immunity. Our work uncovers a dual-component pleiotropic regulation of host inflammation and immunity by MOCCI (C15ORF48) for safeguarding the host during infection and inflammation.

Sequence and structure-based method to predict diacylglycerol lipases in protein sequence

Lipase enzymes play a central role in biotechnology and the food industry. Diacylglyceride lipases (DAG) have received considerable attention due to their physiological significance and potential industrial usage. However, compared to the wide application of triacylglycerol (TAG) lipases, DAG lipases have a limited application due to their low thermostability and specific activity. The molecular basis of substrate specificity of DAG lipases remains elusive, making structure-guided engineering of TAG to DAG lipase difficult. Besides, the number of available DAG lipases is limited compared to TAG lipases. In the current study, we identified structural consensus motifs of DAG lipases that contribute to their DAG specificity on a structural comparison of DAG and TAG lipases. To find potential DAG lipases, sequence motifs and predicted secondary structures were used to screen millions of protein sequences and predict new DAG lipases. In total, 83 new putative DAG lipases were identified. The predicted DAG lipases were validated by expression of randomly chosen putative DAG lipases followed by functional assay for their DAG and TAG specific activity. The reported method is efficient and cost-effective for discovering new DAG lipases used in the food industry after the required characterization to meet potential application needs.

Hydrolysis of organophosphorus by diatom purple acid phosphatase and sequential regulation of cell metabolism

Phosphorus (P) limitation affects phytoplankton growth and population size in aquatic systems, and consequently limits aquatic primary productivity. Plants have evolved a range of metabolic responses to cope with P limitation, such as accumulation of purple acid phosphatases (PAPs) to enhance acquisition of phosphates. However, it remains unknown whether algae have evolved a similar mechanism. In this study, we examined the role of PAPs in the model microalga Phaeodactylum tricornutum. Expression of PAP1 was enhanced in P. tricornutum cells grown on organophosphorus compared to inorganic phosphate. PAP1 overexpression improved cellular growth and biochemical composition in a growth-phase dependent manner. PAP1 promoted growth and photosynthesis during growth phases and reallocated carbon flux towards lipogenesis during the stationary phase. PAP1 was found to be localized in the endoplasmic reticulum and it orchestrated the expression of genes involved in key metabolic pathways and translocation of inorganic P (Pi), thereby improving energy use, reducing equivalents and antioxidant potential. RNAi of PAP1 induced expression of its homolog PAP2, thereby compensating for the Pi scavenging activity of PAP1. Our results demonstrate that PAP1 brings about sequential regulation of metabolism, and provide novel insights into algal phosphorus metabolism and aquatic primary productivity.

Genome-wide identification and expression analysis of detoxification efflux carriers (DTX) genes family under abiotic stresses in flax

The detoxification efflux carriers (DTX)/multidrug and toxic compound extrusion (MATE) transporters encompass an ancient gene family of secondary transporters involved in the process of plant detoxification. A genome-wide analysis of these transporters was carried out in order to better understand the transport of secondary metabolites in flaxseed genome (Linum usitassimum). A total of 73 genes coding for DTX/MATE transporters were identified. Gene structure, protein domain and motif organization were found to be notably conserved over the distinct phylogenetic groups, showing the evolutionary significant role of each class. Gene ontology (GO) annotation revealed a link to transporter activities, response to stimulus and localizations. The presence of various hormone and stress-responsive cis-regulatory elements in promoter regions could be directly correlated with the alteration of their transcripts. Tertiary structure showed conservation for pore size and constrains in the pore, which indicate their involvement in the exclusion of toxic substances from the cell. MicroRNA target analysis revealed that LuDTXs genes were targeted by different classes of miRNA families. Twelve LuDTX genes were chosen for further quantitative real-time polymerase chain reaction analysis in response to cold, salinity and cadmium stress at 0, 6, 12 and 24 hours after treatment. Altogether, the identified members of the DTX gene family, their expression profile, phylogenetic and miRNAs analysis might provide opportunities for future functional validation of this important gene family in flax.

Genome-wide Identification and Characterization of FCS-Like Zinc Finger (FLZ) Family Genes in Maize (Zea mays) and Functional Analysis of ZmFLZ25 in Plant Abscisic Acid Response

FCS-like zinc finger family proteins (FLZs), a class of plant-specific scaffold of SnRK1 complex, are involved in the regulation of various aspects of plant growth and stress responses. Most information of FLZ family genes was obtained from the studies in Arabidopsis thaliana, whereas little is known about the potential functions of FLZs in crop plants. In this study, 37 maize FLZ (ZmFLZ) genes were identified to be asymmetrically distributed on 10 chromosomes and can be divided into three subfamilies. Protein interaction and subcellular localization assays demonstrated that eight typical ZmFLZs interacted and partially co-localized with ZmKIN10, the catalytic α-subunit of the SnRK1 complex in maize leaf mesophyll cells. Expression profile analysis revealed that several ZmFLZs were differentially expressed across various tissues and actively responded to diverse abiotic stresses. In addition, ectopic overexpression of ZmFLZ25 in Arabidopsis conferred hypersensitivity to exogenous abscisic acid (ABA) and triggered higher expression of ABA-induced genes, pointing to the positive regulatory role of ZmFLZ25 in plant ABA signaling, a scenario further evidenced by the interactions between ZmFLZ25 and ABA receptors. In summary, these data provide the most comprehensive information on FLZ family genes in maize, and shed light on the biological function of ZmFLZ25 in plant ABA signaling.

Prunus Hexokinase 3 genes alter primary C-metabolism and promote drought and salt stress tolerance in Arabidopsis transgenic plants

Hexokinases (HXKs) and fructokinases (FRKs) are the only two families of enzymes in plants that have been identified as able to phosphorylate Glucose (Glc) and Fructose (Fru). Glc can only be phosphorylated in plants by HXKs, while Fru can be phosphorylated by either HXKs or FRKs. The various subcellular localizations of HXKs in plants indicate that they are involved in diverse functions, including anther dehiscence and pollen germination, stomatal closure in response to sugar levels, stomatal aperture and reducing transpiration. Its association with modulating programmed cell death, and responses to oxidative stress and pathogen infection (abiotic and biotic stresses) also have been reported. To extend our understanding about the function of HXK-like genes in the response of Prunus rootstocks to abiotic stress, we performed a detailed bioinformatic and functional analysis of hexokinase 3-like genes (HXK3s) from two Prunus rootstock genotypes, 'M.2624' (Prunus cerasifera Ehrh × P. munsoniana W.Wight & Hedrick) and 'M.F12/1' (P. avium L.), which are tolerant and sensitive to hypoxia stress, respectively. A previous large-scale transcriptome sequencing of roots of these rootstocks, showed that this HXK3-like gene that was highly induced in the tolerant genotype under hypoxia conditions. In silico analysis of gene promoters from M.2624 and M.F12/1 genotypes revealed regulatory elements that could explain differential transcriptional profiles of HXK3 genes. Subcellular localization was determinates by both bioinformatic prediction and expression of their protein fused to the green fluorescent protein (GFP) in protoplasts and transgenic plants of Arabidopsis. Both approaches showed that they are expressed in plastids. Metabolomics analysis of Arabidopsis plants ectopically expressing Prunus HXK3 genes revealed that content of several metabolites including phosphorylated sugars (G6P), starch and some metabolites associated with the TCA cycle were affected. These transgenic Arabidopsis plants showed improved tolerance to salt and drought stress under growth chamber conditions. Our results suggest that Prunus HXK3 is a potential candidate for enhancing tolerance to salt and drought stresses in stone fruit trees and other plants.

Genome-scale metabolic model of the diatom Thalassiosira pseudonana highlights the importance of nitrogen and sulfur metabolism in redox balance

Diatoms are unicellular photosynthetic algae known to secrete organic matter that fuels secondary production in the ocean, though our knowledge of how their physiology impacts the composition of dissolved organic matter remains limited. Like all photosynthetic organisms, their use of light for energy and reducing power creates the challenge of avoiding cellular damage. To better understand the interplay between redox balance and organic matter secretion, we reconstructed a genome-scale metabolic model of Thalassiosira pseudonana strain CCMP 1335, a model for diatom molecular biology and physiology, with a 60-year history of studies. The model simulates the metabolic activities of 1,432 genes via a network of 2,792 metabolites produced through 6,079 reactions distributed across six subcellular compartments. Growth was simulated under different steady-state light conditions (5–200 μmol photons m^-2 s^-1) and in a batch culture progressing from exponential growth to nitrate-limitation and nitrogen-starvation. We used the model to examine the dissipation of reductants generated through light-dependent processes and found that when available, nitrate assimilation is an important means of dissipating reductants in the plastid; under nitrate-limiting conditions, sulfate assimilation plays a similar role. The use of either nitrate or sulfate uptake to balance redox reactions leads to the secretion of distinct organic nitrogen and sulfur compounds. Such compounds can be accessed by bacteria in the surface ocean. The model of the diatom Thalassiosira pseudonana provides a mechanistic explanation for the production of ecologically and climatologically relevant compounds that may serve as the basis for intricate, cross-kingdom microbial networks. Diatom metabolism has an important influence on global biogeochemistry; metabolic models of marine microorganisms link genes to ecosystems and may be key to integrating molecular data with models of ocean biogeochemistry.

Single dose of a replication-defective vaccinia virus expressing Zika virus-like particles is protective in mice

Zika virus (ZIKV), a flavivirus transmitted primarily by infected mosquitos, can cause neurological symptoms such as Guillian-Barré syndrome and microcephaly. We developed several vaccinia virus (VACV) vaccine candidates for ZIKV based on replication-inducible VACVs (vINDs) expressing ZIKV pre-membrane (prM) and envelope (E) proteins (vIND-ZIKVs). These vIND-ZIKVs contain elements of the tetracycline operon and replicate only in the presence of tetracyclines. The pool of vaccine candidates was narrowed to one vIND-ZIKV containing a novel mutation in the signal peptide of prM that led to higher expression and secretion of E and production of virus-like particles, which was then tested for safety, immunogenicity, and efficacy in mice. vIND-ZIKV grows to high titers in vitro in the presence of doxycycline (DOX) but is replication-defective in vivo in the absence of DOX, causing no weight loss in mice. C57BL/6 mice vaccinated once with vIND-ZIKV in the absence of DOX (as a replication-defective virus) developed robust levels of E-peptide-specific IFN-γ-secreting splenocytes and anti-E IgG titers, with modest levels of serum-neutralizing antibodies. Vaccinated mice treated with anti-IFNAR1 antibody were completely protected from ZIKV viremia post-challenge after a single dose of vIND-ZIKV. Furthermore, mice with prior immunity to VACV developed moderate anti-E IgG titers that increased after booster vaccination, and were protected from viremia only after two vaccinations with vIND-ZIKV.

The iron-sulfur scaffold protein HCF101 unveils the complexity of organellar evolution in SAR, Haptista and Cryptista

Background

Nbp35-like proteins (Nbp35, Cfd1, HCF101, Ind1, and AbpC) are P-loop NTPases that serve as components of iron-sulfur cluster (FeS) assembly machineries. In eukaryotes, Ind1 is present in mitochondria, and its function is associated with the assembly of FeS clusters in subunits of respiratory Complex I, Nbp35 and Cfd1 are the components of the cytosolic FeS assembly (CIA) pathway, and HCF101 is involved in FeS assembly of photosystem I in plastids of plants (chHCF101). The AbpC protein operates in Bacteria and Archaea. To date, the cellular distribution of these proteins is considered to be highly conserved with only a few exceptions.

Results

We searched for the genes of all members of the Nbp35-like protein family and analyzed their targeting sequences. Nbp35 and Cfd1 were predicted to reside in the cytoplasm with some exceptions of Nbp35 localization to the mitochondria; Ind1was found in the mitochondria, and HCF101 was predicted to reside in plastids (chHCF101) of all photosynthetically active eukaryotes. Surprisingly, we found a second HCF101 paralog in all members of Cryptista, Haptista, and SAR that was predicted to predominantly target mitochondria (mHCF101), whereas Ind1 appeared to be absent in these organisms. We also identified a few exceptions, as apicomplexans possess mHCF101 predicted to localize in the cytosol and Nbp35 in the mitochondria. Our predictions were experimentally confirmed in selected representatives of Apicomplexa (Toxoplasma gondii), Stramenopila (Phaeodactylum tricornutum, Thalassiosira pseudonana), and Ciliophora (Tetrahymena thermophila) by tagging proteins with a transgenic reporter. Phylogenetic analysis suggested that chHCF101 and mHCF101 evolved from a common ancestral HCF101 independently of the Nbp35/Cfd1 and Ind1 proteins. Interestingly, phylogenetic analysis supports rather a lateral gene transfer of ancestral HCF101 from bacteria than its acquisition being associated with either α-proteobacterial or cyanobacterial endosymbionts.

Conclusion

Our searches for Nbp35-like proteins across eukaryotic lineages revealed that SAR, Haptista, and Cryptista possess mitochondrial HCF101. Because plastid localization of HCF101 was only known thus far, the discovery of its mitochondrial paralog explains confusion regarding the presence of HCF101 in organisms that possibly lost secondary plastids (e.g., ciliates, Cryptosporidium) or possess reduced nonphotosynthetic plastids (apicomplexans).

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01777-x.

Iron metabolism strategies in diatoms

Diatoms are one of the most successful group of photosynthetic eukaryotes in the contemporary ocean. They are ubiquitously distributed and are the most abundant primary producers in polar waters. Equally remarkable is their ability to tolerate iron deprivation and respond to periodic iron fertilization. Despite their relatively large cell sizes, diatoms tolerate iron limitation and frequently dominate iron-stimulated phytoplankton blooms, both natural and artificial. Here, we review the main iron use strategies of diatoms, including their ability to assimilate and store a range of iron sources, and the adaptations of their photosynthetic machinery and architecture to iron deprivation. Our synthesis relies on published literature and is complemented by a search of 82 diatom transcriptomes, including information collected from seven representatives of the most abundant diatom genera in the world's oceans.

ATP binding cassette importers in eukaryotic organisms

ATP-binding cassette (ABC) transporters are ubiquitous across all realms of life. Dogma suggests that bacterial ABC transporters include both importers and exporters, whilst eukaryotic members of this family are solely exporters, implying that ABC import function was lost during evolution. This view is being challenged, for example energy-coupling factor (ECF)-type ABC importers appear to fulfil important roles in both algae and plants where they form the ABCI sub-family. Herein we discuss whether bacterial Type I and Type II ABC importers also made the transition into extant eukaryotes. Various studies suggest that Type I importers exist in algae and the liverwort family of primitive non-vascular plants, but not in higher plants. The existence of eukaryotic Type II importers is also supported: a transmembrane protein homologous to vitamin B12 import system transmembrane protein (BtuC), hemin transport system transmembrane protein (HmuU) and high-affinity zinc uptake system membrane protein (ZnuB) is present in the Cyanophora paradoxa genome. This protein has homologs within the genomes of red algae. Furthermore, its candidate nucleotide-binding domain (NBD) shows closest similarity to other bacterial Type II importer NBDs such as BtuD. Functional studies suggest that Type I importers have roles in maintaining sulphate levels in the chloroplast, whilst Type II importers probably act as importers of Mn²⁺ or Zn²⁺ , as inferred by comparisons with bacterial homologs. Possible explanations for the presence of these transporters in simple plants, but not in other eukaryotic organisms, are considered. In order to utilise the existing nomenclature for eukaryotic ABC proteins, we propose that eukaryotic Type I and II importers be classified as ABCJ and ABCK transporters, respectively.

Isolation and characterization of a γ-carbonic anhydrase localized in the mitochondria of Saccharina japonica

Saccharina japonica is an ecologically and economically important seaweed that is dominant in the rocky shores of cold-temperate regions, forms the major component of productive beds, and affects marine environments. S. japonica exhibits a high photosynthetic efficiency in natural seawater with low dissolved CO₂ concentration, thus suggesting the presence of its carbon-concentrating mechanism (CCM). However, the genes, proteins, and pathways involved in the CCM of S. japonica have not been fully identified and characterized. Carbonic anhydrase (CA) is a crucial component of CCM in macroalgae. In this study, the cloning, characterization, and subcellular localization of a specific CA were described. Multisequence alignment and phylogenetic analysis indicated that this CA belonged to the gamma (Sjγ-CA) class. This enzyme has a full-length cDAN of 1370 bp, encodes a protein with 246 amino acids (aa; ca. 25.7 kDa), and contains the mitochondrial transit peptide of 16 aa and LbH_gama_CA_like domain of 159 aa that defined the γ-CA region. The Sjγ-CA was successfully expressed in E. coli BL21 and purified as an active recombinant CA. Immunogold electron microscopy and fluorescence localization illustrated that this enzyme is localized in the mitochondria, and its transcription level is up-regulated by low CO₂ concentration. These findings showed that Sjγ-CA is a possible component of the CCM in S. japonica. This work is the first to report about the mtCA of macroalgae and provides a basis for further analysis on seaweed CCM.

An apoplastic fluid extraction method for the characterization of grapevine leaves proteome and metabolome from a single sample

The analysis of complex biological systems keeps challenging researchers. The main goal of systems biology is to decipher interactions within cells, by integrating datasets from large scale analytical approaches including transcriptomics, proteomics and metabolomics and more specialized 'OMICS' such as epigenomics and lipidomics. Studying different cellular compartments allows a broader understanding of cell dynamics. Plant apoplast, the cellular compartment external to the plasma membrane including the cell wall, is particularly demanding to analyze. Despite our knowledge on apoplast involvement on several processes from cell growth to stress responses, its dynamics is still poorly known due to the lack of efficient extraction processes adequate to each plant system. Analyzing woody plants such as grapevine raises even more challenges. Grapevine is among the most important fruit crops worldwide and a wider characterization of its apoplast is essential for a deeper understanding of its physiology and cellular mechanisms. Here, we describe, for the first time, a vacuum-infiltration-centrifugation method that allows a simultaneous extraction of grapevine apoplastic proteins and metabolites from leaves on a single sample, compatible with high-throughput mass spectrometry analyses. The extracted apoplast from two grapevine cultivars, Vitis vinifera cv 'Trincadeira' and 'Regent', was directly used for proteomics and metabolomics analysis. The proteome was analyzed by nanoLC-MS/MS and more than 700 common proteins were identified, with highly diverse biological functions. The metabolome profile through FT-ICR-MS allowed the identification of 514 unique putative compounds revealing a broad spectrum of molecular classes.

Kinases and protein motifs required for AZI1 plastid localization and trafficking during plant defense induction

The proper subcellular localization of defense factors is an important part of the plant immune system. A key component for systemic resistance, lipid transfer protein (LTP)-like AZI1, is needed for the systemic movement of the priming signal azelaic acid (AZA) and a pool of AZI1 exists at the site of AZA production, the plastid envelope. Moreover, after systemic defense-triggering infections, the proportion of AZI1 localized to plastids increases. However, AZI1 does not possess a classical plastid transit peptide that can explain its localization. Instead, AZI1 uses a bipartite N-terminal signature that allows for its plastid targeting. Furthermore, the kinases MPK3 and MPK6, associated with systemic immunity, promote the accumulation of AZI1 at plastids during priming induction. Our results indicate the existence of a mode of plastid targeting possibly related to defense responses.

Genome-wide identification and expression analysis of the NF-Y transcription factor family in Populus

In the past few years, many studies have reported that the transcription factor Nuclear Factor Y (NF-Y) gene family plays important roles in embryonic development, photosynthesis, flowering time regulation and stress response, in various plants. Although the NF-Y gene family has been systematically studied in many species, little is known about NF-Y genes in Populus. In this study, the NF-Y gene family in the Populus genome was identified and its structural characteristics were described. Fifty-two NF-Y genes were authenticated in the Populus trichocarpa genome and categorized into three subfamilies (NF-YA/B/C) by phylogenetic analysis. Chromosomal localization of these genes revealed that they were distributed randomly across 17 of the 19 chromosomes. Segmental duplication played a vital role in the amplification of Populus NF-Y gene family. Moreover, microsynteny analysis indicated that, among Populus trichocarpa, Arabidopsis thaliana, Vitis vinifera and Carica papaya, NF-Y duplicated regions were more conserved between Populus trichocarpa and Vitis vinifera. Redundant stress-related cis-elements were also found in the promoters of most 13 NF-YA genes and their expression levels varied widely following drought, salt, ABA and cold treatments. Subcellular localization experiments in tobacco showed that PtNF-YA3 was localized in nucleus and cytomembrane, while PtNF-YA4 was only in the nucleus in tobacco. According to the transcriptional activity experiments, neither of them had transcriptional activity in yeast. In summary, a comprehensive analysis of the Populus NF-Y gene family was performed to establish a theoretical basis for further functional studies on this family.

The Cytospora chrysosperma Virulence Effector CcCAP1 Mainly Localizes to the Plant Nucleus To Suppress Plant Immune Responses

Canker disease is caused by the fungus Cytospora chrysosperma and damages a wide range of woody plants, causing major losses to crops and native plants. Plant pathogens secrete virulence-related effectors into host cells during infection to regulate plant immunity and promote colonization. However, the functions of C. chrysosperma effectors remain largely unknown. In this study, we used Agrobacterium tumefaciens-mediated transient expression system in Nicotiana benthamiana and confocal microscopy to investigate the immunoregulation roles and subcellular localization of CcCAP1, a virulence-related effector identified in C. chrysosperma CcCAP1 was significantly induced in the early stages of infection and contains cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 proteins (CAP) superfamily domain with four cysteines. CcCAP1 suppressed the programmed cell death triggered by Bcl-2-associated X protein (BAX) and the elicitin infestin1 (INF1) in transient expression assays with Nicotiana benthamiana The CAP superfamily domain was sufficient for its cell death-inhibiting activity and three of the four cysteines in the CAP superfamily domain were indispensable for its activity. Pathogen challenge assays in N. benthamiana demonstrated that transient expression of CcCAP1 promoted Botrytis cinerea infection and restricted reactive oxygen species accumulation, callose deposition, and defense-related gene expression. In addition, expression of green fluorescent protein-labeled CcCAP1 in N. benthamiana showed that it localized to both the plant nucleus and the cytoplasm, but the nuclear localization was essential for its full immune inhibiting activity. These results suggest that this virulence-related effector of C. chrysosperma modulates plant immunity and functions mainly via its nuclear localization and the CAP domain.IMPORTANCE The data presented in this study provide a key resource for understanding the biology and molecular basis of necrotrophic pathogen responses to Nicotiana benthamiana resistance utilizing effector proteins, and CcCAP1 may be used in future studies to understand effector-triggered susceptibility processes in the Cytospora chrysosperma-poplar interaction system.

The E3-ligase AoUBR1 in N-end rule pathway is involved in the vegetative growth, secretome, and trap formation in Arthrobotrys oligospora

The N-end rule pathway is a regulated protein degradation system. Arthrobotrys oligospora, a typical nematode-trapping fungus, switches its life strategies from saprophytism to carnivorism when capturing free-living nematodes by means of adhesive networks. In this study, a putative E3-ligase AoUBR1 involved in N-end rule pathway was characterized in A. oligospora during vegetative growth and trap formation. Expression of AoUBR1 coding gene was down-regulated during trap formation. Compared with wild type, the AoUBR1 knock-out mutants decreased the vegetative growth, formed less traps, and turned to be sensitive to cold stress, while, AoUBR1 overexpression mutants lost the capacity to produce conidia and also formed less traps. A number of genes differentially expressed by knock-out and overexpression of AoUBR1, which lead to the transcriptional responses associated with plasma membrane, transportation, oxidation, and proteolysis. AoUBR1 knock-out also resulted in the down-regulation of numerous secreted proteins associated with carnivorism and nutrient utilization from nematodes. In addition, AoUBR1 homologs were conserved in nematode-trapping fungi based on the genome searching. Therefore, the results suggested AoUBR1 in A. oligospora and its homologs in other trapping fungi are involved in the lifestyle switch between saprophytism and carnivorism.

Retortamonads from vertebrate hosts share features of anaerobic metabolism and pre-adaptations to parasitism with diplomonads

Although the mitochondria of extant eukaryotes share a single origin, functionally these organelles diversified to a great extent, reflecting lifestyles of the organisms that host them. In anaerobic protists of the group Metamonada, mitochondria are present in reduced forms (also termed hydrogenosomes or mitosomes) and a complete loss of mitochondrion in Monocercomonoides exilis (Metamonada:Preaxostyla) has also been reported. Within metamonads, retortamonads from the gastrointestinal tract of vertebrates form a sister group to parasitic diplomonads (e.g. Giardia and Spironucleus) and have also been hypothesized to completely lack mitochondria. We obtained transcriptomic data from Retortamonas dobelli and R. caviae and searched for enzymes of the core metabolism as well as mitochondrion- and parasitism-related proteins. Our results indicate that retortamonads have a streamlined metabolism lacking pathways for metabolites they are probably capable of obtaining from prey bacteria or their environment, reminiscent of the biochemical arrangement in other metamonads. Retortamonads were surprisingly found do encode homologs of components of Giardia's remarkable ventral disk, as well as homologs of regulatory NEK kinases and secreted lytic enzymes known for involvement in host colonization by Giardia. These can be considered pre-adaptations of these intestinal microorganisms to parasitism. Furthermore, we found traces of the mitochondrial metabolism represented by iron‑sulfur cluster assembly subunits, subunits of mitochondrial translocation and chaperone machinery and, importantly, [FeFe]‑hydrogenases and hydrogenase maturases (HydE, HydF and HydG). Altogether, our results strongly suggest that a remnant mitochondrion is still present.

A novel allele of the Arabidopsis thaliana MACPF protein CAD1 results in deregulated immune signaling

Immune recognition in plants is governed by two major classes of receptors: pattern recognition receptors (PRRs) and nucleotide-binding leucine-rich repeat receptors (NLRs). Located at the cell surface, PRRs bind extracellular ligands originating from microbes (indicative of "non-self") or damaged plant cells (indicative of "infected-self"), and trigger signaling cascades to protect against infection. Located intracellularly, NLRs sense pathogen-induced physiological changes and trigger localized cell death and systemic resistance. Immune responses are under tight regulation in order to maintain homeostasis and promote plant health. In a forward-genetic screen to identify regulators of PRR-mediated immune signaling, we identified a novel allele of the membrane-attack complex and perforin (MACPF)-motif containing protein CONSTITUTIVE ACTIVE DEFENSE 1 (CAD1) resulting from a missense mutation in a conserved N-terminal cysteine. We show that cad1-5 mutants display deregulated immune signaling and symptoms of autoimmunity dependent on the lipase-like protein ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1), suggesting that CAD1 integrity is monitored by the plant immune system. We further demonstrate that CAD1 localizes to both the cytosol and plasma membrane using confocal microscopy and subcellular fractionation. Our results offer new insights into immune homeostasis and provide tools to further decipher the intriguing role of MACPF proteins in plants.

Creatine kinase B controls futile creatine cycling in thermogenic fat

Obesity increases the risk of mortality because of metabolic sequelae such as type 2 diabetes and cardiovascular disease¹. Thermogenesis by adipocytes can counteract obesity and metabolic diseases^2,3. In thermogenic fat, creatine liberates a molar excess of mitochondrial ADP-purportedly via a phosphorylation cycle⁴-to drive thermogenic respiration. However, the proteins that control this futile creatine cycle are unknown. Here we show that creatine kinase B (CKB) is indispensable for thermogenesis resulting from the futile creatine cycle, during which it traffics to mitochondria using an internal mitochondrial targeting sequence. CKB is powerfully induced by thermogenic stimuli in both mouse and human adipocytes. Adipocyte-selective inactivation of Ckb in mice diminishes thermogenic capacity, increases predisposition to obesity, and disrupts glucose homeostasis. CKB is therefore a key effector of the futile creatine cycle.

Genomic rearrangements generate hypervariable mini-chromosomes in host-specific isolates of the blast fungus

Supernumerary mini-chromosomes–a unique type of genomic structural variation–have been implicated in the emergence of virulence traits in plant pathogenic fungi. However, the mechanisms that facilitate the emergence and maintenance of mini-chromosomes across fungi remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), mini-chromosomes have been first described in the early 1990s but, until very recently, have been overlooked in genomic studies. Here we investigated structural variation in four isolates of the blast fungus M. oryzae from different grass hosts and analyzed the sequences of mini-chromosomes in the rice, foxtail millet and goosegrass isolates. The mini-chromosomes of these isolates turned out to be highly diverse with distinct sequence composition. They are enriched in repetitive elements and have lower gene density than core-chromosomes. We identified several virulence-related genes in the mini-chromosome of the rice isolate, including the virulence-related polyketide synthase Ace1 and two variants of the effector gene AVR-Pik. Macrosynteny analyses around these loci revealed structural rearrangements, including inter-chromosomal translocations between core- and mini-chromosomes. Our findings provide evidence that mini-chromosomes emerge from structural rearrangements and segmental duplication of core-chromosomes and might contribute to adaptive evolution of the blast fungus.

The genomes of plant pathogens often exhibit an architecture that facilitates high rates of dynamic rearrangements and genetic diversification in virulence associated regions. These regions, which tend to be gene sparse and repeat rich, are thought to serve as a cradle for adaptive evolution. Supernumerary chromosomes, i.e. chromosomes that are only present in some but not all individuals of a species, are a special type of structural variation that have been observed in plants, animals, and fungi. Here we identified and studied supernumerary mini-chromosomes in the blast fungus Magnaporthe oryzae, a pathogen that causes some of the most destructive plant diseases. We found that rice, foxtail millet and goosegrass isolates of this pathogen contain mini-chromosomes with distinct sequence composition. All mini-chromosomes are rich in repetitive genetic elements and have lower gene densities than core-chromosomes. Further, we identified virulence-related genes on the mini-chromosome of the rice isolate. We observed large-scale genomic rearrangements around these loci, indicative of a role of mini-chromosomes in facilitating genome dynamics. Taken together, our results indicate that mini-chromosomes contribute to genome rearrangements and possibly adaptive evolution of the blast fungus.

Evaluation of CBSX Proteins as Regulators of the Chloroplast Thioredoxin System

The chloroplast-localized cystathionine β-synthase X (CBSX) proteins CBSX1 and CBSX2 have been proposed as modulators of thioredoxins (Trxs). In this study, the contribution of CBSX proteins to the redox regulation of thiol enzymes in the chloroplast Trx system was evaluated both in vitro and in vivo. The in vitro biochemical studies evaluated whether CBSX proteins alter the specificities of classical chloroplastic Trx f and Trx m for their target proteins. However, addition of CBSX proteins did not alter the specificities of Trx f and Trx m for disulfide bond reduction of the photosynthesis-related major thiol enzymes, FBPase, SBPase, and NADP-MDH. In vivo analysis showed that CBSX-deficient mutants grew similarly to wild type plants under continuous normal light conditions and that CBSX deficiency did not affect photo-reduction of photosynthesis-related thiol enzymes by Trx system at several light intensities. Although CBSX proteins have been suggested as modulators in the chloroplast Trx system, our results did not support this model, at least in the cases of FBPase, SBPase, and NADP-MDH in leaves. However, fresh weights of the cbsx2 mutants were decreased under short day. Since Trxs regulate many proteins participating in various metabolic reactions in the chloroplast, CBSX proteins may function to regulate other chloroplast Trx target proteins, or serve as modulators in non-photosynthetic plastids of flowers. As a next stage, further investigations are required to understand the modulation of Trx-dependent redox regulation by plastidal CBSX proteins.

The genetic basis of cytoplasmic male sterility and fertility restoration in wheat

Hybrid wheat varieties give higher yields than conventional lines but are difficult to produce due to a lack of effective control of male fertility in breeding lines. One promising system involves the Rf1 and Rf3 genes that restore fertility of wheat plants carrying Triticum timopheevii-type cytoplasmic male sterility (T-CMS). Here, by genetic mapping and comparative sequence analyses, we identify Rf1 and Rf3 candidates that can restore normal pollen production in transgenic wheat plants carrying T-CMS. We show that Rf1 and Rf3 bind to the mitochondrial orf279 transcript and induce cleavage, preventing expression of the CMS trait. The identification of restorer genes in wheat is an important step towards the development of hybrid wheat varieties based on a CMS-Rf system. The characterisation of their mode of action brings insights into the molecular basis of CMS and fertility restoration in plants.

The Molecular Function of Plant mTERFs as Key Regulators of Organellar Gene Expression

The protein family of mTERFs (mitochondrial transcription termination factors) was initially studied in mammalian and insect mitochondria before the first Arabidopsis mTERF mutant was characterized. More than 10 years of research on the function of plant mTERFs in the flowering plants Arabidopsis thaliana, Zea mays and the green microalga Chlamydomonas reinhardtii has since highlighted that mTERFs are key regulators of organellar gene expression (OGE) in mitochondria and in chloroplasts. Additional functions to be fulfilled by plant mTERFs (e.g. splicing) and the fact that the expression of two organellar genomes had to be facilitated have led to a massive expansion of the plant mTERF portfolio compared to that found in mammals. Plant mTERFs are implicated in all steps of OGE ranging from the modulation of transcription to the maturation of tRNAs and hence translation. Furthermore, being regulators of OGE, mTERFs are required for a successful long-term acclimation to abiotic stress, retrograde signaling and interorganellar communication. Here, I review the recent progress in the elucidation of molecular mTERF functions.

Identification and Characterization of a Periplasmic α-Carbonic Anhydrase (CA) in the Gametophytes of Saccharina japonica (Phaeophyceae)

Periplasmic or external carbonic anhydrases (CAs) have been well accepted as playing a crucial role in the acquisition of dissolved inorganic carbon; however, no cytological evidence or molecular information on these enzymes has been reported in seaweeds to date. In this study, the full-length cDNA sequence coding for a putative periplasmic Sjα-CA2 was cloned from the gametophytes of Saccharina japonica, an industrial brown seaweed. It was 1,728 bp in length and included a 263-bp 5'-untranslated region (UTR), a 577-bp 3'-UTR, and an 888-bp open reading frame encoding a protein precursor consisting of 295 amino acids. The mature protein, after removal of a predicted 28-residue signal peptide, was composed of 267 amino acids with a relative molecular weight of 29.27 kDa. Multisequence alignment and phylogenetic analysis indicated that it was a member of the α-CA family. Enzyme activity assays showed that the recombinant Sjα-CA2 in Escherichia coli possessed CO2 hydration and esterase activities, thus identifying this gene Sjα-CA2 in function. Immunogold electron microscopic observations with the prepared anti-Sjα-CA2 polyclonal antibody illustrated that Sjα-CA2 was located in periplasmic space of the kelp gametophyte cells. Quantitative real-time PCR results revealed that the transcription of Sjα-CA2 was induced by elevated HCO 3 - levels, but it was little changed while the kelp gametophytes were subjected to elevated CO2 concentrations. This study suggests that the periplasmic Sjα-CA2 might play a role in adapting to elevated environmental levels of HCO 3 - by dehydration of HCO 3 - to generate CO2 , which could be readily taken up by S. japonica gametophytes.

Identification and characterization of the LRR repeats in plant LRR-RLKs

Background

Leucine-rich-repeat receptor-like kinases (LRR-RLKs) play central roles in sensing various signals to regulate plant development and environmental responses. The extracellular domains (ECDs) of plant LRR-RLKs contain LRR motifs, consisting of highly conserved residues and variable residues, and are responsible for ligand perception as a receptor or co-receptor. However, there are few comprehensive studies on the ECDs of LRR-RLKs due to the difficulty in effectively identifying the divergent LRR repeats.

Results

In the current study, an efficient LRR motif prediction program, the “Phyto-LRR prediction” program, was developed based on the position-specific scoring matrix algorithm (PSSM) with some optimizations. This program was trained by 16-residue plant-specific LRR-highly conserved segments (HCS) from LRR-RLKs of 17 represented land plant species and a database containing more than 55,000 predicted LRRs based on this program was constructed. Both the prediction tool and database are freely available at http://phytolrr.com/ for website usage and at http://github.com/phytolrr for local usage. The LRR-RLKs were classified into 18 subgroups (SGs) according to the maximum-likelihood phylogenetic analysis of kinase domains (KDs) of the sequences. Based on the database and the SGs, the characteristics of the LRR motifs in the ECDs of the LRR-RLKs were examined, such as the arrangement of the LRRs, the solvent accessibility, the variable residues, and the N-glycosylation sites, revealing a comprehensive profile of the plant LRR-RLK ectodomains.

Conclusion

The “Phyto-LRR prediction” program is effective in predicting the LRR segments in plant LRR-RLKs, which, together with the database, will facilitate the exploration of plant LRR-RLKs functions. Based on the database, comprehensive sequential characteristics of the plant LRR-RLK ectodomains were profiled and analyzed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00344-y.

Identification and development of novel salt-responsive candidate gene based SSRs (cg-SSRs) and MIR gene based SSRs (mir-SSRs) in bread wheat (Triticum aestivum)

Salt stress adversely affects the global wheat production and productivity. To improve salinity tolerance of crops, identification of robust molecular markers is highly imperative for development of salt-tolerant cultivars to mimic yield losses under saline conditions. In this study, we mined 171 salt-responsive genes (including 10 miRNAs) from bread wheat genome using the sequence information of functionally validated salt-responsive rice genes. Salt-stress, tissue and developmental stage-specific expression analysis of RNA-seq datasets revealed the constitutive as well as the inductive response of salt-responsive genes in different tissues of wheat. Fifty-four genotypes were phenotyped for salt stress tolerance. The stress tolerance index of the genotypes ranged from 0.30 to 3.18. In order to understand the genetic diversity, candidate gene based SSRs (cg-SSRs) and MIR gene based SSRs (miR-SSRs) were mined from 171 members of salt-responsive genes of wheat and validated among the contrasting panels of 54 tolerant as well as susceptible wheat genotypes. Among 53 SSR markers screened, 10 cg-SSRs and 8 miR-SSRs were found to be polymorphic. Polymorphic information content between the wheat genotypes ranged from 0.07 to 0.67, indicating the extant of wide genetic variation among the salt tolerant and susceptible genotypes at the DNA level. The genetic diversity analysis based on the allelic data grouped the wheat genotypes into three separate clusters of which single group encompassing most of the salt susceptible genotypes and two of them containing salt tolerance and moderately salt tolerance wheat genotypes were in congruence with penotypic data. Our study showed that both salt-responsive genes and miRNAs based SSRs were more diverse and can be effectively used for diversity analysis. This study reports the first extensive survey on genome-wide analysis, identification, development and validation of salt-responsive cg-SSRs and miR-SSRs in wheat. The information generated in the present study on genetic divergence among genotypes having a differential response to salt will help in the selection of suitable lines as parents for developing salt tolerant cultivars in wheat.

Genome-wide identification, evolution, expression, and alternative splicing profiles of peroxiredoxin genes in cotton

Peroxiredoxin (PRX) is a ubiquitous thioredoxin-dependent peroxidase that can eliminate excessive free radicals produced by stress and protect cells from oxidative damage. PRXs are also involved in reactive oxygen species (ROS)- and redox-dependent signaling by performing redox interactions with other proteins and modify their redox status. At present, PRX family identification, evolution and regulation research has been conducted in some plants; however, systematic research about this family is lacking in cotton. In this study, a total of 44 PRXs were identified in the cotton genome. Phylogenetic and conserved active site analyses showed that the PRXs were divided into six subfamilies according to the conserved site (PxxxTxxC…S…W/F) and conserved cysteinyl residues positions. Segmental duplication and polyploid events were the main methods for PRX family expansion, and the PRXs of diploid G. arboreum were the donors of PRXs in the D subgenomes of allotetraploid G. hirsutum and G. barbadense during the evolution of the PRX family. qRT-PCR analysis confirmed that cis-acting elements play important roles in regulating the expression of PRXs. Alternative splicing events occurred in GhPRX14-D that can increased the complexity of transcripts in G. hirsutum. Subcellular localization showed that most PRX members were located in chloroplasts, the cytoplasmic membrane and the nucleus. Our results provide systematic support for a better understanding of PRXs in cotton and a starting point for further studies of the specific functions of PRXs in cotton.

An 8.22 Mb Assembly and Annotation of the Alpaca (Vicugna pacos) Y Chromosome

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca (Vicugna pacos) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

Enlisting the Ixodes scapularis Embryonic ISE6 Cell Line to Investigate the Neuronal Basis of Tick-Pathogen Interactions

Neuropeptides are small signaling molecules expressed in the tick central nervous system, i.e., the synganglion. The neuronal-like Ixodes scapularis embryonic cell line, ISE6, is an effective tool frequently used for examining tick–pathogen interactions. We detected 37 neuropeptide transcripts in the I. scapularis ISE6 cell line using in silico methods, and six of these neuropeptide genes were used for experimental validation. Among these six neuropeptide genes, the tachykinin-related peptide (TRP) of ISE6 cells varied in transcript expression depending on the infection strain of the tick-borne pathogen, Anaplasma phagocytophilum. The immunocytochemistry of TRP revealed cytoplasmic expression in a prominent ISE6 cell subpopulation. The presence of TRP was also confirmed in A. phagocytophilum-infected ISE6 cells. The in situ hybridization and immunohistochemistry of TRP of I. scapularis synganglion revealed expression in distinct neuronal cells. In addition, TRP immunoreaction was detected in axons exiting the synganglion via peripheral nerves as well as in hemal nerve-associated lateral segmental organs. The characterization of a complete Ixodes neuropeptidome in ISE6 cells may serve as an effective in vitro tool to study how tick-borne pathogens interact with synganglion components that are vital to tick physiology. Therefore, our current study is a potential stepping stone for in vivo experiments to further examine the neuronal basis of tick–pathogen interactions.

The Arabidopsis NOT4A E3 ligase promotes PGR3 expression and regulates chloroplast translation

Chloroplast function requires the coordinated action of nuclear- and chloroplast-derived proteins, including several hundred nuclear-encoded pentatricopeptide repeat (PPR) proteins that regulate plastid mRNA metabolism. Despite their large number and importance, regulatory mechanisms controlling PPR expression are poorly understood. Here we show that the Arabidopsis NOT4A ubiquitin-ligase positively regulates the expression of PROTON GRADIENT REGULATION 3 (PGR3), a PPR protein required for translating several thylakoid-localised photosynthetic components and ribosome subunits within chloroplasts. Loss of NOT4A function leads to a strong depletion of cytochrome b₆f and NAD(P)H dehydrogenase (NDH) complexes, as well as plastid 30 S ribosomes, which reduces mRNA translation and photosynthetic capacity, causing pale-yellow and slow-growth phenotypes. Quantitative transcriptome and proteome analysis of the not4a mutant reveal it lacks PGR3 expression, and that its molecular defects resemble those of a pgr3 mutant. Furthermore, we show that normal plastid function is restored to not4a through transgenic PGR3 expression. Our work identifies NOT4A as crucial for ensuring robust photosynthetic function during development and stress-response, through promoting PGR3 production and chloroplast translation.

Epitope Analysis and Efficacy Evaluation of Phosphatase 2C (PP2C) DNA Vaccine Against Toxoplasma gondii Infection

Toxoplasma gondii infects almost all warm-blooded animals and negatively affects the health of a wide range of these animals, including humans. Protein phosphatase 2C (PP2C) is a T. gondii protein secreted by rhoptry organelles during host cell invasion. However, very little is known about whether this protein can induce protective immunity against T. gondii. In this study, bioinformatics analysis of PP2C revealed some useful information in the context of anti-toxoplasmosis treatments and vaccine research. In addition, the PP2C gene was amplified, and a eukaryotic expression vector (pEGFP-PP2C) was successfully constructed to express PP2C. Finally, the constructed pEGFP-PP2C was injected into mice to evaluate whether it could induce immunoprotection. Compared with the control groups, we found that immunizations with the pEGFP-PP2C plasmid could elicit specific IgG antibodies and cytokines against T. gondii infection. The survival of mice immunized with the pEGFP-PP2C plasmid was significantly prolonged compared with that of the control group mice. Based on the ability of pEGFP-PP2C to induce specific immune responses against T. gondii, we propose that PP2C merits consideration as a potential vaccine candidate against toxoplasmosis.

The IMG/M data management and analysis system v.6.0: new tools and advanced capabilities

The Integrated Microbial Genomes & Microbiomes system (IMG/M: https://img.jgi.doe.gov/m/) contains annotated isolate genome and metagenome datasets sequenced at the DOE's Joint Genome Institute (JGI), submitted by external users, or imported from public sources such as NCBI. IMG v 6.0 includes advanced search functions and a new tool for statistical analysis of mixed sets of genomes and metagenome bins. The new IMG web user interface also has a new Help page with additional documentation and webinar tutorials to help users better understand how to use various IMG functions and tools for their research. New datasets have been processed with the prokaryotic annotation pipeline v.5, which includes extended protein family assignments.

Computational Analysis of African Swine Fever Virus Protein Space for the Design of an Epitope-Based Vaccine Ensemble

African swine fever virus is the etiological agent of African swine fever, a transmissible severe hemorrhagic disease that affects pigs, causing massive economic losses. There is neither a treatment nor a vaccine available, and the only method to control its spread is by extensive culling of pigs. So far, classical vaccine development approaches have not yielded sufficiently good results in terms of concomitant safety and efficacy. Nowadays, thanks to advances in genomic and proteomic techniques, a reverse vaccinology strategy can be explored to design alternative vaccine formulations. In this study, ASFV protein sequences were analyzed using an in-house pipeline based on publicly available immunoinformatic tools to identify epitopes of interest for a prospective vaccine ensemble. These included experimentally validated sequences from the Immune Epitope Database, as well as de novo predicted sequences. Experimentally validated and predicted epitopes were prioritized following a series of criteria that included evolutionary conservation, presence in the virulent and currently circulating variant Georgia 2007/1, and lack of identity to either the pig proteome or putative proteins from pig gut microbiota. Following this strategy, 29 B-cell, 14 CD4⁺ T-cell and 6 CD8⁺ T-cell epitopes were selected, which represent a starting point to investigating the protective capacity of ASFV epitope-based vaccines.

A Lineage-Specific Paralog of Oma1 Evolved into a Gene Family from Which a Suppressor of Male Sterility-Inducing Mitochondria Emerged in Plants

Cytoplasmic male sterility (MS) in plants is caused by MS-inducing mitochondria, which have emerged frequently during plant evolution. Nuclear restorer-of-fertility (Rf)genes can suppress their cognate MS-inducing mitochondria. Whereas many Rfs encode a class of RNA-binding protein, the sugar beet (Caryophyllales) Rf encodes a protein resembling Oma1, which is involved in the quality control of mitochondria. In this study, we investigated the molecular evolution of Oma1 homologs in plants. We analyzed 37 plant genomes and concluded that a single copy is the ancestral state in Caryophyllales. Among the sugar beet Oma1 homologs, the orthologous copy is located in a syntenic region that is preserved in Arabidopsis thaliana. The sugar beet Rf is a complex locus consisting of a small Oma1 homolog family (RF-Oma1 family) unique to sugar beet. The gene arrangement in the vicinity of the locus is seen in some but not all Caryophyllalean plants and is absent from Ar. thaliana. This suggests a segmental duplication rather than a whole-genome duplication as the mechanism of RF-Oma1 evolution. Of thirty-seven positively selected codons in RF-Oma1, twenty-six of these sites are located in predicted transmembrane helices. Phylogenetic network analysis indicated that homologous recombination among the RF-Oma1 members played an important role to generate protein activity related to suppression. Together, our data illustrate how an evolutionarily young Rf has emerged from a lineage-specific paralog. Interestingly, several evolutionary features are shared with the RNA-binding protein type Rfs. Hence, the evolution of the sugar beet Rf is representative of Rf evolution in general.

Dynamics in Secondary Metabolite Gene Clusters in Otherwise Highly Syntenic and Stable Genomes in the Fungal Genus Botrytis

Fungi of the genus Botrytis infect >1,400 plant species and cause losses in many crops. Besides the broad host range pathogen Botrytis cinerea, most other species are restricted to a single host. Long-read technology was used to sequence genomes of eight Botrytis species, mostly pathogenic on Allium species, and the related onion white rot fungus, Sclerotium cepivorum. Most assemblies contained <100 contigs, with the Botrytis aclada genome assembled in 16 gapless chromosomes. The core genome and pan-genome of 16 Botrytis species were defined and the secretome, effector, and secondary metabolite repertoires analyzed. Among those genes, none is shared among all Allium pathogens and absent from non-Allium pathogens. The genome of each of the Allium pathogens contains 8-39 predicted effector genes that are unique for that single species, none stood out as potential determinant for host specificity. Chromosome configurations of common ancestors of the genus Botrytis and family Sclerotiniaceae were reconstructed. The genomes of B. cinerea and B. aclada were highly syntenic with only 19 rearrangements between them. Genomes of Allium pathogens were compared with ten other Botrytis species (nonpathogenic on Allium) and with 25 Leotiomycetes for their repertoire of secondary metabolite gene clusters. The pattern was complex, with several clusters displaying patchy distribution. Two clusters involved in the synthesis of phytotoxic metabolites are at distinct genomic locations in different Botrytis species. We provide evidence that the clusters for botcinic acid production in B. cinerea and Botrytis sinoallii were acquired by horizontal transfer from taxa within the same genus.