A phylogeny of caenorhabditis reveals frequent loss of introns during nematode evolution

Recurrent but Short-Lived Duplications of Centromeric Proteins in Holocentric Caenorhabditis Species

Centromeric histones (CenH3s) are essential for chromosome inheritance during cell division in most eukaryotes. CenH3 genes have rapidly evolved and undergone repeated gene duplications and diversification in many plant and animal species. In Caenorhabditis species, two independent duplications of CenH3 (named hcp-3 for HoloCentric chromosome-binding Protein 3) were previously identified in C. elegans and C. remanei. Using phylogenomic analyses in 32 Caenorhabditis species, we find strict retention of the ancestral hcp-3 gene and 10 independent duplications. Most hcp-3L (hcp-3-like) paralogs are only found in 1-2 species, are expressed in both males and females/hermaphrodites, and encode histone fold domains with 69-100% identity to ancestral hcp-3. We identified novel N-terminal protein motifs, including putative kinetochore protein-interacting motifs and a potential separase cleavage site, which are well conserved across Caenorhabditis HCP-3 proteins. Other N-terminal motifs vary in their retention across paralogs or species, revealing potential subfunctionalization or functional loss following duplication. An N-terminal extension in the hcp-3L gene of C. afra revealed an unprecedented protein fusion, where hcp-3L fused to duplicated segments from hcp-4 (nematode CENP-C). By extending our analyses beyond CenH3, we found gene duplications of six inner and outer kinetochore genes in Caenorhabditis, which appear to have been retained independent of hcp-3 duplications. Our findings suggest that centromeric protein duplications occur frequently in Caenorhabditis nematodes, are selectively retained for short evolutionary periods, then degenerate or are lost entirely. We hypothesize that unique challenges associated with holocentricity in Caenorhabditis may lead to this rapid "revolving door" of kinetochore protein paralogs.

Intron losses and gains in the nematodes

Background

The evolution of spliceosomal introns has been widely studied among various eukaryotic groups. Researchers nearly reached the consensuses on the pattern and the mechanisms of intron losses and gains across eukaryotes. However, according to previous studies that analyzed a few genes or genomes, Nematoda seems to be an eccentric group.

Results

Taking advantage of the recent accumulation of sequenced genomes, we extensively analyzed the intron losses and gains using 104 nematode genomes across all the five Clades of the phylum. Nematodes have a wide range of intron density, from less than one to more than nine per kbp coding sequence. The rates of intron losses and gains exhibit significant heterogeneity both across different nematode lineages and across different evolutionary stages of the same lineage. The frequency of intron losses far exceeds that of intron gains. Five pieces of evidence supporting the model of cDNA-mediated intron loss have been observed in ten Caenorhabditis species, the dominance of the precise intron losses, frequent loss of adjacent introns, high-level expression of the intron-lost genes, preferential losses of short introns, and the preferential losses of introns close to 3′-ends of genes. Like studies in most eukaryotic groups, we cannot find the source sequences for the limited number of intron gains detected in the Caenorhabditis genomes.

Conclusions

These results indicate that nematodes are a typical eukaryotic group rather than an outlier in intron evolution.

Chromosome-level reference genomes for two strains of Caenorhabditis briggsae: an improved platform for comparative genomics

The publication of the Caenorhabditis briggsae reference genome in 2003 enabled the first comparative genomics studies between C. elegans and C. briggsae, shedding light on the evolution of genome content and structure in the Caenorhabditis genus. However, despite being widely used, the currently available C. briggsae reference genome is substantially less complete and structurally accurate than the C. elegans reference genome. Here, we used high-coverage Oxford Nanopore long-read and chromosome conformation capture data to generate chromosome-level reference genomes for two C. briggsae strains: QX1410, a new reference strain closely related to the laboratory AF16 strain, and VX34, a highly divergent strain isolated in China. We also sequenced 99 recombinant inbred lines (RILs) generated from reciprocal crosses between QX1410 and VX34 to create a recombination map and identify chromosomal domains. Additionally, we used both short- and long-read RNA sequencing (RNA-seq) data to generate high-quality gene annotations. By comparing these new reference genomes to the current reference, we reveal that hyper-divergent haplotypes cover large portions of the C. briggsae genome, similar to recent reports in C. elegans and C. tropicalis. We also show that the genomes of selfing Caenorhabditis species have undergone more rearrangement than their outcrossing relatives, which has biased previous estimates of rearrangement rate in Caenorhabditis. These new genomes provide a substantially improved platform for comparative genomics in Caenorhabditis and narrow the gap between the quality of genomic resources available for C. elegans and C. briggsae.

Sex Determination in Nematode Germ Cells

BACKGROUND

Animal germ cells differentiate as sperm or as oocytes. These sexual fates are controlled by complex regulatory pathways to ensure that the proper gametes are made at the appropriate times.

SUMMARY

Nematodes like Caenorhabditis elegans and its close relatives are ideal models for studying how this regulation works, because the XX animals are self-fertile hermaphrodites that produce both sperm and oocytes. In these worms, germ cells use the same signal transduction pathway that functions in somatic cells. This pathway determines the activity of the transcription factor TRA-1, a Gli protein that can repress male genes. However, the pathway is extensively modified in germ cells, largely by the action of translational regulators like the PUF proteins. Many of these modifications play critical roles in allowing the XX hermaphrodites to make sperm in an otherwise female body. Finally, TRA-1 cooperates with chromatin regulators in the germ line to control the activity of fog-1 and fog-3, which are essential for spermatogenesis. FOG-1 and FOG-3 work together to determine germ cell fates by blocking the translation of oogenic transcripts. Key Messages: Although there is great diversity in how germ cell fates are controlled in other animals, many of the key nematode genes are conserved, and the critical role of translational regulators may be universal.

Chromosome-Wide Evolution and Sex Determination in the Three-Sexed Nematode Auanema rhodensis

Trioecy, a mating system in which males, females and hermaphrodites co-exist, is a useful system to investigate the origin and maintenance of alternative mating strategies. In the trioecious nematode Auanema rhodensis, males have one X chromosome (XO), whereas females and hermaphrodites have two (XX). The female vs. hermaphrodite sex determination mechanisms have remained elusive. In this study, RNA-seq analyses show a 20% difference between the L2 hermaphrodite and female gene expression profiles. RNAi experiments targeting the DM (doublesex/mab-3) domain transcription factor dmd-10/11 suggest that the hermaphrodite sexual fate requires the upregulation of this gene. The genetic linkage map (GLM) shows that there is chromosome-wide heterozygosity for the X chromosome in F2 hermaphrodite-derived lines originated from crosses between two parental inbred strains. These results confirm the lack of recombination of the X chromosome in hermaphrodites, as previously reported. We also describe conserved chromosome elements (Nigon elements), which have been mostly maintained throughout the evolution of Rhabditina nematodes. The seven-chromosome karyotype of A. rhodensis, instead of the typical six found in other rhabditine species, derives from fusion/rearrangements events involving three Nigon elements. The A. rhodensis X chromosome is the smallest and most polymorphic with the least proportion of conserved genes. This may reflect its atypical mode of father-to-son transmission and its lack of recombination in hermaphrodites and males. In conclusion, this study provides a framework for studying the evolution of chromosomes in rhabditine nematodes, as well as possible mechanisms for the sex determination in a three-sexed species.

Comparative genomics of 10 new Caenorhabditis species

