Search

The impact of genetic regulatory variation active in early pancreatic development on adult pancreatic disease and traits is not well understood. Here, we generate a panel of 107 fetal-like iPSC-derived pancreatic progenitor cells (iPSC-PPCs) from whole genome-sequenced individuals and identify 4065 genes and 4016 isoforms whose expression and/or alternative splicing are affected by regulatory variation. We integrate eQTLs identified in adult islets and whole pancreas samples, which reveal 1805 eQTL associations that are unique to the fetal-like iPSC-PPCs and 1043 eQTLs that exhibit regulatory plasticity across the fetal-like and adult pancreas tissues. Colocalization with GWAS risk loci for pancreatic diseases and traits show that some putative causal regulatory variants are active only in the fetal-like iPSC-PPCs and likely influence disease by modulating expression of disease-associated genes in early development, while others with regulatory plasticity likely exert their effects in both the fetal and adult pancreas by modulating expression of different disease genes in the two developmental stages.

MOTIVATION: Double-stranded RNAs (dsRNAs) are potent triggers of innate immune responses upon recognition by cytosolic dsRNA sensor proteins. Identification of endogenous dsRNAs helps to better understand the dsRNAome and its relevance to innate immunity related to human diseases. RESULTS: Here, we report dsRID (double-stranded RNA identifier), a machine-learning-based method to predict dsRNA regions in silico, leveraging the power of long-read RNA-sequencing (RNA-seq) and molecular traits of dsRNAs. Using models trained with PacBio long-read RNA-seq data derived from Alzheimers disease (AD) brain, we show that our approach is highly accurate in predicting dsRNA regions in multiple datasets. Applied to an AD cohort sequenced by the ENCODE consortium, we characterize the global dsRNA profile with potentially distinct expression patterns between AD and controls. Together, we show that dsRID provides an effective approach to capture global dsRNA profiles using long-read RNA-seq data. AVAILABILITY AND IMPLEMENTATION: Software implementation of dsRID, and genomic coordinates of regions predicted by dsRID in all samples are available at the GitHub repository: https://github.com/gxiaolab/dsRID.

Different plant species within the grasses were parallel targets of domestication, giving rise to crops with distinct evolutionary histories and traits1. Key traits that distinguish these species are mediated by specialized cell types2. Here we compare the transcriptomes of root cells in three grass species-Zea mays, Sorghum bicolor and Setaria viridis. We show that single-cell and single-nucleus RNA sequencing provide complementary readouts of cell identity in dicots and monocots, warranting a combined analysis. Cell types were mapped across species to identify robust, orthologous marker genes. The comparative cellular analysis shows that the transcriptomes of some cell types diverged more rapidly than those of others-driven, in part, by recruitment of gene modules from other cell types. The data also show that a recent whole-genome duplication provides a rich source of new, highly localized gene expression domains that favour fast-evolving cell types. Together, the cell-by-cell comparative analysis shows how fine-scale cellular profiling can extract conserved modules from a pan transcriptome and provide insight on the evolution of cells that mediate key functions in crops.

Poa pratensis, commonly known as Kentucky bluegrass, is a popular cool-season grass species used as turf in lawns and recreation areas globally. Despite its substantial economic value, a reference genome had not previously been assembled due to the genome's relatively large size and biological complexity that includes apomixis, polyploidy, and interspecific hybridization. We report here a fortuitous de novo assembly and annotation of a P. pratensis genome. Instead of sequencing the genome of a C4 grass, we accidentally sampled and sequenced tissue from a weedy P. pratensis whose stolon was intertwined with that of the C4 grass. The draft assembly consists of 6.09 Gbp with an N50 scaffold length of 65.1 Mbp, and a total of 118 scaffolds, generated using PacBio long reads and Bionano optical map technology. We annotated 256K gene models and found 58% of the genome to be composed of transposable elements. To demonstrate the applicability of the reference genome, we evaluated population structure and estimated genetic diversity in P. pratensis collected from three North American prairies, two in Manitoba, Canada and one in Colorado, USA. Our results support previous studies that found high genetic diversity and population structure within the species. The reference genome and annotation will be an important resource for turfgrass breeding and study of bluegrasses.

The Stellers jay is a familiar bird of western forests from Alaska south to Nicaragua. Here, we report a draft reference assembly for the species generated from PacBio HiFi long-read and Omni-C chromatin-proximity sequencing data as part of the California Conservation Genomics Project (CCGP). Sequenced reads were assembled into 352 scaffolds totaling 1.16 Gb in length. Assembly metrics indicate a highly contiguous and complete assembly with a contig N50 of 7.8 Mb, scaffold N50 of 25.8 Mb, and BUSCO completeness score of 97.2%. Repetitive elements span 16.6% of the genome including nearly 90% of the W chromosome. Compared with high-quality assemblies from other members of the family Corvidae, the Stellers jay genome contains a larger proportion of repetitive elements than 4 crow species (Corvus), but a lower proportion of repetitive elements than the California scrub-jay (Aphelocoma californica). This reference genome will serve as an essential resource for future studies on speciation, local adaptation, phylogeography, and conservation genetics in this species of significant biological interest.

Spray-induced gene silencing (SIGS) is an emerging tool for crop pest protection. It utilizes exogenously applied double-stranded RNA to specifically reduce pest target gene expression using endogenous RNA interference machinery. In this study, SIGS methods were developed and optimized for powdery mildew fungi, which are widespread obligate biotrophic fungi that infect agricultural crops, using the known azole-fungicide target cytochrome P450 51 (CYP51) in the Golovinomyces orontii-Arabidopsis thaliana pathosystem. Additional screening resulted in the identification of conserved gene targets and processes important to powdery mildew proliferation: apoptosis-antagonizing transcription factor in essential cellular metabolism and stress response; lipid catabolism genes lipase a, lipase 1, and acetyl-CoA oxidase in energy production; and genes involved in manipulation of the plant host via abscisic acid metabolism (9-cis-epoxycarotenoid dioxygenase, xanthoxin dehydrogenase, and a putative abscisic acid G-protein coupled receptor) and secretion of the effector protein, effector candidate 2. Powdery mildew is the dominant disease impacting grapes and extensive powdery mildew resistance to applied fungicides has been reported. We therefore developed SIGS for the Erysiphe necator-Vitis vinifera system and tested six successful targets identified using the G. orontii-A. thaliana system. For all targets tested, a similar reduction in powdery mildew disease was observed between systems. This indicates screening of broadly conserved targets in the G. orontii-A. thaliana pathosystem identifies targets and processes for the successful control of other powdery mildew fungi. The efficacy of SIGS on powdery mildew fungi makes SIGS an exciting prospect for commercial powdery mildew control.

