Genome skimming with nanopore sequencing precisely determines global and transposon DNA methylation in vertebrates.

Genome skimming is defined as low-pass sequencing below 0.05× coverage and is typically used for mitochondrial genome recovery and species identification. Long-read nanopore sequencers enable simultaneous reading of both DNA sequence and methylation and can multiplex samples for low-cost genome skimming. Here I present nanopore sequencing as a highly precise platform for global DNA methylation and transposon assessment. At coverage of just 0.001×, or 30 Mb of reads, accuracy is sub-1%. Biological and technical replicates validate high precision. Skimming 40 vertebrate species reveals conserved patterns of global methylation consistent with whole-genome bisulfite sequencing and an average mapping rate >97%. Genome size directly correlates to global DNA methylation, explaining 39% of its variance. Accurate SINE and LINE transposon methylation in both the mouse and primates can be obtained with just 0.0001× coverage, or 3 Mb of reads. Sample multiplexing, field portability, and the low price of this instrument combine to make genome skimming for DNA methylation an accessible method for epigenetic assessment from ecology to epidemiology and for low-resource groups.

De novo sequencing, diploid assembly, and annotation of the black carpenter ant, Camponotus pennsylvanicus, and its symbionts by one person for $1000, using nanopore sequencing.

The black carpenter ant (Camponotus pennsylvanicus) is a pest species found widely throughout North America. From a single individual I used long-read nanopore sequencing to assemble a phased diploid genome of 306 Mb and 60X coverage, with quality assessed by a 97.0% BUSCO score, improving upon other ant assemblies. The mitochondrial genome reveals minor rearrangements from other ants. The reads also allowed assembly of parasitic and symbiont genomes. I include a complete Wolbachia bacterial assembly with a size of 1.2 Mb, as well as a commensal symbiont Blochmannia pennsylvanicus, at 791 kb. DNA methylation and hydroxymethylation were measured at base-pair resolution level from the same reads and confirmed extremely low levels seen in the Formicidae family. There was moderate heterozygosity, with 0.16% of bases being biallelic from the parental haplotypes. Protein prediction yielded 14 415 amino acid sequences with 95.8% BUSCO score and 86% matching to previously known proteins. All assemblies were derived from a single MinION flow cell generating 20 Gb of sequence for a cost of $1047 including consumable reagents. Adding fixed costs for equipment brings the total for an ant-sized genome to less than $5000. All analyses were performed in 1 week on a single desktop computer.

Creating reference animal genomes is typically a large, expensive process. Here I sequenced the genome of the black carpenter ant for only $1000 as a sole researcher in just one week. Along with the nuclear genome, I assembled the mitochondrial genome and two commensal bacteria species living within the ant. Nanopore technology also enabled epigenetic measurements from the same ant and replicated other studies showing very low DNA methylation. The reference genome compared favorably to other ant species in continuity and protein prediction accuracy. This method will allow other low-resource labs to create high quality genome assemblies with a low cost.

Allergen-induced DNA release by the airway epithelium amplifies type 2 immunity.

BACKGROUND: Alternaria alternata and house dust mite exposure evokes IL-33 secretion from the airway epithelium, which functions as an alarmin to stimulate type 2 immunity. Extracellular DNA (eDNA) is also an alarmin that intensifies inflammation in cystic fibrosis, chronic obstructive pulmonary disease, and asthma. OBJECTIVE: We investigated the mechanisms underlying allergen-evoked DNA mobilization and release from the airway epithelium and determined the role of eDNA in type 2 immunity. METHODS: Human bronchial epithelial (hBE) cells were used to characterize allergen-induced DNA mobilization and extracellular release using comet assays to measure DNA fragmentation, Qubit double-stranded DNA assays to measure DNA release, and DNA sequencing to determine eDNA composition. Mice were used to investigate the role of eDNA in type 2 immunity. RESULTS: Alternaria extract rapidly induces mitochondrial and nuclear DNA release from human bronchial epithelial cells, whereas house dust mite extract induces mitochondrial DNA release. Caspase-3 is responsible for nuclear DNA fragmentation and becomes activated after cleavage by furin. Analysis of secreted nuclear DNA showed disproportionally higher amounts of promotor and exon sequences and lower intron and intergenic regions compared to predictions of random DNA fragmentation. In mice, Alternaria-induced type 2 immune responses were blocked by pretreatment with a DNA scavenger. In caspase-3-deficient mice, Alternaria-induced DNA release was suppressed. Furthermore, intranasal administration of mouse genomic DNA with Alternaria amplified secretion of IL-5 and IL-13 into bronchoalveolar lavage fluid while DNA alone had no effect. CONCLUSION: These findings highlight a novel, allergen-induced mechanism of rapid DNA release that amplifies type 2 immunity in airways.

