Search

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.

Cases of Leber congenital amaurosis caused by mutations in CRX (LCA7) exhibit an early form of the disease and show signs of significant photoreceptor dysfunction and eventual loss. To establish a translational in vitro model system to study gene-editing-based therapies, we generated LCA7 retinal organoids harboring a dominant disease-causing mutation in CRX. Our LCA7 retinal organoids develop signs of immature and dysfunctional photoreceptor cells, providing us with a reliable in vitro model to recapitulate LCA7. Furthermore, we performed a proof-of-concept study in which we utilize allele-specific CRISPR/Cas9-based gene editing to knock out mutant CRX and saw moderate rescue of photoreceptor phenotypes in our organoids. This work provides early evidence for an effective approach to treat LCA7, which can be applied more broadly to other dominant genetic diseases.

The human brain is subdivided into distinct anatomical structures, including the neocortex, which in turn encompasses dozens of distinct specialized cortical areas. Early morphogenetic gradients are known to establish early brain regions and cortical areas, but how early patterns result in finer and more discrete spatial differences remains poorly understood1. Here we use single-cell RNA sequencing to profile ten major brain structures and six neocortical areas during peak neurogenesis and early gliogenesis. Within the neocortex, we find that early in the second trimester, a large number of genes are differentially expressed across distinct cortical areas in all cell types, including radial glia, the neural progenitors of the cortex. However, the abundance of areal transcriptomic signatures increases as radial glia differentiate into intermediate progenitor cells and ultimately give rise to excitatory neurons. Using an automated, multiplexed single-molecule fluorescent in situ hybridization approach, we find that laminar gene-expression patterns are highly dynamic across cortical regions. Together, our data suggest that early cortical areal patterning is defined by strong, mutually exclusive frontal and occipital gene-expression signatures, with resulting gradients giving rise to the specification of areas between these two poles throughout successive developmental timepoints.

Retrotransposons mediate gene regulation in important developmental and pathological processes. Here, we characterized the transient retrotransposon induction during preimplantation development of eight mammals. Induced retrotransposons exhibit similar preimplantation profiles across species, conferring gene regulatory activities, particularly through long terminal repeat (LTR) retrotransposon promoters. A mouse-specific MT2B2 retrotransposon promoter generates an N-terminally truncated Cdk2ap1ΔN that peaks in preimplantation embryos and promotes proliferation. In contrast, the canonical Cdk2ap1 peaks in mid-gestation and represses cell proliferation. This MT2B2 promoter, whose deletion abolishes Cdk2ap1ΔN production, reduces cell proliferation and impairs embryo implantation, is developmentally essential. Intriguingly, Cdk2ap1ΔN is evolutionarily conserved in sequence and function yet is driven by different promoters across mammals. The distinct preimplantation Cdk2ap1ΔN expression in each mammalian species correlates with the duration of its preimplantation development. Hence, species-specific transposon promoters can yield evolutionarily conserved, alternative protein isoforms, bestowing them with new functions and species-specific expression to govern essential biological divergence.

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.

Estrogen receptor α (ERα) is a hormone receptor and key driver for over 70% of breast cancers that has been studied for decades as a transcription factor. Unexpectedly, we discover that ERα is a potent non-canonical RNA-binding protein. We show that ERα RNA binding function is uncoupled from its activity to bind DNA and critical for breast cancer progression. Employing genome-wide cross-linking immunoprecipitation (CLIP) sequencing and a functional CRISPRi screen, we find that ERα-associated mRNAs sustain cancer cell fitness and elicit cellular responses to stress. Mechanistically, ERα controls different steps of RNA metabolism. In particular, we demonstrate that ERα RNA binding mediates alternative splicing of XBP1 and translation of the eIF4G2 and MCL1 mRNAs, which facilitates survival upon stress conditions and sustains tamoxifen resistance of cancer cells. ERα is therefore a multifaceted RNA-binding protein, and this activity transforms our knowledge of post-transcriptional regulation underlying cancer development and drug response.

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.

The PD-L1 overexpression is an important event of immune escape and metastasis in triple-negative breast cancer (TNBC), but the molecular mechanism remains to be determined. Interferon gamma (IFNγ) represents a major driving force behind PD-L1 expression in tumor microenvironment, and histone deacetylase 2 (HDAC2) is required for IFN signaling. Here, we investigated the regulation of HDAC2 on the IFNγ-induced PD-L1 expression in TNBC cells. We found the HDAC2 and PD-L1 expression in TNBC was significantly higher than that in non-TNBC, and HDAC2 was positively correlated with PD-L1 expression. HDAC2 promoted PD-L1 induction by upregulating the phosphorylation of JAK1, JAK2, and STAT1, as well as the translocation of STAT1 to the nucleus and the recruitment of STAT1 to the PD-L1 promoter. Meanwhile, HDAC2 was recruited to the PD-L1 promoter by STAT1, and HDAC2 knockout compromised IFNγ-induced upregulation of H3K27, H3K9 acetylation, and the BRD4 recruitment in PD-L1 promoter. In addition, significant inhibition of proliferation, colony formation, migration, and cell cycle of TNBC cells were observed following knockout of HDAC2 in vitro. Furthermore, HDAC2 knockout reduced IFNγ-induced PD-L1 expression, lymphocyte infiltration, and retarded tumor growth and metastasis in the breast cancer mouse models. This study may provide evidence that HDAC2 promotes IFNγ-induced PD-L1 expression, suggesting a way for enhanced antitumor immunity when targeting the HDAC2 in TNBC.