The nematode Caenorhabditis elegans has been central to the understanding of metazoan biology. However, C. elegans is but one species among millions and the significance of this important model organism will only be fully revealed if it is placed in a rich evolutionary context. Global sampling efforts have led to the discovery of over 50 putative species from the genus Caenorhabditis, many of which await formal species description. Here, we present species descriptions for 10 new Caenorhabditis species. We also present draft genome sequences for nine of these new species, along with a transcriptome assembly for one. We exploit these whole-genome data to reconstruct the Caenorhabditis phylogeny and use this phylogenetic tree to dissect the evolution of morphology in the genus. We reveal extensive variation in genome size and investigate the molecular processes that underlie this variation. We show unexpected complexity in the evolutionary history of key developmental pathway genes. These new species and the associated genomic resources will be essential in our attempts to understand the evolutionary origins of the C. elegans model.

From "the Worm" to "the Worms" and Back Again: The Evolutionary Developmental Biology of Nematodes

Since the earliest days of research on nematodes, scientists have noted the developmental and morphological variation that exists within and between species. As various cellular and developmental processes were revealed through intense focus on Caenorhabditis elegans, these comparative studies have expanded. Within the genus Caenorhabditis, they include characterization of intraspecific polymorphisms and comparisons of distinct species, all generally amenable to the same laboratory culture methods and supported by robust genomic and experimental tools. The C. elegans paradigm has also motivated studies with more distantly related nematodes and animals. Combined with improved phylogenies, this work has led to important insights about the evolution of nematode development. First, while many aspects of C. elegans development are representative of Caenorhabditis, and of terrestrial nematodes more generally, others vary in ways both obvious and cryptic. Second, the system has revealed several clear examples of developmental flexibility in achieving a particular trait. This includes developmental system drift, in which the developmental control of homologous traits has diverged in different lineages, and cases of convergent evolution. Overall, the wealth of information and experimental techniques developed in C. elegans is being leveraged to make nematodes a powerful system for evolutionary cellular and developmental biology.

Heterogeneity of intron presence/absence in Olifantiella sp. (Bacillariophyta) contributes to the understanding of intron loss

Although hypotheses have been proposed and developed to interpret the origins and functions of introns, substantial controversies remain about the mechanism of intron evolution. The availability of introns in the intermediate state is quite helpful for resolving this debate. In this study, a new strain of diatom (denominated as DB21-1) was isolated and identified as Olifantiella sp., which possesses multiple types of 18S rDNAs (obtained from genomic DNA; lengths ranged from 2,056 bp to 2,988 bp). Based on alignments between 18S rDNAs and 18S rRNA (obtained from cDNA; 1,783 bp), seven intron insertion sites (IISs) located in the 18S rDNA were identified, each of which displayed the polymorphism of intron presence/absence. Specific primers around each IIS were designed to amplify the introns and the results indicated that introns in the same IIS varied in lengths, while terminal sequences were conserved. Our study showed that the process of intron loss happens via a series of successive steps, and each step could derive corresponding introns under intermediate states. Moreover, these results indicate that the mechanism of genomic deletion that occurs at DNA level can also lead to exact intron loss.

Parallel loss of introns in the ABCB1 gene in angiosperms

BACKGROUND

The presence of non-coding introns is a characteristic feature of most eukaryotic genes. While the size of the introns, number of introns per gene and the number of intron-containing genes can vary greatly between sequenced eukaryotic genomes, the structure of a gene with reference to intron presence and positions is typically conserved in closely related species. Unexpectedly, the ABCB1 (ATP-Binding Cassette Subfamily B Member 1) gene which encodes a P-glycoprotein and underlies dwarfing traits in maize (br2), sorghum (dw3) and pearl millet (d2) displayed considerable variation in intron composition.

RESULTS

An analysis of the ABCB1 gene structure in 80 angiosperms revealed that the number of introns ranged from one to nine. All introns in ABCB1 underwent either a one-time loss (single loss in one lineage/species) or multiple independent losses (parallel loss in two or more lineages/species) with the majority of losses occurring within the grass family. In contrast, the structure of the closest homolog to ABCB1, ABCB19, remained constant in the majority of angiosperms analyzed. Using known phylogenetic relationships within the grasses, we determined the ancestral branch-points where the losses occurred. Intron 7, the longest intron, was lost in only a single species, Mimulus guttatus, following duplication of ABCB1. Semiquantitative PCR showed that the M. guttatus ABCB1 gene copy without intron 7 had significantly lower transcript levels than the gene copy with intron 7. We further demonstrated that intron 7 carried two motifs that were highly conserved across the monocot-dicot divide.

CONCLUSIONS

The ABCB1 gene structure is highly dynamic, while the structure of ABCB19 remained largely conserved through evolution. Precise removal of introns, preferential removal of smaller introns and presence of at least 2 bp of microhomology flanking most introns indicated that intron loss may have predominantly occurred through non-homologous end-joining (NHEJ) repair of double strand breaks. Lack of microhomology in the exon upstream of lost phase I introns was likely due to release of the selective constraint on the penultimate base (3rd base in codon) of the terminal codon by the splicing machinery. In addition to size, the presence of regulatory motifs will make introns recalcitrant to loss.

Trans-splicing enhances translational efficiency in C. elegans

Translational efficiency is subject to extensive regulation. However, the factors influencing such regulation are poorly understood. In Caenorhabditis elegans, 62% of genes are trans-spliced to a specific spliced leader (SL1), which replaces part of the native 5' untranslated region (5' UTR). Given the pivotal role the 5' UTR plays in the regulation of translational efficiency, we hypothesized that SL1 trans-splicing functions to regulate translational efficiency. With genome-wide analysis on Ribo-seq data, polysome profiling experiments, and CRISPR-Cas9-based genetic manipulation of trans-splicing sites, we found four lines of evidence in support of this hypothesis. First, SL1 trans-spliced genes have higher translational efficiencies than non-trans-spliced genes. Second, SL1 trans-spliced genes have higher translational efficiencies than non-trans-spliced orthologous genes in other nematode species. Third, an SL1 trans-spliced isoform has higher translational efficiency than the non-trans-spliced isoform of the same gene. Fourth, deletion of trans-splicing sites of endogenous genes leads to reduced translational efficiency. Importantly, we demonstrated that SL1 trans-splicing plays a key role in enhancing translational efficiencies of essential genes. We further discovered that SL1 trans-splicing likely enhances translational efficiency by shortening the native 5' UTRs, hence reducing the presence of upstream start codons (uAUG) and weakening mRNA secondary structures. Taken together, our study elucidates the global function of trans-splicing in enhancing translational efficiency in nematodes, paving the way for further understanding the genomic mechanisms of translational regulation.

Differential regulation of germ line apoptosis and germ cell differentiation by CPEB family members in C. elegans

Cytoplasmic polyadenylation element binding (CPEB) proteins are evolutionary conserved RNA-binding proteins that control mRNA polyadenylation and translation. Orthologs in humans and other vertebrates are mainly involved in oogenesis. This is also the case for the C. elegans CPEB family member CPB-3, whereas two further CPEB proteins (CPB-1 and FOG-1) are involved in spermatogenesis. Here we describe the characterisation of a new missense allele of cpb-3 and show that loss of cpb-3 function leads to an increase in physiological germ cell death. To better understand the interaction and effect of C. elegans CPEB proteins on processes such as physiological apoptosis, germ cell differentiation, and regulation of gene expression, we characterised changes in the transcriptome and proteome of C. elegans CPEB mutants. Our results show that, despite their sequence similarities CPEB family members tend to have distinct overall effects on gene expression (both at the transcript and protein levels). This observation is consistent with the distinct phenotypes observed in the various CPEB family mutants.

From intronization to intron loss: How the interplay between mRNA-associated processes can shape the architecture and the expression of eukaryotic genes

