Search

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.

DNA methyltransferases interact with their CpG target sites in the context of variable flanking sequences. We investigated DNA methylation by the human DNMT3B catalytic domain using substrate pools containing CpX target sites in randomized flanking context and identified combined effects of CpG recognition and flanking sequence interaction together with complex contact networks involved in balancing the interaction with different flanking sites. DNA methylation rates were more affected by flanking sequences at non-CpG than at CpG sites. We show that T775 has an essential dynamic role in the catalytic mechanism of DNMT3B. Moreover, we identify six amino acid residues in the DNA-binding interface of DNMT3B (N652, N656, N658, K777, N779, and R823), which are involved in the equalization of methylation rates of CpG sites in favored and disfavored sequence contexts by forming compensatory interactions to the flanking residues including a CpG specific contact to an A at the +1 flanking site. Non-CpG flanking preferences of DNMT3B are highly correlated with non-CpG methylation patterns in human cells. Comparison of the flanking sequence preferences of human and mouse DNMT3B revealed subtle differences suggesting a co-evolution of flanking sequence preferences and cellular DNMT targets.

Genetic sequencing, or DNA sequencing, using the Sanger technique has become widely used in the veterinary diagnostic community. This technology plays a role in verification of PCR results and is used to provide the genetic sequence data needed for phylogenetic analysis, epidemiologic studies, and forensic investigations. The Laboratory Technology Committee of the American Association of Veterinary Laboratory Diagnosticians has prepared guidelines for sample preparation, submission to sequencing facilities or instrumentation, quality assessment of nucleic acid sequence data performed, and for generating basic sequencing data and phylogenetic analysis for diagnostic applications. This guidance is aimed at assisting laboratories in providing consistent, high-quality, and reliable sequence data when using Sanger-based genetic sequencing as a component of their laboratory services.

Single-cell RNA-sequencing (scRNA-Seq) is a compelling approach to directly and simultaneously measure cellular composition and state, which can otherwise only be estimated by applying deconvolution methods to bulk RNA-Seq estimates. However, it has not yet become a widely used tool in population-scale analyses, due to its prohibitively high cost. Here we show that given the same budget, the statistical power of cell-type-specific expression quantitative trait loci (eQTL) mapping can be increased through low-coverage per-cell sequencing of more samples rather than high-coverage sequencing of fewer samples. We use simulations starting from one of the largest available real single-cell RNA-Seq data from 120 individuals to also show that multiple experimental designs with different numbers of samples, cells per sample and reads per cell could have similar statistical power, and choosing an appropriate design can yield large cost savings especially when multiplexed workflows are considered. Finally, we provide a practical approach on selecting cost-effective designs for maximizing cell-type-specific eQTL power which is available in the form of a web tool.

Genomes of all characterized higher eukaryotes harbor examples of transposable element (TE) bursts-the rapid amplification of TE copies throughout a genome. Despite their prevalence, understanding how bursts diversify genomes requires the characterization of actively transposing TEs before insertion sites and structural rearrangements have been obscured by selection acting over evolutionary time. In this study, rice recombinant inbred lines (RILs), generated by crossing a bursting accession and the reference Nipponbare accession, were exploited to characterize the spread of the very active Ping/mPing family through a small population and the resulting impact on genome diversity. Comparative sequence analysis of 272 individuals led to the identification of over 14,000 new insertions of the mPing miniature inverted-repeat transposable element (MITE), with no evidence for silencing of the transposase-encoding Ping element. In addition to new insertions, Ping-encoded transposase was found to preferentially catalyze the excision of mPing loci tightly linked to a second mPing insertion. Similarly, structural variations, including deletion of rice exons or regulatory regions, were enriched for those with break points at one or both ends of linked mPing elements. Taken together, these results indicate that structural variations are generated during a TE burst as transposase catalyzes both the high copy numbers needed to distribute linked elements throughout the genome and the DNA cuts at the TE ends known to dramatically increase the frequency of recombination.

The AlkB family of nonheme Fe(II)/2-oxoglutarate-dependent oxygenases are essential regulators of RNA epigenetics by serving as erasers of one-carbon marks on RNA with release of formaldehyde (FA). Two major human AlkB family members, FTO and ALKBH5, both act as oxidative demethylases of N6-methyladenosine (m6A) but furnish different major products, N6-hydroxymethyladenosine (hm6A) and adenosine (A), respectively. Here we identify foundational mechanistic differences between FTO and ALKBH5 that promote these distinct biochemical outcomes. In contrast to FTO, which follows a traditional oxidative N-demethylation pathway to catalyze conversion of m6A to hm6A with subsequent slow release of A and FA, we find that ALKBH5 catalyzes a direct m6A-to-A transformation with rapid FA release. We identify a catalytic R130/K132/Y139 triad within ALKBH5 that facilitates release of FA via an unprecedented covalent-based demethylation mechanism with direct detection of a covalent intermediate. Importantly, a K132Q mutant furnishes an ALKBH5 enzyme with an m6A demethylation profile that resembles that of FTO, establishing the importance of this residue in the proposed covalent mechanism. Finally, we show that ALKBH5 is an endogenous source of FA in the cell by activity-based sensing of FA fluxes perturbed via ALKBH5 knockdown. This work provides a fundamental biochemical rationale for nonredundant roles of these RNA demethylases beyond different substrate preferences and cellular localization, where m6A demethylation by ALKBH5 versus FTO results in release of FA, an endogenous one-carbon unit but potential genotoxin, at different rates in living systems.

The mammary epithelial cell (MEC) system is a bilayered ductal epithelium of luminal and basal cells, maintained by a lineage of stem and progenitor populations. Here, we used integrated single-cell transcriptomics and chromatin accessibility analysis to reconstruct the cell types of the mouse MEC system and their underlying gene regulatory features in an unbiased manner. We define differentiation states within the secretory type of luminal cells, which forms a continuous spectrum of general luminal progenitor and lactation-committed progenitor cells. By integrating single-cell transcriptomics and chromatin accessibility landscapes, we identify cis- and trans-regulatory elements that are differentially activated in the specific epithelial cell types and our newly defined luminal differentiation states. Our work provides a resource to reveal cis/trans-regulatory elements associated with MEC identity and differentiation that will serve as a reference to determine how the chromatin accessibility landscape changes during breast cancer.

