SMART: recent updates, new developments and status in 2020

Systematic Analysis of bHLH Transcription Factors in Cassava Uncovers Their Roles in Postharvest Physiological Deterioration and Cyanogenic Glycosides Biosynthesis

The basic helix-loop-helix (bHLH) proteins are a large superfamily of transcription factors, and play a central role in a wide range of metabolic, physiological, and developmental processes in higher organisms. However, systematic investigation of bHLH gene family in cassava (Manihot esculenta Crantz) has not been reported. In the present study, we performed a genome-wide survey and identified 148 MebHLHs genes were unevenly harbored in 18 chromosomes. Through phylogenetic analyses along with Arabidopsis counterparts, these MebHLHs genes were divided into 19 groups, and each gene contains a similar structure and conserved motifs. Moreover, many cis-acting regulatory elements related to various defense and stress responses showed in MebHLH genes. Interestingly, transcriptome data analyses unveiled 117 MebHLH genes during postharvest physiological deterioration (PPD) process of cassava tuberous roots, while 65 MebHLH genes showed significantly change. Meanwhile, the relative quantitative analysis of 15 MebHLH genes demonstrated that they were sensitive to PPD, suggesting they may involve in PPD process regulation. Cyanogenic glucosides (CGs) biosynthesis during PPD process was increased, silencing of MebHLH72 and MebHLH114 showed that linamarin content was significantly decreased in the leaves. To summarize, the genome-wide identification and expression profiling of MebHLH candidates pave a new avenue for uderstanding their function in PPD and CGs biosynthesis, which will accelerate the improvement of PPD tolerance and decrease CGs content in cassava tuberous roots.

Bioinformatics Approaches for Analyzing Multigene Families Encoding Immune Receptors

Genome sequences are quickly becoming available from a variety of organisms, providing researchers with an abundance of previously inaccessible information and an important source of insight into immune mechanisms. There are a variety of methods to accurately characterize genes from new genome sequences, but immune receptors pose special challenges for these techniques. Immune receptors, particularly those that directly recognize pathogens, often diverge rapidly among species and are commonly found in large, complex multigene families. Because of these characteristics, immune receptors tend to be overlooked or misannotated in large-scale genomic surveys. We describe here a strategy to characterize homologs of immune receptors and to identify putative receptors from newly assembled genome or transcriptome sequences. The description of these protocols is aimed at a typical immunologist and does not rely on substantial a priori knowledge of bioinformatics. The approach is based on using low-stringency sequence searches to identify divergent homologs. For receptors with multiple domains, the intersection of low-stringency searches can be used to identify divergent receptor sequences with high confidence. For multigene families, these predictions can be refined using sequence conservation among gene family paralogs. Assembled genome sequences serve as a critical foundation for subsequent functional characterization and remove long-standing barriers in understanding the evolution of immune recognition systems.

Prediction of Potential Drug Targets and Vaccine Candidates Against Antibiotic-Resistant Pseudomonas aeruginosa

Pseudomonas aeruginosa is one of the leading causes of nosocomial infections, characterized by increasing antibiotic resistance, severity and mortality. Therefore, numerous efforts have been made nowadays to identify new therapeutic targets. This study aimed to find potential drug targets and vaccine candidates in drug-resistant strains of P. aeruginosa. Extensive antibiotic-resistant and carbapenem-resistant strains of P. aeruginosa with complete genome were selected and ten common hypothetical proteins (HPs) containing more than 200 amino acids were obtained. The structural, functional and immunological predictions of these HPs were performed with the utility of bioinformatics approaches. Two common HPs (Gene ID: 2877781645 and 2877781936) among other investigated proteins were revealed as potential candidates for pharmaceutical and vaccine purposes based on structural and physicochemical properties, functional domains, subcellular localizations, signal peptides, toxicity, virulence factor, antigenicity, allergenicity and immunoinformatic predictions. The consequence of this predictive study will assist in novel drug and vaccine design through experimental investigations.

VicRK and CovR polymorphisms in Streptococcus mutans strains associated with cardiovascular infections

Introduction. Streptococcus mutans, a common species of the oral microbiome, expresses virulence genes promoting cariogenic dental biofilms, persistence in the bloodstream and cardiovascular infections.Gap statement. Virulence gene expression is variable among S. mutans strains and controlled by the transcription regulatory systems VicRK and CovR.Aim. This study investigates polymorphisms in the vicRK and covR loci in S. mutans strains isolated from the oral cavity or from the bloodstream, which were shown to differ in expression of covR, vicRK and downstream genes.Methodology. The transcriptional activities of covR, vicR and vicK were compared by RT-qPCR between blood and oral strains after exposure to human serum. PCR-amplified promoter and/or coding regions of covR and vicRK of 18 strains (11 oral and 7 blood) were sequenced and compared to the reference strain UA159.Results. Serum exposure significantly reduced covR and vicR/K transcript levels in most strains (P<0.05), but reductions were higher in oral than in blood strains. Single-nucleotide polymorphisms (SNPs) were detected in covR regulatory and coding regions, but SNPs affecting the CovR effector domain were only present in two blood strains. Although vicR was highly conserved, vicK showed several SNPs, and SNPs affecting VicK regions important for autokinase activity were found in three blood strains.Conclusions. This study reveals transcriptional and structural diversity in covR and vicR/K, and identifies polymorphisms of functional relevance in blood strains, indicating that covR and vicRK might be important loci for S. mutans adaptation to host selective pressures associated with virulence diversity.

In silico analysis and characterization of medicinal mushroom cystathionine beta-synthase as an angiotensin converting enzyme (ACE) inhibitory protein

Angiotensin-converting enzyme (ACE) regulates blood pressure and has been implicated in several conditions including lung injury, fibrosis and Alzheimer's disease. Medicinal mushroom Ganordema lucidum (Reishi) cystathionine beta-synthase (GlCBS) was previously reported to possess ACE inhibitory activities. However, the inhibitory mechanism of CBS protein remains unreported. Therefore, this study integrates in silico sequencing, structural and functional based-analysis, protein modelling, molecular docking and binding affinity calculation to elucidate the inhibitory mechanism of GlCBS and Lignosus rhinocerus (Tiger milk mushroom) CBS protein (LrCBS) towards ACE. In silico analysis indicates that CBSs from both mushrooms share high similarities in terms of physical properties, structural properties and domain distribution. Protein-protein docking analysis revealed that both GlCBS and LrCBS potentially modulate the C-terminal domain of ACE (C-ACE) activity via regulation of chloride activation and/or prevention of substrate entry. GICBS and LrCBS were also shown to interact with ACE at the same region that presumably inhibits the function of ACE.

The Chinese pine genome and methylome unveil key features of conifer evolution

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

Sequence conservation, domain architectures, and phylogenetic distribution of the HD-GYP type c-di-GMP phosphodiesterases

The HD-GYP domain, named after two of its conserved sequence motifs, was first described in 1999 as a specialized version of the widespread HD phosphohydrolase domain that had additional highly conserved amino acid residues. Domain associations of HD-GYP indicated its involvement in bacterial signal transduction and distribution patterns of this domain suggested that it could serve as a hydrolase of the bacterial second messenger c-di-GMP, in addition to or instead of the EAL domain. Subsequent studies confirmed the ability of various HD-GYP domains to hydrolyze c-di-GMP to linear pGpG and/or GMP. Certain HD-GYP-containing proteins hydrolyze another second messenger, cGAMP, and some HD-GYP domains participate in regulatory protein-protein interactions. The recently solved structures of HD-GYP domains from four distinct organisms clarified the mechanisms of c-di-GMP binding and metal-assisted hydrolysis. However, the HD-GYP domain is poorly represented in public domain databases, which causes certain confusion about its phylogenetic distribution, functions, and domain architectures. Here, we present a refined sequence model for the HD-GYP domain and describe the roles of its most conserved residues in metal and/or substrate binding. We also calculate the numbers of HD-GYPs encoded in various genomes and list the most common domain combinations involving HD-GYP, such as the RpfG (REC-HD-GYP), Bd1817 (DUF3391- HD-GYP), and PmGH (GAF-HD-GYP) protein families. We also provide the descriptions of six HD-GYP-associated domains, including four novel integral membrane sensor domains. This work is expected to stimulate studies of diverse HD-GYP-containing proteins, their N-terminal sensor domains and the signals to which they respond. IMPORTANCE The HD-GYP domain forms class II of c-di-GMP phosphodiesterases that control the cellular levels of the universal bacterial second messenger c-di-GMP and therefore affect flagellar and/or twitching motility, cell development, biofilm formation, and, often, virulence. Despite more than 20 years of research, HD-GYP domains are insufficiently characterized; they are often confused with 'classical' HD domains that are involved in various housekeeping activities and may participate in signaling, hydrolyzing (p)ppGpp and c-di-AMP. This work provides an updated description of the HD-GYP domain, including its sequence conservation, phylogenetic distribution, domain architectures, and the most widespread HD-GYP-containing protein families. This work shows that HD-GYP domains are widespread in many environmental bacteria and are predominant c-di-GMP hydrolases in many lineages, including clostridia and deltaproteobacteria.