The mitochondrial genome and Epigenome of the Golden lion Tamarin from fecal DNA using Nanopore adaptive sequencing.

BACKGROUND: The golden lion tamarin (Leontopithecus rosalia) is an endangered Platyrrhine primate endemic to the Atlantic coastal forests of Brazil. Despite ongoing conservation efforts, genetic data on this species remains scarce. Complicating factors include limitations on sample collection and a lack of high-quality reference sequences. Here, we used nanopore adaptive sampling to resequence the L. rosalia mitogenome from feces, a sample which can be collected non-invasively. RESULTS: Adaptive sampling doubled the fraction of both host-derived and mitochondrial sequences compared to sequencing without enrichment. 258x coverage of the L. rosalia mitogenome was achieved in a single flow cell by targeting the unfinished genome of the distantly related emperor tamarin (Saguinus imperator) and the mitogenome of the closely related black lion tamarin (Leontopithecus chrysopygus). The L. rosalia mitogenome has a length of 16,597 bp, sharing 99.68% sequence identity with the L. chrysopygus mitogenome. A total of 38 SNPs between them were identified, with the majority being found in the non-coding D-loop region. DNA methylation and hydroxymethylation were directly detected using a neural network model applied to the raw signal from the MinION sequencer. In contrast to prior reports, DNA methylation was negligible in mitochondria in both CpG and non-CpG contexts. Surprisingly, a quarter of the 642 CpG sites exhibited DNA hydroxymethylation greater than 1% and 44 sites were above 5%, with concentration in the 3' side of several coding regions. CONCLUSIONS: Overall, we report a robust new mitogenome assembly for L. rosalia and direct detection of cytosine base modifications in all contexts.

Epigenetic biomarkers indicate islet cell death in xenotransplantation.

BACKGROUND: Xenotransplantation of porcine islets has emerged in recent decades as a potential treatment for type 1 diabetes (T1D). Current methods of detection, indicative of successful engraftment, occur downstream of actual islet death. Epigenetic biomarkers can be detected in circulating cell-free DNA (cfDNA) to provide an earlier indication of graft dysfunction. AIMS: The present study identified a biomarker of islet death using differential methylation of the insulin gene, INS, originating from ß-cells in porcine islets. MATERIALS & METHODS: Pyrosequencing primers specific for porcine INS were designed to quantify hypomethylation along 12 cysteine-guanine dinucleotide (CpG) sites, including three sites in the cyclic adenosine monophosphate (cAMP) response element (CRE) binding protein 2 (CRE2) binding region of the 5' untranslated region (UTR) and nine sites within intron 2. RESULTS: PCR amplification of bisulfite-converted DNA combined with pyrosequencing data support the conclusion that hypomethylated porcine INS is specific to islet origin. CONCLUSION: Moreover, the results of this study indicate a highly specific epigenetic biomarker, capable of detecting a single islet, supporting the measurement of cfDNA as a biomarker for transplanted islet death. Defining the epigenetic characteristics of porcine-derived islets within cfDNA will be crucial to develop a better understanding of graft survival immunology for transplantation.