Bacterial communities play critical roles across all of Earth's biomes, affecting human health and global ecosystem functioning. They do so under strong constraints exerted by viruses, that is, bacteriophages or 'phages'. Phages can reshape bacterial communities' structure, influence long-term evolution of bacterial populations, and alter host cell metabolism during infection. Metagenomics approaches, that is, shotgun sequencing of environmental DNA or RNA, recently enabled large-scale exploration of phage genomic diversity, yielding several millions of phage genomes now to be further analyzed and characterized. One major challenge however is the lack of direct host information for these phages. Several methods and tools have been proposed to bioinformatically predict the potential host(s) of uncultivated phages based only on genome sequence information. Here we review these different approaches and highlight their distinct strengths and limitations. We also outline complementary experimental assays which are being proposed to validate and refine these bioinformatic predictions.

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.

Platelet aggregation at the site of atherosclerotic vascular injury is the underlying pathophysiology of myocardial infarction and stroke. To build upon prior GWAS, here we report on 16 loci identified through a whole genome sequencing (WGS) approach in 3,855 NHLBI Trans-Omics for Precision Medicine (TOPMed) participants deeply phenotyped for platelet aggregation. We identify the RGS18 locus, which encodes a myeloerythroid lineage-specific regulator of G-protein signaling that co-localizes with expression quantitative trait loci (eQTL) signatures for RGS18 expression in platelets. Gene-based approaches implicate the SVEP1 gene, a known contributor of coronary artery disease risk. Sentinel variants at RGS18 and PEAR1 are associated with thrombosis risk and increased gastrointestinal bleeding risk, respectively. Our WGS findings add to previously identified GWAS loci, provide insights regarding the mechanism(s) by which genetics may influence cardiovascular disease risk, and underscore the importance of rare variant and regulatory approaches to identifying loci contributing to complex phenotypes.

Given the global impact and severity of COVID-19, there is a pressing need for a better understanding of the SARS-CoV-2 genome and mutations. Multi-strain sequence alignments of coronaviruses (CoV) provide important information for interpreting the genome and its variation. We apply a comparative genomics method, ConsHMM, to the multi-strain alignments of CoV to annotate every base of the SARS-CoV-2 genome with conservation states based on sequence alignment patterns among CoV. The learned conservation states show distinct enrichment patterns for genes, protein domains, and other regions of interest. Certain states are strongly enriched or depleted of SARS-CoV-2 mutations, which can be used to predict potentially consequential mutations. We expect the conservation states to be a resource for interpreting the SARS-CoV-2 genome and mutations.

Using the Cap Analysis of Gene Expression (CAGE) technology, the FANTOM5 consortium provided one of the most comprehensive maps of transcription start sites (TSSs) in several species. Strikingly, ~72% of them could not be assigned to a specific gene and initiate at unconventional regions, outside promoters or enhancers. Here, we probe these unassigned TSSs and show that, in all species studied, a significant fraction of CAGE peaks initiate at microsatellites, also called short tandem repeats (STRs). To confirm this transcription, we develop Cap Trap RNA-seq, a technology which combines cap trapping and long read MinION sequencing. We train sequence-based deep learning models able to predict CAGE signal at STRs with high accuracy. These models unveil the importance of STR surrounding sequences not only to distinguish STR classes, but also to predict the level of transcription initiation. Importantly, genetic variants linked to human diseases are preferentially found at STRs with high transcription initiation level, supporting the biological and clinical relevance of transcription initiation at STRs. Together, our results extend the repertoire of non-coding transcription associated with DNA tandem repeats and complexify STR polymorphism.

Functional biomolecules, such as RNA, encapsulated inside a protocellular membrane are believed to have comprised a very early, critical stage in the evolution of life, since membrane vesicles allow selective permeability and create a unit of selection enabling cooperative phenotypes. The biophysical environment inside a protocell would differ fundamentally from bulk solution due to the microscopic confinement. However, the effect of the encapsulated environment on ribozyme evolution has not been previously studied experimentally. Here, we examine the effect of encapsulation inside model protocells on the self-aminoacylation activity of tens of thousands of RNA sequences using a high-throughput sequencing assay. We find that encapsulation of these ribozymes generally increases their activity, giving encapsulated sequences an advantage over nonencapsulated sequences in an amphiphile-rich environment. In addition, highly active ribozymes benefit disproportionately more from encapsulation. The asymmetry in fitness gain broadens the distribution of fitness in the system. Consistent with Fisher's fundamental theorem of natural selection, encapsulation therefore leads to faster adaptation when the RNAs are encapsulated inside a protocell during in vitro selection. Thus, protocells would not only provide a compartmentalization function but also promote activity and evolutionary adaptation during the origin of life.

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.

