Regulatory variation at glypican-3 underlies a major growth QTL in mice

The impact of species-wide gene expression variation on Caenorhabditis elegans complex traits

Phenotypic variation in organism-level traits has been studied in Caenorhabditis elegans wild strains, but the impacts of differences in gene expression and the underlying regulatory mechanisms are largely unknown. Here, we use natural variation in gene expression to connect genetic variants to differences in organismal-level traits, including drug and toxicant responses. We perform transcriptomic analyses on 207 genetically distinct C. elegans wild strains to study natural regulatory variation of gene expression. Using this massive dataset, we perform genome-wide association mappings to investigate the genetic basis underlying gene expression variation and reveal complex genetic architectures. We find a large collection of hotspots enriched for expression quantitative trait loci across the genome. We further use mediation analysis to understand how gene expression variation could underlie organism-level phenotypic variation for a variety of complex traits. These results reveal the natural diversity in gene expression and possible regulatory mechanisms in this keystone model organism, highlighting the promise of using gene expression variation to understand how phenotypic diversity is generated.

Mutational sources of trans-regulatory variation affecting gene expression in Saccharomyces cerevisiae

Heritable variation in a gene’s expression arises from mutations impacting cis- and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene’s expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.

Pleiotropy and epistasis within and between signaling pathways defines the genetic architecture of fungal virulence

Cryptococcal disease is estimated to affect nearly a quarter of a million people annually. Environmental isolates of Cryptococcus deneoformans, which make up 15 to 30% of clinical infections in temperate climates such as Europe, vary in their pathogenicity, ranging from benign to hyper-virulent. Key traits that contribute to virulence, such as the production of the pigment melanin, an extracellular polysaccharide capsule, and the ability to grow at human body temperature have been identified, yet little is known about the genetic basis of variation in such traits. Here we investigate the genetic basis of melanization, capsule size, thermal tolerance, oxidative stress resistance, and antifungal drug sensitivity using quantitative trait locus (QTL) mapping in progeny derived from a cross between two divergent C. deneoformans strains. Using a "function-valued" QTL analysis framework that exploits both time-series information and growth differences across multiple environments, we identified QTL for each of these virulence traits and drug susceptibility. For three QTL we identified the underlying genes and nucleotide differences that govern variation in virulence traits. One of these genes, RIC8, which encodes a regulator of cAMP-PKA signaling, contributes to variation in four virulence traits: melanization, capsule size, thermal tolerance, and resistance to oxidative stress. Two major effect QTL for amphotericin B resistance map to the genes SSK1 and SSK2, which encode key components of the HOG pathway, a fungal-specific signal transduction network that orchestrates cellular responses to osmotic and other stresses. We also discovered complex epistatic interactions within and between genes in the HOG and cAMP-PKA pathways that regulate antifungal drug resistance and resistance to oxidative stress. Our findings advance the understanding of virulence traits among diverse lineages of Cryptococcus, and highlight the role of genetic variation in key stress-responsive signaling pathways as a major contributor to phenotypic variation.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Giant Island Mice Exhibit Widespread Gene Expression Changes in Key Metabolic Organs

Island populations repeatedly evolve extreme body sizes, but the genomic basis of this pattern remains largely unknown. To understand how organisms on islands evolve gigantism, we compared genome-wide patterns of gene expression in Gough Island mice, the largest wild house mice in the world, and mainland mice from the WSB/EiJ wild-derived inbred strain. We used RNA-seq to quantify differential gene expression in three key metabolic organs: gonadal adipose depot, hypothalamus, and liver. Between 4,000 and 8,800 genes were significantly differentially expressed across the evaluated organs, representing between 20% and 50% of detected transcripts, with 20% or more of differentially expressed transcripts in each organ exhibiting expression fold changes of at least 2×. A minimum of 73 candidate genes for extreme size evolution, including Irs1 and Lrp1, were identified by considering differential expression jointly with other data sets: 1) genomic positions of published quantitative trait loci for body weight and growth rate, 2) whole-genome sequencing of 16 wild-caught Gough Island mice that revealed fixed single-nucleotide differences between the strains, and 3) publicly available tissue-specific regulatory elements. Additionally, patterns of differential expression across three time points in the liver revealed that Arid5b potentially regulates hundreds of genes. Functional enrichment analyses pointed to cell cycling, mitochondrial function, signaling pathways, inflammatory response, and nutrient metabolism as potential causes of weight accumulation in Gough Island mice. Collectively, our results indicate that extensive gene regulatory evolution in metabolic organs accompanied the rapid evolution of gigantism during the short time house mice have inhabited Gough Island.

Genetical genomics of growth in a chicken model

BACKGROUND

The genetics underlying body mass and growth are key to understanding a wide range of topics in biology, both evolutionary and developmental. Body mass and growth traits are affected by many genetic variants of small effect. This complicates genetic mapping of growth and body mass. Experimental intercrosses between individuals from divergent populations allows us to map naturally occurring genetic variants for selected traits, such as body mass by linkage mapping. By simultaneously measuring traits and intermediary molecular phenotypes, such as gene expression, one can use integrative genomics to search for potential causative genes.

RESULTS

In this study, we use linkage mapping approach to map growth traits (N = 471) and liver gene expression (N = 130) in an advanced intercross of wild Red Junglefowl and domestic White Leghorn layer chickens. We find 16 loci for growth traits, and 1463 loci for liver gene expression, as measured by microarrays. Of these, the genes TRAK1, OSBPL8, YEATS4, CEP55, and PIP4K2B are identified as strong candidates for growth loci in the chicken. We also show a high degree of sex-specific gene-regulation, with almost every gene expression locus exhibiting sex-interactions. Finally, several trans-regulatory hotspots were found, one of which coincides with a major growth locus.

CONCLUSIONS

These findings not only serve to identify several strong candidates affecting growth, but also show how sex-specificity and local gene-regulation affect growth regulation in the chicken.

Genetics of metabolic syndrome: potential clues from wild-derived inbred mouse strains