Wine cultivars are available to growers in multiple clonal selections with agronomic and enological differences. Phenotypic differences between clones originated from somatic mutations that accrued over thousands of asexual propagation cycles. Genetic diversity between grape cultivars remains unexplored, and tools to discriminate unequivocally clones have been lacking. This study aimed to uncover genetic variations among a group of clonal selections of 4 important Vitis vinifera cultivars: Cabernet sauvignon, Sauvignon blanc, Chardonnay, and Merlot, and use this information to develop genetic markers to discriminate the clones of these cultivars. We sequenced with short-read sequencing technology the genomes of 18 clones, including biological replicates for a total of 46 genomes. Sequences were aligned to their respective cultivars reference genome for variant calling. We used reference genomes of Cabernet sauvignon, Chardonnay, and Merlot and developed a de novo genome assembly of Sauvignon blanc using long-read sequencing. On average, 4 million variants were detected for each clone, with 74.2% being single nucleotide variants and 25.8% being small insertions or deletions (InDel). The frequency of these variants was consistent across all clones. From these variants, we validated 46 clonal markers using high-throughput amplicon sequencing for 77.7% of the evaluated clones, most of them small InDel. These results represent an advance in grapevine genotyping strategies and will benefit the viticulture industry for the characterization and identification of the plant material.

Violacein is a water-insoluble violet pigment produced by various Gram-negative bacteria. The compound and the bacteria that produce it have been gaining attention due to the antimicrobial and proposed antitumour properties of violacein and the possibility that strains producing it may have broad industrial uses. Bacteria that produce violacein have been isolated from diverse environments including fresh and ocean waters, glaciers, tropical soils, trees, fish and the skin of amphibians. We report here the isolation and characterization of six violacein-producing bacterial strains and three non-violacein-producing close relatives, each isolated from either an aquatic environment or moist food materials in northern California, USA. For each isolate, we characterized traditional phenotypes, generated and analysed draft genome sequences, and carried out multiple types of taxonomic, phylogenetic and phylogenomic analyses. Based on these analyses we assign putative identifications to the nine isolates, which include representatives of the genera Chromobacterium, Aquitalea, Iodobacter, Duganella, Massilia and Janthinobacterium. In addition, we discuss the utility of various metrics for taxonomic assignment in these groups including average nucleotide identity, whole genome phylogenetic analysis and extent of recent homologous recombination using the software program PopCOGenT.

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

PurposeThere is limited payer coverage for genome sequencing (GS) relative to exome sequencing (ES) in the U.S. Our objective was to assess payers' considerations for coverage of GS versus coverage of ES and requirements payers have for coverage of GS. The study was conducted by the NIH-funded Clinical Sequencing Evidence-Generating Research Consortium (CSER).MethodsWe conducted semi-structured interviews with representatives of private payer organizations (payers, N = 12) on considerations and evidentiary and other needs for coverage of GS and ES. Data were analyzed using thematic analysis.ResultsWe described four categories of findings and solutions: demonstrated merits of GS versus ES, enhanced methods for evidence generation, consistent laboratory processes/sequencing methods, and enhanced implementation/care delivery. Payers see advantages to GS vs. ES and are open to broader GS coverage but need more proof of these advantages to consider them in coverage decision-making. Next steps include establishing evidence of benefits in specific clinical scenarios, developing quality standards, ensuring transparency of laboratory methods, developing clinical centers of excellence, and incorporating the role of genetic professionals.ConclusionBy comparing coverage considerations for GS and ES, we identified a path forward for coverage of GS. Future research should explicitly address payers' conditions for coverage.

A global international initiative, such as the Earth BioGenome Project (EBP), requires both agreement and coordination on standards to ensure that the collective effort generates rapid progress toward its goals. To this end, the EBP initiated five technical standards committees comprising volunteer members from the global genomics scientific community: Sample Collection and Processing, Sequencing and Assembly, Annotation, Analysis, and IT and Informatics. The current versions of the resulting standards documents are available on the EBP website, with the recognition that opportunities, technologies, and challenges may improve or change in the future, requiring flexibility for the EBP to meet its goals. Here, we describe some highlights from the proposed standards, and areas where additional challenges will need to be met.

BackgroundSeahorses, seadragons, pygmy pipehorses, and pipefishes (Syngnathidae, Syngnathiformes) are among the most recognizable groups of fishes because of their derived morphology, unusual life history, and worldwide distribution. Despite previous phylogenetic studies and recent new species descriptions of syngnathids, the evolutionary relationships among several major groups within this family remain unresolved.ResultsHere, we provide a reconstruction of syngnathid phylogeny based on genome-wide sampling of 1314 ultraconserved elements (UCEs) and expanded taxon sampling to assess the current taxonomy and as a basis for macroevolutionary insights. We sequenced a total of 244 new specimens across 117 species and combined with published UCE data for a total of 183 species of Syngnathidae, about 62% of the described species diversity, to compile the most data-rich phylogeny to date. We estimated divergence times using 14 syngnathiform fossils, including nine fossils with newly proposed phylogenetic affinities, to better characterize current and historical biogeographical patterns, and to reconstruct diversification through time. We present a phylogenetic hypothesis that is well-supported and provides several notable insights into syngnathid evolution. We found nine non-monophyletic genera, evidence for seven cryptic species, five potentially invalid synonyms, and identified a novel sister group to the seahorses, the Indo-Pacific pipefishes Halicampus macrorhynchus and H. punctatus. In addition, the morphologically distinct southwest Pacific seahorse Hippocampus jugumus was recovered as the sister to all other non-pygmy seahorses. As found in many other groups, a high proportion of syngnathid lineages appear to have originated in the Central Indo-Pacific and subsequently dispersed to adjoining regions. Conversely, we also found an unusually high subsequent return of lineages from southern Australasia to the Central Indo-Pacific. Diversification rates rose abrupt

Tobamoviruses are agriculturally relevant viruses that cause crop losses and have infected plants in many regions of the world. These viruses are frequently found in municipal wastewater, likely coming from human diet and industrial waste across wastewater catchment areas. As part of a large wastewater-based epidemiology study across Southern California, we analyzed RNA sequence data from 275 influent wastewater samples obtained from eight wastewater treatment plants with a catchment area of approximately 16 million people from July 2020 to August 2021. We assembled 1,083 high-quality genomes, enumerated viral sequencing reads, and detected thousands of single nucleotide variants from eight common tobamoviruses: bell pepper mottle virus, cucumber green mottle mosaic virus, pepper mild mottle virus, tobacco mild green mosaic virus, tomato brown rugose fruit virus, tomato mosaic virus, tomato mottle mosaic virus, and tropical soda apple mosaic virus. We show that single nucleotide variants had amino acid-altering consequences along with synonymous mutations, which represents potential evolution with functional consequences in genomes of these viruses. Our study shows the importance of wastewater sequencing to monitor the genomic diversity of these plant-infecting viruses, and we suggest that our data could be used to continue tracking the genomic variability of such pathogens. IMPORTANCE Diseases caused by viruses in the genus Tobamovirus cause crop losses around the world. As with other viruses, mutation occurring in the virus's genomes can have functional consequences and may alter viral infectivity. Many of these plant-infecting viruses have been found in wastewater, likely coming from human consumption of infected plants and produce. By sequencing RNA extracted from influent wastewater obtained from eight wastewater treatment plants in Southern California, we assembled high-quality viral genomes and detected thousands of single nucleotide variants from eight tobam

