Mitochondrial DNA Variant in COX1 Subunit Significantly Alters Energy Metabolism of Geographically Divergent Wild Isolates in Caenorhabditis elegans

Mitochondrial DNA (mtDNA) sequence variation can influence the penetrance of complex diseases and climatic adaptation. While studies in geographically defined human populations suggest that mtDNA mutations become fixed when they have conferred metabolic capabilities optimally suited for a specific environment, it has been challenging to definitively assign adaptive functions to specific mtDNA sequence variants in mammals. We investigated whether mtDNA genome variation functionally influences Caenorhabditis elegans wild isolates of distinct mtDNA lineages and geographic origins. We found that, relative to N2 (England) wild-type nematodes, CB4856 wild isolates from a warmer native climate (Hawaii) had a unique p.A12S amino acid substitution in the mtDNA-encoded COX1 core catalytic subunit of mitochondrial complex IV (CIV). Relative to N2, CB4856 worms grown at 20 °C had significantly increased CIV enzyme activity, mitochondrial matrix oxidant burden, and sensitivity to oxidative stress but had significantly reduced lifespan and mitochondrial membrane potential. Interestingly, mitochondrial membrane potential was significantly increased in CB4856 grown at its native temperature of 25 °C. A transmitochondrial cybrid worm strain, chpIR (M, CB4856 > N2), was bred as homoplasmic for the CB4856 mtDNA genome in the N2 nuclear background. The cybrid strain also displayed significantly increased CIV activity, demonstrating that this difference results from the mtDNA-encoded p.A12S variant. However, chpIR (M, CB4856 > N2) worms had significantly reduced median and maximal lifespan relative to CB4856, which may relate to their nuclear–mtDNA genome mismatch. Overall, these data suggest that C. elegans wild isolates of varying geographic origins may adapt to environmental challenges through mtDNA variation to modulate critical aspects of mitochondrial energy metabolism.     

A genome-wide library of CB4856/N2 introgression lines of Caenorhabditis elegans

Recombinant inbred lines (RILs) derived from Caenorhabditis elegans wild-type N2 and CB4856 are increasingly being used for mapping genes underlying complex traits. To speed up mapping and gene discovery, introgression lines (ILs) offer a powerful tool for more efficient QTL identification. We constructed a library of 90 ILs, each carrying a single homozygous CB4856 genomic segment introgressed into the genetic background of N2. The ILs were genotyped by 123 single-nucleotide polymorphism (SNP) markers. The proportion of the CB4856 segments in most lines does not exceed 3%, and together the introgressions cover 96% of the CB4856 genome. The value of the IL library was demonstrated by identifying novel loci underlying natural variation in two ageing-related traits, i.e. lifespan and pharyngeal pumping rate. Bin mapping of lifespan resulted in six QTLs, which all have a lifespan-shortening effect on the CB4856 allele. We found five QTLs for the decrease in pumping rate, of which four colocated with QTLs found for average lifespan. This suggests pleiotropic or closely linked QTL associated with lifespan and pumping rate. Overall, the presented IL library provides a versatile resource toward easier and efficient fine mapping and functional analyses of loci and genes underlying complex traits in C. elegans.     

Rapid high resolution single nucleotide polymorphism-comparative genome hybridization mapping in Caenorhabditis elegans

We have developed a significantly improved and simplified method for high-resolution mapping of phenotypic traits in Caenorhabditis elegans using a combination of single nucleotide polymorphisms (SNPs) and oligo array comparative genome hybridization (array CGH). We designed a custom oligonucleotide array using a subset of confirmed SNPs between the canonical wild-type Bristol strain N2 and the Hawaiian isolate CB4856, populated with densely overlapping 50-mer probes corresponding to both N2 and CB4856 SNP sequences. Using this method a mutation can be mapped to a resolution of ~200 kb in a single genetic cross. Six mutations representing each of the C. elegans chromosomes were detected unambiguously and at high resolution using genomic DNA from populations derived from as few as 100 homozygous mutant segregants of mutant N2/CB4856 heterozygotes. Our method completely dispenses with the PCR, restriction digest, and gel analysis of standard SNP mapping and should be easy to extend to any organism with interbreeding strains. This method will be particularly powerful when applied to difficult or hard-to-map low-penetrance phenotypes. It should also be possible to map polygenic traits using this method.     

Genetic mapping of variation in dauer larvae development in growing populations of Caenorhabditis elegans

