Search

To understand the mitochondrial genome structure of two endangered and long-armed scarab beetles, Propomacrus davidi and Propomacrus bimucronatus, their complete mitogenomes were sequenced for the first time in this study. The complete mitogenomes of P. davidi and P. bimucronatus were 18, 042 bp and 18, 104 bp in length, respectively. The gene orders of their mitogenomes were highly consistent with other Coleopteran species, and the typical ATN was used as the start codon in most protein coding genes. The incomplete stop codon T was used in cox1, cox2, and nad5, and TAN was used as a complete stop codon in most protein coding genes. All predicted tRNAs could form a typical cloverleaf secondary structure, except that trnS1 lacked the dihydrouridine arm. Based on the maximum likelihood and the Bayesian inference methods, phylogenetic trees of 50 species were reconstructed. The results showed that P. davidi, P. bimucronatus, Cheirotonus jansoni and Cheirotonus gestroi clustered in the same branch, and were the most closely related. The results supported that subfamily Euchirinae is a monophyletic group of Scarabaeidae, which was consistent with the morphological classification. These molecular data enriched the complete mitogenome database of Euchirinae, and improved our understanding of the phylogenetic relationship and evolutionary characteristics of these two endangered species. [ABSTRACT FROM AUTHOR]

Sex-determining mechanism is still ambiguous for sea cucumber Holothuria scabra which only manifests gonochorism in gonad. In this study, proteomic analysis was employed to delineate sex-related proteins and genes in gonads of H. scabra, subsequently validated through Quantitative real-time polymerase chain reaction (qRT-PCR). A total of 5,313 proteins were identified via proteome sequencing. Among these, 817 proteins exhibited expression in both the ovary and testis, with 445 proteins displaying up-regulation and 372 proteins showing down-regulation (ovary vs testis). Furthermore, 136 and 69 proteins were identified as ovary-specific and testis-specific Differentially Abundant Proteins (DAPs), respectively. And 9 DAP coding genes which play crucial role in ovary and testis were verified by qRT-PCR. Notably, 24 ovary-bias proteins enriched in ribosome pathway strongly indicated the crucial role of ribosome in ovary. This study serves to furnish novel evidence pertaining to sex differences in H. scabra. [ABSTRACT FROM AUTHOR]

Abstract algebraic concepts such as category are considered cornerstones on which logical consistency relies in any sophisticated study of natural phenomena. However, to the best of our knowledge, in molecular/genetic biology, their application is still severely limited because they capture neither the dynamics nor provide a visual form. The Petri net (PN) has often been used to illustrate visually parallel, asynchronous dynamic events in small data systems. A prototypal hybrid model combining both category theory and extended PNs may instead be indispensable for that purpose. This hybrid model incorporates 1) token-like elements of a group, 2) object-like places of a category, 3) square poles (rather than pentagon poles) that enable unique identifications of single-strand DNA sequences from the shape of its polygonal line, 4) creation/annihilation morphisms that generate/erase tokens, 5) Cartesian products 'Z5×Z2×...' that enable conversions between DNA and RNA sequences, 6) somatic recombinations (VDJ recombinations) for antibodies displayed concretely in category-theoretic form, 7) 'identity protein Δ' translated from a triplet of identity bases 'EEE' as an advanced concept from our previous display of the canonical central dogma, 8) illustrations of an incidence-matrix-like matrix A that includes operators as coordinates, and 9) basic topics concerning the canonical central dogma being displayed concretely using concepts of conventional category theory such as 'adjoint', 'adjoint functor', 'natural transformation', 'Yoneda's lemma' and 'Kan extension'. These ideas provide more advanced tools that expand our previous model concerning nucleic-acid-base sequences. Despite the nascent nature of our methodology, our hybrid model has potential in a variety of applications, illustrated using molecular/genetic sequences, in particular providing a simple dynamic/visual representation. With further improvements, this approach may prove effective in reducing the need for large data-storing systems. [ABSTRACT FROM AUTHOR]

In embryonic development and organogenesis, cells sharing identical genetic codes acquire diverse gene expression states in a highly reproducible spatial distribution, crucial for multicellular formation and quantifiable through positional information. To understand the spontaneous growth of complexity, we constructed a one-dimensional division-decision model, simulating the growth of cells with identical genetic networks from a single cell. Our findings highlight the pivotal role of cell division in providing positional cues, escorting the system toward states rich in information. Moreover, we pinpointed lateral inhibition as a critical mechanism translating spatial contacts into gene expression. Our model demonstrates that the spatial arrangement resulting from cell division, combined with cell lineages, imparts positional information, specifying multiple cell states with increased complexity—illustrated through examples in C.elegans. This study constitutes a foundational step in comprehending developmental intricacies, paving the way for future quantitative formulations to construct synthetic multicellular patterns. Author summary: Embryonic development shapes our bodies from a single cell, determining the placement of the head and tail. But how do cells, all sharing the same genetic code, precisely know what to become? Our mathematical model cracks this code. Envision your information being provided by your neighbors and your 'mom'—the division lineage. Then, appropriate regulatory networks (such as lateral inhibition, the most effective network motif) transform this information into diverse yet robust gene expressions. This math model helps us see the rules behind spontaneously growing complexity, guiding us to create new patterns of cells. It's a big step in understanding how bodies form and opens doors to building cool new structures. [ABSTRACT FROM AUTHOR]

The standard genetic code defines the rules of translation for nearly every life form on Earth. It also determines the amino acid changes accessible via single-nucleotide mutations, thus influencing protein evolvability—the ability of mutation to bring forth adaptive variation in protein function. One of the most striking features of the standard genetic code is its robustness to mutation, yet it remains an open question whether such robustness facilitates or frustrates protein evolvability. To answer this question, we use data from massively parallel sequence-to-function assays to construct and analyze 6 empirical adaptive landscapes under hundreds of thousands of rewired genetic codes, including those of codon compression schemes relevant to protein engineering and synthetic biology. We find that robust genetic codes tend to enhance protein evolvability by rendering smooth adaptive landscapes with few peaks, which are readily accessible from throughout sequence space. However, the standard genetic code is rarely exceptional in this regard, because many alternative codes render smoother landscapes than the standard code. By constructing low-dimensional visualizations of these landscapes, which each comprise more than 16 million mRNA sequences, we show that such alternative codes radically alter the topological features of the network of high-fitness genotypes. Whereas the genetic codes that optimize evolvability depend to some extent on the detailed relationship between amino acid sequence and protein function, we also uncover general design principles for engineering nonstandard genetic codes for enhanced and diminished evolvability, which may facilitate directed protein evolution experiments and the bio-containment of synthetic organisms, respectively. Only "one in a million" possible genetic codes are as robust as the standard genetic code, but does this robustness accelerate or impede adaptive evolution? This study uses experimental data from six massively parallel sequence-to-function assays for four proteins to show that robust genetic codes enhance protein evolvability by producing smooth adaptive landscapes with few peaks. [ABSTRACT FROM AUTHOR]

The relationship between genetic code robustness and protein evolvability is unknown. A new study in PLOS Biology using in silico rewiring of genetic codes and functional protein data identified a positive correlation between code robustness and protein evolvability that is protein-specific. The relationship between genetic code robustness and protein evolvability is unknown. This Primer explores a new PLOS Biology study which uses in silico rewiring of genetic codes and functional protein data to identify a positive correlation between code robustness and protein evolvability that is protein-specific. [ABSTRACT FROM AUTHOR]

In this study, an extensive analysis of microsatellite markers (Single Tandem Repeats—STRs) in Penaeus vannamei was conducted at an advanced level. The markers were thoroughly examined, characterized, and specific markers located within coding regions were identified. Out of a total of 306 STRs, 117 were classified as perfect markers based on their single repeat motif. Among these perfect markers, 62 were found to be associated with predicted coding genes (mRNA), which were involved in various functions such as binding, catalytic activity, ATP-dependent activity, transcription, structural and molecular regulation. To validate the accuracy of the findings, a sample of nine markers was subjected to in vitro testing, which confirmed the presence of polymorphisms within the population. These results suggest the existence of different protein isoforms within the population, indicating the potential of these markers for application in both population and phenotype-genotype association studies. This innovative approach opens up new possibilities for investigating the impact of genomic plasticity in populations of P. vannamei. [ABSTRACT FROM AUTHOR]

Influenza A virus contains regions of its segmented genome associated with ability to package the segments into virions, but many such regions are poorly characterised. We provide detailed predictions of the key locations within these packaging-associated regions, and their structures, by applying a recently-improved pipeline for delineating constrained regions in RNA viruses and applying structural prediction algorithms. We find and characterise other known constrained regions within influenza A genomes, including the region associated with the PA-X frameshift, regions associated with alternative splicing, and constraint around the initiation motif for a truncated PB1 protein, PB1-N92, associated with avian viruses. We further predict the presence of constrained regions that have not previously been described. The extra characterisation our work provides allows investigation of these key regions for drug target potential, and points towards determinants of packaging compatibility between segments. Author summary: Influenza A viruses infect birds and mammals, including humans. Their genetic code is contained in multiple segments; all these segments need to be packaged inside virus particles for effective transmission. Mixing of segments between viruses that infect different organisms can result in pandemics. We have developed a computational pipeline for finding regions within the genetic code where mutation is less often seen; these regions are often important for virus lifecycle. In the case of influenza A, many of the regions we identify are associated with ensuring the segments are efficiently packaged. Other regions are likely to be important in allowing the virus to optimise the ratios of the different proteins it produces. Our work guides how to investigate these regions further for drug target potential. [ABSTRACT FROM AUTHOR]

Processive enzymes like polymerases or ribosomes are often studied in bulk experiments by monitoring time-dependent signals, such as fluorescence time traces. However, due to biomolecular process stochasticity, ensemble signals may lack the distinct features of single-molecule signals. Here, we demonstrate that, under certain conditions, bulk signals from processive reactions can be decomposed to unveil hidden information about individual reaction steps. Using mRNA translation as a case study, we show that decomposing a noisy ensemble signal generated by the translation of mRNAs with more than a few codons is an ill-posed problem, addressable through Tikhonov regularization. We apply our method to the fluorescence signatures of in-vitro translated LepB mRNA and determine codon-position dependent translation rates and corresponding state-specific fluorescence intensities. We find a significant change in fluorescence intensity after the fourth and the fifth peptide bond formation, and show that both codon position and encoded amino acid have an effect on the elongation rate. This demonstrates that our approach enhances the information content extracted from bulk experiments, thereby expanding the range of these time- and cost-efficient methods. Author summary: Ensemble (bulk) in-vitro experiments are less complex and expensive than single-molecule measurements, but their read-outs are blurred due to the stochastic nature of biomolecular processes. Our study investigates how information about individual steps in processive reactions can be deduced from ensemble experiments, despite these inherent challenges. We introduce a least-squares method for decomposing ensemble signals from synchronized enzymatic reactions, allowing us to reveal distinct features of individual reaction steps. Using mRNA translation as a showcase, we propose a systematic approach to assess the feasibility of ensemble signal decomposition. To demonstrate the efficacy and limits of our approach, we decompose simulated fluorescence signatures of in-vitro mRNA translation processes. We then apply our method to experimental data of LepB mRNA translation, where a codon-by-codon analysis reveals that both the codon position and the encoded amino acid impact the translation dynamics. Our approach can be applied to assess bulk signals of other processive enzymatic reactions beyond mRNA translation, provided a comprehensive quantitative understanding of the process is available. The introduced method offers a robust means to extract detailed information from bulk experimental setups, thus helping to bridge the gap between cost-effective bulk experiments and detailed single-molecule studies. [ABSTRACT FROM AUTHOR]

