An intronic enhancer of Bmp6 underlies evolved tooth gain in sticklebacks

Continuous tooth replacement: what can teleost fish teach us?

Most tooth-bearing non-mammalian vertebrates have the capacity to replace their teeth throughout life. This capacity was lost in mammals, which replace their teeth only once at most. Not surprisingly, continuous tooth replacement has attracted much attention. Classical morphological studies (e.g. to analyse patterns of replacement) are now being complemented by molecular studies that investigate the expression of genes involved in tooth formation. This review focuses on ray-finned fish (actinopterygians), which have teeth often distributed throughout the mouth and pharynx, and more specifically on teleost fish, the largest group of extant vertebrates. First we highlight the diversity in tooth distribution and in tooth replacement patterns. Replacement tooth formation can start from a distinct (usually discontinuous and transient) dental lamina, but also in the absence of a successional lamina, e.g. from the surface epithelium of the oropharynx or from the outer dental epithelium of a predecessor tooth. The relationship of a replacement tooth to its predecessor is closely related to whether replacement is the result of a prepattern or occurs on demand. As replacement teeth do not necessarily have the same molecular signature as first-generation teeth, the question of the actual trigger for tooth replacement is discussed. Much emphasis has been laid in the past on the potential role of epithelial stem cells in initiating tooth replacement. The outcome of such studies has been equivocal, possibly related to the taxa investigated, and the permanent or transient nature of the dental lamina. Alternatively, replacement may result from local proliferation of undifferentiated progenitors, stimulated by hitherto unknown, perhaps mesenchymal, factors. So far, the role of the neurovascular link in continuous tooth replacement has been poorly investigated, despite the presence of a rich vascularisation surrounding actinopterygian (as well as chondrichthyan) teeth and despite a complete arrest of tooth replacement after nerve resection. Lastly, tooth replacement is possibly co-opted as a process to expand the number of teeth in a dentition ontogenetically whilst conserving features of the primary dentition. That neither a dental lamina, nor stem cells appear to be required for tooth replacement places teleosts in an advantageous position as models for tooth regeneration in humans, where the dental lamina regresses and epithelial stem cells are considered lost.

Transcriptome Analysis Elucidates the Potential Key Genes Involved in Rib Development in bmp6-Deficient Silver Carp (Hypophthalmichthys molitrix)

Bone morphogenetic protein 6 (BMP-6) is a constituent of the TGF-β superfamily, known for its ability to stimulate bone and cartilage formation. The investigation of bmp6's involvement in the formation of intermuscular bones in fish has garnered significant attention in recent years. The rib cage is an important skeletal structure that plays a protective function for internal organs in fish. However, there has been limited research conducted on the effects of the bmp6 gene on rib development. Silver carp is one of four major fish in China, favoured for its affordability and tender muscle. Nevertheless, the presence of numerous intermuscular bones in silver carp significantly hinders the advancement of its palatability and suitability for processing. This study showcases the effective utilisation of CRISPR/Cas9 technology for the purpose of disrupting the bmp6 gene in silver carp, leading to the creation of chimeras in the P0 generation, marking the first instance of such an achievement. The chimeras exhibited complete viability, normal appearance, and partial intermuscular bones loss, with approximately 30% of them displaying rib bifurcation or bending. Subsequently, a transcriptome analysis on ribs of P0 chimeras and wild-type silver carp was conducted, leading to the identification of 934 genes exhibiting differential expression, of which 483 were found to be up-regulated and 451 were found to be down-regulated. The results of the KEGG analysis revealed that the "NF-kappa B signalling pathway", "Hippo signalling pathway", "osteoclast differentiation", and "haematopoietic cell lineage" exhibited enrichment and displayed a significant correlation with bone development. The up-regulated genes such as tnfα, fos, and ctgf in pathways may facilitate the proliferation and differentiation of osteoclasts, whereas the down-regulation of genes such as tgfb2 and tgfbr1 in pathways may hinder the formation and specialisation of osteoblasts, ultimately resulting in rib abnormalities. This study presents novel findings on the impact of bmp6 gene deletion on the rib development of silver carp, while simultaneously investigating the previously unexplored molecular mechanisms underlying rib defects in fish.

Modulation of tooth regeneration through opposing responses to Wnt and BMP signals in teleosts

Most vertebrate species undergo tooth replacement throughout adult life. This process is marked by the shedding of existing teeth and the regeneration of tooth organs. However, little is known about the genetic circuitry regulating tooth replacement. Here, we tested whether fish orthologs of genes known to regulate mammalian hair regeneration have effects on tooth replacement. Using two fish species that demonstrate distinct modes of tooth regeneration, threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio), we found that transgenic overexpression of four different genes changed tooth replacement rates in the direction predicted by a hair regeneration model: Wnt10a and Grem2a increased tooth replacement rate, whereas Bmp6 and Dkk2 strongly inhibited tooth formation. Thus, similar to known roles in hair regeneration, Wnt and BMP signals promote and inhibit regeneration, respectively. Regulation of total tooth number was separable from regulation of replacement rates. RNA sequencing of stickleback dental tissue showed that Bmp6 overexpression resulted in an upregulation of Wnt inhibitors. Together, these data support a model in which different epithelial organs, such as teeth and hair, share genetic circuitry driving organ regeneration.

Co-option of an Astacin Metalloprotease Is Associated with an Evolutionarily Novel Feeding Morphology in a Predatory Nematode

Animals consume a wide variety of food sources to adapt to different environments. However, the genetic mechanisms underlying the acquisition of evolutionarily novel feeding morphology remain largely unknown. While the nematode Caenorhabditis elegans feeds on bacteria, the satellite species Pristionchus pacificus exhibits predatory feeding behavior toward other nematodes, which is an evolutionarily novel feeding habit. Here, we found that the astacin metalloprotease Ppa-NAS-6 is required for the predatory killing by P. pacificus. Ppa-nas-6 mutants were defective in predation-associated characteristics, specifically the tooth morphogenesis and tooth movement during predation. Comparison of expression patterns and rescue experiments of nas-6 in P. pacificus and C. elegans suggested that alteration of the spatial expression patterns of NAS-6 may be vital for acquiring predation-related traits. Reporter analysis of the Ppa-nas-6 promoter in C. elegans revealed that the alteration in expression patterns was caused by evolutionary changes in cis- and trans-regulatory elements. This study suggests that the co-option of a metalloprotease is involved in an evolutionarily novel feeding morphology.

Membrane-bound transcription factor LRRC4 inhibits glioblastoma cell motility

Membrane-bound transcription factors (MTFs) have been observed in many types of organisms, such as plants, animals and microorganisms. However, the routes of MTF nuclear translocation are not well understood. Here, we reported that LRRC4 is a novel MTF that translocates to the nucleus as a full-length protein via endoplasmic reticulum-Golgi transport, which is different from the previously described nuclear entry mechanism. A ChIP-seq assay showed that LRRC4 target genes were mainly involved in cell motility. We confirmed that LRRC4 bound to the enhancer element of the RAP1GAP gene to activate its transcription and inhibited glioblastoma cell movement by affecting cell contraction and polarization. Furthermore, atomic force microscopy (AFM) confirmed that LRRC4 or RAP1GAP altered cellular biophysical properties, such as the surface morphology, adhesion force and cell stiffness. Thus, we propose that LRRC4 is an MTF with a novel route of nuclear translocation. Our observations demonstrate that LRRC4-null glioblastoma led to disordered RAP1GAP gene expression, which increased cellular movement. Re-expression of LRRC4 enabled it to suppress tumors, and this is a potential for targeted treatment in glioblastoma.

Molecular Mechanisms Underlying Vertebrate Adaptive Evolution: A Systematic Review

