Variability in gene expression underlies incomplete penetrance

The phenotypic differences between individual organisms can often be ascribed to underlying genetic and environmental variation. However, even genetically identical organisms in homogeneous environments vary, indicating that randomness in developmental processes such as gene expression may also generate diversity. To examine the consequences of gene expression variability in multicellular organisms, we studied intestinal specification in the nematode Caenorhabditis elegans in which wild-type cell fate is invariant and controlled by a small transcriptional network. Mutations in elements of this network can have indeterminate effects: some mutant embryos fail to develop intestinal cells, whereas others produce intestinal precursors. By counting transcripts of the genes in this network in individual embryos, we show that the expression of an otherwise redundant gene becomes highly variable in the mutants and that this variation is subjected to a threshold, producing an ON/OFF expression pattern of the master regulatory gene of intestinal differentiation. Our results demonstrate that mutations in developmental networks can expose otherwise buffered stochastic variability in gene expression, leading to pronounced phenotypic variation.National Institutes of Health (U.S.). Pioneer AwardMathematical Sciences Postdoctoral Research Fellowships (DMS-0603392)National Institutes of Health (U.S.). Ruth L. Kirschstein National Research Service Award (5F32GM080966

TbUNC119 and Its Binding Protein Complex Are Essential for Propagation, Motility, and Morphogenesis of Trypanosoma brucei Procyclic Form Cells

Flagellum-mediated motility of Trypanosoma brucei is considered to be essential for the parasite to complete stage development in the tsetse fly vector, while the mechanism by which flagellum-mediated motility is controlled are not fully understood. We thus compared T. brucei whole gene products (amino acid sequence) with Caenorhabditis elegans UNC (uncoordinated) proteins, in order to find uncharacterized motility-related T. brucei genes. Through in silico analysis, we found 88 gene products which were highly similar to C. elegans UNC proteins and categorized them as TbCEUN (T. brucei gene products which have high similarity to C. elegans UNC proteins). Approximately two thirds of the 88 TbCEUN gene products were kinesin-related molecules. A gene product highly similar to C. elegans UNC119 protein was designated as TbUNC119. RNAi-mediated depletion of TbUNC119 showed no apparent phenotype. However, knock-down analysis of both TbUNC119 and its binding protein (TbUNC119BP) which was found by yeast two-hybrid analysis showed characteristic phenotypes, including reduced motility, morphological change (extended cell shape), and cellular apoptosis. Based on the observed phenotypes, possible function of the TbUNC119 and TbUNC119BP is discussed

A Genome-Wide RNAi Screen for Factors Involved in Neuronal Specification in Caenorhabditis elegans

One of the central goals of developmental neurobiology is to describe and understand the multi-tiered molecular events that control the progression of a fertilized egg to a terminally differentiated neuron. In the nematode Caenorhabditis elegans, the progression from egg to terminally differentiated neuron has been visually traced by lineage analysis. For example, the two gustatory neurons ASEL and ASER, a bilaterally symmetric neuron pair that is functionally lateralized, are generated from a fertilized egg through an invariant sequence of 11 cellular cleavages that occur stereotypically along specific cleavage planes. Molecular events that occur along this developmental pathway are only superficially understood. We take here an unbiased, genome-wide approach to identify genes that may act at any stage to ensure the correct differentiation of ASEL. Screening a genome-wide RNAi library that knocks-down 18,179 genes (94% of the genome), we identified 245 genes that affect the development of the ASEL neuron, such that the neuron is either not generated, its fate is converted to that of another cell, or cells from other lineage branches now adopt ASEL fate. We analyze in detail two factors that we identify from this screen: (1) the proneural gene hlh-14, which we find to be bilaterally expressed in the ASEL/R lineages despite their asymmetric lineage origins and which we find is required to generate neurons from several lineage branches including the ASE neurons, and (2) the COMPASS histone methyltransferase complex, which we find to be a critical embryonic inducer of ASEL/R asymmetry, acting upstream of the previously identified miRNA lsy-6. Our study represents the first comprehensive, genome-wide analysis of a single neuronal cell fate decision. The results of this analysis provide a starting point for future studies that will eventually lead to a more complete understanding of how individual neuronal cell types are generated from a single-cell embryo

The Roles and Acting Mechanism of Caenorhabditis elegans DNase II Genes in Apoptotic DNA Degradation and Development

DNase II enzymes are acidic endonucleases that have been implicated in mediating apoptotic DNA degradation, a critical cell death execution event. C. elegans genome contains three DNase II homologues, NUC-1, CRN-6, and CRN-7, but their expression patterns, acting sites, and roles in apoptotic DNA degradation and development are unclear. We have conducted a comprehensive analysis of three C. elegans DNase II genes and found that nuc-1 plays a major role, crn-6 plays an auxiliary role, and crn-7 plays a negligible role in resolving 3′ OH DNA breaks generated in apoptotic cells. Promoter swapping experiments suggest that crn-6 but not crn-7 can partially substitute for nuc-1 in mediating apoptotic DNA degradation and both fail to replace nuc-1 in degrading bacterial DNA in intestine. Despite of their restricted and largely non-overlapping expression patterns, both CRN-6 and NUC-1 can mediate apoptotic DNA degradation in many cells, suggesting that they are likely secreted nucleases that are retaken up by other cells to exert DNA degradation functions. Removal or disruption of NUC-1 secretion signal eliminates NUC-1's ability to mediate DNA degradation across its expression border. Furthermore, blocking cell corpse engulfment does not affect apoptotic DNA degradation mediated by nuc-1, suggesting that NUC-1 acts in apoptotic cells rather than in phagocytes to resolve 3′ OH DNA breaks. Our study illustrates how multiple DNase II nucleases play differential roles in apoptotic DNA degradation and development and reveals an unexpected mode of DNase II action in mediating DNA degradation

The study of Priapulus caudatus reveals conserved molecular patterning underlying different gut morphogenesis in the Ecdysozoa

Background The digestive systems of animals can become highly specialized in response to their exploration and occupation of new ecological niches. Although studies on different animals have revealed commonalities in gut formation, the model systems Caenorhabditis elegans and Drosophila melanogaster, which belong to the invertebrate group Ecdysozoa, exhibit remarkable deviations in how their intestines develop. Their morphological and developmental idiosyncrasies have hindered reconstructions of ancestral gut characters for the Ecdysozoa, and limit comparisons with vertebrate models. In this respect, the phylogenetic position, and slow evolving morphological and molecular characters of marine priapulid worms advance them as a key group to decipher evolutionary events that occurred in the lineages leading to C. elegans and D. melanogaster. Results In the priapulid Priapulus caudatus, the gut consists of an ectodermal foregut and anus, and a mid region of at least partial endodermal origin. The inner gut develops into a 16-cell primordium devoid of visceral musculature, arranged in three mid tetrads and two posterior duplets. The mouth invaginates ventrally and shifts to a terminal anterior position as the ventral anterior ectoderm differentially proliferates. Contraction of the musculature occurs as the head region retracts into the trunk and resolves the definitive larval body plan. Despite obvious developmental differences with C. elegans and D. melanogaster, the expression in P. caudatus of the gut-related candidate genes NK2.1, foxQ2, FGF8/17/18, GATA456, HNF4, wnt1, and evx demonstrate three distinct evolutionarily conserved molecular profiles that correlate with morphologically identified sub-regions of the gut. Conclusions The comparative analysis of priapulid development suggests that a midgut formed by a single endodermal population of vegetal cells, a ventral mouth, and the blastoporal origin of the anus are ancestral features in the Ecdysozoa. Our molecular data on P. caudatus reveal a conserved ecdysozoan gut-patterning program and demonstrates that extreme morphological divergence has not been accompanied by major molecular innovations in transcriptional regulators during digestive system evolution in the Ecdysozoa. Our data help us understand the origins of the ecdysozoan body plan, including those of C. elegans and D. melanogaster, and this is critical for comparisons between these two prominent model systems and their vertebrate counterparts

Stochastic loss and gain of symmetric divisions in the C. elegans epidermis perturbs robustness of stem cell number

Biological systems are subject to inherent stochasticity. Nevertheless, development is remarkably robust, ensuring the consistency of key phenotypic traits such as correct cell numbers in a certain tissue. It is currently unclear which genes modulate phenotypic variability, what their relationship is to core components of developmental gene networks, and what is the developmental basis of variable phenotypes. Here, we start addressing these questions using the robust number of Caenorhabditis elegans epidermal stem cells, known as seam cells, as a readout. We employ genetics, cell lineage tracing, and single molecule imaging to show that mutations in lin-22, a Hes-related basic helix-loop-helix (bHLH) transcription factor, increase seam cell number variability. We show that the increase in phenotypic variability is due to stochastic conversion of normally symmetric cell divisions to asymmetric and vice versa during development, which affect the terminal seam cell number in opposing directions. We demonstrate that LIN-22 acts within the epidermal gene network to antagonise the Wnt signalling pathway. However, lin-22 mutants exhibit cell-to-cell variability in Wnt pathway activation, which correlates with and may drive phenotypic variability. Our study demonstrates the feasibility to study phenotypic trait variance in tractable model organisms using unbiased mutagenesis screens

Cell-type–specific eQTL of primary melanocytes facilitates identification of melanoma susceptibility genes

Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type–specific regulatory landscape of human melanocytes, which give rise to melanoma but account for <5% of typical human skin biopsies, we performed an eQTL analysis in primary melanocyte cultures from 106 newborn males. We identified 597,335 cis-eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans-eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis-regulation of IRF4. Melanocyte eQTLs are enriched in cis-regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes (ZFP90, HEBP1, MSC, CBWD1, and RP11-383H13.1) were associated with melanoma at genome-wide significant P-values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology