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

Microscopic marvels: Decoding the role of micropeptides in innate immunity

The innate immune response is under selection pressures from changing environments and pathogens. While sequence evolution can be studied by comparing rates of amino acid mutations within and between species, how a gene's birth and death contribute to the evolution of immunity is less known. Short open reading frames, once regarded as untranslated or transcriptional noise, can often produce micropeptides of <100 amino acids with a wide array of biological functions. Some micropeptide sequences are well conserved, whereas others have no evolutionary conservation, potentially representing new functional compounds that arise from species-specific adaptations. To date, few reports have described the discovery of novel micropeptides of the innate immune system. The diversity of immune-related micropeptides is a blind spot for gene and functional annotation. Immune-related micropeptides represent a potential reservoir of untapped compounds for understanding and treating disease. This review consolidates what is currently known about the evolution and function of innate immune-related micropeptides to facilitate their investigation.

Seeing the unseen in characterizing RNA editome during rice endosperm development

Rice (Oryza sativa L.) endosperm is essential to provide nutrients for seed germination and determine grain yield. RNA editing, a post-transcriptional modification essential for plant development, unfortunately, is not fully characterized during rice endosperm development. Here, we perform systematic analyses to characterize RNA editome during rice endosperm development. We find that most editing sites are C-to-U CDS-recoding in mitochondria, leading to increased hydrophobic amino acids and changed structures of mitochondrial proteins. Comparative analysis of RNA editome reveals that CDS-recoding sites present higher editing frequencies with lower variabilities and their resultant recoded amino acids tend to exhibit stronger evolutionary conservation across many land plants. Furthermore, we classify mitochondrial genes into three groups, presenting distinct patterns in terms of CDS-recoding events. Besides, we conduct genome-wide screening to detect pentatricopeptide repeat (PPR) proteins and construct PPR-RNA binding profiles, yielding candidate PPR editing factors related to rice endosperm development. Taken together, our findings provide valuable insights for deciphering fundamental mechanisms of rice endosperm development underlying RNA editing machinery.

A comprehensive in silico analysis of putative outer membrane and secretory hydrolases from the pathogenic Leptospira: Possible implications in pathogenesis

Outer surface/membrane and virulent secretory proteins are primarily crucial for pathogenesis. Secreted and outer membrane hydrolases of many pathogens play an important role in attenuating the host immune system. Leptospira expresses many such proteins, and few have been characterized to display various roles, including host immune evasion. However, identification, classification, characterization, and elucidation of the possible role of Leptospira's outer membrane and secretory hydrolases have yet to be explored. In the present study, we used bioinformatics tools to predict exported proteins from the pathogenic Leptospira proteome. Moreover, we focused on secretory and outer membrane putative hydrolases from the exported proteins to generate a deeper understanding. Our analysis yielded four putative outer/secretory hydrolases, LIC_10995, LIC_11183, LIC_11463, and LIC_12988, containing α/β hydrolase fold and displayed similarity with lipase motif. Moreover, their conservation analysis of the predicted hydrolases across the spectrum of different Leptospira species showed high clustering with the pathogenic species. Outer membrane and secretory proteins with lipolytic activity may have a role in pathogenesis. This is the first bioinformatics analysis of secretory and outer membrane α/β hydrolases from leptospiral species. However, experimental studies are indeed required to unravel this possibility.

ADP-glucose pyrophosphorylase gene family in soybean and implications in drought stress tolerance

BACKGROUND

ADP-glucose pyrophosphorylase (AGPase) is the key rate-limiting enzyme in starch biosynthesis pathway, and has been identified as a potential target for manipulation strategies aimed at improving crop yield and quality.

OBJECTIVE

To identify the AGPase gene family members in soybean, and explore the potential implications of GmAGPS2 in drought stress tolerance.

METHODS

The genome-wide identification and sequence analysis of soybean AGPase gene family was carried out by bioinformatics methods. The GmAGP gene expression was analyzed using transcriptome data and quantitative real-time PCR (qRT-PCR). Furthermore, transgenic yeast strains overexpressing GmAGPS2 were generated, and their growth was observed under drought stress.

RESULTS

In this study, we searched for AGPase genes (GmAGP) in the soybean genome and identified a total of 14 GmAGP genes. The GmAGP proteins had a unique conserved NTP_transferase domain and were mainly located in the chloroplast and cytosol. Evolutionarily, the GmAGP proteins can be clustered into two distinct subgroups; within the same subgroup, they displayed a similar distribution pattern of conserved motifs. The GmAGP genes exhibited an uneven distribution on 10 chromosomes, and segmental duplication contributed to AGPase gene family expansion in soybean. The GmAGP genes presented different tissue expression pattern, in which GmAGPL6, GmAGPL9, and GmAGPL10 mainly exhibited tissue-specific expression pattern. The promoter of GmAGP genes had multiple cis-acting elements related to phytohormones and stress responses, and 8 GmAGP genes contained drought-responsive cis-acting elements. qRT‒PCR analysis demonstrated a significant upregulation expression of GmAGPL6, GmAGPL10, and GmAGPS2 in response to drought stress. Further functional analysis indicated that GmAGPS2 gene could improve yeast growth under drought stress conditions and enhance the drought tolerance of yeast.

CONCLUSION

These results will contribute to further elucidation of the function of GmAGP genes, and offer important candidate genes for the genetic improvement of starch and yield-related traits and the breeding of high drought stress tolerance varieties in soybean.

Protein insertion into the inner membrane of mitochondria: routes and mechanisms

The inner membrane of mitochondria contains hundreds of different integral membrane proteins. These proteins transport molecules into and out of the matrix, they carry out multifold catalytic reactions and they promote the biogenesis or degradation of mitochondrial constituents. Most inner membrane proteins are encoded by nuclear genes and synthesized in the cytosol from where they are imported into mitochondria by translocases in the outer and inner membrane. Three different import routes direct proteins into the inner membrane and allow them to acquire their appropriate membrane topology. First, mitochondrial import intermediates can be arrested at the level of the TIM23 inner membrane translocase by a stop-transfer sequence to reach the inner membrane by lateral insertion. Second, proteins can be fully translocated through the TIM23 complex into the matrix from where they insert into the inner membrane in an export-like reaction. Carriers and other polytopic membrane proteins embark on a third insertion pathway: these hydrophobic proteins employ the specialized TIM22 translocase to insert from the intermembrane space (IMS) into the inner membrane. This review article describes these three targeting routes and provides an overview of the machinery that promotes the topogenesis of mitochondrial inner membrane proteins.

Structural, kinetic, and evolutionary peculiarities of HISN3, a plant 5'-ProFAR isomerase

Histidine biosynthesis is essential for the growth and development of plants, where it occurs within chloroplasts. The eleven reactions are catalyzed by eight enzymes, known as HISN1-8, each acting sequentially. Here, we present the crystal structures of a 5'-ProFAR isomerase (HISN3) from the model legume Medicago truncatula bound to its enzymatically synthesized substrate (ProFAR) and product (PrFAR). The active site of MtHISN3 contains a sodium cation that participates in ligand recognition, a feature not observed in bacterial and fungal structures of homologous enzymes. The steady-state kinetics of wild-type MtHISN3 revealed a slightly higher turnover rate compared to its bacterial homologs. Plant HISN3 sequences contain an unusually elongated Lys60-Ser91 fragment, while deletion of the 74-80 region resulted in a 30-fold loss in catalytic efficiency compared to the wild-type. Molecular dynamics simulations suggested that the fragment facilitates product release, thereby contributing to a higher kcat. Moreover, conservation analyses suggested a non-cyanobacterial origin for plant HISN3 enzymes, which is another instance of a non-cyanobacterial enzyme in the plant histidine biosynthetic pathway. Finally, a virtual screening campaign yielded five molecules, with the energy gains ranging between -13.6 and -13.1 kcal/mol, which provide new scaffolds for the future development of herbicides.

Analysis of the Respiratory Activity in the Antarctic Yeast Rhodotorula mucilaginosa M94C9 Reveals the Presence of Respiratory Supercomplexes and Alternative Elements

The respiratory activities of mitochondrial complexes I, II, and IV were analyzed in permeabilized Rhodotorula mucilaginosa cells and isolated mitochondria, and the kinetic parameters K0.5 and Vmax were obtained. No difference in substrate affinities were found between mitochondria and permeabilized cells. The activities of the components of the mitochondrial respiratory chain of the Antarctic yeast R. mucilaginosa M94C9 were identified by in-gel activity and SDS-PAGE. The mitochondria exhibited activity for the classical components of the electron transport chain (Complexes I, II, III, and IV), and supercomplexes were formed by a combination of the respiratory complexes I, III, and IV. Unfortunately, the activities of the monomeric and dimeric forms of the F1F0-ATP synthase were not revealed by the in-gel assay, but the two forms of the ATP synthase were visualized in the SDS-PAGE. Furthermore, two alternative pathways for the oxidation of cytosolic NADH were identified: the alternative NADH dehydrogenase and the glycerol-3-phosphate dehydrogenase. In addition, an NADPH dehydrogenase and a lactate cytochrome b2 dehydrogenase were found. The residual respiratory activity following cyanide addition suggests the presence of an alternative oxidase in cells.