Transcription-coupled processes such as capping, splicing, and cleavage/polyadenylation participate in the journey from genes to proteins. Although they are traditionally thought to serve only as steps in the generation of mature mRNAs, a synthesis of available data indicates that these processes could also act as a driving force for the evolution of eukaryotic genes. A theoretical framework for how mRNA-associated processes may shape gene structure and expression has recently been proposed. Factors that promote splicing and cleavage/polyadenylation in this framework compete for access to overlapping or neighboring signals throughout the transcription cycle. These antagonistic interactions allow mechanisms for intron gain and splice site recognition as well as common trends in eukaryotic gene structure and expression to be coherently integrated. Here, I extend this framework further. Observations that largely (but not exclusively) revolve around the formation of DNA-RNA hybrid structures, called R loops, and promoter directionality are integrated. Additionally, the interplay between splicing factors and cleavage/polyadenylation factors is theorized to also affect the formation of intragenic DNA double-stranded breaks thereby contributing to intron loss. The most notable prediction in this proposition is that RNA molecules can mediate intron loss by serving as a template to repair DNA double-stranded breaks. The framework presented here leverages a vast body of empirical observations, logically extending previous suggestions, and generating verifiable predictions to further substantiate the view that the intracellular environment plays an active role in shaping the structure and the expression of eukaryotic genes.

Mating and male pheromone kill Caenorhabditis males through distinct mechanisms

Differences in longevity between sexes is a mysterious yet general phenomenon across great evolutionary distances. To test the roles of responses to environmental cues and sexual behaviors in longevity regulation, we examined Caenorhabditis male lifespan under solitary, grouped, and mated conditions. We find that neurons and the germline are required for male pheromone-dependent male death. Hermaphrodites with a masculinized nervous system secrete male pheromone and are susceptible to male pheromone killing. Male pheromone-mediated killing is unique to androdioecious Caenorhabditis, and may reduce the number of males in hermaphroditic populations; neither males nor females of gonochoristic species are susceptible to male pheromone killing. By contrast, mating-induced death, which is characterized by germline-dependent shrinking, glycogen loss, and ectopic vitellogenin expression, utilizes distinct molecular pathways and is shared between the sexes and across species. The study of sex- and species-specific regulation of aging reveals deeply conserved mechanisms of longevity and population structure regulation.

DOI:http://dx.doi.org/10.7554/eLife.23493.001

In many animals, different sexes have different life expectancies. This holds true for a roundworm species called Caenorhabditis elegans that has commonly been used to study aging and lifespan. Unlike some related Caenorhabditis roundworm species (which consist of male and female worms), C. elegans worms are predominantly hermaphrodites and reproduce by self-fertilization. C. elegans males are normally rare. However, under stressful conditions the number of males increases to reduce inbreeding and so help the worm population to adapt to the environment.

Investigations into the factors that affect the lifespan of C. elegans have mostly studied hermaphrodites. For example, one recent study showed that mating shortens the lifespan of hermaphrodites. Another study showed that pheromones – hormones that change the behavior of other worms – also shorten hermaphrodite lifespan. The male pheromone is produced by males and sensed by both males and hermaphrodites. But does mating and male pheromone affect the lifespan of male roundworms?

Shi et al. have now studied Caenorhabditis worms of different species and sexes to investigate how sexual behaviors and male pheromone regulate the lifespan of male roundworms. The results of the experiments revealed two distinct mechanisms of male death. Firstly, mating caused the males of many different Caenorhabditis species to shrink and die, and also killed females and hermaphrodites. Secondly, the males of hermaphroditic species – and only these males – could also be killed by male pheromone.

The results suggest that death from mating may be an unavoidable cost of reproducing that is seen across all sexes and species of roundworm. In contrast, death by male pheromone may be a way of culling the male population in hermaphroditic species, for example, after stressful conditions have caused a sudden increase in the number of male worms.

Further work is now needed to investigate the finer details of the mechanisms by which mating and male pheromone cause death. Ultimately, this work in Caenorhabditis could be extended to help us to understand how other animals regulate their lifespan and maintain an optimum ratio of the sexes.

DOI:http://dx.doi.org/10.7554/eLife.23493.002

Evolutionary dynamics of triosephosphate isomerase gene intron location pattern in Metazoa: A new perspective on intron evolution in animals

Intron evolution, including its dynamics in the evolutionary transitions and diversification of eukaryotes, remains elusive. Inadequate taxon sampling due to data shortage, unclear phylogenetic framework, and inappropriate outgroup application might be among the causes. Besides, the integrity of all the introns within a gene was often neglected previously. Taking advantage of the ancient conserved triosephosphate isomerase gene (tim), the relatively robust phylogeny of Metazoa, and choanoflagellates as outgroup, the evolutionary dynamics of tim intron location pattern (ILP) in Metazoa was investigated. From 133 representative species of ten phyla, 30 types of ILPs were identified. A most common one, which harbors the maximum six intron positions, is deduced to be the common ancestral tim ILP of Metazoa, which almost had formed in their protozoan ancestor and was surprisingly retained and passed down till to each ancestors of metazoan phyla. In the subsequent animal diversification, it underwent different evolutionary trajectories: within Deuterostomia, it was almost completely retained only with changes in a few species with relatively recently fast-evolving histories, while within the rapidly radiating Protostomia, besides few but remarkable retention, it usually displayed extensive intron losses and a few gains. Therefore, a common ancestral exon-intron arrangement pattern of an animal gene is definitely discovered; besides the 'intron-rich view' of early animal genes being confirmed, the novel insight that high exon-intron re-arrangements of genes seem to be associated with the relatively recently rapid evolution of lineages/species/genomes but have no correlation with the ancient major evolutionary transitions in animal evolution, is revealed.

"The persistence of memory"-Hermaphroditism in nematodes

Self-fertility has evolved many times in nematodes. This transition often produces an androdioecious species, with XX hermaphrodites and XO males. Although these hermaphrodites resemble females in most respects, early germ cells differentiate as sperm, and late ones as oocytes. The sperm then receive an activation signal, populate the spermathecae, and are stored for later use in self-fertilization. These traits are controlled by complex modifications to the sex-determination and sperm activation pathways, which have arisen independently during the evolution of each hermaphroditic species. This transformation in reproductive strategy then promotes other major changes in the development, evolution, and population structure of these animals. Mol. Reprod. Dev. 84: 144-157, 2017. © 2016 Wiley Periodicals, Inc.

Co-option of alternate sperm activation programs in the evolution of self-fertile nematodes

Self-fertility evolved independently in three species of Caenorhabditis, yet the underlying genetic changes remain unclear. This transition required that XX animals acquire the ability to produce sperm and then signal those sperm to activate and fertilise oocytes. Here, we show that all genes that regulate sperm activation in C. elegans are conserved throughout the genus, even in male/female species. By using gene editing, we show that C. elegans and C. briggsae hermaphrodites use the SPE-8 tyrosine kinase pathway to activate sperm, whereas C. tropicalis hermaphrodites use a TRY-5 serine protease pathway. Finally, our analysis of double mutants shows that these pathways were redundant in ancestral males. Thus, newly evolving hermaphrodites became self-fertile by co-opting either of the two redundant male programs. The existence of these alternatives helps explain the frequent origin of self-fertility in nematode lineages. This work also demonstrates that the new genome-editing techniques allow unprecedented power and precision in evolutionary studies.

Recurrent loss of specific introns during angiosperm evolution

