Julian C. Lui

Biological interpretation of genome-wide association studies using predicted gene functions

The main challenge for gaining biological insights from genetic associations is identifying which genes and pathways explain the associations. Here we present DEPICT, an integrative tool that employs predicted gene functions to systematically prioritize the most likely causal genes at associated loci, highlight enriched pathways and identify tissues/cell types where genes from associated loci are highly expressed. DEPICT is not limited to genes with established functions and prioritizes relevant gene sets for many phenotypes.

An Imprinted Gene Network that Controls Mammalian Somatic Growth is Down-Regulated During Postnatal Growth Deceleration in Multiple Organs

In mammals, somatic growth is rapid in early postnatal life but decelerates with age and eventually halts, thus determining the adult body size of the species. This growth deceleration, which reflects declining proliferation, occurs simultaneously in multiple organs yet appears not to be coordinated by a systemic mechanism. We, therefore, hypothesized that growth deceleration results from a growth-limiting genetic program that is common to multiple tissues. Here, we identified a set of 11 imprinted genes that show down-regulation of mRNA expression with age in multiple organs. For these genes, Igf2, H19, Plagl1, Mest, Peg3, Dlk1, Gtl2, Grb10, Ndn, Cdkn1c, and SLC38a4, the declines show a temporal pattern similar to the decline in growth rate. All 11 genes have been implicated in the control of cell proliferation or somatic growth. Thus, our findings suggest that the declining expression of these genes contributes to coordinate growth deceleration in multiple tissues. We next hypothesized that the coordinate decline in expression of these imprinted genes is caused by altered methylation and consequent silencing of the expressed allele. Contrary to this hypothesis, the methylation status of the promoter regions of Mest, Peg3, and Plagl1 did not change with age. Our findings suggest that a set of growth-regulating imprinted genes is expressed at high levels in multiple tissues in early postnatal life, contributing to rapid somatic growth, but that these genes are subsequently downregulated in multiple tissues simultaneously, contributing to coordinate growth deceleration and cessation, thus imposing a fundamental limit on adult body size.

Mechanisms Limiting Body Growth in Mammals

Recent studies have begun to provide insight into a long-standing mystery in biology-why body growth in animals is rapid in early life but then progressively slows, thus imposing a limit on adult body size. This growth deceleration in mammals is caused by potent suppression of cell proliferation in multiple tissues and is driven primarily by local, rather than systemic, mechanisms. Recent evidence suggests that this progressive decline in proliferation results from a genetic program that occurs in multiple organs and involves the down-regulation of a large set of growth-promoting genes. This program does not appear to be driven simply by time, but rather depends on growth itself, suggesting that the limit on adult body size is imposed by a negative feedback loop. Different organs appear to use different types of information to precisely target their adult size. For example, skeletal and cardiac muscle growth are negatively regulated by myostatin, the concentration of which depends on muscle mass itself. Liver growth appears to be modulated by bile acid flux, a parameter that reflects organ function. In pancreas, organ size appears to be limited by the initial number of progenitor cells, suggesting a mechanism based on cell-cycle counting. Further elucidation of the fundamental mechanisms suppressing juvenile growth is likely to yield important insights into the pathophysiology of childhood growth disorders and of the unrestrained growth of cancer. In addition, improved understanding of these growth-suppressing mechanisms may someday allow their therapeutic suspension in adult tissues to facilitate tissue regeneration.

Short stature, accelerated bone maturation, and early growth cessation due to heterozygous aggrecan mutations.