Prediction of mitochondrial targeting signals and their cleavage sites

In this chapter we survey prediction tools and computational methods for the prediction of amino acid sequence elements which target proteins to the mitochondria. We will primarily focus on the prediction of N-terminal mitochondrial targeting signals (MTSs) and their N-terminal cleavage sites by mitochondrial peptidases. We first give practical details useful for using and installing some prediction tools. Then we describe procedures for preparing datasets of MTS containing proteins for statistical analysis or development of new prediction methods. Following that we lightly survey some of the computational techniques used by prediction tools. Finally, after discussing some caveats regarding the reliability of such methods to predict the effects of mutations on MTS function; we close with a discussion of possible future directions of computer prediction methods related to mitochondrial proteins.

​Fusarium Protein Toolkit: a web-based resource for structural and variant analysis of Fusarium species

BACKGROUND

​​The genus Fusarium poses significant threats to food security and safety worldwide because numerous species of the fungus cause destructive diseases and/or mycotoxin contamination in crops. The adverse effects of climate change are exacerbating some existing threats and causing new problems. These challenges highlight the need for innovative solutions, including the development of advanced tools to identify targets for control strategies.

DESCRIPTION

In response to these challenges, we developed the Fusarium Protein Toolkit (FPT), a web-based tool that allows users to interrogate the structural and variant landscape within the Fusarium pan-genome. The tool displays both AlphaFold and ESMFold-generated protein structure models from six Fusarium species. The structures are accessible through a user-friendly web portal and facilitate comparative analysis, functional annotation inference, and identification of related protein structures. Using a protein language model, FPT predicts the impact of over 270 million coding variants in two of the most agriculturally important species, Fusarium graminearum and F. verticillioides. To facilitate the assessment of naturally occurring genetic variation, FPT provides variant effect scores for proteins in a Fusarium pan-genome based on 22 diverse species. The scores indicate potential functional consequences of amino acid substitutions and are displayed as intuitive heatmaps using the PanEffect framework.

CONCLUSION

FPT fills a knowledge gap by providing previously unavailable tools to assess structural and missense variation in proteins produced by Fusarium. FPT has the potential to deepen our understanding of pathogenic mechanisms in Fusarium, and aid the identification of genetic targets for control strategies that reduce crop diseases and mycotoxin contamination. Such targets are vital to solving the agricultural problems incited by Fusarium, particularly evolving threats resulting from climate change. Thus, FPT has the potential to contribute to improving food security and safety worldwide.

Comparative Genomics of Fungi in Nectriaceae Reveals Their Environmental Adaptation and Conservation Strategies

This study presents the first genome assembly of the freshwater saprobe fungus Neonectria lugdunensis and a comprehensive phylogenomics analysis of the Nectriaceae family, examining genomic traits according to fungal lifestyles. The Nectriaceae family, one of the largest in Hypocreales, includes fungi with significant ecological roles and economic importance as plant pathogens, endophytes, and saprobes. The phylogenomics analysis identified 2684 single-copy orthologs, providing a robust evolutionary framework for the Nectriaceae family. We analyzed the genomic characteristics of 17 Nectriaceae genomes, focusing on their carbohydrate-active enzymes (CAZymes), biosynthetic gene clusters (BGCs), and adaptations to environmental temperatures. Our results highlight the adaptation mechanisms of N. lugdunensis, emphasizing its capabilities for plant litter degradation and enzyme activity in varying temperatures. The comparative genomics of different Nectriaceae lifestyles revealed significant differences in genome size, gene content, repetitive elements, and secondary metabolite production. Endophytes exhibited larger genomes, more effector proteins, and BGCs, while plant pathogens had higher thermo-adapted protein counts, suggesting greater resilience to global warming. In contrast, the freshwater saprobe shows less adaptation to warmer temperatures and is important for conservation goals. This study underscores the importance of understanding fungal genomic adaptations to predict ecosystem impacts and conservation targets in the face of climate change.

Profiling of Phakopsora pachyrhizi transcriptome revealed co-expressed virulence effectors as prospective RNA interference targets for soybean rust management

Soybean rust (SBR), caused by an obligate biotrophic pathogen Phakopsora pachyrhizi, is a devastating disease of soybean worldwide. However, the mechanisms underlying plant invasion by P. pachyrhizi are poorly understood, which hinders the development of effective control strategies for SBR. Here we performed detailed histological characterization on the infection cycle of P. pachyrhizi in soybean and conducted a high-resolution transcriptional dissection of P. pachyrhizi during infection. This revealed P. pachyrhizi infection leads to significant changes in gene expression with 10 co-expressed gene modules, representing dramatic transcriptional shifts in metabolism and signal transduction during different stages throughout the infection cycle. Numerous genes encoding secreted protein are biphasic expressed, and are capable of inhibiting programmed cell death triggered by microbial effectors. Notably, three co-expressed P. pachyrhizi apoplastic effectors (PpAE1, PpAE2, and PpAE3) were found to suppress plant immune responses and were essential for P. pachyrhizi infection. Double-stranded RNA coupled with nanomaterials significantly inhibited SBR infection by targeting PpAE1, PpAE2, and PpAE3, and provided long-lasting protection to soybean against P. pachyrhizi. Together, this study revealed prominent changes in gene expression associated with SBR and identified P. pachyrhizi virulence effectors as promising targets of RNA interference-based soybean protection strategy against SBR.

Arbuscular mycorrhizal fungi alter carbon metabolism and invertase genes expressions of Populus simonii × P. nigra under drought stress

Arbuscular mycorrhizal fungi (AMF) play a crucial role in regulating the allocation of carbon between source and sink tissues in plants and in regulating their stress responses by changing the sucrose biosynthesis, transportation, and catabolism in plants. Invertase, a key enzyme for plant development, participates in the response of plants to drought stress by regulating sucrose metabolism. However, the detailed mechanisms by which INV genes respond to drought stress in mycorrhizal plants remain unclear. This study examined the sugar content, enzyme activity, and expression profiles of INV genes of Populus simonii × P. nigra (PsnINVs) under two inoculation treatments (inoculation or non-inoculation) and two water conditions (well-watered or drought stress). Results showed that under drought stress, AMF up-regulated the expressions of PsnA/NINV1, PsnA/NINV2, PsnA/NINV3, and PsnA/NINV5 in leaves, which may be related to the enhancement of photosynthetic capacity. Additionally, AMF up-regulated the expressions of PsnA/NINV6, PsnA/NINV10, and PsnA/NINV12 in leaves, which may be related to enhancing osmotic regulation ability and drought tolerance.

Strangers in a foreign land: 'Yeastizing' plant enzymes

Expressing plant metabolic pathways in microbial platforms is an efficient, cost-effective solution for producing many desired plant compounds. As eukaryotic organisms, yeasts are often the preferred platform. However, expression of plant enzymes in a yeast frequently leads to failure because the enzymes are poorly adapted to the foreign yeast cellular environment. Here, we first summarize the current engineering approaches for optimizing performance of plant enzymes in yeast. A critical limitation of these approaches is that they are labour-intensive and must be customized for each individual enzyme, which significantly hinders the establishment of plant pathways in cellular factories. In response to this challenge, we propose the development of a cost-effective computational pipeline to redesign plant enzymes for better adaptation to the yeast cellular milieu. This proposition is underpinned by compelling evidence that plant and yeast enzymes exhibit distinct sequence features that are generalizable across enzyme families. Consequently, we introduce a data-driven machine learning framework designed to extract 'yeastizing' rules from natural protein sequence variations, which can be broadly applied to all enzymes. Additionally, we discuss the potential to integrate the machine learning model into a full design-build-test cycle.

The chloroplast pentatricopeptide repeat protein RCN22 regulates tiller number in rice by affecting sugar levels via the TB1-RCN22-RbcL module

As an important yield component, rice tiller number controls panicle number and determines grain yield. Regulation of rice tiller number by chloroplast pentatricopeptide repeat (PPR) proteins has not been reported previously. Here, we report the rice reduced culm number22 (rcn22) mutant, which produces few tillers owing to suppressed tiller bud elongation. Map-based cloning revealed that RCN22 encodes a chloroplast-localized P-type PPR protein. We found that RCN22 specifically binds to the 5' UTR of RbcL mRNA (encoding the large subunit of Rubisco) and enhances its stability. The reduced abundance of RbcL mRNA in rcn22 leads to a lower photosynthetic rate and decreased sugar levels. Consequently, transcript levels of DWARF3 (D3) and TEOSINTE BRANCHED1 (TB1) (which encode negative regulators of tiller bud elongation) are increased, whereas protein levels of the positive regulator DWARF53 (D53) are decreased. Furthermore, high concentrations of sucrose can rescue the tiller bud growth defect of the rcn22 mutant. On the other hand, TB1 directly binds to the RCN22 promoter and downregulates its expression. The tb1/rcn22 double mutant shows a tillering phenotype similar to that of rcn22. Our results suggest that the TB1-RCN22-RbcL module plays a vital role in rice tiller bud elongation by affecting sugar levels.