Early life social and ecological determinants of global DNA methylation in wild spotted hyenas.

Environmental factors early in life can have lasting influence on the development and phenotypes of animals, but the underlying molecular modifications remain poorly understood. We examined cross-sectional associations among early life socioecological factors and global DNA methylation in 293 wild spotted hyenas (Crocuta crocuta) in the Masai Mara National Reserve, Kenya, grouped according to three age classes (cub, subadult and adult). Explanatory variables of interest included annual maternal rank based on outcomes of dyadic agonistic interactions, litter size, wild ungulate prey density and anthropogenic disturbance in the year each hyena was born based on counts of illegal livestock in the Reserve. The dependent variable of interest was global DNA methylation, assessed via the LUminometric Methylation Assay, which provides a percentage methylation value calculated at CCGG sites across the genome. Among cubs, we observed approximately 2.75% higher CCGG methylation in offspring born to high- than low-ranking mothers. Among cubs and subadults, higher anthropogenic disturbance corresponded with greater %CCGG methylation. In both cubs and adults, we found an inverse association between prey density measured before a hyena was 3 months old and %CCGG methylation. Our results suggest that maternal rank, anthropogenic disturbance and prey availability early in life are associated with later life global DNA methylation. Future studies are required to understand the extent to which these DNA methylation patterns relate to adult phenotypes and fitness outcomes.

Functional Conservation in the SIAMESE-RELATED Family of Cyclin-Dependent Kinase Inhibitors in Land Plants.

The best-characterized members of the plant-specific SIAMESE-RELATED (SMR) family of cyclin-dependent kinase inhibitors regulate the transition from the mitotic cell cycle to endoreplication, also known as endoreduplication, an altered version of the cell cycle in which DNA is replicated without cell division. Some other family members are implicated in cell cycle responses to biotic and abiotic stresses. However, the functions of most SMRs remain unknown, and the specific cyclin-dependent kinase complexes inhibited by SMRs are unclear. Here, we demonstrate that a diverse group of SMRs, including an SMR from the bryophyte Physcomitrella patens, can complement an Arabidopsis thaliana siamese (sim) mutant and that both Arabidopsis SIM and P. patens SMR can inhibit CDK activity in vitro. Furthermore, we show that Arabidopsis SIM can bind to and inhibit both CDKA;1 and CDKB1;1. Finally, we show that SMR2 acts to restrict cell proliferation during leaf growth in Arabidopsis and that SIM, SMR1/LGO, and SMR2 play overlapping roles in controlling the transition from cell division to endoreplication during leaf development. These results indicate that differences in SMR function in plant growth and development are primarily due to differences in transcriptional and posttranscriptional regulation, rather than to differences in fundamental biochemical function.

Perinatal bisphenol A exposure promotes dose-dependent alterations of the mouse methylome.

BACKGROUND: Environmental factors during perinatal development may influence developmental plasticity and disease susceptibility via alterations to the epigenome. Developmental exposure to the endocrine active compound, bisphenol A (BPA), has previously been associated with altered methylation at candidate gene loci. Here, we undertake the first genome-wide characterization of DNA methylation profiles in the liver of murine offspring exposed perinatally to multiple doses of BPA through the maternal diet. RESULTS: Using a tiered focusing approach, our strategy proceeds from unbiased broad DNA methylation analysis using methylation-based next generation sequencing technology to in-depth quantitative site-specific CpG methylation determination using the Sequenom EpiTYPER MassARRAY platform to profile liver DNA methylation patterns in offspring maternally exposed to BPA during gestation and lactation to doses ranging from 0 BPA/kg (Ctr), 50 µg BPA/kg (UG), or 50 mg BPA/kg (MG) diet (N = 4 per group). Genome-wide analyses indicate non-monotonic effects of DNA methylation patterns following perinatal exposure to BPA, corroborating previous studies using multiple doses of BPA with non-monotonic outcomes. We observed enrichment of regions of altered methylation (RAMs) within CpG island (CGI) shores, but little evidence of RAM enrichment in CGIs. An analysis of promoter regions identified several hundred novel BPA-associated methylation events, and methylation alterations in the Myh7b and Slc22a12 gene promoters were validated. Using the Comparative Toxicogenomics Database, a number of candidate genes that have previously been associated with BPA-related gene expression changes were identified, and gene set enrichment testing identified epigenetically dysregulated pathways involved in metabolism and stimulus response. CONCLUSIONS: In this study, non-monotonic dose dependent alterations in DNA methylation among BPA-exposed mouse liver samples and their relevant pathways were identified and validated. The comprehensive methylome map presented here provides candidate loci underlying the role of early BPA exposure and later in life health and disease status.