Most eukaryotes harbor two distinct pre-mRNA splicing machineries: the major spliceosome, which removes >99% of introns, and the minor spliceosome, which removes rare, evolutionarily conserved introns. Although hypothesized to serve important regulatory functions, physiologic roles of the minor spliceosome are not well understood. For example, the minor spliceosome component ZRSR2 is subject to recurrent, leukemia-associated mutations, yet functional connections among minor introns, hematopoiesis and cancers are unclear. Here, we identify that impaired minor intron excision via ZRSR2 loss enhances hematopoietic stem cell self-renewal. CRISPR screens mimicking nonsense-mediated decay of minor intron-containing mRNA species converged on LZTR1, a regulator of RAS-related GTPases. LZTR1 minor intron retention was also discovered in the RASopathy Noonan syndrome, due to intronic mutations disrupting splicing and diverse solid tumors. These data uncover minor intron recognition as a regulator of hematopoiesis, noncoding mutations within minor introns as potential cancer drivers and links among ZRSR2 mutations, LZTR1 regulation and leukemias.

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.

Ultraconserved enhancer sequences show perfect conservation between human and rodent genomes, suggesting that their functions are highly sensitive to mutation. However, current models of enhancer function do not sufficiently explain this extreme evolutionary constraint. We subjected 23 ultraconserved enhancers to different levels of mutagenesis, collectively introducing 1,547 mutations, and examined their activities in transgenic mouse reporter assays. Overall, we find that the regulatory properties of ultraconserved enhancers are robust to mutation. Upon mutagenesis, nearly all (19/23, 83%) still functioned as enhancers at one developmental stage, as did most of those tested again later in development (5/9, 56%). Replacement of endogenous enhancers with mutated alleles in mice corroborated results of transgenic assays, including the functional resilience of ultraconserved enhancers to mutation. Our findings show that the currently known activities of ultraconserved enhancers do not necessarily require the perfect conservation observed in evolution and suggest that additional regulatory or other functions contribute to their sequence constraint.

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.

APOBEC mutagenesis, a major driver of cancer evolution, is known for targeting TpC sites in DNA. Recently, we showed that APOBEC3A (A3A) targets DNA hairpin loops. Here, we show that DNA secondary structure is in fact an orthogonal influence on A3A substrate optimality and, surprisingly, can override the TpC sequence preference. VpC (non-TpC) sites in optimal hairpins can outperform TpC sites as mutational hotspots. This expanded understanding of APOBEC mutagenesis illuminates the genomic Twin Paradox, a puzzling pattern of closely spaced mutation hotspots in cancer genomes, in which one is a canonical TpC site but the other is a VpC site, and double mutants are seen only in trans, suggesting a two-hit driver event. Our results clarify this paradox, revealing that both hotspots in these twins are optimal A3A substrates. Our findings reshape the notion of a mutation signature, highlighting the additive roles played by DNA sequence and DNA structure.

Studies of acute myeloid leukemia rely on DNA sequencing and immunophenotyping by flow cytometry as primary tools for disease characterization. However, leukemia tumor heterogeneity complicates integration of DNA variants and immunophenotypes from separate measurements. Here we introduce DAb-seq, a technology for simultaneous capture of DNA genotype and cell surface phenotype from single cells at high throughput, enabling direct profiling of proteogenomic states in tens of thousands of cells. To demonstrate the approach, we analyze the disease of three patients with leukemia over multiple treatment timepoints and disease recurrences. We observe complex genotype-phenotype dynamics that illustrate the subtlety of the disease process and the degree of incongruity between blast cell genotype and phenotype in different clinical scenarios. Our results highlight the importance of combined single-cell DNA and protein measurements to fully characterize the heterogeneity of leukemia.

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.

It has long been known that noncoding genomic regions can be obligate cis elements acted upon in trans by gene products. In viruses, cis elements regulate gene expression, encapsidation, and other maturation processes, but mapping these elements relies on targeted iterative deletion or laborious prospecting for rare spontaneously occurring mutants. Here, we introduce a method to comprehensively map viral cis and trans elements at single-nucleotide resolution by high-throughput random deletion. Variable-size deletions are randomly generated by transposon integration, excision, and exonuclease chewback and then barcoded for tracking via sequencing (i.e., random deletion library sequencing [RanDeL-seq]). Using RanDeL-seq, we generated and screened >23,000 HIV-1 variants to generate a single-base resolution map of HIV-1's cis and trans elements. The resulting landscape recapitulated HIV-1's known cis-acting elements (i.e., long terminal repeat [LTR], Ψ, and Rev response element [RRE]) and, surprisingly, indicated that HIV-1's central DNA flap (i.e., central polypurine tract [cPPT] to central termination sequence [CTS]) is as critical as the LTR, Ψ, and RRE for long-term passage. Strikingly, RanDeL-seq identified a previously unreported ∼300-bp region downstream of RRE extending to splice acceptor 7 that is equally critical for sustained viral passage. RanDeL-seq was also used to construct and screen a library of >90,000 variants of Zika virus (ZIKV). Unexpectedly, RanDeL-seq indicated that ZIKV's cis-acting regions are larger than the untranscribed (UTR) termini, encompassing a large fraction of the nonstructural genes. Collectively, RanDeL-seq provides a versatile framework for generating viral deletion mutants, enabling discovery of replication mechanisms and development of novel antiviral therapeutics, particularly for emerging viral infections.IMPORTANCE Recent studies have renewed interest in developing novel antiviral therapeutics and vaccines based on defective interfering particles (DIPs)-a subset of viral deletion mutants that conditionally replicate. Identifying and engineering DIPs require that viral cis- and trans-acting elements be accurately mapped. Here, we introduce a high-throughput method (random deletion library sequencing [RanDeL-seq]) to comprehensively map cis- and trans-acting elements within a viral genome. RanDeL-seq identified essential cis elements in HIV, including the obligate nature of the once-controversial viral central polypurine tract (cPPT), and identified a new cis region proximal to the Rev responsive element (RRE). RanDeL-seq also identified regions of Zika virus required for replication and packaging. RanDeL-seq is a versatile and comprehensive technique to rapidly map cis and trans regions of a genome.

