Genetically Distinct Behavioral Modules Underlie Natural Variation in Thermal Performance Curves

Thermal reaction norms pervade organismal traits as stereotyped responses to temperature, a fundamental environmental input into sensory and physiological systems. Locomotory behavior represents an especially plastic read-out of animal response, with its dynamic dependence on environmental stimuli presenting a challenge for analysis and for understanding the genomic architecture of heritable variation. Here we characterize behavioral reaction norms as thermal performance curves for the nematode Caenorhabditis briggsae, using a collection of 23 wild isolate genotypes and 153 recombinant inbred lines to quantify the extent of genetic and plastic variation in locomotory behavior to temperature changes. By reducing the dimensionality of the multivariate phenotypic response with a function-valued trait framework, we identified genetically distinct behavioral modules that contribute to the heritable variation in the emergent overall behavioral thermal performance curve. Quantitative trait locus mapping isolated regions on Chromosome II associated with locomotory activity at benign temperatures and Chromosome V loci related to distinct aspects of sensitivity to high temperatures, with each quantitative trait locus explaining up to 28% of trait variation. These findings highlight how behavioral responses to environmental inputs as thermal reaction norms can evolve through independent changes to genetically distinct modular components of such complex phenotypes.

REPRODUCTIVE ISOLATION IN RHABDITIDAE (NEMATODA: SECERNENTEA); MECHANISMS THAT ISOLATE SIX SPECIES OF THREE GENERA

We have attempted interspecific hybridizations among six species of rhabditid nematodes: Caenorhabditis elegans, Caenorhabditis briggsae, Caenorhabditis remanei, Caenorhabditis sp. v, Rhabditis sp., and Pelodera teres. Copulation was observed in all crosses between Caenorhabditis species; however, none resulted in the generation of stable hybrid populations. No copulation was observed in crosses between Caenorhabditis males and Rhabditis or Pelodera females, even when congeneric females were present, suggesting that Caenorhabditis males are able to selectively recognize congeneric females by a short-range stimulus. All pairwise combinations of Caenorhabditis species were isolated to some degree by gametic mechanisms; 7 of 12 combinations were cross infertile and 5 of 12 were cross-fertile but had low brood sizes. In cross-fertile combinations, most hybrid embryos were inviable and arrested prior to gastrulation. Only in crosses of C. briggsae males to C. sp. v females did any hybrids survive embryogenesis. Most of these C. briggsae/C. sp. v hybrids arrested during larval development, and the few that reached adulthood invariably were female. These results are consistent with the presence of at least two lethal factors in the C. briggsae-C. sp. v combination: a maternal lethal factor in the cytoplasm of C. briggsae and a recessive lethal factor on the X chromosome of C. sp. v.

Insights into species divergence and the evolution of hermaphroditism from fertile interspecies hybrids of Caenorhabditis nematodes

The architecture of both phenotypic variation and reproductive isolation are important problems in evolutionary genetics. The nematode genus Caenorhabditis includes both gonochoristic (male/female) and androdioecious (male/hermaprodite) species. However, the natural genetic variants distinguishing reproductive mode remain unknown, and nothing is known about the genetic basis of postzygotic isolation in the genus. Here we describe the hybrid genetics of the first Caenorhabditis species pair capable of producing fertile hybrid progeny, the gonochoristic Caenorhabditis sp. 9 and the androdioecious C. briggsae. Though many interspecies F(1) arrest during embryogenesis, a viable subset develops into fertile females and sterile males. Reciprocal parental crosses reveal asymmetry in male-specific viability, female fertility, and backcross viability. Selfing and spermatogenesis are extremely rare in XX F(1), and almost all hybrid self-progeny are inviable. Consistent with this, F(1) females do not express male-specific molecular germline markers. We also investigated three approaches to producing hybrid hermaphrodites. A dominant mutagenesis screen for self-fertile F(1) hybrids was unsuccessful. Polyploid F(1) hybrids with increased C. briggsae genomic material did show elevated rates of selfing, but selfed progeny were mostly inviable. Finally, the use of backcrosses to render the hybrid genome partial homozygous for C. briggsae alleles did not increase the incidence of selfing or spermatogenesis relative to the F(1) generation. These hybrid animals were genotyped at 23 loci, and significant segregation distortion (biased against C. briggsae) was detected at 13 loci. This, combined with an absence of productive hybrid selfing, prevents formulation of simple hypotheses about the genetic architecture of hermaphroditism. In the near future, this hybrid system will likely be fruitful for understanding the genetics of reproductive isolation in Caenorhabditis.

A toolkit for rapid gene mapping in the nematode Caenorhabditis briggsae