The metabolic syndrome (MetS) is a complex constellation of metabolic abnormalities including obesity, abnormal glucose metabolism, dyslipidemia, and elevated blood pressure that together substantially increase risk for cardiovascular disease and Type 2 diabetes. Both genetic and environmental factors contribute to the development of MetS, but this process is still far from understood. Human studies have revealed only part of the underlying basis. Studies in mice offer many strengths that can complement human studies to help elucidate the etiology and pathophysiology of MetS. Here we review the ways mice can contribute to MetS research. In particular, we focus on the information that can be obtained from studies of the inbred strains, with specific focus on the phenotypes of the wild-derived inbred strains. These are newly derived inbred strains that were created from wild-caught mice. They contain substantial genetic variation that is not present in the classical inbred strains, have phenotypes of relevance for MetS, and various mouse strain resources have been created to facilitate the mining of this new genetic variation. Thus studies using wild-derived inbred strains hold great promise for increasing our understanding of MetS.

Genetics of Rapid and Extreme Size Evolution in Island Mice

Organisms on islands provide a revealing window into the process of adaptation. Populations that colonize islands often evolve substantial differences in body size from their mainland relatives. Although the ecological drivers of this phenomenon have received considerable attention, its genetic basis remains poorly understood. We use house mice (subspecies: Mus musculus domesticus) from remote Gough Island to provide a genetic portrait of rapid and extreme size evolution. In just a few hundred generations, Gough Island mice evolved the largest body size among wild house mice from around the world. Through comparisons with a smaller-bodied wild-derived strain from the same subspecies (WSB/EiJ), we demonstrate that Gough Island mice achieve their exceptional body weight primarily by growing faster during the 6 weeks after birth. We use genetic mapping in large F(2) intercrosses between Gough Island mice and WSB/EiJ to identify 19 quantitative trait loci (QTL) responsible for the evolution of 16-week weight trajectories: 8 QTL for body weight and 11 QTL for growth rate. QTL exhibit modest effects that are mostly additive. We conclude that body size evolution on islands can be genetically complex, even when substantial size changes occur rapidly. In comparisons to published studies of laboratory strains of mice that were artificially selected for divergent body sizes, we discover that the overall genetic profile of size evolution in nature and in the laboratory is similar, but many contributing loci are distinct. Our results underscore the power of genetically characterizing the entire growth trajectory in wild populations and lay the foundation necessary for identifying the mutations responsible for extreme body size evolution in nature.

Quantitative trait loci for energy balance traits in an advanced intercross line derived from mice divergently selected for heat loss

Obesity in human populations, currently a serious health concern, is considered to be the consequence of an energy imbalance in which more energy in calories is consumed than is expended. We used interval mapping techniques to investigate the genetic basis of a number of energy balance traits in an F₁₁ advanced intercross population of mice created from an original intercross of lines selected for increased and decreased heat loss. We uncovered a total of 137 quantitative trait loci (QTLs) for these traits at 41 unique sites on 18 of the 20 chromosomes in the mouse genome, with X-linked QTLs being most prevalent. Two QTLs were found for the selection target of heat loss, one on distal chromosome 1 and another on proximal chromosome 2. The number of QTLs affecting the various traits generally was consistent with previous estimates of heritabilities in the same population, with the most found for two bone mineral traits and the least for feed intake and several body composition traits. QTLs were generally additive in their effects, and some, especially those affecting the body weight traits, were sex-specific. Pleiotropy was extensive within trait groups (body weights, adiposity and organ weight traits, bone traits) and especially between body composition traits adjusted and not adjusted for body weight at sacrifice. Nine QTLs were found for one or more of the adiposity traits, five of which appeared to be unique. The confidence intervals among all QTLs averaged 13.3 Mb, much smaller than usually observed in an F₂ cross, and in some cases this allowed us to make reasonable inferences about candidate genes underlying these QTLs. This study combined QTL mapping with genetic parameter analysis in a large segregating population, and has advanced our understanding of the genetic architecture of complex traits related to obesity.

Preparation and characterization of monoclonal antibody against glypican-3

Glypican-3 (GPC3) has been reported as a novel serum and histochemical marker for hepatocellular carcinoma (HCC) by several groups. As an oncofetal protein, it is expressed abundantly in the fetal liver, inactive in the normal adult liver, and frequently reactivated in HCC. Immunology reagents are urgently needed to proceed with mechanism-related research, clinical validation, and application. In this report, monoclonal antibodies (MAbs) against GPC3 were made from hyperimmune BALB/c mice by injecting 100 μg of purified antigen intraperitoneally. Hybridomas were screened by indirect enzyme-linked immunosorbent assay (ELISA) using purified protein. Finally 13 mouse hybridomas producing MAbs to GPC3 were established. The MAbs obtained were fully characterized using Western blot analysis, immunofluorescence, and immunohistochemistry. The results showed that these antibodies could be used for preliminary application of the next step mechanism-related research and GPC3 expression level analysis.

Overexpression of GPC3 inhibits hepatocellular carcinoma cell proliferation and invasion through induction of apoptosis

Glypican‑3 (GPC3) is a membrane heparan sulfate proteoglycan involved in cell proliferation, differentiation, adhesion, migration and the development of the majority of mesodermal tissues and organs. GPC3 has been found to be important for the occurrence and development of hepatocellular carcinoma (HCC). Therefore, it may be suitable for use as a novel molecular marker for the diagnosis of primary liver cancer. In the present study, the role of GPC3 in the occurrence and development of HCC was determined. GPC3 recombinant vector was transfected into two HCC cell lines, Huh7 and SK‑HEP‑1, to upregulate the expression of GPC3 and examine changes in the biological behavior of the cells. Results indicate that overexpression of GPC3 in Huh7 and SK‑HEP‑1 cells effectively inhibited cell proliferation and cell invasion through induction of apoptosis. However, cotreatment of the cells with insulin‑like growth factor 2 (IGF2) and fibroblast growth factor 2 (FGF2) was found by Annexin V‑PI flow cytometric analysis to significantly inhibit the apoptotic cell death induced by GPC3 overexpression. These observations indicate that GPC3 may act as a negative regulator of IGF2 and FGF2 pathways. Taken together, these results demonstrate that overexpression of GPC3 inhibits the occurrence and development of HCC.