Geosmin synthase ges1 knock-down by siRNA in the dikaryotic fungus Tricholoma vaccinum

Genetic manipulation for generating knock-out experiments is essential in deciphering the precise function of a gene. However, dikaryotic fungi pose the inherent challenge of having two allelic versions of each gene, one in each nucleus. In addition, they often are slow-growing and do not withstand protoplasting, which is why Agrobacterium tumefaciens-mediated transformation has been adapted. To obtain knock-out strains, however, is not feasible with a mere deletion construct transformation and screening for deletions in both nuclear copies. Hence, a convenient method using chemically synthesized dicer substrate interfering RNA (DsiRNA) for posttranscriptional interference of targeted mRNA was developed, based on the fungal dicer/argonaute system inherent in fungi for sequence recognition and degradation. A proof-of-principle using this newly established method for knock-down of the volatile geosmin is presented in the dikaryotic fungus Tricholoma vaccinum that is forming ectomycorrhizal symbiosis with spruce trees. The gene ges1, a terpene synthase, was transcribed with a 50-fold reduction in transcript levels in the knockdown strain. The volatile geosmin was slightly reduced, but not absent in the fungus carrying the knockdown construct pointing at low specificity in other terpene synthases known for that class of enzymes.

Genome-Wide Identification and Expression Analysis of the Plant U-Box Protein Gene Family in Phyllostachys edulis

The U-box gene encodes a ubiquitin ligase that contains a U-box domain. The plant U-box (PUB) protein plays an important role in the plant stress response; however, very few studies have investigated the role of these proteins in Moso bamboo (Phyllostachys edulis). Thus, more research on PUB proteins is necessary to understand the mechanisms of stress tolerance in P. edulis. In this study, we identified 121 members of the PUB family in P. edulis (PePUB), using bioinformatics based on the P. edulis V2 genome build. The U-box genes of P. edulis showed an uneven distribution among the chromosomes. Phylogenetic analysis of the U-box genes between P. edulis and Arabidopsis thaliana suggested that these genes can be classified into eight subgroups (Groups I–VIII) based on their structural and phylogenetic features. All U-box genes and the structure of their encoded proteins were identified in P. edulis. We further investigated the expression pattern of PePUB genes in different tissues, including the leaves, panicles, rhizomes, roots, and shoots. The qRT-PCR results showed that expression of three genes, PePUB15, PePUB92, and PePUB120, was upregulated at low temperatures compared to that at 25°C. The expression levels of two PePUBs, PePUB60 and PePUB120, were upregulated under drought stress. These results suggest that the PePUB genes play an important role in resistance to low temperatures and drought in P. edulis. This research provides new insight into the function, diversity, and characterization of PUB genes in P. edulis and provides a basis for understanding their biological roles and molecular mechanisms.

Genome-Wide Analysis of the Soybean Calmodulin-Binding Protein 60 Family and Identification of GmCBP60A-1 Responses to Drought and Salt Stresses

Calmodulin-binding protein 60 (CBP60) members constitute a plant-specific protein family that plays an important role in plant growth and development. In the soybean genome, nineteen CBP60 members were identified and analyzed for their corresponding sequences and structures to explore their functions. Among GmCBP60A-1, which primarily locates in the cytomembrane, was significantly induced by drought and salt stresses. The overexpression of GmCBP60A-1 enhanced drought and salt tolerance in Arabidopsis, which showed better state in the germination of seeds and the root growth of seedlings. In the soybean hairy roots experiment, the overexpression of GmCBP60A-1 increased proline content, lowered water loss rate and malondialdehyde (MDA) content, all of which likely enhanced the drought and salt tolerance of soybean seedlings. Under stress conditions, drought and salt response-related genes showed significant differences in expression in hairy root soybean plants of GmCBP60A-1-overexpressing and hairy root soybean plants of RNAi. The present study identified GmCBP60A-1 as an important gene in response to salt and drought stresses based on the functional analysis of this gene and its potential underlying mechanisms in soybean stress-tolerance.

Characterization of a novel fructosyltransferase InuCA from Lactobacillus crispatus that attaches to the cell surface by electrostatic interaction

Fructosyltransferases (FTases), a group of carbohydrate-active enzymes, synthesize fructooligosaccharides (FOS) and fructans, which are promising prebiotics for human health. Here we originally identified a novel FTase InuCA from L. crispatus, a dominant species in the vaginal microbiotas of human. InuCA was characterized by a shortest C-terminus and the highest isoelectric point among the reported Lactobacillus FTases. InuCA was an inulosucrase and produced a serial of FOS using sucrose as substrate at a moderate temperature. Surprisingly, the C-terminal deletion mutant synthesized oligosaccharides with fructosyl chain longer than that of the wild type, suggesting that the C-terminal part blocked the binding of long-chain receptor. Moreover, InuCA bound to the cell surface by electrostatic interaction, which was dependent on the environmental pH and represented a distinctive binding mode in FTases. The catalytic and structural properties of InuCA will be contributed to the FTases engineering and the knowledge of the adaptation of L. crispatus in the vaginal environment. Importance L. crispatus is one of the most important species in human vaginal microbiotas and its persistence is strongly negatively correlated with the vaginal diseases. Our research reveals that a novel inulosucrase InuCA is present in L. cirspatus. InuCA keeps the ability to synthesize prebiotic fructo-oligosaccharides, although it lacks a large part of the C-terminal region compared to other FTases. Remarkably, the short C-terminus of InuCA blocks the transfructosylation activity for producing oligosaccharides with longer chain, which is meaningful to the directional modification of FTases and the oligosaccharide products. Besides the catalytic activity, InuCA is anchored on the cell surface dependent on the environmental pH and may be also involved in the adhesion of L. crispatus to the vaginal epithelial cells. Since L. crispatus plays an essential role in the normal vaginal micro-ecosystem, the described work will be helpful to elucidate the functional genes and colonization mechanism of the dominant species.

Comparative transcriptome and metabolome analysis reveal glutathione metabolic network and functional genes underlying blue and red-light mediation in maize seedling leaf

BACKGROUND

Light quality severely affects biosynthesis and metabolism-associated process of glutathione. However, the role of specific light is still unclear on the glutathione metabolism. In this article, comparatively transcriptome and metabolome methods are used to fully understand the blue and red-light conditions working on the glutathione metabolism in maize seedling leaf.

RESULTS

There are 20 differently expressed genes and 4 differently expressed metabolites in KEGG pathway of glutathione metabolism. Among them, 12 genes belong to the glutathione S-transferase family, 3 genes belong to the ascorbate peroxidase gene family and 2 genes belong to the ribonucleoside-diphosphate reductase gene family. Three genes, G6PD, SPDS1, and GPX1 belong to the gene family of glucose 6-phosphate dehydrogenase, spermidine synthase, and glutathione peroxidase, respectively. Four differently expressed metabolites are identified. Three of them, Glutathione disulfide, Glutathione, and l-γ-Glutamyl-L-amino acid are decreased while L-Glutamate is increased. In addition, Through PPI analysis, two annotated genes gst16 and DAAT, and 3 unidentified genes 100381533, pco105094 and umc2770, identified as RPP13-like3, BCAT-like1and GMPS, were obtained. By the analysis of protein sequence and PPI network, we predict that pco105094 and umc2770 were involved in the GSSG-GSH and AsA-GSH cycle in the network of glutathione metabolism.