A number of studies have demonstrated that epigenetic factors regulate plant developmental timing in response to environmental changes. However, we still have an incomplete view of how epigenetic factors can regulate developmental events such as organogenesis, and the transition from cell division to cell expansion, in plants. The small number of cell types and the relatively simple developmental progression required to form the Arabidopsis petal makes it a good model to investigate the molecular mechanisms driving plant organogenesis. In this study, we investigated how the RABBIT EARS (RBE) transcriptional repressor maintains the downregulation of its downstream direct target, TCP5, long after RBE expression dissipates. We showed that RBE recruits the Groucho/Tup1-like corepressor TOPLESS (TPL) to repress TCP5 transcription in petal primordia. This process involves multiple layers of changes such as remodeling of chromatin accessibility, alteration of RNA polymerase activity, and histone modifications, resulting in an epigenetic memory that is maintained through multiple cell divisions. This memory functions to maintain cell divisions during the early phase of petal development, and its attenuation in a cell division-dependent fashion later in development enables the transition from cell division to cell expansion. Overall, this study unveils a novel mechanism by which the memory of an epigenetic state, and its cell-cycle regulated decay, acts as a timer to precisely control organogenesis. Author summary: Epigenetic changes, such as histone modifications or DNA methylation, can modify gene expression without affecting the genetic code. Such alterations can persist for many cell generations, and often are induced by various environmental inputs in both plants and animals. By contrast, little is known of how intrinsic epigenetic changes can act to regulate the temporal progression of development. Here we showed that the expression of RABBIT EARS (RBE) represses the transcription of its target, TCP5, through inducing a variety of epigenetic changes at the TCP5 locus. These changes dissipated over a period of days, allowing for eventual activation of TCP5 expression and a concomitant shift from cell division to cell expansion during Arabidopsis petal development. The gradual loss of repressive epigenetic marks at the TCP5 locus represents an internal timing mechanism to temporally control the development of the Arabidopsis petal. [ABSTRACT FROM AUTHOR]

Understanding air quality requires a comprehensive understanding of its various factors. Most of the association rule techniques focuses on high frequency terms, ignoring the potential importance of low- frequency terms and causing unnecessary storage space waste. Therefore, a dynamic genetic association rule mining algorithm is proposed in this paper, which combines the improved dynamic genetic algorithm with the association rule mining algorithm to realize the importance mining of low- frequency terms. Firstly, in the chromosome coding phase of genetic algorithm, an innovative multi-information coding strategy is proposed, which selectively stores similar values of different levels in one storage unit. It avoids storing all the values at once and facilitates efficient mining of valid rules later. Secondly, by weighting the evaluation indicators such as support, confidence and promotion in association rule mining, a new evaluation index is formed, avoiding the need to set a minimum threshold for high-interest rules. Finally, in order to improve the mining performance of the rules, the dynamic crossover rate and mutation rate are set to improve the search efficiency of the algorithm. In the experimental stage, this paper adopts the 2016 annual air quality data set of Beijing to verify the effectiveness of the unit point multi-information coding strategy in reducing the rule storage air, the effectiveness of mining the rules formed by the low frequency item set, and the effectiveness of combining the rule mining algorithm with the swarm intelligence optimization algorithm in terms of search time and convergence. In the experimental stage, this paper adopts the 2016 annual air quality data set of Beijing to verify the effectiveness of the above three aspects. The unit point multi-information coding strategy reduced the rule space storage consumption by 50%, the new evaluation index can mine more interesting rules whose interest level can be up to 90%, while mining the rules formed by the lower frequency terms, and in terms of search time, we reduced it about 20% compared with some meta-heuristic algorithms, while improving convergence. [ABSTRACT FROM AUTHOR]

Domestic dogs are the primary urban reservoirs of Leishmania infantum, the causative agent of visceral leishmaniasis. In Canine Leishmaniasis (CanL), modulation of the host's immune response may be associated with the expression of small non-coding RNAs called microRNA (miR). miR-194 expression increases in peripheral blood mononuclear cells (PBMCs) of dogs with leishmaniasis with a positive correlation with the parasite load and in silico analysis demonstrated that the TRAF6 gene is the target of miR-194 in PBMCs from diseased dogs. Here, we isolated PBMCs from 5 healthy dogs and 28 dogs with leishmaniasis, naturally infected with L. infantum. To confirm changes in miR-194 and TRAF6 expression, basal expression of miR-194 and gene expression of TRAF6 was measured using qPCR. PBMCs from healthy dogs and dogs with leishmaniasis were transfected with miR-194 scramble, mimic, and inhibitor and cultured at 37° C, 5% CO2 for 48 hours. The expression of possible targets was measured: iNOS, NO, T-bet, GATA3, and FoxP3 were measured using flow cytometry; the production of cytokines IL-1β, IL-4, IL-6, IL-10, TNF-α, IFN-γ, and TGF-β in cell culture supernatants was measured using capture enzyme-linked immunosorbent assays (ELISA). Parasite load was measured using cytometry and qPCR. Functional assays followed by miR-194 inhibitor and IL-1β blockade and assessment of NO production were also performed. Basal miR-194 expression was increased in PBMC from dogs with Leishmaniasis and was negatively correlated with TRAF6 expression. The mimic of miR-194 promoted an increase in parasite load. There were no significant changes in T-bet, GATA3, or FoxP3 expression with miR-194 enhancement or inhibition. Inhibition of miR-194 increased IL-1β and NO in PBMCs from diseased dogs, and blockade of IL-1β following miR-194 inhibition decreased NO levels. These findings suggest that miR-194 is upregulated in PBMCs from dogs with leishmaniasis and increases parasite load, possibly decreasing NO production via IL-1β. These results increase our understanding of the mechanisms of evasion of the immune response by the parasite and the identification of possible therapeutic targets. Author summary: Visceral leishmaniasis (VL) is a zoonosis in tropical and subtropical regions. In the new world, it is caused by the obligate intracellular protozoan Leishmania infantum. In urban areas, domestic dogs are the primary reservoir. Increase in canine leishmaniasis (CanL) contributes to the increase in VL in humans. One mechanism of parasite survival is through regulation of microRNA (miRNA), responsible for translational control of transcripts related to anti-parasite immunity. One example is miR-194, which is increased in the peripheral blood mononuclear cells (PBMCs) of dogs with leishmaniasis. In the present study, we found that increasing miR-194 increases parasite load, and decreasing miR-194 increases IL1β and NO. Subsequently, we found that the increase in NO depends on IL1β. Since NO its directly involved in parasite load control, these findings contribute to understanding how the parasite evades immune responses and paves the way for the discovery of new therapeutic targets. [ABSTRACT FROM AUTHOR]

Staphylococcus aureus is an opportunistic pathogen capable of causing many different human diseases. During colonization and infection, S. aureus will encounter a range of hostile environments, including acidic conditions such as those found on the skin and within macrophages. However, little is known about the mechanisms that S. aureus uses to detect and respond to low pH. Here, we employed a transposon sequencing approach to determine on a genome-wide level the genes required or detrimental for growth at low pH. We identified 31 genes that were essential for the growth of S. aureus at pH 4.5 and confirmed the importance of many of them through follow up experiments using mutant strains inactivated for individual genes. Most of the genes identified code for proteins with functions in cell wall assembly and maintenance. These data suggest that the cell wall has a more important role than previously appreciated in promoting bacterial survival when under acid stress. We also identified several novel processes previously not linked to the acid stress response in S. aureus. These include aerobic respiration and histidine transport, the latter by showing that one of the most important genes, SAUSA300_0846, codes for a previously uncharacterized histidine transporter. We further show that under acid stress, the expression of the histidine transporter gene is increased in WT S. aureus. In a S. aureus SAUSA300_0846 mutant strain expression of the histidine biosynthesis genes is induced under acid stress conditions allowing the bacteria to maintain cytosolic histidine levels. This strain is, however, unable to maintain its cytosolic pH to the same extent as a WT strain, revealing an important function specifically for histidine transport in the acid stress response of S. aureus. Author summary: Staphylococcus aureus is an important human bacterial pathogen that can cause a range of diseases. During infection, the pathogen will encounter a hostile environment within the human host, including acidic conditions such as those found on the skin and within macrophages. The bacterium has developed sophisticated strategies to survive and grow under such harsh conditions. Here we performed a genome wide screen to identify factors that are required by this pathogen to survive under acid stress conditions and identified several novel processes including histidine uptake. Understanding the response of S. aureus to deal with acid stress conditions will help us better manage infections. [ABSTRACT FROM AUTHOR]

The control of reproductive processes in Rhodnius prolixus involves a variety of neuroactive chemicals. Among these, several studies have suggested that the biogenic amine octopamine (OA), might play an active role in these processes. Here, we investigate the molecular profile of the R. prolixus α adrenergic-like OA receptor 1 (RpOAα1-R) and its role in egg production. Comparative molecular analyses confirm that the RpOAα1-R gene codes for a true OAα1 receptor. The RpOAα1-R transcript is highly expressed in tissues associated with egg production, and after a blood meal, which is the stimulus for full egg production in R. prolixus, the RpOAα1-R transcript is upregulated in the ovaries and spermatheca. After RNAi-mediated RpOAα1-R knockdown, an ovarian phenotype characterized by slow egg development is observed. Furthermore, an altered egg phenotype has been characterized with eggs that are deformed. Interestingly, there is no evidence of disruption in vitellogenin (Vg) synthesis by the fat body or uptake by the oocytes. On the other hand, RpOAα1-R downregulation is correlated with defective choriogenesis in the eggs. These results provide critical information concerning the role of OAα1-R in oogenesis in R. prolixus. [ABSTRACT FROM AUTHOR]