piOxi database: a web resource of germline and somatic tissue piRNAs identified by chemical oxidation.

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs that are highly expressed and extensively studied from the germline. piRNAs associate with PIWI proteins to maintain DNA methylation for transposon silencing and transcriptional gene regulation for genomic stability. Mature germline piRNAs have distinct characteristics including a 24- to 32-nucleotide length and a 2'-O-methylation signature at the 3' end. Although recent studies have identified piRNAs in somatic tissues, they remain poorly characterized. For example, we recently demonstrated notable expression of piRNA in the murine soma, and while overall expression was lower than that of the germline, unique characteristics suggested tissue-specific functions of this class. While currently available databases commonly use length and association with PIWI proteins to identify piRNA, few have included a chemical oxidation method that detects piRNA based on its 3' modification. This method leads to reproducible and rigorous data processing when coupled with next-generation sequencing and bioinformatics analysis. Here, we introduce piOxi DB, a user-friendly web resource that provides a comprehensive analysis of piRNA, generated exclusively through sodium periodate treatment of small RNA. The current version of piOxi DB includes 435 749 germline and 9828 somatic piRNA sequences robustly identified from M. musculus, M. fascicularis and H. sapiens. The database provides species- and tissue-specific data that are further analyzed according to chromosome location and correspondence to gene and repetitive elements. piOxi DB is an informative tool to assist broad research applications in the fields of RNA biology, cancer biology, environmental toxicology and beyond. Database URL:  https://pioxidb.dcmb.med.umich.edu/.

The genome of Przewalski's horse (Equus ferus przewalskii).

The Przewalski's horse (Equus ferus przewalskii) is an endangered equid native to the steppes of central Asia. After becoming extinct in the wild multiple conservation efforts convened to preserve the species, including captive breeding programs, reintroduction and monitoring systems, protected lands, and cloning. Availability of a highly contiguous reference genome is essential to support these continued efforts. We used Oxford Nanopore sequencing to produce a scaffold-level 2.5 Gb nuclear assembly and 16,002 bp mitogenome from a captive Przewalski's mare. All assembly drafts were generated from 111 Gb of sequence from a single PromethION R10.4.1 flow cell. The mitogenome contained 37 genes in the standard mammalian configuration and was 99.63% identical to the domestic horse (Equus caballus). The nuclear assembly, EquPr2, contained 2,146 scaffolds with an N50 of 85.1 Mb, 43X mean depth, and BUSCO quality score of 98.92%. EquPr2 successfully improves upon the existing Przewalski's horse reference genome (Burgud), with 25-fold fewer scaffolds, a 166-fold larger N50, and phased pseudohaplotypes. Modified basecalls revealed 79.5% DNA methylation and 2.1% hydroxymethylation globally. Allele-specific methylation analysis between pseudohaplotypes revealed 226 differentially methylated regions (DMRs) in known imprinted genes and loci not previously reported as imprinted. The heterozygosity rate of 0.165% matches previous estimates for the species and compares favorably to other endangered animals. This improved Przewalski's horse assembly will serve as a valuable resource for conservation efforts and comparative genomics investigations.

The genome of Przewalski's horse (Equus ferus przewalskii).

