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

Med-type GATA factors and the evolution of mesendoderm specification in nematodes

In the nematode, C elegans, the divergent GATA-type transcription factors MED-1 and MED-2 are encoded by an unlinked, redundant pair of intronless genes. The med-1,2 genes are among the first to be activated in the embryo and are critical for the specification of the 7-cell stage MS (mesoderm) and E (endoderm) precursor cells. We have previously shown that the binding site recognized by MED-1 is a noncanonical RAGTATAC site that is not expected from the resemblance of its single C4-type zinc finger to those of other known GATA factors, which recognize the consensus HGATAR. To date, no MED-like zinc fingers have been described outside of C. elegans. In order to understand the evolution of these transcription factors, and the evolution of gene networks that specify early cell fates in Caenorhabditis, we have identified med sequence homologs in the related nematodes C. briggsae and C. remanei. While C. briggsae encodes two med-like genes similar to C. elegans, we find evidence for seven distinct med-like genes in C. remanei. Somewhat unexpectedly, the coding regions of all med genes appear to lack introns. We report that the med homologs have similar expression in their respective species. We further show that the C. briggsae homologs, and at least five of the seven C. remanei homologs, can fully complement the embryonic lethal phenotype of a C. elegans med-1,2(-) strain. We conclude that Med function and expression have been conserved over tens of millions of years of evolution, and that there may be a mechanism that selects against the acquisition of introns in these genes. (c) 2005 Elsevier Inc. All rights reserved

Gene expression profiling of gastrocnemius of minimuscle mice.

Few studies have investigated heterogeneity of selection response in replicate lines subjected to equivalent selection. We developed four replicate lines of mice based on high levels of voluntary wheel running (high runner or HR lines) while also maintaining four nonselected control lines. This led to the unexpected discovery of the HR minimuscle (HR(mini)) phenotype, recognized by a 50% reduction in hindlimb muscle mass, which became fixed in 1 of the four HR selected lines. Here, we report genome-wide expression profiling describing transcriptome differences between HR(normal) and HR(mini) medial gastrocnemius. Consistent with the known reduction of type IIB fibers in HR(mini), Myh4 gene expression was −8.82-fold less (P = 0.0001) in HR(mini), which was closely associated with differences in the “calcium signaling” canonical pathway, including structural genes (e.g., Mef2c, twofold greater in HR(mini), P = 0.0003) and myogenic factors (e.g., Myog, 3.8-fold greater in HR(mini), P = 0.0026) associated with slow-type myofibers. The gene that determines the HR(mini) phenotype is known to reside in a 2.6335-Mb interval on mouse chromosome 11 and 7 genes (Myh10, Chrnb1, Acadvl, Senp3, Gabarap, Eif5a, and Clec10a) from this region were differentially expressed. Verification by real-time PCR confirmed 1.5-fold greater (P < 0.05) expression of very long chain acyl-CoA dehydrogenase (Acadvl) in HR(mini). Ten other genes associated with fatty acid metabolism were also upregulated in HR(mini), suggesting differences in the ability to metabolize fatty acids in HR(normal) and HR(mini) muscles. This work provides a resource for understanding differences in muscle phenotypes in populations exhibiting high running capacity