Bacterial micro-compartments (BMC) are complex macromolecular assemblies that participate in varied metabolic processes in about 20% of bacterial species. Most of these organisms carry BMC genetic information organized in operons that often include several paralog genes coding for components of the compartment shell. BMC shell constituents can be classified depending on their oligomerization state as hexamers (BMC-H), pentamers (BMC-P) or trimers (BMC-T). Formation of hetero-oligomers combining different protein homologs is theoretically feasible, something that could ultimately modify BMC shell rigidity or permeability, for instance. Despite that, it remains largely unknown whether hetero-oligomerization is a widespread phenomenon. Here, we demonstrated that the tripartite GFP (tGFP) reporter technology is an appropriate tool that might be exploited for such purposes. Thus, after optimizing parameters such as the size of linkers connecting investigated proteins to GFP10 or GFP11 peptides, the type and strength of promoters, or the impact of placing coding cassettes in the same or different plasmids, homo-oligomerization processes could be successfully monitored for any of the three BMC shell classes. Moreover, the screen perfectly reproduced published data on hetero-association between couples of CcmK homologues from Syn. sp. PCC6803, which were obtained following a different approach. This study paves the way for mid/high throughput screens to characterize the extent of hetero-oligomerization occurrence in BMC-possessing bacteria, and most especially in organisms endowed with several BMC types and carrying numerous shell paralogs. On the other hand, our study also unveiled technology limitations deriving from the low solubility of one of the components of this modified split-GFP approach, the GFP1-9. [ABSTRACT FROM AUTHOR]

In Trypanosoma cruzi DNA is packaged into chromatin by octamers of histone proteins that form nucleosomes. Transcription of protein coding genes in trypanosomes is constitutive producing polycistronic units and gene expression is primarily regulated post-transcriptionally. However, chromatin organization influences DNA dependent processes. Hence, determining nucleosome position is of uppermost importance to understand the peculiarities found in trypanosomes. To map nucleosomes genome-wide in several organisms, digestion of chromatin with micrococcal nuclease followed by deep sequencing has been applied. Nonetheless, the special requirements for cell manipulation and the uniqueness of the chromatin organization in trypanosomes entails a customized analytical approach. In this work, we adjusted this broadly used method to the hybrid reference strain, CL Brener. Particularly, we implemented an exhaustive and thorough computational workflow to overcome the difficulties imposed by this complex genome. We tested the performance of two aligners, Bowtie2 and HISAT2, and discuss their advantages and caveats. Specifically, we highlight the relevance of using the whole genome as a reference instead of the commonly used Esmeraldo-like haplotype to avoid spurious alignments. Additionally, we show that using the whole genome refines the average nucleosome representation, but also the quality of mapping for every region represented. Moreover, we show that the average nucleosome organization around trans-splicing acceptor site described before, is not just an average since the same chromatin pattern is detected for most of the represented regions. In addition, we extended the study to a non-hybrid strain applying the experimental and analytical approach to Sylvio-X10 strain. Furthermore, we provide a source code for the construction of 2D plots and heatmaps which are easy to adapt to any T. cruzi strain. [ABSTRACT FROM AUTHOR]

We analyzed the expression of genes coding for Na+ transporters (OsHKT1.5, OsHKT1.1, OsSOS1, OsSOS2, OsNHX1, OsNHX2), Cl- transporter (OsNRT1, OsCLC, OsCCC1) and gene coding for the transcription factor DREB (OsDREB2) involved in response to desiccation in two cultivars of O. glaberrrima differing in salt-resistance (salt-tolerant cultivar (TOG5307) and salt-sensitive (TOG 5949)) exposed to NaCl, PEG or both agents present simultaneously. Seedlings were grown in iso-osmotic nutrient solution (Ψs = -0.47±0.02 MPa) containing PEG 6,000 12.9% (water stress), NaCl 75 mM (salt stress) and PEG 6.4% + NaCl 37.5 mM (MIX-treatment) during 1 and 7 days. Plants were analyzed for gene expression, mineral nutrients, and photosynthetic-related parameters. Na+ and Cl- accumulations in salt-treated plants were lower in roots and shoots of TOG5307 comparatively to TOG5949 while water content decreased in TOG5307. TOG5307 exhibited tolerance to water stress and maintained higher net photosynthesis and water use efficiency than TOG5949 in response to all treatments, but was less efficient for osmotic adjustment. Dehydration tolerance of TOG5307 involves a higher OsDREB2 expression. TOG5307 also exhibited a higher OsSOS1, OsSOS2, OsNHX1 and OsNHX2 expression than TOG5949 in response to salinity. OsHKT1.5 was slightly induced in the shoot. OsHKT1.1 was recorded in the shoots but remained undetectable in the roots. Chloride and sodium accumulations were strongly reduced in the shoots when PEG was present. Salinity resistance in Oryza glaberrima implies tolerance to dehydration as well as complementary strategies of Na+ exclusion through the SOS system and Na+ tolerance through vacuolar sequestration. [ABSTRACT FROM AUTHOR]

Ribosomes are complex macromolecules assembled from 4 rRNAs and 79 ribosomal proteins (RPs). Their assembly is organized in a highly hierarchical manner, which is thought to avoid dead-end pathways, thereby enabling efficient assembly of ribosomes in the large quantities needed for healthy cellular growth. Moreover, hierarchical assembly also can help ensure that each RP is included in the mature ribosome. Nonetheless, how this hierarchy is achieved remains unknown, beyond the examples that depend on direct RP-RP interactions, which account for only a fraction of the observed dependencies. Using assembly of the small subunit head and a disease-associated mutation in the assembly factor Ltv1 as a model system, we dissect here how the hierarchy in RP binding is constructed. A combination of data from yeast genetics, mass spectrometry, DMS probing and biochemical experiments demonstrate that the LIPHAK-disease-associated Ltv1 mutation leads to global defects in head assembly, which are explained by direct binding of Ltv1 to 5 out of 15 RPs, and indirect effects that affect 4 additional RPs. These indirect effects are mediated by conformational transitions in the nascent subunit that are regulated by Ltv1. Mechanistically, Ltv1 aids the recruitment of some RPs via direct protein-protein interactions, but surprisingly also delays the recruitment of other RPs. Delayed binding of key RPs also delays the acquisition of RNA structure that is stabilized by these proteins. Finally, our data also indicate direct roles for Ltv1 in chaperoning the folding of a key rRNA structural element, the three-helix junction j34-35-38. Thus, Ltv1 plays critical roles in organizing the order of both RP binding to rRNA and rRNA folding, thereby enabling efficient 40S subunit assembly. Author summary: Ribosomes make proteins in all cells. During that process, they are responsible not just for the correct translation of the genetic code into amino acids, but also for selecting the correct mRNAs for translation. Finally, damage in mRNAs requires active translation. Changes in ribosome numbers or ribosome composition can perturb each of these functions and can lead to diseases. Thus, a better understanding of how cells ensure efficient and correct ribosome assembly is paramount. Using a mutation that is associated with a dermatological condition, LIPHAK syndrome, we describe here novel roles for the assembly factor Ltv1 in orchestrating the assembly of RPs into the head of the small ribosomal subunit. Our experiments, a combination of mass spectrometry, RNA structure probing, genetics and biochemical assays, demonstrate that Ltv1 helps recruit some RPs, while delaying the incorporation of others, thereby establishing the hierarchy of RP binding. In addition, our results also demonstrate roles for Ltv1 in modulating rRNA folding. Finally, the data establish the identity of a previously described quality control step that is bypassed by a cancer-associated mutation. [ABSTRACT FROM AUTHOR]

In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria. Author summary: Transgene experiments in Escherichia coli report that codon usage in the first few amino acids after the initiating ATG has a profound influence on the resulting protein level by promoting ribosomal initiation, rather than enabling speedy elongation. For the design of bacterial transgenes and for simple codon-based in silico inference of protein levels, it would be valuable to computationally determine initiation optimality scores for codons for any given species. An attractive approach is to characterize the 5' codon enrichment of highly expressed genes compared with the most lowly expressed genes, just as translational optimality scores of codons (affecting elongation) have been similarly defined employing the full gene body. Using unique resources provided for E. coli we show that, unexpectedly, this doesn't work: the 5' ends of highly expressed genes are enriched, compared to lowly expressed ones, in codons that objectively are associated with low initiation rates. This likely reflects conflicting selection pressures in highly expressed genes which can favour low initiation to promote low noise or high efficiency. While simple prediction of native protein abundance from 5' codon content is then somewhat futile, the objective initiation optimality scores may hold relevance for in silico transgene design for a broad spectrum of bacteria. [ABSTRACT FROM AUTHOR]

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the "universal" genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3'UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids. Author summary: The genetic code is almost universal across life. The vast majority of organisms use the canonical genetic code, which has three stop codons (UAA, UAG, and UGA) and 61 sense codons that code for amino acids. Here, we report the discovery of an unexpected genetic code variant in an uncultured ciliate species from the Oligohymenophorea class, where the canonical stop codons UAA and UAG have been reassigned to code for lysine and glutamic acid, respectively. This is a particularly unusual genetic code reassignment as UAA and UAG differ at the wobble position and their evolution is thought to be coupled. We also report that the remaining stop codon, UGA, is enriched immediately downstream of genes in the same reading frame, suggesting a possible role in minimising deleterious consequences in the event of translational readthrough. Our work documents, for the first time, a genetic code variant where the codons UAA and UAG specify two different amino acids and shows that there are still unexplored genetic code reassignments awaiting discovery. [ABSTRACT FROM AUTHOR]

The unknown pathogenicity of a significant number of variants found in cancer-related genes is attributed to limited epidemiological data, resulting in their classification as variant of uncertain significance (VUS). To date, Breast Cancer gene-2 (BRCA2) has the highest number of VUSs, which has necessitated the development of several robust functional assays to determine their functional significance. Here we report the use of a humanized-mouse embryonic stem cell (mESC) line expressing a single copy of the human BRCA2 for a CRISPR-Cas9-based high-throughput functional assay. As a proof-of-principle, we have saturated 11 codons encoded by BRCA2 exons 3, 18, 19 and all possible single-nucleotide variants in exon 13 and multiplexed these variants for their functional categorization. Specifically, we used a pool of 180-mer single-stranded donor DNA to generate all possible combination of variants. Using a high throughput sequencing-based approach, we show a significant drop in the frequency of non-functional variants, whereas functional variants are enriched in the pool of the cells. We further demonstrate the response of these variants to the DNA-damaging agents, cisplatin and olaparib, allowing us to use cellular survival and drug response as parameters for variant classification. Using this approach, we have categorized 599 BRCA2 variants including 93-single nucleotide variants (SNVs) across the 11 codons, of which 28 are reported in ClinVar. We also functionally categorized 252 SNVs from exon 13 into 188 functional and 60 non-functional variants, demonstrating that saturation genome editing (SGE) coupled with drug sensitivity assays can enhance functional annotation of BRCA2 VUS. Author summary: The exponential rise in genetic sequencing led to the identification of several missense variants. However, its clinical utility is often limited by our ability to determine the functional impact of genetic variants. This has resulted in the identification of many variant of uncertain significance (VUS), particularly in BRCA2. Several functional assays have been developed to ascertain the impact of VUSs on protein function that can be used to determine their pathogenicity. Our mouse ES cell (mESC)-based method includes generation of individual BRCA2 variants in a bacterial artificial chromosome (BAC) by recombineering and assess their expression in mouse ES cells. These cells can then be used to determine how the variants affect mESC viability and sensitivity to DNA damaging agents. Although the assay has high sensitivity and specificity, the process is time consuming and labor-intensive. To overcome these limitations, here we report the development of a CRISPR-Cas9 based high-throughput approach that can be used to classify multiple BRCA2 variants. Specifically, we utilized the BRCA2-independent single-stranded annealing pathway to efficiently knock-in targeted mutations in mESCs using a library of single-stranded DNA donors. Not only do our results support the ClinVar classification of some variants, but we also assign clinical significance to several VUSs, for which no functional data is currently available. [ABSTRACT FROM AUTHOR]