CONCLUSIONS

Compared to red light, blue light remarkably changed the transcription signal transduction and metabolism of glutathione metabolism. Differently expressed genes and metabolic mapped to the glutathione metabolism signaling pathways. In total, we obtained three unidentified genes, and two of them were predicted in current glutathione metabolism network. This result will contribute to the research of glutathione metabolism of maize.

The Role of Outer Membrane Protein(s) Harboring SLH/OprB-Domains in Extracellular Vesicles’ Production in Synechocystis sp. PCC 6803

Cyanobacteria are a group of photosynthetic prokaryotes that contribute to primary production on a global scale. These microorganisms release vesicles to the extracellular environment, spherical nanosized structures, derived essentially from the outer membrane. Even though earlier works in model Gram-negative bacteria have hypothesized that outer membrane stability is crucial in vesicle formation, the mechanisms determining vesicle biogenesis in cyanobacteria remain unknown. Here, we report on the identification of six candidate genes encoding outer membrane proteins harboring SLH/OprB-domains in the genome of the model cyanobacterium Synechocystis sp. PCC 6803. Using a genetics-based approach, one gene was found to encode an essential protein (Slr1841), while the remaining five are not essential for growth under standard conditions. Vesicle production was monitored, and it was found that a mutant in the gene encoding the second most abundant SLH/OprB protein in Synechocystis sp. PCC 6803 outer membrane (Slr1908) produces more vesicles than any of the other tested strains. Moreover, the Slr1908-protein was also found to be important for iron uptake. Altogether, our results suggest that proteins containing the SLH/OprB-domains may have dual biological role, related to micronutrient uptake and to outer membrane stability, which, together or alone, seem to be involved in cyanobacterial vesicle biogenesis.

eggNOG-mapper v2: Functional Annotation, Orthology Assignments, and Domain Prediction at the Metagenomic Scale

Even though automated functional annotation of genes represents a fundamental step in most genomic and metagenomic workflows, it remains challenging at large scales. Here, we describe a major upgrade to eggNOG-mapper, a tool for functional annotation based on precomputed orthology assignments, now optimized for vast (meta)genomic data sets. Improvements in version 2 include a full update of both the genomes and functional databases to those from eggNOG v5, as well as several efficiency enhancements and new features. Most notably, eggNOG-mapper v2 now allows for: 1) de novo gene prediction from raw contigs, 2) built-in pairwise orthology prediction, 3) fast protein domain discovery, and 4) automated GFF decoration. eggNOG-mapper v2 is available as a standalone tool or as an online service at http://eggnog-mapper.embl.de.

eggNOG-mapper v2: Functional Annotation, Orthology Assignments, and Domain Prediction at the Metagenomic Scale

Even though automated functional annotation of genes represents a fundamental step in most genomic and metagenomic workflows, it remains challenging at large scales. Here, we describe a major upgrade to eggNOG-mapper, a tool for functional annotation based on precomputed orthology assignments, now optimized for vast (meta)genomic data sets. Improvements in version 2 include a full update of both the genomes and functional databases to those from eggNOG v5, as well as several efficiency enhancements and new features. Most notably, eggNOG-mapper v2 now allows for: 1) de novo gene prediction from raw contigs, 2) built-in pairwise orthology prediction, 3) fast protein domain discovery, and 4) automated GFF decoration. eggNOG-mapper v2 is available as a standalone tool or as an online service at http://eggnog-mapper.embl.de.

Identification of Genetic Predisposition in Noncirrhotic Portal Hypertension Patients With Multiple Renal Cysts by Integrated Analysis of Whole-Genome and Single-Cell RNA Sequencing

Background and Aims: The multiple renal cysts (MRC) occur in some patients with noncirrhotic portal hypertension (NCPH) could be a subset of ciliopathy. However, the potential genetic influencers and/or determinants in NCPH with MRC are largely unknown. The aim of this study was to explore the potential candidate variants/genes associated with those patients.

Methods: 8,295 cirrhotic patients with portal hypertension were enrolled in cohort 1 and 267 patients affected with NCPH were included in cohort 2. MRC was defined as at least two cysts in both kidneys within a patient detected by ultrasonography or computed tomography. Whole-genome sequencing (WGS) was performed in nine patients (four from cohort 1 and five from cohort 2). Then we integrated WGS and publicly available single-cell RNA sequencing (scRNA-seq) to prioritize potential candidate genes. Genes co-expressed with known pathogenic genes within same cell types were likely associated NCPH with MRC.

Results: The prevalence of MRC in NCPH patients (19.5%, 52/267) was significantly higher than cirrhotic patients (6.2%, 513/8,295). Further, the clinical characteristics of NCPH patients with MRC were distinguishable from cirrhotic patients, including late-onset, more prominent portal hypertension however having preserved liver functions. In the nine whole genome sequenced patients, we identified three patients with early onset harboring compound rare putative pathogenic variants in the known disease gene PKHD1. For the remaining patients, by assessing cilia genes profile in kidney and liver scRNA-seq data, we identified CRB3 was the most co-expressed gene with PKHD1 that highly expressed in ureteric bud cell, kidney stromal cell and hepatoblasts. Moreover, we found a homozygous variant, CRB3 p.P114L, that caused conformational changes in the evolutional conserved domain, which may associate with NCPH with MRC.

Conclusion: ScRNA-seq enables unravelling cell heterogeneity with cell specific gene expression across multiple tissues. With the boosting public accessible scRNA-seq data, we believe our proposed analytical strategy would effectively help disease risk gene identification.

The evolution of N6-methyladenosine regulators in plants

As a reversible modification, N⁶-methyladenosine (m⁶A) plays key roles in series of biological processes. Although it has been found that m⁶A modification is regulated by writers, erasers and readers, their evolutionary processes are still not clearly and systematically described. In the present work, we identified 1592 m⁶A modification regulators from 65 representative plant species and performed the phylogenetic relationships, sequence structure, selection pressure, and codon usage analysis across species. The regulators from different species or subfamilies were distinguishable based on the phylogenetic trees. Although the gene structure was structurally and functionally conserved for each kind of regulators, the unique exon/intron structures and motif organizations were observed among different families. The selection pressure analysis demonstrated that the regulators experienced purifying selection. Interestingly, the selection pressure for the regulators in higher plants was more relaxed, indicating that they might have acquired new functions during evolution. In addition, the different codon usage preferences were observed for the different kinds of m⁶A modification regulators. These results will not only facilitate our understanding of the evolution of m⁶A regulators, but also shed light on how the evolutionary differences affect their functional divergence.

Genome-wide characterization of bZIP gene family identifies potential members involved in flavonoids biosynthesis in Ginkgo biloba L

Ginkgo biloba L. is an ancient relict plant with rich pharmacological activity and nutritional value, and its main physiologically active components are flavonoids and terpene lactones. The bZIP gene family is one of the largest gene families in plants and regulates many processes including pathogen defense, secondary metabolism, stress response, seed maturation, and flower development. In this study, genome-wide distribution of the bZIP transcription factors was screened from G. biloba database in silico analysis. A total of 40 bZIP genes were identified in G. biloba and were divided into 10 subclasses. GbbZIP members in the same group share a similar gene structure, number of introns and exons, and motif distribution. Analysis of tissue expression pattern based on transcriptome indicated that GbbZIP08 and GbbZIP15 were most highly expressed in mature leaf. And the expression level of GbbZIP13 was high in all eight tissues. Correlation analysis and phylogenetic tree analysis suggested that GbbZIP08 and GbbZIP15 might be involved in the flavonoid biosynthesis. The transcriptional levels of 20 GbbZIP genes after SA, MeJA, and low temperature treatment were analyzed by qRT-PCR. The expression level of GbbZIP08 was significantly upregulated under 4°C. Protein–protein interaction network analysis indicated that GbbZIP09 might participate in seed germination by interacting with GbbZIP32. Based on transcriptome and degradome data, we found that 32 out of 117 miRNAs were annotated to 17 miRNA families. The results of this study may provide a theoretical foundation for the functional validation of GbbZIP genes in the future.