Viral genome sequencing has guided our understanding of the spread and extent of genetic diversity of SARS-CoV-2 during the COVID-19 pandemic. SARS-CoV-2 viral genomes are usually sequenced from nasopharyngeal swabs of individual patients to track viral spread. Recently, RT-qPCR of municipal wastewater has been used to quantify the abundance of SARS-CoV-2 in several regions globally. However, metatranscriptomic sequencing of wastewater can be used to profile the viral genetic diversity across infected communities. Here, we sequenced RNA directly from sewage collected by municipal utility districts in the San Francisco Bay Area to generate complete and nearly complete SARS-CoV-2 genomes. The major consensus SARS-CoV-2 genotypes detected in the sewage were identical to clinical genomes from the region. Using a pipeline for single nucleotide variant calling in a metagenomic context, we characterized minor SARS-CoV-2 alleles in the wastewater and detected viral genotypes which were also found within clinical genomes throughout California. Observed wastewater variants were more similar to local California patient-derived genotypes than they were to those from other regions within the United States or globally. Additional variants detected in wastewater have only been identified in genomes from patients sampled outside California, indicating that wastewater sequencing can provide evidence for recent introductions of viral lineages before they are detected by local clinical sequencing. These results demonstrate that epidemiological surveillance through wastewater sequencing can aid in tracking exact viral strains in an epidemic context.

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.

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.

Viruses are integral components of all ecosystems and microbiomes on Earth. Through pervasive infections of their cellular hosts, viruses can reshape microbial community structure and drive global nutrient cycling. Over the past decade, viral sequences identified from genomes and metagenomes have provided an unprecedented view of viral genome diversity in nature. Since 2016, the IMG/VR database has provided access to the largest collection of viral sequences obtained from (meta)genomes. Here, we present the third version of IMG/VR, composed of 18 373 cultivated and 2 314 329 uncultivated viral genomes (UViGs), nearly tripling the total number of sequences compared to the previous version. These clustered into 935 362 viral Operational Taxonomic Units (vOTUs), including 188 930 with two or more members. UViGs in IMG/VR are now reported as single viral contigs, integrated proviruses or genome bins, and are annotated with a new standardized pipeline including genome quality estimation using CheckV, taxonomic classification reflecting the latest ICTV update, and expanded host taxonomy prediction. The new IMG/VR interface enables users to efficiently browse, search, and select UViGs based on genome features and/or sequence similarity. IMG/VR v3 is available at https://img.jgi.doe.gov/vr, and the underlying data are available to download at https://genome.jgi.doe.gov/portal/IMG_VR.

Advances in sequencing technology have allowed researchers to sequence DNA with greater ease and at decreasing costs. Main developments have focused on either sequencing many short sequences or fewer large sequences. Methods for sequencing mid-sized sequences of 600-5,000 bp are currently less efficient. For example, the PacBio Sequel I system yields ~ 100,000-300,000 reads with an accuracy per base pair of 90-99%. We sought to sequence several DNA populations of ~ 870 bp in length with a sequencing accuracy of 99% and to the greatest depth possible. We optimised a simple, robust method to concatenate genes of ~ 870 bp five times and then sequenced the resulting DNA of ~ 5,000 bp by PacBioSMRT long-read sequencing. Our method improved upon previously published concatenation attempts, leading to a greater sequencing depth, high-quality reads and limited sample preparation at little expense. We applied this efficient concatenation protocol to sequence nine DNA populations from a protein engineering study. The improved method is accompanied by a simple and user-friendly analysis pipeline, DeCatCounter, to sequence medium-length sequences efficiently at one-fifth of the cost.