Adaptive evolution is a process in which variation that confers an evolutionary advantage in a specific environmental context arises and is propagated through a population. When investigating this process, researchers have mainly focused on describing advantageous phenotypes or putative advantageous genotypes. A recent increase in molecular data accessibility and technological advances has allowed researchers to go beyond description and to make inferences about the mechanisms underlying adaptive evolution. In this systematic review, we discuss articles from 2016 to 2022 that investigated or reviewed the molecular mechanisms underlying adaptive evolution in vertebrates in response to environmental variation. Regulatory elements within the genome and regulatory proteins involved in either gene expression or cellular pathways have been shown to play key roles in adaptive evolution in response to most of the discussed environmental factors. Gene losses were suggested to be associated with an adaptive response in some contexts. Future adaptive evolution research could benefit from more investigations focused on noncoding regions of the genome, gene regulation mechanisms, and gene losses potentially yielding advantageous phenotypes. Investigating how novel advantageous genotypes are conserved could also contribute to our knowledge of adaptive evolution.

Repeated genetic divergence plays a minor role in repeated phenotypic divergence of lake-stream stickleback

Recent studies have shown that the repeated evolution of similar phenotypes in response to similar ecological conditions (here "parallel evolution") often occurs through mutations in the same genes. However, many previous studies have focused on known candidate genes in a limited number of systems. Thus, the question of how often parallel phenotypic evolution is due to parallel genetic changes remains open. Here, we used quantitative trait locus (QTL) mapping in F2 intercrosses between lake and stream threespine stickleback (Gasterosteus aculeatus) from four independent watersheds on Vancouver Island, Canada to determine whether the same QTL underlie divergence in the same phenotypes across, between, and within watersheds. We find few parallel QTL, even in independent crosses from the same watershed or for phenotypes that have diverged in parallel. These findings suggest that different mutations can lead to similar phenotypes. The low genetic repeatability observed in these lake-stream systems contrasts with the higher genetic repeatability observed in other stickleback systems. We speculate that differences in evolutionary history, gene flow, and/or the strength and direction of selection might explain these differences in genetic parallelism and emphasize that more work is needed to move beyond documenting genetic parallelism to identifying the underlying causes.

Site pleiotropy of a stickleback Bmp6 enhancer

Development and regeneration are orchestrated by gene regulatory networks that operate in part through transcriptional enhancers. Although many enhancers are pleiotropic and are active in multiple tissues, little is known about whether enhancer pleiotropy is due to 1) site pleiotropy, in which individual transcription factor binding sites (TFBS) are required for activity in multiple tissues, or 2) multiple distinct sites that regulate expression in different tissues. Here, we investigated the pleiotropy of an intronic enhancer of the stickleback Bone morphogenetic protein 6 (Bmp6) gene. This enhancer was previously shown to regulate evolved changes in tooth number and tooth regeneration, and is highly pleiotropic, with robust activity in both fins and teeth throughout embryonic, larval, and adult life, and in the heart and kidney in adult fish. We tested the hypothesis that the pleiotropy of this enhancer is due to site pleiotropy of an evolutionarily conserved predicted Foxc1 TFBS. Transgenic analysis and site-directed mutagenesis experiments both deleting and scrambling this predicted Foxc1 TFBS revealed that the binding site is required for enhancer activity in both teeth and fins throughout embryonic, larval, and adult development, and in the heart and kidney in adult fish. Collectively these data support a model where the pleiotropy of this Bmp6 enhancer is due to site pleiotropy and this putative binding site is required for enhancer activity in multiple anatomical sites from the embryo to the adult.

BMP Signaling Pathway in Dentin Development and Diseases

BMP signaling plays an important role in dentin development. BMPs and antagonists regulate odontoblast differentiation and downstream gene expression via canonical Smad and non-canonical Smad signaling pathways. The interaction of BMPs with their receptors leads to the formation of complexes and the transduction of signals to the canonical Smad signaling pathway (for example, BMP ligands, receptors, and Smads) and the non-canonical Smad signaling pathway (for example, MAPKs, p38, Erk, JNK, and PI3K/Akt) to regulate dental mesenchymal stem cell/progenitor proliferation and differentiation during dentin development and homeostasis. Both the canonical Smad and non-canonical Smad signaling pathways converge at transcription factors, such as Dlx3, Osx, Runx2, and others, to promote the differentiation of dental pulp mesenchymal cells into odontoblasts and downregulated gene expressions, such as those of DSPP and DMP1. Dysregulated BMP signaling causes a number of tooth disorders in humans. Mutation or knockout of BMP signaling-associated genes in mice results in dentin defects which enable a better understanding of the BMP signaling networks underlying odontoblast differentiation and dentin formation. This review summarizes the recent advances in our understanding of BMP signaling in odontoblast differentiation and dentin formation. It includes discussion of the expression of BMPs, their receptors, and the implicated downstream genes during dentinogenesis. In addition, the structures of BMPs, BMP receptors, antagonists, and dysregulation of BMP signaling pathways associated with dentin defects are described.