Human herpes simplex virus 2 (HSV-2) is a ubiquitous, slowly evolving DNA virus. HSV-2 has two primary lineages, one found in West and Central Africa and the other found worldwide. Competing hypotheses have been proposed to explain how HSV-2 migrated out-of-Africa (i)HSV-2 followed human migration out-of-Africa 50-100 thousand years ago, or (ii)HSV-2 migrated via the trans-Atlantic slave trade 150-500 years ago. Limited geographic sampling and lack of molecular clock signal has precluded robust comparison. Here, we analyze newly sequenced HSV-2 genomes from Africa to resolve geography and timing of divergence events within HSV-2. Phylogeographic analysis consistently places the ancestor of worldwide dispersal in East Africa, though molecular clock is too slow to be detected using available data. Rates 4.2 × 10-8-5.6 × 10-8 substitutions/site/year, consistent with previous age estimates, suggest a worldwide dispersal 22-29 thousand years ago. Thus, HSV-2 likely migrated with humans from East Africa and dispersed after the Last Glacial Maximum.

Evolutionary biologists have long sought to disentangle phylogenetic relationships among taxa spanning the tree of life, an increasingly important task as anthropogenic influences accelerate population declines and species extinctions, particularly in insects. Phylogenetic analyses are commonly used to identify unique evolutionary lineages, to clarify taxonomic designations of the focal taxa, and to inform conservation decisions. Advances in DNA sequencing techniques have increasingly facilitated the ability of researchers to apply genomic methods to phylogenetic analyses, even for non-model organisms. Stoneflies are non-model insects that are important bioindicators of the quality of freshwater habitats and landscape disturbance as they spend the immature stages of their life cycles in fresh water, and the adult stages in terrestrial environments. Phylogenetic relationships within the stonefly genus Suwallia (Insecta: Plecoptera: Chloroperlidae) are poorly understood, and have never been assessed using molecular data. We used DNA sequence data from genome-wide ultraconserved element loci to generate the first molecular phylogeny for the group and assess its monophyly. We found that Palearctic and Nearctic Suwallia do not form reciprocally monophyletic clades, and that a biogeographic history including dispersal, vicariance, and founder event speciation via jump dispersal best explains the geographic distribution of this group. Our results also strongly suggest that Neaviperla forcipata (Neave, 1929) is nested within Suwallia, and the concept of the genus Suwallia should be revised to include it. Thus, we formally propose a new taxonomic combination wherein Neaviperla forcipata (Neave, 1929) is reclassified as Suwallia forcipata (Neave, 1929). Moreover, some Suwallia species (e.g., S. amoenacolens, S. kerzhneri, S. marginata, S. pallidula, and S. starki) exhibit pronounced cryptic diversity that is worthy of further investigation. These findings provide a first glim

Subviral agents are nucleic acids which lack the features for classification as a virus. Tombusvirus-like associated RNAs (tlaRNAs) are subviral positive-sense, single-stranded RNAs that replicate autonomously, yet depend on a coinfecting virus for encapsidation and transmission. TlaRNAs produce abundant subgenomic RNA (sgRNA) upon infection. Here, we investigate how the well-studied tlaRNA, ST9, produces sgRNA and its function. We found ST9 is a noncoding RNA, due to its lack of protein coding capacity. We used resistance assays with eukaryotic Exoribonuclease-1 (XRN1) to investigate sgRNA production via incomplete degradation of genomic RNA. The ST9 3' untranslated region stalled XRN1 very near the 5' sgRNA end. Thus, the XRN family of enzymes drives sgRNA accumulation in ST9-infected tissue by incomplete degradation of ST9 RNA. This work suggests tlaRNAs are not just parasites of viruses with compatible capsids, but also mutually beneficial partners that influence host cell RNA biology.

The mitochondrial genome (mtDNA) is an important source of disease-causing genetic variability, but existing sequencing methods limit understanding, precluding phased measurement of mutations and clear detection of large sporadic deletions. We adapted a method for amplification-free sequence enrichment using Cas9 cleavage to obtain full length nanopore reads of mtDNA. We then utilized the long reads to phase mutations in a patient with an mtDNA-linked syndrome and demonstrated that this method can map age-induced mtDNA deletions. We believe this method will offer deeper insight into our understanding of mtDNA variation.

Eight new species of Isodiametridae (Acoela) are presented and phylogenetic relationships of the family using 18S and COI sequence data reconstructed. Four species of Baltalimania, two species of Praeaphanostoma and two species of Pseudoposthia are described that all can be distinguished both morphologically and through differences in the DNA sequences as supported by mPTP analysis. The addition of eight new species increases the number of nominal species of Isodiametridae worldwide by 10% and within Sweden by 31%, but also highlights several taxonomic problems. Several species are now confirmed as isodiametrids that are otherwise incongruent with the family's diagnostic characters. Generic keys are presented for Baltalimania and Praeaphanostoma.

Nascent RNA-sequencing tracks transcription at nucleotide resolution. The genomic distribution of engaged transcription complexes, in turn, uncovers functional genomic regions. Here, we provide analytical steps to (1) identify transcribed regulatory elements de novo genome-wide, (2) quantify engaged transcription complexes at enhancers, promoter-proximal regions, divergent transcripts, gene bodies, and termination windows, and (3) measure distribution of transcription machineries and regulatory proteins across functional genomic regions. This protocol tracks engaged transcription complexes across functional genomic regions demonstrated in human K562 erythroleukemia cells. For complete details on the use and execution of this protocol, please refer to Vihervaara et al. (2021)., QC 20220930

Structural variants (SVs) are implicated in the etiology of Mendelian diseases but have been systematically underascertained owing to sequencing technology limitations. Long-read sequencing enables comprehensive detection of SVs, but approaches for prioritization of candidate SVs are needed. Structural variant Annotation and analysis (SvAnna) assesses all classes of SVs and their intersection with transcripts and regulatory sequences, relating predicted effects on gene function with clinical phenotype data. SvAnna places 87% of deleterious SVs in the top ten ranks. The interpretable prioritizations offered by SvAnna will facilitate the widespread adoption of long-read sequencing in diagnostic genomics. SvAnna is available at https://github.com/TheJacksonLaboratory/SvAnn a .

