Specific down-regulation of spermatogenesis genes targeted by 22G RNAs in hybrid sterile males associated with an X-Chromosome introgression

Efficient targeted recombination with CRISPR/Cas9 in hybrids of Caenorhabditis nematodes with suppressed recombination

BACKGROUND

Homology-based recombination (HR) is the cornerstone of genetic mapping. However, a lack of sufficient sequence homology or the presence of a genomic rearrangement prevents HR through crossing, which inhibits genetic mapping in relevant genomic regions. This is particularly true in species hybrids whose genomic sequences are highly divergent along with various genome arrangements, making the mapping of genetic loci, such as hybrid incompatibility (HI) loci, through crossing impractical. We previously mapped tens of HI loci between two nematodes, Caenorhabditis briggsae and C. nigoni, through the repeated backcrossing of GFP-linked C. briggsae fragments into C. nigoni. However, the median introgression size was over 7 Mb, indicating apparent HR suppression and preventing the subsequent cloning of the causative gene underlying a given HI phenotype. Therefore, a robust method that permits recombination independent of sequence homology is desperately desired.

RESULTS

Here, we report a method of highly efficient targeted recombination (TR) induced by CRISPR/Cas9 with dual guide RNAs (gRNAs), which circumvents the HR suppression in hybrids between the two species. We demonstrated that a single gRNA was able to induce efficient TR between highly homologous sequences only in the F1 hybrids but not in the hybrids that carry a GFP-linked C. briggsae fragment in an otherwise C. nigoni background. We achieved highly efficient TR, regardless of sequence homology or genetic background, when dual gRNAs were used that each specifically targeted one parental chromosome. We further showed that dual gRNAs were able to induce efficient TR within genomic regions that had undergone inversion, in which HR-based recombination was expected to be suppressed, supporting the idea that dual-gRNA-induced TR can be achieved through nonhomology-based end joining between two parental chromosomes.

CONCLUSIONS

Recombination suppression can be circumvented through CRISPR/Cas9 with dual gRNAs, regardless of sequence homology or the genetic background of the species hybrid. This method is expected to be applicable to other situations in which recombination is suppressed in interspecies or intrapopulation hybrids.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Highly efficient transgenesis with miniMos in Caenorhabditis briggsae

Caenorhabditis briggsae as a companion species for Caenorhabditis elegans has played an increasingly important role in study of evolution of development and genome and gene regulation. Aided by the isolation of its sister spices, it has recently been established as a model for speciation study. To take full advantage of the species for comparative study, an effective transgenesis method especially those with single-copy insertion is important for functional comparison. Here, we improved a transposon-based transgenesis methodology that had been originally developed in C. elegans but worked marginally in C. briggsae. By incorporation of a heat shock step, the transgenesis efficiency in C. briggsae with a single-copy insertion is comparable to that in C. elegans. We used the method to generate 54 independent insertions mostly consisting of a mCherry tag over the C. briggsae genome. We demonstrated the use of the tags in identifying interacting loci responsible for hybrid male sterility between C. briggsae and Caenorhabditis nigoni when combined with the GFP tags we generated previously. Finally, we demonstrated that C. briggsae tolerates the C. elegans toxin, PEEL-1, but not SUP-35, making the latter a potential negative selection marker against extrachromosomal array.

Genetic exchange with an outcrossing sister species causes severe genome-wide dysregulation in a selfing Caenorhabditis nematode

Different modes of reproduction evolve rapidly, with important consequences for genome composition. Selfing species often occupy a similar niche as their outcrossing sister species with which they are able to mate and produce viable hybrid progeny, raising the question of how they maintain genomic identity. Here, we investigate this issue by using the nematode Caenorhabditis briggsae, which reproduces as a hermaphrodite, and its outcrossing sister species Caenorhabditis nigoni We hypothesize that selfing species might develop some barriers to prevent gene intrusions through gene regulation. We therefore examined gene regulation in the hybrid F2 embryos resulting from reciprocal backcrosses between F1 hybrid progeny and C. nigoni or C. briggsae F2 hybrid embryos with ∼75% of their genome derived from C. briggsae (termed as bB2) were inviable, whereas those with ∼75% of their genome derived from C. nigoni (termed as nB2) were viable. Misregulation of transposable elements, coding genes, and small regulatory RNAs was more widespread in the bB2 compared with the nB2 hybrids, which is a plausible explanation for the differential phenotypes between the two hybrids. Our results show that regulation of the C. briggsae genome is strongly affected by genetic exchanges with its outcrossing sister species, C. nigoni, whereas regulation of the C. nigoni genome is more robust on genetic exchange with C. briggsae The results provide new insights into how selfing species might maintain their identity despite genetic exchanges with closely related outcrossing species.