The RNA polymerase II (Pol II) core promoter is the strategic site of convergence of the signals that lead to the initiation of DNA transcription1-5, but the downstream core promoter in humans has been difficult to understand1-3. Here we analyse the human Pol II core promoter and use machine learning to generate predictive models for the downstream core promoter region (DPR) and the TATA box. We developed a method termed HARPE (high-throughput analysis of randomized promoter elements) to create hundreds of thousands of DPR (or TATA box) variants, each with known transcriptional strength. We then analysed the HARPE data by support vector regression (SVR) to provide comprehensive models for the sequence motifs, and found that the SVR-based approach is more effective than a consensus-based method for predicting transcriptional activity. These results show that the DPR is a functionally important core promoter element that is widely used in human promoters. Notably, there appears to be a duality between the DPR and the TATA box, as many promoters contain one or the other element. More broadly, these findings show that functional DNA motifs can be identified by machine learning analysis of a comprehensive set of sequence variants.

Electrochemical biosensors transduce biochemical events (e.g., DNA hybridization) to electrical signals and can be readily interfaced with electronic instrumentation for portability. Nanostructuring the working electrode enhances sensor performance via augmented effective surface area that increases the capture probability of an analyte. However, increasing the effective surface area via thicker nanostructured electrodes hinders the analyte's permeation into the nanostructured volume and limits its access to deeper electrode surfaces. Here, we use nanoporous gold (np-Au) with various thicknesses and pore morphologies coupled with a methylene blue (MB) reporter-tagged DNA probe for DNA target detection as a model system to study the influence of electrode features on electrochemical sensing performance. Independent of the DNA target concentration, the hybridization current (surrogate for detection sensitivity) increases with the surface enhancement factor (EF), until an EF of ∼5, after which the sensor performance deteriorates. Electrochemical and fluorometric quantification of a desorbed DNA probe suggest that DNA permeation is severely limited for higher EFs. In addition, undesirable capacitive currents disguise the faradaic currents from the MB reporter at larger EFs that require higher square wave voltammetry (SWV) frequencies. Finally, a real-time hybridization study reveals that expanding the effective surface area beyond EFs of ∼5 decreases sensor performance.

CRISPR-Cas12c/d proteins share limited homology with Cas12a and Cas9 bacterial CRISPR RNA (crRNA)-guided nucleases used widely for genome editing and DNA detection. However, Cas12c (C2c3)- and Cas12d (CasY)-catalyzed DNA cleavage and genome editing activities have not been directly observed. We show here that a short-complementarity untranslated RNA (scoutRNA), together with crRNA, is required for Cas12d-catalyzed DNA cutting. The scoutRNA differs in secondary structure from previously described tracrRNAs used by CRISPR-Cas9 and some Cas12 enzymes, and in Cas12d-containing systems, scoutRNA includes a conserved five-nucleotide sequence that is essential for activity. In addition to supporting crRNA-directed DNA recognition, biochemical and cell-based experiments establish scoutRNA as an essential cofactor for Cas12c-catalyzed pre-crRNA maturation. These results define scoutRNA as a third type of transcript encoded by a subset of CRISPR-Cas genomic loci and explain how Cas12c/d systems avoid requirements for host factors including ribonuclease III for bacterial RNA-mediated adaptive immunity.

Massively parallel reporter assays (MPRAs) can simultaneously measure the function of thousands of candidate regulatory sequences (CRSs) in a quantitative manner. In this method, CRSs are cloned upstream of a minimal promoter and reporter gene, alongside a unique barcode, and introduced into cells. If the CRS is a functional regulatory element, it will lead to the transcription of the barcode sequence, which is measured via RNA sequencing and normalized for cellular integration via DNA sequencing of the barcode. This technology has been used to test thousands of sequences and their variants for regulatory activity, to decipher the regulatory code and its evolution, and to develop genetic switches. Lentivirus-based MPRA (lentiMPRA) produces 'in-genome' readouts and enables the use of this technique in hard-to-transfect cells. Here, we provide a detailed protocol for lentiMPRA, along with a user-friendly Nextflow-based computational pipeline-MPRAflow-for quantifying CRS activity from different MPRA designs. The lentiMPRA protocol takes ~2 months, which includes sequencing turnaround time and data processing with MPRAflow.

Blood vascular endothelial cells (BECs) control the immune response by regulating blood flow and immune cell recruitment in lymphoid tissues. However, the diversity of BEC and their origins during immune angiogenesis remain unclear. Here we profile transcriptomes of BEC from peripheral lymph nodes and map phenotypes to the vasculature. We identify multiple subsets, including a medullary venous population whose gene signature predicts a selective role in myeloid cell (vs lymphocyte) recruitment to the medulla, confirmed by videomicroscopy. We define five capillary subsets, including a capillary resident precursor (CRP) that displays stem cell and migratory gene signatures, and contributes to homeostatic BEC turnover and to neogenesis of high endothelium after immunization. Cell alignments show retention of developmental programs along trajectories from CRP to mature venous and arterial populations. Our single cell atlas provides a molecular roadmap of the lymph node blood vasculature and defines subset specialization for leukocyte recruitment and vascular homeostasis.

Many proteins regulate the expression of genes by binding to specific regions encoded in the genome1. Here we introduce a new data set of RNA elements in the human genome that are recognized by RNA-binding proteins (RBPs), generated as part of the Encyclopedia of DNA Elements (ENCODE) project phase III. This class of regulatory elements functions only when transcribed into RNA, as they serve as the binding sites for RBPs that control post-transcriptional processes such as splicing, cleavage and polyadenylation, and the editing, localization, stability and translation of mRNAs. We describe the mapping and characterization of RNA elements recognized by a large collection of human RBPs in K562 and HepG2 cells. Integrative analyses using five assays identify RBP binding sites on RNA and chromatin in vivo, the in vitro binding preferences of RBPs, the function of RBP binding sites and the subcellular localization of RBPs, producing 1,223 replicated data sets for 356 RBPs. We describe the spectrum of RBP binding throughout the transcriptome and the connections between these interactions and various aspects of RNA biology, including RNA stability, splicing regulation and RNA localization. These data expand the catalogue of functional elements encoded in the human genome by the addition of a large set of elements that function at the RNA level by interacting with RBPs.

Sex estimation of skeletons is fundamental to many archaeological studies. Currently, three approaches are available to estimate sex-osteology, genomics, or proteomics, but little is known about the relative reliability of these methods in applied settings. We present matching osteological, shotgun-genomic, and proteomic data to estimate the sex of 55 individuals, each with an independent radiocarbon date between 2,440 and 100 cal BP, from two ancestral Ohlone sites in Central California. Sex estimation was possible in 100% of this burial sample using proteomics, in 91% using genomics, and in 51% using osteology. Agreement between the methods was high, however conflicts did occur. Genomic sex estimates were 100% consistent with proteomic and osteological estimates when DNA reads were above 100,000 total sequences. However, more than half the samples had DNA read numbers below this threshold, producing high rates of conflict with osteological and proteomic data where nine out of twenty conditional DNA sex estimates conflicted with proteomics. While the DNA signal decreased by an order of magnitude in the older burial samples, there was no decrease in proteomic signal. We conclude that proteomics provides an important complement to osteological and shotgun-genomic sex estimation.

