Weiguo Feng

Biomedical Discovery Acceleration, with Applications to Craniofacial Development

The profusion of high-throughput instruments and the explosion of new results in the scientific literature, particularly in molecular biomedicine, is both a blessing and a curse to the bench researcher. Even knowledgeable and experienced scientists can benefit from computational tools that help navigate this vast and rapidly evolving terrain. In this paper, we describe a novel computational approach to this challenge, a knowledge-based system that combines reading, reasoning, and reporting methods to facilitate analysis of experimental data. Reading methods extract information from external resources, either by parsing structured data or using biomedical language processing to extract information from unstructured data, and track knowledge provenance. Reasoning methods enrich the knowledge that results from reading by, for example, noting two genes that are annotated to the same ontology term or database entry. Reasoning is also used to combine all sources into a knowledge network that represents the integration of all sorts of relationships between a pair of genes, and to calculate a combined reliability score. Reporting methods combine the knowledge network with a congruent network constructed from experimental data and visualize the combined network in a tool that facilitates the knowledge-based analysis of that data. An implementation of this approach, called the Hanalyzer, is demonstrated on a large-scale gene expression array dataset relevant to craniofacial development. The use of the tool was critical in the creation of hypotheses regarding the roles of four genes never previously characterized as involved in craniofacial development; each of these hypotheses was validated by further experimental work.

Conditional Deletion of Activating Protein 2α (AP-2α) in the Developing Retina Demonstrates Non-Cell-Autonomous Roles for AP-2α in Optic Cup Development

ABSTRACT Activating protein 2α (AP-2α) is known to be expressed in the retina, and AP-2α-null mice exhibit defects in the developing optic cup, including patterning of the neural retina (NR) and a replacement of the dorsal retinal pigmented epithelium (RPE) with NR. In this study, we analyzed the temporal and spatial retinal expression patterns of AP-2α and created a conditional deletion of AP-2α in the developing retina. AP-2α exhibited a distinct expression pattern in the developing inner nuclear layer of the retina, and colocalization studies indicated that AP-2α was exclusively expressed in postmitotic amacrine cell populations. Targeted deletion of AP-2α in the developing retina did not result in observable retinal defects. Further examination of AP-2α-null mutants revealed that the severity of the RPE defect was variable and, although defects in retinal lamination occur at later embryonic stages, earlier stages showed normal lamination and expression of markers for amacrine and ganglion cells. Together, these data demonstrate that, whereas AP-2α alone does not play an intrinsic role in retinogenesis, it has non-cell-autonomous effects on optic cup development. Additional expression analyses showed that multiple AP-2 proteins are present in the developing retina, which will be important to future studies.

Spatial and Temporal Analysis of Gene Expression during Growth and Fusion of the Mouse Facial Prominences

Orofacial malformations resulting from genetic and/or environmental causes are frequent human birth defects yet their etiology is often unclear because of insufficient information concerning the molecular, cellular and morphogenetic processes responsible for normal facial development. We have, therefore, derived a comprehensive expression dataset for mouse orofacial development, interrogating three distinct regions – the mandibular, maxillary and frontonasal prominences. To capture the dynamic changes in the transcriptome during face formation, we sampled five time points between E10.5–E12.5, spanning the developmental period from establishment of the prominences to their fusion to form the mature facial platform. Seven independent biological replicates were used for each sample ensuring robustness and quality of the dataset. Here, we provide a general overview of the dataset, characterizing aspects of gene expression changes at both the spatial and temporal level. Considerable coordinate regulation occurs across the three prominences during this period of facial growth and morphogenesis, with a switch from expression of genes involved in cell proliferation to those associated with differentiation. An accompanying shift in the expression of polycomb and trithorax genes presumably maintains appropriate patterns of gene expression in precursor or differentiated cells, respectively. Superimposed on the many coordinated changes are prominence-specific differences in the expression of genes encoding transcription factors, extracellular matrix components, and signaling molecules. Thus, the elaboration of each prominence will be driven by particular combinations of transcription factors coupled with specific cell:cell and cell:matrix interactions. The dataset also reveals several prominence-specific genes not previously associated with orofacial development, a subset of which we externally validate. Several of these latter genes are components of bidirectional transcription units that likely share cis-acting sequences with well-characterized genes. Overall, our studies provide a valuable resource for probing orofacial development and a robust dataset for bioinformatic analysis of spatial and temporal gene expression changes during embryogenesis.

Vgll2a is required for neural crest cell survival during zebrafish craniofacial development

Invertebrate and vertebrate vestigial (vg) and vestigial-like (VGLL) genes are involved in embryonic patterning and cell fate determination. These genes encode cofactors that interact with members of the Scalloped/TEAD family of transcription factors and modulate their activity. We have previously shown that, in mice, Vgll2 is differentially expressed in the developing facial prominences. In this study, we show that the zebrafish ortholog vgll2a is expressed in the pharyngeal endoderm and ectoderm surrounding the neural crest derived mesenchyme of the pharyngeal arches. Moreover, both the FGF and retinoic acid (RA) signaling pathways, which are critical components of the hierarchy controlling craniofacial patterning, regulate this domain of vgll2a expression. Consistent with these observations, vgll2a is required within the pharyngeal endoderm for NCC survival and pharyngeal cartilage development. Specifically, knockdown of Vgll2a in zebrafish embryos using Morpholino injection results in increased cell death within the pharyngeal arches, aberrant endodermal pouch morphogenesis, and hypoplastic cranial cartilages. Overall, our data reveal a novel non-cell autonomous role for Vgll2a in development of the NCC-derived vertebrate craniofacial skeleton.