Freshwater Colonization, Adaptation, and Genomic Divergence in Threespine Stickleback

The Threespine Stickleback is ancestrally a marine fish, but many marine populations breed in fresh water (i.e., are anadromous), facilitating their colonization of isolated freshwater habitats a few years after they form. Repeated adaptation to fresh water during at least 10 My and continuing today has led to Threespine Stickleback becoming a premier system to study rapid adaptation. Anadromous and freshwater stickleback breed in sympatry and may hybridize, resulting in introgression of freshwater-adaptive alleles into anadromous populations, where they are maintained at low frequencies as ancient standing genetic variation. Anadromous stickleback have accumulated hundreds of freshwater-adaptive alleles that are disbursed as few loci per marine individual and provide the basis for adaptation when they colonize fresh water. Recent whole-lake experiments in lakes around Cook Inlet, Alaska have revealed how astonishingly rapid and repeatable this process is, with the frequency of 40% of the identified freshwater-adaptive alleles increasing from negligible (∼1%) in the marine founder to ≥50% within ten generations in fresh water, and freshwater phenotypes evolving accordingly. These high rates of genomic and phenotypic evolution imply very intense directional selection on phenotypes of heterozygotes. Sexual recombination rapidly assembles freshwater-adaptive alleles that originated in different founders into multilocus freshwater haplotypes, and regions important for adaptation to freshwater have suppressed recombination that keeps advantageous alleles linked within large haploblocks. These large haploblocks are also older and appear to have accumulated linked advantageous mutations. The contemporary evolution of Threespine Stickleback has provided broadly applicable insights into the mechanisms that facilitate rapid adaptation.

The ectodysplasin-A receptor is a candidate gene for lateral plate number variation in stickleback fish

Variation in lateral plating in stickleback fish represents a classical example of rapid and parallel adaptation in morphology. The underlying genetic architecture involves polymorphism at the ectodysplasin-A gene (EDA). However, lateral plate number is influenced by additional loci that remain poorly characterized. Here, we search for such loci by performing genome-wide differentiation mapping based on pooled whole-genome sequence data from a European stickleback population variable in the extent of lateral plating, while tightly controlling for the phenotypic effect of EDA. This suggests a new candidate locus, the EDA receptor gene (EDAR), for which additional support is obtained by individual-level targeted Sanger sequencing and by comparing allele frequencies among natural populations. Overall, our study illustrates the power of pooled whole-genome sequencing for searching phenotypically relevant loci and opens opportunities for exploring the population genetics and ecological significance of a new candidate locus for stickleback armor evolution.

Evolved Bmp6 enhancer alleles drive spatial shifts in gene expression during tooth development in sticklebacks

Mutations in enhancers have been shown to often underlie natural variation but the evolved differences in enhancer activity can be difficult to identify in vivo. Threespine sticklebacks (Gasterosteus aculeatus) are a robust system for studying enhancer evolution due to abundant natural genetic variation, a diversity of evolved phenotypes between ancestral marine and derived freshwater forms, and the tractability of transgenic techniques. Previous work identified a series of polymorphisms within an intronic enhancer of the Bone morphogenetic protein 6 (Bmp6) gene that are associated with evolved tooth gain, a derived increase in freshwater tooth number that arises late in development. Here, we use a bicistronic reporter construct containing a genetic insulator and a pair of reciprocal two-color transgenic reporter lines to compare enhancer activity of marine and freshwater alleles of this enhancer. In older fish, the two alleles drive partially overlapping expression in both mesenchyme and epithelium of developing teeth, but the freshwater enhancer drives a reduced mesenchymal domain and a larger epithelial domain relative to the marine enhancer. In younger fish, these spatial shifts in enhancer activity are less pronounced. Comparing Bmp6 expression by in situ hybridization in developing teeth of marine and freshwater fish reveals similar evolved spatial shifts in gene expression. Together, these data support a model in which the polymorphisms within this enhancer underlie evolved tooth gain by shifting the spatial expression of Bmp6 during tooth development, and provide a general strategy to identify spatial differences in enhancer activity in vivo.