The Przewalski's horse (Equus ferus przewalskii) is an endangered equid native to the steppes of central Asia. After becoming extinct in the wild, multiple conservation efforts convened to preserve the species including captive breeding programs, reintroduction and monitoring systems, protected lands, and cloning. Availability of a highly contiguous reference genome is essential to support these continued efforts. We used Oxford Nanopore sequencing to produce a scaffold-level 2.5 Gb nuclear assembly and 16,002 bp mitogenome from a captive Przewalski's mare. All assembly drafts were generated from 111 Gb of sequence from a single PromethION R10.4.1 flow cell. The mitogenome contained 37 genes in the standard mammalian configuration and was 99.63% identical to the domestic horse (Equus caballus). The nuclear assembly, EquPr2, contained 2,146 scaffolds with an N50 of 85.1 Mb, 43X mean depth, and BUSCO quality score of 98.92%. EquPr2 successfully improves upon the existing Przewalski's horse reference genome (Burgud), with 25-fold fewer scaffolds, a 166-fold larger N50, and phased pseudohaplotypes. Modified basecalls revealed 79.5% DNA methylation and 2.1% hydroxymethylation globally. Allele-specific methylation analysis between pseudohaplotypes revealed 226 differentially methylated regions (DMRs) in known imprinted genes and loci not previously reported as imprinted. The heterozygosity rate of 0.165% matches previous estimates for the species and compares favorably to other endangered animals. This improved Przewalski's horse assembly will serve as a valuable resource for conservation efforts and comparative genomics investigations.

Phylogenetic and DNA methylation analysis reveal novel regions of variable methylation in the mouse IAP class of transposons.

BACKGROUND: Select retrotransposons in the long terminal repeat (LTR) class exhibit interindividual variation in DNA methylation that is altered by developmental environmental exposures. Yet, neither the full extent of variability at these "metastable epialleles," nor the phylogenetic relationship underlying variable elements is well understood. The murine metastable epialleles, Avy and CabpIAP, result from independent insertions of an intracisternal A particle (IAP) mobile element, and exhibit remarkably similar sequence identity (98.5%). RESULTS: Utilizing the C57BL/6 genome we identified 10802 IAP LTRs overall and a subset of 1388 in a family that includes Avy and CabpIAP. Phylogenetic analysis revealed two duplication and divergence events subdividing this family into three clades. To characterize interindividual variation across clades, liver DNA from 17 isogenic mice was subjected to combined bisulfite and restriction analysis (CoBRA) for 21 separate LTR transposons (7 per clade). The lowest and highest mean methylation values were 59% and 88% respectively, while methylation levels at individual LTRs varied widely, ranging from 9% to 34%. The clade with the most conserved elements had significantly higher mean methylation across LTRs than either of the two diverged clades (p = 0.040 and p = 0.017). Within each mouse, average methylation across all LTRs was not significantly different (71%-74%, p > 0.99). CONCLUSIONS: Combined phylogenetic and DNA methylation analysis allows for the identification of novel regions of variable methylation. This approach increases the number of known metastable epialleles in the mouse, which can serve as biomarkers for environmental modifications to the epigenome.

An Alu-based phylogeny of gibbons (hylobatidae).

