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.

Genomic architecture of 5S rDNA cluster and its variations within and between species.

BACKGROUND: Ribosomal DNAs (rDNAs) are arranged in purely tandem repeats, preventing them from being reliably assembled onto chromosomes during generation of genome assembly. The uncertainty of rDNA genomic structure presents a significant barrier for studying their function and evolution. RESULTS: Here we generate ultra-long Oxford Nanopore Technologies (ONT) and short NGS reads to delineate the architecture and variation of the 5S rDNA cluster in the different strains of C. elegans and C. briggsae. We classify the individual rDNA's repeating units into 25 types based on the unique sequence variations in each unit of C. elegans (N2). We next perform assembly of the cluster by taking advantage of the long reads that carry these units, which led to an assembly of 5S rDNA cluster consisting of up to 167 consecutive 5S rDNA units in the N2 strain. The ordering and copy number of various rDNA units are consistent with the separation time between strains. Surprisingly, we observed a drastically reduced level of variation in the unit composition in the 5S rDNA cluster in the C. elegans CB4856 and C. briggsae AF16 strains than in the C. elegans N2 strain, suggesting that N2, a widely used reference strain, is likely to be defective in maintaining the 5S rDNA cluster stability compared with other wild isolates of C. elegans or C. briggsae. CONCLUSIONS: The results demonstrate that Nanopore DNA sequencing reads are capable of generating assembly of highly repetitive sequences, and rDNA units are highly dynamic both within and between population(s) of the same species in terms of sequence and copy number. The detailed structure and variation of the 5S rDNA units within the rDNA cluster pave the way for functional and evolutionary studies.

Technical challenges in defining RNA modifications.

It is well established that DNA base modifications play a key role in gene regulation during development and in response to environmental stress. This type of epigenetic control of development and environmental responses has been intensively studied over the past few decades. Similar to DNA, various RNA species also undergo modifications that play important roles in, for example, RNA splicing, protein translation, and the avoidance of immune surveillance by host. More than 160 different types of RNA modifications have been identified. In addition to base modifications, RNA modification also involves splicing of pre-mRNAs, leading to as many as tens of transcript isoforms from a single pre-RNA, especially in higher organisms. However, the function, prevalence and distribution of RNA modifications are poorly understood. The lack of a suitable method for the reliable identification of RNA modifications constitutes a significant challenge to studying their functions. This review focuses on the technologies that enable de novo identification of RNA base modifications and the alternatively spliced mRNA transcripts.

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.

Specific Interactions Between Autosome and X Chromosomes Cause Hybrid Male Sterility in Caenorhabditis Species.

Hybrid male progeny from interspecies crosses are more prone to sterility or inviability than hybrid female progeny, and the male sterility and inviability often demonstrate parent-of-origin asymmetry. However, the underlying genetic mechanism of asymmetric sterility or inviability remains elusive. We previously established a genome-wide hybrid incompatibility (HI) landscape between Caenorhabditis briggsae and C. nigoni by phenotyping a large collection of C. nigoni strains each carrying a C. briggsae introgression. In this study, we systematically dissect the genetic mechanism of asymmetric sterility and inviability in both hybrid male and female progeny between the two species. Specifically, we performed reciprocal crosses between C. briggsae and different C. nigoni strains that each carry a GFP-labeled C. briggsae genomic fragment referred to as introgression, and scored the HI phenotypes in the F1 progeny. The aggregated introgressions cover 94.6% of the C. briggsae genome, including 100% of the X chromosome. Surprisingly, we observed that two C. briggsaeX fragments that produce C. nigoni male sterility as an introgression rescued hybrid F1 sterility in males fathered by C. briggsae Subsequent backcrossing analyses indicated that a specific interaction between the X-linked interaction and one autosome introgression is required to rescue the hybrid male sterility. In addition, we identified another two C. briggsae genomic intervals on chromosomes II and IV that can rescue the inviability, but not the sterility, of hybrid F1 males fathered by C. nigoni, suggesting the involvement of differential epistatic interactions in the asymmetric hybrid male fertility and inviability. Importantly, backcrossing of the rescued sterile males with C. nigoni led to the isolation of a 1.1-Mb genomic interval that specifically interacts with an X-linked introgression, which is essential for hybrid male fertility. We further identified three C. briggsae genomic intervals on chromosome I, II, and III that produced inviability in all F1 progeny, dependent on or independent of the parent-of-origin. Taken together, we identified multiple independent interacting loci that are responsible for asymmetric hybrid male and female sterility, and inviability, which lays a foundation for their molecular characterization.

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.

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.