Analysis and prediction of internal mitochondrial targeting signals

Mitochondria consist of several hundreds of proteins, the vast majority of which are synthesized in the cytosol as precursor proteins from where they are targeted to and imported into mitochondria. The transport of proteins into mitochondria relies on specific targeting information encoded within the protein sequence, known as mitochondrial targeting sequences (MTSs). These N-terminal extensions are usually between 8 and 80 residues long and form amphipathic helices with one hydrophobic and one positively charged surface. Receptors on the mitochondrial surface recognize the MTSs and direct precursors through protein-conducting channels in the outer and inner membrane to the mitochondrial matrix, where presequences are often removed by proteases. In addition to these MTSs, many mitochondrial proteins contain internal matrix targeting sequences (iMTSs) which share the same structural features with MTSs. These iMTSs are neither necessary nor sufficient for mitochondrial targeting, however, they help to increase the import-competence of precursor proteins as they bind to the TOM receptors and presumably facilitate the unfolding of precursors on the mitochondrial surface. Prediction algorithms allow the identification of iMTSs in protein sequences. In this chapter, we present iMLP, an agnostic algorithm for the prediction of iMTS propensity profiles. This iMTS prediction tool is provided via an iMLP webservice at http://iMLP.bio.uni-kl.de and is also available as a BioFSharp application that can be executed locally. We describe and explain the usage of this prediction algorithm and how to interpret the results of this valuable tool.

Complete genome of the Medicago anthracnose fungus, Colletotrichum destructivum, reveals a mini-chromosome-like region within a core chromosome

Colletotrichum destructivum (Cd) is a phytopathogenic fungus causing significant economic losses on forage legume crops (Medicago and Trifolium species) worldwide. To gain insights into the genetic basis of fungal virulence and host specificity, we sequenced the genome of an isolate from Medicago sativa using long-read (PacBio) technology. The resulting genome assembly has a total length of 51.7 Mb and comprises ten core chromosomes and two accessory chromosomes, all of which were sequenced from telomere to telomere. A total of 15, 631 gene models were predicted, including genes encoding potentially pathogenicity-related proteins such as candidate-secreted effectors (484), secondary metabolism key enzymes (110) and carbohydrate-active enzymes (619). Synteny analysis revealed extensive structural rearrangements in the genome of Cd relative to the closely related Brassicaceae pathogen, Colletotrichum higginsianum. In addition, a 1.2 Mb species-specific region was detected within the largest core chromosome of Cd that has all the characteristics of fungal accessory chromosomes (transposon-rich, gene-poor, distinct codon usage), providing evidence for exchange between these two genomic compartments. This region was also unique in having undergone extensive intra-chromosomal segmental duplications. Our findings provide insights into the evolution of accessory regions and possible mechanisms for generating genetic diversity in this asexual fungal pathogen.

Expansion of the MutS Gene Family in Plants

The MutS gene family is distributed across the tree of life and is involved in recombination, DNA repair, and protein translation. Multiple evolutionary processes have expanded the set of MutS genes in plants relative to other eukaryotes. Here, we investigate the origins and functions of these plant-specific genes. Land plants, green algae, red algae, and glaucophytes share cyanobacterial-like MutS1 and MutS2 genes that presumably were gained via plastid endosymbiotic gene transfer. MutS1 was subsequently lost in some taxa, including seed plants, whereas MutS2 was duplicated in Viridiplantae (i.e., land plants and green algae) with widespread retention of both resulting paralogs. Viridiplantae also have two anciently duplicated copies of the eukaryotic MSH6 gene (i.e., MSH6 and MSH7) and acquired MSH1 via horizontal gene transfer - potentially from a nucleocytovirus. Despite sharing the same name, "plant MSH1" is not directly related to the gene known as MSH1 in some fungi and animals, which may be an ancestral eukaryotic gene acquired via mitochondrial endosymbiosis and subsequently lost in most eukaryotic lineages. There has been substantial progress in understanding the functions of MSH1 and MSH6/MSH7 in plants, but the roles of the cyanobacterial-like MutS1 and MutS2 genes remain uncharacterized. Known functions of bacterial homologs and predicted protein structures, including fusions to diverse nuclease domains, provide hypotheses about potential molecular mechanisms. Because most plant-specific MutS proteins are targeted to the mitochondria and/or plastids, the expansion of this family appears to have played a large role in shaping plant organelle genetics.

Homologous Delta-12 Fatty Acid Desaturase (FAD2) Genes Affect Gene Expression and Linoleic Acid Levels in Lentinula edodes under Heat Stress

Delta-12 fatty acid desaturases (FAD2s) actively regulate stress responses and cell differentiation in living organisms. In this study, six homologous FAD2 genes were identified based on the genome sequence of Lentinula edodes. Then, the six FAD2 protein sequences were analyzed using bioinformatics tools, including ExPASy ProtParam, SignalP, TMHMM, and TargetP. These analyses were performed to predict the physical and chemical properties, signal peptides, and transmembrane and conserved domains of these proteins. The polypeptide sequences were aligned, and a maximum likelihood phylogenetic tree was constructed using MEGA 7.0 software to elucidate the phylogenetic relationships between homologous FAD2 sequences. The results demonstrated that the FAD2 proteins contained three conserved histidine-rich regions (HXXXH, HXXHH, and HXXHH), which included eight histidine residues. The linoleic acid content and FAD2 enzyme activity were further analyzed, and the levels in the mutagenic heat-tolerant strain 18N44 were lower than those in the wild-type strain 18. Interestingly, the expression levels of the FAD2-2 and FAD2-3 genes under heat stress in strain 18N44 were lower than those in strain 18. These findings indicated that FAD2-2 and FAD2-3 may play major roles in the synthesis of linoleic acid during heat stress.

Chemical Protein Crosslinking-Coupled Mass Spectrometry Reveals Interaction of LHCI with LHCII and LHCSR3 in Chlamydomonas reinhardtii

The photosystem I (PSI) of the green alga Chlamydomonas reinhardtii associates with 10 light-harvesting proteins (LHCIs) to form the PSI-LHCI complex. In the context of state transitions, two LHCII trimers bind to the PSAL, PSAH and PSAO side of PSI to produce the PSI-LHCI-LHCII complex. In this work, we took advantage of chemical crosslinking of proteins in conjunction with mass spectrometry to identify protein-protein interactions between the light-harvesting proteins of PSI and PSII. We detected crosslinks suggesting the binding of LHCBM proteins to the LHCA1-PSAG side of PSI as well as protein-protein interactions of LHCSR3 with LHCA5 and LHCA3. Our data indicate that the binding of LHCII to PSI is more versatile than anticipated and imply that LHCSR3 might be involved in the regulation of excitation energy transfer to the PSI core via LHCA5/LHCA3.

GhCTSF1, a short PPR protein with a conserved role in chloroplast development and photosynthesis, participates in intron splicing of rpoC1 and ycf3-2 transcripts in cotton

Cotton is one of the most important textile fibers worldwide. As crucial agronomic traits, leaves play an essential role in the growth, disease resistance, fiber quality, and yield of cotton plants. Pentatricopeptide repeat (PPR) proteins are a large family of nuclear-encoded proteins involved in organellar or nuclear RNA metabolism. Using a virus-induced gene silencing assay, we found that cotton plants displayed variegated yellow leaf phenotypes with decreased chlorophyll content when expression of the PPR gene GhCTSF1 was silenced. GhCTSF1 encodes a chloroplast-localized protein that contains only two PPR motifs. Disruption of GhCTSF1 substantially reduces the splicing efficiency of rpoC1 intron 1 and ycf3 intron 2. Loss of function of the GhCTSF1 ortholog EMB1417 causes splicing defects in rpoC1 and ycf3-2, leading to impaired chloroplast structure and decreased photosynthetic rates in Arabidopsis. We also found that GhCTSF1 interacts with two splicing factors, GhCRS2 and GhWTF1. Defects in GhCRS2 and GhWTF1 severely affect intron splicing of rpoC1 and ycf3-2 in cotton, leading to defects in chloroplast development and a reduction in photosynthesis. Our results suggest that GhCTSF1 is specifically required for splicing rpoC1 and ycf3-2 in cooperation with GhCRS2 and GhWTF1.

An insulin-like signalling pathway model for Fasciola gigantica

BACKGROUND

