Single-cell RNA-seq of Drosophila miranda testis reveals the evolution and trajectory of germline sex chromosome regulation.

Although sex chromosomes have evolved from autosomes, they often have unusual regulatory regimes that are sex- and cell-type-specific such as dosage compensation (DC) and meiotic sex chromosome inactivation (MSCI). The molecular mechanisms and evolutionary forces driving these unique transcriptional programs are critical for genome evolution but have been, in the case of MSCI in Drosophila, subject to continuous debate. Here, we take advantage of the younger sex chromosomes in D. miranda (XR and the neo-X) to infer how former autosomes acquire sex-chromosome-specific regulatory programs using single-cell and bulk RNA sequencing and ribosome profiling, in a comparative evolutionary context. We show that contrary to mammals and worms, the X down-regulation through germline progression is most consistent with the shutdown of DC instead of MSCI, resulting in half gene dosage at the end of meiosis for all 3 X's. Moreover, lowly expressed germline and meiotic genes on the neo-X are ancestrally lowly expressed, instead of acquired suppression after sex linkage. For the young neo-X, DC is incomplete across all tissue and cell types and this dosage imbalance is rescued by contributions from Y-linked gametologs which produce transcripts that are translated to compensate both gene and protein dosage. We find an excess of previously autosomal testis genes becoming Y-specific, showing that the neo-Y and its masculinization likely resolve sexual antagonism. Multicopy neo-sex genes are predominantly expressed during meiotic stages of spermatogenesis, consistent with their amplification being driven to interfere with mendelian segregation. Altogether, this study reveals germline regulation of evolving sex chromosomes and elucidates the consequences these unique regulatory mechanisms have on the evolution of sex chromosome architecture.

Multi-locus genomic signatures of local adaptation to snow across the landscape in California populations of a willow leaf beetle.

Organisms living in mountains contend with extreme climatic conditions, including short growing seasons and long winters with extensive snow cover. Anthropogenic climate change is driving unprecedented, rapid warming of montane regions across the globe, resulting in reduced winter snowpack. Loss of snow as a thermal buffer may have serious consequences for animals overwintering in soil, yet little is known about how variability in snowpack acts as a selective agent in montane ecosystems. Here, we examine genomic variation in California populations of the leaf beetle Chrysomela aeneicollis, an emerging natural model system for understanding how organisms respond to climate change. We used a genotype-environment association approach to identify genomic signatures of local adaptation to microclimate in populations from three montane regions with variable snowpack and a coastal region with no snow. We found that both winter-associated environmental variation and geographical distance contribute to overall genomic variation across the landscape. We identified non-synonymous variation in novel candidate loci associated with cytoskeletal function, ion transport and membrane stability, cellular processes associated with cold tolerance in other insects. These findings provide intriguing evidence that variation in snowpack imposes selective gradients in montane ecosystems.

Chromosome-level reference genomes of two imperiled desert fishes: spikedace (Meda fulgida) and loach minnow (Tiaroga cobitis).

North American minnows (Cypriniformes: Leuciscidae) comprise a diverse taxonomic group, but many members, particularly those inhabiting deserts, face elevated extinction risks. Despite conservation concerns, leuciscids remain under sampled for reference assemblies relative to other groups of freshwater fishes. Here, we present 2 chromosome-scale reference genome assemblies spikedace (Meda fulgida) and loach minnow (Tiaroga cobitis) using PacBio, Illumina and Omni-C technologies. The complete assembly for spikedace was 882.1 Mb in total length comprised of 83 scaffolds with N50 = 34.8 Mb, L50 = 11, N75 = 32.3 Mb, and L75 = 18. The complete assembly for loach minnow was 1.3 Gb in total length comprised of 550 scaffolds with N50 = 48.6 Mb, L50 = 13, N75 = 42.3 Mb, and L75 = 20. Completeness assessed via Benchmarking Universal Single-Copy Orthologues (BUSCO) metrics using the Actinopterygii BUSCO database showed ∼97% for spikedace and ∼98% for loach minnow of complete BUSCO proportions. Annotation revealed approximately 32.58 and 29.04% of spikedace and loach minnow total genome lengths to be comprised of protein-coding genes, respectively. Comparative genomic analyses of these endangered and co-distributed fishes revealed widespread structural variants, gene family expansions, and evidence of positive selection in both genomes.

Highly-multiplexed and efficient long-amplicon PacBio and Nanopore sequencing of hundreds of full mitochondrial genomes.

