Statistical Inference of Selection and Divergence from a Time-Dependent Poisson Random Field Model

We apply a recently developed time-dependent Poisson random field model to aligned DNA sequences from two related biological species to estimate selection coefficients and divergence time. We use Markov chain Monte Carlo methods to estimate species divergence time and selection coefficients for each locus. The model assumes that the selective effects of non-synonymous mutations are normally distributed across genetic loci but constant within loci, and synonymous mutations are selectively neutral. In contrast with previous models, we do not assume that the individual species are at population equilibrium after divergence. Using a data set of 91 genes in two Drosophila species, D. melanogaster and D. simulans, we estimate the species divergence time (or 1.68 million years, assuming the haploid effective population size years) and a mean selection coefficient per generation . Although the average selection coefficient is positive, the magnitude of the selection is quite small. Results from numerical simulations are also presented as an accuracy check for the time-dependent model

LGR6 Is a High Affinity Receptor of R-Spondins and Potentially Functions as a Tumor Suppressor

BACKGROUND: LGR6 (leucine-rich repeat containing, G protein-coupled receptor 6) is a member of the rhodopsin-like seven transmembrane domain receptor superfamily with the highest homology to LGR4 and LGR5. LGR6 was found as one of the novel genes mutated in colon cancer through total exon sequencing and its promoter region is hypermethylated in 20-50% of colon cancer cases. In the skin, LGR6 marks a population of stem cells that can give rise to all cell lineages. Recently, we and others demonstrated that LGR4 and LGR5 function as receptors of R-spondins to potentiate Wnt/β-catenin signaling. However, the binding affinity and functional response of LGR6 to R-spondins, and the activity of colon cancer mutants of LGR6 have not been determined. PRINCIPAL FINDINGS: We found that LGR6 also binds and responds to R-spondins 1-3 with high affinity to enhance Wnt/β-catenin signaling through increased LRP6 phosphorylation. Similar to LGR4 and LGR5, LGR6 is not coupled to heterotrimeric G proteins or to β-arrestin following R-spondin stimulation. Functional and expression analysis of three somatic mutations identified in colon cancer samples indicates that one mutant fails to bind and respond to R-spondin (loss-of-function), but the other two have no significant effect on receptor function. Overexpression of wild-type LGR6 in HeLa cells leads to increased cell migration following co-treatment with R-spondin1 and Wnt3a when compared to vector control cells or cells overexpressing the loss-of-function mutant. CONCLUSIONS: LGR6 is a high affinity receptor for R-spondins 1-3 and potentially functions as a tumor suppressor despite its positive effect on Wnt/β-catenin signaling

A high resolution genome-wide scan for significant selective sweeps: an application to pooled sequence data in laying chickens

In most studies aimed at localizing footprints of past selection, outliers at tails of the empirical distribution of a given test statistic are assumed to reflect locus-specific selective forces. Significance cutoffs are subjectively determined, rather than being related to a clear set of hypotheses. Here, we define an empirical p-value for the summary statistic by means of a permutation method that uses the observed SNP structure in the real data. To illustrate the methodology, we applied our approach to a panel of 2.9 million autosomal SNPs identified from re-sequencing a pool of 15 individuals from a brown egg layer line. We scanned the genome for local reductions in heterozygosity, suggestive of selective sweeps. We also employed a modified sliding window approach that accounts for gaps in the sequence and increases scanning resolution by moving the overlapping windows by steps of one SNP only, and suggest to call this a "creeping window" strategy. The approach confirmed selective sweeps in the region of previously described candidate genes, i.e. TSHR, PRL, PRLHR, INSR, LEPR, IGF1, and NRAMP1 when used as positive controls. The genome scan revealed 82 distinct regions with strong evidence of selection (genome-wide p-value<0.001), including genes known to be associated with eggshell structure and immune system such as CALB1 and GAL cluster, respectively. A substantial proportion of signals was found in poor gene content regions including the most extreme signal on chromosome 1. The observation of multiple signals in a highly selected layer line of chicken is consistent with the hypothesis that egg production is a complex trait controlled by many genes

The Spatial Distribution of LGR5+ Cells Correlates With Gastric Cancer Progression