BACKGROUND

The nematode C. briggsae serves as a useful model organism for comparative analysis of developmental and behavioral processes. The amenability of C. briggsae to genetic manipulations and the availability of its genome sequence have prompted researchers to study evolutionary changes in gene function and signaling pathways. These studies rely on the availability of forward genetic tools such as mutants and mapping markers.

RESULTS

We have computationally identified more than 30,000 polymorphisms (SNPs and indels) in C. briggsae strains AF16 and HK104. These include 1,363 SNPs that change restriction enzyme recognition sites (snip-SNPs) and 638 indels that range between 7 bp and 2 kb. We established bulk segregant and single animal-based PCR assay conditions and used these to test 107 polymorphisms. A total of 75 polymorphisms, consisting of 14 snip-SNPs and 61 indels, were experimentally confirmed with an overall success rate of 83%. The utility of polymorphisms in genetic studies was demonstrated by successful mapping of 12 mutations, including 5 that were localized to sub-chromosomal regions. Our mapping experiments have also revealed one case of a misassembled contig on chromosome 3.

CONCLUSIONS

We report a comprehensive set of polymorphisms in C. briggsae wild-type strains and demonstrate their use in mapping mutations. We also show that molecular markers can be useful tools to improve the C. briggsae genome sequence assembly. Our polymorphism resource promises to accelerate genetic and functional studies of C. briggsae genes.

Comparison of C. elegans and C. briggsae genome sequences reveals extensive conservation of chromosome organization and synteny

To determine whether the distinctive features of Caenorhabditis elegans chromosomal organization are shared with the C. briggsae genome, we constructed a single nucleotide polymorphism-based genetic map to order and orient the whole genome shotgun assembly along the six C. briggsae chromosomes. Although these species are of the same genus, their most recent common ancestor existed 80-110 million years ago, and thus they are more evolutionarily distant than, for example, human and mouse. We found that, like C. elegans chromosomes, C. briggsae chromosomes exhibit high levels of recombination on the arms along with higher repeat density, a higher fraction of intronic sequence, and a lower fraction of exonic sequence compared with chromosome centers. Despite extensive intrachromosomal rearrangements, 1:1 orthologs tend to remain in the same region of the chromosome, and colinear blocks of orthologs tend to be longer in chromosome centers compared with arms. More strikingly, the two species show an almost complete conservation of synteny, with 1:1 orthologs present on a single chromosome in one species also found on a single chromosome in the other. The conservation of both chromosomal organization and synteny between these two distantly related species suggests roles for chromosome organization in the fitness of an organism that are only poorly understood presently.

INBREEDING AND OUTBREEDING DEPRESSION IN CAENORHABDITIS NEMATODES

The nematode Caenorhabditis elegans reproduces primarily by self-fertilization of hermaphrodites, yet males are present at low frequencies in natural populations (androdioecy). The ancestral state of C. elegans was probably gonochorism (separate males and females), as in its relative C. remanei. Males may be maintained in C. elegans because outcrossed individuals escape inbreeding depression. The level of inbreeding depression is, however, expected to be low in such a highly selfing species, compared with an outcrosser like C. remanei. To investigate these issues, we measured life-history traits in the progeny of inbred versus outcrossed C. elegans and C. remanei individuals derived from recently isolated natural populations. In addition, we maintained inbred lines of C. remanei through 13 generations of full-sibling mating. Highly inbred C. remanei showed dramatic reductions in brood size and relative fitness compared to outcrossed individuals, with evidence of both direct genetic and maternal-effect inbreeding depression. This decline in fitness accumulated over time, causing extinction of nearly 90% of inbred lines, with no evidence of purging of deleterious mutations from the remaining lines. In contrast, pure strains of C. elegans performed better than crosses between strains, indicating outbreeding depression. The results are discussed in relation to the evolution of androdioecy and the effect of mating system on the level of inbreeding depression.

High nucleotide polymorphism and rapid decay of linkage disequilibrium in wild populations of Caenorhabditis remanei

