The developmental basis for scaling of mammalian tooth size.

When evolution leads to differences in body size, organs generally scale along. A well-known example of the tight relationship between organ and body size is the scaling of mammalian molar teeth. To investigate how teeth scale during development and evolution, we compared molar development from initiation through final size in the mouse and the rat. Whereas the linear dimensions of the rat molars are twice that of the mouse molars, their shapes are largely the same. Here, we focus on the first lower molars that are considered the most reliable dental proxy for size-related patterns due to their low within-species variability. We found that scaling of the molars starts early, and that the rat molar is patterned equally as fast but in a larger size than the mouse molar. Using transcriptomics, we discovered that a known regulator of body size, insulin-like growth factor 1 (Igf1), is more highly expressed in the rat molars compared to the mouse molars. Ex vivo and in vivo mouse models demonstrated that modulation of the IGF pathway reproduces several aspects of the observed scaling process. Furthermore, analysis of IGF1-treated mouse molars and computational modeling indicate that IGF signaling scales teeth by simultaneously enhancing growth and by inhibiting the cusp-patterning program, thereby providing a relatively simple mechanism for scaling teeth during development and evolution. Finally, comparative data from shrews to elephants suggest that this scaling mechanism regulates the minimum tooth size possible, as well as the patterning potential of large teeth.

Chromosomal neighbourhoods allow identification of organ specific changes in gene expression.

Although most genes share their chromosomal neighbourhood with other genes, distribution of genes has not been explored in the context of individual organ development; the common focus of developmental biology studies. Because developmental processes are often associated with initially subtle changes in gene expression, here we explored whether neighbouring genes are informative in the identification of differentially expressed genes. First, we quantified the chromosomal neighbourhood patterns of genes having related functional roles in the mammalian genome. Although the majority of protein coding genes have at least five neighbours within 1 Mb window around each gene, very few of these neighbours regulate development of the same organ. Analyses of transcriptomes of developing mouse molar teeth revealed that whereas expression of genes regulating tooth development changes, their neighbouring genes show no marked changes, irrespective of their level of expression. Finally, we test whether inclusion of gene neighbourhood in the analyses of differential expression could provide additional benefits. For the analyses, we developed an algorithm, called DELocal that identifies differentially expressed genes by comparing their expression changes to changes in adjacent genes in their chromosomal regions. Our results show that DELocal removes detection bias towards large changes in expression, thereby allowing identification of even subtle changes in development. Future studies, including the detection of differential expression, may benefit from, and further characterize the significance of gene-gene neighbour relationships.

System-level analyses of keystone genes required for mammalian tooth development.

When a null mutation of a gene causes a complete developmental arrest, the gene is typically considered essential for life. Yet, in most cases, null mutations have more subtle effects on the phenotype. Here we used the phenotypic severity of mutations as a tool to examine system-level dynamics of gene expression. We classify genes required for the normal development of the mouse molar into different categories that range from essential to subtle modification of the phenotype. Collectively, we call these the developmental keystone genes. Transcriptome profiling using microarray and RNAseq analyses of patterning stage mouse molars show highly elevated expression levels for genes essential for the progression of tooth development, a result reminiscent of essential genes in single-cell organisms. Elevated expression levels of progression genes were also detected in developing rat molars, suggesting evolutionary conservation of this system-level dynamics. Single-cell RNAseq analyses of developing mouse molars reveal that even though the size of the expression domain, measured in the number of cells, is the main driver of organ-level expression, progression genes show high cell-level transcript abundances. Progression genes are also upregulated within their pathways, which themselves are highly expressed. In contrast, a high proportion of the genes required for normal tooth patterning are secreted ligands that are expressed in fewer cells than their receptors and intracellular components. Overall, even though expression patterns of individual genes can be highly different, conserved system-level principles of gene expression can be detected using phenotypically defined gene categories.

Vole genomics links determinate and indeterminate growth of teeth.

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Identification of antibodies against HAI-1 and integrin alpha6beta4 as immunohistochemical markers of human villous cytotrophoblast.