Threespine Stickleback: A Model System For Evolutionary Genomics

The repeated adaptation of oceanic threespine sticklebacks to fresh water has made it a premier organism to study parallel evolution. These small fish have multiple distinct ecotypes that display a wide range of diverse phenotypic traits. Ecotypes are easily crossed in the laboratory, and families are large and develop quickly enough for quantitative trait locus analyses, positioning the threespine stickleback as a versatile model organism to address a wide range of biological questions. Extensive genomic resources, including linkage maps, a high-quality reference genome, and developmental genetics tools have led to insights into the genomic basis of adaptation and the identification of genomic changes controlling traits in vertebrates. Recently, threespine sticklebacks have been used as a model system to identify the genomic basis of highly complex traits, such as behavior and host-microbiome and host-parasite interactions. We review the latest findings and new avenues of research that have led the threespine stickleback to be considered a supermodel of evolutionary genomics.

Distinct tooth regeneration systems deploy a conserved battery of genes

BACKGROUND

Vertebrate teeth exhibit a wide range of regenerative systems. Many species, including most mammals, reptiles, and amphibians, form replacement teeth at a histologically distinct location called the successional dental lamina, while other species do not employ such a system. Notably, a 'lamina-less' tooth replacement condition is found in a paraphyletic array of ray-finned fishes, such as stickleback, trout, cod, medaka, and bichir. Furthermore, the position, renewal potential, and latency times appear to vary drastically across different vertebrate tooth regeneration systems. The progenitor cells underlying tooth regeneration thus present highly divergent arrangements and potentials. Given the spectrum of regeneration systems present in vertebrates, it is unclear if morphologically divergent tooth regeneration systems deploy an overlapping battery of genes in their naïve dental tissues.

RESULTS

In the present work, we aimed to determine whether or not tooth progenitor epithelia could be composed of a conserved cell type between vertebrate dentitions with divergent regeneration systems. To address this question, we compared the pharyngeal tooth regeneration processes in two ray-finned fishes: zebrafish (Danio rerio) and threespine stickleback (Gasterosteus aculeatus). These two teleost species diverged approximately 250 million years ago and demonstrate some stark differences in dental morphology and regeneration. Here, we find that the naïve successional dental lamina in zebrafish expresses a battery of nine genes (bmpr1aa, bmp6, cd34, gli1, igfbp5a, lgr4, lgr6, nfatc1, and pitx2), while active Wnt signaling and Lef1 expression occur during early morphogenesis stages of tooth development. We also find that, despite the absence of a histologically distinct successional dental lamina in stickleback tooth fields, the same battery of nine genes (Bmpr1a, Bmp6, CD34, Gli1, Igfbp5a, Lgr4, Lgr6, Nfatc1, and Pitx2) are expressed in the basalmost endodermal cell layer, which is the region most closely associated with replacement tooth germs. Like zebrafish, stickleback replacement tooth germs additionally express Lef1 and exhibit active Wnt signaling. Thus, two fish systems that either have an organized successional dental lamina (zebrafish) or lack a morphologically distinct successional dental lamina (sticklebacks) deploy similar genetic programs during tooth regeneration.

CONCLUSIONS

We propose that the expression domains described here delineate a highly conserved "successional dental epithelium" (SDE). Furthermore, a set of orthologous genes is known to mark hair follicle epithelial stem cells in mice, suggesting that regenerative systems in other epithelial appendages may utilize a related epithelial progenitor cell type, despite the highly derived nature of the resulting functional organs.