Comprehensive Analysis of SRO Gene Family in Sesamum indicum (L.) Reveals Its Association with Abiotic Stress Responses

SIMILAR TO RCD-ONEs (SROs) comprise a small plant-specific gene family which play important roles in regulating numerous growth and developmental processes and responses to environmental stresses. However, knowledge of SROs in sesame (Sesamum indicum L.) is limited. In this study, four SRO genes were identified in the sesame genome. Phylogenetic analysis showed that 64 SROs from 10 plant species were divided into two groups (Group I and II). Transcriptome data revealed different expression patterns of SiSROs over various tissues. Expression analysis showed that Group II SROs, especially SiSRO2b, exhibited a stronger response to various abiotic stresses and phytohormones than those in Group I, implying their crucial roles in response to environmental stimulus and hormone signals. In addition, the co-expression network and protein-protein interaction network indicated that SiSROs are associated with a wide range of stress responses. Moreover, transgenic yeast harboring SiSRO2b showed improved tolerance to salt, osmotic and oxidative stress, indicating SiSRO2b could confer multiple tolerances to transgenic yeast. Taken together, this study not only lays a foundation for further functional dissection of the SiSRO gene family, but also provides valuable gene candidates for genetic improvement of abiotic stress tolerance in sesame.

First Isolation and Molecular Typing of Pathogenic and Intermediate Leptospira Species from Urine of Symptomatic Dogs

Aim of this study was to evaluate, the presence and diversity of Leptospira spp. in blood and urine samples collected from 175 owned-dogs from Sardinia, Italy. After determination of leptospiral infection by microscopic agglutination test (MAT), urine from MAT-positive dogs were examined by real-time polymerase chain reaction (lipL32 rt-PCR) and then isolated by culture. In order to characterize obtained serovars, positive cultures were then subjected to 16S rRNA and secY sequencing, phylogenetic analysis and Multilocus Sequence Typing (MLST). Results showed that seven dogs (4%; 95% CI: 0-55) had Leptospira DNAs in their urine and five strains were isolated from urine cultures. The three different sequence types (ST17, ST198 and ST24) belonging to Leptospira interrogans genomospecies identified by MLST analyses in this study, confirmed that the leptospiral infection was widespread in Sardinian dogs. We also reported the first characterization of a new Leptospira spp. isolated from urine of one dog living in the study area. Whole genome sequencing and phylogenetic analysis, confirmed that this genospecies was closely related to Leptospira hovindhougenii, an intermediate Leptospira spp. with unknown pathogenicity previously isolated from a rat in Denmark. Further studies are required to clarify whether healthy dogs that shed leptospires in their urine could represent a zoonotic risk for humans in this region.

The leptospiral LipL21 and LipL41 proteins exhibit a broad spectrum of interactions with host cell components

Leptospirosis is a globally prevalent zoonotic disease, and is caused by pathogenic spirochetes from the genus Leptospira. LipL21 and LipL41 are lipoproteins expressed strongly on the outer membrane of pathogenic Leptospira spp. Many studies have shown that both proteins are interesting targets for vaccines and diagnosis. However, their role in host-pathogen interactions remains underexplored. Therefore, we evaluated the capacity of LipL21 and LipL41 to bind with glycosaminoglycans (GAGs), the cell receptors and extracellular matrix, and plasma components by ELISA. Both proteins interacted with collagen IV, laminin, E-cadherin, and elastin dose-dependently. A broad-spectrum binding to plasma components was also observed. Only LipL21 interacted with all the GAG components tested, whereas LipL41 presented a concentration-dependent binding only for chondroitin 4 sulfate. Although, both proteins have the ability to interact with fibrinogen, only LipL21 inhibited fibrin clot formation partially. Both proteins exhibited a decrease in plasminogen binding in the presence of amino caproic acid (ACA), a competitive inhibitor of lysine residues, suggesting that their binding occurs via the kringle domains of plasminogen. LipL41, but not LipL21, was able to convert plasminogen to plasmin, and recruit plasminogen from normal human serum, suggesting that the interaction of this protein with plasminogen may occur in physiological conditions. This work provides the first report demonstrating the capacity of LipL21 and LipL41 to interact with a broad range of host components, highlighting their importance in host-Leptospira interactions.

Genome-Wide Identification and Expression Analysis of Potential Antiviral Tripartite Motif Proteins (TRIMs) in Grass Carp (Ctenopharyngodon idella)

Tripartite motif proteins (TRIMs), especially B30.2 domain-containing TRIMs (TRIMs-B30.2), are increasingly well known for their antiviral immune functions in mammals, while antiviral TRIMs are far from being identified in teleosts. In the present study, we identified a total of 42 CiTRIMs from the genome of grass carp, Ctenopharyngodon idella, an important cultured teleost in China, based on hmmsearch and SMART analysis. Among these CiTRIMs, the gene loci of 37 CiTRIMs were located on different chromosomes and shared gene collinearities with homologous counterparts from human and zebrafish genomes. They possessed intact conserved RBCC or RB domain assemblies at their N-termini and eight different domains, including the B30.2 domain, at their C-termini. A total of 19 TRIMs-B30.2 were identified, and most of them were clustered into a large branch of CiTRIMs in the dendrogram. Tissue expression analysis showed that 42 CiTRIMs were universally expressed in various grass carp tissues. A total of 11 significantly differentially expressed CiTRIMs were found in two sets of grass carp transcriptomes during grass carp reovirus (GCRV) infection. Three of them, including Cibtr40, CiTRIM103 and CiTRIM109, which all belonged to TRIMs-B30.2, were associated with the type I interferon response during GCRV infection by weighted network co-expression and gene expression trend analyses, suggesting their involvement in antiviral immunity. These findings may offer useful information for understanding the structure, evolution, and function of TRIMs in teleosts and provide potential antiviral immune molecule markers for grass carp.

Analysis of CcGASA family members in Citrus clementina (Hort. ex Tan.) by a genome-wide approach

The Gibberellic Acid Stimulated Arabidopsis (GASA) proteins were investigated in the study to help understand their possible roles in fruit trees, particularly in Citrus. A total of 18 CcGASA proteins were identified and characterized in Citrus clementina via a genome-wide approach. It was shown that the CcGASA proteins structurally shared a conserved GASA domain but varied considerably in primary sequences and motif compositions. Thus, they could be classified into three major phylogenetic groups, G1~G3, and two groups, G1 and G3 could be further classified into subgroups. The cis- elements on all CcGASA promoters were identified and categorized, and the associated transcription factors were predicted. In addition, the possible interactions between the CcGASA proteins and other proteins were predicted. All the clues suggested that these genes should be involved in defense against biotic and abiotic stresses and in growth and development. The notion was further supported by gene expression analysis that showed these genes were more or less responsive to the treatments of plant hormones (GA₃, SA, ABA and IAA), and infections of citrus canker pathogen Xanthomonas citri. It was noted that both the segmental and the tandem duplications had played a role in the expansion of the CcGASA gene family in Citrus. Our results showed that the members of the CcGASA gene family should have structurally and functionally diverged to different degrees, and hence, the representative group members should be individually investigated to dissect their specific roles.

Structural characterization of a galectin isolated from the marine sponge Chondrilla caribensis with leishmanicidal potential

BACKGROUND

Solving primary structure of lectins leads to an understanding of the physiological roles within an organism and its biotechnological potential. Only eight sponge lectins have had their primary structure fully determined.

METHODS

The primary structure of CCL, Chondrilla caribensis lectin, was determined by tandem mass spectrometry. The three-dimensional structure was predicted and the protein-carbohydrate interaction analysed by molecular docking. Furthermore, the anti-leishmanial activity was observed by assays with Leishmania infantum.

RESULTS