The genetic architecture of growth rate in juvenile Takifugu species

Closely related species have often evolved dramatic differences in body size. Takifugu rubripes (fugu) is a large marine pufferfish whose genome has been sequenced, whereas T. niphobles is the smallest species among Takifugu. We show that, unsurprisingly, the juvenile growth rate of T. rubripes is higher than that of T. niphobles in a laboratory setting. We produced F(2) progenies of their F(1) hybrids and found one quantitative trait locus (QTL) significantly associated with variation in juvenile body size. This QTL region (3.5 Mb) contains no known genes directly related to growth phenotype (such as IGFs) except Fgf21, which inhibits growth hormone signaling in mouse. The QTL in Takifugu spp. is distinct from the region previously known to control body size variations in stickleback or tilapia. Our results suggest that in the fish tested herein, genomic regions underlying body size evolution might have different genetic origins. They also suggest that many diverse traits in Takifugu spp. are amenable to genetic mapping.

LRP1 mediates Hedgehog-induced endocytosis of the GPC3-Hedgehog complex

Glypican-3 (GPC3) is a heparan sulfate (HS) proteoglycan that is bound to the cell membrane through a glycosylphosphatidylinositol link. This glypican regulates embryonic growth by inhibiting the hedgehog (Hh) signaling pathway. GPC3 binds Hh and competes with Patched (Ptc), the Hh receptor, for Hh binding. The interaction of Hh with GPC3 triggers the endocytosis and degradation of the GPC3-Hh complex with the consequent reduction of Hh available for binding to Ptc. Currently, the molecular mechanisms by which the GPC3-Hh complex is internalized remains unknown. Here we show that the low-density-lipoprotein receptor-related protein-1 (LRP1) mediates the Hh-induced endocytosis of the GPC3-Hh complex, and that this endocytosis is necessary for the Hh-inhibitory activity of GPC3. Furthermore, we demonstrate that GPC3 binds through its HS chains to LRP1, and that this interaction causes the removal of GPC3 from the lipid rafts domains.

Mapping genetic loci that interact with myostatin to affect growth traits

Myostatin, or GDF8, is an inhibitor of skeletal muscle growth. A non-functional myostatin mutation leads to a double muscling phenotype in some species, for example, mice, cattle and humans. Previous studies have indicated that there are loci in the genome that interact with myostatin to control backfat depth and other complex traits. We now report a quantitative trait loci (QTL) mapping study designed to identify loci that interact with myostatin to impact growth traits in mice. Body weight and average daily gain traits were collected on F2 progeny derived from a myostatin-null C57BL/6 strain by M16i cross. In all, 44 main effect QTL were detected above a 5% genome-wide significance threshold when an interval mapping method was used. An additional 37 QTL were identified to significantly interact with myostatin, sex or reciprocal cross. A total of 12 of these QTL interacted with myostatin genotype. These results provide a foundation for the further fine mapping of genome regions that harbor loci that interact with myostatin.

Intersubspecific subcongenic mouse strain analysis reveals closely linked QTLs with opposite effects on body weight

A previous genome-wide QTL study revealed many QTLs affecting postnatal body weight and growth in an intersubspecific backcross mouse population between the C57BL/6J (B6) strain and wild Mus musculus castaneus mice captured in the Philippines. Subsequently, several closely linked QTLs for body composition traits were revealed in an F(2) intercross population between B6 and B6.Cg-Pbwg1, a congenic strain on the B6 genetic background carrying the growth QTL Pbwg1 on proximal chromosome 2. However, no QTL affecting body weight has been duplicated in the F(2) population, except for mapping an overdominant QTL that causes heterosis of body weight. In this study, we developed 17 intersubspecific subcongenic strains with overlapping and nonoverlapping castaneus regions from the B6.Cg-Pbwg1 congenic strain in order to search for and genetically dissect QTLs affecting body weight into distinct closely linked loci. Phenotypic comparisons of several developed subcongenic strains with the B6 strain revealed that two closely linked but distinct QTLs that regulate body weight, named Pbwg1.11 and Pbwg1.12, are located on an 8.9-Mb region between D2Mit270 and D2Mit472 and on the next 3.6-Mb region between D2Mit205 and D2Mit182, respectively. Further analyses using F(2) segregating populations obtained from intercrosses between B6 and each of the two selected subcongenic strains confirmed the presence of these two body weight QTLs. Pbwg1.11 had an additive effect on body weight at 6, 10, and 13 weeks of age, and its castaneus allele decreased it. In contrast, the castaneus allele at Pbwg1.12 acted in a dominant fashion and surprisingly increased body weight at 6, 10, and 13 weeks of age despite the body weight of wild castaneus mice being 60% of that of B6 mice. These findings illustrate the complex genetic nature of body weight regulation and support the importance of subcongenic mouse analysis to dissect closely linked loci.

The genetic structure of domestic rabbits

Understanding the genetic structure of domestic species provides a window into the process of domestication and motivates the design of studies aimed at making links between genotype and phenotype. Rabbits exhibit exceptional phenotypic diversity, are of great commercial value, and serve as important animal models in biomedical research. Here, we provide the first comprehensive survey of nucleotide polymorphism and linkage disequilibrium (LD) within and among rabbit breeds. We resequenced 16 genomic regions in population samples of both wild and domestic rabbits and additional 35 fragments in 150 rabbits representing six commonly used breeds. Patterns of genetic variation suggest a single origin of domestication in wild populations from France, supporting historical records that place rabbit domestication in French monasteries. Levels of nucleotide diversity both within and among breeds were ~0.2%, but only 60% of the diversity present in wild populations from France was captured by domestic rabbits. Despite the recent origin of most breeds, levels of population differentiation were high (F(ST) = 17.9%), but the majority of polymorphisms were shared and thus transferable among breeds. Coalescent simulations suggest that domestication began with a small founding population of less than 1,200 individuals. Taking into account the complex demographic history of domestication with two successive bottlenecks, two loci showed deviations that were consistent with artificial selection, including GPC4, which is known to be associated with growth rates in humans. Levels of diversity were not significantly different between autosomal and X-linked loci, providing no evidence for differential contributions of males and females to the domesticated gene pool. The structure of LD differed substantially within and among breeds. Within breeds, LD extends over large genomic distances. Markers separated by 400 kb typically showed r(2) higher than 0.2, and some LD extended up to 3,200 kb. Much less LD was found among breeds. This advantageous LD structure holds great promise for reducing the interval of association in future mapping studies.