Syncytiotrophoblast and invasive extravillous trophoblast arise from a common stem cell, namely villous cytotrophoblast, but have very different characteristics. The study of the differentiation process relies on the availability of suitable markers for these different cell types of developing placenta. In this work, we have produced monoclonal antibodies that are specific to human villous cytotrophoblast. Monoclonal antibody (MAb) MG2 was specific to villous cytotrophoblast across gestation, and recognizes hepatocyte growth factor activator inhibitor type 1. MAb MD10 stained villous cytotrophoblast across gestation and also some endothelial cells, particularly in the second or third trimester. MAb MD10 recognizes human integrin alpha6beta4. As a test for specificity, the novel MAbs were also used for staining of frozen tissue from human colon carcinoma. The results show that the two antibodies can be used as tools to study human villous cytotrophoblasts and also human tumors. The MG2 antibody seems most specific and promising for the study of various aspects of human villous cytotrophoblast.

Mice lacking a Myc enhancer that includes human SNP rs6983267 are resistant to intestinal tumors.

Multiple cancer-associated single-nucleotide polymorphisms (SNPs) have been mapped to conserved sequences within a 500-kilobase region upstream of the MYC oncogene on human chromosome 8q24. These SNPs may affect cancer development through altered regulation of MYC expression, but this hypothesis has been difficult to confirm. We generated mice deficient in Myc-335, a putative MYC regulatory element that contains rs6983267, a SNP accounting for more human cancer-related morbidity than any other genetic variant or mutation. In Myc-335 null mice, Myc transcripts were expressed in the intestinal crypts in a pattern similar to that in wild-type mice but at modestly reduced levels. The mutant mice displayed no overt phenotype but were markedly resistant to intestinal tumorigenesis induced by the APCmin mutation. These results establish that a cancer-associated SNP identified in human genome-wide association studies has a functional effect in vivo.

The common colorectal cancer predisposition SNP rs6983267 at chromosome 8q24 confers potential to enhanced Wnt signaling.

Homozygosity for the G allele of rs6983267 at 8q24 increases colorectal cancer (CRC) risk approximately 1.5 fold. We report here that the risk allele G shows copy number increase during CRC development. Our computer algorithm, Enhancer Element Locator (EEL), identified an enhancer element that contains rs6983267. The element drove expression of a reporter gene in a pattern that is consistent with regulation by the key CRC pathway Wnt. rs6983267 affects a binding site for the Wnt-regulated transcription factor TCF4, with the risk allele G showing stronger binding in vitro and in vivo. Genome-wide ChIP assay revealed the element as the strongest TCF4 binding site within 1 Mb of MYC. An unambiguous correlation between rs6983267 genotype and MYC expression was not detected, and additional work is required to scrutinize all possible targets of the enhancer. Our work provides evidence that the common CRC predisposition associated with 8q24 arises from enhanced responsiveness to Wnt signaling.

High-throughput assay for determining specificity and affinity of protein-DNA binding interactions.

Limited information exists for the binding specificities of many important transcription factors. To address this, we have previously developed a microwell-based assay for directly measuring the affinity of DNA-protein binding interactions. We describe here the detailed protocol for determining sequence specificities of DNA-binding proteins using this assay. The described method is rapid; after preparation of the reagents, the assay can be run in a single day, and its throughput can be increased further by automation. The method is quantitative but requires prior knowledge of one high-affinity binding site for the protein of interest. The protocol can be adapted for determining the effect of protein modifications and protein-protein interactions on DNA-binding specificity, and for engineering proteins with new DNA-binding specificities. In addition, the method is suitable for high-throughput screening to identify proteins or small molecules that modulate protein-DNA binding interactions.

Genome-wide prediction of mammalian enhancers based on analysis of transcription-factor binding affinity.

Understanding the regulation of human gene expression requires knowledge of the "second genetic code," which consists of the binding specificities of transcription factors (TFs) and the combinatorial code by which TF binding sites are assembled to form tissue-specific enhancer elements. Using a novel high-throughput method, we determined the DNA binding specificities of GLIs 1-3, Tcf4, and c-Ets1, which mediate transcriptional responses to the Hedgehog (Hh), Wnt, and Ras/MAPK signaling pathways. To identify mammalian enhancer elements regulated by these pathways on a genomic scale, we developed a computational tool, enhancer element locator (EEL). We show that EEL can be used to identify Hh and Wnt target genes and to predict activated TFs based on changes in gene expression. Predictions validated in transgenic mouse embryos revealed the presence of multiple tissue-specific enhancers in mouse c-Myc and N-Myc genes, which has implications for organ-specific growth control and tumor-type specificity of oncogenes.