In the nematode Caenorhabditis elegans, the appropriate induction of dauer larvae development within growing populations is likely to be a primary determinant of genotypic fitness. The underlying genetic architecture of natural genetic variation in dauer formation has, however, not been thoroughly investigated. Here, we report extensive natural genetic variation in dauer larvae development within growing populations across multiple wild isolates. Moreover, bin mapping of introgression lines (ILs) derived from the genetically divergent isolates N2 and CB4856 reveals 10 quantitative trait loci (QTLs) affecting dauer formation. Comparison of individual ILs to N2 identifies an additional eight QTLs, and sequential IL analysis reveals six more QTLs. Our results also show that a behavioural, laboratory-derived, mutation controlled by the neuropeptide Y receptor homolog npr-1 can affect dauer larvae development in growing populations. These findings illustrate the complex genetic architecture of variation in dauer larvae formation in C. elegans and may help to understand how the control of variation in dauer larvae development has evolved.     

Recombinational Landscape and Population Genomics of Caenorhabditis elegans

Recombination rate and linkage disequilibrium, the latter a function of population genomic processes, are the critical parameters for mapping by linkage and association, and their patterns in Caenorhabditis elegans are poorly understood. We performed high-density SNP genotyping on a large panel of recombinant inbred advanced intercross lines (RIAILs) of C. elegans to characterize the landscape of recombination and, on a panel of wild strains, to characterize population genomic patterns. We confirmed that C. elegans autosomes exhibit discrete domains of nearly constant recombination rate, and we show, for the first time, that the pattern holds for the X chromosome as well. The terminal domains of each chromosome, spanning about 7% of the genome, exhibit effectively no recombination. The RIAILs exhibit a 5.3-fold expansion of the genetic map. With median marker spacing of 61 kb, they are a powerful resource for mapping quantitative trait loci in C. elegans. Among 125 wild isolates, we identified only 41 distinct haplotypes. The patterns of genotypic similarity suggest that some presumed wild strains are laboratory contaminants. The Hawaiian strain, CB4856, exhibits genetic isolation from the remainder of the global population, whose members exhibit ample evidence of intercrossing and recombining. The population effective recombination rate, estimated from the pattern of linkage disequilibrium, is correlated with the estimated meiotic recombination rate, but its magnitude implies that the effective rate of outcrossing is extremely low, corroborating reports of selection against recombinant genotypes. Despite the low population, effective recombination rate and extensive linkage disequilibrium among chromosomes, which are techniques that account for background levels of genomic similarity, permit association mapping in wild C. elegans strains.     

A Caenorhabditis elegans Wild Type Defies the Temperature–Size Rule Owing to a Single Nucleotide Polymorphism in tra-3

Ectotherms rely for their body heat on surrounding temperatures. A key question in biology is why most ectotherms mature at a larger size at lower temperatures, a phenomenon known as the temperature–size rule. Since temperature affects virtually all processes in a living organism, current theories to explain this phenomenon are diverse and complex and assert often from opposing assumptions. Although widely studied, the molecular genetic control of the temperature–size rule is unknown. We found that the Caenorhabditis elegans wild-type N2 complied with the temperature–size rule, whereas wild-type CB4856 defied it. Using a candidate gene approach based on an N2 × CB4856 recombinant inbred panel in combination with mutant analysis, complementation, and transgenic studies, we show that a single nucleotide polymorphism in tra-3 leads to mutation F96L in the encoded calpain-like protease. This mutation attenuates the ability of CB4856 to grow larger at low temperature. Homology modelling predicts that F96L reduces TRA-3 activity by destabilizing the DII-A domain. The data show that size adaptation of ectotherms to temperature changes may be less complex than previously thought because a subtle wild-type polymorphism modulates the temperature responsiveness of body size. These findings provide a novel step toward the molecular understanding of the temperature–size rule, which has puzzled biologists for decades.     

Mapping determinants of gene expression plasticity by genetical genomics in C. elegans

Recent genetical genomics studies have provided intimate views on gene regulatory networks. Gene expression variations between genetically different individuals have been mapped to the causal regulatory regions, termed expression quantitative trait loci. Whether the environment-induced plastic response of gene expression also shows heritable difference has not yet been studied. Here we show that differential expression induced by temperatures of 16 °C and 24 °C has a strong genetic component in Caenorhabditis elegans recombinant inbred strains derived from a cross between strains CB4856 (Hawaii) and N2 (Bristol). No less than 59% of 308 trans-acting genes showed a significant eQTL-by-environment interaction, here termed plasticity quantitative trait loci. In contrast, only 8% of an estimated 188 cis-acting genes showed such interaction. This indicates that heritable differences in plastic responses of gene expression are largely regulated in trans. This regulation is spread over many different regulators. However, for one group of trans-genes we found prominent evidence for a common master regulator: a transband of 66 coregulated genes appeared at 24 °C. Our results suggest widespread genetic variation of differential expression responses to environmental impacts and demonstrate the potential of genetical genomics for mapping the molecular determinants of phenotypic plasticity.     