The glypican 3-hosted murine mir717 gene: sequence conservation, seed region polymorphisms and putative targets

Mir717 (mmu-mir-717) was first reported in mouse and resides in the intron 3 of glypican 3 (Gpc3) gene. Our present study revealed that this microRNA (miRNA) gene is conserved among 26 mammalian species and harbors polymorphic sites within the mature seed region in mice. Our finding represents a rare four layer genomic overlap consisting of growth associated quantitative trait locus (QTL), body mass associated Gpc3 gene, highly conserved miRNA gene and mature miRNA seed single nucleotide polymorphism (SNP) identified in the lean mouse strain 129/Sv. Additionally, genes potentially targeted by Mir717 include 91 genes associated with obesity and related phenotypes in mammals. Our analysis provides a basis for further experiments to causally connect the identified SNP and Mir717 gene itself to obesity regulation. Furthermore, our bioinformatics analysis now enables functional annotation of Mir717 orthologs in other species, thus determining the effect of its target genes on fat-related traits.

A wild derived quantitative trait locus on mouse chromosome 2 prevents obesity

Background

The genetic architecture of multifactorial traits such as obesity has been poorly understood. Quantitative trait locus (QTL) analysis is widely used to localize loci affecting multifactorial traits on chromosomal regions. However, large confidence intervals and small phenotypic effects of identified QTLs and closely linked loci are impeding the identification of causative genes that underlie the QTLs. Here we developed five subcongenic mouse strains with overlapping and non-overlapping wild-derived genomic regions from an F2 intercross of a previously developed congenic strain, B6.Cg-Pbwg1, and its genetic background strain, C57BL/6J (B6). The subcongenic strains developed were phenotyped on low-fat standard chow and a high-fat diet to fine-map a previously identified obesity QTL. Microarray analysis was performed with Affymetrix GeneChips to search for candidate genes of the QTL.

Results

The obesity QTL was physically mapped to an 8.8-Mb region of mouse chromosome 2. The wild-derived allele significantly decreased white fat pad weight, body weight and serum levels of glucose and triglyceride. It was also resistant to the high-fat diet. Among 29 genes residing within the 8.8-Mb region, Gpd2, Upp2, Acvr1c, March7 and Rbms1 showed great differential expression in livers and/or gonadal fat pads between B6.Cg-Pbwg1 and B6 mice.

Conclusions

The wild-derived QTL allele prevented obesity in both mice fed a low-fat standard diet and mice fed a high-fat diet. This finding will pave the way for identification of causative genes for obesity. A further understanding of this unique QTL effect at genetic and molecular levels may lead to the discovery of new biological and pathologic pathways associated with obesity.

Advances in the relationship between glypican-3 and primary hepatic carcinoma

Glypican-3 (GPC3) is a membrane heparin sulfate proteoglycan that is expressed abundantly in the fetal liver, inactive in the normal adult liver and frequently reactivated in hepatocellular carcinoma (HCC). Serum soluble GPC3 is a new diagnostic and prognostic biomarker for HCC. GPC3 is also a potential target for targeted therapy of HCC. The expression of GPC3 is upregulated at the early stage of HCC, which is associated with sulfatase-2, zinc-fingers and homeoboxes 2 (ZHX2) and alpha-fetoprotein regulator 2 (AFR2). GPC3 overexpression can activate integrin, insulin-like growth factor and Wnt signaling to promote HCC development.

[Development, optimization and application of the expression analysis platform based on multiplex quantitative RT-PCR using fluorescent universal primers]

A multiplex quantitative RT-PCR technology with a universal fluorescent primer was established. This technology employs a chimeric-primer-induced-universal-primer amplification method that ensures target genes amplified in a constant ratio. This technique was cost-effective, moderate-throughput, and reliable in quantification of gene expression. It is complementary to cDNA chip, which has low quantitative accuracy , and Real-time quantitative PCR with low throughput, through improving the entire process of expression profiling analysis. Eleven genes within a QTL segment regulating mouse puberty onset on chromosome X were investigated to construct and optimize the method. The sensitivity of detection (102 copies) was determined, the concentration ratio of universal primer and chimeric forward primers (1:1) was optimized, and the accuracy and repeatability were validated. The method of Touchdown PCR with addition of universal primers significantly improved amplification of genes expressed in low abundance. After testing the expression profile of 11 genes in hypothalamus and testis in two mouse strains C3H/HeJ and C57BL/6J at the age of 15 d, one gene named PHF6 was found differentially expressed for further function analysis.

Glypican-1 controls brain size through regulation of fibroblast growth factor signaling in early neurogenesis

BACKGROUND

Cell surface heparan sulfate proteoglycans (HSPGs) act as co-receptors for multiple families of growth factors that regulate animal cell proliferation, differentiation and patterning. Elimination of heparan sulfate during brain development is known to produce severe structural abnormalities. Here we investigate the developmental role played by one particular HSPG, glypican-1 (Gpc1), which is especially abundant on neuronal cell membranes, and is the major HSPG of the adult rodent brain.

RESULTS

Mice with a null mutation in Gpc1 were generated and found to be viable and fertile. The major phenotype associated with Gpc1 loss is a highly significant reduction in brain size, with only subtle effects on brain patterning (confined to the anterior cerebellum). The brain size difference emerges very early during neurogenesis (between embryonic days 8.5 and 9.5), and remains roughly constant throughout development and adulthood. By examining markers of different signaling pathways, and the differentiation behaviors of cells in the early embryonic brain, we infer that Gpc1(-/-) phenotypes most likely result from a transient reduction in fibroblast growth factor (FGF) signaling. Through the analysis of compound mutants, we provide strong evidence that Fgf17 is the FGF family member through which Gpc1 controls brain size.

