Bringing immersive science to undergraduate laboratory courses using CRISPR gene knockouts in frogs and butterflies

Adaptive introgression of a visual preference gene

Visual preferences are important drivers of mate choice and sexual selection, but little is known of how they evolve at the genetic level. In this study, we took advantage of the diversity of bright warning patterns displayed by Heliconius butterflies, which are also used during mate choice. Combining behavioral, population genomic, and expression analyses, we show that two Heliconius species have evolved the same preferences for red patterns by exchanging genetic material through hybridization. Neural expression of regucalcin1 correlates with visual preference across populations, and disruption of regucalcin1 with CRISPR-Cas9 impairs courtship toward conspecific females, providing a direct link between gene and behavior. Our results support a role for hybridization during behavioral evolution and show how visually guided behaviors contributing to adaptation and speciation are encoded within the genome.

Cis-regulatory modes of Ultrabithorax inactivation in butterfly forewings

Hox gene clusters encode transcription factors that drive regional specialization during animal development: for example the Hox factor Ubx is expressed in the insect metathoracic (T3) wing appendages and differentiates them from T2 mesothoracic identities. Hox transcriptional regulation requires silencing activities that prevent spurious activation and regulatory crosstalks in the wrong tissues, but this has seldom been studied in insects other than Drosophila, which shows a derived Hox dislocation into two genomic clusters that disjoined Antennapedia (Antp) and Ultrabithorax (Ubx). Here, we investigated how Ubx is restricted to the hindwing in butterflies, amidst a contiguous Hox cluster. By analysing Hi-C and ATAC-seq data in the butterfly Junonia coenia, we show that a Topologically Associated Domain (TAD) maintains a hindwing-enriched profile of chromatin opening around Ubx. This TAD is bordered by a Boundary Element (BE) that separates it from a region of joined wing activity around the Antp locus. CRISPR mutational perturbation of this BE releases ectopic Ubx expression in forewings, inducing homeotic clones with hindwing identities. Further mutational interrogation of two non-coding RNA encoding regions and one putative cis-regulatory module within the Ubx TAD cause rare homeotic transformations in both directions, indicating the presence of both activating and repressing chromatin features. We also describe a series of spontaneous forewing homeotic phenotypes obtained in Heliconius butterflies, and discuss their possible mutational basis. By leveraging the extensive wing specialization found in butterflies, our initial exploration of Ubx regulation demonstrates the existence of silencing and insulating sequences that prevent its spurious expression in forewings.

Frizzled2 receives WntA signaling during butterfly wing pattern formation

Butterfly color patterns provide visible and biodiverse phenotypic readouts of the patterning processes. Although the secreted ligand WntA has been shown to instruct the color pattern formation in butterflies, its mode of reception remains elusive. Butterfly genomes encode four homologs of the Frizzled-family of Wnt receptors. Here, we show that CRISPR mosaic knockouts of frizzled2 (fz2) phenocopy the color pattern effects of WntA loss of function in multiple nymphalids. Whereas WntA mosaic clones result in intermediate patterns of reduced size, fz2 clones are cell-autonomous, consistent with a morphogen function. Shifts in expression of WntA and fz2 in WntA crispant pupae show that they are under positive and negative feedback, respectively. Fz1 is required for Wnt-independent planar cell polarity in the wing epithelium. Fz3 and Fz4 show phenotypes consistent with Wnt competitive-antagonist functions in vein formation (Fz3 and Fz4), wing margin specification (Fz3), and color patterning in the Discalis and Marginal Band Systems (Fz4). Overall, these data show that the WntA/Frizzled2 morphogen-receptor pair forms a signaling axis that instructs butterfly color patterning and shed light on the functional diversity of insect Frizzled receptors.

Functional genomic tools for emerging model species

Most studies in the field of ecology and evolution aiming to connect genotype to phenotype rarely validate identified loci using functional tools. Recent developments in RNA interference (RNAi) and clustered regularly interspaced palindromic repeats (CRISPR)-Cas genome editing have dramatically increased the feasibility of functional validation. However, these methods come with specific challenges when applied to emerging model organisms, including limited spatial control of gene silencing, low knock-in efficiencies, and low throughput of functional validation. Moreover, many functional studies to date do not recapitulate ecologically relevant variation, and this limits their scope for deeper insights into evolutionary processes. We therefore argue that increased use of gene editing by allelic replacement through homology-directed repair (HDR) would greatly benefit the field of ecology and evolution.

CRISPR in butterflies: An undergraduate lab experience to inactivate wing patterning genes during development

CRISPR is a technique increasingly used in the laboratory for both fundamental and applied research. We designed and implemented a lab experience for undergraduates to carry out CRISPR technology in the lab, and knockout the wing patterning genes optix and WntA in Vanessa cardui butterflies. Students obtained spectacular phenotypic mutants of butterfly wings color and patterns, awakening curiosity about how genomes encode morphology. In addition, students successfully used molecular techniques to genotype and screen wild-type caterpillar larvae and butterflies for CRISPR edits in genes. Student feedback suggests that they experienced a meaningful process of scientific inquiry by carrying out the whole CRISPR workflow process, from the design and delivery of CRISPR components through microinjection of butterfly eggs, the rearing of live animals through their complete life cycle, and molecular and phenotypic analyses of the resulting mutants. We discuss our experience using CRISP genome editing experiments in butterflies to expose students to hands-on research experiences probing gene-to-phenotype relationships in a charismatic and live organism.

Optimization of CRISPR/Cas9-mediated gene disruption in Xenopus laevis using a phenotypic image analysis technique

