Search

Abstract Background Epigenome‐wide association studies have identified multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. This study aimed to test the hypothesis that an alcohol consumption epigenetic risk score (ERS) is associated with blood pressure (BP) traits. Results We implemented an ERS based on a previously reported epigenetic signature of 144 alcohol-associated CpGs in meta-analysis of participants of European ancestry. We found a one-unit increment of ERS was associated with eleven drinks of alcohol consumed per day, on average, across several cohorts (p

Abstract Background Dyslipidemia, which is characterized by an unfavorable lipid profile, is a key risk factor for cardiovascular disease (CVD). Understanding the relationships between epigenetic aging and lipid levels may help guide early prevention and treatment efforts for dyslipidemia. Methods We used weighted linear regression to cross-sectionally investigate the associations between five measures of epigenetic age acceleration estimated from whole blood DNA methylation (HorvathAge Acceleration, HannumAge Acceleration, PhenoAge Acceleration, GrimAge Acceleration, and DunedinPACE) and four blood lipid measures (total cholesterol (TC), LDL-C, HDL-C, and triglycerides (TG)) in 3,813 participants (mean age = 70 years) from the Health and Retirement Study (HRS). As a sensitivity analysis, we examined the same associations in participants who fasted prior to the blood draw (n = 2,531) and in participants who did not take lipid-lowering medication (n = 1,869). Using interaction models, we also examined whether demographic factors including age, sex, and educational attainment modified the relationships between epigenetic age acceleration and blood lipids. Results After adjusting for age, race/ethnicity, sex, fasting status, and lipid-lowering medication use, greater epigenetic age acceleration was associated with lower TC, HDL-C, and LDL-C, and higher TG (p

Abstract In silico identification of viral anti-CRISPR proteins (Acrs) has relied largely on the guilt-by-association method using known Acrs or anti-CRISPR associated proteins (Acas) as the bait. However, the low number and limited spread of the characterized archaeal Acrs and Aca hinders our ability to identify Acrs using guilt-by-association. Here, based on the observation that the few characterized archaeal Acrs and Aca are transcribed immediately post viral infection, we hypothesize that these genes, and many other unidentified anti-defense genes (ADG), are under the control of conserved regulatory sequences including a strong promoter, which can be used to predict anti-defense genes in archaeal viruses. Using this consensus sequence based method, we identify 354 potential ADGs in 57 archaeal viruses and 6 metagenome-assembled genomes. Experimental validation identified a CRISPR subtype I-A inhibitor and the first virally encoded inhibitor of an archaeal toxin-antitoxin based immune system. We also identify regulatory proteins potentially akin to Acas that can facilitate further identification of ADGs combined with the guilt-by-association approach. These results demonstrate the potential of regulatory sequence analysis for extensive identification of ADGs in viruses of archaea and bacteria.

Abstract Thermus thermophilus bacteriophage P23-45 encodes a giant 5,002-residue tail tape measure protein (TMP) that defines the length of its extraordinarily long tail. Here, we show that the N-terminal portion of P23-45 TMP is an unusual RNA polymerase (RNAP) homologous to cellular RNAPs. The TMP-fused virion RNAP transcribes pre-early phage genes, including a gene that encodes another, non-virion RNAP, that transcribes early and some middle phage genes. We report the crystal structures of both P23-45 RNAPs. The non-virion RNAP has a crab-claw-like architecture. By contrast, the virion RNAP adopts a unique flat structure without a clamp. Structure and sequence comparisons of the P23-45 RNAPs with other RNAPs suggest that, despite the extensive functional differences, the two P23-45 RNAPs originate from an ancient gene duplication in an ancestral phage. Our findings demonstrate striking adaptability of RNAPs that can be attained within a single virus species.

A comprehensive census of McrBC systems, among the most common forms of prokaryotic Type IV restriction systems, followed by phylogenetic analysis, reveals their enormous abundance in diverse prokaryotes and a plethora of genomic associations. We focus on a previously uncharacterized branch, which we denote coiled-coil nuclease tandems (CoCoNuTs) for their salient features: the presence of extensive coiled-coil structures and tandem nucleases. The CoCoNuTs alone show extraordinary variety, with three distinct types and multiple subtypes. All CoCoNuTs contain domains predicted to interact with translation system components, such as OB-folds resembling the SmpB protein that binds bacterial transfer-messenger RNA (tmRNA), YTH-like domains that might recognize methylated tmRNA, tRNA, or rRNA, and RNA-binding Hsp70 chaperone homologs, along with RNases, such as HEPN domains, all suggesting that the CoCoNuTs target RNA. Many CoCoNuTs might additionally target DNA, via McrC nuclease homologs. Additional restriction systems, such as Type I RM, BREX, and Druantia Type III, are frequently encoded in the same predicted superoperons. In many of these superoperons, CoCoNuTs are likely regulated by cyclic nucleotides, possibly, RNA fragments with cyclic termini, that bind associated CARF (CRISPR-Associated Rossmann Fold) domains. We hypothesize that the CoCoNuTs, together with the ancillary restriction factors, employ an echeloned defense strategy analogous to that of Type III CRISPR-Cas systems, in which an immune response eliminating virus DNA and/or RNA is launched first, but then, if it fails, an abortive infection response leading to PCD/dormancy via host RNA cleavage takes over.

ABSTRACT Endosomal sorting complexes required for transport (ESCRT) play key roles in protein sorting between membrane-bounded compartments of eukaryotic cells. Homologs of many ESCRT components are identifiable in various groups of archaea, especially in Asgardarchaeota, the archaeal phylum that is currently considered to include the closest relatives of eukaryotes, but not in bacteria. We performed a comprehensive search for ESCRT protein homologs in archaea and reconstructed ESCRT evolution using the phylogenetic tree of Vps4 ATPase (ESCRT IV) as a scaffold and using sensitive protein sequence analysis and comparison of structural models to identify previously unknown ESCRT proteins. Several distinct groups of ESCRT systems in archaea outside of Asgard were identified, including proteins structurally similar to ESCRT-I and ESCRT-II, and several other domains involved in protein sorting in eukaryotes, suggesting an early origin of these components. Additionally, distant homologs of CdvA proteins were identified in Thermoproteales which are likely components of the uncharacterized cell division system in these archaea. We propose an evolutionary scenario for the origin of eukaryotic and Asgard ESCRT complexes from ancestral building blocks, namely, the Vps4 ATPase, ESCRT-III components, wH (winged helix-turn-helix fold) and possibly also coiled-coil, and Vps28-like domains. The last archaeal common ancestor likely encompassed a complex ESCRT system that was involved in protein sorting. Subsequent evolution involved either simplification, as in the TACK superphylum, where ESCRT was co-opted for cell division, or complexification as in Asgardarchaeota. In Asgardarchaeota, the connection between ESCRT and the ubiquitin system that was previously considered a eukaryotic signature was already established.IMPORTANCEAll eukaryotic cells possess complex intracellular membrane organization. Endosomal sorting complexes required for transport (ESCRT) play a central role in membrane remodeling which is essential for cellular functionality in eukaryotes. Recently, it has been shown that Asgard archaea, the archaeal phylum that includes the closest known relatives of eukaryotes, encode homologs of many components of the ESCRT systems. We employed protein sequence and structure comparisons to reconstruct the evolution of ESCRT systems in archaea and identified several previously unknown homologs of ESCRT subunits, some of which can be predicted to participate in cell division. The results of this reconstruction indicate that the last archaeal common ancestor already encoded a complex ESCRT system that was involved in protein sorting. In Asgard archaea, ESCRT systems evolved toward greater complexity, and in particular, the connection between ESCRT and the ubiquitin system that was previously considered a eukaryotic signature was established.

Introduction: Rate-limiting enzymes (RLEs) are innate slow points in metabolic pathways, and many function in bio-processes related to nutrient sensing. Many RLEs carry causal mutations relevant to inherited metabolic disorders. Because the activity of RLEs in cardiovascular health is poorly characterized, our objective was to assess their involvement in cardiometabolic health and disease and where altered biophysical and biochemical functions can promote disease. Methods: A dataset of 380 human RLEs was compared to protein and gene datasets for factors likely to contribute to cardiometabolic disease, including proteins showing significant age-related altered expression in blood and genetic loci with variants that associate with common cardiometabolic phenotypes. The biochemical reactions catalyzed by RLEs were evaluated for metabolites enriched in RLE subsets associating with various cardiometabolic phenotypes. Most significance tests were based on Z-score enrichment converted to p values with a normal distribution function. Results: Of 380 RLEs analyzed, 112 function in mitochondria, and 53 are assigned to inherited metabolic disorders. There was a depletion of RLE proteins known as aging biomarkers. At the gene level, RLEs were assessed for common genetic variants that associated with important cardiometabolic traits of LDL-cholesterol or any of the five outcomes pertinent to metabolic syndrome. This revealed several RLEs with links to cardiometabolic traits, from a minimum of 26 for HDL-cholesterol to a maximum of 45 for plasma glucose. Analysis of these GWAS-linked RLEs for enrichment of the molecular constituents of the catalyzed reactions disclosed a number of significant phenotype-metabolite links. These included blood pressure with acetate (p = 2.2 × 10−4) and NADP+ (p = 0.0091), plasma HDL-cholesterol and triglyceride with diacylglycerol (p = 2.6 × 10−5, 6.4 × 10−5, respectively) and diolein (p = 2.2 × 10−6, 5.9 × 10−6), and waist circumference with d-glucosamine-6-phosphate (p = 1.8 × 10−4). Conclusion: In the context of cardiometabolic health, aging, and disease, these results highlight key diet-derived metabolites that are central to specific rate-limited processes that are linked to cardiometabolic health. These metabolites include acetate and diacylglycerol, pertinent to blood pressure and triglycerides, respectively, as well as diacylglycerol and HDL-cholesterol.

ABSTRACT The identification of microbial genes essential for survival as those with lethal knockout phenotype (LKP) is a common strategy for functional interrogation of genomes. However, interpretation of the LKP is complicated because a substantial fraction of the genes with this phenotype remains poorly functionally characterized. Furthermore, many genes can exhibit LKP not because their products perform essential cellular functions but because their knockout activates the toxicity of other genes (conditionally essential genes). We analyzed the sets of LKP genes for two archaea, Methanococcus maripaludis and Sulfolobus islandicus, using a variety of computational approaches aiming to differentiate between essential and conditionally essential genes and to predict at least a general function for as many of the proteins encoded by these genes as possible. This analysis allowed us to predict the functions of several LKP genes including previously uncharacterized subunit of the GINS protein complex with an essential function in genome replication and of the KEOPS complex that is responsible for an essential tRNA modification as well as GRP protease implicated in protein quality control. Additionally, several novel antitoxins (conditionally essential genes) were predicted, and this prediction was experimentally validated by showing that the deletion of these genes together with the adjacent genes apparently encoding the cognate toxins caused no growth defect. We applied principal component analysis based on sequence and comparative genomic features showing that this approach can separate essential genes from conditionally essential ones and used it to predict essential genes in other archaeal genomes.IMPORTANCEOnly a relatively small fraction of the genes in any bacterium or archaeon is essential for survival as demonstrated by the lethal effect of their disruption. The identification of essential genes and their functions is crucial for understanding fundamental cell biology. However, many of the genes with a lethal knockout phenotype remain poorly functionally characterized, and furthermore, many genes can exhibit this phenotype not because their products perform essential cellular functions but because their knockout activates the toxicity of other genes. We applied state-of-the-art computational methods to predict the functions of a number of uncharacterized genes with the lethal knockout phenotype in two archaeal species and developed a computational approach to predict genes involved in essential functions. These findings advance the current understanding of key functionalities of archaeal cells.

Abstract Several groups of bacteria have complex life cycles involving cellular differentiation and multicellular structures. For example, actinobacteria of the genus Streptomyces form multicellular vegetative hyphae, aerial hyphae, and spores. However, similar life cycles have not yet been described for archaea. Here, we show that several haloarchaea of the family Halobacteriaceae display a life cycle resembling that of Streptomyces bacteria. Strain YIM 93972 (isolated from a salt marsh) undergoes cellular differentiation into mycelia and spores. Other closely related strains are also able to form mycelia, and comparative genomic analyses point to gene signatures (apparent gain or loss of certain genes) that are shared by members of this clade within the Halobacteriaceae. Genomic, transcriptomic and proteomic analyses of non-differentiating mutants suggest that a Cdc48-family ATPase might be involved in cellular differentiation in strain YIM 93972. Additionally, a gene encoding a putative oligopeptide transporter from YIM 93972 can restore the ability to form hyphae in a Streptomyces coelicolor mutant that carries a deletion in a homologous gene cluster (bldKA-bldKE), suggesting functional equivalence. We propose strain YIM 93972 as representative of a new species in a new genus within the family Halobacteriaceae, for which the name Actinoarchaeum halophilum gen. nov., sp. nov. is herewith proposed. Our demonstration of a complex life cycle in a group of haloarchaea adds a new dimension to our understanding of the biological diversity and environmental adaptation of archaea.

Malnutrition is one of the major factors of bone and cartilage disorders. Pacific cod (Gadus macrocephalus) processing waste is a cheap and highly promising source of bioactive substances, including collagen-derived peptides and amino acids, for bone and cartilage structure stabilization. The addition of these substances to a functional drink is one of the ways to achieve their fast intestinal absorption. Collagen hydrolysate was obtained via enzymatic hydrolysis, ultrafiltration, freeze-drying, and grinding to powder. The lyophilized hydrolysate was a light gray powder with high protein content (>90%), including collagen (about 85% of total protein) and a complete set of essential and non-essential amino acids. The hydrolysate had no observed adverse effect on human mesenchymal stem cell morphology, viability, or proliferation. The hydrolysate was applicable as a protein food supply or a structure-forming food component due to the presence of collagen fiber fragments. An isotonic fitness drink (osmolality 298.1 ± 2.1 mOsm/L) containing hydrolysate and vitamin C as a cofactor in collagen biosynthesis was prepared. The addition of the hydrolysate did not adversely affect its organoleptic parameters. The production of such functional foods and drinks is one of the beneficial ways of fish processing waste utilization.

The female patient with extracapillary pauci-immune AncA associated glomerulonephritis having antibodies against myeloperoxidase and dialysis-requiring renal failure is described. The computed tomography revealed the mucous membrane thickening in the maxillary sinus. The pulse dose glucocorticosteroid and cyclophosphamide therapy contributed to the decrease in the levels of antibodies against myeloperoxidase and improvement of the nitrogen excretion function of the kidneys.

Abstract Objectives A previously developed decision model to prioritize surgical procedures in times of scarce surgical capacity used quality of life (QoL) primarily derived from experts in one center. These estimates are key input of the model, and might be more context-dependent than the other input parameters (age, survival). The aim of this study was to validate our model by replicating these QoL estimates. Methods The original study estimated QoL of patients in need of commonly performed procedures in live expert-panel meetings. This study replicated this procedure using a web-based Delphi approach in a different hospital. The new QoL scores were compared with the original scores using mixed effects linear regression. The ranking of surgical procedures based on combined QoL values from the validation and original study was compared to the ranking based solely on the original QoL values. Results The overall mean difference in QoL estimates between the validation study and the original study was − 0.11 (95% CI: -0.12 - -0.10). The model output (DALY/month delay) based on QoL data from both studies was similar to the model output based on the original data only: The Spearman’s correlation coefficient between the ranking of all procedures before and after including the new QoL estimates was 0.988. Discussion Even though the new QoL estimates were systematically lower than the values from the original study, the ranking for urgency based on health loss per unit of time delay of procedures was consistent. This underscores the robustness and generalizability of the decision model for prioritization of surgical procedures.

Hemangioblastomas are benign tumors that develop from the vessels of the central nervous system and can be a manifestation of autosomal dominant von Hippel-Lindau disease. Statistically, they account for 1.5‒2.5% of all intracranial tumors and 2‒15% of spinal cord tumor lesions. There are very few publications on the intramedullary localization of these neoplasms.The patient, 45 years old, a serviceman, presented with complaints of headache, slight unsteadiness of gait, as well as slight weakness in the right extremities, more pronounced in the right upper extremity, periodic numbness of the upper extremities, which progressed and made further service impossible. On neurological examination: pupils D=S, light reflexes were brisk, eye movement was fully preserved and horizontal nystagmus. BNI - PS - I, BNI - NS - I. HB - I. GR - I. Swallowing and phonation were fully preserved. There was a slight hemiparesis on the right. Hemihypesthesia on the right was more prominent in the upper extremity. Ataxia of mixed genesis. Pelvic organs function was preserved. Periodic constipation for up to 7 days. Magnetic resonance imaging revealed a multifocal brain lesion. Supratentorially, a cystic mass measuring 2.60×2.12×2.14 cm with a solid component up to 1.5 cm in the diameter was detected in the of the thickened corpus callosum. Intramedullary cystic-solid lesion of the medulla oblongata extending to the cervical spinal cord with conventional dimensions of the solid component 1.76×1.23×1.57 cm and the cystic component was 1.52×1.62×1.22 cm. Magnetic resonance imaging of the cervicothoracic region of the spinal cord revealed significant hydromyelitic expansion of the central canal from the C2 level to the Th3 level (up to 10mm in the diameter). Endovascular embolization of the neoplasm with a liquid embolic agent (Phil) and following microsurgical en bloc tumor resection were performed.Hemangioblastomas with an intramedullary location are extremely difficult and risky for the surgical removal. The presence of a cystic component in the hemangioblastoma structure or perifocally gives a chance to remove such a neoplasm avoiding risk of functional deterioration. Preoperative endovascular obliteration of hemangioblastoma vascularity is considered as an effective measure, although it is associated with a risk of cerebral ischemia in corresponding brain structures. Another crucial issues of intramedullary hemangioblastoma surgery includes multimodal intraoperative neuromonitoring and maximum possible intraoperative magnification to minimize injury of adjacent brain structures. Present practice shows that multiple intramedullary cysts in the spinal cord commonly regressing after neoplasm removal and leads to improvement of neurological deficit within a relatively short period of time.Conclusions. Hemangioblastomas with intramedullary growth and perifocal cyst can be surgically removed with a good functional result and comprehensive approach includes preoperative selective angiography with endovascular embolization, multimodal intraoperative neuromonitoring and appropriate microsurgical technique.

Abstract Background The burden of the COVID-19 pandemic resulted in a reduction of available health care capacity for regular care. To guide prioritisation of semielective surgery in times of scarcity, we previously developed a decision model to quantify the expected health loss due to delay of surgery, in an academic hospital setting. The aim of this study is to validate our decision model in a nonacademic setting and include additional elective surgical procedures. Methods In this study, we used the previously published three-state cohort state-transition model, to evaluate the health effects of surgery postponement for 28 surgical procedures commonly performed in nonacademic hospitals. Scientific literature and national registries yielded nearly all input parameters, except for the quality of life (QoL) estimates which were obtained from experts using the Delphi method. Two expert panels, one from a single nonacademic hospital and one from different nonacademic hospitals in the Netherlands, were invited to estimate QoL weights. We compared estimated model results (disability adjusted life years (DALY)/month of surgical delay) based on the QoL estimates from the two panels by calculating the mean difference and the correlation between the ranks of the different surgical procedures. The eventual model was based on the combined QoL estimates from both panels. Results Pacemaker implantation was associated with the most DALY/month of surgical delay (0.054 DALY/month, 95% CI: 0.025–0.103) and hemithyreoidectomy with the least DALY/month (0.006 DALY/month, 95% CI: 0.002–0.009). The overall mean difference of QoL estimates between the two panels was 0.005 (95% CI -0.014–0.004). The correlation between ranks was 0.983 (p

The cryo-EM structure of the IscB–ωRNA–target DNA complex provide insights into the ωRNA-guided DNA cleavage mechanism by IscB and the evolution of the type II CRISPR-Cas9 effector complexes

Abstract Background Evolutionary rate is a key characteristic of gene families that is linked to the functional importance of the respective genes as well as specific biological functions of the proteins they encode. Accurate estimation of evolutionary rates is a challenging task that requires precise phylogenetic analysis. Here we present an easy to estimate protein family level measure of sequence variability based on alignment column homogeneity in multiple alignments of protein sequences from Clade-Specific Clusters of Orthologous Genes (csCOGs). Results We report genome-wide estimates of variability for 8 diverse groups of bacteria and archaea and investigate the connection between variability and various genomic and biological features. The variability estimates are based on homogeneity distributions across amino acid sequence alignments and can be obtained for multiple groups of genomes at minimal computational expense. About half of the variance in variability values can be explained by the analyzed features, with the greatest contribution coming from the extent of gene paralogy in the given csCOG. The correlation between variability and paralogy appears to originate, primarily, not from gene duplication, but from acquisition of distant paralogs and xenologs, introducing sequence variants that are more divergent than those that could have evolved in situ during the lifetime of the given group of organisms. Both high-variability and low-variability csCOGs were identified in all functional categories, but as expected, proteins encoded by integrated mobile elements as well as proteins involved in defense functions and cell motility are, on average, more variable than proteins with housekeeping functions. Additionally, using linear discriminant analysis, we found that variability and fraction of genomes carrying a given gene are the two variables that provide the best prediction of gene essentiality as compared to the results of transposon mutagenesis in Sulfolobus islandicus. Conclusions Variability, a measure of sequence diversity within an alignment relative to the overall diversity within a group of organisms, offers a convenient proxy for evolutionary rate estimates and is informative with respect to prediction of functional properties of proteins. In particular, variability is a strong predictor of gene essentiality for the respective organisms and indicative of sub- or neofunctionalization of paralogs.

Abstract Background Bacteria and archaea produce an enormous diversity of modified peptides that are involved in various forms of inter-microbial conflicts or communication. A vast class of such peptides are Ribosomally synthesized, Postranslationally modified Peptides (RiPPs), and a major group of RiPPs are graspetides, so named after ATP-grasp ligases that catalyze the formation of lactam and lactone linkages in these peptides. The diversity of graspetides, the multiple proteins encoded in the respective Biosynthetic Gene Clusters (BGCs) and their evolution have not been studied in full detail. In this work, we attempt a comprehensive analysis of the graspetide-encoding BGCs and report a variety of novel graspetide groups as well as ancillary proteins implicated in graspetide biosynthesis and expression. Results We compiled a comprehensive, manually curated set of graspetides that includes 174 families including 115 new families with distinct patterns of amino acids implicated in macrocyclization and further modification, roughly tripling the known graspetide diversity. We derived signature motifs for the leader regions of graspetide precursors that could be used to facilitate graspetide prediction. Graspetide biosynthetic gene clusters and specific precursors were identified in bacterial divisions not previously known to encode RiPPs, in particular, the parasitic and symbiotic bacteria of the Candidate phyla radiation. We identified Bacteroides-specific biosynthetic gene clusters (BGC) that include remarkable diversity of graspetides encoded in the same loci which predicted to be modified by the same ATP-grasp ligase. We studied in details evolution of recently characterized chryseoviridin BGCs and showed that duplication and horizonal gene exchange both contribute to the diversification of the graspetides during evolution. Conclusions We demonstrate previously unsuspected diversity of graspetide sequences, even those associated with closely related ATP-grasp enzymes. Several previously unnoticed families of proteins associated with graspetide biosynthetic gene clusters are identified. The results of this work substantially expand the known diversity of RiPPs and can be harnessed to further advance approaches for their identification.

Neuraminidase (NA)-based immunity could reduce the harmful impact of novel antigenic variants of influenza viruses. The detection of neuraminidase-inhibiting (NI) antibodies in parallel with anti-hemagglutinin (HA) antibodies may enhance research on the immunogenicity and duration of antibody responses to influenza vaccines. To assess anti-NA antibodies after vaccination with seasonal inactivated influenza vaccines, we used the enzyme-linked lectin assay, and anti-HA antibodies were detected in the hemagglutination inhibition assay. The dynamics of the anti-NA antibody response differed depending on the virus subtype: antibodies to A/H3N2 virus neuraminidase increased later than antibodies to A/H1N1pdm09 subtype neuraminidase and persisted longer. In contrast to HA antibodies, the fold increase in antibody titers to NA after vaccination poorly depended on the preexisting level. At the same time, NA antibody levels after vaccination directly correlated with titers before vaccination. A difference was found in response to NA antigen between split and subunit-adjuvanted vaccines and in NA functional activity in the vaccine formulations.

Many microorganisms feed on the tissue and recalcitrant bone materials from dead animals, however little is known about the collaborative effort and characteristics of their enzymes. In this study, microbial metagenomes from symbionts of the marine bone-dwelling worm Osedax mucofloris, and from microbial biofilms growing on experimentally deployed bone surfaces were screened for specialized bone-degrading enzymes. A total of 2,043 taxonomically (closest match within 40 phyla) and functionally (1 proteolytic and 9 glycohydrolytic activities) diverse and non-redundant sequences (median pairwise identity of 23.6%) encoding such enzymes were retrieved. The taxonomic assignation and the median identity of 72.2% to homologous proteins reflect microbial and functional novelty associated to a specialized bone-degrading marine community. Binning suggests that only one generalist hosting all ten targeted activities, working in synergy with multiple specialists hosting a few or individual activities. Collagenases were the most abundant enzyme class, representing 48% of the total hits. A total of 47 diverse enzymes, representing 8 hydrolytic activities, were produced in Escherichia coli, whereof 13 were soluble and active. The biochemical analyses revealed a wide range of optimal pH (4.0–7.0), optimal temperature (5–65 °C), and of accepted substrates, specific to each microbial enzyme. This versatility may contribute to a high environmental plasticity of bone-degrading marine consortia that can be confronted to diverse habitats and bone materials. Through bone-meal degradation tests, we further demonstrated that some of these enzymes, particularly those from Flavobacteriaceae and Marinifilaceae, may be an asset for development of new value chains in the biorefinery industry.

CRISPR-Cas is a host adaptive immunity system and viruses harbor diverse anti-CRISPR proteins (Acrs). Here, the authors develop a random forest machine-learning approach to predict Acrs, identifying 2500 candidate Acr families, which expand the current repertoire of predicted Acrs by two orders of magnitude.

Abstract Macrostomum lignano is a free-living flatworm that is emerging as an attractive experimental animal for research on a broad range of biological questions. One feature setting it apart from other flatworms is the successful establishment of transgenesis methods, facilitated by a steady supply of eggs in the form of single-cell zygotes that can be readily manipulated. This, in combination with the transparency of the animal and its small size, creates practical advantages for imaging and fluorescence-activated cell sorting in studies related to stem cell biology and regeneration. M. lignano can regenerate most of its body parts, including the germline, thanks to the neoblasts, which represent the flatworm stem cell system. Interestingly, neoblasts seem to have a high capacity of cellular maintenance, as M. lignano can survive up to 210 Gy of γ-irradiation, and partially offset the negative consequence of ageing. As a non-self-fertilizing simultaneous hermaphrodite that reproduces in a sexual manner, M. lignano is also used to study sexual selection and other evolutionary aspects of sexual reproduction. Work over the past several years has led to the development of molecular resources and tools, including high-quality genome and transcriptome assemblies, transcriptional profiling of the germline and somatic neoblasts, gene knockdown, and in situ hybridization. The increasingly detailed characterization of this animal has also resulted in novel research questions, such as bio-adhesion based on its adhesion-release glands and genome evolution due to its recent whole-genome duplication.

ABSTRACT Denham Harman’s oxidative damage theory identifies superoxide (O2•−) radicals as central agents of aging and radiation injury, with Mn2+-dependent superoxide dismutase (MnSOD) as the principal O2•−-scavenger. However, in the radiation-resistant nematode Caenorhabditis elegans, the mitochondrial antioxidant enzyme MnSOD is dispensable for longevity, and in the model bacterium Deinococcus radiodurans, it is dispensable for radiation resistance. Many radiation-resistant organisms accumulate small-molecule Mn2+-antioxidant complexes well-known for their catalytic ability to scavenge O2•−, along with MnSOD, as exemplified by D. radiodurans. Here, we report experiments that relate the MnSOD and Mn-antioxidant content to aging and oxidative stress resistances and which indicate that C. elegans, like D. radiodurans, may rely on Mn-antioxidant complexes as the primary defense against reactive oxygen species (ROS). Wild-type and ΔMnSOD D. radiodurans and C. elegans were monitored for gamma radiation sensitivities over their life spans while gauging Mn2+-antioxidant content by electron paramagnetic resonance (EPR) spectroscopy, a powerful new approach to determining the in vivo Mn-antioxidant content of cells as they age. As with D. radiodurans, MnSOD is dispensable for radiation survivability in C. elegans, which hyperaccumulates Mn-antioxidants exceptionally protective of proteins. Unexpectedly, ΔMnSOD mutants of both the nematodes and bacteria exhibited increased gamma radiation survival compared to the wild-type. In contrast, the loss of MnSOD renders radiation-resistant bacteria sensitive to atmospheric oxygen during desiccation. Our results support the concept that the disparate responses to oxidative stress are explained by the accumulation of Mn-antioxidant complexes which protect, complement, and can even supplant MnSOD. IMPORTANCE The current theory of cellular defense against oxidative damage identifies antioxidant enzymes as primary defenders against ROS, with MnSOD being the preeminent superoxide (O2•−) scavenger. However, MnSOD is shown to be dispensable both for radiation resistance and longevity in model organisms, the bacterium Deinococcus radiodurans and the nematode Caenorhabditis elegans. Measured by electron paramagnetic resonance (EPR) spectroscopy, small-molecule Mn-antioxidant content was shown to decline in unison with age-related decreases in cell proliferation and radioresistance, which again are independent of MnSOD presence. Most notably, the Mn-antioxidant content of C. elegans drops precipitously in the last third of its life span, which links with reports that the steady-state level of oxidized proteins increases exponentially during the last third of the life span in animals. This leads us to propose that global responses to oxidative stress must be understood through an extended theory that includes small-molecule Mn-antioxidants as potent O2•−-scavengers that complement, and can even supplant, MnSOD.

The principal biological function of bacterial and archaeal CRISPR systems is RNA-guided adaptive immunity against viruses and other mobile genetic elements (MGEs). These systems show remarkable evolutionary plasticity and functional versatility at multiple levels, including both the defense mechanisms that lead to direct, specific elimination of the target DNA or RNA and those that cause programmed cell death (PCD) or induction of dormancy. This flexibility is also evident in the recruitment of CRISPR systems for nondefense functions. Defective CRISPR systems or individual CRISPR components have been recruited by transposons for RNA-guided transposition, by plasmids for interplasmid competition, and by viruses for antidefense and interviral conflicts. Additionally, multiple highly derived CRISPR variants of yet unknown functions have been discovered. A major route of innovation in CRISPR evolution is the repurposing of diverged repeat variants encoded outside CRISPR arrays for various structural and regulatory functions. The evolutionary plasticity and functional versatility of CRISPR systems are striking manifestations of the ubiquitous interplay between defense and “normal” cellular functions. The CRISPR systems show remarkable functional versatility beyond their principal function as an adaptive immune mechanism. This Essay discusses how derived CRISPR systems have been recruited by transposons on multiple occasions and mediate RNA-guided transposition; derived CRISPR RNAs are frequently recruited for regulatory functions.

Sweet dessert watermelon (Citrullus lanatus) is one of the most important vegetable crops consumed throughout the world. The chemical composition of watermelon provides both high nutritional value and various health benefits. The present manuscript introduces a catalog of 1,679 small molecules occurring in the watermelon and their cheminformatics analysis for diverse features. In this catalog, the phytochemicals are associated with the literature describing their presence in the watermelon plant, and when possible, concentration values in various plant parts (flesh, seeds, leaves, roots, rind). Also cataloged are the chemical classes, molecular weight and formula, chemical structure, and certain physical and chemical properties for each phytochemical. In our view, knowing precisely what is in what we eat, as this catalog does for watermelon, supports both the rationale for certain controlled feeding studies in the field of precision nutrition, and plant breeding efforts for the development of new varieties with enhanced concentrations of specific phytochemicals. Additionally, improved and comprehensive collections of natural products accessible to the public will be especially useful to researchers in nutrition, cheminformatics, bioinformatics, and drug development, among other disciplines.

Molecular mechanisms involved in biological conflicts and self vs nonself recognition in archaea remain poorly characterized. We apply phylogenomic analysis to identify a hypervariable gene module that is widespread among Thermococcales. These loci consist of an upstream gene coding for a large protein containing several immunoglobulin (Ig) domains and unique combinations of downstream genes, some of which also contain Ig domains. In the large Ig domain containing protein, the C-terminal Ig domain sequence is hypervariable, apparently, as a result of recombination between genes from different Thermococcales. To reflect the hypervariability, we denote this gene module VARTIG (VARiable Thermococcales IG). The overall organization of the VARTIG modules is similar to the organization of Polymorphic Toxin Systems (PTS). Archaeal genomes outside Thermococcales encode a variety of Ig domain proteins, but no counterparts to VARTIG and no Ig domains with comparable levels of variability. The specific functions of VARTIG remain unknown but the identified features of this system imply three testable hypotheses: (i) involvement in inter-microbial conflicts analogous to PTS, (ii) role in innate immunity analogous to the vertebrate complement system, and (iii) function in self vs nonself discrimination analogous to the vertebrate Major Histocompatibility Complex. The latter two hypotheses seem to be of particular interest given the apparent analogy to the vertebrate immunity.

The Cas12b family of CRISPR nucleases has been underutilized in mammalian cells due to the high temperature requirement of known members. Here the authors engineer BhCas12b to overcome this limitation for robust and specific genome editing applications in human cells.

Abstract Background The Cas4 family endonuclease is a component of the adaptation module in many variants of CRISPR-Cas adaptive immunity systems. Unlike most of the other Cas proteins, Cas4 is often encoded outside CRISPR-cas loci (solo-Cas4) and is also found in mobile genetic elements (MGE-Cas4). Results As part of our ongoing investigation of CRISPR-Cas evolution, we explored the phylogenomics of the Cas4 family. About 90% of the archaeal genomes encode Cas4 compared to only about 20% of the bacterial genomes. Many archaea encode both the CRISPR-associated form (CAS-Cas4) and solo-Cas4, whereas in bacteria, this combination is extremely rare. The solo-cas4 genes are over-represented in environmental bacteria and archaea with small genomes that typically lack CRISPR-Cas, suggesting that Cas4 could perform uncharacterized defense or repair functions in these microbes. Phylogenomic analysis indicates that both the CRISPR-associated cas4 genes are often transferred horizontally but almost exclusively, as part of the adaptation module. The evolutionary integrity of the adaptation module sharply contrasts the rampant shuffling of CRISPR-cas modules whereby a given variant of the adaptation module can combine with virtually any effector module. The solo-cas4 genes evolve primarily via vertical inheritance and are subject only to occasional horizontal transfer. The selection pressure on cas4 genes does not substantially differ between CAS-Cas4 and solo-cas4, and is close to the genomic median. Thus, cas4 genes, similarly to cas1 and cas2, evolve similarly to ‘regular’ microbial genes involved in various cellular functions, showing no evidence of direct involvement in virus-host arms races. A notable feature of the Cas4 family evolution is the frequent recruitment of cas4 genes by various mobile genetic elements (MGE), particularly, archaeal viruses. The functions of Cas4 in these elements are unknown and potentially might involve anti-defense roles. Conclusions Unlike most of the other Cas proteins, Cas4 family members are as often encoded by stand-alone genes as they are incorporated in CRISPR-Cas systems. In addition, cas4 genes were repeatedly recruited by MGE, perhaps, for anti-defense functions. Experimental characterization of the solo and MGE-encoded Cas4 nucleases is expected to reveal currently uncharacterized defense and anti-defense systems and their interactions with CRISPR-Cas systems.

Abstract Background Minimally invasive specific biomarkers of neurodegenerative diseases (NDs) would facilitate patient selection and disease progression monitoring. We describe the assessment of circulating brain-enriched microRNAs as potential biomarkers for Alzheimer’s disease (AD), frontotemporal dementia (FTD), Parkinson’s disease (PD), and amyotrophic lateral sclerosis (ALS). Methods In this case-control study, the plasma samples were collected from 250 research participants with a clinical diagnosis of AD, FTD, PD, and ALS, as well as from age- and sex-matched control subjects (n = 50 for each group), recruited from 2003 to 2015 at the University of Pennsylvania Health System, including the Alzheimer’s Disease Center, the Parkinson’s Disease and Movement Disorders Center, the Frontotemporal Degeneration Center, and the Amyotrophic Lateral Sclerosis Clinic. Each group was randomly divided into training and confirmation sets of equal size. To evaluate the potential of circulating microRNAs enriched in specific brain regions affected by NDs and present in synapses as biomarkers of NDs, the levels of 37 brain-enriched and inflammation-associated microRNAs in the plasma of all participants were measured using individual qRT-PCR. A “microRNA pair” approach was used for data normalization. Results MicroRNA pairs and their combinations (classifiers) capable of differentiating NDs from control and from each other were defined using independently and jointly analyzed training and confirmation datasets. AD, PD, FTD, and ALS are differentiated from control with accuracy of 0.89, 0.90, 0.88, and 0.83 (AUCs, 0.96, 0.96, 0.94, and 0.93), respectively; NDs are differentiated from each other with accuracy ranging from 0.77 (AUC, 0.87) for AD vs. FTD to 0.93 (AUC, 0.98) for AD vs. ALS. The data further indicate sex dependence of some microRNA markers. The average increase in accuracy in distinguishing ND from control for all and male/female groups is 0.06; the largest increase is for ALS, from 0.83 for all participants to 0.92/0.98 for male/female participants. Conclusions The work presented here suggests the possibility of developing microRNA-based diagnostics for detection and differentiation of NDs. Larger multicenter clinical studies are needed to further evaluate circulating brain-enriched microRNAs as biomarkers for NDs and to investigate their association with other ND biomarkers in clinical trial settings.

In the absence of complete genomes, the metabolic capabilities of uncultured ARMAN-like archaea have been uncertain. Here, Golyshina et al. apply an enrichment culture technique and find that the ungapped genome of the ARMAN-like archaeon Mia14 has lost key metabolic pathways, suggesting dependence on the host archaeon Cuniculiplasma divulgatum.

Abstract The free-living flatworm Macrostomum lignano is a model organism for evolutionary and developmental biology studies. Recently, an unusual karyotypic diversity was revealed in this species. Specifically, worms are either ‘normal’ 2n = 8, or they are aneuploid with one or two additional large chromosome(s) (i.e. 2n = 9 or 2n = 10, respectively). Aneuploid worms did not show visible behavioral or morphological abnormalities and were successful in reproduction. In this study, we generated microdissected DNA probes from chromosome 1 (further called MLI1), chromosome 2 (MLI2), and a pair of similar-sized smaller chromosomes (MLI3, MLI4). FISH using these probes revealed that MLI1 consists of contiguous regions homologous to MLI2-MLI4, suggesting that MLI1 arose due to the whole genome duplication and subsequent fusion of one full chromosome set into one large metacentric chromosome. Therefore, one presumably full haploid genome was packed into MLI1, leading to hidden tetraploidy in the M. lignano genome. The study of Macrostomum sp. 8 — a sibling species of M. lignano — revealed that it usually has one additional pair of large chromosomes (2n = 10) showing a high homology to MLI1, thus suggesting hidden hexaploidy in its genome. Possible evolutionary scenarios for the emergence of the M. lignano and Macrostomum sp. 8 genomes are discussed.

Abstract The genetic platforms of Deinococcus species remain the only systems in which massive ionizing radiation (IR)-induced genome damage can be investigated in vivo at exposures commensurate with cellular survival. We report the whole genome sequence of the extremely IR-resistant rod-shaped bacterium Deinococcus ficus KS 0460 and its phenotypic characterization. Deinococcus ficus KS 0460 has been studied since 1987, first under the name Deinobacter grandis, then Deinococcus grandis. The D. ficus KS 0460 genome consists of a 4.019 Mbp sequence (69.7% GC content and 3894 predicted genes) divided into six genome partitions, five of which are confirmed to be circular. Circularity was determined manually by mate pair linkage. Approximately 76% of the predicted proteins contained identifiable Pfam domains and 72% were assigned to COGs. Of all D. ficus KS 0460 proteins, 79% and 70% had homologues in Deinococcus radiodurans ATCC BAA-816 and Deinococcus geothermalis DSM 11300, respectively. The most striking differences between D. ficus KS 0460 and D. radiodurans BAA-816 identified by the comparison of the KEGG pathways were as follows: (i) D. ficus lacks nine enzymes of purine degradation present in D. radiodurans, and (ii) D. ficus contains eight enzymes involved in nitrogen metabolism, including nitrate and nitrite reductases, that D. radiodurans lacks. Moreover, genes previously considered to be important to IR resistance are missing in D. ficus KS 0460, namely, for the Mn-transporter nramp, and proteins DdrF, DdrJ and DdrK, all of which are also missing in Deinococcus deserti. Otherwise, D. ficus KS 0460 exemplifies the Deinococcus lineage.

Abstract We report a comprehensive computational dissection of the domain architecture of the SAMD9 family proteins that are involved in antivirus and antitumor response in humans. We show that the SAMD9 protein family is represented in most animals and also, unexpectedly, in bacteria, in particular actinomycetes. From the N to C terminus, the core SAMD9 family architecture includes DNA/RNA-binding AlbA domain, a variant Sir2-like domain, a STAND-like P-loop NTPase, an array of TPR repeats and an OB-fold domain with predicted RNA-binding properties. Vertebrate SAMD9 family proteins contain the eponymous SAM domain capable of polymerization, whereas some family members from other animals instead contain homotypic adaptor domains of the DEATH superfamily, known as dedicated components of apoptosis networks. Such complex domain architecture is reminiscent of the STAND superfamily NTPases that are involved in various signaling processes, including programmed cell death, in both eukaryotes and prokaryotes. These findings suggest that SAMD9 is a hub of a novel, evolutionarily conserved defense network that remains to be characterized. Reviewers This article was reviewed by Igor B. Zhulin and Mensur Dlakic.

Abstract Background Evolution of bacterial and archaeal genomes is a highly dynamic process that involves intensive loss of genes as well as gene gain via horizontal transfer, with a lesser contribution from gene duplication. The rates of these processes can be estimated by comparing genomes that are linked by an evolutionary tree. These estimated rates of genome dynamics events substantially differ for different functional classes of genes. The genes involved in defense against viruses and other invading DNA are among those that are gained and lost at the highest rates. Results We employed a stochastic birth-and-death model to obtain maximum likelihood estimates of the rates of gain and loss of defense genes in 35 groups of closely related bacterial genomes and one group of archaeal genomes. We find that on average, the defense genes experience 1.4 fold higher flux than the rest of microbial genes. This excessive flux of defense genes over the genomic mean is consistent across diverse microbial groups. The few exceptions include intracellular parasites with small, degraded genomes that possess few defense systems which are more stable than in other microbes. Generally, defense genes follow the previously established pattern of genome dynamics, with gene family loss being about 3 times more common than gain and an order of magnitude more common than expansion or contraction of gene families. Case by case analysis of the evolutionary dynamics of defense genes indicates frequent multiple events in the same locus and widespread involvement of mobile elements in the gain and loss of defense genes. Conclusions Evolution of microbial defense systems is highly dynamic but, notwithstanding the host-parasite arms race, generally follows the same trends that have been established for the rest of the genes. Apart from the paucity and the low flux of defense genes in parasitic bacteria with deteriorating genomes, there is no clear connection between the evolutionary regime of defense systems and microbial life style.

ABSTRACT Dephospho-coenzyme A (dephospho-CoA) kinase (DPCK) catalyzes the ATP-dependent phosphorylation of dephospho-CoA, the final step in coenzyme A (CoA) biosynthesis. DPCK has been identified and characterized in bacteria and eukaryotes but not in archaea. The hyperthermophilic archaeon Thermococcus kodakarensis encodes two homologs of bacterial DPCK and the DPCK domain of eukaryotic CoA synthase, TK1334 and TK2192. We purified the recombinant TK1334 and TK2192 proteins and found that they lacked DPCK activity. Bioinformatic analyses showed that, in several archaea, the uncharacterized gene from arCOG04076 protein is fused with the gene for phosphopantetheine adenylyltransferase (PPAT), which catalyzes the reaction upstream of the DPCK reaction in CoA biosynthesis. This observation suggested that members of arCOG04076, both fused to PPAT and standalone, could be the missing archaeal DPCKs. We purified the recombinant TK1697 protein, a standalone member of arCOG04076 from T. kodakarensis, and demonstrated its GTP-dependent DPCK activity. Disruption of the TK1697 resulted in CoA auxotrophy, indicating that TK1697 encodes a DPCK that contributes to CoA biosynthesis in T. kodakarensis. TK1697 homologs are widely distributed in archaea, suggesting that the arCOG04076 protein represents a novel family of DPCK that is not homologous to bacterial and eukaryotic DPCKs but is distantly related to bacterial and eukaryotic thiamine pyrophosphokinases. We also constructed and characterized gene disruption strains of TK0517 and TK2128, homologs of bifunctional phosphopantothenoylcysteine synthetase-phosphopantothenoylcysteine decarboxylase and PPAT, respectively. Both strains displayed CoA auxotrophy, indicating their contribution to CoA biosynthesis. Taken together with previous studies, the results experimentally validate the entire CoA biosynthesis pathway in T. kodakarensis. IMPORTANCE CoA is utilized in a wide range of metabolic pathways, and its biosynthesis is essential for all life. Pathways for CoA biosynthesis in bacteria and eukaryotes have been established. In archaea, however, the enzyme that catalyzes the final step in CoA biosynthesis, dephospho-CoA kinase (DPCK), had not been identified. In the present study, bioinformatic analyses identified a candidate for the DPCK in archaea, which was biochemically and genetically confirmed in the hyperthermophilic archaeon Thermococcus kodakarensis. Genetic analyses on genes presumed to encode bifunctional phosphopantothenoylcysteine synthetase-phosphopantothenoylcysteine decarboxylase and phosphopantetheine adenylyltransferase confirmed their involvement in CoA biosynthesis. Taken together with previous studies, the results reveal the entire pathway for CoA biosynthesis in a single archaeon and provide insight into the different mechanisms of CoA biosynthesis and their distribution in nature.

ABSTRACT Numerous, diverse, highly variable defense and offense genetic systems are encoded in most bacterial genomes and are involved in various forms of conflict among competing microbes or their eukaryotic hosts. Here we focus on the offense and self-versus-nonself discrimination systems encoded by archaeal genomes that so far have remained largely uncharacterized and unannotated. Specifically, we analyze archaeal genomic loci encoding polymorphic and related toxin systems and ribosomally synthesized antimicrobial peptides. Using sensitive methods for sequence comparison and the “guilt by association” approach, we identified such systems in 141 archaeal genomes. These toxins can be classified into four major groups based on the structure of the components involved in the toxin delivery. The toxin domains are often shared between and within each system. We revisit halocin families and substantially expand the halocin C8 family, which was identified in diverse archaeal genomes and also certain bacteria. Finally, we employ features of protein sequences and genomic locus organization characteristic of archaeocins and polymorphic toxins to identify candidates for analogous but not necessarily homologous systems among uncharacterized protein families. This work confidently predicts that more than 1,600 archaeal proteins, currently annotated as “hypothetical” in public databases, are components of conflict and self-versus-nonself discrimination systems. IMPORTANCE Diverse and highly variable systems involved in biological conflicts and self-versus-nonself discrimination are ubiquitous in bacteria but much less studied in archaea. We performed comprehensive comparative genomic analyses of the archaeal systems that share components with analogous bacterial systems and propose an approach to identify new systems that could be involved in these functions. We predict polymorphic toxin systems in 141 archaeal genomes and identify new, archaea-specific toxin and immunity protein families. These systems are widely represented in archaea and are predicted to play major roles in interactions between species and in intermicrobial conflicts. This work is expected to stimulate experimental research to advance the understanding of poorly characterized major aspects of archaeal biology.

BackgroundEpigenome-wide association studies have identified multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. This study aimed to test the hypothesis that an alcohol consumption epigenetic risk score (ERS) is associated with blood pressure (BP) traits.ResultsWe implemented an ERS based on a previously reported epigenetic signature of 144 alcohol-associated CpGs in meta-analysis of participants of European ancestry. We found a one-unit increment of ERS was associated with eleven drinks of alcohol consumed per day, on average, across several cohorts (p

ABSTRACT All organisms must adapt to ever-changing environmental conditions and accordingly have evolved diverse signal transduction systems. In bacteria, the most abundant networks are built around the two-component signal transduction systems that include histidine kinases and receiver domains. In contrast, eukaryotic signal transduction is dominated by serine/threonine/tyrosine protein kinases. Both of these systems are also found in archaea, but they are not as common and diversified as their bacterial and eukaryotic counterparts, suggesting the possibility that archaea have evolved other, still uncharacterized signal transduction networks. Here we propose a role for KaiC family ATPases, known to be key components of the circadian clock in cyanobacteria, in archaeal signal transduction. The KaiC family is notably expanded in most archaeal genomes, and although most of these ATPases remain poorly characterized, members of the KaiC family have been shown to control archaellum assembly and have been found to be a stable component of the gas vesicle system in Halobacteria. Computational analyses described here suggest that KaiC-like ATPases and their homologues with inactivated ATPase domains are involved in many other archaeal signal transduction pathways and comprise major hubs of complex regulatory networks. We predict numerous input and output domains that are linked to KaiC-like proteins, including putative homologues of eukaryotic DEATH domains that could function as adapters in archaeal signaling networks. We further address the relationships of the archaeal family of KaiC homologues to the bona fide KaiC of cyanobacteria and implications for the existence of a KaiC-based circadian clock apparatus in archaea. IMPORTANCE Little is currently known about signal transduction pathways in Archaea. Recent studies indicate that KaiC-like ATPases, known as key components of the circadian clock apparatus in cyanobacteria, are involved in the regulation of archaellum assembly and, likely, type IV pili and the gas vesicle system in Archaea. We performed comprehensive comparative genomic analyses of the KaiC family. A vast protein interaction network was revealed, with KaiC family proteins as hubs for numerous input and output components, many of which are shared with two-component signal transduction systems. Putative KaiC-based signal transduction systems are predicted to regulate the activities of membrane-associated complexes and individual proteins, such as signal recognition particle and membrane transporters, and also could be important for oxidative stress response regulation. KaiC-centered signal transduction networks are predicted to play major roles in archaeal physiology, and this work is expected to stimulate their experimental characterization.

With the continuously accelerating genome sequencing from diverse groups of archaea and bacteria, accurate identification of gene orthology and availability of readily expandable clusters of orthologous genes are essential for the functional annotation of new genomes. We report an update of the collection of archaeal Clusters of Orthologous Genes (arCOGs) to cover, on average, 91% of the protein-coding genes in 168 archaeal genomes. The new arCOGs were constructed using refined algorithms for orthology identification combined with extensive manual curation, including incorporation of the results of several completed and ongoing research projects in archaeal genomics. A new level of classification is introduced, superclusters that untie two or more arCOGs and more completely reflect gene family evolution than individual, disconnected arCOGs. Assessment of the current archaeal genome annotation in public databases indicates that consistent use of arCOGs can significantly improve the annotation quality. In addition to their utility for genome annotation, arCOGs also are a platform for phylogenomic analysis. We explore this aspect of arCOGs by performing a phylogenomic study of the Thermococci that are traditionally viewed as the basal branch of the Euryarchaeota. The results of phylogenomic analysis that involved both comparison of multiple phylogenetic trees and a search for putative derived shared characters by using phyletic patterns extracted from the arCOGs reveal a likely evolutionary relationship between the Thermococci, Methanococci, and Methanobacteria. The arCOGs are expected to be instrumental for a comprehensive phylogenomic study of the archaea.

ABSTRACT Clustered regularly interspaced short palindromic repeats and CRISPR-associated protein (CRISPR-Cas) systems store the memory of past encounters with foreign DNA in unique spacers that are inserted between direct repeats in CRISPR arrays. For only a small fraction of the spacers, homologous sequences, called protospacers, are detectable in viral, plasmid, and microbial genomes. The rest of the spacers remain the CRISPR “dark matter.” We performed a comprehensive analysis of the spacers from all CRISPR-cas loci identified in bacterial and archaeal genomes, and we found that, depending on the CRISPR-Cas subtype and the prokaryotic phylum, protospacers were detectable for 1% to about 19% of the spacers (~7% global average). Among the detected protospacers, the majority, typically 80 to 90%, originated from viral genomes, including proviruses, and among the rest, the most common source was genes that are integrated into microbial chromosomes but are involved in plasmid conjugation or replication. Thus, almost all spacers with identifiable protospacers target mobile genetic elements (MGE). The GC content, as well as dinucleotide and tetranucleotide compositions, of microbial genomes, their spacer complements, and the cognate viral genomes showed a nearly perfect correlation and were almost identical. Given the near absence of self-targeting spacers, these findings are most compatible with the possibility that the spacers, including the dark matter, are derived almost completely from the species-specific microbial mobilomes. IMPORTANCE The principal function of CRISPR-Cas systems is thought to be protection of bacteria and archaea against viruses and other parasitic genetic elements. The CRISPR defense function is mediated by sequences from parasitic elements, known as spacers, that are inserted into CRISPR arrays and then transcribed and employed as guides to identify and inactivate the cognate parasitic genomes. However, only a small fraction of the CRISPR spacers match any sequences in the current databases, and of these, only a minority correspond to known parasitic elements. We show that nearly all spacers with matches originate from viral or plasmid genomes that are either free or have been integrated into the host genome. We further demonstrate that spacers with no matches have the same properties as those of identifiable origins, strongly suggesting that all spacers originate from mobile elements.

ABSTRACT Cas1 integrase is the key enzyme of the clustered regularly interspaced short palindromic repeat (CRISPR)-Cas adaptation module that mediates acquisition of spacers derived from foreign DNA by CRISPR arrays. In diverse bacteria, the cas1 gene is fused (or adjacent) to a gene encoding a reverse transcriptase (RT) related to group II intron RTs. An RT-Cas1 fusion protein has been recently shown to enable acquisition of CRISPR spacers from RNA. Phylogenetic analysis of the CRISPR-associated RTs demonstrates monophyly of the RT-Cas1 fusion, and coevolution of the RT and Cas1 domains. Nearly all such RTs are present within type III CRISPR-Cas loci, but their phylogeny does not parallel the CRISPR-Cas type classification, indicating that RT-Cas1 is an autonomous functional module that is disseminated by horizontal gene transfer and can function with diverse type III systems. To compare the sequence pools sampled by RT-Cas1-associated and RT-lacking CRISPR-Cas systems, we obtained samples of a commercially grown cyanobacterium—Arthrospira platensis. Sequencing of the CRISPR arrays uncovered a highly diverse population of spacers. Spacer diversity was particularly striking for the RT-Cas1-containing type III-B system, where no saturation was evident even with millions of sequences analyzed. In contrast, analysis of the RT-lacking type III-D system yielded a highly diverse pool but reached a point where fewer novel spacers were recovered as sequencing depth was increased. Matches could be identified for a small fraction of the non-RT-Cas1-associated spacers, and for only a single RT-Cas1-associated spacer. Thus, the principal source(s) of the spacers, particularly the hypervariable spacer repertoire of the RT-associated arrays, remains unknown. IMPORTANCE While the majority of CRISPR-Cas immune systems adapt to foreign genetic elements by capturing segments of invasive DNA, some systems carry reverse transcriptases (RTs) that enable adaptation to RNA molecules. From analysis of available bacterial sequence data, we find evidence that RT-based RNA adaptation machinery has been able to join with CRISPR-Cas immune systems in many, diverse bacterial species. To investigate whether the abilities to adapt to DNA and RNA molecules are utilized for defense against distinct classes of invaders in nature, we sequenced CRISPR arrays from samples of commercial-scale open-air cultures of Arthrospira platensis, a cyanobacterium that contains both RT-lacking and RT-containing CRISPR-Cas systems. We uncovered a diverse pool of naturally occurring immune memories, with the RT-lacking locus acquiring a number of segments matching known viral or bacterial genes, while the RT-containing locus has acquired spacers from a distinct sequence pool for which the source remains enigmatic.

Clustering of functionally related genes in operons allows for coregulated gene expression in prokaryotes. This is advantageous when equal amounts of gene products are required. Production of protein complexes with an uneven stoichiometry, however, requires tuning mechanisms to generate subunits in appropriate relative quantities. Using comparative genomic analysis, we show that differential translation is a key determinant of modulated expression of genes clustered in operons and that codon bias generally is the best in silico indicator of unequal protein production. Variable ribosome density profiles of polycistronic transcripts correlate strongly with differential translation patterns. In addition, we provide experimental evidence that de novo initiation of translation can occur at intercistronic sites, allowing for differential translation of any gene irrespective of its position on a polycistronic messenger. Thus, modulation of translation efficiency appears to be a universal mode of control in bacteria and archaea that allows for differential production of operon-encoded proteins.

Over the past decade, the free-living flatworm Macrostomum lignano has been successfully used in many areas of biology, including embryology, stem cells, sexual selection, bioadhesion and aging. The increased use of this powerful laboratory model, including the establishment of genomic resources and tools, makes it essential to have a detailed description of the chromosome organization of this species, previously suggested to have a karyotype with 2n = 8 and one pair of large and three pairs of small metacentric chromosomes. We performed cytogenetic analyses for chromosomes of one commonly used inbred line of M. lignano (called DV1) and uncovered unexpected chromosome number variation in the form of aneuploidies of the largest chromosomes. These results prompted us to perform karyotypic studies in individual specimens of this and other lines of M. lignano reared under laboratory conditions, as well as in freshly field-collected specimens from different natural populations. Our analyses revealed a high frequency of aneuploids and in some cases other numerical and structural chromosome abnormalities in laboratory-reared lines of M. lignano, and some cases of aneuploidy were also found in freshly field-collected specimens. Moreover, karyological analyses were performed in specimens of three further species: Macrostomum sp. 8 (a close relative of M. lignano), M. spirale and M. hystrix. Macrostomum sp. 8 showed a karyotype that was similar to that of M. lignano, with tetrasomy for its largest chromosome being the most common karyotype, while the other two species showed a simpler karyotype that is more typical of the genus Macrostomum. These findings suggest that M. lignano and Macrostomum sp. 8 can be used as new models for studying processes of partial genome duplication in genome evolution.

ABSTRACT The Trojan horse Escherichia coli antibiotic microcin C (McC) consists of a heptapeptide attached to adenosine through a phosphoramidate linkage. McC is synthesized by the MccB enzyme, which terminally adenylates the ribosomally synthesized heptapeptide precursor MccA. The peptide part is responsible for McC uptake; it is degraded inside the cell to release a toxic nonhydrolyzable aspartyl-adenylate. Bionformatic analysis reveals that diverse bacterial genomes encoding mccB homologues also contain adjacent short open reading frames that may encode MccA-like adenylation substrates. Using chemically synthesized predicted peptide substrates and recombinant cognate MccB protein homologs, adenylated products were obtained in vitro for predicted MccA peptide-MccB enzyme pairs from Helicobacter pylori, Streptococcus thermophilus, Lactococcus johnsonii, Bartonella washoensis, Yersinia pseudotuberculosis, and Synechococcus sp. Some adenylated products were shown to inhibit the growth of E. coli by targeting aspartyl-tRNA synthetase, the target of McC. IMPORTANCE Our results prove that McC-like adenylated peptides are widespread and are encoded by both Gram-negative and Gram-positive bacteria and by cyanobacteria, opening ways for analyses of physiological functions of these compounds and for creation of microcin C-like antibiotics targeting various bacteria.

The elaborate eukaryotic DNA replication machinery evolved from the archaeal ancestors that themselves show considerable complexity. Here we discuss the comparative genomic and phylogenetic analysis of the core replication enzymes, the DNA polymerases, in archaea and their relationships with the eukaryotic polymerases. In archaea, there are three groups of family B DNA polymerases, historically known as PolB1, PolB2 and PolB3. All three groups appear to descend from the last common ancestors of the extant archaea but their subsequent evolutionary trajectories seem to have been widely different. Although PolB3 is present in all archaea, with the exception of Thaumarchaeota, and appears to be directly involved in lagging strand replication, the evolution of this gene does not follow the archaeal phylogeny, conceivably due to multiple horizontal transfers and/or dramatic differences in evolutionary rates. In contrast, PolB1 is missing in Euryarchaeota but otherwise seems to have evolved vertically. The third archaeal group of family B polymerases, PolB2, includes primarily proteins in which the catalytic centers of the polymerase and exonuclease domains are disrupted and accordingly the enzymes appear to be inactivated. The members of the PolB2 group are scattered across archaea and might be involved in repair or regulation of replication along with inactivated members of the RadA family ATPases and an additional, uncharacterized protein that are encoded within the same predicted operon. In addition to the family B polymerases, all archaea, with the exception of the Crenarchaeota, encode enzymes of a distinct family D the origin of which is unclear. We examine multiple considerations that appear compatible with the possibility that family D polymerases are highly derived homologs of family B. The eukaryotic DNA polymerases show a highly complex relationship with their archaeal ancestors including contributions of proteins and domains from both the family B and the fam

Many neurodegenerative diseases, such as Alzheimer’s disease, Parkinson disease, vascular and frontotemporal dementias, as well as other chronic neurological pathologies, are characterized by slow development with a long asymptomatic period followed by a stage with mild clinical symptoms. As a consequence, these serious pathologies are diagnosed late in the course of a disease, when massive death of neurons has already occurred and effective therapeutic intervention is problematic. Thus, the development of screening tests capable of detecting neurodegenerative diseases during early, preferably asymptomatic, stages is a high unmet need. Since such tests are to be used for screening of large populations, they should be non-invasive and relatively inexpensive. Further, while subjects identified by screening tests can be further tested with more invasive and expensive methods, e.g. analysis of cerebrospinal fluid or imaging techniques, to be of practical utility screening tests should have high sensitivity and specificity. In this review, we discuss advantages and disadvantages of various approaches to developing screening tests based on analysis of circulating cell-free miRNA. Applications of circulating miRNA-based tests for diagnosis of acute and chronic brain pathologies, for research of normal brain aging, and for disease and treatment monitoring are also discussed.