The spread of extended-spectrum beta-lactamases (ESBLs) in nosocomial and community-acquired enterobacteria is an important challenge for clinicians due to the limited therapeutic options for infections that are caused by these organisms. Here, we developed a panel of ESBL coding genes, evaluated the abundance and prevalence of ESBL encoding genes in patients undergoing H. pylori eradication therapy, and summarized the effects of eradication therapy on functional profiles of the gut microbiome. To assess the repertoire of known beta lactamase (BL) genes, they were divided into clusters according to their evolutionary relation. Primers were designed for amplification of cluster marker regions, and the efficiency of this amplification panel was assessed in 120 fecal samples acquired from 60 patients undergoing H. pylori eradication therapy. In addition, fecal samples from an additional 30 patients were used to validate the detection efficiency of the developed ESBL panel. The presence for majority of targeted clusters was confirmed by NGS of amplification products. Metagenomic sequencing revealed that the abundance of ESBL genes within the pool of microorganisms was very low. The global relative abundances of the ESBL-coding gene clusters did not differ significantly among treatment states. However, at the level of each cluster, classical ESBL producers such as Klebsiella sp. for blaOXY (p = 0.0076), Acinetobacter sp. for blaADC (p = 0.02297) and others, differed significantly with a tendency to decrease compared to the pre- and post-eradication states. Only 13 clusters were common across all three datasets, suggesting a patient-specific distribution profile of ESBL-coding genes. The number of AMR genes detected in the post-eradication state was higher than that in the pre-eradication state, which could be attributed, at least in part, to the therapy. This study demonstrated that the ESBL screening panel was effective in targeting ESBL-coding gene clusters from bacterial DNA and that minor differences exist in the abundance and prevalence of ESBL-coding gene levels before and after eradication therapy. [ABSTRACT FROM AUTHOR]

The majority of rare diseases are genetic, and regardless of advanced high-throughput genomics-based investigations, 60% of patients remain undiagnosed. A major factor limiting our ability to identify disease-causing alterations is a poor understanding of the morbid and normal human genome. A major genomic contributor of which function and distribution remain largely unstudied are the transposable elements (TE), which constitute 50% of our genome. Here we aim to resolve this knowledge gap and increase the diagnostic yield of rare disease patients investigated with clinical genome sequencing. To this end we characterized TE insertions in 1000 Swedish individuals from the SweGen dataset and 2504 individuals from the 1000 Genomes Project (1KGP), creating seven population-specific TE insertion databases. Of note, 66% of TE insertions in SweGen were present at >1% in the 1KGP databases, proving that most insertions are common across populations. Focusing on the rare TE insertions, we show that even though ~0.7% of those insertions affect protein coding genes, they rarely affect known disease casing genes (<0.1%). Finally, we applied a TE insertion identification workflow on two clinical cases where disease causing TE insertions were suspected and could verify the presence of pathogenic TE insertions in both. Altogether we demonstrate the importance of TE insertion detection and highlight possible clinical implications in rare disease diagnostics. [ABSTRACT FROM AUTHOR]

Astragalus is the largest flowering plant genus. We assembled the plastid genomes of four Astragalus species (Astragalus iranicus, A. macropelmatus, A. mesoleios, A. odoratus) using next-generation sequencing and analyzed their plastomes including genome organization, codon usage, nucleotide diversity, prediction of RNA editing and etc. The total length of the newly sequenced Astragalus plastomes ranged from 121,050 bp to 123,622 bp, with 110 genes comprising 76 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes. Comparative analysis of the chloroplast genomes of Astragalus revealed several hypervariable regions comprising three non-coding sites (trnQ(UUG)–accD, rps7 –trnV(GAC) and trnR(ACG)–trnN(GUU)) and four protein-coding genes (ycf1, ycf2, accD and clpP), which have potential as molecular markers. Positive selection signatures were found in five genes in Astragalus species including rps11, rps15, accD, clpP and ycf1. The newly sequenced species, A. macropelmatus, has an approximately 13-kb inversion in IR region. Phylogenetic analysis based on 75 protein-coding gene sequences confirmed that Astragalus form a monophyletic clade within the tribe Galegeae and Oxytropis is sister group to the Coluteoid clade. The results of this study may helpful in elucidating the chloroplast genome structure, understanding the evolutionary dynamics at genus Astragalus and IRLC levels and investigating the phylogenetic relationships. Moreover, the newly plastid genomes sequenced have been increased the plastome data resources on Astragalus that can be useful in further phylogenomic studies. [ABSTRACT FROM AUTHOR]

One of the main characteristics of prokaryotic genomes is the ratio in which guanine-cytosine bases are used in their DNA sequences. This is known as the genomic GC content and varies widely, from values below 20% to values greater than 74%. It has been demonstrated that the genomic GC content varies in accordance with the phylogenetic distribution of organisms and influences the amino acid composition of their corresponding proteomes. This bias is particularly important for amino acids that are coded by GC content-rich codons such as alanine, glycine, and proline, as well as amino acids that are coded by AT-rich codons, such as lysine, asparagine, and isoleucine. In our study, we extend these results by considering the effect of the genomic GC content on the secondary structure of proteins. On a set of 192 representative prokaryotic genomes and proteome sequences, we identified through a bioinformatic study that the composition of the secondary structures of the proteomes varies in relation to the genomic GC content; random coils increase as the genomic GC content increases, while alpha-helices and beta-sheets present an inverse relationship. In addition, we found that the tendency of an amino acid to form part of a secondary structure of proteins is not ubiquitous, as previously expected, but varies according to the genomic GC content. Finally, we discovered that for some specific groups of orthologous proteins, the GC content of genes biases the composition of secondary structures of the proteins for which they code. [ABSTRACT FROM AUTHOR]

Adenosine-to-inosine RNA editing is essential to prevent undesired immune activation. This diverse process alters the genetic content of the RNA and may recode proteins, change splice sites and miRNA targets, and mimic genomic mutations. Recent studies have associated or implicated aberrant editing with pathological conditions, including cancer, autoimmune diseases, and neurological and psychiatric conditions. RNA editing patterns in cardiovascular tissues have not been investigated systematically so far, and little is known about its potential role in cardiac diseases. Some hints suggest robust editing in this system, including the fact that ADARB1 (ADAR2), the main coding-sequence editor, is most highly expressed in these tissues. Here we characterized RNA editing in the heart and arteries and examined a contributory role to the development of atherosclerosis and two structural heart diseases -Ischemic and Dilated Cardiomyopathies. Analyzing hundreds of RNA-seq samples taken from the heart and arteries of cardiac patients and controls, we find that global editing, alongside inflammatory gene expression, is increased in patients with atherosclerosis, cardiomyopathies, and heart failure. We describe a single recoding editing site and suggest it as a target for focused research. This recoding editing site in the IGFBP7 gene is one of the only evolutionary conserved sites between mammals, and we found it exhibits consistently increased levels of editing in these patients. Our findings reveal that RNA editing is abundant in arteries and is elevated in several key cardiovascular conditions. They thus provide a roadmap for basic and translational research of RNA as a mediator of atherosclerosis and non-genetic cardiomyopathies. Author summary: The human genetic code is highly preserved, yet RNA editing is a process that changes it in RNA molecules. This may alter the resultant proteins or the properties of the RNA strands themselves, leading to changes in their special structure and affecting their interaction with the immune system. Unedited RNA may be highly immunogenic, and studies of recent years demonstrated that dysregulated editing causes an inflammatory response. We show that human RNA editing is especially active in the arteries and investigate its role in cardiovascular diseases such as atherosclerosis and cardiomyopathies. These diseases have an underlying inflammatory component that eventually leads to the loss of normal structure and function. We show that RNA editing is increased in patients and draw a line connecting it with a concomitant increase in inflammatory pathways. We also detected specific coding editing sites, mainly in the IGFBP7 gene, where aberrant editing is present, resulting in an altered protein product. We delineate these changes and suggest they may contribute to or are a marker of the underlying pathology. Our findings indicate that RNA editing is a key player in several cardiovascular diseases. [ABSTRACT FROM AUTHOR]

Transcription factors read the genome, fundamentally connecting DNA sequence to gene expression across diverse cell types. Determining how, where, and when TFs bind chromatin will advance our understanding of gene regulatory networks and cellular behavior. The 2017 ENCODE-DREAM in vivo Transcription-Factor Binding Site (TFBS) Prediction Challenge highlighted the value of chromatin accessibility data to TFBS prediction, establishing state-of-the-art methods for TFBS prediction from DNase-seq. However, the more recent Assay-for-Transposase-Accessible-Chromatin (ATAC)-seq has surpassed DNase-seq as the most widely-used chromatin accessibility profiling method. Furthermore, ATAC-seq is the only such technique available at single-cell resolution from standard commercial platforms. While ATAC-seq datasets grow exponentially, suboptimal motif scanning is unfortunately the most common method for TFBS prediction from ATAC-seq. To enable community access to state-of-the-art TFBS prediction from ATAC-seq, we (1) curated an extensive benchmark dataset (127 TFs) for ATAC-seq model training and (2) built "maxATAC", a suite of user-friendly, deep neural network models for genome-wide TFBS prediction from ATAC-seq in any cell type. With models available for 127 human TFs, maxATAC is the largest collection of high-performance TFBS prediction models for ATAC-seq. maxATAC performance extends to primary cells and single-cell ATAC-seq, enabling improved TFBS prediction in vivo. We demonstrate maxATAC's capabilities by identifying TFBS associated with allele-dependent chromatin accessibility at atopic dermatitis genetic risk loci. Author summary: Proteins called transcription factors interpret the genome, reading both DNA sequence and chromatin state, to orchestrate gene expression across the diversity of human cell types. In any given cell type, most chromatin is "inaccessible", and only those parts of the genetic code needed or likely to be needed soon are "accessible" for transcription factor binding to affect gene expression and cellular behavior. Hundreds of transcription factors are expressed in a given cell type and context (e.g., age, disease), and knowledge of their context-specific DNA binding sites is key to uncovering how transcription factors regulate cellular behaviors in health or disease. However, experimentally profiling the >1,600 human transcription factors across all cell types and contexts is infeasible. We built a suite of computational models "maxATAC" to predict transcription factor binding from a measurement of accessible chromatin, ATAC-seq. Importantly, ATAC-seq is feasible even at single-cell resolution. Thus, this data type, in combination with maxATAC, can be used to infer transcription factor binding sites in directly-relevant cell types isolated from physiological and disease settings, enabling insights into disease mechanisms, including how genetic variants and cellular context impact transcription factor binding, gene expression patterns and disease risk. [ABSTRACT FROM AUTHOR]