Of different types of noncoding RNAs, microRNAs (miRNAs) have arguably been in the spotlight over the last decade. As post-transcriptional regulators of gene expression, miRNAs play key roles in various cellular pathways, including both development and response to a/biotic stress, such as drought and diseases. Having high-quality reference genome sequences enabled identification and annotation of miRNAs in several plant species, where miRNA sequences are highly conserved. As computational miRNA identification and annotation processes are mostly error-prone processes, homology-based predictions increase prediction accuracy. We developed and have improved the miRNA annotation pipeline, SUmir, in the last decade, which has been used for several plant genomes since then. This study presents a fully automated, new miRNA pipeline, mirMachine (miRNA Machine), by (i) adding an additional filtering step on the secondary structure predictions, (ii) making it fully automated, and (iii) introducing new options to predict either known miRNA based on homology or novel miRNAs based on small RNA sequencing reads using the previous pipeline. The new miRNA pipeline, mirMachine, was tested using The Arabidopsis Information Resource, TAIR10, release of the Arabidopsis genome and the International Wheat Genome Sequencing Consortium (IWGSC) wheat reference genome v2.

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 (accession no. EU082088), which may affect protein function and additional drug sensitivity. CONCLUSION: Based on the findings of the molecular study, it can be assumed that the Bangladeshi isolates and all Chinese isolates of DPV may have a common ancestry.

Research on RNA function and therapeutic potential is dominated by the use of chemoengineered RNA mimics. Recent efforts have led to the establishment of novel technologies for the production of recombinant or bioengineered RNA molecules, which should better recapitulate the structures, functions and safety profiles of natural RNAs because both are produced and folded in living cells. Herein, we describe a robust approach for reproducible fermentation production of bioengineered RNA agents (BERAs) carrying warhead miRNAs, siRNAs, aptamers, or other forms of small RNAs, based upon an optimal hybrid tRNA/pre-miRNA carrier. Target BERA/sRNAs are readily purified by fast protein liquid chromatography (FPLC) to a high degree of homogeneity (>97%). This approach offers a consistent high-level expression (>30% of total bacterial RNAs) and large-scale production of ready-to-use BERAs (multiple to tens milligrams from 1 L bacterial culture).

The functional annotation of livestock genomes is crucial for understanding the molecular mechanisms that underpin complex traits of economic importance, adaptive evolution and comparative genomics. Here, we provide the most comprehensive catalogue to date of regulatory elements in the pig (Sus scrofa) by integrating 223 epigenomic and transcriptomic data sets, representing 14 biologically important tissues. We systematically describe the dynamic epigenetic landscape across tissues by functionally annotating 15 different chromatin states and defining their tissue-specific regulatory activities. We demonstrate that genomic variants associated with complex traits and adaptive evolution in pig are significantly enriched in active promoters and enhancers. Furthermore, we reveal distinct tissue-specific regulatory selection between Asian and European pig domestication processes. Compared with human and mouse epigenomes, we show that porcine regulatory elements are more conserved in DNA sequence, under both rapid and slow evolution, than those under neutral evolution across pig, mouse, and human. Finally, we provide biological insights on tissue-specific regulatory conservation, and by integrating 47 human genome-wide association studies, we demonstrate that, depending on the traits, mouse or pig might be more appropriate biomedical models for different complex traits and diseases.

Aberrant alternative splicing (AS) is a hallmark of cancer and a potential target for novel anti-cancer therapeutics. Breast cancer-associated AS events are known to be linked to disease progression, metastasis, and survival of breast cancer patients. To identify altered AS programs occurring in metastatic breast cancer, we perform a global analysis of AS events by using RNA-mediated oligonucleotide annealing, selection, and ligation coupled with next-generation sequencing (RASL-seq). We demonstrate that, relative to low-metastatic, high-metastatic breast cancer cells show different AS choices in genes related to cancer progression. Supporting a global reshape of cancer-related splicing profiles in metastatic breast cancer we found an enrichment of RNA-binding motifs recognized by several splicing regulators, which have aberrant expression levels or activity during breast cancer progression, including SRSF1. Among SRSF1-regulated targets we found DCUN1D5, a gene for which skipping of exon 4 in its pre-mRNA introduces a premature termination codon (PTC), thus generating an unstable transcript degraded by nonsense-mediated mRNA decay (NMD). Significantly, distinct breast cancer subtypes show different DCUN1D5 isoform ratios with metastatic breast cancer expressing the highest level of the NMD-insensitive DCUN1D5 mRNA, thus showing high DCUN1D5 expression levels, which are ultimately associated with poor overall and relapse-free survival in breast cancer patients. Collectively, our results reveal global AS features of metastatic breast tumors, which open new possibilities for the treatment of these aggressive tumor types.

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.

Leafy head is a unique type of plant architecture found in some vegetable crops, with leaves bending inward to form a compact head. The genetic and molecular mechanisms underlying leafy head in vegetables remain poorly understood. We genetically fine-mapped and cloned a major quantitative trait locus controlling heading in lettuce. The candidate gene (LsKN1) is a homolog of knotted 1 (KN1) from Zea mays Complementation and CRISPR/Cas9 knockout experiments confirmed the role of LsKN1 in heading. In heading lettuce, there is a CACTA-like transposon inserted into the first exon of LsKN1 (LsKN1▽). The transposon sequences act as a promoter rather than an enhancer and drive high expression of LsKN1▽. The enhanced expression of LsKN1▽ is necessary but not sufficient for heading in lettuce. Data from ChIP-sequencing, electrophoretic mobility shift assays, and dual luciferase assays indicate that the LsKN1▽ protein binds the promoter of LsAS1 and down-regulates its expression to alter leaf dorsoventrality. This study provides insight into plant leaf development and will be useful for studies on heading in other vegetable crops.