BACKGROUND: Mitochondrial genome sequences have become critical to the study of biodiversity. Genome skimming and other short-read based methods are the most common approaches, but they are not well-suited to scale up to multiplexing hundreds of samples. Here, we report on a new approach to sequence hundreds to thousands of complete mitochondrial genomes in parallel using long-amplicon sequencing. We amplified the mitochondrial genome of 677 specimens in two partially overlapping amplicons and implemented an asymmetric PCR-based indexing approach to multiplex 1,159 long amplicons together on a single PacBio SMRT Sequel II cell. We also tested this method on Oxford Nanopore Technologies (ONT) MinION R9.4 to assess if this method could be applied to other long-read technologies. We implemented several optimizations that make this method significantly more efficient than alternative mitochondrial genome sequencing methods. RESULTS: With the PacBio sequencing data we recovered at least one of the two fragments for 96% of samples (~ 80-90%) with mean coverage ~ 1,500x. The ONT data recovered less than 50% of input fragments likely due to low throughput and the design of the Barcoded Universal Primers which were optimized for PacBio sequencing. We compared a single mitochondrial gene alignment to half and full mitochondrial genomes and found, as expected, increased tree support with longer alignments, though whole mitochondrial genomes were not significantly better than half mitochondrial genomes. CONCLUSIONS: This method can effectively capture thousands of long amplicons in a single run and be used to build more robust phylogenies quickly and effectively. We provide several recommendations for future users depending on the evolutionary scale of their system. A natural extension of this method is to collect multi-locus datasets consisting of mitochondrial genomes and several long nuclear loci at once.

A chromosome-scale genome assembly and evaluation of mtDNA variation in the willow leaf beetle Chrysomela aeneicollis.

The leaf beetle Chrysomela aeneicollis has a broad geographic range across Western North America but is restricted to cool habitats at high elevations along the west coast. Central California populations occur only at high altitudes (2,700-3,500 m) where they are limited by reduced oxygen supply and recent drought conditions that are associated with climate change. Here, we report a chromosome-scale genome assembly alongside a complete mitochondrial genome and characterize differences among mitochondrial genomes along a latitudinal gradient over which beetles show substantial population structure and adaptation to fluctuating temperatures. Our scaffolded genome assembly consists of 21 linkage groups; one of which we identified as the X chromosome based on female/male whole genome sequencing coverage and orthology with Tribolium castaneum. We identified repetitive sequences in the genome and found them to be broadly distributed across all linkage groups. Using a reference transcriptome, we annotated a total of 12,586 protein-coding genes. We also describe differences in putative secondary structures of mitochondrial RNA molecules, which may generate functional differences important in adaptation to harsh abiotic conditions. We document substitutions at mitochondrial tRNA molecules and substitutions and insertions in the 16S rRNA region that could affect intermolecular interactions with products from the nuclear genome. This first chromosome-level reference genome will enable genomic research in this important model organism for understanding the biological impacts of climate change on montane insects.

A Reference Genome Assembly of Hybrid-Derived California Wild Radish (Raphanus sativus × raphanistrum).

For agriculturally important plants, pollination and herbivory are 2 ecological factors that play into the success of crop yields. Each is also important in natural environments where invasive plants and their effect on species interactions may alter the native ecology. The California Wild Radish (Raphanus sativus × raphanistrum), a hybrid derived from an agriculturally important crop and a nonnative cultivar, is common in California. Remarkably, it has recently replaced wild populations of both progenitor species. Experiments on phenotypic variation for petal color and antiherbivore defenses suggest both pairs of polymorphisms are maintained as a result of pollinator- and herbivore-mediated natural selection. This species provides an opportunity to understand how natural selection shapes the evolution of ecologically important traits when traits are constrained by 2 opposing forces. Here we provide the genome assembly of the California Wild Radish displaying improvement to currently existing genomes for agronomically important crucifers. This genome sequence provides the tools to dissect the genomic architecture of traits related to herbivory and pollination using natural variation in the wild as well as the ability to infer demographic and selective history in the context of hybridization. Study systems like these will improve our understanding and predictions of evolutionary change for correlated traits.

Neo-sex chromosome evolution shapes sex-dependent asymmetrical introgression barrier.