The insulin/insulin-like signalling (IIS) pathway is common in mammals and invertebrates, and the IIS pathway is unknown in Fasciola gigantica. In the present study, the IIS pathway was reconstructed in F. gigantica. We defined the components involved in the IIS pathway and investigated the transcription profiles of these genes for all developmental stages of F. gigantica. In addition, the presence of these components in excretory and secretory products (ESPs) was predicted via signal peptide annotation.

RESULTS

The core components of the IIS pathway were detected in F. gigantica. Among these proteins, one ligand (FgILP) and one insulin-like molecule binding protein (FgIGFBP) were analysed. Interestingly, three receptors (FgIR-1/FgIR-2/FgIR-3) were detected, and a novel receptor, FgIR-3, was screened, suggesting novel functions. Fg14-3-3ζ, Fgirs, and Fgpp2a exhibited increased transcription in 42-day-old juveniles and 70-day-old juveniles, while Fgilp, Fgigfb, Fgsgk-1, Fgakt-1, Fgir-3, Fgpten, and Fgaap-1 exhibited increased transcription in metacercariae. FgILP, FgIGFBP, FgIR-2, FgIR-3, and two transcription factors (FgHSF-1 and FgSKN-1) were predicted to be present in FgESPs, indicating their exogenous roles.

CONCLUSIONS

This study helps to elucidate the signal transduction pathway of IIS in F. gigantica, which will aid in understanding the interaction between flukes and hosts, as well as in understanding fluke developmental regulation, and will also lay a foundation for further characterisation of the IIS pathways of trematodes.

Melatonin-Regulated Chaperone Binding Protein Plays a Key Role in Cadmium Stress Tolerance in Rice, Revealed by the Functional Characterization of a Novel Serotonin N-Acetyltransferase 3 (SNAT3) in Rice

The study of the mechanisms by which melatonin protects against cadmium (Cd) toxicity in plants is still in its infancy, particularly at the molecular level. In this study, the gene encoding a novel serotonin N-acetyltransferase 3 (SNAT3) in rice, a pivotal enzyme in the melatonin biosynthetic pathway, was cloned. Rice (Oryza sativa) OsSNAT3 is the first identified plant ortholog of archaeon Thermoplasma volcanium SNAT. The purified recombinant OsSNAT3 catalyzed the conversion of serotonin and 5-methoxytryptamine to N-acetylserotonin and melatonin, respectively. The suppression of OsSNAT3 by RNAi led to a decline in endogenous melatonin levels followed by a reduction in Cd tolerance in transgenic RNAi rice lines. In addition, the expression levels of genes encoding the endoplasmic reticulum (ER) chaperones BiP3, BiP4, and BiP5 were much lower in RNAi lines than in the wild type. In transgenic rice plants overexpressing OsSNAT3 (SNAT3-OE), however, melatonin levels were higher than in wild-type plants. SNAT3-OE plants also tolerated Cd stress, as indicated by seedling growth, malondialdehyde, and chlorophyll levels. BiP4 expression was much higher in the SNAT3-OE lines than in the wild type. These results indicate that melatonin engineering could help crops withstand Cd stress, resulting in high yields in Cd-contaminated fields.

Properties and Functional Analysis of Two Chorismate Mutases from Maritime Pine

Through the shikimate pathway, a massive metabolic flux connects the central carbon metabolism with the synthesis of chorismate, the common precursor of the aromatic amino acids phenylalanine, tyrosine, and tryptophan, as well as other compounds, including salicylate or folate. The alternative metabolic channeling of chorismate involves a key branch-point, finely regulated by aromatic amino acid levels. Chorismate mutase catalyzes the conversion of chorismate to prephenate, a precursor of phenylalanine and tyrosine and thus a vast repertoire of fundamental derived compounds, such as flavonoids or lignin. The regulation of this enzyme has been addressed in several plant species, but no study has included conifers or other gymnosperms, despite the importance of the phenolic metabolism for these plants in processes such as lignification and wood formation. Here, we show that maritime pine (Pinus pinaster Aiton) has two genes that encode for chorismate mutase, PpCM1 and PpCM2. Our investigations reveal that these genes encode plastidial isoenzymes displaying activities enhanced by tryptophan and repressed by phenylalanine and tyrosine. Using phylogenetic studies, we have provided new insights into the possible evolutionary origin of the cytosolic chorismate mutases in angiosperms involved in the synthesis of phenylalanine outside the plastid. Studies based on different platforms of gene expression and co-expression analysis have allowed us to propose that PpCM2 plays a central role in the phenylalanine synthesis pathway associated with lignification.

Tardigrade secretory proteins protect biological structures from desiccation

Tardigrades, microscopic animals that survive a broad range of environmental stresses, express a unique set of proteins termed tardigrade-specific intrinsically disordered proteins (TDPs). TDPs are often expressed at high levels in tardigrades upon desiccation, and appear to mediate stress adaptation. Here, we focus on the proteins belonging to the secreted family of tardigrade proteins termed secretory-abundant heat soluble ("SAHS") proteins, and investigate their ability to protect diverse biological structures. Recombinantly expressed SAHS proteins prevent desiccated liposomes from fusion, and enhance desiccation tolerance of E. coli and Rhizobium tropici upon extracellular application. Molecular dynamics simulation and comparative structural analysis suggest a model by which SAHS proteins may undergo a structural transition upon desiccation, in which removal of water and solutes from a large internal cavity in SAHS proteins destabilizes the beta-sheet structure. These results highlight the potential application of SAHS proteins as stabilizing molecules for preservation of cells.

Root Primordium Defective 1 Encodes an Essential PORR Protein Required for the Splicing of Mitochondria-Encoded Group II Introns and for Respiratory Complex I Biogenesis

Cellular respiration involves complex organellar metabolic activities that are pivotal for plant growth and development. Mitochondria contain their own genetic system (mitogenome, mtDNA), which encodes key elements of the respiratory machinery. Plant mtDNAs are notably larger than their counterparts in Animalia, with complex genome organization and gene expression characteristics. The maturation of the plant mitochondrial transcripts involves extensive RNA editing, trimming and splicing events. These essential processing steps rely on the activities of numerous nuclear-encoded cofactors, which may also play key regulatory roles in mitochondrial biogenesis and function and hence in plant physiology. Proteins that harbor the plant organelle RNA recognition (PORR) domain are represented in a small gene family in plants. Several PORR members, including WTF1, WTF9 and LEFKOTHEA, are known to act in the splicing of organellar group II introns in angiosperms. The AT4G33495 gene locus encodes an essential PORR protein in Arabidopsis, termed ROOT PRIMORDIUM DEFECTIVE 1 (RPD1). A null mutation of At.RPD1 causes arrest in early embryogenesis, while the missense mutant lines, rpd1.1 and rpd1.2, exhibit a strong impairment in root development and retarded growth phenotypes, especially under high-temperature conditions. Here, we further show that RPD1 functions in the splicing of introns that reside in the coding regions of various complex I (CI) subunits (i.e. nad2, nad4, nad5 and nad7), as well as in the maturation of the ribosomal rps3 pre-RNA in Arabidopsis mitochondria. The altered growth and developmental phenotypes and modified respiration activities are tightly correlated with respiratory chain CI defects in rpd1 mutants.

Multi-epitopevaccines, from design to expression; an in silico approach

The development of vaccines against a wide range of infectious diseases and pathogens often relies on multi-epitope strategies that can effectively stimulate both humoral and cellular immunity. Immunoinformatics tools play a pivotal role in designing such vaccines, enhancing immune response potential, and minimizing the risk of failure. This review presents a comprehensive overview of practical tools for epitope prediction and the associated immune responses. These immunoinformatics tools facilitate the selection of epitopes based on parameters such as antigenicity, absence of toxic and allergenic sequences, secondary and tertiary structures, sequence conservation, and population coverage. The chosen epitopes can be tailored for B-cells or T-cells, both of which require further assessments covered in this study. We offer a range of suitable linkers that effectively separate cytotoxic T lymphocyte and helper T lymphocyte epitopes while preserving their functionality. Additionally, we identify various adjuvants for specific purposes. We delve into the evaluation of MHC-epitope interactions, MHC clusters, and the simulation of final constructs through molecular docking techniques. We provide diverse linkers and adjuvants optimized for epitope functions to bolster immune responses through epitope attachment. By leveraging these comprehensive tools, the development of multi-epitope vaccines holds the promise of robust immunity and a significant reduction in experimental costs.

Combined sensory, volatilome and transcriptome analyses identify a limonene terpene synthase as a major contributor to the characteristic aroma of a Coffea arabica L. specialty coffee

BACKGROUND

The fruity aromatic bouquet of coffee has attracted recent interest to differentiate high value market produce as specialty coffee. Although the volatile compounds present in green and roasted coffee beans have been extensively described, no study has yet linked varietal molecular differences to the greater abundance of specific substances and support the aroma specificity of specialty coffees.

RESULTS