The amino acid sequence consists of 142 amino acids with a calculated molecular mass of 15,443 Da. The lectin has a galectin-like domain architecture. As observed in other sponge galectins, the signature sequence of a highly conserved domain was also identified in CCL with some modifications. CCL exhibits a typical galectin structure consisting of a β-sandwich. Molecular docking showed that the amino acids interacting with CCL ligands at the monosaccharide binding site are mostly the same as those conserved in this family of lectins. Through its interaction with L. infantum glycans, CCL was able to inhibit the development of this parasite. CCL also induced apoptosis after eliciting ROS production and altering the membrane integrity of Leishmania infantum promastigote.

CONCLUSIONS

CCL joins the restricted group of sponge lectins with determined primary structure and very high biotechnological potential owing to its promising results against pathogens that cause Leishmaniasis.

GENERAL SIGNIFICANCE

As the determination of primary structure is important for biological studies, now CCL can become a sponge galectin with an exciting future in the field of human health.

ATG8-Interacting Motif: Evolution and Function in Selective Autophagy of Targeting Biological Processes

Autophagy is an evolutionarily conserved vacuolar process functioning in the degradation of cellular components for reuse. In plants, autophagy is generally activated upon stress and its regulation is executed by numbers of AuTophaGy-related genes (ATGs), of which the ATG8 plays a dual role in both biogenesis of autophagosomes and recruitment of ATG8-interacting motif (AIM) anchored selective autophagy receptors (SARs). Such motif is either termed as AIM or ubiquitin-interacting motif (UIM), corresponding to the LC3-interacting region (LIR)/AIM docking site (LDS) or the UIM docking site (UDS) of ATG8, respectively. To date, dozens of AIM or UIM containing SARs have been characterized. However, the knowledge of these motifs is still obscured. In this review, we intend to summarize the current understanding of SAR proteins and discuss the conservation and diversification of the AIMs/UIMs, expectantly providing new insights into the evolution of them in various biological processes in plants.

PRGdb 4.0: an updated database dedicated to genes involved in plant disease resistance process