The CRISPR/Cas9 method has become popular for gene disruption experiments in Xenopus laevis. However, the experimental conditions that influence the efficiency of CRISPR/Cas9 remain unclear. To that end, we developed an image analysis technique for the semi-quantitative evaluation of the pigment phenotype resulting from the disruption of tyrosinase genes in X. laevis using a CRISPR/Cas9 approach, and then examined the effects of varying five experimental parameters (timing of the CRISPR reagent injection into developing embryos; amount of Cas9 mRNA in the injection reagent; total injection volume per embryo; number of injection sites per embryo; and the culture temperature of the injected embryos) on the gene disruption efficiency. The results of this systematic analysis suggest that the highest possible efficiency of target gene disruption can be achieved by injecting a total of 20 nL of the CRISPR reagent containing 1500 pg of Cas9 mRNA or 4 ng of Cas9 protein into two separate locations (10 nL each) of one-cell stage embryos cultured at 22°C. This study also highlights the importance of balancing the experimental parameters for increasing gene disruption efficiency and provides valuable insights into the optimal conditions for applying the CRISPR/Cas9 system to new experimental organisms.

Implementing CRISPR-Cas9 Yeast Practicals into Biology Curricula

CRISPR-Cas9 is a genome-editing technique that has been widely adopted thanks to its simplicity, efficiency, and broad application potential. Due to its advantages and pervasive use, there have been attempts to include this method in the existing curricula for students majoring in various disciplines of biology. In this perspective, we summarize the existing CRISPR-Cas courses that harness a well-established model organism: baker's yeast, Saccharomyces cerevisiae. As an example, we present a detailed description of a fully hands-on, flexible, robust, and cost-efficient practical CRISPR-Cas9 course, where students participate in yeast genome editing at every stage-from the bioinformatic design of single-guide RNA, through molecular cloning and yeast transformation, to the final confirmation of the introduced mutation. Finally, we emphasize that in addition to providing experimental skills and theoretical knowledge, the practical courses on CRISPR-Cas represent ideal platforms for discussing the ethical implications of the democratization of biology.

High-Quality Genome Assembly and Comprehensive Transcriptome of the Painted Lady Butterfly Vanessa cardui

The painted lady butterfly, Vanessa cardui, has the longest migration routes, the widest hostplant diversity, and one of the most complex wing patterns of any insect. Due to minimal culturing requirements, easily characterized wing pattern elements, and technical feasibility of CRISPR/Cas9 genome editing, V. cardui is emerging as a functional genomics model for diverse research programs. Here, we report a high-quality, annotated genome assembly of the V. cardui genome, generated using 84× coverage of PacBio long-read data, which we assembled into 205 contigs with a total length of 425.4 Mb (N50 = 10.3 Mb). The genome was very complete (single-copy complete Benchmarking Universal Single-Copy Orthologs [BUSCO] 97%), with contigs assembled into presumptive chromosomes using synteny analyses. Our annotation used embryonic, larval, and pupal transcriptomes, and 20 transcriptomes across five different wing developmental stages. Gene annotations showed a high level of accuracy and completeness, with 14,437 predicted protein-coding genes. This annotated genome assembly constitutes an important resource for diverse functional genomic studies ranging from the developmental genetic basis of butterfly color pattern, to coevolution with diverse hostplants.

Porcupine/Wntless-dependent trafficking of the conserved WntA ligand in butterflies

Wnt ligands are key signaling molecules in animals, but little is known about the evolutionary dynamics and mode of action of the WntA orthologs, which are not present in the vertebrates or in Drosophila. Here we show that the WntA subfamily evolved at the base of the Bilateria + Cnidaria clade, and conserved the thumb region and Ser209 acylation site present in most other Wnts, suggesting WntA requires the core Wnt secretory pathway. WntA proteins are distinguishable from other Wnts by a synapomorphic Iso/Val/Ala216 amino-acid residue that replaces the otherwise ubiquitous Thr216 position. WntA embryonic expression is conserved between beetles and butterflies, suggesting functionality, but the WntA gene was lost three times within arthropods, in podoplean copepods, in the cyclorrhaphan fly radiation, and in ensiferan crickets and katydids. Finally, CRISPR mosaic knockouts (KOs) of porcupine and wntless phenocopied the pattern-specific effects of WntA KOs in the wings of Vanessa cardui butterflies. These results highlight the molecular conservation of the WntA protein across invertebrates, and imply it functions as a typical Wnt ligand that is acylated and secreted through the Porcupine/Wntless secretory pathway.

Ultrabithorax Is a Micromanager of Hindwing Identity in Butterflies and Moths

The forewings and hindwings of butterflies and moths (Lepidoptera) are differentiated from each other, with segment-specific morphologies and color patterns that mediate a wide range of functions in flight, signaling, and protection. The Hox geneUltrabithorax(Ubx) is a master selector gene that differentiates metathoracic from mesothoracic identities across winged insects, and previous work has shown this role extends to at least some of the color patterns from the butterfly hindwing. Here we used CRISPR targeted mutagenesis to generateUbxloss-of-function somatic mutations in two nymphalid butterflies (Junonia coenia,Vanessa cardui) and a pyralid moth (Plodia interpunctella). The resulting mosaic clones yielded hindwing-to-forewing transformations, showingUbxis necessary for specifying many aspects of hindwing-specific identities, including scale morphologies, color patterns, and wing venation and structure. These homeotic phenotypes showed cell-autonomous, sharp transitions between mutant and non-mutant scales, except for clones that encroached into the border ocelli (eyespots) and resulted in composite and non-autonomous effects on eyespot ring determination. In the pyralid moth, homeotic clones converted the folding and depigmented hindwing into rigid and pigmented composites, affected the wing-coupling frenulum, and induced ectopic scent-scales in male androconia. These data confirmUbxis a master selector of lepidopteran hindwing identity and suggest it acts on many gene regulatory networks involved in wing development and patterning.