This study compared four Arabica genotypes including one, Geisha Especial, suggested to generate specialty coffee. Formal sensory evaluations of coffee beverages stressed the importance of coffee genotype in aroma perception and that Geisha Especial-made coffee stood out by having fine fruity, and floral, aromas and a more balanced acidity. Comparative SPME-GC-MS analyses of green and roasted bean volatile compounds indicated that those of Geisha Especial differed by having greater amounts of limonene and 3-methylbutanoic acid in agreement with the coffee cup aroma perception. A search for gene ontology differences of ripening beans transcriptomes of the four varieties revealed that they differed by metabolic processes linked to terpene biosynthesis due to the greater gene expression of prenyl-pyrophosphate biosynthetic genes and terpene synthases. Only one terpene synthase (CaTPS10-like) had an expression pattern that paralleled limonene loss during the final stage of berry ripening and limonene content in the studied four varieties beans. Its functional expression in tobacco leaves confirmed its functioning as a limonene synthase.

CONCLUSIONS

Taken together, these data indicate that coffee variety genotypic specificities may influence ripe berry chemotype and final coffee aroma unicity. For the specialty coffee variety Geisha Especial, greater expression of terpene biosynthetic genes including CaTPS10-like, a limonene synthase, resulted in the greater abundance of limonene in green beans, roasted beans and a unique citrus note of the coffee drink.

Codon usage bias of secretory protein in Fusarium oxysporum f. sp. cubense tropical race 4

Banana Fusarium oxysporum f. sp. cubense tropical race 4 (Foc-TR4) is a highly destructive pathogen that infects nearly all major banana cultivars and has a tendency to spread further. Secreted proteins play a crucial role in the process of Fusarium wilt infection in bananas. In this study, we analyzed the codon usage bias (CUB) of the Foc-TR4 classical secretory protein genome for the first time and observed a strong bias toward codons ending with C. We found that 572 out of the 14,543 amino acid sequences in the Foc-TR4 genome exhibited characteristics of classical secretory proteins. The CUB was largely influenced by selection optimization pressure, as indicated by the ENC value and neutral plot analysis. Among the identified codons, such as UCC and CCC, 11 were found to be optimal for Foc-TR4 gene expression. Codons with higher GC content and a C base in the third position showed greater selectivity. The CUB in the secretory proteins encoded by Foc-TR4 provides insights into their evolutionary patterns, contributing to the development and screening of novel and effective antifungal drugs.

Discovery of Three Bemisia tabaci Effectors and Their Effect on Gene Expression in Planta

Bemisia tabaci (whitefly) is a polyphagous agroeconomic pest species complex. Two members of this species complex, Mediterranean (MED) and Middle-East-Asia Minor 1 (MEAM1), have a worldwide distribution and have been shown to manipulate plant defenses through effectors. In this study, we used three different strategies to identify three MEAM1 proteins that can act as effectors. Effector B1 was identified using a bioinformatics-driven effector-mining strategy, whereas effectors S1 and P1 were identified in the saliva of whiteflies collected from artificial diet and in phloem exudate of tomato on which nymphs were feeding, respectively. These three effectors were B. tabaci specific and able to increase whitefly fecundity when transiently expressed in tobacco plants (Nicotiana tabacum). Moreover, they reduced growth of Pseudomonas syringae pv. tabaci in Nicotiana benthamiana. All three effectors changed gene expression in planta, and B1 and S1 also changed phytohormone levels. Gene ontology and KEGG pathway enrichment analysis pinpointed plant-pathogen interaction and photosynthesis as the main enriched pathways for all three effectors. Our data thus show the discovery and validation of three new B. tabaci MEAM1 effectors that increase whitefly fecundity and modulate plant immunity. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Mutation of negative regulatory gene CEHC1 encoding an FBXO3 protein results in normoxic expression of HYDA genes in Chlamydomonas reinhardtii

Oxygen is known to prevent hydrogen production in Chlamydomonas, both by inhibiting the hydrogenase enzyme and by preventing the accumulation of HYDA-encoding transcripts. We developed a screen for mutants showing constitutive accumulation of HYDA1 transcripts in the presence of oxygen. A reporter gene required for ciliary motility, placed under the control of the HYDA1 promoter, conferred motility only in hypoxic conditions. By selecting for mutants able to swim even in the presence of oxygen we obtained strains that express the reporter gene constitutively. One mutant identified a gene encoding an F-box only protein 3 (FBXO3), known to participate in ubiquitylation and proteasomal degradation pathways in other eukaryotes. Transcriptome profiles revealed that the mutation, termed cehc1-1 , leads to constitutive expression of HYDA1 and other genes regulated by hypoxia, and of many genes known to be targets of CRR1, a transcription factor in the nutritional copper signaling pathway. CRR1 was required for the constitutive expression of the HYDA1 reporter gene in cehc1-1 mutants. The CRR1 protein, which is normally degraded in Cu-supplemented cells, was stabilized in cehc1-1 cells, supporting the conclusion that CEHC1 acts to facilitate the degradation of CRR1. Our results reveal a novel negative regulator in the CRR1 pathway and possibly other pathways leading to complex metabolic changes associated with response to hypoxia.

PPR596 Is Required for nad2 Intron Splicing and Complex I Biogenesis in Arabidopsis

Mitochondria are essential organelles that generate energy via oxidative phosphorylation. Plant mitochondrial genome encodes some of the respiratory complex subunits, and these transcripts require accurate processing, including C-to-U RNA editing and intron splicing. Pentatricopeptide repeats (PPR) proteins are involved in various organellar RNA processing events. PPR596, a P-type PPR protein, was previously identified to function in the C-to-U editing of mitochondrial rps3 transcripts in Arabidopsis. Here, we demonstrate that PPR596 functions in the cis-splicing of nad2 intron 3 in mitochondria. Loss of the PPR596 function affects the editing at rps3eU1344SS, impairs nad2 intron 3 splicing and reduces the mitochondrial complex I's assembly and activity, while inducing alternative oxidase (AOX) gene expression. This defect in nad2 intron splicing provides a plausible explanation for the slow growth of the ppr595 mutants. Although a few P-type PPR proteins are involved in RNA C-to-U editing, our results suggest that the primary function of PPR596 is intron splicing.

Targeting imidazole-glycerol phosphate dehydratase in plants: novel approach for structural and functional studies, and inhibitor blueprinting

The histidine biosynthetic pathway (HBP) is targeted for herbicide design with preliminary success only regarding imidazole-glycerol phosphate dehydratase (IGPD, EC 4.2.1.19), or HISN5, as referred to in plants. HISN5 catalyzes the sixth step of the HBP, in which imidazole-glycerol phosphate (IGP) is dehydrated to imidazole-acetol phosphate. In this work, we present high-resolution cryoEM and crystal structures of Medicago truncatula HISN5 (MtHISN5) in complexes with an inactive IGP diastereoisomer and with various other ligands. MtHISN5 can serve as a new model for plant HISN5 structural studies, as it enables resolving protein-ligand interactions at high (2.2 Å) resolution using cryoEM. We identified ligand-binding hotspots and characterized the features of plant HISN5 enzymes in the context of the HISN5-targeted inhibitor design. Virtual screening performed against millions of small molecules not only revealed candidate molecules but also identified linkers for fragments that were experimentally confirmed to bind. Based on experimental and computational approaches, this study provides guidelines for designing symmetric HISN5 inhibitors that can reach two neighboring active sites. Finally, we conducted analyses of sequence similarity networks revealing that plant HISN5 enzymes derive from cyanobacteria. We also adopted a new approach to measure MtHISN5 enzymatic activity using isothermal titration calorimetry and enzymatically synthesized IGP.

Genome-wide identification, stress- and hormone-responsive expression characteristics, and regulatory pattern analysis of Scutellaria baicalensis SbSPLs

SQUAMOSA PROMOTER BINDING PROTEIN-LIKEs (SPLs) encode plant-specific transcription factors that regulate plant growth and development, stress response, and metabolite accumulation. However, there is limited information on Scutellaria baicalensis SPLs. In this study, 14 SbSPLs were identified and divided into 8 groups based on phylogenetic relationships. SbSPLs in the same group had similar structures. Abscisic acid-responsive (ABRE) and MYB binding site (MBS) cis-acting elements were found in the promoters of 8 and 6 SbSPLs. Segmental duplications and transposable duplications were the main causes of SbSPL expansion. Expression analysis based on transcriptional profiling showed that SbSPL1, SbSPL10, and SbSPL13 were highly expressed in roots, stems, and flowers, respectively. Expression analysis based on quantitative real-time polymerase chain reaction (RT‒qPCR) showed that most SbSPLs responded to low temperature, drought, abscisic acid (ABA) and salicylic acid (SA), among which the expression levels of SbSPL7/9/10/12 were significantly upregulated in response to abiotic stress. These results indicate that SbSPLs are involved in the growth, development and stress response of S. baicalensis. In addition, 8 Sba-miR156/157 s were identified, and SbSPL1-5 was a potential target of Sba-miR156/157 s. The results of target gene prediction and coexpression analysis together indicated that SbSPLs may be involved in the regulation of L-phenylalanine (L-Phe), lignin and jasmonic acid (JA) biosynthesis. In summary, the identification and characterization of the SbSPL gene family lays the foundation for functional research and provides a reference for improved breeding of S. baicalensis stress resistance and quality traits.