Regulation of bacteriophage gene expression involves repressor proteins that bind and downregulate early lytic promoters. A large group of mycobacteriophages code for repressors that are unusual in also terminating transcription elongation at numerous binding sites (stoperators) distributed across the phage genome. Here we provide the X-ray crystal structure of a mycobacteriophage immunity repressor bound to DNA, which reveals the binding of a monomer to an asymmetric DNA sequence using two independent DNA binding domains. The structure is supported by small-angle X-ray scattering, DNA binding, molecular dynamics, and in vivo immunity assays. We propose a model for how dual DNA binding domains facilitate regulation of both transcription initiation and elongation, while enabling evolution of other superinfection immune specificities.

BackgroundThe eukaryotic parasite Plasmodium falciparum causes millions of malarial infections annually while drug resistance to common anti-malarials is further confounding eradication efforts. Translation is an attractive therapeutic target that will benefit from a deeper mechanistic understanding. As the rate limiting step of translation, initiation is a primary driver of translational efficiency. It is a complex process regulated by both cis and trans acting factors, providing numerous potential targets. Relative to model organisms and humans, P. falciparum mRNAs feature unusual 5' untranslated regions suggesting cis-acting sequence complexity in this parasite may act to tune levels of protein synthesis through their effects on translational efficiency.MethodsHere, in vitro translation is deployed to compare the role of cis-acting regulatory sequences in P. falciparum and humans. Using parasite mRNAs with high or low translational efficiency, the presence, position, and termination status of upstream "AUG"s, in addition to the base composition of the 5' untranslated regions, were characterized.ResultsThe density of upstream "AUG"s differed significantly among the most and least efficiently translated genes in P. falciparum, as did the average "GC" content of the 5' untranslated regions. Using exemplars from highly translated and poorly translated mRNAs, multiple putative upstream elements were interrogated for impact on translational efficiency. Upstream "AUG"s were found to repress translation to varying degrees, depending on their position and context, while combinations of upstream "AUG"s had non-additive effects. The base composition of the 5' untranslated regions also impacted translation, but to a lesser degree. Surprisingly, the effects of cis-acting sequences were remarkably conserved between P. falciparum and humans.ConclusionsWhile translational regulation is inherently complex, this work contributes toward a more comprehensive unde

The major biological methyl donor, S-adenosylmethionine (adoMet) synthesis occurs mainly in the liver. Methionine adenosyltransferase 1A (MAT1A) and glycine N-methyltransferase (GNMT) are two key enzymes involved in the functional implications of that variation. We collected 42 RNA-seq data from paired hepatocellular carcinoma (HCC) and its adjacent normal liver tissue from the Cancer Genome Atlas (TCGA). There was no mutation found in MAT1A or GNMT RNA in the 42 HCC patients. The 11,799 genes were annotated in the RNA-Seq data, and their expression levels were used to investigate the phenotypes of low MAT1A and low GNMT by Gene Set Enrichment Analysis (GSEA). The REACTOME_TRANSLATION gene set was enriched and visualized in a heatmap along with corresponding differences in gene expression between low MAT1A versus high MAT1A and low GNMT versus high GNMT. We identified 43 genes of the REACTOME_TRANSLATION gene set that are powerful prognosis factors in HCC. The significantly predicted genes were referred into eukaryotic translation initiation (EIF3B, EIF3K), eukaryotic translation elongation (EEF1D), and ribosomal proteins (RPs). Cell models expressing various MAT1A and GNMT proved that simultaneous restoring the expression of MAT1A and GNMT decreased cell proliferation, invasion, as well as the REACTOME_TRANSLATION gene EEF1D, consistent with a better prognosis in human HCC. We demonstrated new findings that downregulation or defect in MAT1A and GNMT genes can enrich the protein-associated translation process that may account for poor HCC prognosis. This is the first study demonstrated that MAT1A and GNMT, the 2 key enzymes involved in methionine cycle, could attenuate the function of ribosome translation. We propose a potential novel mechanism by which the diminished GNMT and MAT1A expression may confer poor prognosis for HCC.

HIV-1-infected cells that persist despite antiretroviral therapy (ART) are frequently considered "transcriptionally silent," but active viral gene expression may occur in some cells, challenging the concept of viral latency. Applying an assay for profiling the transcriptional activity and the chromosomal locations of individual proviruses, we describe a global genomic and epigenetic map of transcriptionally active and silent proviral species and evaluate their longitudinal evolution in persons receiving suppressive ART. Using genome-wide epigenetic reference data, we show that proviral transcriptional activity is associated with activating epigenetic chromatin features in linear proximity of integration sites and in their inter- and intrachromosomal contact regions. Transcriptionally active proviruses were actively selected against during prolonged ART; however, this pattern was violated by large clones of virally infected cells that may outcompete negative selection forces through elevated intrinsic proliferative activity. Our results suggest that transcriptionally active proviruses are dynamically evolving under selection pressure by host factors.

Small RNAs (sRNAs) are known to regulate pathogenic plant-microbe interactions. Emerging evidence from the study of these model systems suggests that microRNAs (miRNAs) can be translocated between microbes and plants to facilitate symbiosis. The roles of sRNAs in mutualistic mycorrhizal fungal interactions, however, are largely unknown. In this study, we characterized miRNAs encoded by the ectomycorrhizal fungus Pisolithus microcarpus and investigated their expression during mutualistic interaction with Eucalyptus grandis. Using sRNA sequencing data and in situ miRNA detection, a novel fungal miRNA, Pmic_miR-8, was found to be transported into E. grandis roots after interaction with P. microcarpus Further characterization experiments demonstrate that inhibition of Pmic_miR-8 negatively impacts the maintenance of mycorrhizal roots in E. grandis, while supplementation of Pmic_miR-8 led to deeper integration of the fungus into plant tissues. Target prediction and experimental testing suggest that Pmic_miR-8 may target the host NB-ARC domain containing transcripts, suggesting a potential role for this miRNA in subverting host signaling to stabilize the symbiotic interaction. Altogether, we provide evidence of previously undescribed cross-kingdom sRNA transfer from ectomycorrhizal fungi to plant roots, shedding light onto the involvement of miRNAs during the developmental process of mutualistic symbioses.