MotivationGenetic variation in regulatory elements can alter transcription factor (TF) binding by mutating a TF binding motif, which in turn may affect the activity of the regulatory elements. However, it is unclear which motifs are prone to impact transcriptional regulation if mutated. Current motif analysis tools either prioritize TFs based on motif enrichment without linking to a function or are limited in their applications due to the assumption of linearity between motifs and their functional effects.ResultsWe present MAGGIE (Motif Alteration Genome-wide to Globally Investigate Elements), a novel method for identifying motifs mediating TF binding and function. By leveraging measurements from diverse genotypes, MAGGIE uses a statistical approach to link mutations of a motif to changes of an epigenomic feature without assuming a linear relationship. We benchmark MAGGIE across various applications using both simulated and biological datasets and demonstrate its improvement in sensitivity and specificity compared with the state-of-the-art motif analysis approaches. We use MAGGIE to gain novel insights into the divergent functions of distinct NF-κB factors in pro-inflammatory macrophages, revealing the association of p65-p50 co-binding with transcriptional activation and the association of p50 binding lacking p65 with transcriptional repression.Availability and implementationThe Python package for MAGGIE is freely available at https://github.com/zeyang-shen/maggie. The accession number for the NF-κB ChIP-seq data generated for this study is Gene Expression Omnibus: GSE144070.Supplementary informationSupplementary data are available at Bioinformatics online.

Recognition of highly degenerate mammalian splice sites by the core spliceosomal machinery is regulated by several protein factors that predominantly bind exonic splicing motifs. These are postulated to be single-stranded in order to be functional, yet knowledge of secondary structural features that regulate the exposure of exonic splicing motifs across the transcriptome is not currently available. Using transcriptome-wide RNA structural information we show that retained introns in mouse are commonly flanked by a short (≲70 nucleotide), highly base-paired segment upstream and a predominantly single-stranded exonic segment downstream. Splicing assays with select pre-mRNA substrates demonstrate that loops immediately upstream of the introns contain pre-mRNA-specific splicing enhancers, the substitution or hybridization of which impedes splicing. Additionally, the exonic segments flanking the retained introns appeared to be more enriched in a previously identified set of hexameric exonic splicing enhancer (ESE) sequences compared to their spliced counterparts, suggesting that base-pairing in the exonic segments upstream of retained introns could be a means for occlusion of ESEs. The upstream exonic loops of the test substrate promoted recruitment of splicing factors and consequent pre-mRNA structural remodeling, leading up to assembly of the early spliceosome. These results suggest that disruption of exonic stem-loop structures immediately upstream (but not downstream) of the introns regulate alternative splicing events, likely through modulating accessibility of splicing factors.

Meningiomas are the most common primary intracranial tumors, but meningioma metastases are rare. Accordingly, the clinical workup, diagnostic testing, and molecular classification of metastatic meningioma is incompletely understood. Here, we present a case report of multiply recurrent meningioma complicated by liver metastasis. We discuss the patient presentation, imaging findings, and conventional histopathologic characterization of both the intracranial lesion and the metastatic focus. Further, we perform multiplatform molecular profiling, comprised of DNA methylation arrays and RNA-sequencing, of six stereotactically-guided samples from the intracranial meningioma and a single ultrasound-guided liver metastasis biopsy. Our results show that DNA methylation clusters distinguish the liver metastasis from the intracranial meningioma samples, and identify a small focus of hepatocyte contamination with the liver biopsy. Nonetheless, DNA methylation-based classification accurately identifies the liver metastasis as a meningioma with high confidence. We also find that clustering of RNA-sequencing results distinguishes the liver metastasis from the intracranial meningiomas samples, but that differential gene expression classification is confounded by hepatocyte-specific gene expression programs in the liver metastasis. In sum, this case report sheds light on the comparative biology of intracranial and metastatic meningioma. Furthermore, our results support methylation-based classification as a robust method of diagnosing metastatic lesions, underscore the broad utility of DNA methylation array profiling in diagnostic pathology, and caution against the routine use of bulk RNA-sequencing for identifying tumor signatures in heterogeneous metastatic lesions.

Anthocyanin pigments furnish a powerful visual output of the stress and metabolic status of Arabidopsis thaliana plants. Essential for pigment accumulation is TRANSPARENT TESTA19 (TT19), a glutathione S-transferase proposed to bind and stabilize anthocyanins, participating in their vacuolar sequestration, a function conserved across the flowering plants. Here, we report the identification of genetic suppressors that result in anthocyanin accumulation in the absence of TT19. We show that mutations in RDR6, SGS3, or DCL4 suppress the anthocyanin defect of tt19 by pushing carbon towards flavonoid biosynthesis. This effect is not unique to tt19 and extends to at least one other anthocyanin pathway gene mutant. This synergy between mutations in components of the RDR6-SGS3-DCL4 siRNA system and the flavonoid pathway reveals genetic/epigenetic mechanisms regulating metabolic fluxes.

Pediatric cataract is clinically and genetically heterogeneous and is the most common cause of childhood blindness worldwide. In this study, we aimed to identify disease-causing variants in three large British families and one isolated case with autosomal dominant congenital cataract, using whole exome sequencing. We identified four different heterozygous variants, three in the large families and one in the isolated case. Family A, with a novel missense variant (c.178G>C, p.Gly60Arg) in GJA8 with lamellar cataract; family B, with a recurrent variant in GJA8 (c.262C>T, p.Pro88Ser) associated with nuclear cataract; and family C, with a novel variant in GJA3 (c.771dupC, p.Ser258GlnfsTer68) causing a lamellar phenotype. Individual D had a novel variant in GJA3 (c.82G>T, p.Val28Leu) associated with congenital cataract. Each sequence variant was found to co-segregate with disease. Here, we report three novel and one recurrent disease-causing sequence variant in the gap junctional protein encoding genes causing autosomal dominant congenital cataract. Our study further extends the mutation spectrum of these genes and further facilitates clinical diagnosis. A recurrent p.P88S variant in GJA8 causing isolated nuclear cataract provides evidence of further phenotypic heterogeneity associated with this variant.

Members of the conserved Argonaute (Ago) protein family provide defence against invading nucleic acids in eukaryotes in the process of RNA interference. Many prokaryotes also contain Ago proteins that are predicted to be active nucleases; however, their functional activities in host cells remain poorly understood. Here, we characterize the in vitro and in vivo properties of the SeAgo protein from the mesophilic cyanobacterium Synechococcus elongatus. We show that SeAgo is a DNA-guided nuclease preferentially acting on single-stranded DNA targets, with non-specific guide-independent activity observed for double-stranded substrates. The SeAgo gene is steadily expressed in S. elongatus; however, its deletion or overexpression does not change the kinetics of cell growth. When purified from its host cells or from heterologous E. coli, SeAgo is loaded with small guide DNAs whose formation depends on the endonuclease activity of the argonaute protein. SeAgo co-purifies with SSB proteins suggesting that they may also be involved in DNA processing. The SeAgo-associated small DNAs are derived from diverse genomic locations, with certain enrichment for the proposed sites of chromosomal replication initiation and termination, but show no preference for an endogenous plasmid. Therefore, promiscuous genome sampling by SeAgo does not have great effects on cell physiology and plasmid maintenance.