Tfap2a-dependent changes in mouse facial morphology result in clefting that can be ameliorated by a reduction in Fgf8 gene dosage

Failure of facial prominence fusion causes cleft lip and palate (CL/P), a common human birth defect. Several potential mechanisms can be envisioned that would result in CL/P, including failure of prominence growth and/or alignment as well as a failure of fusion of the juxtaposed epithelial seams. Here, using geometric morphometrics, we analyzed facial outgrowth and shape change over time in a novel mouse model exhibiting fully penetrant bilateral CL/P. This robust model is based upon mutations in Tfap2a, the gene encoding transcription factor AP-2α, which has been implicated in both syndromic and non-syndromic human CL/P. Our findings indicate that aberrant morphology and subsequent misalignment of the facial prominences underlies the inability of the mutant prominences to fuse. Exencephaly also occured in some of the Tfap2a mutants and we observed additional morphometric differences that indicate an influence of neural tube closure defects on facial shape. Molecular analysis of the CL/P model indicates that Fgf signaling is misregulated in the face, and that reducing Fgf8 gene dosage can attenuate the clefting pathology by generating compensatory changes. Furthermore, mutations in either Tfap2a or Fgf8 increase variance in facial shape, but the combination of these mutations restores variance to normal levels. The alterations in variance provide a potential mechanistic link between clefting and the evolution and diversity of facial morphology. Overall, our findings suggest that CL/P can result from small gene-expression changes that alter the shape of the facial prominences and uncouple their coordinated morphogenesis, which is necessary for normal fusion.

Identification and analysis of a conserved Tcfap2a intronic enhancer element required for expression in facial and limb bud mesenchyme.

ABSTRACT Tcfap2a, the gene encoding the mouse AP-2α transcription factor, is required for normal development of multiple structures during embryogenesis, including the face and limbs. Using comparative sequence analysis and transgenic-mouse experiments we have identified an intronic enhancer within this gene that directs expression to the face and limb mesenchyme. There are two conserved sequence blocks within this intron, and the larger of these directs tissue-specific activity and is found in all vertebrate Tcfap2a genes analyzed. To assess the role of the enhancer in regulating endogenous mouse Tcfap2a expression, we have deleted this cis-regulatory sequence from the genome. Loss of this element severely impairs Tcfap2a expression in the limb bud mesenchyme but generates only a modest reduction in the facial mesenchyme. The reduction in Tcfap2a transcription is accompanied by altered patterning of the forelimb, resulting in postaxial polydactyly. These results indicate that the major role for this enhancer resides within the limb bud, and it serves to maintain a level of Tcfap2a expression that limits the size of the hand plate and the associated number of digit primordia. The potential role of this cis-acting sequence in modeling the size and shape of the face and limbs during evolution is discussed.

Leveraging existing biological knowledge in the identification of candidate genes for facial dysmorphology.

BackgroundIn response to the frequently overwhelming output of high-throughput microarray experiments, we propose a methodology to facilitate interpretation of biological data in the context of existing knowledge. Through the probabilistic integration of explicit and implicit data sources a functional interaction network can be constructed. Each edge connecting two proteins is weighted by a confidence value capturing the strength and reliability of support for that interaction given the combined data sources. The resulting network is examined in conjunction with expression data to identify groups of genes with significant temporal or tissue specific patterns. In contrast to unstructured gene lists, these networks often represent coherent functional groupings.ResultsBy linking from shared functional categorizations to primary biological resources we apply this method to craniofacial microarray data, generating biologically testable hypotheses and identifying candidate genes for craniofacial development.ConclusionThe novel methodology presented here illustrates how the effective integration of pre-existing biological knowledge and high-throughput experimental data drives biological discovery and hypothesis generation.

Systems biology of facial development: contributions of ectoderm and mesenchyme

The rapid increase in gene-centric biological knowledge coupled with analytic approaches for genomewide data integration provides an opportunity to develop systems-level understanding of facial development. Experimental analyses have demonstrated the importance of signaling between the surface ectoderm and the underlying mesenchyme are coordinating facial patterning. However, current transcriptome data from the developing vertebrate face is dominated by the mesenchymal component, and the contributions of the ectoderm are not easily identified. We have generated transcriptome datasets from critical periods of mouse face formation that enable gene expression to be analyzed with respect to time, prominence, and tissue layer. Notably, by separating the ectoderm and mesenchyme we considerably improved the sensitivity compared to data obtained from whole prominences, with more genes detected over a wider dynamic range. From these data we generated a detailed description of ectoderm-specific developmental programs, including pan-ectodermal programs, prominence- specific programs and their temporal dynamics. The genes and pathways represented in these programs provide mechanistic insights into several aspects of ectodermal development. We also used these data to identify co-expression modules specific to facial development. We then used 14 co-expression modules enriched for genes involved in orofacial clefts to make specific mechanistic predictions about genes involved in tongue specification, in nasal process patterning and in jaw development. Our multidimensional gene expression dataset is a unique resource for systems analysis of the developing face; our co-expression modules are a resource for predicting functions of poorly annotated genes, or for predicting roles for genes that have yet to be studied in the context of facial development; and our analytic approaches provide a paradigm for analysis of other complex developmental programs.