Gibbons (Hylobatidae) are small, arboreal apes indigenous to Southeast Asia that diverged from other apes ∼15-18 Ma. Extant lineages radiated rapidly 6-10 Ma and are organized into four genera (Hylobates, Hoolock, Symphalangus, and Nomascus) consisting of 12-19 species. The use of short interspersed elements (SINEs) as phylogenetic markers has seen recent popularity due to several desirable characteristics: the ancestral state of a locus is known to be the absence of an element, rare potentially homoplasious events are relatively easy to resolve, and samples can be quickly and inexpensively genotyped. During radiation of primates, one particular family of SINEs, the Alu family, has proliferated in primate genomes. Nomascus leucogenys (northern white-cheeked gibbon) sequences were analyzed for repetitive content with RepeatMasker using a custom library. The sequences containing Alu elements identified as members of a gibbon-specific subfamily were then compared with orthologous positions in other primate genomes. A primate phylogenetic panel consisting of 18 primate species, including 13 gibbon species representing all four extant genera, was assayed for all loci, and a total of 125 gibbon-specific Alu insertions were identified. The resulting amplification patterns were used to generate a phylogenetic tree. We demonstrate significant support for Symphalangus as the most basal lineage within the family. Our findings also place Nomascus as a derived lineage, sister to Hoolock, with the Nomascus-Hoolock clade sister to Hylobates. Further, our analysis groups N. leucogenys and Nomascus siki as sister taxa to the exclusion of the other Nomascus species assayed. This study represents the first use of SINEs to determine the genus level phylogenetic relationships within the family Hylobatidae. These relationships have been resolved with robust support at most internal nodes, demonstrating the utility of SINE-based phylogenetic analysis. We postulate that hybridization and rapid radiation may have contributed to the complex and contradictory findings of the previous studies. Our findings will aid in the conservation of these threatened primates and inform future studies of the biogeographical history and distribution of modern gibbon species.

Rex1/Zfp42 as an epigenetic regulator for genomic imprinting.

Zfp42/Rex1 (reduced expression gene 1) is a well-known stem-cell marker that has been duplicated from YY1 in the eutherian lineage. In the current study, we characterized the in vivo roles of Rex1 using a mutant mouse line disrupting its transcription. In contrast to the ubiquitous expression of YY1, Rex1 is expressed only during spermatogenesis and early embryogenesis and also in a very limited area of the placenta. Yet, the gene dosage of Rex1 is very critical for the survival of the late-stage embryos and neonates. This delayed phenotypic consequence suggests potential roles for Rex1 in establishing and maintaining unknown epigenetic modifications. Consistently, Rex1-null blastocysts display hypermethylation in the differentially methylated regions (DMRs) of Peg3 and Gnas imprinted domains, which are known to contain YY1 binding sites. Further analyses confirmed in vivo binding of Rex1 only to the unmethylated allele of these two regions. Thus, Rex1 may function as a protector for these DMRs against DNA methylation. Overall, the functional connection of Rex1 to genomic imprinting represents another case where newly made genes have co-evolved with lineage-specific phenomena.

Genome Skimming with Nanopore Sequencing Precisely Determines Global and Transposon DNA Methylation in Vertebrates.

Genome skimming is defined as low-pass sequencing below 0.05X coverage and is typically used for mitochondrial genome recovery and species identification. Long read nanopore sequencers enable simultaneous reading of both DNA sequence and methylation and can multiplex samples for low-cost genome skimming. Here I present nanopore sequencing as a highly precise platform for global DNA methylation and transposon assessment. At coverage of just 0.001X, or 30 Mb of reads, accuracy is sub-1%. Biological and technical replicates validate high precision. Skimming 40 vertebrate species reveals conserved patterns of global methylation consistent with whole genome bisulfite sequencing and an average mapping rate above 97%. Genome size directly correlates to global DNA methylation, explaining 44% of its variance. Accurate SINE and LINE transposon methylation in both mouse and primates can be obtained with just 0.0001X coverage, or 3 Mb of reads. Sample multiplexing, field portability, and the low price of this instrument combine to make genome skimming for DNA methylation an accessible method for epigenetic assessment from ecology to epidemiology, and by low resource groups.

Intergenerational arsenic exposure on the mouse epigenome and metabolic physiology.