MotivationRNA-sequencing (RNA-seq) enables global identification of RNA-editing sites in biological systems and disease. A salient step in many studies is to identify editing sites that statistically associate with treatment (e.g. case versus control) or covary with biological factors, such as age. However, RNA-seq has technical features that incumbent tests (e.g. t-test and linear regression) do not consider, which can lead to false positives and false negatives.ResultsIn this study, we demonstrate the limitations of currently used tests and introduce the method, RNA-editing tests (REDITs), a suite of tests that employ beta-binomial models to identify differential RNA editing. The tests in REDITs have higher sensitivity than other tests, while also maintaining the type I error (false positive) rate at the nominal level. Applied to the GTEx dataset, we unveil RNA-editing changes associated with age and gender, and differential recoding profiles between brain regions.Availability and implementationREDITs are implemented as functions in R and freely available for download at https://github.com/gxiaolab/REDITs. The repository also provides a code example for leveraging parallelization using multiple cores.

CRISPR-Cas systems are adaptive immune systems that protect bacteria from bacteriophage (phage) infection1. To provide immunity, RNA-guided protein surveillance complexes recognize foreign nucleic acids, triggering their destruction by Cas nucleases2. While the essential requirements for immune activity are well understood, the physiological cues that regulate CRISPR-Cas expression are not. Here, a forward genetic screen identifies a two-component system (KinB-AlgB), previously characterized in the regulation of Pseudomonas aeruginosa alginate biosynthesis3,4, as a regulator of the expression and activity of the P. aeruginosa Type I-F CRISPR-Cas system. Downstream of KinB-AlgB, activators of alginate production AlgU (a σE orthologue) and AlgR repress CRISPR-Cas activity during planktonic and surface-associated growth5. AmrZ, another alginate regulator6, is triggered to repress CRISPR-Cas immunity upon surface association. Pseudomonas phages and plasmids have taken advantage of this regulatory scheme and carry hijacked homologs of AmrZ that repress CRISPR-Cas expression and activity. This suggests that while CRISPR-Cas regulation may be important to limit self-toxicity, endogenous repressive pathways represent a vulnerability for parasite manipulation.

Developmental progression depends on temporally defined changes in gene expression mediated by transient exposure of lineage intermediates to signals in the progenitor niche. To determine whether cell-intrinsic epigenetic mechanisms contribute to signal-induced transcriptional responses, here we manipulate the signalling environment and activity of the histone demethylase LSD1 during differentiation of hESC-gut tube intermediates into pancreatic endocrine cells. We identify a transient requirement for LSD1 in endocrine cell differentiation spanning a short time-window early in pancreas development, a phenotype we reproduced in mice. Examination of enhancer and transcriptome landscapes revealed that LSD1 silences transiently active retinoic acid (RA)-induced enhancers and their target genes. Furthermore, prolonged RA exposure phenocopies LSD1 inhibition, suggesting that LSD1 regulates endocrine cell differentiation by limiting the duration of RA signalling. Our findings identify LSD1-mediated enhancer silencing as a cell-intrinsic epigenetic feedback mechanism by which the duration of the transcriptional response to a developmental signal is limited.

Achromatopsia (ACHM) is an autosomal recessive disease that results in severe visual loss. Symptoms of ACHM include impaired visual acuity, nystagmus, and photoaversion starting from infancy; furthermore, ACHM is associated with bilateral foveal hypoplasia and absent or severely reduced cone photoreceptor function on electroretinography. Here, we performed genetic sequencing in 3 patients from 2 families with ACHM, identifying and functionally characterizing 2 mutations in the activating transcription factor 6 (ATF6) gene. We identified a homozygous deletion covering exons 8-14 of the ATF6 gene from 2 siblings from the same family. In another patient from a different family, we identified a heterozygous deletion covering exons 2 and 3 of the ATF6 gene found in trans with a previously identified ATF6 c.970C>T (p.Arg324Cys) ACHM disease allele. Recombinant ATF6 proteins bearing these exon deletions showed markedly impaired transcriptional activity by qPCR and RNA-Seq analysis compared with WT-ATF6. Finally, RNAscope revealed that ATF6 and the related ATF6B transcripts were expressed in cones as well as in all retinal layers in normal human retina. Overall, our data identify loss-of-function ATF6 disease alleles that cause human foveal disease.

Mucopolysaccharidosis (MPS) is a metabolic storage disorder caused by the deficiency of any lysosomal enzyme required for the breakdown of glycosaminoglycans. A 15-month-old Boston Terrier presented with clinical signs consistent with lysosomal storage disease including corneal opacities, multifocal central nervous system disease and progressively worsening clinical course. Diagnosis was confirmed at necropsy based on histopathologic evaluation of multiple organs demonstrating accumulation of mucopolysaccharides. Whole genome sequencing was used to uncover a frame-shift insertion affecting the alpha-L-iduronidase (IDUA) gene (c.19_20insCGGCCCCC), a mutation confirmed in another Boston Terrier presented 2 years later with a similar clinical picture. Both dogs were homozygous for the IDUA mutation and shared coat colors not recognized as normal for the breed by the American Kennel Club. In contrast, the mutation was not detected in 120 unrelated Boston Terriers as well as 202 dogs from other breeds. Recent inbreeding to select for recessive and unusual coat colors may have concentrated this relatively rare allele in the breed. The identification of the variant enables ante-mortem diagnosis of similar cases and selective breeding to avoid the spread of this disease in the breed. Boston Terriers carrying this variant represent a promising model for MPS I with neurological abnormalities in humans.

Gene expression levels vary across developmental stage, cell type, and region in the brain. Genomic variants also contribute to the variation in expression, and some neuropsychiatric disorder loci may exert their effects through this mechanism. To investigate these relationships, we present BrainVar, a unique resource of paired whole-genome and bulk tissue RNA sequencing from the dorsolateral prefrontal cortex of 176 individuals across prenatal and postnatal development. Here we identify common variants that alter gene expression (expression quantitative trait loci [eQTLs]) constantly across development or predominantly during prenatal or postnatal stages. Both "constant" and "temporal-predominant" eQTLs are enriched for loci associated with neuropsychiatric traits and disorders and colocalize with specific variants. Expression levels of more than 12,000 genes rise or fall in a concerted late-fetal transition, with the transitional genes enriched for cell-type-specific genes and neuropsychiatric risk loci, underscoring the importance of cataloging developmental trajectories in understanding cortical physiology and pathology.