In this study we tested the prevalence, histoanatomical distribution and tumour biological significance of the Wnt target protein and cancer stem cell marker LGR5 in tumours of the human gastrointestinal tract. Differential expression of LGR5 was studied on transcriptional (real-time polymerase chain reaction) and translational level (immunohistochemistry) in malignant and corresponding non-malignant tissues of 127 patients comprising six different primary tumour sites, i.e. oesophagus, stomach, liver, pancreas, colon and rectum. The clinico-pathological significance of LGR5 expression was studied in 100 patients with gastric carcinoma (GC). Non-neoplastic tissue usually harboured only very few scattered LGR5+ cells. The corresponding carcinomas of the oesophagus, stomach, liver, pancreas, colon and rectum showed significantly more LGR5+ cells as well as significantly higher levels of LGR5-mRNA compared with the corresponding non-neoplastic tissue. Double staining experiments revealed a coexpression of LGR5 with the putative stem cell markers CD44, Musashi-1 and ADAM17. Next we tested the hypothesis that the sequential changes of gastric carcinogenesis, i.e. chronic atrophic gastritis, intestinal metaplasia and invasive carcinoma, are associated with a reallocation of the LGR5+ cells. Interestingly, the spatial distribution of LGR5 changed: in non-neoplastic stomach mucosa, LGR5+ cells were found predominantly in the mucous neck region; in intestinal metaplasia LGR5+ cells were localized at the crypt base, and in GC LGR5+ cells were present at the luminal surface, the tumour centre and the invasion front. The expression of LGR5 in the tumour centre and invasion front of GC correlated significantly with the local tumour growth (T-category) and the nodal spread (N-category). Furthermore, patients with LGR5+ GCs had a shorter median survival (28.0±8.6 months) than patients with LGR5− GCs (54.5±6.3 months). Our results show that LGR5 is differentially expressed in gastrointestinal cancers and that the spatial histoanatomical distribution of LGR5+ cells has to be considered when their tumour biological significance is sought

Multiple Roles of Integrin-Linked Kinase in Epidermal Development, Maturation and Pigmentation Revealed by Molecular Profiling

Integrin-linked kinase (ILK) is an important scaffold protein that mediates a variety of cellular responses to integrin stimulation by extracellular matrix proteins. Mice with epidermis-restricted inactivation of the Ilk gene exhibit pleiotropic phenotypic defects, including impaired hair follicle morphogenesis, reduced epidermal adhesion to the basement membrane, compromised epidermal integrity, as well as wasting and failure to thrive leading to perinatal death. To better understand the underlying molecular mechanisms that cause such a broad range of alterations, we investigated the impact of Ilk gene inactivation on the epidermis transcriptome. Microarray analysis showed over 700 differentially regulated mRNAs encoding proteins involved in multiple aspects of epidermal function, including keratinocyte differentiation and barrier formation, inflammation, regeneration after injury, and fundamental epidermal developmental pathways. These studies also revealed potential effects on genes not previously implicated in ILK functions, including those important for melanocyte and melanoblast development and function, regulation of cytoskeletal dynamics, and homeobox genes. This study shows that ILK is a critical regulator of multiple aspects of epidermal function and homeostasis, and reveals the previously unreported involvement of ILK not only in epidermal differentiation and barrier formation, but also in melanocyte genesis and function

Wnt5a Regulates Midbrain Dopaminergic Axon Growth and Guidance

During development, precise temporal and spatial gradients are responsible for guiding axons to their appropriate targets. Within the developing ventral midbrain (VM) the cues that guide dopaminergic (DA) axons to their forebrain targets remain to be fully elucidated. Wnts are morphogens that have been identified as axon guidance molecules. Several Wnts are expressed in the VM where they regulate the birth of DA neurons. Here, we describe that a precise temporo-spatial expression of Wnt5a accompanies the development of nigrostriatal projections by VM DA neurons. In mice at E11.5, Wnt5a is expressed in the VM where it was found to promote DA neurite and axonal growth in VM primary cultures. By E14.5, when DA axons are approaching their striatal target, Wnt5a causes DA neurite retraction in primary cultures. Co-culture of VM explants with Wnt5a-overexpressing cell aggregates revealed that Wnt5a is capable of repelling DA neurites. Antagonism experiments revealed that the effects of Wnt5a are mediated by the Frizzled receptors and by the small GTPase, Rac1 (a component of the non-canonical Wnt planar cell polarity pathway). Moreover, the effects were specific as they could be blocked by Wnt5a antibody, sFRPs and RYK-Fc. The importance of Wnt5a in DA axon morphogenesis was further verified in Wnt5a−/− mice, where fasciculation of the medial forebrain bundle (MFB) as well as the density of DA neurites in the MFB and striatal terminals were disrupted. Thus, our results identify a novel role of Wnt5a in DA axon growth and guidance

Non-equivalence of Wnt and R-spondin ligands during Lgr5+ intestinal stem-cell self-renewal