Coccidioidomycosis is a common fungal disease that is endemic to arid and semi-arid regions of both American continents. Coccidioides immitis and Coccidioides posadasii are the etiological agents of the disease, also known as Valley Fever. For several decades, the C. posadasii strain Silveira has been used widely in vaccine studies, is the source strain for production of diagnostic antigens, and is a widely used experimental strain for functional studies. In 2009, the genome was sequenced using Sanger sequencing technology, and a draft assembly and annotation were made available. In this study, the genome of the Silveira strain was sequenced using single molecule real-time sequencing PacBio technology, assembled into chromosomal-level contigs, genotyped, and the genome was reannotated using sophisticated and curated in silico tools. This high-quality genome sequencing effort has improved our understanding of chromosomal structure, gene set annotation, and lays the groundwork for identification of structural variants (e.g. transversions, translocations, and copy number variants), assessment of gene gain and loss, and comparison of transposable elements in future phylogenetic and population genomics studies.

The cellular response to stress is an important determinant of disease pathogenesis. Uncovering the molecular fingerprints of distinct stress responses may identify novel biomarkers and key signaling pathways for different diseases. Emerging evidence shows that transfer RNA-derived small RNAs (tDRs) play pivotal roles in stress responses. However, RNA modifications present on tDRs are barriers to accurately quantifying tDRs using traditional small RNA sequencing. Here, AlkB-facilitated methylation sequencing is used to generate a comprehensive landscape of cellular and extracellular tDR abundances in various cell types during different stress responses. Extracellular tDRs are found to have distinct fragmentation signatures from intracellular tDRs and these tDR signatures are better indicators of different stress responses than miRNAs. These distinct extracellular tDR fragmentation patterns and signatures are also observed in plasma from patients on cardiopulmonary bypass. It is additionally demonstrated that angiogenin and RNASE1 are themselves regulated by stressors and contribute to the stress-modulated abundance of sub-populations of cellular and extracellular tDRs. Finally, a sub-population of extracellular tDRs is identified for which AGO2 appears to be required for their expression. Together, these findings provide a detailed profile of stress-responsive tDRs and provide insight about tDR biogenesis and stability in response to cellular stressors.

Over the last 25 years, biology has entered the genomic era and is becoming a science of 'big data'. Most interpretations of genomic analyses rely on accurate functional annotations of the proteins encoded by more than 500 000 genomes sequenced to date. By different estimates, only half the predicted sequenced proteins carry an accurate functional annotation, and this percentage varies drastically between different organismal lineages. Such a large gap in knowledge hampers all aspects of biological enterprise and, thereby, is standing in the way of genomic biology reaching its full potential. A brainstorming meeting to address this issue funded by the National Science Foundation was held during 3-4 February 2022. Bringing together data scientists, biocurators, computational biologists and experimentalists within the same venue allowed for a comprehensive assessment of the current state of functional annotations of protein families. Further, major issues that were obstructing the field were identified and discussed, which ultimately allowed for the proposal of solutions on how to move forward.

The giant ciliate Stentor coeruleus is a classical model system for studying regeneration and morphogenesis in a single cell. The anterior of the cell is marked by an array of cilia, known as the oral apparatus, which can be induced to shed and regenerate in a series of reproducible morphological steps, previously shown to require transcription. If a cell is cut in half, each half regenerates an intact cell. We used RNA sequencing (RNAseq) to assay the dynamic changes in Stentor's transcriptome during regeneration, after both oral apparatus shedding and bisection, allowing us to identify distinct temporal waves of gene expression including kinases, RNA -binding proteins, centriole biogenesis factors, and orthologs of human ciliopathy genes. By comparing transcriptional profiles of different regeneration events, we identified distinct modules of gene expression corresponding to oral apparatus regeneration, posterior holdfast regeneration, and recovery after wounding. By measuring gene expression after blocking translation, we show that the sequential waves of gene expression involve a cascade mechanism in which later waves of expression are triggered by translation products of early-expressed genes. Among the early-expressed genes, we identified an E2F transcription factor and the RNA-binding protein Pumilio as potential regulators of regeneration based on the expression pattern of their predicted target genes. RNAi-mediated knockdown experiments indicate that Pumilio is required for regenerating oral structures of the correct size. E2F is involved in the completion of regeneration but is dispensable for earlier steps. This work allows us to classify regeneration genes into groups based on their potential role for regeneration in distinct cell regeneration paradigms, and provides insight into how a single cell can coordinate complex morphogenetic pathways to regenerate missing structures.

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Structural variants (SVs) are implicated in the etiology of Mendelian diseases but have been systematically underascertained owing to sequencing technology limitations. Long-read sequencing enables comprehensive detection of SVs, but approaches for prioritization of candidate SVs are needed. Structural variant Annotation and analysis (SvAnna) assesses all classes of SVs and their intersection with transcripts and regulatory sequences, relating predicted effects on gene function with clinical phenotype data. SvAnna places 87% of deleterious SVs in the top ten ranks. The interpretable prioritizations offered by SvAnna will facilitate the widespread adoption of long-read sequencing in diagnostic genomics. SvAnna is available at https://github.com/TheJacksonLaboratory/SvAnn a .

Regulation of bacteriophage gene expression involves repressor proteins that bind and downregulate early lytic promoters. A large group of mycobacteriophages code for repressors that are unusual in also terminating transcription elongation at numerous binding sites (stoperators) distributed across the phage genome. Here we provide the X-ray crystal structure of a mycobacteriophage immunity repressor bound to DNA, which reveals the binding of a monomer to an asymmetric DNA sequence using two independent DNA binding domains. The structure is supported by small-angle X-ray scattering, DNA binding, molecular dynamics, and in vivo immunity assays. We propose a model for how dual DNA binding domains facilitate regulation of both transcription initiation and elongation, while enabling evolution of other superinfection immune specificities.

Coccidioidomycosis is a common fungal disease that is endemic to arid and semi-arid regions of both American continents. Coccidioides immitis and Coccidioides posadasii are the etiological agents of the disease, also known as Valley Fever. For several decades, the C. posadasii strain Silveira has been used widely in vaccine studies, is the source strain for production of diagnostic antigens, and is a widely used experimental strain for functional studies. In 2009, the genome was sequenced using Sanger sequencing technology, and a draft assembly and annotation were made available. In this study, the genome of the Silveira strain was sequenced using single molecule real-time sequencing PacBio technology, assembled into chromosomal-level contigs, genotyped, and the genome was reannotated using sophisticated and curated in silico tools. This high-quality genome sequencing effort has improved our understanding of chromosomal structure, gene set annotation, and lays the groundwork for identification of structural variants (e.g. transversions, translocations, and copy number variants), assessment of gene gain and loss, and comparison of transposable elements in future phylogenetic and population genomics studies.