As ribosomes translate the genetic code, they can encounter a variety of obstacles that hinder their progress. If ribosomes stall for prolonged times, cells suffer due to the loss of translating ribosomes and the accumulation of aberrant protein products. Thus to protect cells, stalled ribosomes experience a series of reactions to relieve the stall and degrade the offending mRNA, a process known as No-Go mRNA Decay (NGD). While much of the machinery for NGD is known, the precise ordering of events and factors along this pathway has not been tested. Here, we deploy C. elegans to unravel the coordinated events comprising NGD. Utilizing a novel reporter and forward and reverse genetics, we identify the machinery required for NGD. Our subsequent molecular analyses define a functional requirement for ubiquitination on at least two ribosomal proteins (eS10 and uS10), and we show that ribosomes lacking ubiquitination sites on eS10 and uS10 fail to perform NGD in vivo. We show that the nuclease NONU-1 acts after the ubiquitin ligase ZNF-598, and discover a novel requirement for the ribosome rescue factors HBS-1/PELO-1 in mRNA decay via NONU-1. Taken together, our work demonstrates mechanisms by which ribosomes signal to effectors of mRNA repression, and we delineate links between repressive factors working toward a well-defined NGD pathway. Author summary: Ribosomes are large molecular machines entrusted with the essential task of reading mRNAs and building proteins. It is critical that ribosomes proceed with their work undeterred, so as to avoid traffic jams between ribosomes on mRNAs; however, they sometimes experience challenges which make traffic jams unavoidable. In these cases, cells have evolved machinery to clean up stalled ribosomes and remove offending mRNAs. Much of this machinery is now identified, but the relationships between factors has yet to be tested. Here, in C. elegans we identify the machinery that responds to stalls, and we order factors relative to one another. We find that ribosomes must be tagged by ubiquitin to signal mRNA decay, as ribosomes lacking ubiquitin sites fail to elicit mRNA decay. Furthermore, we find that these ubiquitin signals recruit an enzyme which cuts mRNA, and we present evidence that ribosome rescue factors are required for the mRNA decay reaction as well. Overall, our work provides new connections between factors orchestrating the elimination of problematic mRNAs that stall ribosomes. [ABSTRACT FROM AUTHOR]

The fast rate of viral mutations of SARS CoV-2 result in decrease in the efficacy of the vaccines that have been developed before the emergence of these mutations. Thus, it is believed that using additional measures to combat the virus is not only advisable but also beneficial. Two antiviral drugs were authorized for emergency use by the FDA, namely Pfizer's two-drug regimen sold under the brand name Paxlovid, and Merck's drug Lagevrio. Pfizer's two-drug combination consists of nirmatrelvir, a protease inhibitor that blocks coronavirus ability to multiply and another antiviral, ritonavir, that lowers the rate of drug clearance to boost the longevity and activity of the protease inhibitor. Merck's drug Lagevrio (molnupiravir) is a nucleoside analogue with a mechanism of action that aims to introduce errors into the genetic code of the virus. We believe the armament against the virus can be augmented by the addition of another class of enzyme inhibitors that are required for viral survival and its ability to replicate. Enzymes like nsp14 and nsp10/16 methyltransferases (MTases) represent another class of drug targets since they are required for viral RNA translation and evading the host immune system. In this communication, we have successfully verified that the MTase-Glo, which is universal and homogeneous MTase assay can be used to screen for inhibitors of the two pivotal enzymes nsp14 and nsp16 of SARS CoV-2. Furthermore, we have carried out extensive studies on those enzymes using different RNA substrates and tested their activity using various inhibitors and verified the utility of this assay for use in drug screening programs. We anticipate our work will be pursued further to screen for large libraries to discover new and selective inhibitors for the viral enzymes particularly that these enzymes are structurally different from their mammalian counterparts. [ABSTRACT FROM AUTHOR]

A field isolate (Reo/SDWF /Pheasant/17608/20) of avian orthoreovirus (ARV), isolated from a flock of game-pheasants in Weifang, Shandong Province, was genetically characterized being a field variant or novel strain in our recent research studies in conducting whole genome sequencing by using Next-Generation Sequencing (NGS) technique on Illumina MiSeq platform. Among a total of 870,197 35-151-mer sequencing reads, 297,711 reads (34.21%) were identified as ARV sequences. The de novo assembly of the ARV reads resulted in generation of 10 ARV-related contigs with the average sequencing coverage from 1390× to 1977× according to 10 ARV genome segments. The complete genomes of this pheasant-origin ARV (Reo/SDWF /Pheasant/17608/20) were 23,495 bp in length and consist of 10 dsRNA segments ranged from 1192 bp (S4) to 3958 bp (L1) encoding 12 viral proteins. Sequence comparison between the SDWF17608 and classic ARV reference strains revealed that 58.1–100% nucleotide (nt) identities and 51.4–100% amino acid (aa) identities were in genome segment coding genes. The 10 RNA segments had conversed termini at 5' (5'-GCUUUU) and 3' (UCAUC-3') side, which were identical to the most published ARV strains. Phylogenetic analysis revealed that this pheasant ARV field variant was closely related with chicken ARV strains in 7 genome segment genes, but it possessed significant sequence divergence in M1, M3 and S2 segments. These findings suggested that this pheasant-origin field variant was a divergent ARV strain and was likely originated from reassortments between different chicken ARV strains. [ABSTRACT FROM AUTHOR]

The luciferin sulfokinase (coelenterazine sulfotransferase) of Renilla was previously reported to activate the storage form, luciferyl sulfate (coelenterazine sulfate) to luciferin (coelenterazine), the substrate for the luciferase bioluminescence reaction. The gene coding for the coelenterazine sulfotransferase has not been identified. Here we used a combined proteomic/transcriptomic approach to identify and clone the sulfotransferase cDNA. Multiple isoforms of coelenterazine sulfotransferase were identified from the anthozoan Renilla muelleri by intersecting its transcriptome with the LC-MS/MS derived peptide sequences of coelenterazine sulfotransferase purified from Renilla. Two of the isoforms were expressed in E. coli, purified, and partially characterized. The encoded enzymes display sulfotransferase activity that is comparable to that of the native sulfotransferase isolated from Renilla reniformis that was reported in 1970. The bioluminescent assay for sensitive detection of 3'-phosphoadenosine 5'-phosphate (PAP) using the recombinant sulfotransferase is demonstrated. [ABSTRACT FROM AUTHOR]

Multi-marker metabarcoding is increasingly being used to generate biodiversity information across different domains of life from microbes to fungi to animals such as for molecular ecology and biomonitoring applications in different sectors from academic research to regulatory agencies and industry. Current popular bioinformatic pipelines support microbial and fungal marker analysis, while ad hoc methods are often used to process animal metabarcode markers from the same study. MetaWorks provides a harmonized processing environment, pipeline, and taxonomic assignment approach for demultiplexed Illumina reads for all biota using a wide range of metabarcoding markers such as 16S, ITS, and COI. A Conda environment is provided to quickly gather most of the programs and dependencies for the pipeline. Several workflows are provided such as: taxonomically assigning exact sequence variants, provides an option to generate operational taxonomic units, and facilitates single-read processing. Pipelines are automated using Snakemake to minimize user intervention and facilitate scalability. All pipelines use the RDP classifier to provide taxonomic assignments with confidence measures. We extend the functionality of the RDP classifier for taxonomically assigning 16S (bacteria), ITS (fungi), and 28S (fungi), to also support COI (eukaryotes), rbcL (eukaryotes, land plants, diatoms), 12S (fish, vertebrates), 18S (eukaryotes, diatoms) and ITS (fungi, plants). MetaWorks properly handles ITS by trimming flanking conserved rRNA gene regions as well as protein coding genes by providing two options for removing obvious pseudogenes. MetaWorks can be downloaded from https://github.com/terrimporter/MetaWorks and quickstart instructions, pipeline details, and a tutorial for new users can be found at https://terrimporter.github.io/MetaWorksSite. [ABSTRACT FROM AUTHOR]

Messenger RNA (mRNA) translation can lead to higher rates of mRNA decay, suggesting the ribosome plays a role in mRNA destruction. Furthermore, mRNA features, such as codon identities, which are directly probed by the ribosome, correlate with mRNA decay rates. Many amino acids are encoded by synonymous codons, some of which are decoded by more abundant tRNAs leading to more optimal translation and increased mRNA stability. Variable translation rates for synonymous codons can lead to ribosomal collisions as ribosomes transit regions with suboptimal codons, and ribosomal collisions can promote mRNA decay. In addition to different translation rates, the presence of certain codons can also lead to higher or lower rates of amino acid misincorporation which could potentially lead to protein misfolding if a substituted amino acid fails to make critical contacts in a structure. Here, we test whether Geneticin—G418, an aminoglycoside antibiotic known to promote amino acid misincorporation—affects mRNA stability. We observe that G418 decreases firefly luciferase mRNA stability in an in vitro translation system and also reduces mRNA stability in mouse embryonic stem cells (mESCs). G418-sensitive mRNAs are enriched for certain optimal codons that contain G or C in the wobble position, arguing that G418 blunts the stabilizing effects of codon optimality. [ABSTRACT FROM AUTHOR]