Polypeptide analysis of genotypic variants of occluded Heliothis spp. baculoviruses

The structural polypeptides of 12 baculovirus isolates which included nuclear polyhedrosis viruses (NPVs) and granulosis viruses (GVs) obtained from four different species of the insect genus Heliothis collected in different geographical regions of the world were characterized using sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The matrix proteins were compared according to their molecular weights and peptide profiles produced after limited proteolysis. Examination of the matrix and virion polypeptide profiles revealed three major polypeptide phenotypes which corresponded to the three baculovirus morphological groups; singly embedded nuclear polyhedrosis viruses (SNPVs), multiply embedded nuclear polyhedrosis viruses (MNPVs), and granulosis viruses (GVs). Enveloped nucleocapsid polypeptide profiles of isolates within each NPV phenotype differed in only one polypeptide whereas the two GV isolates differed by as many as five polypeptides. Nucleocapsid polypeptide profiles of isolates within each of the NPV subgroups were identical while those profiles from the GV nucleocapsids differed slightly in molecular weight of one polypeptide.     

Remarkably Divergent Regions Punctuate the Genome Assembly of the Caenorhabditis elegans Hawaiian Strain CB4856

The Hawaiian strain (CB4856) of Caenorhabditis elegans is one of the most divergent from the canonical laboratory strain N2 and has been widely used in developmental, population, and evolutionary studies. To enhance the utility of the strain, we have generated a draft sequence of the CB4856 genome, exploiting a variety of resources and strategies. When compared against the N2 reference, the CB4856 genome has 327,050 single nucleotide variants (SNVs) and 79,529 insertion–deletion events that result in a total of 3.3 Mb of N2 sequence missing from CB4856 and 1.4 Mb of sequence present in CB4856 but not present in N2. As previously reported, the density of SNVs varies along the chromosomes, with the arms of chromosomes showing greater average variation than the centers. In addition, we find 61 regions totaling 2.8 Mb, distributed across all six chromosomes, which have a greatly elevated SNV density, ranging from 2 to 16% SNVs. A survey of other wild isolates show that the two alternative haplotypes for each region are widely distributed, suggesting they have been maintained by balancing selection over long evolutionary times. These divergent regions contain an abundance of genes from large rapidly evolving families encoding F-box, MATH, BATH, seven-transmembrane G-coupled receptors, and nuclear hormone receptors, suggesting that they provide selective advantages in natural environments. The draft sequence makes available a comprehensive catalog of sequence differences between the CB4856 and N2 strains that will facilitate the molecular dissection of their phenotypic differences. Our work also emphasizes the importance of going beyond simple alignment of reads to a reference genome when assessing differences between genomes.     

Geometric effects in gas vesicle buckling under ultrasound

Acoustic reporter genes based on gas vesicles (GVs) have enabled the use of ultrasound to noninvasively visualize cellular function in vivo. The specific detection of GV signals relative to background acoustic scattering in tissues is facilitated by nonlinear ultrasound imaging techniques taking advantage of the sonomechanical buckling of GVs. However, the effect of geometry on the buckling behavior of GVs under exposure to ultrasound has not been studied. To understand such geometric effects, we developed computational models of GVs of various lengths and diameters and used finite element simulations to predict their threshold buckling pressures and postbuckling deformations. We demonstrated that the GV diameter has an inverse cubic relation to the threshold buckling pressure, whereas length has no substantial effect. To complement these simulations, we experimentally probed the effect of geometry on the mechanical properties of GVs and the corresponding nonlinear ultrasound signals. The results of these experiments corroborate our computational predictions. This study provides fundamental insights into how geometry affects the sonomechanical properties of GVs, which, in turn, can inform further engineering of these nanostructures for high-contrast, nonlinear ultrasound imaging.     

Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