The common ancestor of the self-fertilizing nematodes Caenorhabditis elegans and C. briggsae must have reproduced by obligate outcrossing, like most species in this genus. However, we have only a limited understanding about how genetic variation is patterned in such male-female (gonochoristic) Caenorhabditis species. Here, we report results from surveying nucleotide variation of six nuclear loci in a broad geographic sample of wild isolates of the gonochoristic C. remanei. We find high levels of diversity in this species, with silent-site diversity averaging 4.7%, implying an effective population size close to 1 million. Additionally, the pattern of polymorphisms reveals little evidence for population structure or deviation from neutral expectations, suggesting that the sampled C. remanei populations approximate panmixis and demographic equilibrium. Combined with the observation that linkage disequilibrium between pairs of polymorphic sites decays rapidly with distance, this suggests that C. remanei will provide an excellent system for identifying the genetic targets of natural selection from deviant patterns of polymorphism and linkage disequilibrium. The patterns revealed in this obligately outcrossing species may provide a useful model of the evolutionary circumstances in C. elegans' gonochoristic progenitor. This will be especially important if self-fertilization evolved recently in C. elegans history, because most of the evolutionary time separating C. elegans from its known relatives would have occurred in a state of obligate outcrossing.

Reproductive isolation in Caenorhabditis: terminal phenotypes of hybrid embryos

Several interspecific combinations of the "elegans" group of Caenorhabditis species are cross-fertile. Most F1 hybrids from these crosses arrest during embryogenesis. Developmental defects observed in hybrid embryos include defects in gastrulation initiation, defects in embryonic compaction, and defects in embryonic elongation. These reproductive barriers have arisen multiple times in the evolution of Caenorhabditis.

TGFbeta-like signaling and spicule development in Caenorhabditis elegans

A TGFbeta-like signal is required for spicule development in Caenorhabditis elegans males. This signal appears to originate in the male-specific musculature and is required for the migrations of cells within the proctodeum. The migrations of these cells form cellular molds, the spicule traces, in which the cuticle of the spicules is secreted. Mutations in daf-4, sma-2, sma-3, and sma-4, which disrupt TGFbeta-like signaling, result in aberrant migrations and morphologically abnormal spicules. daf-4, and hence the TGFbeta-like signal, is required prior to or during cell migrations. Therefore, the TGFbeta-like signal may act to prime the migrating cells or as a guidance cue. Mutations in lin-31 result in identical cell migration and spicule morphology defects. Thus, lin-31, which encodes a "winged helix" protein (Miller et al., Genes Dev. 7, 933-947, 1993), may be a component of this TGFbeta-like signaling pathway.

Caenorhabditis elegans genes sma-2, sma-3, and sma-4 define a conserved family of transforming growth factor beta pathway components

Although transforming growth factor beta (TGF-beta) superfamily ligands play critical roles in diverse developmental processes, how cells transduce signals from these ligands is still poorly understood. Cell surface receptors for these ligands have been identified, but their cytoplasmic targets are unknown. We have identified three Caenorhabditis elegans genes, sma-2, sma-3, and sma-4, that have mutant phenotypes similar to those of the TGF-beta-like receptor gene daf-4, indicating that they are required for daf-4-mediated developmental processes. We show that sma-2 functions in the same cells as daf-4, consistent with a role in transducing signals from the receptor. These three genes define a protein family, the dwarfins, that includes the Mad gene product, which participates in the decapentaplegic TGF-beta-like pathway in Drosophila [Sekelsky, J. J., Newfeld, S. J., Raftery, L. A., Chartoff, E. H. & Gelbart, W. M. (1995) Genetics 139, 1347-1358]. The identification of homologous components of these pathways in distantly related organisms suggests that dwarfins may be universally required for TGF-beta-like signal transduction. In fact, we have isolated highly conserved dwarfins from vertebrates, indicating that these components are not idiosyncratic to invertebrates. These analyses suggest that dwarfins are conserved cytoplasmic signal transducers.

Pattern formation in the nematode epidermis: determination of the arrangement of peripheral sense organs in the C. elegans male tail

The developmental process that determines the arrangement of ray sensilla in the Caenorhabditis elegans male tail has been studied. It is shown that the adult arrangement of rays is determined by the placement of ray cells at specific sites in the epidermis of the last larval (L4) stage. Placement of ray cells at specific epidermal sites results from the generation of neurons and support cells in the epidermis near to their final positions, and the subsequent refinement of these positions by an active mechanism involving specific cellular associations. Positions of ray cells and adjacent epidermal cells have been studied during ray development by means of indirect immunofluorescence staining with an antibody to a cell junctional antigen. Mutations are described in six genes that alter the adult arrangement of the rays, frequently resulting in fusion of rays. Changes in the adult pattern of rays in mutants appear to result from prior changes in the epidermal positions of ray cells, and for two mutants it is suggested that this may be due to the inappropriate clustering of processes from neurons and support cells of adjacent rays. Development of the wild-type arrangement of rays appears to require the specification of molecular differences between the rays that affect the specificity of their cellular associations.

Differential DNA amplification and copy number control in the hypotrichous ciliate Euplotes crassus

