Browsing by Author "Weisz, David"

Now showing 1 - 10 of 10
  • Thumbnail Image
    Item
    A Chromosome-length Assembly of the Black Petaltail (Tanypteryx hageni) Dragonfly
    (Oxford University Press, 2023) Tolman, Ethan R; Beatty, Christopher D; Bush, Jonas; Kohli, Manpreet; Moreno, Carlos M; Ware, Jessica L; Weber, K Scott; Khan, Ruqayya; Maheshwari, Chirag; Weisz, David; Dudchenko, Olga; Aiden, Erez Lieberman; Frandsen, Paul B; Center for Theoretical Biological Physics
    We present a chromosome-length genome assembly and annotation of the Black Petaltail dragonfly (Tanypteryx hageni). This habitat specialist diverged from its sister species over 70 million years ago, and separated from the most closely related Odonata with a reference genome 150 million years ago. Using PacBio HiFi reads and Hi-C data for scaffolding we produce one of the most high-quality Odonata genomes to date. A scaffold N50 of 206.6 Mb and a single copy BUSCO score of 96.2% indicate high contiguity and completeness.
  • Thumbnail Image
    Item
    A rapid, low-cost, and highly sensitive SARS-CoV-2 diagnostic based on whole-genome sequencing
    (Public Library of Science, 2023) Adastra, Per A.; Durand, Neva C.; Mitra, Namita; Pulido, Saul Godinez; Mahajan, Ragini; Blackburn, Alyssa; Colaric, Zane L.; Theisen, Joshua W. M.; Weisz, David; Dudchenko, Olga; Gnirke, Andreas; Rao, Suhas S. P.; Kaur, Parwinder; Aiden, Erez Lieberman; Aiden, Aviva Presser; Center for Theoretical Biological Physics
    Early detection of SARS-CoV-2 infection is key to managing the current global pandemic, as evidence shows the virus is most contagious on or before symptom onset. Here, we introduce a low-cost, high-throughput method for diagnosing and studying SARS-CoV-2 infection. Dubbed Pathogen-Oriented Low-Cost Assembly & Re-Sequencing (POLAR), this method amplifies the entirety of the SARS-CoV-2 genome. This contrasts with typical RT-PCR-based diagnostic tests, which amplify only a few loci. To achieve this goal, we combine a SARS-CoV-2 enrichment method developed by the ARTIC Network (https://artic.network/) with short-read DNA sequencing and de novo genome assembly. Using this method, we can reliably (>95% accuracy) detect SARS-CoV-2 at a concentration of 84 genome equivalents per milliliter (GE/mL). The vast majority of diagnostic methods meeting our analytical criteria that are currently authorized for use by the United States Food and Drug Administration with the Coronavirus Disease 2019 (COVID-19) Emergency Use Authorization require higher concentrations of the virus to achieve this degree of sensitivity and specificity. In addition, we can reliably assemble the SARS-CoV-2 genome in the sample, often with no gaps and perfect accuracy given sufficient viral load. The genotypic data in these genome assemblies enable the more effective analysis of disease spread than is possible with an ordinary binary diagnostic. These data can also help identify vaccine and drug targets. Finally, we show that the diagnoses obtained using POLAR of positive and negative clinical nasal mid-turbinate swab samples 100% match those obtained in a clinical diagnostic lab using the Center for Disease Control’s 2019-Novel Coronavirus test. Using POLAR, a single person can manually process 192 samples over an 8-hour experiment at the cost of ~$36 per patient (as of December 7th, 2022), enabling a 24-hour turnaround with sequencing and data analysis time. We anticipate that further testing and refinement will allow greater sensitivity using this approach.
  • Thumbnail Image
    Item
    Chromosome size affects sequence divergence between species through the interplay of recombination and selection
    (Wiley, 2022) Tigano, Anna; Khan, Ruqayya; Omer, Arina D.; Weisz, David; Dudchenko, Olga; Multani, Asha S.; Pathak, Sen; Behringer, Richard R.; Aiden, Erez L.; Fisher, Heidi; MacManes, Matthew D.; Center for Theoretical and Biological Physics
    The structure of the genome shapes the distribution of genetic diversity and sequence divergence. To investigate how the relationship between chromosome size and recombination rate affects sequence divergence between species, we combined empirical analyses and evolutionary simulations. We estimated pairwise sequence divergence among 15 species from three different mammalian clades—Peromyscus rodents, Mus mice, and great apes—from chromosome-level genome assemblies. We found a strong significant negative correlation between chromosome size and sequence divergence in all species comparisons within the Peromyscus and great apes clades but not the Mus clade, suggesting that the dramatic chromosomal rearrangements among Mus species may have masked the ancestral genomic landscape of divergence in many comparisons. Our evolutionary simulations showed that the main factor determining differences in divergence among chromosomes of different sizes is the interplay of recombination rate and selection, with greater variation in larger populations than in smaller ones. In ancestral populations, shorter chromosomes harbor greater nucleotide diversity. As ancestral populations diverge, diversity present at the onset of the split contributes to greater sequence divergence in shorter chromosomes among daughter species. The combination of empirical data and evolutionary simulations revealed that chromosomal rearrangements, demography, and divergence times may also affect the relationship between chromosome size and divergence, thus deepening our understanding of the role of genome structure in the evolution of species divergence.
  • Thumbnail Image
    Item
    Chromosome-length genome assembly and linkage map of a critically endangered Australian bird: the helmeted honeyeater
    (Oxford University Press, 2022) Robledo-Ruiz, Diana A.; Gan, Han Ming; Kaur, Parwinder; Dudchenko, Olga; Weisz, David; Khan, Ruqayya; Aiden, Erez Lieberman; Osipova, Ekaterina; Hiller, Michael; Morales, Hernán E.; Magrath, Michael J.L.; Clarke, Rohan H.; Sunnucks, Paul; Pavlova, Alexandra; Center for Theoretical Biological Physics
    The helmeted honeyeater (Lichenostomus melanops cassidix) is a Critically Endangered bird endemic to Victoria, Australia. To aid its conservation, the population is the subject of genetic rescue. To understand, monitor, and modulate the effects of genetic rescue on the helmeted honeyeater genome, a chromosome-length genome and a high-density linkage map are required.We used a combination of Illumina, Oxford Nanopore, and Hi-C sequencing technologies to assemble a chromosome-length genome of the helmeted honeyeater, comprising 906 scaffolds, with length of 1.1 Gb and scaffold N50 of 63.8 Mb. Annotation comprised 57,181 gene models. Using a pedigree of 257 birds and 53,111 single-nucleotide polymorphisms, we obtained high-density linkage and recombination maps for 25 autosomes and Z chromosome. The total sex-averaged linkage map was 1,347 cM long, with the male map being 6.7% longer than the female map. Recombination maps revealed sexually dimorphic recombination rates (overall higher in males), with average recombination rate of 1.8 cM/Mb. Comparative analyses revealed high synteny of the helmeted honeyeater genome with that of 3 passerine species (e.g., 32 Hi-C scaffolds mapped to 30 zebra finch autosomes and Z chromosome). The genome assembly and linkage map suggest that the helmeted honeyeater exhibits a fission of chromosome 1A into 2 chromosomes relative to zebra finch. PSMC analysis showed a ∼15-fold decline in effective population size to ∼60,000 from mid- to late Pleistocene.The annotated chromosome-length genome and high-density linkage map provide rich resources for evolutionary studies and will be fundamental in guiding conservation efforts for the helmeted honeyeater.
  • Thumbnail Image
    Item
    A Chromosome-Length Reference Genome for the Endangered Pacific Pocket Mouse Reveals Recent Inbreeding in a Historically Large Population
    (Oxford University Press, 2022) Wilder, Aryn P; Dudchenko, Olga; Curry, Caitlin; Korody, Marisa; Turbek, Sheela P; Daly, Mark; Misuraca, Ann; Wang, Gaojianyong; Khan, Ruqayya; Weisz, David; Fronczek, Julie; Aiden, Erez Lieberman; Houck, Marlys L; Shier, Debra M; Ryder, Oliver A; Steiner, Cynthia C; Center for Theoretical Biological Physics
    High-quality reference genomes are fundamental tools for understanding population history, and can provide estimates of genetic and demographic parameters relevant to the conservation of biodiversity. The federally endangered Pacific pocket mouse (PPM), which persists in three small, isolated populations in southern California, is a promising model for studying how demographic history shapes genetic diversity, and how diversity in turn may influence extinction risk. To facilitate these studies in PPM, we combined PacBio HiFi long reads with Omni-C and Hi-C data to generate a de novo genome assembly, and annotated the genome using RNAseq. The assembly comprised 28 chromosome-length scaffolds (N50 = 72.6 MB) and the complete mitochondrial genome, and included a long heterochromatic region on chromosome 18 not represented in the previously available short-read assembly. Heterozygosity was highly variable across the genome of the reference individual, with 18% of windows falling in runs of homozygosity (ROH) >1 MB, and nearly 9% in tracts spanning >5 MB. Yet outside of ROH, heterozygosity was relatively high (0.0027), and historical Ne estimates were large. These patterns of genetic variation suggest recent inbreeding in a formerly large population. Currently the most contiguous assembly for a heteromyid rodent, this reference genome provides insight into the past and recent demographic history of the population, and will be a critical tool for management and future studies of outbreeding depression, inbreeding depression, and genetic load.
  • Thumbnail Image
    Item
    Emx2 underlies the development and evolution of marsupial gliding membranes
    (Springer Nature, 2024) Moreno, Jorge A.; Dudchenko, Olga; Feigin, Charles Y.; Mereby, Sarah A.; Chen, Zhuoxin; Ramos, Raul; Almet, Axel A.; Sen, Harsha; Brack, Benjamin J.; Johnson, Matthew R.; Li, Sha; Wang, Wei; Gaska, Jenna M.; Ploss, Alexander; Weisz, David; Omer, Arina D.; Yao, Weijie; Colaric, Zane; Kaur, Parwinder; Leger, Judy St; Nie, Qing; Mena, Alexandria; Flanagan, Joseph P.; Keller, Greta; Sanger, Thomas; Ostrow, Bruce; Plikus, Maksim V.; Kvon, Evgeny Z.; Aiden, Erez Lieberman; Mallarino, Ricardo; Center for Theoretical Biological Physics
    Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium—an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.
  • Thumbnail Image
    Item
    Improved high quality sand fly assemblies enabled by ultra low input long read sequencing
    (Springer Nature, 2024) Huang, Michelle; Kingan, Sarah; Shoue, Douglas; Nguyen, Oanh; Froenicke, Lutz; Galvin, Brendan; Lambert, Christine; Khan, Ruqayya; Maheshwari, Chirag; Weisz, David; Maslen, Gareth; Davison, Helen; Aiden, Erez Lieberman; Korlach, Jonas; Dudchenko, Olga; McDowell, Mary Ann; Richards, Stephen; Center for Theoretical Biological Physics
    Phlebotomine sand flies are the vectors of leishmaniasis, a neglected tropical disease. High-quality reference genomes are an important tool for understanding the biology and eco-evolutionary dynamics underpinning disease epidemiology. Previous leishmaniasis vector reference sequences were limited by sequencing technologies available at the time and inadequate for high-resolution genomic inquiry. Here, we present updated reference assemblies of two sand flies, Phlebotomus papatasi and Lutzomyia longipalpis. These chromosome-level assemblies were generated using an ultra-low input library protocol, PacBio HiFi long reads, and Hi-C technology. The new P. papatasi reference has a final assembly span of 351.6 Mb and contig and scaffold N50s of 926 kb and 111.8 Mb, respectively. The new Lu. longipalpis reference has a final assembly span of 147.8 Mb and contig and scaffold N50s of 1.09 Mb and 40.6 Mb, respectively. Benchmarking Universal Single-Copy Orthologue (BUSCO) assessments indicated 94.5% and 95.6% complete single copy insecta orthologs for P. papatasi and Lu. longipalpis. These improved assemblies will serve as an invaluable resource for future genomic work on phlebotomine sandflies.
  • Thumbnail Image
    Item
    MCPH1 inhibits Condensin II during interphase by regulating its SMC2-Kleisin interface
    (eLife Sciences Publications Ltd, 2021) Houlard, Martin; Cutts, Erin E.; Shamim, Muhammad S.; Godwin, Jonathan; Weisz, David; Presser Aiden, Aviva; Aiden, Erez Lieberman; Schermelleh, Lothar; Vannini, Alessandro; Nasmyth, Kim; Center for Theoretical Biological Physics
    Dramatic change in chromosomal DNA morphology between interphase and mitosis is a defining features of the eukaryotic cell cycle. Two types of enzymes, namely cohesin and condensin confer the topology of chromosomal DNA by extruding DNA loops. While condensin normally configures chromosomes exclusively during mitosis, cohesin does so during interphase. The processivity of cohesin’s loop extrusion during interphase is limited by a regulatory factor called WAPL, which induces cohesin to dissociate from chromosomes via a mechanism that requires dissociation of its kleisin from the neck of SMC3. We show here that a related mechanism may be responsible for blocking condensin II from acting during interphase. Cells derived from patients affected by microcephaly caused by mutations in the MCPH1 gene undergo premature chromosome condensation. We show that deletion of Mcph1 in mouse embryonic stem cells unleashes an activity of condensin II that triggers formation of compact chromosomes in G1 and G2 phases, accompanied by enhanced mixing of A and B chromatin compartments, and this occurs even in the absence of CDK1 activity. Crucially, inhibition of condensin II by MCPH1 depends on the binding of a short linear motif within MCPH1 to condensin II’s NCAPG2 subunit. MCPH1’s ability to block condensin II’s association with chromatin is abrogated by the fusion of SMC2 with NCAPH2, hence may work by a mechanism similar to cohesin. Remarkably, in the absence of both WAPL and MCPH1, cohesin and condensin II transform chromosomal DNAs of G2 cells into chromosomes with a solenoidal axis.
  • Thumbnail Image
    Item
    Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample
    (Elsevier, 2024) Sandoval-Velasco, Marcela; Dudchenko, Olga; Rodríguez, Juan Antonio; Pérez Estrada, Cynthia; Dehasque, Marianne; Fontsere, Claudia; Mak, Sarah S. T.; Khan, Ruqayya; Contessoto, Vinícius G.; Oliveira Junior, Antonio B.; Kalluchi, Achyuth; Zubillaga Herrera, Bernardo J.; Jeong, Jiyun; Roy, Renata P.; Christopher, Ishawnia; Weisz, David; Omer, Arina D.; Batra, Sanjit S.; Shamim, Muhammad S.; Durand, Neva C.; O’Connell, Brendan; Roca, Alfred L.; Plikus, Maksim V.; Kusliy, Mariya A.; Romanenko, Svetlana A.; Lemskaya, Natalya A.; Serdyukova, Natalya A.; Modina, Svetlana A.; Perelman, Polina L.; Kizilova, Elena A.; Baiborodin, Sergei I.; Rubtsov, Nikolai B.; Machol, Gur; Rath, Krisha; Mahajan, Ragini; Kaur, Parwinder; Gnirke, Andreas; Garcia-Treviño, Isabel; Coke, Rob; Flanagan, Joseph P.; Pletch, Kelcie; Ruiz-Herrera, Aurora; Plotnikov, Valerii; Pavlov, Innokentiy S.; Pavlova, Naryya I.; Protopopov, Albert V.; Di Pierro, Michele; Graphodatsky, Alexander S.; Lander, Eric S.; Rowley, M. Jordan; Wolynes, Peter G.; Onuchic, José N.; Dalén, Love; Marti-Renom, Marc A.; Gilbert, M. Thomas P.; Aiden, Erez Lieberman; Center for Theoretical Biological Physics
    Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth’s death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.
  • Thumbnail Image
    Item
    Topologically Associated Domains Delineate Susceptibility to Somatic Hypermutation
    (Elsevier, 2019) Senigl, Filip; Maman, Yaakov; Dinesh, Ravi K.; Alinikula, Jukka; Seth, Rashu B.; Pecnova, Lubomira; Omer, Arina D.; Rao, Suhas S. P.; Weisz, David; Buerstedde, Jean-Marie; Aiden, Erez Lieberman; Casellas, Rafael; Hejnar, Jiri; Schatz, David G.; Center for Theoretical Biological Physics
    Somatic hypermutation (SHM) introduces point mutations into immunoglobulin (Ig) genes but also causes mutations in other parts of the genome. We have used lentiviral SHM reporter vectors to identify regions of the genome that are susceptible (“hot”) and resistant (“cold”) to SHM, revealing that SHM susceptibility and resistance are often properties of entire topologically associated domains (TADs). Comparison of hot and cold TADs reveals that while levels of transcription are equivalent, hot TADs are enriched for the cohesin loader NIPBL, super-enhancers, markers of paused/stalled RNA polymerase 2, and multiple important B cell transcription factors. We demonstrate that at least some hot TADs contain enhancers that possess SHM targeting activity and that insertion of a strong Ig SHM-targeting element into a cold TAD renders it hot. Our findings lead to a model for SHM susceptibility involving the cooperative action of cis-acting SHM targeting elements and the dynamic and architectural properties of TADs.