Sex chromosomes play a special role in the evolution of reproductive barriers between species. Here we describe conflicting roles of nascent sex chromosomes on patterns of introgression in an experimental hybrid swarm. Drosophila nasuta and Drosophila albomicans are recently diverged, fully fertile sister species that have different sex chromosome systems. The fusion between an autosome (Muller CD) with the ancestral X and Y gave rise to neo-sex chromosomes in D. albomicans, while Muller CD remains unfused in D. nasuta. We found that a large block containing overlapping inversions on the neo-sex chromosome stood out as the strongest barrier to introgression. Intriguingly, the neo-sex chromosome introgression barrier is asymmetrical and sex-dependent. Female hybrids showed significant D. albomicans­biased introgression on Muller CD (neo-X excess), while males showed heterosis with excessive (neo-X, D. nasuta Muller CD) genotypes. We used a population genetic model to dissect the interplay of sex chromosome drive, heterospecific pairing incompatibility between the neo-sex chromosomes and unfused Muller CD, neo-Y disadvantage, and neo-X advantage in generating the observed sex chromosome genotypes in females and males. We show that moderate neo-Y disadvantage and D. albomicans specific meiotic drive are required to observe female-specific D. albomicans­biased introgression in this system, together with pairing incompatibility and neo-X advantage. In conclusion, this hybrid swarm between a young species pair sheds light onto the multifaceted roles of neo-sex chromosomes in a sex-dependent asymmetrical introgression barrier at a species boundary.

Transposable element accumulation drives size differences among polymorphic Y Chromosomes in Drosophila.

Y Chromosomes of many species are gene poor and show low levels of nucleotide variation, yet they often display high amounts of structural diversity. Dobzhansky cataloged several morphologically distinct Y Chromosomes in Drosophila pseudoobscura that differ in size and shape, but the molecular causes of their large size differences are unclear. Here we use cytogenetics and long-read sequencing to study the sequence content of polymorphic Y Chromosomes in D. pseudoobscura We show that Y Chromosomes differ almost twofold in size, ranging from 30 to 60 Mb. Most of this size difference is caused by a handful of active transposable elements (TEs) that have recently expanded on the largest Y Chromosome, with different elements being responsible for Y expansion on differently sized D. pseudoobscura Y's. We show that Y Chromosomes differ in their heterochromatin enrichment and expression of Y-enriched TEs, and also influence expression of dozens of autosomal and X-linked genes. The same helitron element that showed the most drastic amplification on the largest Y in D. pseudoobscura independently amplified on a polymorphic large Y Chromosome in Drosophila affinis, suggesting that some TEs are inherently more prone to become deregulated on Y Chromosomes.

Dynamics and Impacts of Transposable Element Proliferation in the Drosophila nasuta Species Group Radiation.

Transposable element (TE) mobilization is a constant threat to genome integrity. Eukaryotic organisms have evolved robust defensive mechanisms to suppress their activity, yet TEs can escape suppression and proliferate, creating strong selective pressure for host defense to adapt. This genomic conflict fuels a never-ending arms race that drives the rapid evolution of TEs and recurrent positive selection of genes involved in host defense; the latter has been shown to contribute to postzygotic hybrid incompatibility. However, how TE proliferation impacts genome and regulatory divergence remains poorly understood. Here, we report the highly complete and contiguous (N50 = 33.8-38.0 Mb) genome assemblies of seven closely related Drosophila species that belong to the nasuta species group-a poorly studied group of flies that radiated in the last 2 My. We constructed a high-quality de novo TE library and gathered germline RNA-seq data, which allowed us to comprehensively annotate and compare TE insertion patterns between the species, and infer the evolutionary forces controlling their spread. We find a strong negative association between TE insertion frequency and expression of genes nearby; this likely reflects survivor bias from reduced fitness impact of TEs inserting near lowly expressed, nonessential genes, with limited TE-induced epigenetic silencing. Phylogenetic analyses of insertions of 147 TE families reveal that 53% of them show recent amplification in at least one species. The most highly amplified TE is a nonautonomous DNA element (Drosophila INterspersed Element; DINE) which has gone through multiple bouts of expansions with thousands of full-length copies littered throughout each genome. Across all TEs, we find that TEs expansions are significantly associated with high expression in the expanded species consistent with suppression escape. Thus, whereas horizontal transfer followed by the invasion of a naïve genome has been highlighted to explain the long-term survival of TEs, our analysis suggests that evasion of host suppression of resident TEs is a major strategy to persist over evolutionary times. Altogether, our results shed light on the heterogenous and context-dependent nature in which TEs affect gene regulation and the dynamics of rampant TE proliferation amidst a recently radiated species group.