CONCLUSION

These data add to a growing literature that implicates the glypican family of HSPGs in organ size control. They also argue that, among heparan sulfate-dependent signaling molecules, FGFs are disproportionately sensitive to loss of HSPGs. Finally, because heterozygous Gpc1 mutant mice were found to have brain sizes half-way between homozygous and wild type, the data imply that endogenous HSPG levels quantitatively control growth factor signaling, a finding that is both novel and relevant to the general question of how the activities of co-receptors are exploited during development.

Overgrowth of a mouse model of Simpson-Golabi-Behmel syndrome is partly mediated by Indian hedgehog

Loss-of-function mutations of Glypican 3 (Gpc3) cause the Simpson-Golabi-Behmel overgrowth syndrome (SGBS), and developmental overgrowth is observed in Gpc3-null mice, a mouse model for SGBS. We recently reported that GPC3 inhibits Hedgehog (Hh) signalling by inducing its endocytosis and degradation. Here, we show that the developmental overgrowth observed in Gpc3-null mice is, at least in part, a consequence of the hyperactivation of the Hh pathway. We bred Gpc3-null mice with mice that are Hh signalling-deficient owing to the lack of Indian Hh (Ihh), one of the three mammalian Hhs. We found that the Gpc3-null mice showed a 29.9% overgrowth in an Ihh wild-type background, whereas an Ihh-null background partly rescues the overgrowth caused by the lack of Gpc3 as the double mutants were 19.8% bigger than the Ihh-null mice. Consistent with the role of GPC3 in Hh endocytosis and degradation, the Gpc3-null mice show increased levels of Ihh protein and signalling, but similar levels of Ihh messenger RNA.

Glypican-3 inhibits Hedgehog signaling during development by competing with patched for Hedgehog binding

Loss-of-function mutations in glypican-3 (GPC3), one of the six mammalian glypicans, causes the Simpson-Golabi-Behmel overgrowth syndrome (SGBS), and GPC3 null mice display developmental overgrowth. Because the Hedgehog signaling pathway positively regulates body size, we hypothesized that GPC3 acts as an inhibitor of Hedgehog activity during development. Here, we show that GPC3 null embryos display increased Hedgehog signaling and that GPC3 inhibits Hedgehog activity in cultured mouse embryonic fibroblasts. In addition, we report that GPC3 interacts with high affinity with Hedgehog but not with its receptor, Patched, and that GPC3 competes with Patched for Hedgehog binding. Furthermore, GPC3 induces Hedgehog endocytosis and degradation. Surprisingly, the heparan sulfate chains of GPC3 are not required for its interaction with Hedgehog. We conclude that GPC3 acts as a negative regulator of Hedgehog signaling during mammalian development and that the overgrowth observed in SGBS patients is, at least in part, the consequence of hyperactivation of the Hedgehog signaling pathway.

Heritability in the genomics era--concepts and misconceptions

Heritability allows a comparison of the relative importance of genes and environment to the variation of traits within and across populations. The concept of heritability and its definition as an estimable, dimensionless population parameter was introduced by Sewall Wright and Ronald Fisher nearly a century ago. Despite continuous misunderstandings and controversies over its use and application, heritability remains key to the response to selection in evolutionary biology and agriculture, and to the prediction of disease risk in medicine. Recent reports of substantial heritability for gene expression and new estimation methods using marker data highlight the relevance of heritability in the genomics era.

Comparative genomics for detecting human disease genes

Originally, comparative genomics was geared toward defining the synteny of genes between species. As the human genome project accelerated, there was an increase in the number of tools and means to make comparisons culminating in having the genomic sequence for a large number of organisms spanning the evolutionary tree. With this level of resolution and a long history of comparative biology and comparative genetics, it is now possible to use comparative genomics to build or select better animal models and to facilitate gene discovery. Comparative genomics takes advantage of the functional genetic information from other organisms, (vertebrates and invertebrates), to apply it to the study of human physiology and disease. It allows for the identification of genes and regulatory regions, and for acquiring knowledge about gene function. In this chapter, the current state of comparative genomics and the available tools are discussed in the context of developing animal model systems that reflect the clinical picture.

Epistasis between hyperglycemic QTLs revealed in a double congenic of the OLETF rat

Glucose homeostasis is believed to be regulated by multiple genetic components, in addition to numerous external factors. It is therefore crucial to dissect and understand what roles each causative gene plays in maintaining proper glucose metabolism. In OLETF (Otsuka Long-Evans Tokushima Fatty) rat, a model of polygenic type 2 diabetes, at least 14 quantitative trait loci (QTLs) influencing plasma glucose levels were identified. In congenic strains some of the OLETF allelic variants were shown to increase glucose levels. In this study the focus was on two of the hyperglycemic loci, Nidd1/of and Nidd2/of. Congenic rats possessing OLETF genome fragment at either locus showed similar levels of mild hyperglycemia. A newly established double congenic rat showed a further aggravation of hyperglycemia. The Nidd1/of locus was also shown to function in the reduction of plasma leptin levels and fat weights, while the Nidd2/of locus led to increased plasma insulin and fat weights. Interestingly, both plasma leptin and fat weights reverted to the control levels in the double congenic rat. These results indicate that there is an epistatic interaction between the two loci. However, it is unlikely that the abnormal level of enhanced glucose homeostasis is mediated, at least not directly, by leptin or fat mass.

Dissection of a genetic locus influencing renal function in the rat and its concordance with kidney disease loci on human chromosome 1q21