Numerous instances of presence/absence variations for introns have been documented in eukaryotes, and some cases of recurrent loss of the same intron have been suggested. However, there has been no comprehensive or phylogenetically deep analysis of recurrent intron loss. Of 883 cases of intron presence/absence variation that we detected in five sequenced grass genomes, 93 were confirmed as recurrent losses and the rest could be explained by single losses (652) or single gains (118). No case of recurrent intron gain was observed. Deep phylogenetic analysis often indicated that apparent intron gains were actually numerous independent losses of the same intron. Recurrent loss exhibited extreme non-randomness, in that some introns were removed independently in many lineages. The two larger genomes, maize and sorghum, were found to have a higher rate of both recurrent loss and overall loss and/or gain than foxtail millet, rice or Brachypodium. Adjacent introns and small introns were found to be preferentially lost. Intron loss genes exhibited a high frequency of germ line or early embryogenesis expression. In addition, flanking exon A+T-richness and intron TG/CG ratios were higher in retained introns. This last result suggests that epigenetic status, as evidenced by a loss of methylated CG dinucleotides, may play a role in the process of intron loss. This study provides the first comprehensive analysis of recurrent intron loss, makes a series of novel findings on the patterns of recurrent intron loss during the evolution of the grass family, and provides insight into the molecular mechanism(s) underlying intron loss.

The spliceosomal introns are nucleotide sequences that interrupt coding regions of eukaryotic genes and are removed by RNA splicing after transcription. Recent studies have reported several examples of possible recurrent intron loss or gain, i.e., introns that are independently removed from or inserted into the identical sites more than once in an investigated phylogeny. However, the frequency, evolutionary patterns or other characteristics of recurrent intron turnover remain unknown. We provide results for the first comprehensive analysis of recurrent intron turnover within a plant family and show that recurrent intron loss represents a considerable portion of all intron losses identified and intron loss events far outnumber intron gain events. We also demonstrate that recurrent intron loss is non-random, affecting only a small number of introns that are repeatedly lost, and that different lineages show significantly different rates of intron loss. Our results suggest a possible role of DNA methylation in the process of intron loss. Moreover, this study provides strong support for the model of intron loss by reverse transcriptase mediated conversion of genes by their processed mRNA transcripts.

A mitochondrial intron in a verongid sponge

We discovered for the first time a mitochondrial intron in a non-tetillid demosponge, which sheds new light on the interpretation of mitochondrial intron evolution among non-bilaterian animals and has consequences for phylogenetic and DNA barcoding studies. The newly discovered class 1 intron of Aplysinella rhax (Verongida) CO1 has an ORF for a putative LAGLIDADG-type and resembles other sponge and cnidarian mitochondrial introns. Our analysis of the Aplysinella rhax intron underlines that the patchy distribution of introns in sponges is caused by a combination of horizontal and vertical transmission. Further implications for CO1 phylogenetic and barcoding projects are discussed.

Reverse transcriptase and intron number evolution

BACKGROUND

Introns are universal in eukaryotic genomes and play important roles in transcriptional regulation, mRNA export to the cytoplasm, nonsense-mediated decay as both a regulatory and a splicing quality control mechanism, R-loop avoidance, alternative splicing, chromatin structure, and evolution by exon-shuffling.

METHODS

Sixteen complete fungal genomes were used 13 of which were sequenced and annotated by JGI. Ustilago maydis, Cryptococcus neoformans, and Coprinus cinereus (also named Coprinopsis cinerea) were from the Broad Institute. Gene models from JGI-annotated genomes were taken from the GeneCatalog track that contained the best representative gene models. Varying fractions of the GeneCatalog were manually curated by external users. For clarity, we used the JGI unique database identifier.

RESULTS

The last common ancestor of eukaryotes (LECA) has an estimated 6.4 coding exons per gene (EPG) and evolved into the diverse eukaryotic life forms, which is recapitulated by the development of a stem cell. We found a parallel between the simulated reverse transcriptase (RT)-mediated intron loss and the comparative analysis of 16 fungal genomes that spanned a wide range of intron density. Although footprints of RT (RTF) were dynamic, relative intron location (RIL) to the 5'-end of mRNA faithfully traced RT-mediated intron loss and revealed 7.7 EPG for LECA. The mode of exon length distribution was conserved in simulated intron loss, which was exemplified by the shared mode of 75 nt between fungal and Chlamydomonas genomes. The dominant ancient exon length was corroborated by the average exon length of the most intron-rich genes in fungal genomes and consistent with ancient protein modules being ~25 aa. Combined with the conservation of a protein length of 400 aa, the earliest ancestor of eukaryotes could have 16 EPG. During earlier evolution, Ascomycota's ancestor had significantly more 3'-biased RT-mediated intron loss that was followed by dramatic RTF loss. There was a down trend of EPG from more conserved to less conserved genes. Moreover, species-specific genes have higher exon-densities, shorter exons, and longer introns when compared to genes conserved at the phylum level. However, intron length in species-specific genes became shorter than that of genes conserved in all species after genomes experiencing drastic intron loss. The estimated EPG from the most frequent exon length is more than double that from the RIL method.

CONCLUSIONS

This implies significant intron loss during the very early period of eukaryotic evolution. De novo gene-birth contributes to shorter exons, longer introns, and higher exon-density in species-specific genes relative to conserved genes.

Analysis of CACTA transposases reveals intron loss as major factor influencing their exon/intron structure in monocotyledonous and eudicotyledonous hosts

Background

CACTA elements are DNA transposons and are found in numerous organisms. Despite their low activity, several thousand copies can be identified in many genomes. CACTA elements transpose using a ‘cut-and-paste’ mechanism, which is facilitated by a DDE transposase. DDE transposases from CACTA elements contain, despite their conserved function, different exon numbers among various CACTA families. While earlier studies analyzed the ancestral history of the DDE transposases, no studies have examined exon loss and gain with a view of mechanisms that could drive the changes.

Results

We analyzed 64 transposases from different CACTA families among monocotyledonous and eudicotyledonous host species. The annotation of the exon/intron boundaries showed a range from one to six exons. A robust multiple sequence alignment of the 64 transposases based on their protein sequences was created and used for phylogenetic analysis, which revealed eight different clades. We observed that the exon numbers in CACTA transposases are not specific for a host genome. We found that ancient CACTA lineages diverged before the divergence of monocotyledons and eudicotyledons. Most exon/intron boundaries were found in three distinct regions among all the transposases, grouping 63 conserved intron/exon boundaries.

Conclusions

We propose a model for the ancestral CACTA transposase gene, which consists of four exons, that predates the divergence of the monocotyledons and eudicotyledons. Based on this model, we propose pathways of intron loss or gain to explain the observed variation in exon numbers. While intron loss appears to have prevailed, a putative case of intron gain was nevertheless observed.

The evolutionary origins and consequences of self-fertility in nematodes

Self-fertile hermaphrodites have evolved from male/female ancestors in many nematode species, and this transition occurred on three independent occasions in the genus Caenorhabditis. Genetic analyses in Caenorhabditis show that the origin of hermaphrodites required two types of changes: alterations to the sex-determination pathway that allowed otherwise female animals to make sperm during larval development, and the production of signals from the gonad that caused these sperm to activate and fertilize oocytes. Comparisons of C. elegans and C. briggsae hermaphrodites show that the ancestral sex-determination pathway has been altered in multiple unique ways. Some of these changes must have precipitated the production of sperm in XX animals, and others were modifying mutations that increased the efficiency of hermaphroditic reproduction. Reverse genetic experiments show that XX animals acquired the ability to activate sperm by co-opting one of the two redundant pathways that normally work in males. Finally, the adoption of a hermaphroditic lifestyle had profound effects on ecological and sexual interactions and genomic organization. Thus, nematode mating systems are ideal for elucidating the origin of novel traits, and studying the influence of developmental processes on evolutionary change.

Dependence of the sperm/oocyte decision on the nucleosome remodeling factor complex was acquired during recent Caenorhabditis briggsae evolution