Comparative transcriptomics between Drosophila mojavensis and D. arizonae reveals transgressive gene expression and underexpression of spermatogenesis-related genes in hybrid testes

Interspecific hybridization is a stressful condition that can lead to sterility and/or inviability through improper gene regulation in Drosophila species with a high divergence time. However, the extent of these abnormalities in hybrids of recently diverging species is not well known. Some studies have shown that in Drosophila, the mechanisms of postzygotic isolation may evolve more rapidly in males than in females and that the degree of viability and sterility is associated with the genetic distance between species. Here, we used transcriptomic comparisons between two Drosophila mojavensis subspecies and D. arizonae (repleta group, Drosophila) and identified greater differential gene expression in testes than in ovaries. We tested the hypothesis that the severity of the interspecies hybrid phenotype is associated with the degree of gene misregulation. We showed limited gene misregulation in fertile females and an increase in the amount of misregulation in males with more severe sterile phenotypes (motile vs. amotile sperm). In addition, for these hybrids, we identified candidate genes that were mostly associated with spermatogenesis dysfunction.

Widespread misregulation of inter-species hybrid transcriptomes due to sex-specific and sex-chromosome regulatory evolution

When gene regulatory networks diverge between species, their dysfunctional expression in inter-species hybrid individuals can create genetic incompatibilities that generate the developmental defects responsible for intrinsic post-zygotic reproductive isolation. Both cis- and trans-acting regulatory divergence can be hastened by directional selection through adaptation, sexual selection, and inter-sexual conflict, in addition to cryptic evolution under stabilizing selection. Dysfunctional sex-biased gene expression, in particular, may provide an important source of sexually-dimorphic genetic incompatibilities. Here, we characterize and compare male and female/hermaphrodite transcriptome profiles for sibling nematode species Caenorhabditis briggsae and C. nigoni, along with allele-specific expression in their F1 hybrids, to deconvolve features of expression divergence and regulatory dysfunction. Despite evidence of widespread stabilizing selection on gene expression, misexpression of sex-biased genes pervades F1 hybrids of both sexes. This finding implicates greater fragility of male genetic networks to produce dysfunctional organismal phenotypes. Spermatogenesis genes are especially prone to high divergence in both expression and coding sequences, consistent with a "faster male" model for Haldane's rule and elevated sterility of hybrid males. Moreover, underdominant expression pervades male-biased genes compared to female-biased and sex-neutral genes and an excess of cis-trans compensatory regulatory divergence for X-linked genes underscores a "large-X effect" for hybrid male expression dysfunction. Extensive regulatory divergence in sex determination pathway genes likely contributes to demasculinization of XX hybrids. The evolution of genetic incompatibilities due to regulatory versus coding sequence divergence, however, are expected to arise in an uncorrelated fashion. This study identifies important differences between the sexes in how regulatory networks diverge to contribute to sex-biases in how genetic incompatibilities manifest during the speciation process.

Speciation and the developmental alarm clock

New species arise as the genomes of populations diverge. The developmental 'alarm clock' of speciation sounds off when sufficient divergence in genetic control of development leads hybrid individuals to infertility or inviability, the world awoken to the dawn of new species with intrinsic post-zygotic reproductive isolation. Some developmental stages will be more prone to hybrid dysfunction due to how molecular evolution interacts with the ontogenetic timing of gene expression. Considering the ontogeny of hybrid incompatibilities provides a profitable connection between 'evo-devo' and speciation genetics to better link macroevolutionary pattern, microevolutionary process, and molecular mechanisms. Here, we explore speciation alongside development, emphasizing their mutual dependence on genetic network features, fitness landscapes, and developmental system drift. We assess models for how ontogenetic timing of reproductive isolation can be predictable. Experiments and theory within this synthetic perspective can help identify new rules of speciation as well as rules in the molecular evolution of development.