Assessing sequence variation, haplotype analysis and molecular characterisation of aspartate kinase2 (ask2) gene regulating methionine biosynthesis in diverse maize inbreds

Traditional maize grain is deficient in methionine, an essential amino acid required for proper growth and development in humans and poultry birds. Thus, development of high methionine maize (HMM) assumes great significance in alleviating malnutrition through sustainable and cost-effective approach. Of various genetic loci, aspartate kinase2 (ask2) gene plays a pivotal role in regulating methionine accumulation in maize. Here, we sequenced the entire ask2 gene of 5394 bp with 13 exons in five wild and five mutant maize inbreds to understand variation at nucleotide level. Sequence analysis revealed that an SNP in exon-13 caused thymine to adenine transversion giving rise to a favourable mutant allele associated with leucine to glutamine substitution in mutant ASK2 protein. Gene-based diversity analysis with 11 InDel markers grouped 48 diverse inbreds into three major clusters with an average genetic dissimilarity of 0.570 (range, 0.0-0.9). The average major allele frequency, gene diversity and PIC are 0.693, 0.408 and 0.341, respectively. A total of 45 haplotypes of the ask2 gene were identified among the maize inbreds. Evolutionary relationship analysis performed among 22 orthologues grouped them into five major clusters. The number of exons varied from 7 to 17, with length varying from 12 to 495 bp among orthologues. ASK2 protein with 565 amino acids was predicted to be in homo-dimeric state with lysine and tartaric acid as binding ligands. Amino acid kinase and ACT domains were found to be conserved in maize and orthologues. The study depicted the presence of enough genetic diversity in ask2 gene in maize, and development of HMM can be accelerated through introgression of favourable allele of ask2 into the parental lines of elite hybrids using molecular breeding.

The PPO family in Nicotiana tabacum is an important regulator to participate in pollination

Polyphenol oxidases (PPOs) are type-3 copper enzymes and are involved in many biological processes. However, the potential functions of PPOs in pollination are not fully understood. In this work, we have screened 13 PPO members in Nicotiana. tabacum (named NtPPO1-13, NtPPOs) to explore their characteristics and functions in pollination. The results show that NtPPOs are closely related to PPOs in Solanaceae and share conserved domains except NtPPO4. Generally, NtPPOs are diversely expressed in different tissues and are distributed in pistil and male gametes. Specifically, NtPPO9 and NtPPO10 are highly expressed in the pistil and mature anther. In addition, the expression levels and enzyme activities of NtPPOs are increased after N. tabacum self-pollination. Knockdown of NtPPOs would affect pollen growth after pollination, and the purines and flavonoid compounds are accumulated in self-pollinated pistil. Altogether, our findings demonstrate that NtPPOs potentially play a role in the pollen tube growth after pollination through purines and flavonoid compounds, and will provide new insights into the role of PPOs in plant reproduction.

Encyclopedia of Family A DNA Polymerases Localized in Organelles: Evolutionary Contribution of Bacteria Including the Proto-Mitochondrion

DNA polymerases synthesize DNA from deoxyribonucleotides in a semiconservative manner and serve as the core of DNA replication and repair machinery. In eukaryotic cells, there are 2 genome-containing organelles, mitochondria, and plastids, which were derived from an alphaproteobacterium and a cyanobacterium, respectively. Except for rare cases of genome-lacking mitochondria and plastids, both organelles must be served by nucleus-encoded DNA polymerases that localize and work in them to maintain their genomes. The evolution of organellar DNA polymerases has yet to be fully understood because of 2 unsettled issues. First, the diversity of organellar DNA polymerases has not been elucidated in the full spectrum of eukaryotes. Second, it is unclear when the DNA polymerases that were used originally in the endosymbiotic bacteria giving rise to mitochondria and plastids were discarded, as the organellar DNA polymerases known to date show no phylogenetic affinity to those of the extant alphaproteobacteria or cyanobacteria. In this study, we identified from diverse eukaryotes 134 family A DNA polymerase sequences, which were classified into 10 novel types, and explored their evolutionary origins. The subcellular localizations of selected DNA polymerases were further examined experimentally. The results presented here suggest that the diversity of organellar DNA polymerases has been shaped by multiple transfers of the PolI gene from phylogenetically broad bacteria, and their occurrence in eukaryotes was additionally impacted by secondary plastid endosymbioses. Finally, we propose that the last eukaryotic common ancestor may have possessed 2 mitochondrial DNA polymerases, POP, and a candidate of the direct descendant of the proto-mitochondrial DNA polymerase I, rdxPolA, identified in this study.

Genome-wide identification of heat shock protein gene family and their responses to pathogen challenge in Trachinotus ovatus

Heat Shock Proteins (HSPs) are a widely distributed family of proteins produced in response to heat and other stresses. To develop a deeper understanding of the mechanisms governing expression of HSPs in the bony fish Trachinotus ovatus, we carried out a whole genome analysis and identified 43 HSP genes. Based on their phylogenetic relationships with Danio rerio, Seriola dumerili, and Seriola lalandi, they were divided into four subfamilies: HSP20, HSP60, HSP70, and HSP90. We performed an analysis of the predicted physicochemical properties and subcellular localization of proteins encoded by these genes. The chromosomal localization results showed that the HSP genes are distributed across 20 chromosomes of T. ovatus.These genes were found to be expressed in different tissues, and they showed differential expression in the immune response against Streptococcus agalactiae. However, there was no significant differential expression in the different skin tissue locations of T. ovatus after infection by Cryptocaryon irritans Brown. This study provides basic information for further research on the evolution and structure and function of HSPs in teleosts.

Multi-omics analyses reveal the importance of chromoplast plastoglobules in carotenoid accumulation in citrus fruit

Chromoplasts act as a metabolic sink for carotenoids, in which plastoglobules serve as versatile lipoprotein particles. PGs in chloroplasts have been characterized. However, the features of PGs from non-photosynthetic plastids are poorly understood. We found that the development of chromoplast plastoglobules (CPGs) in globular and crystalloid chromoplasts of citrus is associated with alterations in carotenoid storage. Using Nycodenz density gradient ultracentrifugation, an efficient protocol for isolating highly purified CPGs from sweet orange (Citrus sinensis) pulp was established. Forty-four proteins were defined as likely comprise the core proteome of CPGs using comparative proteomics analysis. Lipidome analysis of different chromoplast microcompartments revealed that the nonpolar microenvironment within CPGs was modified by 35 triacylglycerides, two sitosterol esters, and one stigmasterol ester. Manipulation of the CPG-localized gene CsELT1 (esterase/lipase/thioesterase) in citrus calli resulted in increased lipids and carotenoids, which is further evidence that the nonpolar microenvironment of CPGs contributes to carotenoid accumulation and storage in the chromoplasts. This multi-feature analysis of CPGs sheds new light on the role of chromoplasts in carotenoid metabolism, paving the way for manipulating carotenoid content in citrus fruit and other crops.

Association Analysis Provides Insights into Plant Mitonuclear Interactions

Cytonuclear interaction refers to the complex and ongoing process of coevolution between nuclear and organelle genomes, which are responsible for cellular respiration, photosynthesis, lipid metabolism, etc. and play a significant role in adaptation and speciation. There have been a large number of studies to detect signatures of cytonuclear interactions. However, identification of the specific nuclear and organelle genetic polymorphisms that are involved in these interactions within a species remains relatively rare. The recent surge in whole genome sequencing has provided us an opportunity to explore cytonuclear interaction from a population perspective. In this study, we analyzed a total of 3,439 genomes from 7 species to identify signals of cytonuclear interactions by association (linkage disequilibrium) analysis of variants in both the mitochondrial and nuclear genomes across flowering plants. We also investigated examples of nuclear loci identified based on these association signals using subcellular localization assays, gene editing, and transcriptome sequencing. Our study provides a novel perspective on the investigation of cytonuclear coevolution, thereby enriching our understanding of plant fitness and offspring sterility.

Diurnal Rhythms in the Red Seaweed Gracilariopsis chorda are Characterized by Unique Regulatory Networks of Carbon Metabolism