The advent of single-molecule probing techniques has revolutionized the biomedical and life science fields and has spurred the development of a new class of labs-on-chip based on powerful biosensors. Nanopores represent one of the most recent and most promising single molecule sensing paradigms that is seeing increased chip-scale integration for improved convenience and performance. Due to their physical structure, nanopores are highly sensitive, require low sample volume, and offer label-free, amplification-free, high-throughput real-time detection and identification of biomolecules. Over the last 25 years, nanopores have been extensively employed to detect a variety of biomolecules with a growing range of applicatons ranging from nucleic acid sequencing to ultrasensitive diagnostics to single-molecule biophysics. Nanopores, in particular those in solid-state membranes, also have the potential for integration with other technologies such as optics, plasmonics, microfluidics, and optofluidics to perform more complex tasks for an ever-expanding demand. A number of breakthrough results using integrated nanopore platforms have already been reported, and more can be expected as nanopores remain the focus of innovative research and are finding their way into commercial instruments. This review provides an overview of different aspects and challenges of nanopore technology with a focus on chip-scale integration of solid-state nanopores for biosensing and bioanalytical applications.

RNA guanine quadruplexes (rG4) assume important roles in post-transcriptional regulations of gene expression, which are often modulated by rG4-binding proteins. Hence, understanding the biological functions of rG4s requires the identification and functional characterizations of rG4-recognition proteins. By employing a bioinformatic approach based on the analysis of overlap between peaks obtained from rG4-seq analysis and those detected in >230 eCLIP-seq datasets for RNA-binding proteins generated from the ENCODE project, we identified a large number of candidate rG4-binding proteins. We showed that one of these proteins, G3BP1, is able to bind directly to rG4 structures with high affinity and selectivity, where the binding entails its C-terminal RGG domain and is further enhanced by its RRM domain. Additionally, our seCLIP-Seq data revealed that pyridostatin, a small-molecule rG4 ligand, could displace G3BP1 from mRNA in cells, with the most pronounced effects being observed for the 3-untranslated regions (3-UTR) of mRNAs. Moreover, luciferase reporter assay results showed that G3BP1 positively regulates mRNA stability through its binding with rG4 structures. Together, we identified a number of candidate rG4-binding proteins and validated that G3BP1 can bind directly with rG4 structures and regulate the stabilities of mRNAs.