Grand Challenges in Comparative Tooth Biology

Teeth are a model system for integrating developmental genomics, functional morphology, and evolution. We are at the cusp of being able to address many open issues in comparative tooth biology and we outline several of these newly tractable and exciting research directions. Like never before, technological advances and methodological approaches are allowing us to investigate the developmental machinery of vertebrates and discover both conserved and excitingly novel mechanisms of diversification. Additionally, studies of the great diversity of soft tissues, replacement teeth, and non-trophic functions of teeth are providing new insights into dental diversity. Finally, we highlight several emerging model groups of organisms that are at the forefront of increasing our appreciation of the mechanisms underlying tooth diversification.

An osmolality/salinity-responsive enhancer 1 (OSRE1) in intron 1 promotes salinity induction of tilapia glutamine synthetase

Euryhaline tilapia (Oreochromis mossambicus) are fish that tolerate a wide salinity range from fresh water to > 3× seawater. Even though the physiological effector mechanisms of osmoregulation that maintain plasma homeostasis in fresh water and seawater fish are well known, the corresponding molecular mechanisms that control switching between hyper- (fresh water) and hypo-osmoregulation (seawater) remain mostly elusive. In this study we show that hyperosmotic induction of glutamine synthetase represents a prominent part of this switch. Proteomics analysis of the O. mossambicus OmB cell line revealed that glutamine synthetase is transcriptionally regulated by hyperosmolality. Therefore, the 5' regulatory sequence of O. mossambicus glutamine synthetase was investigated. Using an enhancer trapping assay, we discovered a novel osmosensitive mechanism by which intron 1 positively mediates glutamine synthetase transcription. Intron 1 includes a single, functional copy of an osmoresponsive element, osmolality/salinity-responsive enhancer 1 (OSRE1). Unlike for conventional enhancers, the hyperosmotic induction of glutamine synthetase by intron 1 is position dependent. But irrespective of intron 1 position, OSRE1 deletion from intron 1 abolishes hyperosmotic enhancer activity. These findings indicate that proper intron 1 positioning and the presence of an OSRE1 in intron 1 are required for precise enhancement of hyperosmotic glutamine synthetase expression.

Transcriptomic analysis of dead end knockout testis reveals germ cell and gonadal somatic factors in Atlantic salmon

Background

Sustainability challenges are currently hampering an increase in salmon production. Using sterile salmon can solve problems with precocious puberty and genetic introgression from farmed escapees to wild populations. Recently sterile salmon was produced by knocking out the germ cell-specific dead end (dnd). Several approaches may be applied to inhibit Dnd function, including gene knockout, knockdown or immunization. Since it is challenging to develop a successful treatment against a gene product already existing in the body, alternative targets are being explored. Germ cells are surrounded by, and dependent on, gonadal somatic cells. Targeting genes essential for the survival of gonadal somatic cells may be good alternative targets for sterility treatments. Our aim was to identify and characterize novel germ cell and gonadal somatic factors in Atlantic salmon.

Results

We have for the first time analysed RNA-sequencing data from germ cell-free (GCF)/dnd knockout and wild type (WT) salmon testis and searched for genes preferentially expressed in either germ cells or gonadal somatic cells. To exclude genes with extra-gonadal expression, our dataset was merged with available multi-tissue transcriptome data. We identified 389 gonad specific genes, of which 194 were preferentially expressed within germ cells, and 11 were confined to gonadal somatic cells. Interestingly, 5 of the 11 gonadal somatic transcripts represented genes encoding secreted TGF-β factors; gsdf, inha, nodal and two bmp6-like genes, all representative vaccine targets. Of these, gsdf and inha had the highest transcript levels. Expression of gsdf and inha was further confirmed to be gonad specific, and their spatial expression was restricted to granulosa and Sertoli cells of the ovary and testis, respectively. Finally, we show that inha expression increases with puberty in both ovary and testis tissue, while gsdf expression does not change or decreases during puberty in ovary and testis tissue, respectively.

Conclusions