Axon injury is a hallmark of many neurodegenerative diseases, often resulting in neuronal cell death and functional impairment. Dual leucine zipper kinase (DLK) has emerged as a key mediator of this process. However, while DLK inhibition is robustly protective in a wide range of neurodegenerative disease models, it also inhibits axonal regeneration. Indeed, there are no genetic perturbations that are known to both improve long-term survival and promote regeneration. To identify such a neuroprotective target, we conducted a set of complementary high-throughput screens using a protein kinase inhibitor library in human stem cell-derived retinal ganglion cells (hRGCs). Overlapping compounds that promoted both neuroprotection and neurite outgrowth were bioinformatically deconvoluted to identify specific kinases that regulated neuronal death and axon regeneration. This work identified the role of germinal cell kinase four (GCK-IV) kinases in cell death and additionally revealed their unexpected activity in suppressing axon regeneration. Using an adeno-associated virus (AAV) approach, coupled with genome editing, we validated that GCK-IV kinase knockout improves neuronal survival, comparable to that of DLK knockout, while simultaneously promoting axon regeneration. Finally, we also found that GCK-IV kinase inhibition also prevented the attrition of RGCs in developing retinal organoid cultures without compromising axon outgrowth, addressing a major issue in the field of stem cell-derived retinas. Together, these results demonstrate a role for the GCK-IV kinases in dissociating the cell death and axonal outgrowth in neurons and their druggability provides for therapeutic options for neurodegenerative diseases.

The abnormal PI3K/AKT/mTOR pathway is one of the most common genomic abnormalities in breast cancers including triple-negative breast cancer (TNBC), and pharmacologic inhibition of these aberrations has shown activity in TNBC patients. Here, we designed and identified a small-molecule Comp34 that suppresses both AKT and mTOR protein expression and exhibits robust cytotoxicity towards TNBC cells but not nontumorigenic normal breast epithelial cells. Mechanically, long noncoding RNA (lncRNA) AL354740.1-204 (also named as NUDT3-AS4) acts as a microRNA sponge to compete with AKT1/mTOR mRNAs for binding to miR-99s, leading to decrease in degradation of AKT1/mTOR mRNAs and subsequent increase in AKT1/mTOR protein expression. Inhibition of lncRNA-NUDT3-AS4 and suppression of the NUDT3-AS4/miR-99s association contribute to Comp34-affected biologic pathways. In addition, Comp34 alone is effective in cells with secondary resistance to rapamycin, the best-known inhibitor of mTOR, and displays a greater in vivo antitumor efficacy and lower toxicity than rapamycin in TNBC xenografted models. In conclusion, NUDT3-AS4 may play a proproliferative role in TNBC and be considered a relevant therapeutic target, and Comp34 presents promising activity as a single agent to inhibit TNBC through regulation of NUDT3-AS4 and miR-99s.

Fluid intake is an essential innate behaviour that is mainly caused by two distinct types of thirst1-3. Increased blood osmolality induces osmotic thirst that drives animals to consume pure water. Conversely, the loss of body fluid induces hypovolaemic thirst, in which animals seek both water and minerals (salts) to recover blood volume. Circumventricular organs in the lamina terminalis are critical sites for sensing both types of thirst-inducing stimulus4-6. However, how different thirst modalities are encoded in the brain remains unknown. Here we employed stimulus-to-cell-type mapping using single-cell RNA sequencing to identify the cellular substrates that underlie distinct types of thirst. These studies revealed diverse types of excitatory and inhibitory neuron in each circumventricular organ structure. We show that unique combinations of these neuron types are activated under osmotic and hypovolaemic stresses. These results elucidate the cellular logic that underlies distinct thirst modalities. Furthermore, optogenetic gain of function in thirst-modality-specific cell types recapitulated water-specific and non-specific fluid appetite caused by the two distinct dipsogenic stimuli. Together, these results show that thirst is a multimodal physiological state, and that different thirst states are mediated by specific neuron types in the mammalian brain.

Annotations of evolutionary sequence constraint based on multi-species genome alignments and genome-wide maps of epigenomic marks and transcription factor binding provide important complementary information for understanding the human genome and genetic variation. Here we developed the Constrained Non-Exonic Predictor (CNEP) to quantify the evidence of each base in the genome being in an evolutionarily constrained non-exonic element from an input of over 60,000 epigenomic and transcription factor binding features. We find that the CNEP score outperforms baseline and related existing scores at predicting evolutionarily constrained non-exonic bases from such data. However, a subset of them are still not well predicted by CNEP. We developed a complementary Conservation Signature Score by CNEP (CSS-CNEP) that is predictive of those bases. We further characterize the nature of constrained non-exonic bases with low CNEP scores using additional types of information. CNEP and CSS-CNEP are resources for analyzing constrained non-exonic bases in the genome.