BACKGROUND: Duck viral enteritis, commonly known as duck plague (DP), is an acute and contagious fatal disease in ducks, geese, and swans caused by the DP virus (DPV). It poses a serious threat to the growth of duck farming in the Haor (wetland) areas of Bangladesh. AIM: This study aimed to detect the circulating DPV by molecular characterization, followed by phylogenetic analysis, targeting the UL30 gene in infected ducks from five Haor districts in Bangladesh and to observe the variation in the genome sequence between the field virus and vaccine strain of DPV. METHODS: A total of 150 samples (liver, 50; intestine, 50; and oropharyngeal tissue, 50) were collected from DP-suspected sick/dead ducks from 50 affected farms in Kishoreganj, Netrokona, B. Baria, Habiganj, and Sunamganj districts in Bangladesh. For the identification of DPV in collected samples, polymerase chain reaction (PCR) was utilized. Nucleotide sequences of the amplified UL30 gene were compared with those of other DPV strains available in GenBank. RESULTS: Of the 150 samples, 90 (60%) were found to be positive for DPV, as confirmed by PCR. Organ-wise prevalence was higher in the liver (72%), followed by the intestine (64%) and oropharyngeal tissue (44%). Regarding areas, the highest and lowest prevalence in the liver and intestine was observed in Habiganj and B. Baria, respectively, whereas the highest and lowest prevalence in the oropharyngeal tissue was observed in B. Baria and Habiganj, respectively. Two isolates, BAU/KA/DPV(B1)/2014 from Kishoreganj and BAU/KA/DPV(B4)/2014 from Sunamganj were sequenced, and phylogenetic analysis revealed that these isolates are evolutionarily closely related to Chinese isolates of DPV. Additionally, the isolates of DPV BAU/KA/DPV(B1)/2014 and BAU/KA/DPV(B4)/2014 showed the highest (98%) similarity to each other. The nucleotide sequence of the isolate BAU/KA/DPV(B1)/2014 exhibited higher nucleotide variability (246 nucleotides) than that of the vaccine strain (ac

Egg-laying mammals (monotremes) are the only extant mammalian outgroup to therians (marsupial and eutherian animals) and provide key insights into mammalian evolution1,2. Here we generate and analyse reference genomes of the platypus (Ornithorhynchus anatinus) and echidna (Tachyglossus aculeatus), which represent the only two extant monotreme lineages. The nearly complete platypus genome assembly has anchored almost the entire genome onto chromosomes, markedly improving the genome continuity and gene annotation. Together with our echidna sequence, the genomes of the two species allow us to detect the ancestral and lineage-specific genomic changes that shape both monotreme and mammalian evolution. We provide evidence that the monotreme sex chromosome complex originated from an ancestral chromosome ring configuration. The formation of such a unique chromosome complex may have been facilitated by the unusually extensive interactions between the multi-X and multi-Y chromosomes that are shared by the autosomal homologues in humans. Further comparative genomic analyses unravel marked differences between monotremes and therians in haptoglobin genes, lactation genes and chemosensory receptor genes for smell and taste that underlie the ecological adaptation of monotremes.

Small molecules that target the spliceosome SF3B complex are potent inhibitors of cancer cell growth. The compounds affect an early stage of spliceosome assembly when U2 snRNP first engages the branch point sequence of an intron. Employing an inactive herboxidiene analog (iHB) as a competitor, we investigated factors that influence inhibitor interactions with SF3B to interfere with pre-mRNA splicing in vitro. Order-of-addition experiments show that inhibitor interactions are long lasting and affected by both temperature and the presence of ATP. Our data are also consistent with the model that not all SF3B conformations observed in structural studies are conducive to productive inhibitor interactions. Notably, SF3B inhibitors do not impact an ATP-dependent rearrangement in U2 snRNP that exposes the branch binding sequence for base pairing. We also report extended structure-activity relationship analysis of the splicing inhibitor herboxidiene. We identified features of the tetrahydropyran ring that mediate its interactions with SF3B and its ability to interfere with splicing. In the context of recent structures of SF3B bound to inhibitor, our results lead us to extend the model for early spliceosome assembly and inhibitor mechanism. We postulate that interactions between a carboxylic acid substituent of herboxidiene and positively charged SF3B1 side chains in the inhibitor binding channel are needed to maintain inhibitor occupancy while counteracting the SF3B transition to a closed state that is required for stable U2 snRNP interactions with the intron.

Dysregulated mTORC1 signaling alters a wide range of cellular processes, contributing to metabolic disorders and cancer. Defining the molecular details of downstream effectors is thus critical for uncovering selective therapeutic targets. We report that mTORC1 and its downstream kinase S6K enhance eIF4A/4B-mediated translation of Wilms' tumor 1-associated protein (WTAP), an adaptor for the N6-methyladenosine (m6A) RNA methyltransferase complex. This regulation is mediated by 5' UTR of WTAP mRNA that is targeted by eIF4A/4B. Single-nucleotide-resolution m6A mapping revealed that MAX dimerization protein 2 (MXD2) mRNA contains m6A, and increased m6A modification enhances its degradation. WTAP induces cMyc-MAX association by suppressing MXD2 expression, which promotes cMyc transcriptional activity and proliferation of mTORC1-activated cancer cells. These results elucidate a mechanism whereby mTORC1 stimulates oncogenic signaling via m6A RNA modification and illuminates the WTAP-MXD2-cMyc axis as a potential therapeutic target for mTORC1-driven cancers.

BackgroundTranscriptomic data has demonstrated utility to advance the study of physiological diversity and organisms' responses to environmental stressors. However, a lack of genomic resources and challenges associated with collecting high-quality RNA can limit its application for many wild populations. Minimally invasive blood sampling combined with de novo transcriptomic approaches has great potential to alleviate these barriers. Here, we advance these goals for marine turtles by generating high quality de novo blood transcriptome assemblies to characterize functional diversity and compare global transcriptional profiles between tissues, species, and foraging aggregations.ResultsWe generated high quality blood transcriptome assemblies for hawksbill (Eretmochelys imbricata), loggerhead (Caretta caretta), green (Chelonia mydas), and leatherback (Dermochelys coriacea) turtles. The functional diversity in assembled blood transcriptomes was comparable to those from more traditionally sampled tissues. A total of 31.3% of orthogroups identified were present in all four species, representing a core set of conserved genes expressed in blood and shared across marine turtle species. We observed strong species-specific expression of these genes, as well as distinct transcriptomic profiles between green turtle foraging aggregations that inhabit areas of greater or lesser anthropogenic disturbance.ConclusionsObtaining global gene expression data through non-lethal, minimally invasive sampling can greatly expand the applications of RNA-sequencing in protected long-lived species such as marine turtles. The distinct differences in gene expression signatures between species and foraging aggregations provide insight into the functional genomics underlying the diversity in this ancient vertebrate lineage. The transcriptomic resources generated here can be used in further studies examining the evolutionary ecology and anthropogenic impacts on marine turtles.

Spermatogonial stem cells (SSCs) are essential for male fertility. Here, we report that mouse SSC generation is driven by a transcription factor (TF) cascade controlled by the homeobox protein, RHOX10, which acts by driving the differentiation of SSC precursors called pro-spermatogonia (ProSG). We identify genes regulated by RHOX10 in ProSG in vivo and define direct RHOX10-target genes using several approaches, including a rapid temporal induction assay: iSLAMseq. Together, these approaches identify temporal waves of RHOX10 direct targets, as well as RHOX10 secondary-target genes. Many of the RHOX10-regulated genes encode proteins with known roles in SSCs. Using an in vitro ProSG differentiation assay, we find that RHOX10 promotes mouse ProSG differentiation through a conserved transcriptional cascade involving the key germ-cell TFs DMRT1 and ZBTB16. Our study gives important insights into germ cell development and provides a blueprint for how to define TF cascades.

BackgroundExtracellular vesicles (EVs) are regulators of cell-cell interactions and mediators of horizontal transfer of bioactive molecules between cells. EV-mediated cell-cell interactions play roles in physiological and pathophysiological processes, which maybe modulated by exposure to pathogens and cocaine use. However, the effect of pathogens and cocaine use on EV composition and function are not fully understood.ResultsHere, we used systems biology and multi-omics analysis to show that HIV infection (HIV +) and cocaine (COC) use (COC +) promote the release of semen-derived EVs (SEV) with dysregulated extracellular proteome (exProtein), miRNAome (exmiR), and exmiR networks. Integrating SEV proteome and miRNAome revealed a significant decrease in the enrichment of disease-associated, brain-enriched, and HIV-associated miR-128-3p (miR-128) in HIV + COC + SEV with a concomitant increase in miR-128 targets-PEAK1 and RND3/RhoE. Using two-dimensional-substrate single cell haptotaxis, we observed that in the presence of HIV + COC + SEV, contact guidance provided by the extracellular matrix (ECM, collagen type 1) network facilitated far-ranging haptotactic cues that guided monocytes over longer distances. Functionalizing SEV with a miR-128 mimic revealed that the strategic changes in monocyte haptotaxis are in large part the result of SEV-associated miR-128.ConclusionsWe propose that compositionally and functionally distinct HIV + COC + and HIV-COC- SEVs and their exmiR networks may provide cells relevant but divergent haptotactic guidance in the absence of chemotactic cues, under both physiological and pathophysiological conditions.

BackgroundThe all-electronic Single Molecule Break Junction (SMBJ) method is an emerging alternative to traditional polymerase chain reaction (PCR) techniques for genetic sequencing and identification. Existing work indicates that the current spectra recorded from SMBJ experimentations contain unique signatures to identify known sequences from a dataset. However, the spectra are typically extremely noisy due to the stochastic and complex interactions between the substrate, sample, environment, and the measuring system, necessitating hundreds or thousands of experimentations to obtain reliable and accurate results.ResultsThis article presents a DNA sequence identification system based on the current spectra of ten short strand sequences, including a pair that differs by a single mismatch. By employing a gradient boosted tree classifier model trained on conductance histograms, we demonstrate that extremely high accuracy, ranging from approximately 96 % for molecules differing by a single mismatch to 99.5 % otherwise, is possible. Further, such accuracy metrics are achievable in near real-time with just twenty or thirty SMBJ measurements instead of hundreds or thousands. We also demonstrate that a tandem classifier architecture, where the first stage is a multiclass classifier and the second stage is a binary classifier, can be employed to boost the single mismatched pair's identification accuracy to 99.5 %.ConclusionsA monolithic classifier, or more generally, a multistage classifier with model specific parameters that depend on experimental current spectra can be used to successfully identify DNA strands.

Although infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has pleiotropic and systemic effects in some individuals1-3, many others experience milder symptoms. Here, to gain a more comprehensive understanding of the distinction between severe and mild phenotypes in the pathology of coronavirus disease 2019 (COVID-19) and its origins, we performed a whole-blood-preserving single-cell analysis protocol to integrate contributions from all major immune cell types of the blood-including neutrophils, monocytes, platelets, lymphocytes and the contents of the serum. Patients with mild COVID-19 exhibit a coordinated pattern of expression of interferon-stimulated genes (ISGs)3 across every cell population, whereas these ISG-expressing cells are systemically absent in patients with severe disease. Paradoxically, individuals with severe COVID-19 produce very high titres of anti-SARS-CoV-2 antibodies and have a lower viral load compared to individuals with mild disease. Examination of the serum from patients with severe COVID-19 shows that these patients uniquely produce antibodies that functionally block the production of the ISG-expressing cells associated with mild disease, by activating conserved signalling circuits that dampen cellular responses to interferons. Overzealous antibody responses pit the immune system against itself in many patients with COVID-19, and perhaps also in individuals with other viral infections. Our findings reveal potential targets for immunotherapies in patients with severe COVID-19 to re-engage viral defence.

The application of next-generation sequencing extends from microbial identification to epidemiologic insight and antimicrobial resistance prediction. Despite this potential, the roadblock for clinical laboratories lies in implementation and validation of such complex technology and data analysis. Herein, a validation study used whole-genome sequencing (WGS) for pan-bacterial identification (ID) in a clinical laboratory, and assessed its clinical relevance. A diverse set of 125 bacterial isolates, including a subset without genus (25) and/or species (10) ID, were analyzed by de novo assembly and reference genome mapping. The 16S rRNA, rpoB, and groEL genes were used for ID. Using WGS, 100% (89 of 89) and 89% (79 of 89) of isolates were identified at the genus and species-levels, respectively. WGS also provided improved results for 71% (25/35) isolates originally reported with genus-only or descriptive IDs. Chart review identified cases in which improved genus and/or species-level ID by WGS may have had a positive impact on patient care. Reasons included the use of an ineffective antibiotic owing to unclear ID, use of antibiotics when not clinically indicated, and help with an outbreak investigation. The implementation of next-generation sequencing in a clinical microbiology setting is a challenging but necessary task. This study provides a model for the validation and implementation of bacterial ID by WGS in such a setting.

Expression of a few master transcription factors can reprogram the epigenetic landscape and three-dimensional chromatin topology of differentiated cells and achieve pluripotency. During reprogramming, thousands of long-range chromatin contacts are altered, and changes in promoter association with enhancers dramatically influence transcription. Molecular participants at these sites have been identified, but how this re-organization might be orchestrated is not known. Biomolecular condensation is implicated in subcellular organization, including the recruitment of RNA polymerase in transcriptional activation. Here, we show that reprogramming factor KLF4 undergoes biomolecular condensation even in the absence of its intrinsically disordered region. Liquid-liquid condensation of the isolated KLF4 DNA binding domain with a DNA fragment from the NANOG proximal promoter is enhanced by CpG methylation of a KLF4 cognate binding site. We propose KLF4-mediated condensation as one mechanism for selectively organizing and re-organizing the genome based on the local sequence and epigenetic state.

Under the Baltimore nucleic acid-based virus classification scheme, the herpesvirus human cytomegalovirus (HCMV) is a Class I virus, meaning that it contains a double-stranded DNA genome-and no RNA. Here, we report sub-particle cryoEM reconstructions of HCMV virions at 2.9 Å resolution revealing structures resembling non-coding transfer RNAs (tRNAs) associated with the virion's capsid-bound tegument protein, pp150. Through deep sequencing, we show that these RNA sequences match human tRNAs, and we built atomic models using the most abundant tRNA species. Based on our models, tRNA recruitment is mediated by the electrostatic interactions between tRNA phosphate groups and the helix-loop-helix motif of HCMV pp150. The specificity of these interactions may explain the absence of such tRNA densities in murine cytomegalovirus and other human herpesviruses.

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences that provides a single access point to 44 RNA resources and >18 million ncRNA sequences from a wide range of organisms and RNA types. RNAcentral now also includes secondary (2D) structure information for >13 million sequences, making RNAcentral the world's largest RNA 2D structure database. The 2D diagrams are displayed using R2DT, a new 2D structure visualization method that uses consistent, reproducible and recognizable layouts for related RNAs. The sequence similarity search has been updated with a faster interface featuring facets for filtering search results by RNA type, organism, source database or any keyword. This sequence search tool is available as a reusable web component, and has been integrated into several RNAcentral member databases, including Rfam, miRBase and snoDB. To allow for a more fine-grained assignment of RNA types and subtypes, all RNAcentral sequences have been annotated with Sequence Ontology terms. The RNAcentral database continues to grow and provide a central data resource for the RNA community. RNAcentral is freely available at https://rnacentral.org.

Ribosome biogenesis, which takes place mainly in the nucleolus, involves coordinated expression of pre-ribosomal RNAs (pre-rRNAs) and ribosomal proteins, pre-rRNA processing, and subunit assembly with the aid of numerous assembly factors. Our previous study showed that the Arabidopsis thaliana protein arginine methyltransferase AtPRMT3 regulates pre-rRNA processing; however, the underlying molecular mechanism remains unknown. Here, we report that AtPRMT3 interacts with Ribosomal Protein S2 (RPS2), facilitating processing of the 90S/Small Subunit (SSU) processome and repressing nucleolar stress. We isolated an intragenic suppressor of atprmt3-2, which rescues the developmental defects of atprmt3-2 while produces a putative truncated AtPRMT3 protein bearing the entire N-terminus but lacking an intact enzymatic activity domain We further identified RPS2 as an interacting partner of AtPRMT3, and found that loss-of-function rps2a2b mutants were phenotypically reminiscent of atprmt3, showing pleiotropic developmental defects and aberrant pre-rRNA processing. RPS2B binds directly to pre-rRNAs in the nucleus, and such binding is enhanced in atprmt3-2. Consistently, multiple components of the 90S/SSU processome were more enriched by RPS2B in atprmt3-2, which accounts for early pre-rRNA processing defects and results in nucleolar stress. Collectively, our study uncovered a novel mechanism by which AtPRMT3 cooperates with RPS2B to facilitate the dynamic assembly/disassembly of the 90S/SSU processome during ribosome biogenesis and repress nucleolar stress.

Diversity-generating retroelements (DGRs) vary protein sequences to the greatest extent known in the natural world. These elements are encoded by constituents of the human microbiome and the microbial 'dark matter'. Variation occurs through adenine-mutagenesis, in which genetic information in RNA is reverse transcribed faithfully to cDNA for all template bases but adenine. We investigated the determinants of adenine-mutagenesis in the prototypical Bordetella bacteriophage DGR through an in vitro system composed of the reverse transcriptase bRT, Avd protein, and a specific RNA. We found that the catalytic efficiency for correct incorporation during reverse transcription by the bRT-Avd complex was strikingly low for all template bases, with the lowest occurring for adenine. Misincorporation across a template adenine was only somewhat lower in efficiency than correct incorporation. We found that the C6, but not the N1 or C2, purine substituent was a key determinant of adenine-mutagenesis. bRT-Avd was insensitive to the C6 amine of adenine but recognized the C6 carbonyl of guanine. We also identified two bRT amino acids predicted to nonspecifically contact incoming dNTPs, R74 and I181, as promoters of adenine-mutagenesis. Our results suggest that the overall low catalytic efficiency of bRT-Avd is intimately tied to its ability to carry out adenine-mutagenesis.