Origins and evolution of extreme life span in Pacific Ocean rockfishes.

Pacific Ocean rockfishes (genus Sebastes) exhibit extreme variation in life span, with some species being among the most long-lived extant vertebrates. We de novo assembled the genomes of 88 rockfish species and from these identified repeated signatures of positive selection in DNA repair pathways in long-lived taxa and 137 longevity-associated genes with direct effects on life span through insulin signaling and with pleiotropic effects through size and environmental adaptations. A genome-wide screen of structural variation reveals copy number expansions in the immune modulatory butyrophilin gene family in long-lived species. The evolution of different rockfish life histories is coupled to genetic diversity and reshapes the mutational spectrum driving segregating CpG→TpG variants in long-lived species. These analyses highlight the genetic innovations that underlie life history trait adaptations and, in turn, how they shape genomic diversity.

Phosphoglycerate Mutase 1 Activates DNA Damage Repair via Regulation of WIP1 Activity.

The metabolic enzyme phosphoglycerate mutase 1 (PGAM1) is overexpressed in several types of cancer, suggesting an additional function beyond its established role in the glycolytic pathway. We here report that PGAM1 is overexpressed in gliomas where it increases the efficiency of the DNA damage response (DDR) pathway by cytoplasmic binding of WIP1 phosphatase, thereby preventing WIP1 nuclear translocation and subsequent dephosphorylation of the ATM signaling pathway. Silencing of PGAM1 expression in glioma cells consequently decreases formation of γ-H2AX foci, increases apoptosis, and decreases clonogenicity following irradiation (IR) and temozolomide (TMZ) treatment. Furthermore, mice intracranially implanted with PGAM1-knockdown cells have significantly improved survival after treatment with IR and TMZ. These effects are counteracted by exogenous expression of two kinase-dead PGAM1 mutants, H186R and Y92F, indicating an important non-enzymatic function of PGAM1. Our findings identify PGAM1 as a potential therapeutic target in gliomas.

Chromosome-Level Assembly of Drosophila bifasciata Reveals Important Karyotypic Transition of the X Chromosome.

The Drosophila obscura species group is one of the most studied clades of Drosophila and harbors multiple distinct karyotypes. Here we present a de novo genome assembly and annotation of D. bifasciata, a species which represents an important subgroup for which no high-quality chromosome-level genome assembly currently exists. We combined long-read sequencing (Nanopore) and Hi-C scaffolding to achieve a highly contiguous genome assembly approximately 193 Mb in size, with repetitive elements constituting 30.1% of the total length. Drosophila bifasciata harbors four large metacentric chromosomes and the small dot, and our assembly contains each chromosome in a single scaffold, including the highly repetitive pericentromeres, which were largely composed of Jockey and Gypsy transposable elements. We annotated a total of 12,821 protein-coding genes and comparisons of synteny with D. athabasca orthologs show that the large metacentric pericentromeric regions of multiple chromosomes are conserved between these species. Importantly, Muller A (X chromosome) was found to be metacentric in D. bifasciata and the pericentromeric region appears homologous to the pericentromeric region of the fused Muller A-AD (XL and XR) of pseudoobscura/affinis subgroup species. Our finding suggests a metacentric ancestral X fused to a telocentric Muller D and created the large neo-X (Muller A-AD) chromosome ∼15 MYA. We also confirm the fusion of Muller C and D in D. bifasciata and show that it likely involved a centromere-centromere fusion.

Dynamic turnover of centromeres drives karyotype evolution in Drosophila.

Centromeres are the basic unit for chromosome inheritance, but their evolutionary dynamics is poorly understood. We generate high-quality reference genomes for multiple Drosophila obscura group species to reconstruct karyotype evolution. All chromosomes in this lineage were ancestrally telocentric and the creation of metacentric chromosomes in some species was driven by de novo seeding of new centromeres at ancestrally gene-rich regions, independently of chromosomal rearrangements. The emergence of centromeres resulted in a drastic size increase due to repeat accumulation, and dozens of genes previously located in euchromatin are now embedded in pericentromeric heterochromatin. Metacentric chromosomes secondarily became telocentric in the pseudoobscura subgroup through centromere repositioning and a pericentric inversion. The former (peri)centric sequences left behind shrunk dramatically in size after their inactivation, yet contain remnants of their evolutionary past, including increased repeat-content and heterochromatic environment. Centromere movements are accompanied by rapid turnover of the major satellite DNA detected in (peri)centromeric regions.