Inorganic arsenic (iAs) is one of the largest toxic exposures to impact humanity worldwide. Exposure to iAs during pregnancy may disrupt the proper remodeling of the epigenome of F1 developing offspring and potentially their F2 grand-offspring via disruption of fetal primordial germ cells (PGCs). There is a limited understanding between the correlation of disease phenotype and methylation profile within offspring of both generations and whether it persists to adulthood. Our study aims to understand the intergenerational effects of in utero iAs exposure on the epigenetic profile and onset of disease phenotypes within F1 and F2 adult offspring, despite the lifelong absence of direct arsenic exposure within these generations. We exposed F0 female mice (C57BL6/J) to the following doses of iAs in drinking water 2 weeks before pregnancy until the birth of the F1 offspring: 1, 10, 245, and 2300 ppb. We found sex- and dose-specific changes in weight and body composition that persist from early time to adulthood within both generations. Fasting blood glucose challenge suggests iAs exposure causes dysregulation of glucose metabolism, revealing generational, exposure, and sex-specific differences. Toward understanding the mechanism, genome-wide DNA methylation data highlights exposure-specific patterns in liver, finding dysregulation within genes associated with cancer, T2D, and obesity. We also identified regions containing persistently differentially methylated CpG sites between F1 and F2 generations. Our results indicate the F1 developing embryos and their PGCs, which will result in F2 progeny, retain epigenetic damage established during the prenatal period and are associated with adult metabolic dysfunction.

Chromosome-level, nanopore-only genome and allele-specific DNA methylation of Pallas's cat, Otocolobus manul.

Pallas's cat, or the manul cat (Otocolobus manul), is a small felid native to the grasslands and steppes of central Asia. Population strongholds in Mongolia and China face growing challenges from climate change, habitat fragmentation, poaching, and other sources. These threats, combined with O. manul's zoo collection popularity and value in evolutionary biology, necessitate improvement of species genomic resources. We used standalone nanopore sequencing to assemble a 2.5 Gb, 61-contig nuclear assembly and 17097 bp mitogenome for O. manul. The primary nuclear assembly had 56× sequencing coverage, a contig N50 of 118 Mb, and a 94.7% BUSCO completeness score for Carnivora-specific genes. High genome collinearity within Felidae permitted alignment-based scaffolding onto the fishing cat (Prionailurus viverrinus) reference genome. Manul contigs spanned all 19 felid chromosomes with an inferred total gap length of less than 400 kilobases. Modified basecalling and variant phasing produced an alternate pseudohaplotype assembly and allele-specific DNA methylation calls; 61 differentially methylated regions were identified between haplotypes. Nearest features included classical imprinted genes, non-coding RNAs, and putative novel imprinted loci. The assembled mitogenome successfully resolved existing discordance between Felinae nuclear and mtDNA phylogenies. All assembly drafts were generated from 158 Gb of sequence using seven minION flow cells.

Metagenomic surveillance for bacterial tick-borne pathogens using nanopore adaptive sampling.

Technological and computational advancements in the fields of genomics and bioinformatics are providing exciting new opportunities for pathogen discovery and genomic surveillance. In particular, single-molecule nucleotide sequence data originating from Oxford Nanopore Technologies (ONT) sequencing platforms can be bioinformatically leveraged, in real-time, for enhanced biosurveillance of a vast array of zoonoses. The recently released nanopore adaptive sampling (NAS) strategy facilitates immediate mapping of individual nucleotide molecules to a given reference as each molecule is being sequenced. User-defined thresholds then allow for the retention or rejection of specific molecules, informed by the real-time reference mapping results, as they are physically passing through a given sequencing nanopore. Here, we show how NAS can be used to selectively sequence DNA of multiple bacterial tick-borne pathogens circulating in wild populations of the blacklegged tick vector, Ixodes scapularis.

Nanopore adaptive sampling for targeted mitochondrial genome sequencing and bloodmeal identification in hematophagous insects.