Previously, we conducted a genome scan on a population derived from the Dahl salt-sensitive hypertensive (S) and the spontaneously hypertensive rat (SHR) using urinary albumin excretion (UAE) as our primary measure of renal function. We identified 10 quantitative trait loci (QTL) linked to several renal and/or cardiovascular traits. In particular, linkage and subsequent congenic strain analysis demonstrated that the loci on chromosome 2 had a large and significant effect on UAE compared with the S rat. The present work sought to characterize the chromosome 2 congenic strain [S.SHR] by conducting a time-course analysis (week 4-20), including evaluating additional renal parameters, histology, electron microscopy, and gene expression/ pathway analysis. Throughout the time course the congenic strain consistently maintained a threefold reduction in UAE compared with S rats and was supported by the histological findings of significantly reduced glomerular, tubular and interstitial changes. Gene expression/pathway analysis performed at week 4, 12, and 20 revealed that pathways involved in cellular assembly and organization, cellular movement, and immune response were controlled differently between the S and congenic. When all the data are considered, the chromosome 2 congenic appears to attenuate renal damage primarily through an altered fibrotic response. Recombinant progeny testing was employed to reduce the QTL to approximately 1.5 cM containing several interesting candidate genes. The concordance of this rat QTL with renal disease loci on human chromosome 1q21 demonstrate that elucidating the causative gene and mechanism of the rat QTL may be of particular importance for understanding kidney disease in humans.

Genome clashes in hybrids: insights from gene expression

In interspecific hybrids, novel phenotypes often emerge from the interaction of two divergent genomes. Interactions between the two transcriptional networks are assumed to contribute to these unpredicted new phenotypes by inducing novel patterns of gene expression. Here we provide a review of the recent literature on the accumulation of regulatory incompatibilities. We review specific examples of regulatory incompatibilities reported at particular loci as well as genome-scale surveys of gene expression in interspecific hybrids. Finally, we consider and preview novel technologies that could help decipher how divergent transcriptional networks interact in hybrids between species.

Identification of positional candidate genes for body weight and adiposity in subcongenic mice

We previously constructed a congenic mouse, B6.S-D2Mit194-D2Mit311 (B6.S-2) with 27 Mb of SPRET/Ei donor DNA on distal chromosome 2 in a C57BL/6J background that captured an obesity quantitative trait locus (QTL). Mice homozygous for SPRET/Ei alleles at the donor region had decreased body weight and obesity-related phenotypes (Diament AL, Farahani P, Chiu S, Fisler J, Warden CH. Mamm Genome 15: 452-459, 2004). In this study, we constructed five overlapping subcongenics with smaller SPRET/Ei donor regions to fine map the underlying gene(s). One of the five subcongenic lines derived from the B6.S-2 founding congenic, B6.S-2A, captured the body weight and adiposity phenotypes in a donor region with a maximum size of 7.4 Mb. Homozygous SPRET/Ei donor alleles in both the founding congenic and the derived B6.S-2A subcongenic exhibited significant decreases in body weight, multiple fat pad weights, and adiposity index (total fat pad weight divided by body weight). Interval-specific microarray analysis in four tissues for donor region genes from the founding B6.S-2 congenic identified several differentially expressed genes mapping to the B6.S-2A subcongenic donor region, including prohormone convertase 2 (PC2; gene name: Pcsk2). Quantitative real-time PCR confirmed a modest decrease of PC2 expression in brains of mice homozygous for SPRET/Ei donor alleles. Analysis of the relative levels of mRNA for B6 and SPRET/Ei in heterozygous congenic mice showed differentially higher expression of the C57BL/6J allele over the SPRET/Ei allele, indicating a cis regulation of differential expression. Using subcongenic mapping, we successfully narrowed a body weight and obesity QTL interval and identified PC2 as a positional candidate gene.

A growth QTL (Pbwg1) region of mouse chromosome 2 contains closely linked loci affecting growth and body composition

Previous QTL studies have identified 24 QTLs for body weight and growth from 3 to 10 weeks after birth in an intersubspecific backcross mouse population between C57BL/6J and wild Mus musculus castaneus that has 60% of the body size of C57BL/6J. The castaneus allele at the most potent QTL (Pbwg1) on proximal chromosome 2 retards growth. In this study we have developed a congenic strain with a 44.1-Mb interval containing the castaneus allele at Pbwg1 by recurrent backcrossing to C57BL/6J. The congenic mouse developed was characterized by significantly higher body weight gain between 1 and 3 weeks of age and lower weight of white fat pads at 10 weeks of age than C57BL/6J. However, no clear difference in body weight at 1-10 weeks of age was observed between congenic and C57BL/6J strains. QTL analysis with 269 F(2) mice between the two strains did not identify any QTLs for body weight at 1, 3, 6, and 10 weeks of age, but it discovered eight closely linked QTLs affecting body weight gain from 1 to 3 weeks of age, lean body weight, weight of white fat pads, and body length within the Pbwg1 region. The castaneus alleles at all fat pad QTLs reduced the phenotypes, whereas at the remaining growth and body composition QTLs, they increased the trait values. These results illustrate that Pbwg1, which initially appeared to be a single locus, was resolved into several loci with opposite effects on the composition traits of overall body weight. This gives a reason for the loss of the Pbwg1 effect found in the original backcross population.

Fine mapping reveals sex bias in quantitative trait loci affecting growth, skeletal size and obesity-related traits on mouse chromosomes 2 and 11

Previous speed congenic analysis has suggested that the expression of growth and obesity quantitative trait loci (QTL) on distal mouse chromosomes (MMU) 2 and 11, segregating between the CAST/EiJ (CAST) and C57BL/6J-hg/hg (HG) strains, is dependent on sex. To confirm, fine map, and further evaluate QTL x sex interactions, we constructed congenic by recipient F2 crosses for the HG.CAST-(D2Mit329-D2Mit457)N(6) (HG2D) and HG.CAST-(D11Mit260-D11Mit255)N(6) (HG11) congenic strains. Over 700 F2 mice were densely genotyped and phenotyped for a panel of 40 body and organ weight, skeletal length, and obesity-related traits at 9 weeks of age. Linkage analysis revealed 20 QTL affecting a representative subset of phenotypes in HG2DF2 and HG11F2 mice. The effect of sex was quantified by comparing two linear models: the first model included sex as an additive covariate and the second incorporated sex as an additive and an interactive covariate. Of the 20 QTL, 8 were sex biased, sex specific, or sex antagonistic. Most traits were regulated by single QTL; however, two closely linked loci were identified for five traits in HG2DF2 mice. Additionally, the confidence intervals for most QTL were significantly reduced relative to the original mapping results, setting the stage for quantitative trait gene (QTG) discovery. These results highlight the importance of assessing the contribution of sex in complex trait analyses.