Direct full-length RNA sequencing reveals unexpected transcriptome complexity during Caenorhabditis elegans development

Massively parallel sequencing of the polyadenylated RNAs has played a key role in delineating transcriptome complexity, including alternative use of an exon, promoter, 5′ or 3′ splice site or polyadenylation site, and RNA modification. However, reads derived from the current RNA-seq technologies are usually short and deprived of information on modification, compromising their potential in defining transcriptome complexity. Here, we applied a direct RNA sequencing method with ultralong reads using Oxford Nanopore Technologies to study the transcriptome complexity in Caenorhabditis elegans. We generated approximately six million reads using native poly(A)-tailed mRNAs from three developmental stages, with average read lengths ranging from 900 to 1100 nt. Around half of the reads represent full-length transcripts. To utilize the full-length transcripts in defining transcriptome complexity, we devised a method to classify the long reads as the same as existing transcripts or as a novel transcript using sequence mapping tracks rather than existing intron/exon structures, which allowed us to identify roughly 57,000 novel isoforms and recover at least 26,000 out of the 33,500 existing isoforms. The sets of genes with differential expression versus differential isoform usage over development are largely different, implying a fine-tuned regulation at isoform level. We also observed an unexpected increase in putative RNA modification in all bases in the coding region relative to the UTR, suggesting their possible roles in translation. The RNA reads and the method for read classification are expected to deliver new insights into RNA processing and modification and their underlying biology in the future.

Males, Outcrossing, and Sexual Selection in Caenorhabditis Nematodes

Males of Caenorhabditis elegans provide a crucial practical tool in the laboratory, but, as the rarer and more finicky sex, have not enjoyed the same depth of research attention as hermaphrodites. Males, however, have attracted the attention of evolutionary biologists who are exploiting the C. elegans system to test longstanding hypotheses about sexual selection, sexual conflict, transitions in reproductive mode, and genome evolution, as well as to make new discoveries about Caenorhabditis organismal biology. Here, we review the evolutionary concepts and data informed by study of males of C. elegans and other Caenorhabditis We give special attention to the important role of sperm cells as a mediator of inter-male competition and male-female conflict that has led to drastic trait divergence across species, despite exceptional phenotypic conservation in many other morphological features. We discuss the evolutionary forces important in the origins of reproductive mode transitions from males being common (gonochorism: females and males) to rare (androdioecy: hermaphrodites and males) and the factors that modulate male frequency in extant androdioecious populations, including the potential influence of selective interference, host-pathogen coevolution, and mutation accumulation. Further, we summarize the consequences of males being common vs rare for adaptation and for trait divergence, trait degradation, and trait dimorphism between the sexes, as well as for molecular evolution of the genome, at both micro-evolutionary and macro-evolutionary timescales. We conclude that C. elegans male biology remains underexploited and that future studies leveraging its extensive experimental resources are poised to discover novel biology and to inform profound questions about animal function and evolution.

Comparative mitochondrial genomics reveals a possible role of a recent duplication of NADH dehydrogenase subunit 5 in gene regulation

Mitochondrial genome (mtDNA) carries not only well-conserved protein coding, tRNA and rRNA genes, but also highly variable non-coding regions (NCRs). However, the NCRs show poor conservation across species, making their function and evolution elusive. Identification and functional characterization of NCRs across species would be critical for addressing these questions. To this end, we devised a computational pipeline and performed de novo assembly and annotation of mtDNA from 19 Caenorhabditis species using next-generation sequencing (NGS) data. The mtDNAs for 14 out of the 19 species are reported for the first time. Comparison of the 19 genomes reveals species-specific sampling of partial displacement-loop (D-loop) sequence as a novel NCR inserted into a unique tRNA cluster, suggesting an important role of the D-loop and the tRNA cluster in shaping NCR evolution. Intriguingly, RNA-Seq analysis suggests that a novel NCR resulting from a recent duplication of NADH dehydrogenase subunit 5 (ND5) could be utilized as a 3' UTR for up-regulation of its upstream gene. The expression analysis shows a species- and sex-specific expression of mitochondrial genes encoded by mtDNA and nucleus, respectively. Our analyses provide important insights into the function and evolution of mitochondrial NCRs and pave the way for further studying the function and evolution of mitochondrial genome.

Genetic basis to hybrid inviability is more complex than hybrid male sterility in Caenorhabditis nematodes