The canonical Wnt/β-catenin signaling pathway governs diverse developmental, homeostatic and pathologic processes. Palmitoylated Wnt ligands engage cell surface Frizzled (Fzd) receptors and Lrp5/6 co-receptors enabling β-catenin nuclear translocation and Tcf/Lef-dependent gene transactivation1–3. Mutations in Wnt downstream signaling components have revealed diverse functions presumptively attributed to Wnt ligands themselves, although direct attribution remains elusive, as complicated by redundancy between 19 mammalian Wnts and 10 Fzds1 and Wnt hydrophobicity2,3. For example, individual Wnt ligand mutations have not revealed homeostatic phenotypes in the intestinal epithelium4, an archetypal canonical Wnt pathway-dependent rapidly self-renewing tissue whose regeneration is fueled by proliferative crypt Lgr5+ intestinal stem cells (ISCs)5–9. R-spondin ligands (Rspo1–4) engage distinct Lgr4-6 and Rnf43/Znrf3 receptor classes10–13, markedly potentiate canonical Wnt/β-catenin signaling and induce intestinal organoid growth in vitro and Lgr5+ ISCs in vivo8,14–17. However, the interchangeability, functional cooperation and relative contributions of Wnt versus Rspo ligands to in vivo canonical Wnt signaling and ISC biology remain unknown. Here, we deconstructed functional roles of Wnt versus Rspo ligands in the intestinal crypt stem cell niche. We demonstrate that the default fate of Lgr5+ ISCs is lineage commitment, escape from which requires both Rspo and Wnt ligands. However, gain-of-function studies using Rspo versus a novel non-lipidated Wnt analog reveal qualitatively distinct, non-interchangeable roles for these ligands in ISCs. Wnts are insufficient to induce Lgr5+ ISC self-renewal, but rather confer a basal competency by maintaining Rspo receptor expression that enables Rspo to actively drive and specify the extent of stem cell expansion. This functionally non-equivalent yet cooperative interplay between Wnt and Rspo ligands establishes a molecular precedent for regulation of mammalian stem cells by distinct priming and self-renewal factors, with broad implications for precision control of tissue regeneration

The RSPO–LGR4/5–ZNRF3/RNF43 module controls liver zonation and size

LGR4/5 receptors and their cognate RSPO ligands potentiate Wnt/β-catenin signalling and promote proliferation and tissue homeostasis in epithelial stem cell compartments. In the liver, metabolic zonation requires a Wnt/β-catenin signalling gradient, but the instructive mechanism controlling its spatiotemporal regulation is not known. We have now identified the RSPO-LGR4/5-ZNRF3/RNF43 module as a master regulator of Wnt/β-catenin-mediated metabolic liver zonation. Liver-specific LGR4/5 loss of function (LOF) or RSPO blockade disrupted hepatic Wnt/β-catenin signalling and zonation. Conversely, pathway activation in ZNRF3/RNF43 LOF mice or with recombinant RSPO1 protein expanded the hepatic Wnt/β-catenin signalling gradient in a reversible and LGR4/5-dependent manner. Recombinant RSPO1 protein increased liver size and improved liver regeneration, whereas LGR4/5 LOF caused the opposite effects, resulting in hypoplastic livers. Furthermore, we show that LGR4(+) hepatocytes throughout the lobule contribute to liver homeostasis without zonal dominance. Taken together, our results indicate that the RSPO-LGR4/5-ZNRF3/RNF43 module controls metabolic liver zonation and is a hepatic growth/size rheostat during development, homeostasis and regeneration

Holoprosencephaly

Holoprosencephaly (HPE) is a complex brain malformation resulting from incomplete cleavage of the prosencephalon, occurring between the 18th and the 28th day of gestation and affecting both the forebrain and the face. It is estimated to occur in 1/16,000 live births and 1/250 conceptuses. Three ranges of increasing severity are described: lobar, semi-lobar and alobar HPE. Another milder subtype of HPE called middle interhemispheric variant (MIHF) or syntelencephaly is also reported. In most of the cases, facial anomalies are observed in HPE, like cyclopia, proboscis, median or bilateral cleft lip/palate in severe forms, ocular hypotelorism or solitary median maxillary central incisor in minor forms. These latter midline defects can occur without the cerebral malformations and then are called microforms. Children with HPE have many medical problems: developmental delay and feeding difficulties, epilepsy, instability of temperature, heart rate and respiration. Endocrine disorders like diabetes insipidus, adrenal hypoplasia, hypogonadism, thyroid hypoplasia and growth hormone deficiency are frequent. To date, seven genes have been positively implicated in HPE: Sonic hedgehog (SHH), ZIC2, SIX3, TGIF, PTCH, GLI2 and TDGF1. A molecular diagnosis can be performed by gene sequencing and allele quantification for the four main genes SHH, ZIC2, SIX3 and TGIF. Major rearrangements of the subtelomeres can also be identified by multiplex ligation-dependent probe amplification (MLPA). Nevertheless, in about 70% of cases, the molecular basis of the disease remains unknown, suggesting the existence of several other candidate genes or environmental factors. Consequently, a "multiple-hit hypothesis" of genetic and/or environmental factors (like maternal diabetes) has been proposed to account for the extreme clinical variability. In a practical approach, prenatal diagnosis is based on ultrasound and magnetic resonance imaging (MRI) rather than on molecular diagnosis. Treatment is symptomatic and supportive, and requires a multidisciplinary management. Child outcome depends on the HPE severity and the medical and neurological complications associated. Severely affected children have a very poor prognosis. Mildly affected children may exhibit few symptoms and may live a normal life