Ancient and recent positive selection transformed opioid cis-regulation in humans

Changes in the cis-regulation of neural genes likely contributed to the evolution of our species' unique attributes, but evidence of a role for natural selection has been lacking. We found that positive natural selection altered the cis-regulation of human prodynorphin, the precursor molecule for a suite of endogenous opioids and neuropeptides with critical roles in regulating perception, behavior, and memory. Independent lines of phylogenetic and population genetic evidence support a history of selective sweeps driving the evolution of the human prodynorphin promoter. In experimental assays of chimpanzee-human hybrid promoters, the selected sequence increases transcriptional inducibility. The evidence for a change in the response of the brain's natural opioids to inductive stimuli points to potential human-specific characteristics favored during evolution. In addition, the pattern of linked nucleotide and microsatellite variation among and within modern human populations suggests that recent selection, subsequent to the fixation of the human-specific mutations and the peopling of the globe, has favored different prodynorphin cis-regulatory alleles in different parts of the world.

RNA world - the dark matter of evolutionary genomics

For a long time, molecular evolutionary biologists have been focused on DNA and proteins, whereas RNA has lived in the shadow of its famous chemical cousins as a mere intermediary. Although this perspective has begun to change since genome-wide transcriptional profiling was successfully extended to evolutionary biology, it still echoes in evolutionary literature. In this mini-review, new developments of RNA biochemistry and transcriptomics are brought to the attention of evolutionary biologists. In particular, the unexpected abundance and functional significance of noncoding RNAs is briefly reviewed. Noncoding RNAs control a remarkable range of biological pathways and processes, all with obvious fitness consequences, such as initiation of translation, mRNA abundance, transposon jumping, chromosome architecture, stem cell maintenance, development of brain and muscles, insulin secretion, cancerogenesis and plant resistance to viral infections.

Marker-assisted selection based on a multi-trait economic index in chicken: experimental results and simulation

A method proposed herein allows simultaneous selection for several production traits, taking into consideration their marginal economic values (i.e. the economic value of a trait's additional unit). This economic index-marker assisted selection (EI-MAS) method is based on the calculation of the predicted economic breeding value (BV), using information on DNA markers that have previously been found to be associated with relevant quantitative trait loci. Based on the proposed method, results with real birds showed that sire progeny performance was significantly correlated with expected performance (r = 0.61-0.76; P = 0.03-0.01). Simulation analysis using a computer program written specifically for this purpose suggested that the relative advantage of EI-MAS would be large for traits with low heritability values. As expected, the response to EI-MAS was higher when the map distance between the marker and the quantitative trait gene was small, and vice versa. A large number of distantly located markers, spread 10 cM apart, yielded higher response to selection than a small number of closely located markers spread 3 cM apart. Additionally, the response to EI-MAS was higher when a large number (ca.150) of progeny was used for the prediction equation.

Phenotypic characterisation of extreme growth-selected mouse lines: An important prerequisite for future QTL analysis

A unique set of seven mouse lines, long-term selected for high growth, from different laboratories around the world has been comprehensively compared to evaluate these resources for future QTL and gene mapping for growth traits. The heaviest line (DUH) was 40% (males) to 44% (females) heavier than the smallest line (ROH) at birth, and 105% (males) to 114% (females) heavier at 98 d. Body conformation (body length and width, body areas), body composition (dry matter, fat, fatty acid composition, organ weights), and skeletal muscle cellularity also differed substantially. DUH was more than 20% longer (12.3 cm) compared to the shortest line ROH (9.7 cm). DAH (22.5%) had the highest percentage of gonadal fat and the leanest was BEH (7.7%). Line BEH (0.49 g) showed the highest weight for the left M. rectus femoris, which was 2.1 times higher, compared to ROH (0.23 g). These results suggest that different alleles, and possibly different physiological pathways, have contributed to the selection response in the different lines. Therefore these selection lines are an important tool with which to identify the genetic and physiological basis of growth as they may contain many, if not all, growth promoting alleles.

Allele-specific KRT1 expression is a complex trait

The differential expression of alleles occurs commonly in humans and is likely an important genetic factor underlying heritable differences in phenotypic traits. Understanding the molecular basis of allelic expression differences is thus an important challenge. Although many genes have been shown to display differential allelic expression, this is the first study to examine in detail the cumulative effects of multiple cis-regulatory polymorphisms responsible for allele-specific expression differences. We have used a variety of experimental approaches to identify and characterize cis-regulatory polymorphisms responsible for the extreme allele-specific expression differences of keratin-1 (KRT1) in human white blood cells. The combined data from our analyses provide strong evidence that the KRT1 allelic expression differences result from the haplotypic combinations and interactions of five cis-regulatory single nucleotide polymorphisms (SNPs) whose alleles differ in their affinity to bind transcription factors and modulate KRT1 promoter activity. Two of these cis-regulatory SNPs bind transcriptional activators with the alleles on the high-expressing KRT1 haplotype pattern having a higher affinity than the alleles on the low-expressing haplotype pattern. In contrast, the other three cis-regulatory SNPs bind transcriptional inhibitors with the alleles on the low-expressing haplotype pattern having a higher affinity than the alleles on the high-expressing haplotype pattern. Our study provides important new insights into the degree of complexity that the cis-regulatory sequences responsible for allele-specific transcriptional regulation have. These data suggest that allelic expression differences result from the cumulative contribution of multiple DNA sequence polymorphisms, with each having a small effect, and that allele-specific expression can thus be viewed as a complex trait.

Despite the fact that all humans share nearly identical DNA sequences, individuals exhibit tremendous variation in heritable traits, such as height, weight, and skin texture. Recent evidence suggests that expression level differences between different copies (alleles) of a gene contribute to these observed differences in heritable traits. Currently, the mechanisms underlying allele-expression level differences are poorly understood. In this report the authors identified and characterized a set of five single nucleotide polymorphisms (SNPs) contributing to extreme expression differences between keratin-1 (KRT1) alleles in humans. Each of the five SNPs is found in a different regulatory sequence in the proximity of KRT1. The SNPs cause different copies of the five regulatory sequences to differ in their affinities to bind transcription factors controlling KRT1 expression. The extreme KRT1 allele-expression level differences result from the cumulative contributions of these five SNPs which are tightly linked and inherited in two common fixed sets, a low- and a high-expressing set. The study provides important new insights into the complexities of the mechanisms underlying allele-expression level differences. These complexities may explain the difficulties researchers frequently encounter when trying to discover the “causative SNP” in an interval identified as associated with an inherited trait in a genetic study.

PAPPA2, an enzyme that cleaves an insulin-like growth-factor-binding protein, is a candidate gene for a quantitative trait locus affecting body size in mice

Identifying genes responsible for quantitative variation remains a major challenge. We previously identified a quantitative trait locus (QTL) affecting body size that segregated between two inbred strains of mice. By fine mapping, we have refined the location of this QTL to a genomic region containing only four protein-coding genes. One of these genes, PAPPA2, is a strong candidate because it codes for an enzyme that cleaves insulin-like growth-factor-binding protein 5 (IGFBP-5), an important stimulator of bone formation. Among littermates that segregate only for the four-gene region, we show that the QTL has a significant effect on the circulating levels of IGFBP-5 and IGFBP-3 (the latter subject to limited degradation by PAPPA2), but not on levels of IGFBP-2 and IGFBP-4, which are not cleaved by PAPPA2. There are 14 nonsynonymous SNPs among QTL alleles, which may affect the activity of the translated protein. The refinement of the target region to four genes and the finding that the QTL affects IGFBP-5 levels suggest that PAPPA2 may be involved with normal postnatal growth. Our mapping results also illustrate the potentially fractal nature of QTL: as we mapped our QTL with increasing resolution, what appeared to be a single QTL resolved into three closely linked QTL (previous work), and then one of these was further dissected into two in this study.

Association of the OLR1 Gene with Milk Composition in Holstein Dairy Cattle

Oxidized low-density lipoprotein receptor (OLR1) is the major protein that binds, internalizes, and degrades oxidized low-density lipoprotein. The role of OLR1 in lipid metabolism and the results of previous whole-genome scan studies prompted the investigation of OLR1 as a candidate gene affecting milk composition traits. Direct cDNA and genomic sequencing of OLR1 revealed 2 single nucleotide polymorphisms (SNP) in exon 4, 5 SNP in intron 4, and 1 in the 3' untranslated region (UTR). Four intragenic haplotypes comprising SNP positions 7,160, 7,161, 7,278, 7,381, 7,409, 7,438, 7,512, and 8,232 were inferred. Haplotype analysis showed that one of the haplotypes was associated with a significant increase in fat yield and fat percentage. Single SNP analysis showed that allele C of SNP 8,232 (in the 3'-UTR) had significant effects on fat yield and fat percentage, whereas SNP 7,160 and 7,161 (in exon 4) had no significant effects. Both single SNP and haplotype analyses indicate that SNP 8,232 in the 3'-UTR is associated with milk fat yield and percentage and it may be in linkage disequilibrium with the functional polymorphism. To provide support for the hypothesis that SNP 8,232 is responsible for OLR1 expression, OLR1 expression levels in individuals bearing different genotypes were assessed. It was found that OLR1 expression was reduced in genotype AA individuals compared with CC and AC individuals, suggesting that A at position 8,232 may be the nucleotide causing decreased OLR1 expression. The 3'-UTR polymorphism found in this study might control translation or stability of OLR1 mRNA.