Hybrid male sterility often evolves before female sterility or inviability of hybrids, implying that the accumulation of divergence between separated lineages should lead hybrid male sterility to have a more polygenic basis. However, experimental evidence is mixed. Here, we use the nematodes Caenorhabditis remanei and C. latens to characterize the underlying genetic basis of asymmetric hybrid male sterility and hybrid inviability. We demonstrate that hybrid male sterility is consistent with a simple genetic basis, involving a single X-autosome incompatibility. We also show that hybrid inviability involves more genomic compartments, involving diverse nuclear-nuclear incompatibilities, a mito-nuclear incompatibility, and maternal effects. These findings demonstrate that male sensitivity to genetic perturbation may be genetically simple compared to hybrid inviability in Caenorhabditis and motivates tests of generality for the genetic architecture of hybrid incompatibility across the breadth of phylogeny.

Genomic basis of recombination suppression in the hybrid between Caenorhabditis briggsae and C. nigoni

DNA recombination is required for effective segregation and diversification of genomes and for the successful completion of meiosis. Recent studies in various species hybrids have demonstrated a genetic link between DNA recombination and speciation. Consistent with this, we observed a striking suppression of recombination in the hybrids between two nematodes, the hermaphroditic Caenorhabditis briggsae and the gonochoristic C. nigoni. To unravel the molecular basis underlying the recombination suppression in their hybrids, we generated a C. nigoni genome with chromosome-level contiguity and produced an improved C. briggsae genome with resolved gaps up to 2.8 Mb. The genome alignment reveals not only high sequence divergences but also pervasive intra- and inter-chromosomal sequence re-arrangements between the two species, which are plausible culprits for the observed suppression. Comparison of recombination boundary sequences suggests that recombination in the hybrid requires extensive sequence homology, which is rarely seen between the two genomes. The new genomes and genomic libraries form invaluable resources for studying genome evolution, hybrid incompatibilities and sex evolution for this pair of model species.

Rapid genome shrinkage in a self-fertile nematode reveals sperm competition proteins

To reveal impacts of sexual mode on genome content, we compared chromosome-scale assemblies of the outcrossing nematode Caenorhabditis nigoni to its self-fertile sibling species, C. briggsaeC. nigoni's genome resembles that of outcrossing relatives but encodes 31% more protein-coding genes than C. briggsaeC. nigoni genes lacking C. briggsae orthologs were disproportionately small and male-biased in expression. These include the male secreted short (mss) gene family, which encodes sperm surface glycoproteins conserved only in outcrossing species. Sperm from mss-null males of outcrossing C. remanei failed to compete with wild-type sperm, despite normal fertility in noncompetitive mating. Restoring mss to C. briggsae males was sufficient to enhance sperm competitiveness. Thus, sex has a pervasive influence on genome content that can be used to identify sperm competition factors.

X exceptionalism in Caenorhabditis speciation

Speciation genetics research in diverse organisms shows the X-chromosome to be exceptional in how it contributes to "rules" of speciation. Until recently, however, the nematode phylum has been nearly silent on this issue, despite the model organism Caenorhabditis elegans having touched most other topics in biology. Studies of speciation with Caenorhabditis accelerated with the recent discovery of species pairs showing partial interfertility. The resulting genetic analyses of reproductive isolation in nematodes demonstrate key roles for the X-chromosome in hybrid male sterility and inviability, opening up new understanding of the genetic causes of Haldane's rule, Darwin's corollary to Haldane's rule, and enabling tests of the large-X effect hypothesis. Studies to date implicate improper chromatin regulation of the X-chromosome by small RNA pathways as integral to hybrid male dysfunction. Sexual transitions in reproductive mode to self-fertilizing hermaphroditism inject distinctive molecular evolutionary features into the speciation process for some species. Caenorhabditis also provides unique opportunities for analysis in a system with XO sex determination that lacks a Y-chromosome, sex chromosome-dependent sperm competition differences and mechanisms of gametic isolation, exceptional accessibility to the development process and rapid experimental evolution. As genetic analysis of reproductive isolation matures with investigation of multiple pairs of Caenorhabditis species and new species discovery, nematodes will provide a powerful complement to more established study organisms for deciphering the genetic basis of and rules to speciation.