The Plant Resistance Genes database (PRGdb; http://prgdb.org/prgdb4/) has been greatly expanded, keeping pace with the increasing amount of available knowledge and data (sequenced proteomes, cloned genes, public analysis data, etc.). The easy-to-use style of the database website has been maintained, while an updated prediction tool, more data and a new section have been added. This new section will contain plant resistance transcriptomic experiments, providing additional easy-to-access experimental information. DRAGO3, the tool for automatic annotation and prediction of plant resistance genes behind PRGdb, has been improved in both accuracy and sensitivity, leading to more reliable predictions. PRGdb offers 199 reference resistance genes and 586.652 putative resistance genes from 182 sequenced proteomes. Compared to the previous release, PRGdb 4.0 has increased the number of reference resistance genes from 153 to 199, the number of putative resistance genes from 177K from 76 proteomes to 586K from 182 sequenced proteomes. A new section has been created that collects plant-pathogen transcriptomic data for five species of agricultural interest. Thereby, with these improvements and data expansions, PRGdb 4.0 aims to serve as a reference to the plant scientific community and breeders worldwide, helping to further study plant resistance mechanisms that contribute to fighting pathogens.

Comprehensive Analysis of the SUV Gene Family in Allopolyploid Brassica napus and Its Diploid Ancestors

SUV (the Suppressor of variegation [Su(var)] homologs and related) gene family is a subgroup of the SET gene family. According to the SRA domain and WIYLD domain distributions, it can be divided into two categories, namely SUVH (the Suppressor of variegation [Su(var)] homologs) and SUVR (the Suppressor of variegation [Su(var)] related). In this study, 139 SUV genes were identified in allopolyploid Brassica napus and its diploid ancestors, and their evolutionary relationships, protein properties, gene structures, motif distributions, transposable elements, cis-acting elements and gene expression patterns were analyzed. Our results showed that the SUV gene family of B. napus was amplified during allopolyploidization, in which the segmental duplication and TRD played critical roles. After the separation of Brassica and Arabidopsis lineages, orthologous gene analysis showed that many SUV genes were lost during the evolutionary process in B. rapa, B. oleracea and B. napus. The analysis of the gene and protein structures and expression patterns of 30 orthologous gene pairs which may have evolutionary relationships showed that most of them were conserved in gene structures and protein motifs, but only four gene pairs had the same expression patterns.

Centromere-Specific Retrotransposons and Very-Long-Chain Fatty Acid Biosynthesis in the Genome of Yellowhorn (Xanthoceras sorbifolium, Sapindaceae), an Oil-Producing Tree With Significant Drought Resistance

In-depth genome characterization is still lacking for most of biofuel crops, especially for centromeres, which play a fundamental role during nuclear division and in the maintenance of genome stability. This study applied long-read sequencing technologies to assemble a highly contiguous genome for yellowhorn (Xanthoceras sorbifolium), an oil-producing tree, and conducted extensive comparative analyses to understand centromere structure and evolution, and fatty acid biosynthesis. We produced a reference-level genome of yellowhorn, ∼470 Mb in length with ∼95% of contigs anchored onto 15 chromosomes. Genome annotation identified 22,049 protein-coding genes and 65.7% of the genome sequence as repetitive elements. Long terminal repeat retrotransposons (LTR-RTs) account for ∼30% of the yellowhorn genome, which is maintained by a moderate birth rate and a low removal rate. We identified the centromeric regions on each chromosome and found enrichment of centromere-specific retrotransposons of LINE1 and Gypsy in these regions, which have evolved recently (∼0.7 MYA). We compared the genomes of three cultivars and found frequent inversions. We analyzed the transcriptomes from different tissues and identified the candidate genes involved in very-long-chain fatty acid biosynthesis and their expression profiles. Collinear block analysis showed that yellowhorn shared the gamma (γ) hexaploidy event with Vitis vinifera but did not undergo any further whole-genome duplication. This study provides excellent genomic resources for understanding centromere structure and evolution and for functional studies in this important oil-producing plant.

Synthaser: a CD-Search enabled Python toolkit for analysing domain architecture of fungal secondary metabolite megasynth(et)ases

BACKGROUND

Fungi are prolific producers of secondary metabolites (SMs), which are bioactive small molecules with important applications in medicine, agriculture and other industries. The backbones of a large proportion of fungal SMs are generated through the action of large, multi-domain megasynth(et)ases such as polyketide synthases (PKSs) and nonribosomal peptide synthetases (NRPSs). The structure of these backbones is determined by the domain architecture of the corresponding megasynth(et)ase, and thus accurate annotation and classification of these architectures is an important step in linking SMs to their biosynthetic origins in the genome.

RESULTS

Here we report synthaser, a Python package leveraging the NCBI's conserved domain search tool for remote prediction and classification of fungal megasynth(et)ase domain architectures. Synthaser is capable of batch sequence analysis, and produces rich textual output and interactive visualisations which allow for quick assessment of the megasynth(et)ase diversity of a fungal genome. Synthaser uses a hierarchical rule-based classification system, which can be extensively customised by the user through a web application ( http://gamcil.github.io/synthaser ). We show that synthaser provides more accurate domain architecture predictions than comparable tools which rely on curated profile hidden Markov model (pHMM)-based approaches; the utilisation of the NCBI conserved domain database also allows for significantly greater flexibility compared to pHMM approaches. In addition, we demonstrate how synthaser can be applied to large scale genome mining pipelines through the construction of an Aspergillus PKS similarity network.

CONCLUSIONS

Synthaser is an easy to use tool that represents a significant upgrade to previous domain architecture analysis tools. It is freely available under a MIT license from PyPI ( https://pypi.org/project/synthaser ) and GitHub ( https://github.com/gamcil/synthaser ).

Hymenoptera Genome Database: new genomes and annotation datasets for improved go enrichment and orthologue analyses

We report an update of the Hymenoptera Genome Database (HGD; http://HymenopteraGenome.org), a genomic database of hymenopteran insect species. The number of species represented in HGD has nearly tripled, with fifty-eight hymenopteran species, including twenty bees, twenty-three ants, eleven wasps and four sawflies. With a reorganized website, HGD continues to provide the HymenopteraMine genomic data mining warehouse and JBrowse/Apollo genome browsers integrated with BLAST. We have computed Gene Ontology (GO) annotations for all species, greatly enhancing the GO annotation data gathered from UniProt with more than a ten-fold increase in the number of GO-annotated genes. We have also generated orthology datasets that encompass all HGD species and provide orthologue clusters for fourteen taxonomic groups. The new GO annotation and orthology data are available for searching in HymenopteraMine, and as bulk file downloads.

Analysis on Drug-Resistance-Associated Mutations among Multidrug-Resistant Mycobacterium tuberculosis Isolates in China

As the causative bacteria of tuberculosis, Mycobacteriumtuberculosis (M. tb) is aggravated by the emergence of its multidrug-resistant isolates in China. Mutations of six of the most frequently reported resistant genes (rpoB, katG, inhA, embB, gyrA, and rpsL) were detected for rifampicin (RIF), isoniazid (INH), ethambutol (EMB), ofloxacin (OFX), and streptomycin (STR) in this study. The amino acid missense mutations (MMs) and their corresponding single nucleotide polymorphism mutations for all drug-resistant (DR) isolates are described in detail. All isolates were divided into non-extensively drug-resistant (Non-XDR) and preXDR/XDR groups. No statistical differences were detected among MMs and linked MMs (LMs) between the two groups, except for rpsL 88 (p = 0.037). In the preXDR/XDR group, the occurrence of MMs in rpoB, katG, and inhA developed phenotypic resistance and MMs of rpoB 531, katG 315, rpsL 43, and rpsL 88 could develop high levels of DR. It is necessary to carry out epidemiological investigations of DR gene mutations in the local region, and thus provide necessary data to support the design of new technologies for rapid detection of resistant M. tb and the optimization of detection targets.

The Methods of Digging for "Gold" within the Salt: Characterization of Halophilic Prokaryotes and Identification of Their Valuable Biological Products Using Sequencing and Genome Mining Tools

Halophiles, the salt-loving organisms, have been investigated for at least a hundred years. They are found in all three domains of life, namely Archaea, Bacteria, and Eukarya, and occur in saline and hypersaline environments worldwide. They are already a valuable source of various biomolecules for biotechnological, pharmaceutical, cosmetological and industrial applications. In the present era of multidrug-resistant bacteria, cancer expansion, and extreme environmental pollution, the demand for new, effective compounds is higher and more urgent than ever before. Thus, the unique metabolism of halophilic microorganisms, their low nutritional requirements and their ability to adapt to harsh conditions (high salinity, high pressure and UV radiation, low oxygen concentration, hydrophobic conditions, extreme temperatures and pH, toxic compounds and heavy metals) make them promising candidates as a fruitful source of bioactive compounds. The main aim of this review is to highlight the nucleic acid sequencing experimental strategies used in halophile studies in concert with the presentation of recent examples of bioproducts and functions discovered in silico in the halophile's genomes. We point out methodological gaps and solutions based on in silico methods that are helpful in the identification of valuable bioproducts synthesized by halophiles. We also show the potential of an increasing number of publicly available genomic and metagenomic data for halophilic organisms that can be analysed to identify such new bioproducts and their producers.

PQN-59 antagonizes microRNA-mediated repression during post-embryonic temporal patterning and modulates translation and stress granule formation in C. elegans

microRNAs (miRNAs) are potent regulators of gene expression that function in a variety of developmental and physiological processes by dampening the expression of their target genes at a post-transcriptional level. In many gene regulatory networks (GRNs), miRNAs function in a switch-like manner whereby their expression and activity elicit a transition from one stable pattern of gene expression to a distinct, equally stable pattern required to define a nascent cell fate. While the importance of miRNAs that function in this capacity are clear, we have less of an understanding of the cellular factors and mechanisms that ensure the robustness of this form of regulatory bistability. In a screen to identify suppressors of temporal patterning phenotypes that result from ineffective miRNA-mediated target repression, we identified pqn-59, an ortholog of human UBAP2L, as a novel factor that antagonizes the activities of multiple heterochronic miRNAs. Specifically, we find that depletion of pqn-59 can restore normal development in animals with reduced lin-4 and let-7-family miRNA activity. Importantly, inactivation of pqn-59 is not sufficient to bypass the requirement of these regulatory RNAs within the heterochronic GRN. The pqn-59 gene encodes an abundant, cytoplasmically-localized, unstructured protein that harbors three essential "prion-like" domains. These domains exhibit LLPS properties in vitro and normally function to limit PQN-59 diffusion in the cytoplasm in vivo. Like human UBAP2L, PQN-59's localization becomes highly dynamic during stress conditions where it re-distributes to cytoplasmic stress granules and is important for their formation. Proteomic analysis of PQN-59 complexes from embryonic extracts indicates that PQN-59 and human UBAP2L interact with orthologous cellular components involved in RNA metabolism and promoting protein translation and that PQN-59 additionally interacts with proteins involved in transcription and intracellular transport. Finally, we demonstrate that pqn-59 depletion reduces protein translation and also results in the stabilization of several mature miRNAs (including those involved in temporal patterning). These data suggest that PQN-59 may ensure the bistability of some GRNs that require miRNA functions by promoting miRNA turnover and, like UBAP2L, enhancing protein translation.

Characterization of large yellow croaker (Larimichthys crocea) osteoprotegerin and its role in the innate immune response against to Vibrio alginolyticus

Osteoprotegerin (OPG) is a member of the tumor necrosis factor receptor superfamily, contributing to inflammation, apoptosis, and differentiation. However, the function of OPG in the host immune system of teleosts remains unclear. Here, we cloned the cDNA of the LcOPG gene from large yellow croaker. LcOPG mRNA was expressed in all analyzed tissues and was upregulated by Vibrio alginolyticus infection in immune tissues and monocytes/macrophages (MO/MФ). Subsequently, the LcOPG protein was expressed and purified using a prokaryotic expression system. Recombinant LcOPG protein (rLcOPG) treatment suppressed V. alginolyticus-induced pro-inflammatory cytokine and enhanced V. alginolyticus-induced anti-inflammatory cytokine mRNA expression. Furthermore, rLcOPG decreased V. alginolyticus-induced MO/MФ apoptosis. Therefore, the results indicate that LcOPG might play a role in the immune response of V. alginolyticus-infected large yellow croaker.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Sensory perception in bacterial cyclic diguanylate signal transduction

Cyclic diguanylate (c-di-GMP) signal transduction systems provide bacteria with the ability to sense changing cell status or environmental conditions and then execute suitable physiological and social behaviors in response. In this review, we provide a comprehensive census of the stimuli and receptors that are linked to the modulation of intracellular c-di-GMP. Emerging evidence indicates that c-di-GMP networks sense light, surfaces, energy, redox potential, respiratory electron acceptors, temperature, and structurally diverse biotic and abiotic chemicals. Bioinformatic analysis of sensory domains in diguanylate cyclases and c-di-GMP-specific phosphodiesterases as well as the receptor complexes associated with them reveals that these functions are linked to a diverse repertoire of protein domain families. We describe the principles of stimulus perception learned from studying these modular sensory devices, illustrate how they are assembled in varied combinations with output domains, and summarize a system for classifying these sensor proteins based on their complexity. Biological information processing via c-di-GMP signal transduction not only is fundamental to bacterial survival in dynamic environments but also is being used to engineer gene expression circuitry and synthetic proteins with à la carte biochemical functionalities.

RIP3 Associates with RIP1, TRIF, MAVS, and Also IRF3/7 in Host Innate Immune Signaling in Large Yellow Croaker Larimichthys crocea

Receptor-interacting protein 3 (RIP3) has been demonstrated to be a key regulator not only in cell death pathways including apoptosis and necroptosis but also in inflammation and host immune responses. In this study, a RIP3 ortholog named Lc-RIP3 is identified in large yellow croaker (Larimichthys crocea). The open reading frame (ORF) of Lc-RIP3 is 1524 bp long and encodes a protein of 507 amino acids (aa). The deduced Lc-RIP3 protein has an N-terminal kinase domain and a C-terminal RHIM domain, and the genome organization of Lc-RIP3 is conserved in teleosts with 12 exons and 11 introns but is different from that in mammals, which comprises 10 exons and 9 introns. Confocal microscopy revealed that Lc-RIP3 is a cytosolic protein. The expression analysis at the mRNA level indicated that Lc-RIP3 is ubiquitously distributed in various tissues/organs, and could be up-regulated under poly I:C, LPS, PGN, and Pseudomonas plecoglossicida stimulation in vivo. Notably, Lc-RIP3 could induce NF-κB but not IRF3 activation. In addition, Lc-RIP3 co-expression with Lc-TRIF, Lc-MAVS, or Lc-IRF3 significantly abolishes the activation of NF-κB but enhances the induction of IRF3 activity. Moreover, NF-κB activity could be up-regulated when Lc-RIP3 is co-expressed with Lc-RIP1 or Lc-IRF7. These results collectively indicate that Lc-RIP3 acts as an important regulator in host innate immune signaling in teleosts.

Extending the Horizon of Homology Detection with Coevolution-based Structure Prediction

Traditional sequence analysis algorithms fail to identify distant homologies when they lie beyond a detection horizon. In this review, we discuss how co-evolution-based contact and distance prediction methods are pushing back this homology detection horizon, thereby yielding new functional insights and experimentally testable hypotheses. Based on correlated substitutions, these methods divine three-dimensional constraints among amino acids in protein sequences that were previously devoid of all annotated domains and repeats. The new algorithms discern hidden structure in an otherwise featureless sequence landscape. Their revelatory impact promises to be as profound as the use, by archaeologists, of ground-penetrating radar to discern long-hidden, subterranean structures. As examples of this, we describe how triplicated structures reflecting longin domains in MON1A-like proteins, or UVR-like repeats in DISC1, emerge from their predicted contact and distance maps. These methods also help to resolve structures that do not conform to a "beads-on-a-string" model of protein domains. In one such example, we describe CFAP298 whose ubiquitin-like domain was previously challenging to perceive owing to a large sequence insertion within it. More generally, the new algorithms permit an easier appreciation of domain families and folds whose evolution involved structural insertion or rearrangement. As we exemplify with α1-antitrypsin, coevolution-based predicted contacts may also yield insights into protein dynamics and conformational change. This new combination of structure prediction (using innovative co-evolution based methods) and homology inference (using more traditional sequence analysis approaches) shows great promise for bringing into view a sea of evolutionary relationships that had hitherto lain far beyond the horizon of homology detection.

An ancient RAB5 governs the formation of additional vacuoles and cell shape in petunia petals

Homologous ("canonical") RAB5 proteins regulate endosomal trafficking to lysosomes in animals and to the central vacuole in plants. Epidermal petal cells contain small vacuoles (vacuolinos) that serve as intermediate stations for proteins on their way to the central vacuole. Here, we show that transcription factors required for vacuolino formation in petunia induce expression of RAB5a. RAB5a defines a previously unrecognized clade of canonical RAB5s that is evolutionarily and functionally distinct from ARA7-type RAB5s, which act in trafficking to the vacuole. Loss of RAB5a reduces cell height and abolishes vacuolino formation, which cannot be rescued by the ARA7 homologs, whereas constitutive RAB5a (over)expression alters the conical cell shape and promotes homotypic vacuolino fusion, resulting in oversized vacuolinos. These findings provide a rare example of how gene duplication and neofunctionalization increased the complexity of membrane trafficking during evolution and suggest a mechanism by which cells may form multiple vacuoles with distinct content and function.

A novel decrystallizing protein CxEXL22 from Arthrobotrys sp. CX1 capable of synergistically hydrolyzing cellulose with cellulases

A novel expansin-like protein (CxEXL22) has been identified and characterized from newly isolated Arthrobotrys sp. CX1 that can cause cellulose decrystallization. Unlike previously reported expansin-like proteins from microbes, CxEXL22 has a parallel β-sheet domain at the N terminal, containing many hydrophobic residues to form the hydrophobic surface as part of the groove. The direct phylogenetic relationship implied the genetic transfers occurred from nematode to nematicidal fungal Arthrobotrys sp. CX1. CxEXL22 showed strong activity for the hydrolysis of hydrogen bonds between cellulose molecules, especially when highly crystalline cellulose was used as substrate. The hydrolysis efficiency of Avicel was increased 7.9-fold after pretreating with CxEXL22. The rupture characterization of crystalline region indicated that CxEXL22 strongly binds cellulose and breaks up hydrogen bonds in the crystalline regions of cellulose to split cellulose chains, causing significant depolymerization to expose much more microfibrils and enhances cellulose accessibility.

The first chromosome-level genome assembly of a green lacewing Chrysopa pallens and its implication for biological control

Many lacewing species (Insecta: Neuroptera) are important predators of pests with great potential in biological control. So far, there is no chromosome‐level published genome available for Neuroptera. Here we report a high‐quality chromosome‐level reference genome for a green lacewing species Chrysopa pallens (Neuroptera: Chrysopidae), which is one of the most important insect natural enemies used in pest biocontrol. The genome was sequenced using a combination of PacBio and Hi‐C technologies and assembled into seven chromosomes with a total size of 517.21 Mb, occupying 96.07% of the genome sequence. A total of 12,840 protein‐coding genes were identified and approximately 206.21 Mb of repeated sequences were annotated. Phylogenetic analyses indicated that C. pallens diverged from its common ancestor with Tribolium castaneum (Coleoptera) approximately 300 million years ago. The gene families involved in digestion, detoxification, chemoreception, carbohydrate metabolism, immunity, nerves and development were significantly expanded, revealing the potential genomic basis for the polyphagia of C. pallens and its role as an excellent biocontrol agent. This high‐quality genome of C. pallens will provide an important genomic resource for future population genetics, evolutionary and phylogenetic investigations of Chrysopidae as well as comparative genomic studies of Neuropterida and other insects.

The non-canonical functions of telomerase reverse transcriptase gene GlTert on regulating fungal growth, oxidative stress, and ganoderic acid biosynthesis in Ganoderma lucidum

The telomerase reverse transcriptase (TERT) is the core catalytic subunit of telomerase. Its canonical function is synthesizing telomeric repeats to maintain telomere length and chromosomal stability. Accumulating evidence suggests that TERT has other important fundamental functions in addition to its catalytic telomere repeat synthesis activity. However, the non-canonical roles of TERT independent of its enzymatic activity are not clear in filamentous fungi. In the present study, we characterized the GlTert gene in Ganoderma lucidum. The non-canonical roles of GlTert were explored using GlTert-silenced strains (Terti8 and Terti25) obtained by RNA interference. Silencing GlTert delayed the fungal growth, decreased the length between hyphal branches, and induced fungal resistance to oxidative stress in G. ludicum. Further examination revealed that the intracellular ROS (reactive oxygen species) levels were increased while the enzyme activities of the antioxidant systems (superoxide dismutase, catalase, glutathione peroxidase, and ascorbate peroxidase) were decreased in GlTert-silenced strains. In addition, silencing GlTert decreased the ganoderic acid (GA) biosynthesis of G. lucidum. Taken together, our results indicate that GlTert plays a fundamental function on fungal growth, oxidative stress, and GA biosynthesis in G. lucidum, providing new insights for the canonical functions of TERT in filamentous fungi. KEY POINTS: • GlTert affected fungal growth and hyphal branching of G. lucidum. • Silencing GlTert increased the intracellular ROS levels of G. lucidum. • GlTert regulated GA biosynthesis of G. lucidum.

Unraveling Toxoplasma gondii GT1 Strain Virulence and New Protein-Coding Genes with Proteogenomic Analyses

Toxoplasma gondii is one of the most widespread parasites of great relevance to planetary health. It infects approximately one-third of the world population. T. gondii establishes itself in warm-blooded animals and causes adverse health outcomes, particularly in immunocompromised patients. T. gondii is also widely used as a model organism to study other related apicomplexan parasites, which requires a deeper understanding of its molecular biology. Type I strains (GT1 and RH) of T. gondii are considered the most virulent forms. The whole-genome sequencing of T. gondii annotated 8460 predicted gene models in the parasite. To this end, the proteogenomics technology allows harnessing of mass spectrometry (MS)-derived proteomic data to unravel new protein-coding genes, not to mention validation and correction of the existing gene models. In this study using the proteogenomic approach, we report the identification of 31 novel protein-coding genes while reannotating 88 existing gene models. Notably, the genome annotations were corrected for genes, such as SAG5C, GRA6, ROP4, ROP5, and ROP26. The associated proteins are known to play important roles in host-parasite interactions, particularly in relation to parasite virulence, suppression of host immune response, and distinctively pertinent for the survival of the parasite inside the host system. These new findings offer new insights, informing planetary health broadly and the knowledge base on T. gondii virulence specifically. The proteogenomics approach also provides a concrete example to study related apicomplexan organisms of relevance to planetary health, and so as to develop new diagnostics and therapeutics against toxoplasmosis and related diseases.

Genome-Wide Identification and Analysis of the MADS-Box Gene Family in American Beautyberry (Callicarpa americana)

The MADS-box gene family encodes a number of transcription factors that play key roles in various plant growth and development processes from response to environmental cues to cell differentiation and organ identity, especially the floral organogenesis, as in the prominent ABCDE model of flower development. Recently, the genome of American beautyberry (Callicarpa americana) has been sequenced. It is a shrub native to the southern region of United States with edible purple-colored berries; it is a member of the Lamiaceae family, a family of medical and agricultural importance. Seventy-eight MADS-box genes were identified from 17 chromosomes of the C. americana assembled genome. Peptide sequences blast and analysis of phylogenetic relationships with MADS-box genes of Sesame indicum, Solanum lycopersicum, Arabidopsis thaliana, and Amborella trichopoda were performed. Genes were separated into 32 type I and 46 type II MADS-box genes. C. americana MADS-box genes were clustered into four groups: MIKC^C, MIKC*, Mα-type, and Mγ-type, while the Mβ-type group was absent. Analysis of the gene structure revealed that from 1 to 15 exons exist in C. americana MADS-box genes. The number of exons in type II MADS-box genes (5–15) greatly exceeded the number in type I genes (1–9). The motif distribution analysis of the two types of MADS-box genes showed that type II MADS-box genes contained more motifs than type I genes. These results suggested that C. americana MADS-box genes type II had more complex structures and might have more diverse functions. The role of MIKC-type MADS-box genes in flower and fruit development was highlighted when the expression profile was analyzed in different organs transcriptomes. This study is the first genome-wide analysis of the C. americana MADS-box gene family, and the results will further support any functional and evolutionary studies of C. americana MADS-box genes and serve as a reference for related studies of other plants in the medically important Lamiaceae family.

Expression of immunoglobulin constant domain genes in neurons of the mouse central nervous system

General consensus states that immunoglobulins are exclusively expressed by B lymphocytes to form the first line of defense against common pathogens. Here, we provide compelling evidence for the expression of two heavy chain immunoglobulin genes in subpopulations of neurons in the mouse brain and spinal cord. RNA isolated from excitatory and inhibitory neurons through ribosome affinity purification revealed Ighg3 and Ighm transcripts encoding for the constant (Fc), but not the variable regions of IgG3 and IgM. Because, in the absence of the variable immunoglobulin regions, these transcripts lack the canonical transcription initiation site used in lymphocytes, we screened for alternative 5' transcription start sites and identified a novel 5' exon adjacent to a proposed promoter element. Immunohistochemical, Western blot, and in silico analyses strongly support that these neuronal transcripts are translated into proteins containing four Immunoglobulin domains. Our data thus demonstrate the expression of two Fc-encoding genes Ighg3 and Ighm in spinal and supraspinal neurons of the murine CNS and suggest a hitherto unknown function of the encoded proteins.

Two Variants of KIT Causing White Patterning in Stock-Type Horses

Over 30 polymorphisms in the KIT Proto-Oncogene Receptor Tyrosine Kinase (KIT) gene have been implicated in white spotting patterns ranging from small areas to full dermal depigmentation in the horse. We performed a candidate-gene exon sequencing approach on KIT and MITF, 2 known causatives of white spotting patterns, within 2 families of horses of unknown white spotting. Family 1 (Fam1, N = 5) consisted of a Quarter Horse stallion and 4 offspring with white spotting pattern ranging from legs, lower ventral, and head regions with jagged borders, to almost complete white. The second family (Fam2, N = 7) consisted of 6 half-sibling American Paint Horse/Quarter Horse and their dam, demonstrating unpigmented limbs with belly spots and an extensive white patterning on the face. This approach resulted in 2 variants significantly associated with familial phenotypes, where Fam1 variant is an indel leading to a frameshift mutation, and Fam2 a non-synonymous SNP. We validated the variants within an unrelated population of horses (Fam2 variant, P = 0.00271944) as well as for protein functional impact with ExPASy, Protter, Phyre2, SMART, PROVEAN, SIFT, and I-TASSER, confirming the reported associations. Fam1 associated variant, deemed W31, alters the protein sequence, leading to an early stop codon truncating the normal amino acid sequence from 972 to just 115 amino acids. Fam2 associated variant, deemed W32, may have a subtle impact on receptor function or could be in linkage with a non-coding or regulatory change creating the mild spotting pattern observed in this family.

PI3K Plays an Essential Role in Planarian Regeneration and Tissue Maintenance

Phosphatidylinositol 3-kinase (PI3K) signaling plays a central role in various biological processes, and its abnormality leads to a broad spectrum of human diseases, such as cancer, fibrosis, and immunological disorders. However, the mechanisms by which PI3K signaling regulates the behavior of stem cells during regeneration are poorly understood. Planarian flatworms possess abundant adult stem cells (called neoblasts) allowing them to develop remarkable regenerative capabilities, thus the animals represent an ideal model for studying stem cells and regenerative medicine in vivo. In this study, the spatiotemporal expression pattern of Djpi3k, a PI3K ortholog in the planarian Dugesia japonica, was investigated and suggests its potential role in wound response and tissue regeneration. A loss-of-function study was conducted using small molecules and RNA interference technique, providing evidence that PI3K signaling is required for blastema regrowth and cilia maintenance during planarian regeneration and homeostasis. Interestingly, the mitotic and apoptotic responses to amputation are substantially abated in PI3K inhibitor-treated regenerating animals, while knockdown of Djpi3k alleviates the mitotic response and postpones the peak of apoptotic cell death, which may contribute to the varying degrees of regenerative defects induced by the pharmacological and genetic approaches. These observations reveal novel roles for PI3K signaling in the regulation of the cellular responses to amputation during planarian regeneration and provide insights for investigating the disease-related genes in the regeneration-competent organism in vivo.

Whole-Genome Duplication and Purifying Selection Contributes to the Functional Redundancy of Auxin Response Factor (ARF) Genes in Foxtail Millet (Setaria italica L.)

Auxin response factors (ARFs) play crucial roles in auxin-mediated response, whereas molecular genetics of ARF genes was seldom investigated in Setaria italica, an important crop and C₄ model plant. In the present study, genome-wide evolutionary analysis of ARFs was performed in S. italica. Twenty-four SiARF genes were identified and unevenly distributed on eight of the nine chromosomes in S. italica. Duplication mode exploration implied that 13 SiARF proteins were originated from whole-genome duplication and suffered purifying selection. Phylogeny reconstruction of SiARFs by maximum likelihood and neighbor-joining trees revealed SiARFs could be divided into four clades. SiARFs clustered within the same clade shared similar gene structure and protein domain composition, implying functional redundancy. Moreover, amino acid composition of the middle regions was conserved in SiARFs belonged to the same clade. SiARFs were categorized into either activators or repressors according to the enrichment of specific amino acids. Intrinsic disorder was featured in the middle regions of ARF activators. Finally, expression profiles of SiARFs under hormone and abiotic stress treatment not only revealed their potential function in stress response but also indicate their functional redundancy. Overall, our results provide insights into evolutionary aspects of SiARFs and benefit for further functional characterization.