CONTEXT Many children with idiopathic short stature have a delayed bone age. Idiopathic short stature with advanced bone age is far less common. OBJECTIVE The aim was to identify underlying genetic causes of short stature with advanced bone age. SETTING AND DESIGN We used whole-exome sequencing to study three families with autosomal-dominant short stature, advanced bone age, and premature growth cessation. RESULTS Affected individuals presented with short stature [adult heights -2.3 to -4.2 standard deviation scores (SDS)] with histories of early growth cessation or childhood short stature (height SDS -1.9 to -3.5 SDS), advancement of bone age, and normal endocrine evaluations. Whole-exome sequencing identified novel heterozygous variants in ACAN, which encodes aggrecan, a proteoglycan in the extracellular matrix of growth plate and other cartilaginous tissues. The variants were present in all affected, but in no unaffected, family members. In Family 1, a novel frameshift mutation in exon 3 (c.272delA) was identified, which is predicted to cause early truncation of the aggrecan protein. In Family 2, a base-pair substitution was found in a highly conserved location within a splice donor site (c.2026+1G>A), which is also likely to alter the amino acid sequence of a large portion of the protein. In Family 3, a missense variant (c.7064T>C) in exon 14 affects a highly conserved residue (L2355P) and is strongly predicted to perturb protein function. CONCLUSIONS Our study demonstrates that heterozygous mutations in ACAN can cause a mild skeletal dysplasia, which presents clinically as short stature with advanced bone age. The accelerating effect on skeletal maturation has not previously been noted in the few prior reports of human ACAN mutations. Our findings thus expand the spectrum of ACAN defects and provide a new molecular genetic etiology for the unusual child who presents with short stature and accelerated skeletal maturation.

Synthesizing genome-wide association studies and expression microarray reveals novel genes that act in the human growth plate to modulate height

Previous meta-analysis of genome-wide association (GWA) studies has identified 180 loci that influence adult height. However, each GWA locus typically comprises a set of contiguous genes, only one of which presumably modulates height. We reasoned that many of the causative genes within these loci influence height because they are expressed in and function in the growth plate, a cartilaginous structure that causes bone elongation and thus determines stature. Therefore, we used expression microarray studies of mouse and rat growth plate, human disease databases and a mouse knockout phenotype database to identify genes within the GWAS loci that are likely required for normal growth plate function. Each of these approaches identified significantly more genes within the GWA height loci than at random genomic locations (P < 0.0001 each), supporting the validity of the approach. The combined analysis strongly implicates 78 genes in growth plate function, including multiple genes that participate in PTHrP-IHH, BMP and CNP signaling, and many genes that have not previously been implicated in the growth plate. Thus, this analysis reveals a large number of novel genes that regulate human growth plate chondrogenesis and thereby contribute to the normal variations in human adult height. The analytic approach developed for this study may be applied to GWA studies for other common polygenic traits and diseases, thus providing a new general strategy to identify causative genes within GWA loci and to translate genetic associations into mechanistic biological insights.

Copy number variants in patients with short stature

Height is a highly heritable and classic polygenic trait. Recent genome-wide association studies (GWAS) have revealed that at least 180 genetic variants influence adult height. However, these variants explain only about 10% of the phenotypic variation in height. Genetic analysis of short individuals can lead to the discovery of novel rare gene defects with a large effect on growth. In an effort to identify novel genes associated with short stature, genome-wide analysis for copy number variants (CNVs), using single-nucleotide polymorphism arrays, in 162 patients (149 families) with short stature was performed. Segregation analysis was performed if possible, and genes in CNVs were compared with information from GWAS, gene expression in rodents’ growth plates and published information. CNVs were detected in 40 families. In six families, a known cause of short stature was found (SHOX deletion or duplication, IGF1R deletion), in two combined with a de novo potentially pathogenic CNV. Thirty-three families had one or more potentially pathogenic CNVs (n=40). In 24 of these families, segregation analysis could be performed, identifying three de novo CNVs and nine CNVs segregating with short stature. Four were located near loci associated with height in GWAS (ADAMTS17, TULP4, PRKG2/BMP3 and PAPPA). Besides six CNVs known to be causative for short stature, 40 CNVs with possible pathogenicity were identified. Segregation studies and bioinformatics analysis suggested various potential candidate genes.

Novel microcephalic primordial dwarfism disorder associated with variants in the centrosomal protein ninein.