We describe an RT-PCR protocol that allows high-resolution mapping of splicing products and isoforms using fluorescently labeled primers. Each species contains one fluorescent group allowing a direct comparison of the different isoforms despite size differences. A custom-size ladder enables the precise determination of cDNA lengths and discrimination of isoforms differing by less than five nucleotides on polyacrylamide gels. This protocol also allows the detection of products from in vitro splicing reactions, circumventing the need to use radiolabeled transcripts. For complete details on the use and execution of this protocol, please refer to Gabunilas and Chanfreau (2016).

Exoribonuclease-resistant RNAs (xrRNAs) are discrete elements that block the progression of 5' to 3' exoribonucleases using specifically folded RNA structures. A recently discovered class of xrRNA is widespread in several genera of plant-infecting viruses, within both noncoding and protein-coding subgenomic RNAs. The structure of one such xrRNA from a dianthovirus revealed three-dimensional details of the resistant fold but did not answer all questions regarding the conservation and diversity of this xrRNA class. Here, we present the crystal structure of a representative polerovirus xrRNA that contains sequence elements that diverge from the previously solved structure. This new structure rationalizes previously unexplained sequence conservation patterns and shows interactions not present in the first structure. Together, the structures of these xrRNAs from dianthovirus and polerovirus genera support the idea that these plant virus xrRNAs fold through a defined pathway that includes a programmed intermediate conformation. This work deepens our knowledge of the structure-function relationship of xrRNAs and shows how evolution can craft similar RNA folds from divergent sequences.

BackgroundWe have previously reported 18S and 25S ribosomal RNA molecules in Candida albicans resistant to processive 5' → 3' exonuclease, appearing as cells approached stationary growth phase. Initial analysis pointed to extra phosphate(s) at their 5'- end raising the possibility that they were newly transcribed. Here we report on additional experiments exploring this possibility and try to establish which of the RNA polymerases may be transcribing them.ResultsOligo-ligation and primer extension again showed the presence of extra phosphate at the 5'-end of the reported processing sites for both 18S and 25S ribosomal RNA components. Inhibition of Pol I with BMH-21 increased the presence of the molecules. Quantitation with an Agilent Bioanalyzer showed that resistant 18S and 25S molecules are primarily produced in the nucleus. Utilizing an RNA cap specific antibody, a signal could be detected on these molecules via immunoblotting; such signal could be eliminated by decapping reaction. Both the cap specific antibody and eIF4E cap-binding protein, increased fold enrichment upon quantitative amplification. Antibodies specific for the RNA Polymerase II c-terminal domain and TFIIB initiator factor showed the presence of Pol II on DNA sequences for both 18S and 25S molecules in chromatin precipitation and qPCR assays. Rapamycin inhibition of TOR complex also resulted in an increase of resistant 18S and 25S molecules.ConclusionsThese data raise the possibility of a role for RNA Polymerase II in the production of 18S and 25S molecules and indicate that efforts for more direct proof may be worthwhile. If definitively proven it will establish an additional role for RNA Polymerase II in ribosomal production.

CRISPR-Cas technologies have enabled programmable gene editing in eukaryotes and prokaryotes. However, the leading Cas9 and Cas12a enzymes are limited in their ability to make large deletions. Here, we used the processive nuclease Cas3, together with a minimal Type I-C Cascade-based system for targeted genome engineering in bacteria. DNA cleavage guided by a single CRISPR RNA generated large deletions (7-424 kilobases) in Pseudomonas aeruginosa with near-100% efficiency, while Cas9 yielded small deletions and point mutations. Cas3 generated bidirectional deletions originating from the programmed site, which was exploited to reduce the P. aeruginosa genome by 837 kb (13.5%). Large deletion boundaries were efficiently specified by a homology-directed repair template during editing with Cascade-Cas3, but not Cas9. A transferable 'all-in-one' vector was functional in Escherichia coli, Pseudomonas syringae and Klebsiella pneumoniae, and endogenous CRISPR-Cas use was enhanced with an 'anti-anti-CRISPR' strategy. P. aeruginosa Type I-C Cascade-Cas3 (PaeCas3c) facilitates rapid strain manipulation with applications in synthetic biology, genome minimization and the removal of large genomic regions.

DNA methyltransferase DNMT3A is essential for establishment of mammalian DNA methylation during development. The R882H DNMT3A is a hotspot mutation in acute myeloid leukemia (AML) causing aberrant DNA methylation. However, how this mutation affects the structure and function of DNMT3A remains unclear. Here we report structural characterization of wild-type and R882H-mutated DNMT3A in complex with DNA substrates with different sequence contexts. A loop from the target recognition domain (TRD loop) recognizes the CpG dinucleotides in a +1 flanking site-dependent manner. The R882H mutation reduces the DNA binding at the homodimeric interface, as well as the molecular link between the homodimeric interface and TRD loop, leading to enhanced dynamics of TRD loop. Consistently, in vitro methylation analyses indicate that the R882H mutation compromises the enzymatic activity, CpG specificity and flanking sequence preference of DNMT3A. Together, this study uncovers multiple defects of DNMT3A caused by the R882H mutation in AML.