While splicing changes caused by somatic mutations in SF3B1 are known, identifying full-length isoform changes may better elucidate the functional consequences of these mutations. We report nanopore sequencing of full-length cDNA from CLL samples with and without SF3B1 mutation, as well as normal B cell samples, giving a total of 149 million pass reads. We present FLAIR (Full-Length Alternative Isoform analysis of RNA), a computational workflow to identify high-confidence transcripts, perform differential splicing event analysis, and differential isoform analysis. Using nanopore reads, we demonstrate differential 3' splice site changes associated with SF3B1 mutation, agreeing with previous studies. We also observe a strong downregulation of intron retention events associated with SF3B1 mutation. Full-length transcript analysis links multiple alternative splicing events together and allows for better estimates of the abundance of productive versus unproductive isoforms. Our work demonstrates the potential utility of nanopore sequencing for cancer and splicing research.

Higher expression of the human DNA repair enzyme MUTYH has previously been shown to be strongly associated with reduced survival in a panel of 24 human lymphoblastoid cell lines exposed to the alkylating agent N-methyl-N'-nitro-N-nitrosoguanidine (MNNG). The molecular mechanism of MUTYH-enhanced MNNG cytotoxicity is unclear, because MUTYH has a well-established role in the repair of oxidative DNA lesions. Here, we show in mouse embryonic fibroblasts (MEFs) that this MNNG-dependent phenotype does not involve oxidative DNA damage and occurs independently of both O6-methyl guanine adduct cytotoxicity and MUTYH-dependent glycosylase activity. We found that blocking of abasic (AP) sites abolishes higher survival of Mutyh-deficient (Mutyh-/-) MEFs, but this blockade had no additive cytotoxicity in WT MEFs, suggesting the cytotoxicity is due to MUTYH interactions with MNNG-induced AP sites. We found that recombinant mouse MUTYH tightly binds AP sites opposite all four canonical undamaged bases and stimulated apurinic/apyrimidinic endonuclease 1 (APE1)-mediated DNA incision. Consistent with these observations, we found that stable expression of WT, but not catalytically-inactive MUTYH, enhances MNNG cytotoxicity in Mutyh-/- MEFs and that MUTYH expression enhances MNNG-induced genomic strand breaks. Taken together, these results suggest that MUTYH enhances the rapid accumulation of AP-site intermediates by interacting with APE1, implicating MUTYH as a factor that modulates the delicate process of base-excision repair independently of its glycosylase activity.

The probability of point mutations is expected to be highly influenced by the flanking nucleotides that surround them, known as the sequence context. This phenomenon may be mainly attributed to the enzyme that modifies or mutates the genetic material, because most enzymes tend to have specific sequence contexts that dictate their activity. Here, we develop a statistical model that allows for the detection and evaluation of the effects of different sequence contexts on mutation rates from deep population sequencing data. This task is computationally challenging, as the complexity of the model increases exponentially as the context size increases. We established our novel Bayesian method based on sparse model selection methods, with the leading assumption that the number of actual sequence contexts that directly influence mutation rates is minuscule compared with the number of possible sequence contexts. We show that our method is highly accurate on simulated data using pentanucleotide contexts, even when accounting for noisy data. We next analyze empirical population sequencing data from polioviruses and HIV-1 and detect a significant enrichment in sequence contexts associated with deamination by the cellular deaminases ADAR 1/2 and APOBEC3G, respectively. In the current era, where next-generation sequencing data are highly abundant, our approach can be used on any population sequencing data to reveal context-dependent base alterations and may assist in the discovery of novel mutable sites or editing sites.

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

A single enzyme active site that catalyzes multiple reactions is a well-established biochemical theme, but how one nuclease site cleaves both DNA strands of a double helix has not been well understood. In analyzing site-specific DNA cleavage by the mammalian RAG1-RAG2 recombinase, which initiates V(D)J recombination, we find that the active site is reconfigured for the two consecutive reactions and the DNA double helix adopts drastically different structures. For initial nicking of the DNA, a locally unwound and unpaired DNA duplex forms a zipper via alternating interstrand base stacking, rather than melting as generally thought. The second strand cleavage and formation of a hairpin-DNA product requires a global scissor-like movement of protein and DNA, delivering the scissile phosphate into the rearranged active site.

PurposeExome sequencing (ES) has the potential to improve management of congenital anomalies and neurodevelopmental disorders in fetuses, infants, and children. US payers are key stakeholders in patient access to ES. We examined how payers view insurance coverage and clinical utility of pediatric and prenatal ES.MethodsWe employed the framework approach of qualitative research to conduct this study. The study cohort represented 14 payers collectively covering 170,000,000 enrollees.ResultsSeventy-one percent of payers covered pediatric ES despite perceived insufficient evidence because they saw merit in available interventions or in ending the diagnostic odyssey. None covered prenatal ES, because they saw no merit. For pediatric ES, 50% agreed with expanded aspects of clinical utility (e.g., information utility), and 21% considered them sufficient for coverage. For prenatal ES, payers saw little utility until in utero interventions become available.ConclusionThe perceived merit of ES is becoming a factor in payers' coverage for serious diseases with available interventions, even when evidence is perceived insufficient. Payers' views on ES's clinical utility are expanding to include informational utility, aligning with the views of patients and other stakeholders. Our findings inform clinical research, patient advocacy, and policy-making, allowing them to be more relevant to payers.

BackgroundIn cystic fibrosis (CF), the spectrum and frequency of CFTR variants differ by geography and race/ethnicity. CFTR variants in White patients are well-described compared with Latino patients. No studies of CFTR variants have been done in patients with CF in the Dominican Republic or Puerto Rico.MethodsCFTR was sequenced in 61 Dominican Republican patients and 21 Puerto Rican patients with CF and greater than ​​​​60 mmol/L sweat chloride. The spectrum of CFTR variants was identified and the proportion of patients with 0, 1, or 2 CFTR variants identified was determined. The functional effects of identified CFTR variants were investigated using clinical annotation databases and computational prediction tools.ResultsOur study found 10% of Dominican patients had two CFTR variants identified compared with 81% of Puerto Rican patients. No CFTR variants were identified in 69% of Dominican patients and 10% of Puerto Rican patients. In Dominican patients, there were 19 identified CFTR variants, accounting for 25 out of 122 disease alleles (20%). In Puerto Rican patients, there were 16 identified CFTR variants, accounting for 36 out of 42 disease alleles (86%) in Puerto Rican patients. Thirty CFTR variants were identified overall. The most frequent variants for Dominican patients were p.Phe508del and p.Ala559Thr and for Puerto Rican patients were p.Phe508del, p.Arg1066Cys, p.Arg334Trp, and p.I507del.ConclusionsIn this first description of the CFTR variants in patients with CF from the Dominican Republic and Puerto Rico, there was a low detection rate of two CFTR variants after full sequencing with the majority of patients from the Dominican Republic without identified variants.