This study contributes with transcriptome data on salmon testis tissue with and without germ cells. We provide a list of novel and known germ cell- and gonad somatic specific transcripts, and show that the expression of two highly active gonadal somatic secreted TGF-β factors, gsdf and inha, are located within granulosa and Sertoli cells.

Dissecting Trait Variation across Species Barriers

Dissecting the basis of naturally occurring trait variation is one of the central goals of modern genetics. For eukaryotes, classic methods for this purpose rely on screens of recombinants from matings between distinct parents. These tools cannot be used in studies of species that cannot mate to form recombinant progeny in the first place. However, new approaches are coming online to shuffle the genomes of otherwise incompatible species. With them, geneticists can elucidate how evolution built a new trait, even if it happened millions of years ago, in a lineage that is now reproductively cutoff from its closest relatives.

Searching for Sympatric Speciation in the Genomic Era

Sympatric speciation illustrates how natural and sexual selection may create new species in isolation without geographic barriers. However, recent genomic reanalyses of classic examples of sympatric speciation reveal complex histories of secondary gene flow from outgroups into the radiation. In contrast, the rich theoretical literature on this process distinguishes among a diverse range of models based on simple genetic histories and different types of reproductive isolating barriers. Thus, there is a need to revisit how to connect theoretical models of sympatric speciation and their predictions to empirical case studies in the face of widespread gene flow. Here, theoretical differences among different types of sympatric speciation and speciation-with-gene-flow models are reviewed and summarized, and genomic analyses are proposed for distinguishing which models apply to case studies based on the timing and function of adaptive introgression. Investigating whether secondary gene flow contributed to reproductive isolation is necessary to test whether predictions of theory are ultimately borne out in nature.

Efficient CRISPR-Cas9 editing of major evolutionary loci in sticklebacks

BACKGROUND

Stickleback fish are widely used to study the genetic and ecological basis of phenotypic evolution. Although several major loci have now been identified that contribute to evolutionary differences between wild populations, further study of the phenotypes associated with particular genes and mutations has been limited by the difficulty of generating targeted mutations at precise locations in the stickleback genome.

APPROACH AND AIMS

We compared different methods of expressing single-guide RNAs (sgRNAs) and Cas9 activity in fertilized stickleback eggs. We used an easily scored pigmentation gene (SLC24A5) to screen for molecular lesions, phenotypic effects, and possible germline transmission of newly induced alleles. We then used the optimized CRISPR methods to target two major evolutionary loci in sticklebacks, KITLG and EDA. We hypothesized that coding region mutations in the KITLG gene would alter body pigmentation and possibly sex determination, and that mutations in the EDA gene would disrupt the formation of most armor plates, fin rays, spines, teeth, and gill rakers.

RESULTS

Targeted deletions were successfully induced at each target locus by co-injecting one-cell stage stickleback embryos with either Cas9 mRNA or Cas9 protein, together with sgRNAs designed to protein-coding exons. Founder animals were typically mosaic for multiple mutations, which they transmitted through the germline at overall rates of 21 to 100%. We found that the copy of KITLG on the X chromosome (KITLGX) has diverged from the KITLG on the Y chromosome (KITLGY). Predicted loss-of-function mutations in the KITLGX gene dramatically altered pigmentation in both external skin and internal organ, but the same was not true for KITLGY mutations. Predicted loss-of-function mutations in either the KITLGX or KITLGY genes did not lead to sex reversal or prevent fertility. Homozygous loss-of-function mutations in the EDA gene led to complete loss of armor plates, severe reduction or loss of most soft rays in the dorsal, anal, and caudal fins, and severe reductions in tooth and gill raker number. In contrast, long dorsal and pelvic spines remained intact in EDA mutant animals, suggesting that common co-segregation of plate loss and spine reduction in wild populations is unlikely to be due to pleiotropic effects of EDA mutations.

CONCLUSION

CRISPR-Cas9 approaches can be used to induce germline mutations in key evolutionary loci in sticklebacks. Targeted coding region mutations confirm an important role for KITLG and EDA in skin pigmentation and armor plate reduction, respectively. They also provide new information about the functions of these genes in other body structures.