BACKGROUND: Blood-feeding insects are important vectors for an array of zoonotic pathogens. While previous efforts toward generating molecular resources have largely focused on major vectors of global medical and veterinary importance, molecular data across a large number of hematophagous insect taxa remain limited. Advancements in long-read sequencing technologies and associated bioinformatic pipelines provide new opportunities for targeted sequencing of insect mitochondrial (mt) genomes. For engorged hematophagous insects, such technologies can be leveraged for both insect mitogenome genome assembly and identification of vertebrate blood-meal sources. METHODS: We used nanopore adaptive sampling (NAS) to sequence genomic DNA from four species of field-collected, blood-engorged mosquitoes (Aedes and Culex spp.) and one deer fly (Chrysops sp.). NAS was used for bioinformatical enrichment of mtDNA reads of hematophagous insects and potential vertebrate blood-meal hosts using publically available mt genomes as references. We also performed an experimental control to compare results of traditional non-NAS nanopore sequencing to the mt genome enrichment by the NAS method. RESULTS: Complete mitogenomes were assembled and annotated for all five species sequenced with NAS: Aedes trivittatus, Aedes vexans, Culex restuans, Culex territans and the deer fly, Chrysops niger. In comparison to data generated during our non-NAS control experiment, NAS yielded a substantially higher proportion of reference-mapped mtDNA reads, greatly streamlining downstream mitogenome assembly and annotation. The NAS-assembled mitogenomes ranged in length from 15,582 to 16,045 bp, contained between 78.1% and 79.0% A + T content and shared the anticipated arrangement of 13 protein-coding genes, two ribosomal RNAs, and 22 transfer RNAs. Maximum likelihood phylogenies were generated to further characterize each insect species. Additionally, vertebrate blood-meal analysis was successful in three samples sequenced, with mtDNA-based phylogenetic analyses revealing that blood-meal sources for Chrysops niger, Culex restuans and Aedes trivittatus were human, house sparrow (Passer domesticus) and eastern cottontail rabbit (Sylvilagus floridanus), respectively. CONCLUSIONS: Our findings show that NAS has dual utility to simultaneously molecularly identify hematophagous insects and their blood-meal hosts. Moreover, our data indicate NAS can facilitate a wide array of mitogenomic systematic studies through novel 'phylogenetic capture' methods. We conclude that the NAS approach has great potential for broadly improving genomic resources used to identify blood-feeding insects, answer phylogenetic questions and elucidate complex pathways for the transmission of vector-borne pathogens.

Role of the endocannabinoid system in fragile X syndrome: potential mechanisms for benefit from cannabidiol treatment.

Multiple lines of evidence suggest a central role for the endocannabinoid system (ECS) in the neuronal development and cognitive function and in the pathogenesis of fragile X syndrome (FXS). This review describes the ECS, its role in the central nervous system, how it is dysregulated in FXS, and the potential role of cannabidiol as a treatment for FXS. FXS is caused by deficiency or absence of the fragile X messenger ribonucleoprotein 1 (FMR1) protein, FMRP, typically due to the presence of >200 cytosine, guanine, guanine sequence repeats leading to methylation of the FMR1 gene promoter. The absence of FMRP, following FMR1 gene-silencing, disrupts ECS signaling, which has been implicated in FXS pathogenesis. The ECS facilitates synaptic homeostasis and plasticity through the cannabinoid receptor 1, CB1, on presynaptic terminals, resulting in feedback inhibition of neuronal signaling. ECS-mediated feedback inhibition and synaptic plasticity are thought to be disrupted in FXS, leading to overstimulation, desensitization, and internalization of presynaptic CB1 receptors. Cannabidiol may help restore synaptic homeostasis by acting as a negative allosteric modulator of CB1, thereby attenuating the receptor overstimulation, desensitization, and internalization. Moreover, cannabidiol affects DNA methylation, serotonin 5HT1A signal transduction, gamma-aminobutyric acid receptor signaling, and dopamine D2 and D3 receptor signaling, which may contribute to beneficial effects in patients with FXS. Consistent with these proposed mechanisms of action of cannabidiol in FXS, in the CONNECT-FX trial the transdermal cannabidiol gel, ZYN002, was associated with improvements in measures of social avoidance, irritability, and social interaction, particularly in patients who are most affected, showing ≥90% methylation of the FMR1 gene.