Over the last 20 years, studies of Caenorhabditis elegans natural diversity have demonstrated the power of quantitative genetic approaches to reveal the evolutionary, ecological, and genetic factors that shape traits. These studies complement the use of the laboratory-adapted strain N2 and enable additional discoveries not possible using only one genetic background. In this chapter, we describe how to perform quantitative genetic studies in Caenorhabditis, with an emphasis on C. elegans. These approaches use correlations between genotype and phenotype across populations of genetically diverse individuals to discover the genetic causes of phenotypic variation. We present methods that use linkage, near-isogenic lines, association, and bulk-segregant mapping, and we describe the advantages and disadvantages of each approach. The power of C. elegans quantitative genetic mapping is best shown in the ability to connect phenotypic differences to specific genes and variants. We will present methods to narrow genomic regions to candidate genes and then tests to identify the gene or variant involved in a quantitative trait. The same features that make C. elegans a preeminent experimental model animal contribute to its exceptional value as a tool to understand natural phenotypic variation.     

The antimicrobial susceptibility,biofilm formation and genotypic profiles of Staphylococcus haemolyticus from bloodstream infections

We analysed the antimicrobial susceptibility, biofilm formation and genotypic profiles of 27 isolates of Staphylococcus haemolyticus obtained from the blood of 19 patients admitted to a hospital in Rio de Janeiro, Brazil. Our analysis revealed a clinical significance of 36.8% and a multi-resistance rate of 92.6% among these isolates. All but one isolate carried the mecA gene. The staphylococcal cassette chromosome mec type I was the most prevalent mec element detected (67%). Nevertheless, the isolates showed clonal diversity based on pulsed-field gel electrophoresis analysis. The ability to form biofilms was detected in 66% of the isolates studied. Surprisingly, no icaAD genes were found among the biofilm-producing isolates.     

Microsporidia Are Natural Intracellular Parasites of the Nematode Caenorhabditis elegans

For decades the soil nematode Caenorhabditis elegans has been an important model system for biology, but little is known about its natural ecology. Recently, C. elegans has become the focus of studies of innate immunity and several pathogens have been shown to cause lethal intestinal infections in C. elegans. However none of these pathogens has been shown to invade nematode intestinal cells, and no pathogen has been isolated from wild-caught C. elegans. Here we describe an intracellular pathogen isolated from wild-caught C. elegans that we show is a new species of microsporidia. Microsporidia comprise a large class of eukaryotic intracellular parasites that are medically and agriculturally important, but poorly understood. We show that microsporidian infection of the C. elegans intestine proceeds through distinct stages and is transmitted horizontally. Disruption of a conserved cytoskeletal structure in the intestine called the terminal web correlates with the release of microsporidian spores from infected cells, and appears to be part of a novel mechanism by which intracellular pathogens exit from infected cells. Unlike in bacterial intestinal infections, the p38 MAPK and insulin/insulin-like growth factor (IGF) signaling pathways do not appear to play substantial roles in resistance to microsporidian infection in C. elegans. We found microsporidia in multiple wild-caught isolates of Caenorhabditis nematodes from diverse geographic locations. These results indicate that microsporidia are common parasites of C. elegans in the wild. In addition, the interaction between C. elegans and its natural microsporidian parasites provides a system in which to dissect intracellular intestinal infection in vivo and insight into the diversity of pathogenic mechanisms used by intracellular microbes.     

Shifts in the Distribution of Mass Densities Is a Signature of Caloric Restriction in Caenorhabditis elegans

Although the starvation response of the model multicellular organism Caenorhabditis elegans is a subject of much research, there is no convenient phenotypic readout of caloric restriction that can be applicable to large numbers of worms. This paper describes the distribution of mass densities of populations of C. elegans, from larval stages up to day one of adulthood, using isopycnic centrifugation, and finds that density is a convenient, if complex, phenotypic readout in C. elegans. The density of worms in synchronized populations of wildtype N2 C. elegans grown under standard solid-phase culture conditions was normally distributed, with distributions peaked sharply at a mean of 1.091 g/cm³ for L1, L2 and L3 larvae, 1.087 g/cm³ for L4 larvae, 1.081 g/cm³ for newly molted adults, and 1.074 g/cm³ at 24 hours of adulthood. The density of adult worms under starvation stress fell well outside this range, falling to a mean value of 1.054 g/cm³ after eight hours of starvation. This decrease in density correlated with the consumption of stored glycogen in the food-deprived worms. The density of the worms increased when deprived of food for longer durations, corresponding to a shift in the response of the worms: worms sacrifice their bodies by retaining larvae, which consume the adults from within. Density-based screens with the drug Ivermectin on worms cultured on single plates resulted in a clear bimodal (double-peaked) distribution of densities corresponding to drug exposed and non-exposed worms. Thus, measurements of changes in density could be used to conduct screens on the effects of drugs on several populations of worms cultured on single plates.     

Unisexual and Heterosexual Meiotic Reproduction Generate Aneuploidy and Phenotypic Diversity De Novo in the Yeast Cryptococcus neoformans