The human obesity gene map: the 2005 update

This paper presents the 12th update of the human obesity gene map, which incorporates published results up to the end of October 2005. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, transgenic and knockout murine models relevant to obesity, quantitative trait loci (QTL) from animal cross-breeding experiments, association studies with candidate genes, and linkages from genome scans is reviewed. As of October 2005, 176 human obesity cases due to single-gene mutations in 11 different genes have been reported, 50 loci related to Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. There are 244 genes that, when mutated or expressed as transgenes in the mouse, result in phenotypes that affect body weight and adiposity. The number of QTLs reported from animal models currently reaches 408. The number of human obesity QTLs derived from genome scans continues to grow, and we now have 253 QTLs for obesity-related phenotypes from 61 genome-wide scans. A total of 52 genomic regions harbor QTLs supported by two or more studies. The number of studies reporting associations between DNA sequence variation in specific genes and obesity phenotypes has also increased considerably, with 426 findings of positive associations with 127 candidate genes. A promising observation is that 22 genes are each supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. The electronic version of the map with links to useful publications and relevant sites can be found at http://obesitygene.pbrc.edu.

Analysis of allelic differential expression in human white blood cells

Allelic variation of gene expression is common in humans, and is of interest because of its potential contribution to variation in heritable traits. To identify human genes with allelic expression differences, we genotype DNA and examine mRNA isolated from the white blood cells of 12 unrelated individuals using oligonucleotide arrays containing 8406 exonic SNPs. Of the exonic SNPs, 1983, located in 1389 genes, are both expressed in the white blood cells and heterozygous in at least one of the 12 individuals, and thus can be examined for differential allelic expression. Of the 1389 genes, 731 (53%) show allele expression differences in at least one individual. To gain insight into the regulatory mechanisms governing allelic expression differences, we analyze a set of 60 genes containing exonic SNPs that are heterozygous in three or more samples, and for which all heterozygotes display differential expression. We find three patterns of allelic expression, suggesting different underlying regulatory mechanisms. Exonic SNPs in three of the 60 genes are monoallelically expressed in the human white blood cells, and when examined in families show expression of only the maternal copy, consistent with regulation by imprinting. Approximately one-third of the genes have the same allele expressed more highly in all heterozygotes, suggesting that their regulation is predominantly influenced by cis-elements in strong linkage disequilibrium with the assayed exonic SNP. The remaining two-thirds of the genes have different alleles expressed more highly in different heterozygotes, suggesting that their expression differences are influenced by factors not in strong linkage disequilibrium with the assayed exonic SNP.