The major families of chromatin remodelers have been conserved throughout eukaryotic evolution. Because they play broad, pleiotropic roles in gene regulation, it was not known if their functions could change rapidly. Here, we show that major alterations in the use of chromatin remodelers are possible, because the nucleosome remodeling factor (NURF) complex has acquired a unique role in the sperm/oocyte decision of the nematode Caenorhabditis briggsae. First, lowering the activity of C. briggsae NURF-1 or ISW-1, the core components of the NURF complex, causes germ cells to become oocytes rather than sperm. This observation is based on the analysis of weak alleles and null mutations that were induced with TALENs and on RNA interference. Second, qRT-polymerase chain reaction data show that the C. briggsae NURF complex promotes the expression of Cbr-fog-1 and Cbr-fog-3, two genes that control the sperm/oocyte decision. This regulation occurs in the third larval stage and affects the expression of later spermatogenesis genes. Third, double mutants reveal that the NURF complex and the transcription factor TRA-1 act independently on Cbr-fog-1 and Cbr-fog-3. TRA-1 binds both promoters, and computer analyses predict that these binding sites are buried in nucleosomes, so we suggest that the NURF complex alters chromatin structure to allow TRA-1 access to Cbr-fog-1 and Cbr-fog-3. Finally, lowering NURF activity by mutation or RNA interference does not affect this trait in other nematodes, including the sister species C. nigoni, so it must have evolved recently. We conclude that altered chromatin remodeling could play an important role in evolutionary change.

Rogue sperm indicate sexually antagonistic coevolution in nematodes

Intense reproductive competition often continues long after animals finish mating. In many species, sperm from one male compete with those from others to find and fertilize oocytes. Since this competition occurs inside the female reproductive tract, she often influences the outcome through physical or chemical factors, leading to cryptic female choice. Finally, traits that help males compete with each other are sometimes harmful to females, and female countermeasures may thwart the interests of males, which can lead to an arms race between the sexes known as sexually antagonistic coevolution. New studies from Caenorhabditis nematodes suggest that males compete with each other by producing sperm that migrate aggressively and that these sperm may be more likely to win access to oocytes. However, one byproduct of this competition appears to be an increased probability that these sperm will go astray, invading the ovary, prematurely activating oocytes, and sometimes crossing basement membranes and leaving the gonad altogether. These harmful effects are sometimes observed in crosses between animals of the same species but are most easily detected in interspecies crosses, leading to dramatically lowered fitness, presumably because the competitiveness of the sperm and the associated female countermeasures are not precisely matched. This mismatch is most obvious in crosses involving individuals from androdioecious species (which have both hermaphrodites and males), as predicted by the lower levels of sperm competition these species experience. These results suggest a striking example of sexually antagonistic coevolution and dramatically expand the value of nematodes as a laboratory system for studying postcopulatory interactions.

PCR survey of 50 introns in animals: cross-amplification of homologous EPIC loci in eight non-bilaterian, protostome and deuterostome phyla

Exon Primed Intron Crossing (EPIC) markers provide molecular tools that are susceptible to be variable within species while remaining amplifiable by PCR using potentially universal primers. In this study we tested the possibility of obtaining PCR products from 50 EPIC markers on 23 species belonging to seven different phyla (Porifera, Cnidaria, Arthropoda, Nematoda, Mollusca, Annelida, Echinodermata) using 70 new primer pairs. A previous study had identified and tested those loci in a dozen species, including another phylum, Urochordata (Chenuil et al., 2010). Results were contrasted among species. The best results were achieved with the oyster (Mollusca) where 28 loci provided amplicons susceptible to contain an intron according to their size. This was however not the case with the other mollusk Crepidula fornicata, which seems to have undergone a reduction in intron number or intron size. In the Porifera, 13 loci appeared susceptible to contain an intron, a surprisingly high number for this phylum considering its phylogenetic distance with genomic data used to design the primers. For two cnidarian species, numerous loci (24) were obtained. Ecdysozoan phyla (arthropods and nematodes) proved less successful than others as expected considering reports of their rapid rate of genome evolution and the worst results were obtained for several arthropods. Some general patterns among phyla arose, and we discuss how the results of this EPIC survey may give new insights into genome evolution of the study species. This work confirms that this set of EPIC loci provides an easy-to-use toolbox to identify genetic markers potentially useful for population genetics, phylogeography or phylogenetic studies for a large panel of metazoan species. We then argue that obtaining diploid sequence genotypes for these loci became simple and affordable owing to Next-Generation Sequencing development. Species surveyed in this study belong to several genera (Acanthaster, Alvinocaris, Aplysina, Aurelia, Crepidula, Eunicella, Hediste, Hemimysis, Litoditis, Lophelia, Mesopodopsis, Mya, Ophiocten, Ophioderma, Ostrea, Pelagia, Platynereis, Rhizostoma, Rimicaris), two of them, belonging to the family Vesicomydae and Eunicidae, could not be determined at the genus level.

Evolutionary change within a bipotential switch shaped the sperm/oocyte decision in hermaphroditic nematodes

A subset of transcription factors like Gli2 and Oct1 are bipotential--they can activate or repress the same target, in response to changing signals from upstream genes. Some previous studies implied that the sex-determination protein TRA-1 might also be bipotential; here we confirm this hypothesis by identifying a co-factor, and use it to explore how the structure of a bipotential switch changes during evolution. First, null mutants reveal that C. briggsae TRR-1 is required for spermatogenesis, RNA interference implies that it works as part of the Tip60 Histone Acetyl Transferase complex, and RT-PCR data show that it promotes the expression of Cbr-fog-3, a gene needed for spermatogenesis. Second, epistasis tests reveal that TRR-1 works through TRA-1, both to activate Cbr-fog-3 and to control the sperm/oocyte decision. Since previous studies showed that TRA-1 can repress fog-3 as well, these observations demonstrate that it is bipotential. Third, TRR-1 also regulates the development of the male tail. Since Cbr-tra-2 Cbr-trr-1 double mutants resemble Cbr-tra-1 null mutants, these two regulatory branches control all tra-1 activity. Fourth, striking differences in the relationship between these two branches of the switch have arisen during recent evolution. C. briggsae trr-1 null mutants prevent hermaphrodite spermatogenesis, but not Cbr-fem null mutants, which disrupt the other half of the switch. On the other hand, C. elegans fem null mutants prevent spermatogenesis, but not Cel-trr-1 mutants. However, synthetic interactions confirm that both halves of the switch exist in each species. Thus, the relationship between the two halves of a bipotential switch can shift rapidly during evolution, so that the same phenotype is produce by alternative, complementary mechanisms.

Mechanisms of intron loss and gain in the fission yeast Schizosaccharomyces

The fission yeast, Schizosaccharomyces pombe, is an important model species with a low intron density. Previous studies showed extensive intron losses during its evolution. To test the models of intron loss and gain in fission yeasts, we conducted a comparative genomic analysis in four Schizosaccharomyces species. Both intronization and de-intronization were observed, although both were at a low frequency. A de-intronization event was caused by a degenerative mutation in the branch site. Four cases of imprecise intron losses were identified, indicating that genomic deletion is not a negligible mechanism of intron loss. Most intron losses were precise deletions of introns, and were significantly biased to the 3′ sides of genes. Adjacent introns tended to be lost simultaneously. These observations indicated that the main force shaping the exon-intron structures of fission yeasts was precise intron losses mediated by reverse transcriptase. We found two cases of intron gains caused by tandem genomic duplication, but failed to identify the mechanisms for the majority of the intron gain events observed. In addition, we found that intron-lost and intron-gained genes had certain similar features, such as similar Gene Ontology categories and expression levels.

Splicing-related features of introns serve to propel evolution