Patterns of non-uniform usage of synonymous codons vary across genes in an organism and between species across all domains of life. This codon usage bias (CUB) is due to a combination of non-adaptive (e.g. mutation biases) and adaptive (e.g. natural selection for translation efficiency/accuracy) evolutionary forces. Most models quantify the effects of mutation bias and selection on CUB assuming uniform mutational and other non-adaptive forces across the genome. However, non-adaptive nucleotide biases can vary within a genome due to processes such as biased gene conversion (BGC), potentially obfuscating signals of selection on codon usage. Moreover, genome-wide estimates of non-adaptive nucleotide biases are lacking for non-model organisms. We combine an unsupervised learning method with a population genetics model of synonymous coding sequence evolution to assess the impact of intragenomic variation in non-adaptive nucleotide bias on quantification of natural selection on synonymous codon usage across 49 Saccharomycotina yeasts. We find that in the absence of a priori information, unsupervised learning can be used to identify genes evolving under different non-adaptive nucleotide biases. We find that the impact of intragenomic variation in non-adaptive nucleotide bias varies widely, even among closely-related species. We show that the overall strength and direction of translational selection can be underestimated by failing to account for intragenomic variation in non-adaptive nucleotide biases. Interestingly, genes falling into clusters identified by machine learning are also physically clustered across chromosomes. Our results indicate the need for more nuanced models of sequence evolution that systematically incorporate the effects of variable non-adaptive nucleotide biases on codon frequencies. Author summary: Codon usage bias (CUB), or the unequal usage of codons of the same amino acid (i.e. synonymous codons), has been observed in species across all domains of life. CUB is known to be shaped by both non-adaptive (e.g. mutation biases) and adaptive (e.g. natural selection for translation efficiency/accuracy) evolution. A key challenge for researchers is disentangling the role of various processes shaping codon usage, often for the purpose of identifying codons favored by natural selection, sometimes referred to as "optimal" or "preferred" codons. Despite large variation in non-adaptive nucleotide biases within a genome, most methods to quantify natural selection typically ignore this variation for the sake of simplicity. Here, we combine a population genetics model with unsupervised machine learning to identify genes evolving under different non-adaptive nucleotide biases across 49 budding yeasts species. We find that ignoring for variation in non-adaptive nucleotide biases can obfuscate signals of selection on codon usage. Our results indicate the need for more nuanced models of coding sequence evolution. [ABSTRACT FROM AUTHOR]

Central nervous system-expressed long non-coding RNAs (lncRNAs) are often located in the genome close to protein coding genes involved in transcriptional control. Such lncRNA-protein coding gene pairs are frequently temporally and spatially co-expressed in the nervous system and are predicted to act together to regulate neuronal development and function. Although some of these lncRNAs also bind and modulate the activity of the encoded transcription factors, the regulatory mechanisms controlling co-expression of neighbouring lncRNA-protein coding genes remain unclear. Here, we used high resolution NG Capture-C to map the cis-regulatory interaction landscape of the key neuro-developmental Paupar-Pax6 lncRNA-mRNA locus. The results define chromatin architecture changes associated with high Paupar-Pax6 expression in neurons and identify both promoter selective as well as shared cis-regulatory-promoter interactions involved in regulating Paupar-Pax6 co-expression. We discovered that the TCF7L2 transcription factor, a regulator of chromatin architecture and major effector of the Wnt signalling pathway, binds to a subset of these candidate cis-regulatory elements to coordinate Paupar and Pax6 co-expression. We describe distinct roles for Paupar in Pax6 expression control and show that the Paupar DNA locus contains a TCF7L2 bound transcriptional silencer whilst the Paupar transcript can act as an activator of Pax6. Our work provides important insights into the chromatin interactions, signalling pathways and transcription factors controlling co-expression of adjacent lncRNAs and protein coding genes in the brain. Author summary: Long non-coding RNA (lncRNA) genes are often co-expressed in the brain with neighbouring protein coding genes involved in gene expression control. Such lncRNA-protein coding gene pairs are predicted to work together to regulate neuronal development and function. Despite this, the regulatory mechanisms controlling their co-expression is poorly understood. Here, we identify the chromatin interactions and DNA sequences controlling expression of the Paupar-Pax6 lncRNA-protein coding locus in the brain. We show that the Wnt signalling pathway, a key regulator of development and disease, acts through the TCF7L2 transcription factor to co-ordinate their expression. We find that Paupar DNA contains sequences that silence Pax6 whilst the Paupar RNA can act as an activator of Pax6. Our work generates new insights into the complex regulatory relationship controlling the function of neighbouring lncRNA-protein coding genes and is important for understanding development of the brain. [ABSTRACT FROM AUTHOR]

Prunus persica is one of the main stone fruit crops in Crimea and southern Russia. The P. persica genome has recently been sequenced and annotated in good quality. However, for a deeper assessment of the peach genome, it is necessary to include in the research other cultivars that are in the collection of the Nikitsky Botanical Garden. The cultivars of the Nikitsky Botanical Garden are unique and differ from Western European and American ones, as they are derived from cultivars and forms originating from Central Asian, North Caucasian, Transcaucasian and Eastern European countries. In this paper, we present the assembly of the P. persica cv. 'Sovetskiy' genome obtained using Oxford Nanopore long reads and Illumina short reads by hybrid assembly methods. The assembled genome of P. persica cv. 'Sovetskiy' is 206.26 MB in 226 scaffolds, with N50 24 Mb, including 8 chromosomes. It contains 27140 coding genes, 26973 (99.38%) of which are annotated in at least one functional database. More than 36.05% of the genome regions were identified as repeating elements. [ABSTRACT FROM AUTHOR]

Avian influenza virus (AIV) represents a major concern with productive implications in poultry systems but it is also a zoonotic agent that possesses an intrinsic pandemic risk. AIV is an enveloped, negative-sense and single-stranded RNA virus with a segmented genome. The eight genomic segments, comprising the whole genome, encode for eleven proteins. Within these proteins, Hemagglutinin (HA) and Neuraminidase (NA) are the most relevant for studies of evolution and pathogenesis considering their role in viral replication, and have also been used for classification purposes. Migratory birds are the main hosts and play a pivotal role in viral evolution and dissemination due to their migratory routes that comprise large regions worldwide. Altogether, viral and reservoir factors contribute to the emergence of avian influenza viruses with novel features and pathogenic potentials. The study aimed to conduct surveillance of AIVs in wild birds from Peru. A multi-site screening of feces of migratory birds was performed to isolate viruses and to characterize the whole genome sequences, especially the genes coding for HA and NA proteins. Four-hundred-twenty-one (421) fecal samples, collected between March 2019 and March 2020 in Lima, were obtained from 21 species of wild birds. From these, we isolated five AIV from whimbrel, kelp gull, Franklin's gulls and Mallard, which were of low pathogenicity, including four subtypes as H6N8, H13N6, H6N2 and H2N6. Genetic analysis of HA and NA genes revealed novel features in these viruses and phylogenetic analysis exhibited a close relationship with those identified in North America (US and Canada). Furthermore, H2N6 isolate presented a NA sequence with higher genetic relationship to Chilean isolates. These results highlight that the geographical factor is of major relevance in the evolution of AIV, suggesting that AIV circulating in Peru might represent a new site for the emergence of reassortant AIVs. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) are small non-coding RNAs that regulate protein-coding gene expression primarily found in plants and animals. Fungi produce microRNA-like RNAs (milRNAs) that are structurally similar to miRNAs and functionally important in various biological processes. The fungus Fusarium oxysporum f. sp. cubense (Foc) is the causal agent of Banana Fusarium vascular wilt that threatens global banana production. It remains uncharacterized about the biosynthesis and functions of milRNAs in Foc. In this study, we investigated the biological function of milRNAs contributing to Foc pathogenesis. Within 24 hours post infecting the host, the Argonaute coding gene FoQDE2, and two Dicer coding genes FoDCL1 and FoDCL2, all of which are involved in milRNA biosynthesis, were significantly induced. FoQDE2 deletion mutant exhibited decreased virulence, suggesting the involvement of milRNA biosynthesis in the Foc pathogenesis. By small RNA sequencing, we identified 364 small RNA-producing loci in the Foc genome, 25 of which were significantly down-regulated in the FoQDE2 deletion mutant, from which milR-87 was verified as a FoQDE2-depedent milRNA based on qRT-PCR and Northern blot analysis. Compared to the wild-type, the deletion mutant of milR-87 was significantly reduced in virulence, while overexpression of milR-87 enhanced disease severity, confirming that milR-87 is crucial for Foc virulence in the infection process. We furthermore identified FOIG_15013 (a glycosyl hydrolase-coding gene) as the direct target of milR-87 based on the expression of FOIG_15013-GFP fusion protein. The FOIG_15013 deletion mutant displayed similar phenotypes as the overexpression of milR-87, with a dramatic increase in the growth, conidiation and virulence. Transient expression of FOIG_15013 in Nicotiana benthamiana leaves activates the host defense responses. Collectively, this study documents the involvement of milRNAs in the manifestation of the devastating fungal disease in banana, and demonstrates the importance of milRNAs in the pathogenesis and other biological processes. Further analyses of the biosynthesis and expression regulation of fungal milRNAs may offer a novel strategy to combat devastating fungal diseases. Author summary: The fungus Fusarium oxysporum f. sp. cubense (Foc) is the causal agent of Banana Fusarium vascular wilt that threatens global banana production. However, knowledge about pathogenesis of Foc is limited. In particular, pathogenic regulatory mechanism of the microRNA like small RNAs (milRNAs) found in Foc is unknown. Here, we found that FoQDE2, an Argonaute coding gene, and two Dicer coding genes FoDCL1 and FoDCL2, which are involved in milRNA biosynthesis, are significantly induced during the early infection stage of Foc. The results suggested that the milRNAs biosynthesis mediated by these genes may play an active role in the infection process of Foc. Based on this assumption, we subsequently found a FoQDE2-dependent milRNA (milR-87) and identified its target gene. Functional analysis showed that FoQDE2, milR-87 and its target gene were involved in the pathogenicity of Foc in different degree. The studies help us gain insight into the pathogenesis with FoQDE2, milR-87, and its target gene as central axis in Foc. The identified pathogenicity-involved milRNA provides an active target for developing novel and efficient biocontrol agents against Banana Fusarium wilt. [ABSTRACT FROM AUTHOR]