Cellular and physiological cycles are driven by endogenous pacemakers, the diurnal and circadian rhythms. Key functions such as cell cycle progression and cellular metabolism are under rhythmic regulation, thereby maintaining physiological homeostasis. The photoreceptors phytochrome and cryptochrome, in response to light cues, are central input pathways for physiological cycles in most photosynthetic organisms. However, among Archaeplastida, red algae are the only taxa that lack phytochromes. Current knowledge about oscillatory rhythms is primarily derived from model species such as Arabidopsis thaliana and Chlamydomonas reinhardtii in the Viridiplantae, whereas little is known about these processes in other clades of the Archaeplastida, such as the red algae (Rhodophyta). We used genome-wide expression profiling of the red seaweed Gracilariopsis chorda and identified 3,098 rhythmic genes. Here, we characterized possible cryptochrome-based regulation and photosynthetic/cytosolic carbon metabolism in this species. We found a large family of cryptochrome genes in G. chorda that display rhythmic expression over the diurnal cycle and may compensate for the lack of phytochromes in this species. The input pathway gates regulatory networks of carbon metabolism which results in a compact and efficient energy metabolism during daylight hours. The system in G. chorda is distinct from energy metabolism in most plants, which activates in the dark. The green lineage, in particular, land plants, balance water loss and CO2 capture in terrestrial environments. In contrast, red seaweeds maintain a reduced set of photoreceptors and a compact cytosolic carbon metabolism to thrive in the harsh abiotic conditions typical of intertidal zones.

Distinct structural motifs are necessary for targeting and import of Tim17 in Trypanosoma brucei mitochondrion

Nuclear-encoded mitochondrial proteins are correctly translocated to their proper sub-mitochondrial destination using location-specific mitochondrial targeting signals and via multi-protein import machineries (translocases) in the outer and inner mitochondrial membranes (TOM and TIMs, respectively). However, targeting signals of multi-pass Tims are less defined. Here, we report the characterization of the targeting signals of Trypanosoma brucei Tim17 (TbTim17), an essential component of the most divergent TIM complex. TbTim17 possesses a characteristic secondary structure including four predicted transmembrane (TM) domains in the center with hydrophilic N- and C-termini. After examining mitochondrial localization of various deletion and site-directed mutants of TbTim17 in T. brucei using subcellular fractionation and confocal microscopy, we located at least two internal targeting signals (ITS): (i) within TM1 (31-50 AAs) and (ii) TM4 + loop 3 (120-136 AAs). Both signals are required for proper targeting and integration of TbTim17 in the membrane. Furthermore, a positively charged residue (K122) is critical for mitochondrial localization of TbTim17. This is the first report of characterizing the ITS for a multipass inner membrane protein in a divergent eukaryote, like T. brucei.IMPORTANCEAfrican trypanosomiasis (AT) is a deadly disease in human and domestic animals, caused by the parasitic protozoan Trypanosoma brucei. Therefore, AT is not only a concern for human health but also for economic development in the vast area of sub-Saharan Africa. T. brucei possesses a single mitochondrion per cell that imports hundreds of nuclear-encoded mitochondrial proteins for its functions. T. brucei Tim17 (TbTim17), an essential component of the TbTIM17 complex, is a nuclear-encoded protein; thus, it is necessary to be imported from the cytosol to form the TbTIM17 complex. Here, we demonstrated that the internal targeting signals within the transmembrane 1 (TM1) and TM4 with loop 3, and residue K122 are required collectively for import and integration of TbTim17 in the T. brucei mitochondrion. This information could be utilized to block TbTim17 function and parasite growth.

Transcriptomic profiling of different developmental stages reveals parasitic strategies of Wohlfahrtia magnifica, a myiasis-causing flesh fly

BACKGROUND

Wohlfahrtia magnifica is an obligatory parasite that causes myiasis in several warm-blooded vertebrates. Adult females deposit the first-stage larvae directly onto wounds or natural body orifices (e.g., genitalia) of the host, from where they quickly colonize the host tissue and feed on it for development. The infestation of W. magnifica can lead to health issues, welfare concerns, and substantial economic losses. To date, little is known about the molecular mechanisms of the W. magnifica-causing myiasis.

RESULTS

In this study, we collected parasitic-stage larvae of W. magnifica from wounds of naturally infested Bactrian camels, as well as pupae and adult flies reared in vitro from the wound-collected larvae, for investigating the gene expression profiles of the different developmental stages of W. magnifica, with a particular focus on examining gene families closely related to the parasitism of the wound-collected larvae. As key proteins related to the parasite-host interaction, 2049 excretory/secretory (ES) proteins were identified in W. magnifica through the integration of multiple bioinformatics approaches. Functional analysis indicates that these ES proteins are primarily involved in cuticle development, peptidase activity, immune response, and metabolic processes. The global investigation of gene expression at different developmental stages using pairwise comparisons and weighted correlation network analysis (WGCNA) showed that the upregulated genes during second-stage larvae were related to cuticle development, peptidase activity, and RNA transcription and translation; during third-stage larvae to peptidase inhibitor activity and nutrient reservoir activity; during pupae to cell and tissue morphogenesis and cell and tissue development; and during adult flies to signal perception, many of them involved in light perception, and adult behavior, e.g., feeding, mating, and locomotion. Specifically, the expression level analysis of the likely parasitism-related genes in parasitic wound-collected larvae revealed a significant upregulation of 88 peptidase genes (including 47 serine peptidase genes), 110 cuticle protein genes, and 21 heat shock protein (hsp) genes. Interestingly, the expression of 2 antimicrobial peptide (AMP) genes, including 1 defensin and 1 diptericin, was also upregulated in the parasitic larvae.

CONCLUSIONS

We identified ES proteins in W. magnifica and investigated their functional distribution. In addition, gene expression profiles at different developmental stages of W. magnifica were examined. Specifically, we focused on gene families closely related to parasitism of wound-collected larvae. These findings shed light on the molecular mechanisms underlying the life cycle of the myiasis-causing fly, especially during the parasitic larval stages, and provide guidance for the development of control measures against W. magnifica.

Expression analysis of genes encoding extracellular leucine-rich repeat proteins in Arabidopsis thaliana

Leucine-rich repeat (LRR)-containing proteins have been identified in diverse species, including plants. The diverse intracellular and extracellular LRR variants are responsible for numerous biological processes. We analyzed the expression patterns of Arabidopsis thaliana extracellular LRR (AtExLRR) genes, 10 receptor-like proteins, and 4 additional genes expressing the LRR-containing protein by a promoter: β-glucuronidase (GUS) study. According to in silico expression studies, several AtExLRR genes were expressed in a tissue- or stage-specific and abiotic/hormone stress-responsive manner, indicating their potential participation in specific biological processes. Based on the promoter: GUS assay, AtExLRRs were expressed in different cells and organs. A quantitative real-time PCR investigation revealed that the expressions of AtExLRR3 and AtExLRR9 were distinct under various abiotic stress conditions. This study investigated the potential roles of extracellular LRR proteins in plant growth, development, and response to various abiotic stresses.

Improving Signal and Transit Peptide Predictions Using AlphaFold2-predicted Protein Structures