The role of spliceosomal intronic structures played in evolution has only begun to be elucidated. Comparative genomic analyses of fungal snoRNA sequences, which are often contained within introns and/or exons, revealed that about one-third of snoRNA-associated introns in three major snoRNA gene clusters manifested polymorphisms, likely resulting from intron loss and gain events during fungi evolution. Genomic deletions can clearly be observed as one mechanism underlying intron and exon loss, as well as generation of complex introns where several introns lie in juxtaposition without intercalating exons. Strikingly, by tracking conserved snoRNAs in introns, we found that some introns had moved from one position to another by excision from donor sites and insertion into target sties elsewhere in the genome without needing transposon structures. This study revealed the origin of many newly gained introns. Moreover, our analyses suggested that intron-containing sequences were more prone to sustainable structural changes than DNA sequences without introns due to intron's ability to jump within the genome via unknown mechanisms. We propose that splicing-related structural features of introns serve as an additional motor to propel evolution.

Frequency of intron loss correlates with processed pseudogene abundance: a novel strategy to test the reverse transcriptase model of intron loss

Background

Although intron loss in evolution has been described, the mechanism involved is still unclear. Three models have been proposed, the reverse transcriptase (RT) model, genomic deletion model and double-strand-break repair model. The RT model, also termed mRNA-mediated intron loss, suggests that cDNA molecules reverse transcribed from spliced mRNA recombine with genomic DNA causing intron loss. Many studies have attempted to test this model based on its predictions, such as simultaneous loss of adjacent introns, 3'-side bias of intron loss, and germline expression of intron-lost genes. Evidence either supporting or opposing the model has been reported. The mechanism of intron loss proposed in the RT model shares the process of reverse transcription with the formation of processed pseudogenes. If the RT model is correct, genes that have produced more processed pseudogenes are more likely to undergo intron loss.

Results

In the present study, we observed that the frequency of intron loss is correlated with processed pseudogene abundance by analyzing a new dataset of intron loss obtained in mice and rats. Furthermore, we found that mRNA molecules of intron-lost genes are mostly translated on free cytoplasmic ribosomes, a feature shared by mRNA molecules of the parental genes of processed pseudogenes and long interspersed elements. This feature is likely convenient for intron-lost gene mRNA molecules to be reverse transcribed. Analyses of adjacent intron loss, 3'-side bias of intron loss, and germline expression of intron-lost genes also support the RT model.

Conclusions

Compared with previous evidence, the correlation between the abundance of processed pseudogenes and intron loss frequency more directly supports the RT model of intron loss. Exploring such a correlation is a new strategy to test the RT model in organisms with abundant processed pseudogenes.

Modeling the evolution dynamics of exon-intron structure with a general random fragmentation process

Background

Most eukaryotic genes are interrupted by spliceosomal introns. The evolution of exon-intron structure remains mysterious despite rapid advance in genome sequencing technique. In this work, a novel approach is taken based on the assumptions that the evolution of exon-intron structure is a stochastic process, and that the characteristics of this process can be understood by examining its historical outcome, the present-day size distribution of internal translated exons (exon). Through the combination of simulation and modeling the size distribution of exons in different species, we propose a general random fragmentation process (GRFP) to characterize the evolution dynamics of exon-intron structure. This model accurately predicts the probability that an exon will be split by a new intron and the distribution of novel insertions along the length of the exon.

Results

As the first observation from this model, we show that the chance for an exon to obtain an intron is proportional to its size to the 3rd power. We also show that such size dependence is nearly constant across gene, with the exception of the exons adjacent to the 5^′ UTR. As the second conclusion from the model, we show that intron insertion loci follow a normal distribution with a mean of 0.5 (center of the exon) and a standard deviation of 0.11. Finally, we show that intron insertions within a gene are independent of each other for vertebrates, but are more negatively correlated for non-vertebrate. We use simulation to demonstrate that the negative correlation might result from significant intron loss during evolution, which could be explained by selection against multi-intron genes in these organisms.

Conclusions

The GRFP model suggests that intron gain is dynamic with a higher chance for longer exons; introns are inserted into exons randomly with the highest probability at the center of the exon. GRFP estimates that there are 78 introns in every 10 kb coding sequences for vertebrate genomes, agreeing with empirical observations. GRFP also estimates that there are significant intron losses in the evolution of non-vertebrate genomes, with extreme cases of around 57% intron loss in Drosophila melanogaster, 28% in Caenorhabditis elegans, and 24% in Oryza sativa.

Gene make-up: rapid and massive intron gains after horizontal transfer of a bacterial α-amylase gene to Basidiomycetes

BACKGROUND

Increasing genome data show that introns, a hallmark of eukaryotes, already existed at a high density in the last common ancestor of extant eukaryotes. However, intron content is highly variable among species. The tempo of intron gains and losses has been irregular and several factors may explain why some genomes are intron-poor whereas other are intron-rich.

RESULTS

We studied the dynamics of intron gains and losses in an α-amylase gene, whose product breaks down starch and other polysaccharides. It was transferred from an Actinobacterium to an ancestor of Agaricomycotina. This gene underwent further duplications in several species. The results indicate a high rate of intron insertions soon after the gene settled in the fungal genome. A number of these oldest introns, regularly scattered along the gene, remained conserved. Subsequent gains and losses were lineage dependent, with a majority of losses. Moreover, a few species exhibited a high number of both specific intron gains and losses in recent periods. There was little sequence conservation around insertion sites, then probably little information for splicing, whereas splicing sites, inside introns, showed typical and conserved patterns. There was little variation of intron size.

CONCLUSIONS

Since most Basidiomycetes have intron-rich genomes and this richness was ancestral in Fungi, long before the transfer event, we suggest that the new gene was shaped to comply with requirements of the splicing machinery, such as short exon and intron sizes, in order to be correctly processed.

Control of oocyte growth and meiotic maturation in Caenorhabditis elegans

In sexually reproducing animals, oocytes arrest at diplotene or diakinesis and resume meiosis (meiotic maturation) in response to hormones. Chromosome segregation errors in female meiosis I are the leading cause of human birth defects, and age-related changes in the hormonal environment of the ovary are a suggested cause. Caenorhabditis elegans is emerging as a genetic paradigm for studying hormonal control of meiotic maturation. The meiotic maturation processes in C. elegans and mammals share a number of biological and molecular similarities. Major sperm protein (MSP) and luteinizing hormone (LH), though unrelated in sequence, both trigger meiotic resumption using somatic Gα(s)-adenylate cyclase pathways and soma-germline gap-junctional communication. At a molecular level, the oocyte responses apparently involve the control of conserved protein kinase pathways and post-transcriptional gene regulation in the oocyte. At a cellular level, the responses include cortical cytoskeletal rearrangement, nuclear envelope breakdown, assembly of the acentriolar meiotic spindle, chromosome segregation, and likely changes important for fertilization and the oocyte-to-embryo transition. This chapter focuses on signaling mechanisms required for oocyte growth and meiotic maturation in C. elegans and discusses how these mechanisms coordinate the completion of meiosis and the oocyte-to-embryo transition.

Near intron pairs and the metazoan tree