The alpha subunit of the cytoplasmic Phenylalanyl tRNA synthetase (α-PheRS, FARSA in humans) displays cell growth and proliferation activities and its elevated levels can induce cell fate changes and tumor-like phenotypes that are neither dependent on the canonical function of charging tRNAPhe with phenylalanine nor on stimulating general translation. In intestinal stem cells of Drosophila midguts, α-PheRS levels are naturally slightly elevated and human FARSA mRNA levels are elevated in multiple cancers. In the Drosophila midgut model, elevated α-PheRS levels caused the accumulation of many additional proliferating cells resembling intestinal stem cells (ISCs) and enteroblasts (EBs). This phenotype partially resembles the tumor-like phenotype described as Notch RNAi phenotype for the same cells. Genetic interactions between α-PheRS and Notch suggest that their activities neutralize each other and that elevated α-PheRS levels attenuate Notch signaling when Notch induces differentiation into enterocytes, type II neuroblast stem cell proliferation, or transcription of a Notch reporter. These non-canonical functions all map to the N-terminal part of α-PheRS which accumulates naturally in the intestine. This truncated version of α-PheRS (α-S) also localizes to nuclei and displays weak sequence similarity to the Notch intracellular domain (NICD), suggesting that α-S might compete with the NICD for binding to a common target. Supporting this hypothesis, the tryptophan (W) residue reported to be key for the interaction between the NICD and the Su(H) BTD domain is not only conserved in α-PheRS and α-S, but also essential for attenuating Notch signaling. Author summary: Aminoacyl tRNA synthetases charge tRNAs with their cognate amino acid to ensure proper decoding of the genetic code during translation. Independent of its aminoacylation function, the alpha subunit of Drosophila cytoplasmic Phenylalanyl tRNA synthetase (α-PheRS, FARSA in humans) has an additional activity that promotes growth and proliferation. Here we describe that elevated α-PheRS levels also induce cell fate changes and tumorous phenotypes in Drosophila midguts. Excessive proliferating cells with stem and progenitor cell characteristics accumulate and the composition of the terminally differentiated cells changes, too. This phenotype together with observed genetic interactions between α-PheRS and Notch levels show that α-PheRS counteracts Notch signaling in many different tissues and developmental stages. This novel activity of α-PheRS maps to its N-terminal part, which is naturally produced. The fragment contains a DNA binding domain, translocates into nuclei, and displays essential similarities to a Notch domain that binds to the downstream transcription factor. This suggests that it might be competing with Notch for binding to a common target. Not only because Notch plays important roles in many tumors, but also because FARSA mRNA levels are considerably upregulated in many tumors, this novel activity deserves more attention for cancer research. [ABSTRACT FROM AUTHOR]

Pack-TYPE transposable elements (TEs) are a group of non-autonomous DNA transposons found in plants. These elements can efficiently capture and shuffle coding DNA across the host genome, accelerating the evolution of genes. Despite their relevance for plant genome plasticity, the detection and study of Pack-TYPE TEs are challenging due to the high similarity these elements have with genes. Here, we produced an automated annotation pipeline designed to study Pack-TYPE elements and used it to successfully annotate and analyse more than 10,000 new Pack-TYPE TEs in the rice and maize genomes. Our analysis indicates that Pack-TYPE TEs are an abundant and heterogeneous group of elements. We found that these elements are associated with all main superfamilies of Class II DNA transposons in plants and likely share a similar mechanism to capture new chromosomal DNA sequences. Furthermore, we report examples of the direct contribution of these TEs to coding genes, suggesting a generalised and extensive role of Pack-TYPE TEs in plant genome evolution. Author summary: Transposable Elements (TEs) are genetic DNA sequences able to move across the genome, and their transposition activity is associated with genome plasticity and gene evolution. However, most of these elements exhibit "selfish" behaviour, meaning that they mainly transpose their own DNA sequence and only exceptionally might rearrange the DNA of coding genes. Pack-TYPE TEs, found in plants, represent an important exception, and they can efficiently capture and shuffle DNA sequences captured from the genome, accelerating the evolution of genes. We provide here the first automatic pipeline designed explicitly for the annotation of Pack-TYPE TEs. We used our approach to systematically investigate Pack-TYPE TEs in the rice and maize reference genomes, and annotated thousands of new elements in these species. We demonstrate that Pack-TYPE elements are abundant in plants and we report several examples of coding genes originated as a consequence of the mobilization of these elements. [ABSTRACT FROM AUTHOR]

Gene overlap occurs when two or more genes are encoded by the same nucleotides. This phenomenon is found in all taxonomic domains, but is particularly common in viruses, where it may increase the information content of compact genomes or influence the creation of new genes. Here we report a global comparative study of overlapping open reading frames (OvRFs) of 12,609 virus reference genomes in the NCBI database. We retrieved metadata associated with all annotated open reading frames (ORFs) in each genome record to calculate the number, length, and frameshift of OvRFs. Our results show that while the number of OvRFs increases with genome length, they tend to be shorter in longer genomes. The majority of overlaps involve +2 frameshifts, predominantly found in dsDNA viruses. Antisense overlaps in which one of the ORFs was encoded in the same frame on the opposite strand (−0) tend to be longer. Next, we develop a new graph-based representation of the distribution of overlaps among the ORFs of genomes in a given virus family. In the absence of an unambiguous partition of ORFs by homology at this taxonomic level, we used an alignment-free k-mer based approach to cluster protein coding sequences by similarity. We connect these clusters with two types of directed edges to indicate (1) that constituent ORFs are adjacent in one or more genomes, and (2) that these ORFs overlap. These adjacency graphs not only provide a natural visualization scheme, but also a novel statistical framework for analyzing the effects of gene- and genome-level attributes on the frequencies of overlaps. Author summary: Gene overlap occurs when the same part of a genome encodes two or more genes. This phenomenon is found in all biological domains of life, but it is particularly common in viruses, where it may play a role in making viral genomes more compact. To understand the prevalence of overlapping genes in viruses, we analyzed over 12,000 genomes of every known type of virus for which this genetic information is available. Although overlaps are more abundant in viruses with larger genomes, for instance, they are also significantly shorter. Overlaps in which one of the genes is read in the opposite direction (−0 overlaps) tend to be longer, which may be an emergent property of the universal genetic code. We developed a new computational method to analyze and visualize the distribution of overlaps among genomes belonging to a group (family) of viruses as a network. This approach enabled us to identify distinct patterns in the organization of genomes within virus families; for example, gene overlap in the coronavirus family tends to involve non-essential genes outside of the "core" of the network of genes. [ABSTRACT FROM AUTHOR]

The introduction of frameshifting non-3n indels enables the identification of gene-trait associations. However, it has been hypothesised that recovery of the original reading frame owing to usage of non-canonical splice forms could cause rescue. To date there is very little evidence for organism-level rescue by such a mechanism and it is unknown how commonly indels induce, or are otherwise associated with, frame-restoring splice forms. We perform CRISPR/Cas9 editing of randomly selected loci in rice to investigate these issues. We find that the majority of loci have a frame-restoring isoform. Importantly, three quarters of these isoforms are not seen in the absence of the indels, consistent with indels commonly inducing novel isoforms. This is supported by analysis in the context of NMD knockdowns. We consider in detail the two top rescue candidates, in wax deficient anther 1 (wda1) and brittle culm (bc10), finding that organismal-level rescue in both cases is strong but owing to different splice modification routes. More generally, however, as frame-restoring isoforms are low abundance and possibly too disruptive, such rescue we suggest to be the rare exception, not the rule. Nonetheless, assuming that indels commonly induce frame-restoring isoforms, these results emphasize the need to examine RNA level effects of non-3n indels and suggest that multiple non-3n indels in any given gene are advisable to probe a gene's trait associations. Author summary: As protein coding genes are read in units of three (codons), insertions or deletions (indels) that are not a multiple of three long (non 3n) are expected to be especially harmful. Whether they are is important both for interpreting the results of non-3n indel experiments to probe a gene's functional importance and for diagnostics. Particularly enigmatic are incidences where some non-3n changes in a gene compromise phenotypes while other seemingly comparable ones do not. One explanation for the latter is that a non-3n indel might be rescued via a frame-restoring splice form. Here we examine this hypothesis by inducing non-3n indels in many genes in rice and find that many non-3n indels are associated with a splice form that restores the reading frame. In the majority of these cases the indel appears to induce the potential rescuing splice form. We examine two top hit cases in detail and show functional rescue by splice modification. More generally, the frame-restoring forms are, however, low abundance and probably result in compromised proteins. We conclude then that splice mediated rescue is possible, but probably uncommon. Nonetheless it should not be overlooked in experimental design and interpretation. [ABSTRACT FROM AUTHOR]

Viruses are highly evolvable, but what traits endow this property? The high mutation rates of viruses certainly play a role, but factors that act above the genetic code, like protein thermostability, are also expected to contribute. We studied how the thermostability of a model virus, bacteriophage λ, affects its ability to evolve to use a new receptor, a key evolutionary transition that can cause host-range evolution. Using directed evolution and synthetic biology techniques we generated a library of host-recognition protein variants with altered stabilities and then tested their capacity to evolve to use a new receptor. Variants fell within three stability classes: stable, unstable, and catastrophically unstable. The most evolvable were the two unstable variants, whereas seven of eight stable variants were significantly less evolvable, and the two catastrophically unstable variants could not grow. The slowly evolving stable variants were delayed because they required an additional destabilizing mutation. These results are particularly noteworthy because they contradict a widely supported contention that thermostabilizing mutations enhance evolvability of proteins by increasing mutational robustness. Our work suggests that the relationship between thermostability and evolvability is more complex than previously thought, provides evidence for a new molecular model of host-range expansion evolution, and identifies instability as a potential predictor of viral host-range evolution. Author summary: Understanding how viruses evolve to infect new hosts is critical for predicting host shifts as well as tuning host-range in phage therapy applications. Yet a mechanistic understanding of the molecular steps required to shift hosts has not been achieved. For this study we examined the evolutionary potential of different strains of a model virus, bacteriophage λ, to gain the ability to use a new receptor, a key step in host shifts. We discovered that λ variants with destabilized host-recognition proteins were more likely to evolve the necessary mutations to use the new receptor than stabilized variants. However, destabilization was only beneficial to a certain point and variants with overly unstable proteins lost all function. These results led us to propose a new molecular model for receptor use evolution in λ; 1) destabilizing mutations evolve that provide protein structural flexibility that allows new protein conformations to form that are able to interact with the new receptors, and 2) mutations evolve that alter the binding surface chemical properties to assist interactions with the new receptor. Our work with a model virus-host system, points to the potential use of viral stability as a phenotypic indicator of the capacity for virus host-range evolution. [ABSTRACT FROM AUTHOR]

During cell devision, maintaining the epigenetic information encoded in histone modification patterns is crucial for survival and identity of cells. The faithful inheritance of the histone marks from the parental to the daughter strands is a puzzle, given that each strand gets only half of the parental nucleosomes. Mapping DNA replication and reconstruction of modifications to equivalent problems in communication of information, we ask how well enzymes can recover the parental modifications, if they were ideal computing machines. Studying a parameter regime where realistic enzymes can function, our analysis predicts that enzymes may implement a critical threshold filling algorithm which fills unmodified regions of length at most k. This algorithm, motivated from communication theory, is derived from the maximum à posteriori probability (MAP) decoding which identifies the most probable modification sequence based on available observations. Simulations using our method produce modification patterns similar to what has been observed in recent experiments. We also show that our results can be naturally extended to explain inheritance of spatially distinct antagonistic modifications. Author summary: Chromatin is essentially the DNA that is folded and packaged with the help of proteins. While the nucleotide sequence in the DNA codes genetic information, the packaging of the DNA into chromatin encodes extra layer of information. This epigenetic code regulates reading of the genetic code and provides identity to cells—whether the cell is a skin cell or a brain cell, for example. In this work we examine a long standing puzzle, that is, how the epigenetic code in the form of histone modification patterns may get inherited, when a cell divides. Using theoretical arguments, we present an algorithm that enzymes could be executing so that the epigenetic code can be inherited with minimal error, soon after DNA replication. [ABSTRACT FROM AUTHOR]