Genome editing followed by reproductive cloning was previously used to produce two hornless dairy bulls. We crossed one genome-edited dairy bull, homozygous for the dominant PC Celtic POLLED allele, with horned cows (pp) and obtained six heterozygous (PCp) polled calves. The calves had no horns and were otherwise healthy and phenotypically unremarkable. We conducted whole-genome sequencing of all animals using an Illumina HiSeq4000 to achieve ~20× coverage. Bioinformatics analyses revealed the bull was a compound heterozygote, carrying one naturally occurring PC Celtic POLLED allele and an allele containing an additional introgression of the homology-directed repair donor plasmid along with the PC Celtic allele. These alleles segregated in the offspring of this bull, and inheritance of either allele produced polled calves. No other unintended genomic alterations were observed. These data can be used to inform conversations in the scientific community, with regulatory authorities and with the public around 'intentional genomic alterations' and future regulatory actions regarding genome-edited animals.

Nuclear factor κB (NF-κB) RelA is the potent transcriptional activator of inflammatory response genes. We stringently defined a list of direct RelA target genes by integrating physical (chromatin immunoprecipitation sequencing [ChIP-seq]) and functional (RNA sequencing [RNA-seq] in knockouts) datasets. We then dissected each gene's regulatory strategy by testing RelA variants in a primary-cell genetic-complementation assay. All endogenous target genes require RelA to make DNA-base-specific contacts, and none are activatable by the DNA binding domain alone. However, endogenous target genes differ widely in how they employ the two transactivation domains. Through model-aided analysis of the dynamic time-course data, we reveal the gene-specific synergy and redundancy of TA1 and TA2. Given that post-translational modifications control TA1 activity and intrinsic affinity for coactivators determines TA2 activity, the differential TA logics suggests context-dependent versus context-independent control of endogenous RelA-target genes. Although some inflammatory initiators appear to require co-stimulatory TA1 activation, inflammatory resolvers are a part of the NF-κB RelA core response.

DNA from formalin-preserved tissue could unlock a vast repository of genetic information stored in museums worldwide. However, formaldehyde crosslinks proteins and DNA, and prevents ready amplification and DNA sequencing. Formaldehyde acylation also fragments the DNA. Treatment with proteinase K proteolyzes crosslinked proteins to rescue the DNA, though the process is quite slow. To reduce processing time and improve rescue efficiency, we applied the mechanical energy of a vortex fluidic device (VFD) to drive the catalytic activity of proteinase K and recover DNA from American lobster tissue (Homarus americanus) fixed in 3.7% formalin for >1-year. A scan of VFD rotational speeds identified the optimal rotational speed for recovery of PCR-amplifiable DNA and while 500+ base pairs were sequenced, shorter read lengths were more consistently obtained. This VFD-based method also effectively recovered DNA from formalin-preserved samples. The results provide a roadmap for exploring DNA from millions of historical and even extinct species.

CRISPR-based gene drives can spread through wild populations by biasing their own transmission above the 50% value predicted by Mendelian inheritance. These technologies offer population-engineering solutions for combating vector-borne diseases, managing crop pests, and supporting ecosystem conservation efforts. Current technologies raise safety concerns for unintended gene propagation. Herein, we address such concerns by splitting the drive components, Cas9 and gRNAs, into separate alleles to form a trans-complementing split-gene-drive (tGD) and demonstrate its ability to promote super-Mendelian inheritance of the separate transgenes. This dual-component configuration allows for combinatorial transgene optimization and increases safety by restricting escape concerns to experimentation windows. We employ the tGD and a small-molecule-controlled version to investigate the biology of component inheritance and resistant allele formation, and to study the effects of maternal inheritance and impaired homology on efficiency. Lastly, mathematical modeling of tGD spread within populations reveals potential advantages for improving current gene-drive technologies for field population modification.

Sequences encoding Olduvai protein domains (formerly DUF1220) show the greatest human lineage-specific increase in copy number of any coding region in the genome and have been associated, in a dosage-dependent manner, with brain size, cognitive aptitude, autism, and schizophrenia. Tandem intragenic duplications of a three-domain block, termed the Olduvai triplet, in four NBPF genes in the chromosomal 1q21.1-0.2 region, are primarily responsible for the striking human-specific copy number increase. Interestingly, most of the Olduvai triplets are adjacent to, and transcriptionally coregulated with, three human-specific NOTCH2NL genes that have been shown to promote cortical neurogenesis. Until now, the underlying genomic events that drove the Olduvai hyperamplification in humans have remained unexplained. Here, we show that the presence or absence of an alternative first exon of the Olduvai triplet perfectly discriminates between amplified (58/58) and unamplified (0/12) triplets. We provide sequence and breakpoint analyses that suggest the alternative exon was produced by an nonallelic homologous recombination-based mechanism involving the duplicative transposition of an existing Olduvai exon found in the CON3 domain, which typically occurs at the C-terminal end of NBPF genes. We also provide suggestive in vitro evidence that the alternative exon may promote instability through a putative G-quadraplex (pG4)-based mechanism. Lastly, we use single-molecule optical mapping to characterize the intragenic structural variation observed in NBPF genes in 154 unrelated individuals and 52 related individuals from 16 families and show that the presence of pG4-containing Olduvai triplets is strongly correlated with high levels of Olduvai copy number variation. These results suggest that the same driver of genomic instability that allowed the evolutionarily recent, rapid, and extreme human-specific Olduvai expansion remains highly active in the human genome.

Previous efforts to evaluate the detection of human papilloma viral (HPV) DNA in whole saliva as a diagnostic measure for HPV-associated oropharyngeal cancer (HPV-OPC) have not shown sufficient clinical performance. We hypothesize that salivary exosomes are packaged with HPV-associated biomarkers, and efficient enrichment of salivary exosomes through isolation can enhance diagnostic and prognostic performance for HPV-OPC. In this study, an acoustofluidic (the fusion of acoustics and microfluidics) platform was developed to perform size-based isolation of salivary exosomes. These data showed that this platform is capable of consistently isolating exosomes from saliva samples, regardless of viscosity variation and collection method. Compared with the current gold standard, differential centrifugation, droplet digital RT-PCR analysis showed that the average yield of salivary exosomal small RNA from the acoustofluidic platform is 15 times higher. With this high-yield exosome isolation platform, we show that HPV16 DNA could be detected in isolated exosomes from the saliva of HPV-associated OPC patients at 80% concordance with tissues/biopsies positive for HPV16. Overall, these data demonstrated that the acoustofluidic platform can achieve high-purity and high-yield salivary exosome isolation for downstream salivary exosome-based liquid biopsy applications. Additionally, HPV16 DNA sequences in HPV-OPC patients are packaged in salivary exosomes and their isolation will enhance the detection of HPV16 DNA.