Aneuploidy is known to be deleterious and underlies several common human diseases, including cancer and genetic disorders such as trisomy 21 in Down''s syndrome. In contrast, aneuploidy can also be advantageous and in fungi confers antifungal drug resistance and enables rapid adaptive evolution. We report here that sexual reproduction generates phenotypic and genotypic diversity in the human pathogenic yeast Cryptococcus neoformans, which is globally distributed and commonly infects individuals with compromised immunity, such as HIV/AIDS patients, causing life-threatening meningoencephalitis. C. neoformans has a defined a-α opposite sexual cycle; however, >99% of isolates are of the α mating type. Interestingly, α cells can undergo α-α unisexual reproduction, even involving genotypically identical cells. A central question is: Why would cells mate with themselves given that sex is costly and typically serves to admix preexisting genetic diversity from genetically divergent parents? In this study, we demonstrate that α-α unisexual reproduction frequently generates phenotypic diversity, and the majority of these variant progeny are aneuploid. Aneuploidy is responsible for the observed phenotypic changes, as chromosome loss restoring euploidy results in a wild-type phenotype. Other genetic changes, including diploidization, chromosome length polymorphisms, SNPs, and indels, were also generated. Phenotypic/genotypic changes were not observed following asexual mitotic reproduction. Aneuploidy was also detected in progeny from a-α opposite-sex congenic mating; thus, both homothallic and heterothallic sexual reproduction can generate phenotypic diversity de novo. Our study suggests that the ability to undergo unisexual reproduction may be an evolutionary strategy for eukaryotic microbial pathogens, enabling de novo genotypic and phenotypic plasticity and facilitating rapid adaptation to novel environments.     

In Silico Molecular Comparisons of C. elegans and Mammalian Pharmacology Identify Distinct Targets That Regulate Feeding

Phenotypic screens can identify molecules that are at once penetrant and active on the integrated circuitry of a whole cell or organism. These advantages are offset by the need to identify the targets underlying the phenotypes. Additionally, logistical considerations limit screening for certain physiological and behavioral phenotypes to organisms such as zebrafish and C. elegans. This further raises the challenge of elucidating whether compound-target relationships found in model organisms are preserved in humans. To address these challenges we searched for compounds that affect feeding behavior in C. elegans and sought to identify their molecular mechanisms of action. Here, we applied predictive chemoinformatics to small molecules previously identified in a C. elegans phenotypic screen likely to be enriched for feeding regulatory compounds. Based on the predictions, 16 of these compounds were tested in vitro against 20 mammalian targets. Of these, nine were active, with affinities ranging from 9 nM to 10 µM. Four of these nine compounds were found to alter feeding. We then verified the in vitro findings in vivo through genetic knockdowns, the use of previously characterized compounds with high affinity for the four targets, and chemical genetic epistasis, which is the effect of combined chemical and genetic perturbations on a phenotype relative to that of each perturbation in isolation. Our findings reveal four previously unrecognized pathways that regulate feeding in C. elegans with strong parallels in mammals. Together, our study addresses three inherent challenges in phenotypic screening: the identification of the molecular targets from a phenotypic screen, the confirmation of the in vivo relevance of these targets, and the evolutionary conservation and relevance of these targets to their human orthologs.     

Heterologous expression of cyanobacterial gas vesicle proteins in Saccharomyces cerevisiae

Given the potential applications of gas vesicles (GVs) in multiple fields including antigen-displaying and imaging, heterologous reconstitution of synthetic GVs is an attractive and interesting study that has translational potential. Here, we attempted to express and assemble GV proteins (GVPs) into GVs using the model eukaryotic organism Saccharomyces cerevisiae. We first selected and expressed two core structural proteins, GvpA and GvpC from cyanobacteria Anabaena flos-aquae and Planktothrix rubescens, respectively. We then optimized the protein production conditions and validated GV assembly in the context of GV shapes. We found that when two copies of anaA were integrated into the genome, the chromosomal expression of AnaA resulted in GV production regardless of GvpC expression. Next, we co-expressed chaperone-RFP with the GFP-AnaA to aid the AnaA aggregation. The co-expression of individual chaperones (Hsp42, Sis1, Hsp104, and GvpN) with AnaA led to the formation of larger inclusions and enhanced the sequestration of AnaA into the perivacuolar site. To our knowledge, this represents the first study on reconstitution of GVs in S. cerevisiae. Our results could provide insights into optimizing conditions for heterologous protein production as well as the reconstitution of other synthetic microcompartments in yeast.