Non-Mendelian assortment of homologous autosomes of different sizes in males is the ancestral state in the Caenorhabditis lineage

Organismal genome sizes vary by six orders of magnitude and appear positively correlated with organismal size and complexity. Neutral models have been proposed to explain the broad patterns of genome size variation based on organism population sizes. In the Caenorhabditis genus, hermaphrodite genomes are smaller than those of gonochoristic species. One possible driving force for this genome size difference could be non-random chromosome segregation. In Caenorhabditis elegans, chromosome assortment is non-independent and violates Mendel's second law. In males, the shorter homologue of a heterozygous autosome pair preferentially co-segregates with the X chromosome while the longer one preferentially co-segregates with the nullo-X (O) chromosome in a process we call "skew". Since hermaphrodites preferentially receive the shorter chromosomes and can start populations independently, their genome size would be predicted to decrease over evolutionary time. If skew is an important driver for genome size reduction in hermaphroditic Caenorhabditis species, then it should be present in all congeneric species. In this study, we tested this hypothesis and found that skew is present in all eight examined species. Our results suggest that skew is likely the ancestral state in this genus. More speculatively, skew may drive genome size patterns in hermaphroditic species in other nematodes.

Beyond speciation genes: an overview of genome stability in evolution and speciation

Genome stability ensures individual fitness and reliable transmission of genetic information. Hybridization between diverging lineages can trigger genome instability, highlighting its potential role in post-zygotic reproductive isolation. We argue that genome instability is not merely one of several types of hybrid incompatibility, but rather that genome stability is one of the very first and most fundamental traits that can break down when two diverged genomes are combined. Future work will reveal how frequent and predictable genome instability is in hybrids, how it affects hybrid fitness, and whether it is a direct cause or consequence of speciation.

Comparative proteome analysis between C . briggsae embryos and larvae reveals a role of chromatin modification proteins in embryonic cell division

Caenorhabditis briggsae has emerged as a model for comparative biology against model organism C. elegans. Most of its cell fate specifications are completed during embryogenesis whereas its cell growth is achieved mainly in larval stages. The molecular mechanism underlying the drastic developmental changes is poorly understood. To gain insights into the molecular changes between the two stages, we compared the proteomes between the two stages using iTRAQ. We identified a total of 2,791 proteins in the C. briggsae embryos and larvae, 247 of which undergo up- or down-regulation between the two stages. The proteins that are upregulated in the larval stages are enriched in the Gene Ontology categories of energy production, protein translation, and cytoskeleton; whereas those upregulated in the embryonic stage are enriched in the categories of chromatin dynamics and posttranslational modification, suggesting a more active chromatin modification in the embryos than in the larva. Perturbation of a subset of chromatin modifiers followed by cell lineage analysis suggests their roles in controlling cell division pace. Taken together, we demonstrate a general molecular switch from chromatin modification to metabolism during the transition from C. briggsae embryonic to its larval stages using iTRAQ approach. The switch might be conserved across metazoans.

Revisiting Suppression of Interspecies Hybrid Male Lethality in Caenorhabditis Nematodes

Within the nematode genus Caenorhabditis, Caenorhabditis briggsae and C. nigoni are among the most closely related species known. They differ in sexual mode, with C. nigoni retaining the ancestral XO male-XX female outcrossing system, while C. briggsae recently evolved self-fertility and an XX-biased sex ratio. Wild-type C. briggsae and C. nigoni can produce fertile hybrid XX female progeny, but XO progeny are either 100% inviable (when C. briggsae is the mother) or viable but sterile (when C. nigoni is the mother). A recent study provided evidence suggesting that loss of the Cbr-him-8 meiotic regulator in C. briggsae hermaphrodites allowed them to produce viable and fertile hybrid XO male progeny when mated to C. nigoni Because such males would be useful for a variety of genetic experiments, we sought to verify this result. Preliminary crosses with wild-type C. briggsae hermaphrodites occasionally produced fertile males, but they could not be confirmed to be interspecies hybrids. Using an RNA interference (RNAi) protocol that eliminates any possibility of self-progeny in Cbr-him-8 hermaphrodites, we found sterile males bearing the C. nigoni X chromosome, but no fertile males bearing the C. briggsae X, as in wild-type crosses. Our results suggest that the apparent rescue of XO hybrid viability and fertility is due to incomplete purging of self-sperm prior to mating.