During macronuclear development in hypotrichous ciliated protozoans, several thousand macronuclear DNA molecules are amplified several-hundred fold. We investigated the regulation of this amplification by determining the copy numbers of three different macronuclear DNA molecules in the hypotrichous ciliate Euplotes crassus. Two of the macronuclear DNA molecules were present in approximately 1,000 copies per cell, while the third was present in approximately 6,500 copies per cell. These reiteration levels were achieved either during macronuclear development, or shortly thereafter, and were maintained during vegetative growth. The most abundant macronuclear DNA molecule is present as a single-copy sequence in the micronuclear genome. Thus, its high copy number results from differential amplification. These results indicate that DNA amplification during macronuclear development is regulated individually for each macronuclear DNA molecule.

Properties of a class of genes required for ray morphogenesis in Caenorhabditis elegans

We have identified eight mutations that define at least five terminal differentiation genes (ram genes) whose products are required during the extension of the male-specific ray sensilla in Caenorhabditis elegans. ram gene mutations result in morphological abnormalities in the sensory rays but do not appear to interfere with ray functions. A similar ray morphology phenotype was observed in males harboring mutations in three previously defined genes, dpy-11, dpy-18 and sqt-1, that also affect body shape. One of these genes, sqt-1, is known to encode a collagen. Mutations in different ram genes failed to complement, from which we infer that their gene products functionally interact. For one ram gene, failure to complement was shown to result from haploinsufficiency. Intergenic noncomplementation did not extend to the body morphology genes. The temperature-sensitive periods of both ram and body morphology mutations corresponded to the period of development in which ray extension occurs. We propose that ram gene products act together in a critical interaction between the rays and the cuticle required for wild-type ray morphology.

Micronuclear genome organization in Euplotes crassus: a transposonlike element is removed during macronuclear development

After mating, hypotrichous ciliated protozoa transform a set of their micronuclear chromosomes into thousands of short, linear DNA molecules that form the macronuclear genome. To examine micronuclear genome organization in the hypotrich Euplotes crassus, we have analyzed two cloned segments of micronuclear DNA as well as the macronuclear DNA molecules that are derived from them. E. crassus was found to display a number of features characteristic of other hypotrich genomes, including (i) clustering and close spacing of the precursors of macronuclear DNA molecules, (ii) the frequent occurrence of internal eliminated sequences within macronuclear precursors, (iii) overlapping macronuclear precursors, (iv) lack of telomeric repeats at the ends of macronuclear precursors, and (v) alternative processing of the micronuclear chromosome to yield multiple macronuclear DNA molecules. In addition, a moderately repetitive, transposonlike element that interrupts the precursors of two macronuclear DNA molecules has been identified and characterized. This transposonlike element, designated Tec1, is shown to be reproducibly removed from one of the macronuclear precursors during independent episodes of macronuclear development.

Characterization of chromosome fragmentation in two protozoans and identification of a candidate fragmentation sequence in Euplotes crassus

Following the sexual cycle, hypotrichous ciliated protozoans fragment a set of their micronuclear chromosomes to generate the thousands of short, linear DNA molecules present in the transcriptionally active macronucleus. We have used a hybrid selection procedure to examine macronuclear DNA molecules for subtelomeric length heterogeneity to determine whether chromosome fragmentation occurs at unique or multiple sites. The results suggest that multiple, but closely spaced, chromosome fragmentation sites are used by Oxytricha nova. In contrast, Euplotes crassus uses unique chromosome fragmentation sites in a reproducible manner to generate the ends of macronuclear DNA molecules. Additional studies compared DNA sequences in the vicinity of chromosome fragmentation sites in an attempt to define cis-acting sequences that direct the fragmentation process. A conserved sequence was found near chromosome fragmentation sites in E. crassus. The location of the conserved sequence suggests that chromosome fragmentation involves staggered cuts of the micronuclear DNA molecules.

Genetic characterization and use of a restriction fragment length variant in the hypotrichous ciliate Euplotes crassus

Two forms of a macronuclear DNA molecule differing in the presence or absence of a restriction endonuclease recognition site have been detected in the hypotrichous ciliate Euplotes crassus. Through a series of genetic crosses the two forms were shown to be allelic, being derived from a single micronuclear genetic locus. This restriction fragment length variant (RFLV) was used as a genetic marker to determine that the migratory and stationary pronuclei generated during mating can be genetically non-identical. In addition, the RFLV was used to investigate the efficiency of processing of the alternate alleles during macronuclear development and their subsequent transmission during vegetative growth. Little or no bias in the processing and/or amplification of the two alleles was observed during macronuclear development. During vegetative growth, however, changes in the relative amounts of the two alleles were observed.