Gene structure data can substantially advance our understanding of metazoan evolution and deliver an independent approach to resolve conflicts among existing hypotheses. Here, we used changes of spliceosomal intron positions as novel phylogenetic marker to reconstruct the animal tree. This kind of data is inferred from orthologous genes containing mutually exclusive introns at pairs of sequence positions in close proximity, so-called near intron pairs (NIPs). NIP data were collected for 48 species and utilized as binary genome-level characters in maximum parsimony (MP) analyses to reconstruct deep metazoan phylogeny. All groupings that were obtained with more than 80% bootstrap support are consistent with currently supported phylogenetic hypotheses. This includes monophyletic Chordata, Vertebrata, Nematoda, Platyhelminthes and Trochozoa. Several other clades such as Deuterostomia, Protostomia, Arthropoda, Ecdysozoa, Spiralia, and Eumetazoa, however, failed to be recovered due to a few problematic taxa such as the mite Ixodesand the warty comb jelly Mnemiopsis. The corresponding unexpected branchings can be explained by the paucity of synapomorphic changes of intron positions shared between some genomes, by the sensitivity of MP analyses to long-branch attraction (LBA), and by the very unequal evolutionary rates of intron loss and intron gain during evolution of the different subclades of metazoans. In addition, we obtained an assemblage of Cnidaria, Porifera, and Placozoa as sister group of Bilateria+Ctenophora with medium support, a disputable, but remarkable result. We conclude that NIPs can be used as phylogenetic characters also within a broader phylogenetic context, given that they have emerged regularly during evolution irrespective of the large variation of intron density across metazoan genomes.

Structural and functional characterization of ribosomal protein gene introns in sponges

Ribosomal protein genes (RPGs) are a powerful tool for studying intron evolution. They exist in all three domains of life and are much conserved. Accumulating genomic data suggest that RPG introns in many organisms abound with non-protein-coding-RNAs (ncRNAs). These ancient ncRNAs are small nucleolar RNAs (snoRNAs) essential for ribosome assembly. They are also mobile genetic elements and therefore probably important in diversification and enrichment of transcriptomes through various mechanisms such as intron/exon gain/loss. snoRNAs in basal metazoans are poorly characterized. We examined 449 RPG introns, in total, from four demosponges: Amphimedon queenslandica, Suberites domuncula, Suberites ficus and Suberites pagurorum and showed that RPG introns from A. queenslandica share position conservancy and some structural similarity with "higher" metazoans. Moreover, our study indicates that mobile element insertions play an important role in the evolution of their size. In four sponges 51 snoRNAs were identified. The analysis showed discrepancies between the snoRNA pools of orthologous RPG introns between S. domuncula and A. queenslandica. Furthermore, these two sponges show as much conservancy of RPG intron positions between each other as between themselves and human. Sponges from the Suberites genus show consistency in RPG intron position conservation. However, significant differences in some of the orthologous RPG introns of closely related sponges were observed. This indicates that RPG introns are dynamic even on these shorter evolutionary time scales.

Evolutionary patterns of RNA-based gene duplicates in Caenorhabditis nematodes coincide with their genomic features

BACKGROUND

RNA-based gene duplicates (retrocopies) played pivotal roles in many physiological processes. Nowadays, functional retrocopies have been systematically identified in several mammals, fruit flies, plants, zebrafish and other chordates, etc. However, studies about this kind of duplication in Caenorhabditis nematodes have not been reported.

FINDINGS

We identified 43, 48, 43, 9, and 42 retrocopies, of which 6, 15, 18, 3, and 13 formed chimeric genes in C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei, respectively. At least 5 chimeric types exist in Caenorhabditis species, of which retrocopy recruiting both N and C terminus is the commonest one. Evidences from different analyses demonstrate many retrocopies and almost all chimeric genes may be functional in these species. About half of retrocopies in each species has coordinates in other species, and we suggest that retrocopies in closely related species may be helpful in identifying retrocopies for one certain species.

CONCLUSIONS

A number of retrocopies and chimeric genes exist in Caenorhabditis genomes, and some of them may be functional. The evolutionary patterns of these genes may correlate with their genomic features, such as the activity of retroelements, the high rate of mutation and deletion rate, and a large proportion of genes subject to trans-splicing.

SPE-44 implements sperm cell fate

The sperm/oocyte decision in the hermaphrodite germline of Caenorhabditis elegans provides a powerful model for the characterization of stem cell fate specification and differentiation. The germline sex determination program that governs gamete fate has been well studied, but direct mediators of cell-type-specific transcription are largely unknown. We report the identification of spe-44 as a critical regulator of sperm gene expression. Deletion of spe-44 causes sperm-specific defects in cytokinesis, cell cycle progression, and organelle assembly resulting in sterility. Expression of spe-44 correlates precisely with spermatogenesis and is regulated by the germline sex determination pathway. spe-44 is required for the appropriate expression of several hundred sperm-enriched genes. The SPE-44 protein is restricted to the sperm-producing germline, where it localizes to the autosomes (which contain sperm genes) but is excluded from the transcriptionally silent X chromosome (which does not). The orthologous gene in other Caenorhabditis species is similarly expressed in a sex-biased manner, and the protein likewise exhibits autosome-specific localization in developing sperm, strongly suggestive of an evolutionarily conserved role in sperm gene expression. Our analysis represents the first identification of a transcriptional regulator whose primary function is the control of gamete-type-specific transcription in this system.

Stem cells give rise to the variety of specialized cell types within an organism. The decision to adopt a particular cell fate, a process known as specification or determination, requires the coordinated expression of all of the genes needed for that specialized cell to develop and function properly. Understanding the mechanisms that govern these patterns of gene expression is critical to our understanding of stem cell fate specification. We study this process in a nematode species that makes both sperm and eggs from the same stem cell population. We have identified a gene, named spe-44, that is required for the proper expression of sperm genes (but not egg genes). Mutation in spe-44 produces sterile sperm with developmental defects. spe-44 is controlled by factors that govern the sperm/egg decision, and its function in controlling sperm gene expression appears to be conserved in other nematode species.

Context-dependent function of a conserved translational regulatory module

The modification of transcriptional regulation is a well-documented evolutionary mechanism in both plants and animals, but post-transcriptional controls have received less attention. The derived hermaphrodite of C. elegans has regulated spermatogenesis in an otherwise female body. The PUF family RNA-binding proteins FBF-1 and FBF-2 limit XX spermatogenesis by repressing the male-promoting proteins FEM-3 and GLD-1. Here, we examine the function of PUF homologs from other Caenorhabditis species, with emphasis on C. briggsae, which evolved selfing convergently. C. briggsae lacks a bona fide fbf-1/2 ortholog, but two members of the related PUF-2 subfamily, Cbr-puf-2 and Cbr-puf-1.2, do have a redundant germline sex determination role. Surprisingly, this is to promote, rather than limit, hermaphrodite spermatogenesis. We provide genetic, molecular and biochemical evidence that Cbr-puf-2 and Cbr-puf-1.2 repress Cbr-gld-1 by a conserved mechanism. However, Cbr-gld-1 acts to limit, rather than promote, XX spermatogenesis. As with gld-1, no sex determination function for fbf or puf-2 orthologs is observed in gonochoristic Caenorhabditis. These results indicate that PUF family genes were co-opted for sex determination in each hermaphrodite via their long-standing association with gld-1, and that their precise sex-determining roles depend on the species-specific context in which they act. Finally, we document non-redundant roles for Cbr-puf-2 in embryonic and early larval development, the latter role being essential. Thus, recently duplicated PUF paralogs have already acquired distinct functions.

Independent recruitments of a translational regulator in the evolution of self-fertile nematodes