Many proteins contain cleavable signal or transit peptides that direct them to their final subcellular locations. Such peptides are usually predicted from sequence alone using methods such as TargetP 2.0 and SignalP 6.0. While these methods are usually very accurate, we show here that an analysis of a protein's AlphaFold2-predicted structure can often be used to identify false positive predictions. We start by showing that when given a protein's full-length sequence, AlphaFold2 builds experimentally annotated signal and transit peptides in orientations that point away from the main body of the protein. This indicates that AlphaFold2 correctly identifies that a signal is not destined to be part of the mature protein's structure and suggests, as a corollary, that predicted signals that AlphaFold2 folds with high confidence into the main body of the protein are likely to be false positives. To explore this idea, we analyzed predicted signal peptides in 48 proteomes made available in DeepMind's AlphaFold2 database (https://alphafold.ebi.ac.uk). Applying TargetP 2.0 and SignalP 6.0 to the 561,562 proteins in the database results in 95,236 being predicted to contain a cleavable signal or transit peptide. In 95.1% of these cases, the AlphaFold2 structure of the full-length protein is fully consistent with the prediction of TargetP 2.0 or SignalP 6.0. In the remaining 4.9% of cases where the AlphaFold2 structure does not appear consistent with the prediction, the signal is often only predicted with low confidence. The potential false positives identified here may be useful for training even more accurate signal prediction methods.

Novel Plastid Genome Characteristics in Fugacium kawagutii and the Trend of Accelerated Evolution of Plastid Proteins in Dinoflagellates

Typical (peridinin-containing) dinoflagellates possess plastid genomes composed of small plasmids named "minicircles". Despite the ecological importance of dinoflagellate photosynthesis in corals and marine ecosystems, the structural characteristics, replication dynamics, and evolutionary forcing of dinoflagellate plastid genomes remain poorly understood. Here, we sequenced the plastid genome of the symbiodiniacean species Fugacium kawagutii and conducted comparative analyses. We identified psbT-coding minicircles, features previously not found in Symbiodiniaceae. The copy number of F. kawagutii minicircles showed a strong diel dynamics, changing between 3.89 and 34.3 copies/cell and peaking in mid-light period. We found that F. kawagutii minicircles are the shortest among all dinoflagellates examined to date. Besides, the core regions of the minicircles are highly conserved within genus in Symbiodiniaceae. Furthermore, the codon usage bias of the plastid genomes in Heterocapsaceae, Amphidiniaceae, and Prorocentraceae species are greatly influenced by selection pressure, and in Pyrocystaceae, Symbiodiniaceae, Peridiniaceae, and Ceratiaceae species are influenced by both natural selection pressure and mutation pressure, indicating a family-level distinction in codon usage evolution in dinoflagellates. Phylogenetic analysis using 12 plastid-encoded proteins and five nucleus-encoded plastid proteins revealed accelerated evolution trend of both plastid- and nucleus-encoded plastid proteins in peridinin- and fucoxanthin-dinoflagellate plastids compared to plastid proteins of nondinoflagellate algae. These findings shed new light on the structure and evolution of plastid genomes in dinoflagellates, which will facilitate further studies on the evolutionary forcing and function of the diverse dinoflagellate plastids. The accelerated evolution documented here suggests plastid-encoded sequences are potentially useful for resolving closely related dinoflagellates.

The uterine secretory cycle: recurring physiology of endometrial outputs that setup the uterine luminal microenvironment

Conserved in female reproduction across all mammalian species is the estrous cycle and its regulation by the hypothalamic-pituitary-gonadal (HPG) axis, a collective of intersected hormonal events that are crucial for ensuring uterine fertility. Nonetheless, knowledge of the direct mediators that synchronously shape the uterine microenvironment for successive yet distinct events, such as the transit of sperm and support for progressive stages of preimplantation embryo development, remain principally deficient. Toward understanding the timed endometrial outputs that permit luminal events as directed by the estrous cycle, we used Bovidae as a model system to uniquely surface sample and study temporal shifts to in vivo endometrial transcripts that encode for proteins destined to be secreted. The results revealed the full quantitative profile of endometrial components that shape the uterine luminal microenvironment at distinct phases of the estrous cycle (estrus, metestrus, diestrus, and proestrus). In interpreting this comprehensive log of stage-specific endometrial secretions, we define the "uterine secretory cycle" and extract a predictive understanding of recurring physiological actions regulated within the uterine lumen in anticipation of sperm and preimplantation embryonic stages. This repetitive microenvironmental preparedness to sequentially provide operative support was a stable intrinsic framework, with only limited responses to sperm or embryos if encountered in the lumen within the cyclic time period. In uncovering the secretory cycle and unraveling realistic biological processes, we present novel foundational knowledge of terminal effectors controlled by the HPG axis to direct a recurring sequence of vital functions within the uterine lumen.NEW & NOTEWORTHY This study unravels the recurring sequence of changes within the uterus that supports vital functions (sperm transit and development of preimplantation embryonic stages) during the reproductive cycle in female Ruminantia. These data present new systems knowledge in uterine reproductive physiology crucial for setting up in vitro biomimicry and artificial environments for assisted reproduction technologies for a range of mammalian species.

Reconstruction of the TGF-β signaling pathway of Fasciola gigantica

In the present study, we reconstructed the transforming growth factor beta (TGF-β) signaling pathway for Fasciola gigantica, which is a neglected tropical pathogen. We defined the components involved in the TGF-β signaling pathway and investigated the transcription profiles of these genes for all developmental stages of F. gigantica. In addition, the presence of these components in excretory and secretory products (FgESP) was predicted via signal peptide annotation. The core components of the TGF-β signaling pathway have been detected in F. gigantica; classical and nonclassical single transduction pathways were constructed. Four ligands have been detected, which may mediate the TGF-β signaling pathway and BMP signaling pathway. Two ligand-binding type II receptors were detected, and inhibitory Smad7 was not detected. TLP, BMP-3, BMP-1, and ActRIb showed higher transcription in 42-day juvenile and 70-day juvenile, while ActRIIa, Smad1, ActRIIb, Smad8, KAT2B, and PP2A showed higher transcription in egg. TLM, Ski, Smad6, BMPRI, p70S6K, Smad2, Smad3, TgfβRI, Smad4, and p300 showed higher transcription in metacercariae. Four ligands, 2 receptors and 3 Smads are predicted to be present in the FgESP, suggesting their potential extrinsic function. This study should help to understand signal transduction in the TGF-β signaling pathway in F. gigantica. In addition, this study helps to illustrate the complex mechanisms involved in developmental processes and F. gigantica - host interaction and paves the way for further characterization of the signaling pathway in trematodes.

Engineering of CYP82Y1, a cytochrome P450 monooxygenase: a key enzyme in noscapine biosynthesis in opium poppy

Protein engineering provides a powerful base for the circumvention of challenges tied with characteristics accountable for enzyme functions. CYP82Y1 introduces a hydroxyl group (-OH) into C1 of N-methylcanadine as the substrate to yield 1-hydroxy-N-methylcanadine. This chemical process has been found to be the gateway to noscapine biosynthesis. Owning to the importance of CYP82Y1 in this biosynthetic pathway, it has been selected as a target for enzyme engineering. The insertion of tags to the N- and C-terminal of CYP82Y1 was assessed for their efficiencies for improvement of the physiological performances of CYP82Y1. Although these attempts achieved some positive results, further strategies are required to dramatically enhance the CYP82Y1 activity. Here methods that have been adopted to achieve a functionally improved CYP82Y1 will be reviewed. In addition, the possibility of recruitment of other techniques having not yet been implemented in CYP82Y1 engineering, including the substitution of the residues located in the substrate recognition site, formation of the synthetic fusion proteins, and construction of the artificial lipid-based scaffold will be discussed. Given the fact that the pace of noscapine synthesis is constrained by the CYP82Y1-catalyzing step, the methods proposed here are capable of accelerating the rate of reaction performed by CYP82Y1 through improving its properties, resulting in the enhancement of noscapine accumulation.

Gleaning Euglenozoa-specific DNA polymerases in public single-cell transcriptome data

Multiple genes encoding family A DNA polymerases (famA DNAPs), which are evolutionary relatives of DNA polymerase I (PolI) in bacteria and phages, have been found in eukaryotic genomes, and many of these proteins are used mainly in organelles. Among members of the phylum Euglenozoa, distinct types of famA DNAP, PolIA, PolIBCD+, POP, and eugPolA, have been found. It is intriguing how the suite of famA DNAPs had been established during the evolution of Euglenozoa, but the DNAP data have not been sampled from the taxa that sufficiently represent the diversity of this phylum. In particular, little sequence data were available for basal branching species in Euglenozoa until recently. Thanks to the single-cell transcriptome data from symbiontids and phagotrophic euglenids, we have an opportunity to cover the "hole" in the repertory of famA DNAPs in the deep branches in Euglenozoa. The current study identified 16 new famA DNAP sequences in the transcriptome data from 33 phagotrophic euglenids and two symbiontids, respectively. Based on the new famA DNAP sequences, the updated diversity and evolution of famA DNAPs in Euglenozoa are discussed.

Genome-wide survey and evolutionary history of the pectin methylesterase (PME) gene family in the Dothideomycetes class of fungi

Once deposited in the plant cell wall, pectin undergoes demethylesterification by endogenous pectin methylesterases (PMEs), which play various roles in growth and development, including defense against pathogen attacks. Pathogen PMEs can alter pectin's methylesterification pattern, increasing its susceptibility to degradation by other fungal pectinases and thus playing a critical role as virulence factors during early infection stages. To investigate the evolutionary history of PMEs in the Dothideomycetes class of fungi, we obtained genomic data from 15 orders (79 species) and added genomic data from 61 isolates of Corynespora cassiicola. Our analyses involved maximum likelihood phylogenies, gene genealogies, and selection analyses. Additionally, we measured PME gene expression levels of C. cassiicola using soybean as a host through RT-qPCR assays. We recovered 145 putative effector PMEs and 57 putative non-effector PMEs from across the Dothideomycetes. The PME gene family exhibits a small size (up to 5 members per genome) and comprises three major clades. The evolutionary patterns of the PME1 and PME2 clades were largely shaped by duplications and recurring gene retention events, while biased gene loss characterized the small-sized PME3 clade. The presence of five members in the PME gene family of C. cassiicola suggests that the family may play a key role in the evolutionary success of C. cassiicola as a polyphagous plant pathogen. The haplogroups Cc_PME1.1 and Cc_PME1.2 exhibited an accelerated rate of evolution, whereas Cc_PME2.1, Cc_PME2.2, and Cc_PME2.3 seem to be under strong purifying selective constraints. All five PME genes were expressed during infection of soybean leaves, with the highest levels during from six to eight days post-inoculation. The highest relative expression level was measured for CC_29_g7533, a member of the Cc_PME2.3 clade, while the remaining four genes had relatively lower levels of expression.