Accurate methods for tracking individuals are crucial to the success of fisheries and aquaculture management. Management of migratory salmonid populations, which are important for the health of many economies, ecosystems, and indigenous cultures, is particularly dependent on data gathered from tagged fish. However, the physical tagging methods currently used have many challenges including cost, variable marker retention, and information limited to tagged individuals. Genetic tracking methods combat many of the problems associated with physical tags, but have their own challenges including high cost, potentially difficult marker design, and incompatibility of markers across species. Here we show the feasibility of a new genotyping method for parent-based tagging (PBT), where individuals are tracked through the inherent genetic relationships with their parents. We found that Rapture sequencing, a combination of restriction-site associated DNA and capture sequencing, provides sufficient data for parentage assignment. Additionally, the same capture bait set, which targets specific restriction-site associated DNA loci, can be used for both Rainbow Trout Oncorhynchus mykiss and Chinook Salmon Oncorhynchus tshawytscha. We input 248 single nucleotide polymorphisms from 1,121 samples to parentage assignment software and compared parent-offspring relationships of the spawning pairs recorded in a hatchery. Interestingly, our results suggest sperm contamination during hatchery spawning occurred in the production of 14% of offspring, further confirming the need for genetic tagging in accurately tracking individuals. PBT with Rapture successfully assigned progeny to parents with a 98.86% accuracy with sufficient genetic data. Cost for this pilot study was approximately $3 USD per sample. As costs vary based on the number of markers used and individuals sequenced, we expect that when implemented at a large-scale, per sample costs could be further decreased. We conclude that Rapture PBT provides a cost-effective and accurate alternative to the physical coded wire tags, and other genetic-based methods.

Estrogen contributes to females' strong antibody response to microbial vaccines and proneness to autoimmunity, particularly antibody-mediated systemic autoimmunity, in females. We have hypothesized that this is due to estrogen-mediated potentiation of class switch DNA recombination (CSR) and somatic hypermutation (SHM). As we have shown, estrogen boosts AID expression, which is critical for both CSR and SHM, through upregulation of HoxC4, which together with NF-κB critically mediates Aicda (AID gene) promoter activation. We contend here that additional regulation of Aicda expression by estrogen occurs through epigenetic mechanisms. As we have shown, histone deacetylase inhibitors (HDIs) short-chain fatty acid (SCFA) butyrate and propionate as well as the pharmacologic HDI valproic acid upregulate miRNAs that silence AID expression, thereby modulating specific antibody responses in C57BL/6 mice and autoantibody responses in lupus-prone MRL/Fas lpr/lpr mice. Here, using constitutive knockout Esr1 -/- mice and B cells as well as conditional knockout Aicda cre/cre Esr1 flox/flox mice and B cells, we showed that the HDI-mediated downregulation of Aicda expression as well as the maturation of antibody and autoantibody responses is reversed by estrogen and enhanced by deletion of ERα or E2 inhibition. Estrogen's reversion of HDI-mediated inhibition of Aicda and CSR in antibody and autoantibody responses occurred through downregulation of B cell miR-26a, which, as we showed, targets Aicda mRNA 3'UTR. miR-26a was significantly upregulated by HDIs. Accordingly, enforced expression of miR-26a reduced Aicda expression and CSR, while miR-26a-sponges (competitive inhibitors of miR-26a) increased Aicda expression and CSR. Thus, our findings show that estrogen reverses the HDI-mediated downregulation of AID expression and CSR through selective modulation of miR-26a. They also provide mechanistic insights into the immunomodulatory activity of this hormone and a proof-of-principle for using combined ER inhibitor-HDI as a potential therapeutic approach.

A variety of traits are necessary for bacterial colonization of the interior of plant hosts, including well-studied virulence effectors as well as other phenotypes contributing to bacterial growth and survival within the apoplast. High-throughput methods such as transposon sequencing (TnSeq) are powerful tools to identify such genes in bacterial pathogens. However, there is little information as to the distinctiveness of traits required for bacterial colonization of different hosts. Here, we utilize randomly barcoded TnSeq (RB-TnSeq) to identify the genes that contribute to the ability of Pseudomonas syringae strain B728a to grow within common bean (Phaseolus vulgaris), lima bean (Phaseolus lunatus), and pepper (Capsicum annuum); species representing two different plant families. The magnitude of contribution of most genes to apoplastic fitness in each of the plant hosts was similar. However, 50 genes significantly differed in their fitness contributions to growth within these species. These genes encoded proteins in various functional categories including polysaccharide synthesis and transport, amino acid metabolism and transport, cofactor metabolism, and phytotoxin synthesis and transport. Six genes that encoded unannotated, hypothetical proteins also contributed differentially to growth in these hosts. The genetic repertoire of a relatively promiscuous pathogen such as P. syringae may thus be shaped, at least in part, by the conditional contribution of some fitness determinants.