Pleiotropic developmental regulators have been repeatedly linked to the evolution of anatomical novelties. Known mechanisms include cis-regulatory DNA changes that alter regulator transcription patterns or modify target-gene linkages. Here, we examine the role of another form of regulation, translational control, in the repeated evolution of self-fertile hermaphroditism in Caenorhabditis nematodes. Caenorhabditis elegans hermaphrodites initiate spermatogenesis in an otherwise female body through translational repression of the gene tra-2. This repression is mediated by GLD-1, an RNA-binding protein also required for oocyte meiosis and differentiation. By contrast, we show that in the convergently hermaphroditic Caenorhabditis briggsae, GLD-1 acts to promote oogenesis. The opposite functions of gld-1 in these species are not gene-intrinsic, but instead result from the unique contexts for its action that evolved in each. In C. elegans, GLD-1 became essential for promoting XX spermatogenesis via changes in the tra-2 mRNA and evolution of the species-specific protein FOG-2. C. briggsae GLD-1 became an essential repressor of sperm-promoting genes, including Cbr-puf-8, and did not evolve a strong association with tra-2. Despite its variable roles in sex determination, the function of gld-1 in female meiotic progression is ancient and conserved. This conserved role may explain why gld-1 is repeatedly recruited to regulate hermaphroditism. We conclude that, as with transcription factors, spatially localized translational regulators play important roles in the evolution of anatomical novelties.

TRY-5 is a sperm-activating protease in Caenorhabditis elegans seminal fluid

Seminal fluid proteins have been shown to play important roles in male reproductive success, but the mechanisms for this regulation remain largely unknown. In Caenorhabditis elegans, sperm differentiate from immature spermatids into mature, motile spermatozoa during a process termed sperm activation. For C. elegans males, sperm activation occurs during insemination of the hermaphrodite and is thought to be mediated by seminal fluid, but the molecular nature of this activity has not been previously identified. Here we show that TRY-5 is a seminal fluid protease that is required in C. elegans for male-mediated sperm activation. We observed that TRY-5::GFP is expressed in the male somatic gonad and is transferred along with sperm to hermaphrodites during mating. In the absence of TRY-5, male seminal fluid loses its potency to transactivate hermaphrodite sperm. However, TRY-5 is not required for either hermaphrodite or male fertility, suggesting that hermaphrodite sperm are normally activated by a distinct hermaphrodite-specific activator to which male sperm are also competent to respond. Within males, TRY-5::GFP localization within the seminal vesicle is antagonized by the protease inhibitor SWM-1. Together, these data suggest that TRY-5 functions as an extracellular activator of C. elegans sperm. The presence of TRY-5 within the seminal fluid couples the timing of sperm activation to that of transfer of sperm into the hermaphrodite uterus, where motility must be rapidly acquired. Our results provide insight into how C. elegans has adopted sex-specific regulation of sperm motility to accommodate its male-hermaphrodite mode of reproduction.

A bias caused by ectopic development produces sexually dimorphic sperm in nematodes

Self-fertile hermaphrodites have evolved independently several times in the genus Caenorhabditis [1, 2]. These XX hermaphrodites make smaller sperm than males [3, 4], which they use to fertilize their own oocytes. Because larger sperm outcompete smaller sperm in nematodes [3-5], it had been assumed that this dimorphism evolved in response to sperm competition. However, we show that it was instead caused by a developmental bias. When we transformed females of the species Caenorhabditis remanei into hermaphrodites [6], their sperm were significantly smaller than those of males. Because this species never makes hermaphrodites in the wild, this dimorphism cannot be due to selection. Instead, analyses of the related nematode Caenorhabditis elegans suggest that this dimorphism might reflect the development of sperm within the distinct physiological environment of the hermaphrodite gonad. These results reveal a new mechanism for some types of developmental bias-the effects of a novel physical location alter the development of ectopic cells in predictable ways.

The role of reverse transcriptase in intron gain and loss mechanisms

Intron density is highly variable across eukaryotic species. It seems that different lineages have experienced considerably different levels of intron gain and loss events, but the reasons for this are not well known. A large number of mechanisms for intron loss and gain have been suggested, and most of them have at least some level of indirect support. We therefore figured out that the variability in intron density can be a reflection of the fact that different mechanisms are active in different lineages. Quite a number of these putative mechanisms, both for intron loss and for intron gain, postulate that the enzyme reverse transcriptase (RT) has a key role in the process. In this paper, we lay out three predictions whose approval or falsification gives indication for the involvement of RT in intron gain and loss processes. Testing these predictions requires data on the intron gain and loss rates of individual genes along different branches of the eukaryotic phylogenetic tree. So far, such rates could not be computed, and hence, these predictions could not be rigorously evaluated. Here, we use a maximum likelihood algorithm that we have devised in the past, Evolutionary Reconstruction by Expectation Maximization, which allows the estimation of such rates. Using this algorithm, we computed the intron loss and gain rates of more than 300 genes in each branch of the phylogenetic tree of 19 eukaryotic species. Based on that we found only little support for RT activity in intron gain. In contrast, we suggest that RT-mediated intron loss is a mechanism that is very efficient in removing introns, and thus, its levels of activity may be a major determinant of intron number. Moreover, we found that intron gain and loss rates are negatively correlated in intron-poor species but are positively correlated for intron-rich species. One explanation to this is that intron gain and loss mechanisms in intron-rich species (like metazoans) share a common mechanistic component, albeit not a RT.

Phylogenetic distribution of intron positions in alpha-amylase genes of bilateria suggests numerous gains and losses

Most eukaryotes have at least some genes interrupted by introns. While it is well accepted that introns were already present at moderate density in the last eukaryote common ancestor, the conspicuous diversity of intron density among genomes suggests a complex evolutionary history, with marked differences between phyla. The question of the rates of intron gains and loss in the course of evolution and factors influencing them remains controversial. We have investigated a single gene family, alpha-amylase, in 55 species covering a variety of animal phyla. Comparison of intron positions across phyla suggests a complex history, with a likely ancestral intronless gene undergoing frequent intron loss and gain, leading to extant intron/exon structures that are highly variable, even among species from the same phylum. Because introns are known to play no regulatory role in this gene and there is no alternative splicing, the structural differences may be interpreted more easily: intron positions, sizes, losses or gains may be more likely related to factors linked to splicing mechanisms and requirements, and to recognition of introns and exons, or to more extrinsic factors, such as life cycle and population size. We have shown that intron losses outnumbered gains in recent periods, but that "resets" of intron positions occurred at the origin of several phyla, including vertebrates. Rates of gain and loss appear to be positively correlated. No phase preference was found. We also found evidence for parallel gains and for intron sliding. Presence of introns at given positions was correlated to a strong protosplice consensus sequence AG/G, which was much weaker in the absence of intron. In contrast, recent intron insertions were not associated with a specific sequence. In animal Amy genes, population size and generation time seem to have played only minor roles in shaping gene structures.

In silico characterization and comparative genomic analysis of the Culex quinquefasciatus glutathione S-transferase (GST) supergene family

The glutathione S-transferases (GSTs) are phase II class of detoxification enzymes that are involved both directly and indirectly in insecticide resistance mechanisms. The Culex quinquefasciatus GST superfamily was analyzed by utilizing the public domain Culex genome sequence. In total, 35 cytosolic (seven classes) and 5 microsomal putatively active GSTs were retrieved, classified, and annotated. The study revealed the presence of three unclassified GSTs. Of 35 cytosolic GSTs, 65% contributed by insect specific Delta-Epsilon classes. Gene cluster analysis revealed that most of the genes of Delta, Epsilon, and Theta classes were organized into gene clusters. The gene organization analysis revealed the dominance of phase "0" introns in the Culex GST family. The studies on intron loss and gain events revealed that the Delta GSTs have experienced a higher number of loss and gains during their evolution. A positive correlation was observed between the phylogenetic relationship of members of the GST superfamily and their corresponding exon-intron organization. In addition, the genes within the gene clusters revealed the monophyletic phylogenetic relationship implying the importance of gene duplication events in the gene families' evolution. Finally, the comparative genomic analysis has shown a complex evolutionary scenario associated with the GST supergene family evolution in insects.