Coding Region Research Articles

Highly accurate classification and discovery of microbial protein-coding gene functions using FunGeneTyper: an extensible deep learning framework.

High-throughput DNA sequencing technologies decode tremendous amounts of microbial protein-coding gene sequences. However, accurately assigning protein functions to novel gene sequences remain a challenge. To this end, we developed FunGeneTyper, an extensible framework with two new deep learning models (i.e., FunTrans and FunRep), structured databases, and supporting resources for achieving highly accurate (Accuracy > 0.99, F1-score > 0.97) and fine-grained classification of antibiotic resistance genes (ARGs) and virulence factor genes. Using an experimentally confirmed dataset of ARGs comprising remote homologous sequences as the test set, our framework achieves by-far-the-best performance in the discovery of new ARGs from human gut (F1-score: 0.6948), wastewater (0.6072), and soil (0.5445) microbiomes, beating the state-of-the-art bioinformatics tools and sequence alignment-based (F1-score: 0.0556-0.5065) and domain-based (F1-score: 0.2630-0.5224) annotation approaches. Furthermore, our framework is implemented as a lightweight, privacy-preserving, and plug-and-play neural network module, facilitating its versatility and accessibility to developers and users worldwide. We anticipate widespread utilization of FunGeneTyper (https://github.com/emblab-westlake/FunGeneTyper) for precise classification of protein-coding gene functions and the discovery of numerous valuable enzymes. This advancement will have a significant impact on various fields, including microbiome research, biotechnology, metagenomics, and bioinformatics.

Knockdown of NFIC Promotes Bovine Myoblast Proliferation through the CENPF/CDK1 Axis.

As cellular transcription factors and DNA replicators, nuclear factor I (NFI) family members play an important role in mammalian development. However, there is still a lack of research on the muscle regeneration of NFI family members in cattle. In this study, the analysis of NFI family factors was conducted on their characterization, phylogenetics, and functional domains. We found that NFI family members were relatively conserved among different species, but there was heterogeneity in amino acid sequences, DNA coding sequences, and functional domain among members. Furthermore, among NFI family factors, we observed that NFIC exhibited highly expression in bovine muscle tissues, particularly influencing the expression of proliferation marker genes in myoblasts. To investigate the influence of NFIC on myoblast proliferation, we knocked down NFIC (si-NFIC) and found that the proliferation of myoblasts was significantly promoted. In terms of regulation mechanism, we identified that si-NFIC could counteract the inhibitory effect of the cell cycle inhibitor RO-3306. Interestingly, CENPF, as the downstream target gene of NFIC, could affect the expression of CDK1, CCNB1, and actively regulate the cell cycle pathway and cell proliferation. In addition, when CENPF was knocked down, the phosphorylation of p53 and the expression of Bax were increased, but the expression of Bcl2 was inhibited. Our findings mainly highlight the mechanism by which NFIC acts on the CENPF/CDK1 axis to regulate the proliferation of bovine myoblasts.

Mapping medically relevant RNA isoform diversity in the aged human frontal cortex with deep long-read RNA-seq.

Determining whether the RNA isoforms from medically relevant genes have distinct functions could facilitate direct targeting of RNA isoforms for disease treatment. Here, as a step toward this goal for neurological diseases, we sequenced 12 postmortem, aged human frontal cortices (6 Alzheimer disease cases and 6 controls; 50% female) using one Oxford Nanopore PromethION flow cell per sample. We identified 1,917 medically relevant genes expressing multiple isoforms in the frontal cortex where 1,018 had multiple isoforms with different protein-coding sequences. Of these 1,018 genes, 57 are implicated in brain-related diseases including major depression, schizophrenia, Parkinson's disease and Alzheimer disease. Our study also uncovered 53 new RNA isoforms in medically relevant genes, including several where the new isoform was one of the most highly expressed for that gene. We also reported on five mitochondrially encoded, spliced RNA isoforms. We found 99 differentially expressed RNA isoforms between cases with Alzheimer disease and controls.

Metformin exposure during pregnancy and lactation affects offspring's long-term body weight and adipose tissue mass independent of the maternal metabolic state

The increasing prevalence of obesity, type 2 diabetes mellitus (T2DM), and gestational diabetes (GDM) among pregnant women has risen dramatically worldwide. The antihyperglycemic drug metformin is the most common drug for T2DM treatment in non-pregnant individuals; nevertheless, it is increasingly being used for diabetes-complicated pregnancies. Studies on the long-term metabolic effects of this drug in offspring remain scarce. This work aimed to determine the effect of metformin exposure during pregnancy and lactation on the offspring of a model of diet-induced maternal hyperglycemia.Cohorts of pregnant mice were fed a 46% fat diet (HFD) or a control standard diet (SD). A group of dams were exposed to metformin during pregnancy and lactation. After weaning, the offspring were fed SD for 8 weeks and then challenged with a 46% HFD after puberty for 12 weeks. Irrespective of the maternal diet, offspring of metformin-exposed mothers had a lower body weight and reduced inguinal white adipose tissue (iWAT) mass after HFD challenge. This was associated with increased expression of Pparg, Fabp4, Glut4, Srebp1, and Fasn in the iWAT during adulthood in the metabolically impaired dams exposed to metformin, suggesting increased adipogenesis and de novo lipogenesis. Increased expression of Fasn associated with decreased methylation levels at its promoter and proximal coding region in the iWAT was found. These results suggest that metformin modulates gene expression levels by epigenetic mechanisms in maternal metabolic-impaired conditions.

Two telomere-to-telomere gapless genomes reveal insights into Capsicum evolution and capsaicinoid biosynthesis

Chili pepper (Capsicum) is known for its unique fruit pungency due to the presence of capsaicinoids. The evolutionary history of capsaicinoid biosynthesis and the mechanism of their tissue specificity remain obscure due to the lack of high-quality Capsicum genomes. Here, we report two telomere-to-telomere (T2T) gap-free genomes of C. annuum and its wild nonpungent relative C. rhomboideum to investigate the evolution of fruit pungency in chili peppers. We precisely delineate Capsicum centromeres, which lack high-copy tandem repeats but are extensively invaded by CRM retrotransposons. Through phylogenomic analyses, we estimate the evolutionary timing of capsaicinoid biosynthesis. We reveal disrupted coding and regulatory regions of key biosynthesis genes in nonpungent species. We also find conserved placenta-specific accessible chromatin regions, which likely allow for tissue-specific biosynthetic gene coregulation and capsaicinoid accumulation. These T2T genomic resources will accelerate chili pepper genetic improvement and help to understand Capsicum genome evolution.

A screen for Plasmodium falciparum sporozoite surface protein binding to human hepatocyte surface receptors identifies novel host–pathogen interactions

BackgroundSporozoite invasion of hepatocytes is an essential step in the Plasmodium life-cycle and has similarities, at the cellular level, to merozoite invasion of erythrocytes. In the case of the Plasmodium blood-stage, efforts to identify host–pathogen protein–protein interactions have yielded important insights including vaccine candidates. In the case of sporozoite-hepatocyte invasion, the host–pathogen protein–protein interactions involved are poorly understood.MethodsTo gain a better understanding of the protein–protein interaction between the sporozoite ligands and host receptors, a systematic screen was performed. The previous Plasmodium falciparum and human surface protein ectodomain libraries were substantially extended, resulting in the creation of new libraries comprising 88 P. falciparum sporozoite protein coding sequences and 182 sequences encoding human hepatocyte surface proteins. Having expressed recombinant proteins from these sequences, a plate-based assay was used, capable of detecting low affinity interactions between recombinant proteins, modified for enhanced throughput, to screen the proteins for interactions. The novel interactions identified in the screen were characterized biochemically, and their essential role in parasite invasion was further elucidated using antibodies and genetically manipulated Plasmodium parasites.ResultsA total of 7540 sporozoite-hepatocyte protein pairs were tested under conditions capable of detecting interactions of at least 1.2 µM KD. An interaction between the human fibroblast growth factor receptor 4 (FGFR4) and the P. falciparum protein Pf34 is identified and reported here, characterizing its affinity and demonstrating the blockade of the interaction by reagents, including a monoclonal antibody. Furthermore, further interactions between Pf34 and a second P. falciparum rhoptry neck protein, PfRON6, and between human low-density lipoprotein receptor (LDLR) and the P. falciparum protein PIESP15 are identified.Conditional genetic deletion confirmed the essentiality of PfRON6 in the blood-stage, consistent with the important role of this protein in parasite lifecycle. Pf34 was refractory to attempted genetic modification. Antibodies to Pf34 abrogated the interaction and had a modest effect upon sporozoite invasion into primary human hepatocytes.ConclusionPf34 and PfRON6 may be members of a functionally important invasion complex which could be a target for future interventions. The modified interaction screening assay, protein expression libraries and P. falciparum mutant parasites reported here may be a useful tool for protein interaction discovery and antigen candidate screening which could be of wider value to the scientific community.

Leaderless foot-and-mouth disease virus serotype O did not cause clinical disease and failed to establish a persistent infection in cattle

ABSTRACT The foot-and-mouth disease virus (FMDV) Leader proteinase Lpro inhibits host mRNA translation and blocks the interferon response which promotes viral survival. Lpro is not required for viral replication in vitro but serotype A FMDV lacking Lpro has been shown to be attenuated in cattle and pigs. However, it is not known, whether leaderless viruses can cause persistent infection in vivo after simulated natural infection and whether the attenuated phenotype is the same in other serotypes. We have generated an FMDV O/FRA/1/2001 variant lacking most of the Lpro coding region (ΔLb). Cattle were inoculated intranasopharyngeally and observed for 35 days to determine if O FRA/1/2001 ΔLb is attenuated during the acute phase of infection and whether it can maintain a persistent infection in the upper respiratory tract. We found that although this leaderless virus can replicate in vitro in different cell lines, it is unable to establish an acute infection with vesicular lesions and viral shedding nor is it able to persistently infect bovine pharyngeal tissues.

TEX13B is essential for metabolic reprogramming during germ cell differentiation.

What is the functional significance of Tex13b in male germ cell development and differentiation? Tex13b regulates male germ cell differentiation by metabolic reprogramming during spermatogenesis. Studies in mice and humans suggest that TEX13B is a transcription factor and is exclusively expressed in germ cells. We sequenced the coding regions of TEX13B in 628 infertile men and 427 ethnically matched fertile control men. Further, to identify the molecular function of Tex13b, we created a Tex13b knockout and conditional overexpression system in GC-1spg (hereafter, GC-1) cells. Our recent exome sequencing study identified novel candidate genes for male infertility. TEX13B was found to be one of the potential candidates, hence we explored the role of TEX13B in male infertility within a large infertile case-control cohort. We performed functional analyses of Tex13b in a GC-1 cell line using CRISPR-Cas9. We differentially labelled the cell proteins by stable isotope labelling of amino acids in cell culture (SILAC) and performed mass spectrometry-based whole-cell proteomics to identify the differential protein regulation in knockout cells compared to wild-type cells. We found that Tex13b knockout leads to downregulation of the OXPHOS complexes and upregulation of glycolysis genes, which was further validated by western blotting. These results were further confirmed by respirometry analysis in Tex13b knockout cells. Further, we also performed a conditional overexpression of TEX13B in GC-1 cells and studied the expression of OXPHOS complex proteins by western blotting. We identified a rare variant, rs775429506 (p.Gly237Glu), exclusively in two non-obstructive-azoospermia (NOA) men, that may genetically predispose these men for infertility. Further, we demonstrated that Tex13b functions in the transcription regulation of OXPHOS complexes. N/A. We examined the function of Tex13b in GC-1 in vitro by knocking out and conditional overexpression, for understanding the function of Tex13b in germ cells. Unfortunately, this could not be replicated in either an animal model or in patient-derived tissue due to the non-availability of an animal model or patient's testis biopsies. This study identified that Tex13b plays an important role in male germ cell development and differentiation. The findings of this study would be useful for screening infertile males with spermatogenic failure and counselling them before the implementation of assisted reproduction technique(s). Funding was provided by the Council of Scientific and Industrial Research (CSIR) under the network project (BSC0101 and MLP0113) and SERB, the Department of Science and Technology, Government of India (J C Bose Fellowship: JCB/2019/000027). The authors do not have any competing interest.

A novel image encryption algorithm based on hyperchaotic system with cross-feedback structure and diffusive DNA coding operations

A novel image encryption algorithm based on hyperchaotic system with cross-feedback structure and diffusive DNA coding operations

Comparative genomic analyses provide insight into the pathogenicity of three Pseudomonas syringae pv. actinidiae strains from Anhui Province, China

BackgroundPseudomonas syringae pv. actinidiae (Psa) is an important bacterial plant pathogen that causes severe damage to the kiwifruit industry worldwide. Three Psa strains were recently obtained from different kiwifruit orchards in Anhui Province, China. The present study mainly focused on the variations in virulence and genome characteristics of these strains based on the pathogenicity assays and comparative genomic analyses.ResultsThree strains were identified as biovar 3 (Psa3), along with strain QSY6 showing higher virulence than JZY2 and YXH1 in pathogenicity assays. The whole genome assembly revealed that each of the three strains had a circular chromosome and a complete plasmid. The chromosome sizes ranged from 6.5 to 6.6 Mb with a GC content of approximately 58.39 to 58.46%, and a predicted number of protein-coding sequences ranging from 5,884 to 6,019. The three strains clustered tightly with 8 Psa3 reference strains in terms of average nucleotide identity (ANI), whole-genome-based phylogenetic analysis, and pangenome analysis, while they were evolutionarily distinct from other biovars (Psa1 and Psa5). Variations were observed in the repertoire of effectors of the type III secretion system among all 15 strains. Moreover, synteny analysis of the three sequenced strains revealed eight genomic regions containing 308 genes exclusively present in the highly virulent strain QSY6. Further investigation of these genes showed that 16 virulence-related genes highlight several key factors, such as effector delivery systems (type III secretion systems) and adherence (type IV pilus), which might be crucial for the virulence of QSY6.ConclusionThree Psa strains were identified and showed variant virulence in kiwifruit plant. Complete genome sequences and comparative genomic analyses further provided a theoretical basis for the potential pathogenic factors responsible for kiwifruit bacterial canker.

Knockdown of PHOX2B in the retrotrapezoid nucleus reduces the central CO2 chemoreflex in rats

PHOX2B is a transcription factor essential for the development of different classes of neurons in the central and peripheral nervous system. Heterozygous mutations in the PHOX2B coding region are responsible for the occurrence of Congenital Central Hypoventilation Syndrome (CCHS), a rare neurological disorder characterised by inadequate chemosensitivity and life-threatening sleep-related hypoventilation. Animal studies suggest that chemoreflex defects are caused in part by the improper development or function of PHOX2B expressing neurons in the retrotrapezoid nucleus (RTN), a central hub for CO2 chemosensitivity. Although the function of PHOX2B in rodents during development is well established, its role in the adult respiratory network remains unknown. In this study, we investigated whether reduction in PHOX2B expression in chemosensitive neuromedin-B (NMB) expressing neurons in the RTN altered respiratory function. Four weeks following local RTN injection of a lentiviral vector expressing the short hairpin RNA (shRNA) targeting Phox2b mRNA, a reduction of PHOX2B expression was observed in Nmb neurons compared to both naive rats and rats injected with the non-target shRNA. PHOX2B knockdown did not affect breathing in room air or under hypoxia, but ventilation was significantly impaired during hypercapnia. PHOX2B knockdown did not alter Nmb expression but it was associated with reduced expression of both Task2 and Gpr4, two CO2/pH sensors in the RTN. We conclude that PHOX2B in the adult brain has an important role in CO2 chemoreception and reduced PHOX2B expression in CCHS beyond the developmental period may contribute to the impaired central chemoreflex function.

Acquisition of different transcriptional shear mRNA and biological function of porcine interleukin 18 binding protein in PRRSV infection.

Interleukin-18 binding protein (IL-18BP), a natural regulator molecule of the pro-inflammatory cytokine interleukin-18 (IL-18), plays an important role in regulating the expression of the cellular immunity factor interferon-γ (IFN-γ). In a previous RNA-seq analysis of porcine alveolar macrophages (PAM) infected with the TIM and TJ strains of porcine reproductive and respiratory syndrome virus (PRRSV), we unexpectedly found that the mRNA expression of porcine interleukin 18-binding protein (pIL-18BP) in PAM cells infected with the TJM strain was significantly higher than that infected with the TJ strain. Studies have shown that human interleukin-18 binding protein (hIL-18bp) plays an important role in regulating cellular immunity in the course of the disease. However, there is a research gap on pIL-18BP. At the same time, PRRSV infection in pigs triggers weak cellular immune response problems. To explore the expression and the role of pIL-18BP in the cellular immune response induced by PRRSV, we strived to acquire the pIL-18BP gene from PAM or peripheral blood mononuclear cell (PBMC) with RT-PCR and sequencing. Furthermore, pIL-18BP and pIL-18 were both expressed prokaryotically and eukaryotically. The colocalization and interaction based on recombinant pIL-18BP and pIL-18 on cells were confirmed in vitro. Finally, the expression of pIL-18BP, pIL-18, and pIFN-γ was explored in pigs with different PRRSV infection states to interpret the biological function of pIL-18BP in vivo. The results showed there were five shear mutants of pIL-18BP. The mutant with the longest coding region was selected for subsequent functional validation. First, it was demonstrated that TJM-induced pIL-18BP mRNA expression was higher than that of TJ. A direct interaction between pIL-18BP and pIL-18 was confirmed through fluorescence colocalization, bimolecular fluorescent complimentary (BIFC), and co-immunoprecipitation (CO-IP). pIL-18BP also can regulate pIFN-γ mRNA expression. Finally, the expression of pIL-18BP, pIL-18, and pIFN-γ was explored in different PRRSV infection states. Surprisingly, both mRNA and protein expression of pIL-18 were suppressed. These findings fill the gap in understanding the roles played by pIL-18BP in PRRSV infection and provide a foundation for further research.IMPORTANCEPRRSV-infected pigs elicit a weak cellular immune response and the mechanisms of cellular immune regulation induced by PRRSV have not yet been fully elucidated. In this study, we investigated the role of pIL-18BP in PRRSV-induced immune response referring to the regulation of human IL-18BP to human interferon-gamma (hIFN-γ). This is expected to be used as a method to enhance the cellular immune response induced by the PRRSV vaccine. Here, we mined five transcripts of the pIL-18BP gene and demonstrated that it interacts with pIL-18 and regulates pIFN-γ mRNA expression. Surprisingly, we also found that both mRNA and protein expression of pIL-18 were suppressed under different PRRSV strains of infection status. These results have led to a renewed understanding of the roles of pIL-18BP and pIL-18 in cellular immunity induced by PRRSV infection, which has important implications for the prevention and control of PRRS.

Knockdown of PHOX2B in the retrotrapezoid nucleus reduces the central CO2 chemoreflex in rats.

PHOX2B is a transcription factor essential for the development of different classes of neurons in the central and peripheral nervous system. Heterozygous mutations in the PHOX2B coding region are responsible for the occurrence of Congenital Central Hypoventilation Syndrome (CCHS), a rare neurological disorder characterised by inadequate chemosensitivity and life-threatening sleep-related hypoventilation. Animal studies suggest that chemoreflex defects are caused in part by the improper development or function of PHOX2B expressing neurons in the retrotrapezoid nucleus (RTN), a central hub for CO2 chemosensitivity. Although the function of PHOX2B in rodents during development is well established, its role in the adult respiratory network remains unknown. In this study, we investigated whether reduction in PHOX2B expression in chemosensitive neuromedin-B (NMB) expressing neurons in the RTN altered respiratory function. Four weeks following local RTN injection of a lentiviral vector expressing the short hairpin RNA (shRNA) targeting Phox2b mRNA, a reduction of PHOX2B expression was observed in Nmb neurons compared to both naive rats and rats injected with the non-target shRNA. PHOX2B knockdown did not affect breathing in room air or under hypoxia, but ventilation was significantly impaired during hypercapnia. PHOX2B knockdown did not alter Nmb expression but it was associated with reduced expression of both Task2 and Gpr4, two CO2/pH sensors in the RTN. We conclude that PHOX2B in the adult brain has an important role in CO2 chemoreception and reduced PHOX2B expression in CCHS beyond the developmental period may contribute to the impaired central chemoreflex function.

Ribosome decision graphs for the representation of eukaryotic RNA translation complexity.

The application of ribosome profiling has revealed an unexpected abundance of translation in addition to that responsible for the synthesis of previously annotated protein-coding regions. Multiple short sequences have been found to be translated within single RNA molecules, within both annotated protein-coding and noncoding regions. The biological significance of this translation is a matter of intensive investigation. However, current schematic or annotation-based representations of mRNA translation generally do not account for the apparent multitude of translated regions within the same molecules. They also do not take into account the stochasticity of the process that allows alternative translations of the same RNA molecules by different ribosomes. There is a need for formal representations of mRNA complexity that would enable the analysis of quantitative information on translation and more accurate models for predicting the phenotypic effects of genetic variants affecting translation. To address this, we developed a conceptually novel abstraction that we term ribosome decision graphs (RDGs). RDGs represent translation as multiple ribosome paths through untranslated and translated mRNA segments. We termed the latter "translons." Nondeterministic events, such as initiation, reinitiation, selenocysteine insertion, or ribosomal frameshifting, are then represented as branching points. This representation allows for an adequate representation of eukaryotic translation complexity and focuses on locations critical for translation regulation. We show how RDGs can be used for depicting translated regions and for analyzing genetic variation and quantitative genome-wide data on translation for characterization of regulatory modulators of translation.

Antisense RNA regulates glutamine synthetase in a heterocyst-forming cyanobacterium.

Glutamine synthetase (GS) is a key enzyme involved in nitrogen assimilation and the maintenance of C/N balance, and it is strictly regulated in all bacteria. In cyanobacteria, GS expression is controlled by nitrogen control A (NtcA) transcription factor, which operates global nitrogen regulation in these photosynthetic organisms. Furthermore, posttranslational regulation of GS is operated by protein-protein interaction with GS inactivating factors (IFs). In this study, we describe an additional regulatory mechanism involving an antisense RNA. In Nostoc sp. PCC 7120, the gifA gene (encoding GS inactivating factor IF7) is transcribed downstream of the GS (glnA) gene, from the opposite strand, and the gifA mRNA extends into the glnA coding sequence in antisense orientation. Therefore, the dual RNA transcript that encodes gifA constitutes two functional regions: a 5' protein-coding region, encoding IF7, and a 3' untranslated region that acts as an antisense to glnA. By increasing the levels of such antisense RNA either in cis or in trans, we demonstrate that the amount of GS activity can be modulated by the presence of the antisense RNA. The tail-to-tail disposition of the glnA and gifA genes observed in many cyanobacterial strains from the Nostocales clade suggests the prevalence of such antisense RNA-mediated regulation of GS in this group of cyanobacteria.

Precise correction of a spectrum of β-thalassemia mutations in coding and non-coding regions by base editors

β-thalassemia/HbE results from mutations in the β-globin locus that impede the production of functional adult hemoglobin. Base editors (BEs) could facilitate the correction of the point mutations with minimal or no indel creation, but its efficiency and bystander editing for the correction of β-thalassemia mutations in coding and non-coding regions remains unexplored. Here, we screened BE variants in HUDEP-2 cells for their ability to correct a spectrum of β-thalassemia mutations that were integrated into the genome as fragments ofHBB. The identified targets were introduced into their endogenous genomic location using BEs and Cas9/homology-directed repair (HDR) to create cellular models with β-thalassemia/HbE. These β-thalassemia/HbE models were then used to assess the efficiency of correction in the native locus and functional β-globin restoration. Most bystander edits produced near target sites did not interfere with adult hemoglobin expression and are not predicted to be pathogenic. Further, the effectiveness of BE was validated for the correction of the pathogenic HbE variant in severe β0/βE-thalassaemia patient cells. Overall, our study establishes a novel platform to screen and select optimal BE tools for therapeutic genome editing by demonstrating the precise, efficient, and scarless correction of pathogenic point mutations spanning multiple regions of HBB including the promoter, intron, and exons.

Complete genome sequence of Paraburkholderia terrae strain KU-46, a 2,4-dinitrophenol-degrading bacterium.

Paraburkholderia terrae strain KU-46 has been studied for its capability to degrade 2,4-dinitrophenol. Here, we present the complete 10,833,180bp genome of this microorganism, comprising five circular chromosomes housing 9,797 protein-coding sequences. The genes responsible for 2,4-dinitrophenol and 4-nitrophenol degradation are located on chromosome 2.

Molecular Evolution of RAMOSA1 (RA1) in Land Plants.

RAMOSA1 (RA1) is a Cys2-His2-type (C2H2) zinc finger transcription factor that controls plant meristem fate and identity and has played an important role in maize domestication. Despite its importance, the origin of RA1 is unknown, and the evolution in plants is only partially understood. In this paper, we present a well-resolved phylogeny based on 73 amino acid sequences from 48 embryophyte species. The recovered tree topology indicates that, during grass evolution, RA1 arose from two consecutive SUPERMAN duplications, resulting in three distinct grass sequence lineages: RA1-like A, RA1-like B, and RA1; however, most of these copies have unknown functions. Our findings indicate that RA1 and RA1-like play roles in the nucleus despite lacking a traditional nuclear localization signal. Here, we report that copies diversified their coding region and, with it, their protein structure, suggesting different patterns of DNA binding and protein-protein interaction. In addition, each of the retained copies diversified regulatory elements along their promoter regions, indicating differences in their upstream regulation. Taken together, the evidence indicates that the RA1 and RA1-like gene families in grasses underwent subfunctionalization and neofunctionalization enabled by gene duplication.

The complete chloroplast genome of Durio zibethinus L. cultivar Ri6 (Helicteroideae, Malvaceae)

Durian, a member of the Malvaceae family, is famous for its delicious fruits, which have strong scents and are rich in nutrients. In this study, we sequenced and characterized the complete chloroplast genome of Durio zibethinus L. 1774 cultivar Ri6, a popular durian cultivar in Vietnam, using the Illumina Hiseq platform. The results showed a circular chloroplast genome composed of a large single copy of 96,115 bp, a small single copy of 20,819 bp, and two inverted repeat regions of 24,185 bp. This genome consisted of 79 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. The overall GC content of this genome was 35.7%. Phylogenetic analysis inferred from 78 protein-coding regions revealed monophyly of Durio species and a close relationship between D. zibethinus cultivar Ri6 and cultivar Mongthong. This study provides essential information for further studies examining genetic population, breedings, and species identification among Durio taxa and cultivars.

Increased ANKRD1 Levels in Early Senescence Mediated by RBMS1-Elicited ANKRD1 mRNA Stabilization.

Cellular senescence is a dynamic biological process triggered by sublethal cell damage and driven by specific changes in gene expression programs. We recently identified ANKRD1 (ankyrin repeat domain 1) as a protein strongly elevated after triggering senescence in fibroblasts. Here, we set out to investigate the mechanisms driving the elevated production of ANKRD1 in the early stages of senescence. Our results indicated that the rise in ANKRD1 levels after triggering senescence using etoposide (Eto) was the result of moderate increases in transcription and translation, and robust mRNA stabilization. Antisense oligomer (ASO) pulldown followed by mass spectrometry revealed a specific interaction of the RNA-binding protein RBMS1 with ANKRD1 mRNA that was confirmed by ribonucleoprotein immunoprecipitation analysis. RBMS1 abundance decreased in the nucleus and increased in the cytoplasm during Eto-induced senescence; in agreement with the hypothesis that RBMS1 may participate in post-transcriptional stabilization of ANKRD1 mRNA, silencing RBMS1 reduced, while overexpressing RBMS1 enhanced ANKRD1 mRNA half-life after Eto treatment. A segment proximal to the ANKRD1 coding region was identified as binding RBMS1 and conferring RBMS1-dependent increased expression of a heterologous reporter. We propose that RBMS1 increases expression of ANKRD1 during the early stages of senescence by stabilizing ANKRD1 mRNA.