Pseudomonas aeruginosa is an important opportunistic pathogen capable of causing variety of infections in humans. The type III secretion system (T3SS) is a critical virulence determinant of P. aeruginosa in the host infections. Expression of the T3SS is regulated by ExsA, a master regulator that activates the expression of all known T3SS genes. Expression of the exsA gene is controlled at both transcriptional and posttranscriptional levels. Here, we screened a P. aeruginosa transposon (Tn5) insertional mutant library and found rplI, a gene coding for the ribosomal large subunit protein L9, to be a repressor for the T3SS gene expression. Combining real-time quantitative PCR (qPCR), western blotting and lacZ fusion assays, we show that RplI controls the expression of exsA at the posttranscriptional level. Further genetic experiments demonstrated that RplI mediated control of the exsA translation involves 5' untranslated region (5' UTR). A ribosome immunoprecipitation assay and qPCR revealed higher amounts of a 24 nt fragment from exsA mRNA being associated with ribosomes in the ΔrplI mutant. An interaction between RplI and exsA mRNA harboring its 24 nt, but not 12 nt, 5' UTR was confirmed by RNA Gel Mobility Shift and Microscale Thermophoresis assays. Overall, this study identifies the ribosomal large subunit protein L9 as a novel T3SS repressor that inhibits ExsA translation in P. aeruginosa. Author summary: Ribosomes provide all living organisms the capacity to synthesize proteins. The production of many ribosomal proteins is often controlled by an autoregulatory feedback mechanism. P. aeruginosa is an opportunistic human pathogen and its type III secretion system (T3SS) is a critical virulence determinant in host infections. In this study, by screening a Tn5 mutant library, we identified rplI, encoding ribosomal large subunit protein L9, as a novel repressor for the T3SS. Further exploring the regulatory mechanism, we found that the RplI protein interacts with the 5' UTR (5' untranslated region) of exsA, a gene coding for transcriptional activator of the T3SS. Such an interaction likely blocks ribosome loading on the exsA 5' UTR, inhibiting the initiation of exsA translation. The significance of this work is in the identification of a novel repressor for the T3SS and elucidation of its molecular mechanism. Furthermore, this work provides evidence for individual ribosomal protein regulating mRNA translation beyond its autogenous feedback control. [ABSTRACT FROM AUTHOR]

Studies of drug resistance in the protozoan parasites of the genus Leishmania have been helpful in revealing biochemical pathways as potential drug targets. The chlorinated glutamine analogue acivicin has shown good activity against Leishmania cells and was shown to target several enzymes containing amidotransferase domains. We selected a Leishmania tarentolae clone for acivicin resistance. The genome of this resistant strain was sequenced and the gene coding for the amidotransferase domain-containing GMP synthase was found to be amplified. Episomal expression of this gene in wild-type L. tarentolae revealed a modest role in acivicin resistance. The most prominent defect observed in the resistant mutant was reduced uptake of glutamate, and through competition experiments we determined that glutamate and acivicin, but not glutamine, share the same transporter. Several amino acid transporters (AATs) were either deleted or mutated in the resistant cells. Some contributed to the acivicin resistance phenotype although none corresponded to the main glutamate transporter. Through sequence analysis one AAT on chromosome 22 corresponded to the main glutamate transporter. Episomal expression of the gene coding for this transporter in the resistant mutant restored glutamate transport and acivicin susceptibility. Its genetic knockout led to reduced glutamate transport and acivicin resistance. We propose that acivicin binds covalently to this transporter and as such leads to decreased transport of glutamate and acivicin thus leading to acivicin resistance. Author summary: Studies of drug resistance in the protozoan parasites of the genus Leishmania have been helpful in revealing biochemical pathways as potential drug targets. Here we report on the characterization at the genomics and metabolomics levels of a L. tarentolae strain made resistant to acivicin, an analogue of glutamine with activity against this parasite. We found that resistance to acivicin is accompanied by a reduced uptake and intracellular levels of glutamate and that both are expected to share the same transporter. Through gene overexpression and disruption studies we identified the main amino acid transporter responsible for glutamate uptake. [ABSTRACT FROM AUTHOR]

Protein posttranslational modifications (PTMs) regulate the biological processes of human diseases by genetic code expansion and cellular pathophysiology regulation; however, system-wide changes in PTM levels in the intracerebral hemorrhage (ICH) brain remain poorly understood. Succinylation refers to a major PTM during the regulation of multiple biological processes. In this study, according to the methods of quantitative succinyllysine proteomics based on high-resolution mass spectrometry, we investigated ICH-associated brain protein succinyllysine modifications and obtained 3,680 succinylated sites and quantified around 3,530 sites. Among them, 25 succinyllysine sites on 23 proteins were upregulated (hypersuccinylated), whereas 13 succinyllysine sites on 12 proteins were downregulated (hyposuccinylated) following ICH. The cell component enrichment analysis of these succinylproteins with significant changes showed that 58.3% of the hyposuccinylated proteins were observed in the mitochondria, while the hyper-succinylproteins located in mitochondria decreased in the percentage to about 35% in ICH brains with a concomitant increase in the percentage of cytoplasm to 30.4%. Further bioinformatic analysis showed that the succinylproteins were mostly mitochondria and synapse-related subcellular located and involved in many pathophysiological processes, like metabolism, synapse working, and ferroptosis. Moreover, the integrative analysis of our succinylproteomics data and previously published transcriptome data showed that the mRNAs matched by most differentially succinylated proteins were especially highly expressed in neurons, endothelial cells, and astrocytes. Our study uncovers some succinylation-affected processes and pathways in response to ICH brains and gives us novel insights into understanding pathophysiological processes of brain injury caused by ICH. [ABSTRACT FROM AUTHOR]

Mutations in the PINK1 and PRKN genes are the most common cause of early-onset familial Parkinson disease. These genes code for the PINK1 and Parkin proteins, respectively, which are involved in the degradation of dysfunctional mitochondria through mitophagy. An early step in PINK1 –Parkin mediated mitophagy is the ubiquitination of the mitofusin proteins MFN1 and -2. The ubiquitination of MFN1 and -2 in patient samples may therefore serve as a biomarker to determine the functional effects of PINK1 and PRKN mutations, and to screen idiopathic patients for potential mitophagy defects. We aimed to characterise the expression of the PINK1 –Parkin mitophagy machinery in peripheral blood mononuclear cells (PBMCs) and assess if these cells could serve as a platform to evaluate mitophagy via analysis of MFN1 and -2 ubiquitination. Mitophagy was induced through mitochondrial depolarisation by treatment with the protonophore CCCP and ubiquitinated MFN proteins were analysed by western blotting. In addition, PINK1 and PRKN mRNA and protein expression levels were characterised with reverse transcriptase quantitative PCR and western blotting, respectively. Whilst CCCP treatment led to MFN ubiquitination in primary fibroblasts, SH-SY5Y neuroblastoma cells and Jurkat leukaemic cells, treatment of PBMCs did not induce ubiquitination of MFN. PRKN mRNA and protein was readily detectable in PBMCs at comparable levels to those observed in Jurkat and fibroblast cells. In contrast, PINK1 protein was undetectable and PINK1 mRNA levels were remarkably low in control PBMCs. Our findings suggest that the PINK1 –Parkin mitophagy signalling pathway is not functional in PBMCs. Therefore, PBMCs are not a suitable biosample for analysis of mitophagy function in Parkinson disease patients. [ABSTRACT FROM AUTHOR]

We sequenced maxicircles from T. cruzi strains representative of the species evolutionary diversity by using long-read sequencing, which allowed us to uncollapse their repetitive regions, finding that their real lengths range from 35 to 50 kb. T. cruzi maxicircles have a common architecture composed of four regions: coding region (CR), AT-rich region, short (SR) and long repeats (LR). Distribution of genes, both in order and in strand orientation are conserved, being the main differences the presence of deletions affecting genes coding for NADH dehydrogenase subunits, reinforcing biochemical findings that indicate that complex I is not functional in T. cruzi. Moreover, the presence of complete minicircles into maxicircles of some strains lead us to think about the origin of minicircles. Finally, a careful phylogenetic analysis was conducted using coding regions of maxicircles from up to 29 strains, and 1108 single copy nuclear genes from all of the DTUs, clearly establishing that taxonomically T. cruzi is a complex of species composed by group 1 that contains clades A (TcI), B (TcIII)and D (TcIV), and group 2 (1 and 2 do not coincide with groups I and II described decades ago) containing clade C (TcII), being all hybrid strains of the BC type. Three variants of maxicircles exist in T. cruzi: a, b and c, in correspondence with clades A, B, and C from mitochondrial phylogenies. While A and C carry maxicircles a and c respectively, both clades B and D carry b maxicircle variant; hybrid strains also carry the b- variant. We then propose a new nomenclature that is self-descriptive and makes use of both the phylogenetic relationships and the maxicircle variants present in T. cruzi. [ABSTRACT FROM AUTHOR]

Human immunogenetic variation in the form of HLA and KIR types has been shown to be strongly associated with a multitude of immune-related phenotypes. However, association studies involving immunogenetic loci most commonly involve simple analyses of classical HLA allelic diversity, resulting in limitations regarding the interpretability and reproducibility of results. We here present MiDAS, a comprehensive R package for immunogenetic data transformation and statistical analysis. MiDAS recodes input data in the form of HLA alleles and KIR types into biologically meaningful variables, allowing HLA amino acid fine mapping, analyses of HLA evolutionary divergence as well as experimentally validated HLA-KIR interactions. Further, MiDAS enables comprehensive statistical association analysis workflows with phenotypes of diverse measurement scales. MiDAS thus closes the gap between the inference of immunogenetic variation and its efficient utilization to make relevant discoveries related to immune and disease biology. It is freely available under a MIT license. Author summary: Genetic association studies of complex traits often yield highly significant associations in genetic loci coding for HLA or KIR genes, which have central functions for immune responses. Although the roles of these genes, for example in antigen presentation or immune signaling, are well established, their extreme degree of variability makes it challenging to infer mechanistic hypotheses. Starting with HLA or KIR typing data, our software tool MiDAS facilitates statistical association testing, but also recodes and groups immunogenetic information according to function and validated biological interactions. For instance, we can test for association on the level of the actual amino acid sequence, investigating whether an association is strongest for amino acids that determine whether or not a given antigen can be presented by a HLA protein. We can also group HLA alleles according to their interaction with specific KIR on Natural Killer (NK) cells, and a significant association of such interactions might implicate NK cells in our phenotype of interest. In summary, MiDAS offers straightforward workflows for the analysis of immunogenetic data from discovery to functional fine-mapping. [ABSTRACT FROM AUTHOR]