BackgroundMitochondria support critical cellular functions, such as energy production through oxidative phosphorylation, regulation of reactive oxygen species, apoptosis, and calcium homeostasis.ObjectiveGiven the heightened level of cellular activity in patients with asthma, we sought to determine whether mitochondrial DNA (mtDNA) copy number measured in peripheral blood differed between individuals with and without asthma.MethodsWhole genome sequence data was generated as part of the Trans-Omics for Precision Medicine (TOPMed) Program on participants from the Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-ethnicity (SAPPHIRE) and the Study of African Americans, Asthma, Genes, & Environment II (SAGE II). We restricted our analysis to individuals who self-identified as African American (3,651 asthma cases and 1,344 controls). Mitochondrial copy number was estimated using the sequencing read depth ratio for the mitochondrial and nuclear genomes. Respiratory complex expression was assessed using RNA-sequencing.ResultsAverage mitochondrial copy number was significantly higher among individuals with asthma when compared with controls (SAPPHIRE: 218.60 vs. 200.47, P

Signal transducer and activator of transcription-3 (STAT3) is an oncogenic transcription factor implicated in carcinogenesis, tumor progression, and drug resistance in head and neck squamous cell carcinoma (HNSCC). A decoy oligonucleotide targeting STAT3 offers a promising anti-tumor strategy, but achieving targeted tumor delivery of the decoy with systemic administration poses a significant challenge. We previously showed the potential for STAT3 decoy-loaded microbubbles, in conjunction with ultrasound targeted microbubble cavitation (UTMC), to decrease tumor growth in murine squamous cell carcinoma. As a next step towards clinical translation, we sought to determine the anti-tumor efficacy of our STAT3 decoy delivery platform against human HNSCC and the effect of higher STAT3 decoy microbubble loading on tumor cell inhibition. STAT3 decoy was loaded on cationic lipid microbubbles (STAT3-MB) or loaded on liposome-conjugated lipid microbubbles to form STAT3-loaded liposome-microbubble complexes (STAT3-LPX). UTMC treatment efficacy with these two formulations was evaluated in vitro using viability and apoptosis assays in CAL33 (human HNSCC) cells. Anti-cancer efficacy in vivo was performed in a CAL33 tumor murine xenograft model. UTMC with STAT3-MB caused significantly lower CAL33 cell viability compared to UTMC with STAT3-LPX (56.8±8.4% vs 84.5±8.8%, respectively, p

Recent advances in polymerase engineering have made it possible to isolate aptamers from libraries of synthetic genetic polymers (XNAs) with backbone structures that are distinct from those found in nature. However, nearly all of the XNA aptamers produced thus far have been generated against protein targets, raising significant questions about the ability of XNA aptamers to recognize small molecule targets. Here, we report the evolution of an ATP-binding aptamer composed entirely of α-L-threose nucleic acid (TNA). A chemically synthesized version of the best aptamer sequence shows high affinity to ATP and strong specificity against other naturally occurring ribonucleotide triphosphates. Unlike its DNA and RNA counterparts that are susceptible to nuclease digestion, the ATP-binding TNA aptamer exhibits high biological stability against hydrolytic enzymes that rapidly degrade DNA and RNA. Based on these findings, we suggest that TNA aptamers could find widespread use as molecular recognition elements in diagnostic and therapeutic applications that require high biological stability.

DNA methylation maintenance by DNMT1 is an essential process in mammals but molecular mechanisms connecting DNA methylation patterns and enzyme activity remain elusive. Here, we systematically analyzed the specificity of DNMT1, revealing a pronounced influence of the DNA sequences flanking the target CpG site on DNMT1 activity. We determined DNMT1 structures in complex with preferred DNA substrates revealing that DNMT1 employs flanking sequence-dependent base flipping mechanisms, with large structural rearrangements of the DNA correlating with low catalytic activity. Moreover, flanking sequences influence the conformational dynamics of the active site and cofactor binding pocket. Importantly, we show that the flanking sequence preferences of DNMT1 highly correlate with genomic methylation in human and mouse cells, and 5-azacytidine triggered DNA demethylation is more pronounced at CpG sites with flanks disfavored by DNMT1. Overall, our findings uncover the intricate interplay between CpG-flanking sequence, DNMT1-mediated base flipping and the dynamic landscape of DNA methylation.

As members of the plant microbiota, arbuscular mycorrhizal fungi (AMF, Glomeromycotina) symbiotically colonize plant roots. AMF also possess their own microbiota, hosting some uncultivable endobacteria. Ongoing research has revealed the genetics underlying plant responses to colonization by AMF, but the fungal side of the relationship remains in the dark. Here, we sequenced the genome of Gigaspora margarita, a member of the Gigasporaceae in an early diverging group of the Glomeromycotina. In contrast to other AMF, G. margarita may host distinct endobacterial populations and possesses the largest fungal genome so far annotated (773.104 Mbp), with more than 64% transposable elements. Other unique traits of the G. margarita genome include the expansion of genes for inorganic phosphate metabolism, the presence of genes for production of secondary metabolites and a considerable number of potential horizontal gene transfer events. The sequencing of G. margarita genome reveals the importance of its immune system, shedding light on the evolutionary pathways that allowed early diverging fungi to interact with both plants and bacteria.

The functionality of long noncoding RNAs (lncRNAs) is disputed. In general, lncRNAs are under weak selective pressures, suggesting that the majority of lncRNAs may be nonfunctional. However, although some surveys showed negligible phenotypic effects upon lncRNA perturbation, key biological roles were demonstrated for individual lncRNAs. Most lncRNAs with proven functions were implicated in gene expression regulation, in pathways related to cellular pluripotency, differentiation, and organ morphogenesis, suggesting that functional lncRNAs may be more abundant in embryonic development, rather than in adult organs. To test this hypothesis, we perform a multidimensional comparative transcriptomics analysis, across five developmental time points (two embryonic stages, newborn, adult, and aged individuals), four organs (brain, kidney, liver, and testes), and three species (mouse, rat, and chicken). We find that, overwhelmingly, lncRNAs are preferentially expressed in adult and aged testes, consistent with the presence of permissive transcription during spermatogenesis. LncRNAs are often differentially expressed among developmental stages and are less abundant in embryos and newborns compared with adult individuals, in agreement with a requirement for tighter expression control and less tolerance for noisy transcription early in development. For differentially expressed lncRNAs, we find that the patterns of expression variation among developmental stages are generally conserved between mouse and rat. Moreover, lncRNAs expressed above noise levels in somatic organs and during development show higher evolutionary conservation, in particular, at their promoter regions. Thus, we show that functionally constrained lncRNA loci are enriched in developing organs, and we suggest that many of these loci may function in an RNA-independent manner.

The California Arbovirus Surveillance Program was initiated over 50 years ago to track endemic encephalitides and was enhanced in 2000 to include West Nile virus (WNV) infections in humans, mosquitoes, sentinel chickens, dead birds and horses. This comprehensive statewide program is a function of strong partnerships among the California Department of Public Health (CDPH), the University of California, and local vector control and public health agencies. This manuscript summarizes WNV surveillance data in California since WNV was first detected in 2003 in southern California. From 2003 through 2018, 6,909 human cases of WNV disease, inclusive of 326 deaths, were reported to CDPH, as well as 730 asymptomatic WNV infections identified during screening of blood and organ donors. Of these, 4,073 (59.0%) were reported as West Nile neuroinvasive disease. California's WNV disease burden comprised 15% of all cases that were reported to the U.S. Centers for Disease Control and Prevention during this time, more than any other state. Additionally, 1,299 equine WNV cases were identified, along with detections of WNV in 23,322 dead birds, 31,695 mosquito pools, and 7,340 sentinel chickens. Annual enzootic detection of WNV typically preceded detection in humans and prompted enhanced intervention to reduce the risk of WNV transmission. Peak WNV activity occurred from July through October in the Central Valley and southern California. Less than five percent of WNV activity occurred in other regions of the state or outside of this time. WNV continues to be a major threat to public and wild avian health in California, particularly in southern California and the Central Valley during summer and early fall months. Local and state public health partners must continue statewide human and mosquito surveillance and facilitate effective mosquito control and bite prevention measures.