Zebrafish Sensitivity to Botulinum Neurotoxins.

Botulinum neurotoxins (BoNT) are the most potent known toxins. The mouse LD50 assay is the gold standard for testing BoNT potency, but is not sensitive enough to detect the extremely low levels of neurotoxin that may be present in the serum of sensitive animal species that are showing the effects of BoNT toxicity, such as channel catfish affected by visceral toxicosis of catfish. Since zebrafish are an important animal model for diverse biomedical and basic research, they are readily available and have defined genetic lines that facilitate reproducibility. This makes them attractive for use as an alternative bioassay organism. The utility of zebrafish as a bioassay model organism for BoNT was investigated. The 96 h median immobilizing doses of BoNT/A, BoNT/C, BoNT/E, and BoNT/F for adult male Tübingen strain zebrafish (0.32 g mean weight) at 25 °C were 16.31, 124.6, 4.7, and 0.61 picograms (pg)/fish, respectively. These findings support the use of the zebrafish-based bioassays for evaluating the presence of BoNT/A, BoNT/E, and BoNT/F. Evaluating the basis of the relatively high resistance of zebrafish to BoNT/C and the extreme sensitivity to BoNT/F may reveal unique functional patterns to the action of these neurotoxins.

Zebrafish (Danio rerio) bioassay for visceral toxicosis of catfish and botulinum neurotoxin serotype E.

Visceral toxicosis of catfish (VTC), a sporadic disease of cultured channel catfish (Ictalurus punctatus) often with high mortality, is caused by botulinum neurotoxin serotype E (BoNT/E). Presumptive diagnosis of VTC is based on characteristic clinical signs and lesions, and the production of these signs and mortality after sera from affected fish is administered to sentinel catfish. The diagnosis is confirmed if the toxicity is neutralized with BoNT/E antitoxin. Because small catfish are often unavailable, the utility of adult zebrafish (Danio rerio) was evaluated in BoNT/E and VTC bioassays. Channel catfish and zebrafish susceptibilities were compared using trypsin-activated BoNT/E in a 96-hr trial by intracoelomically administering 0, 1.87, 3.7, 7.5, 15, or 30 pg of toxin per gram of body weight (g-bw) of fish. All of the zebrafish died at the 7.5 pg/g-bw and higher, while the catfish died at the 15 pg/g-bw dose and higher. To test the bioassay, sera from VTC-affected fish or control sera were intracoelomically injected at a dose of 10 µl per zebrafish and 20 µl/g-bw for channel catfish. At 96 hr post-injection, 78% of the zebrafish and 50% of the catfish receiving VTC sera died, while no control fish died. When the VTC sera were preincubated with BoNT/E antitoxin, they became nontoxic to zebrafish. Histology of zebrafish injected with either VTC serum or BoNT/E demonstrated renal necrosis. Normal catfish serum was toxic to larval zebrafish in immersion exposures, abrogating their utility in VTC bioassays. The results demonstrate bioassays using adult zebrafish for detecting BoNT/E and VTC are sensitive and practical.

Determination of the median lethal dose of botulinum serotype E in channel catfish fingerlings.

The median lethal dose of botulinum serotype E in 5.3-g channel catfish Ictalurus punctatus fingerlings was determined. Five tanks (five fish/tank) were assigned to each of the following treatment groups: 70, 50, 35, 25, or 15 pg of purified botulinum serotype E. Fish were injected intracoelomically and observed for 96 h. Administration of the toxin resulted in initial hyperactivity followed by erratic swimming, paresis, and death. The cumulative mortality by treatment group was 100% at 70 pg, 96% at 50 pg, 100% at 35 pg, 88% at 25 pg, and 56% at 15 pg. The median lethal dose was calculated as 13.7 pg/fish (equivalent to a 0.81 median lethal dose for mice Mus musculus) using a logistic regression model. All fish were necropsied; lesions included exophthalmia, ascites, splenic congestion, intussusception of the intestines, congested spleen, and blanching of the intestinal tract. The resultant clinical signs and lesions were similar to those noted in the syndrome of visceral toxicosis of catfish. This study indicates that channel catfish are more sensitive to the effects of botulinum serotype E than laboratory mice, and the signs and lesions of visceral toxicosis of catfish were replicated by injecting catfish with the toxin.