CONTEXT Microcephalic primordial dwarfism (MPD) is a rare, severe form of human growth failure in which growth restriction is evident in utero and continues into postnatal life. Single causative gene defects have been identified in a number of patients with MPD, and all involve genes fundamental to cellular processes including centrosome functions. OBJECTIVE The objective of the study was to find the genetic etiology of a novel presentation of MPD. DESIGN The design of the study was whole-exome sequencing performed on two affected sisters in a single family. Molecular and functional studies of a candidate gene were performed using patient-derived primary fibroblasts and a zebrafish morpholino oligonucleotides knockdown model. PATIENTS Two sisters presented with a novel subtype of MPD, including severe intellectual disabilities. MAIN OUTCOME MEASURES NIN, encoding Ninein, a centrosomal protein critically involved in asymmetric cell division, was identified as a candidate gene, and functional impacts in fibroblasts and zebrafish were studied. RESULTS From 34,606 genomic variants, two very rare missense variants in NIN were identified. Both probands were compound heterozygotes. In the zebrafish, ninein knockdown led to specific and novel defects in the specification and morphogenesis of the anterior neuroectoderm, resulting in a deformity of the developing cranium with a small, squared skull highly reminiscent of the human phenotype. CONCLUSION We identified a novel clinical subtype of MPD in two sisters who have rare variants in NIN. We show, for the first time, that reduction of ninein function in the developing zebrafish leads to specific deficiencies of brain and skull development, offering a developmental basis for the myriad phenotypes in our patients.

Coordinated postnatal down-regulation of multiple growth-promoting genes: evidence for a genetic program limiting organ growth

Children grow, but adults do not. The cessation of growth in multiple organs is the end result of a progressive decline in cell proliferation beginning in early life. The mechanisms responsible for this growth deceleration are largely unknown. Using expression microarray and real‐time PCR, we identified a common program of gene expression in lung, kidney, and liver during growth deceleration in juvenile rats. Gene ontology analyses and siRNA‐mediated knockdown in vitro indicated that many of the down‐regulated genes are growth promoting. Down‐regulated genes in the program showed declining histone H3K4 trimethylation with age, implicating underlying epigenetic mechanisms. To investigate the physiological processes driving the genetic program, a tryptophan‐deficient diet was used to temporarily inhibit juvenile growth in newborn rats for 4 wk. Afterward, microarray analysis showed that the genetic program had been delayed, implying that it is driven by body growth itself rather than age. Taken together, the findings suggest that growth in early life induces progressive down‐regulation of a large set of proliferation‐stimulating genes, causing organ growth to slow and eventually cease.—Lui, J. C., Forcinito, P., Chang, M., Chen, W., Barnes, K. M., Baron, J. Coordinated postnatal down‐regulation of multiple growth‐promoting genes: evidence for a genetic program limiting organ growth. FASEB J. 24, 3083–3092 (2010). www.fasebj.org

Organization of the Indian hedgehog – parathyroid hormone-related protein system in the postnatal growth plate

In embryonic growth cartilage, Indian hedgehog (Ihh) and parathyroid hormone-related protein (PTHrP) participate in a negative feedback loop that regulates chondrocyte differentiation. Postnatally, this region undergoes major structural and functional changes. To explore the organization of the Ihh–PTHrP system in postnatal growth plate, we microdissected growth plates of 7-day-old rats into their constituent zones and assessed expression of genes participating in the h–PTHrP feedback loop. Ihh, Patched 1, Smoothened, Gli1, Gli2, Gli3, and Pthr1 were expressed in regions analogous to the expression domains in embryonic growth cartilage. However, PTHrP was expressed in resting zone cartilage, a site that differs from the embryonic source, the periarticular cells. We then used mice in which lacZ has replaced coding sequences of Gli1 and thus serves as a marker for active hedgehog signaling. At 1, 4, 8, and 12 weeks of age, lacZ expression was detected in a pattern analogous to that of embryonic cartilage. The findings support the hypothesis that the embryonic Ihh–PTHrP feedback loop is maintained in the postnatal growth plate except that the source of PTHrP has shifted to a more proximal location in the resting zone.

Growth plate senescence and catch-up growth.

Longitudinal bone growth is rapid in prenatal and early postnatal life, but then slows with age and eventually ceases. This growth deceleration is caused primarily by a decrease in chondrocyte proliferation, and is associated with other structural, functional, and molecular changes collectively termed growth plate senescence. Current evidence suggests that growth plate senescence occurs because the progenitor chondrocytes in the resting zone have a limited replicative capacity which is gradually exhausted with increasing cell division. In addition, recent experimental findings from laboratory and clinical studies suggest that growth plate senescence explains the phenomenon of catch-up growth. Growth-inhibiting conditions such as glucocorticoid excess and hypothyroidism delay the program of growth plate senescence. Consequently, growth plates are less senescent after these conditions resolve and therefore grow more rapidly than is normal for age, resulting in catch-up growth.