BackgroundThe syphilis epidemic continues to cause substantial morbidity and mortality worldwide, particularly in low- and middle-income countries, despite several recent disease control initiatives. Though our understanding of the pathogenesis of this disease and the biology of the syphilis agent, Treponema pallidum subsp. pallidum has improved over the last two decades, further research is necessary to improve clinical diagnosis and disease management protocols. Additionally, such research efforts could contribute to the identification of possible targets for the development of an effective vaccine to stem syphilis spread.MethodsThis study will recruit two cohorts of participants with active syphilis infection, one with de novo infection, one with repeat infection. Whole blood specimens will be collected from each study participant at baseline, 4, 12, 24, 36, and 48 weeks, to track specific markers of their immunological response, as well as to compare humoral reactivity to Treponema pallidum antigens between the two groups. Additionally, we will use serum specimens to look for unique cytokine patterns in participants with early syphilis. Oral and blood samples, as well as samples from any syphilitic lesions present, will also be collected to sequence any Treponema pallidum DNA found.DiscussionBy furthering our understanding of syphilis pathogenesis and human host immune response to Treponema pallidum, we will provide important data that will help in development of new point-of-care tests that could better identify active infection, leading to improved syphilis diagnosis and management. Findings could also contribute to vaccine development efforts.

BackgroundIntermediate filament proteins that construct the nuclear lamina of a cell include the Lamin A/C proteins encoded by the LMNA gene, and are implicated in fundamental processes such as nuclear structure, gene expression, and signal transduction. LMNA mutations predominantly affect mesoderm-derived cell lineages in diseases collectively termed as laminopathies that include dilated cardiomyopathy with conduction defects, different forms of muscular dystrophies, and premature aging syndromes as Hutchinson-Gilford Progeria Syndrome. At present, our understanding of the molecular mechanisms regulating tissue-specific manifestations of laminopathies are still limited.MethodsTo gain deeper insight into the molecular mechanism of a novel LMNA splice-site mutation (c.357-2A > G) in an affected family with cardiac disease, we conducted deep RNA sequencing and pathway analysis for nine fibroblast samples obtained from three patients with cardiomyopathy, three unaffected family members, and three unrelated, unaffected individuals. We validated our findings by quantitative PCR and protein studies.ResultsWe identified eight significantly differentially expressed genes between the mutant and non-mutant fibroblasts, that included downregulated insulin growth factor binding factor protein 5 (IGFBP5) in patient samples. Pathway analysis showed involvement of the ERK/MAPK signaling pathway consistent with previous studies. We found no significant differences in gene expression for Lamin A/C and B-type lamins between the groups. In mutant fibroblasts, RNA-seq confirmed that only the LMNA wild type allele predominately was expressed, and Western Blot showed normal Lamin A/C protein levels.ConclusionsIGFBP5 may contribute in maintaining signaling pathway homeostasis, which may lead to the absence of notable molecular and structural abnormalities in unaffected tissues such as fibroblasts. Compensatory mechanisms from other nuclear membrane proteins were not found. Our results also demonstrate that only one copy of the wild type allele is sufficient for normal levels of Lamin A/C protein to maintain physiological function in an unaffected cell type. This suggests that affected cell types such as cardiac tissues may be more sensitive to haploinsufficiency of Lamin A/C. These results provide insight into the molecular mechanism of disease with a possible explanation for the tissue specificity of LMNA-related dilated cardiomyopathy.

BackgroundCell-free DNA's (cfDNA) use as a biomarker in cancer is challenging due to genetic heterogeneity of malignancies and rarity of tumor-derived molecules. Here we describe and demonstrate a novel machine-learning guided panel design strategy for improving the detection of tumor variants in cfDNA. Using this approach, we first generated a model to classify and score candidate variants for inclusion on a prostate cancer targeted sequencing panel. We then used this panel to screen tumor variants from prostate cancer patients with localized disease in both in silico and hybrid capture settings.MethodsWhole Genome Sequence (WGS) data from 550 prostate tumors was analyzed to build a targeted sequencing panel of single point and small (

Since the pioneering work of Ned Seeman in the early 1980s, the use of the DNA molecule as a construction material experienced a rapid growth and led to the establishment of a new field of science, nowadays called structural DNA nanotechnology. Here, the self-recognition properties of DNA are employed to build micrometer-large molecular objects with nanometer-sized features, thus bridging the nano- to the microscopic world in a programmable fashion. Distinct design strategies and experimental procedures have been developed over the years, enabling the realization of extremely sophisticated structures with a level of control that approaches that of natural macromolecular assemblies. Nevertheless, our understanding of the building process, i.e., what defines the route that goes from the initial mixture of DNA strands to the final intertwined superstructure, is, in some cases, still limited. In this review, we describe the main structural and energetic features of DNA nanoconstructs, from the simple Holliday junction to more complicated DNA architectures, and present the theoretical frameworks that have been formulated until now to explain their self-assembly. Deeper insights into the underlying principles of DNA self-assembly may certainly help us to overcome current experimental challenges and foster the development of original strategies inspired to dissipative and evolutive assembly processes occurring in nature.