Homeostatic Control of Sebaceous Glands by Innate Lymphoid Cells Regulates Commensal Bacteria Equilibrium.

Immune cells and epithelium form sophisticated barrier systems in symbiotic relationships with microbiota. Evidence suggests that immune cells can sense microbes through intact barriers, but regulation of microbial commensalism remain largely unexplored. Here, we uncovered spatial compartmentalization of skin-resident innate lymphoid cells (ILCs) and modulation of sebaceous glands by a subset of RORγt+ ILCs residing within hair follicles in close proximity to sebaceous glands. Their persistence in skin required IL-7 and thymic stromal lymphopoietin, and localization was dependent on the chemokine receptor CCR6. ILC subsets expressed TNF receptor ligands, which limited sebocyte growth by repressing Notch signaling pathway. Consequently, loss of ILCs resulted in sebaceous hyperplasia with increased production of antimicrobial lipids and restricted commensalism of Gram-positive bacterial communities. Thus, epithelia-derived signals maintain skin-resident ILCs that regulate microbial commensalism through sebaceous gland-mediated tuning of the barrier surface, highlighting an immune-epithelia circuitry that facilitates host-microbe symbiosis.

Genetic variants regulating immune cell levels in health and disease.

The complex network of specialized cells and molecules in the immune system has evolved to defend against pathogens, but inadvertent immune system attacks on "self" result in autoimmune disease. Both genetic regulation of immune cell levels and their relationships with autoimmunity are largely undetermined. Here, we report genetic contributions to quantitative levels of 95 cell types encompassing 272 immune traits, in a cohort of 1,629 individuals from four clustered Sardinian villages. We first estimated trait heritability, showing that it can be substantial, accounting for up to 87% of the variance (mean 41%). Next, by assessing ∼8.2 million variants that we identified and confirmed in an extended set of 2,870 individuals, 23 independent variants at 13 loci associated with at least one trait. Notably, variants at three loci (HLA, IL2RA, and SH2B3/ATXN2) overlap with known autoimmune disease associations. These results connect specific cellular phenotypes to specific genetic variants, helping to explicate their involvement in disease.

GWAS of genetic factors affecting white blood cell morphological parameters in Sardinians uncovers influence of chromosome 11 innate immunity gene cluster on eosinophil morphology.

Few genome-wide association studies (GWAS) analyzing genetic regulation of morphological traits of white blood cells have been reported. We carried out a GWAS of 12 morphological traits in 869 individuals from the general population of Sardinia, Italy. These traits, included measures of cell volume, conductivity and light scatter in four white-cell populations (eosinophils, lymphocytes, monocytes, neutrophils). This analysis yielded seven statistically significant signals, four of which were novel (four novel, PRG2, P2RX3, two of CDK6). Five signals were replicated in the independent INTERVAL cohort of 11 822 individuals. The most interesting signal with large effect size on eosinophil scatter (P-value = 8.33 x 10-32, beta = -1.651, se = 0.1351) falls within the innate immunity cluster on chromosome 11, and is located in the PRG2 gene. Computational analyses revealed that a rare, Sardinian-specific PRG2:p.Ser148Pro mutation modifies PRG2 amino acid contacts and protein dynamics in a manner that could possibly explain the changes observed in eosinophil morphology. Our discoveries shed light on genetics of morphological traits. For the first time, we describe such large effect size on eosinophils morphology that is relatively frequent in Sardinian population.

Risk of severe COVID-19 infection among adults with prior exposure to children.

Susceptibility and severity of COVID-19 infection vary widely. Prior exposure to endemic coronaviruses, common in young children, may protect against SARS-CoV-2. We evaluated risk of severe COVID-19 among adults with and without exposure to young children in a large, integrated healthcare system. Adults with children 0-5 years were matched 1:1 to adults with children 6-11 years, 12-18 years, and those without children based upon a COVID-19 propensity score and risk factors for severe COVID-19. COVID-19 infections, hospitalizations, and need for intensive care unit (ICU) were assessed in 3,126,427 adults, of whom 24% (N = 743,814) had children 18 years or younger, and 8.8% (N = 274,316) had a youngest child 0-5 years. After 1:1 matching, propensity for COVID-19 infection and risk factors for severe COVID-19 were well balanced between groups. Rates of COVID-19 infection were slightly higher for adults with exposure to older children (incident risk ratio, 1.09, 95% confidence interval, [1.05-1.12] and IRR 1.09 [1.05-1.13] for adults with children 6-11 and 12-18, respectively), compared to those with children 0-5 years, although no difference in rates of COVID-19 illness requiring hospitalization or ICU admission was observed. However, adults without exposure to children had lower rates of COVID-19 infection (IRR 0.85, [0.83-0.87]) but significantly higher rates of COVID-19 hospitalization (IRR 1.49, [1.29-1.73]) and hospitalization requiring ICU admission (IRR 1.76, [1.19-2.58]) compared to those with children aged 0-5. In a large, real-world population, exposure to young children was associated with less severe COVID-19 illness. Endemic coronavirus cross-immunity may play a role in protection against severe COVID-19.

A fast linkage method for population GWAS cohorts with related individuals.

Linkage analysis, a class of methods for detecting co-segregation of genomic segments and traits in families, was used to map disease-causing genes for decades before genotyping arrays and dense SNP genotyping enabled genome-wide association studies in population samples. Population samples often contain related individuals, but the segregation of alleles within families is rarely used because traditional linkage methods are computationally inefficient for larger datasets. Here, we describe Population Linkage, a novel application of Haseman-Elston regression as a method of moments estimator of variance components and their standard errors. We achieve additional computational efficiency by using modern methods for detection of IBD segments and variance component estimation, efficient preprocessing of input data, and minimizing redundant numerical calculations. We also refined variance component models to account for the biases in population-scale methods for IBD segment detection. We ran Population Linkage on four blood lipid traits in over 70,000 individuals from the HUNT and SardiNIA studies, successfully detecting 25 known genetic signals. One notable linkage signal that appeared in both was for low-density lipoprotein (LDL) cholesterol levels in the region near the gene APOE (LOD = 29.3, variance explained = 4.1%). This is the region where the missense variants rs7412 and rs429358, which together make up the ε2, ε3, and ε4 alleles each account for 2.4% and 0.8% of variation in circulating LDL cholesterol. Our results show the potential for linkage analysis and other large-scale applications of method of moments variance components estimation.

Actively transcribed rDNA and distal junction (DJ) sequence are involved in association of NORs with nucleoli.

Although they are organelles without a limiting membrane, nucleoli have an exclusive structure, built upon the rDNA-rich acrocentric short arms of five human chromosomes (nucleolar organizer regions or NORs). This has raised the question: what are the structural features of a chromosome required for its inclusion in a nucleolus? Previous work has suggested that sequences adjacent to the tandemly repeated rDNA repeat units (DJ, distal junction sequence) may be involved, and we have extended such studies by addressing several issues related to the requirements for the association of NORs with nucleoli. We exploited both a set of somatic cell hybrids containing individual human acrocentric chromosomes and a set of Human Artificial Chromosomes (HACs) carrying different parts of a NOR, including an rDNA unit or DJ or PJ (proximal junction) sequence. Association of NORs with nucleoli was increased when constituent rDNA was transcribed and may be also affected by the status of heterochromatin blocks formed next to the rDNA arrays. Furthermore, our data suggest that a relatively small size DJ region, highly conserved in evolution, is also involved, along with the rDNA repeats, in the localization of p-arms of acrocentric chromosomes in nucleoli. Thus, we infer a cooperative action of rDNA sequence-stimulated by its activity-and sequences distal to rDNA contributing to incorporation into nucleoli. Analysis of NOR sequences also identified LncRNA_038958 in the DJ, a candidate transcript with the region of the suggested promoter that is located close to the DJ/rDNA boundary and contains CTCF binding sites. This LncRNA may affect RNA Polymerase I and/or nucleolar activity. Our findings provide the basis for future studies to determine which RNAs and proteins interact critically with NOR sequences to organize the higher-order structure of nucleoli and their function in normal cells and pathological states.

Evidence of Polygenic Adaptation in Sardinia at Height-Associated Loci Ascertained from the Biobank Japan.

Adult height is one of the earliest putative examples of polygenic adaptation in humans. However, this conclusion was recently challenged because residual uncorrected stratification from large-scale consortium studies was considered responsible for the previously noted genetic difference. It thus remains an open question whether height loci exhibit signals of polygenic adaptation in any human population. We re-examined this question, focusing on one of the shortest European populations, the Sardinians, in addition to mainland European populations. We utilized height-associated loci from the Biobank Japan (BBJ) dataset to further alleviate concerns of biased ascertainment of GWAS loci and showed that the Sardinians remain significantly shorter than expected under neutrality (∼0.22 standard deviation shorter than Utah residents with ancestry from northern and western Europe [CEU] on the basis of polygenic height scores, p = 3.89 × 10-4). We also found the trajectory of polygenic height scores between the Sardinian and the British populations diverged over at least the last 10,000 years (p = 0.0082), consistent with a signature of polygenic adaptation driven primarily by the Sardinian population. Although the polygenic score-based analysis showed a much subtler signature in mainland European populations, we found a clear and robust adaptive signature in the UK population by using a haplotype-based statistic, the trait singleton density score (tSDS), driven by the height-increasing alleles (p = 9.1 × 10-4). In summary, by ascertaining height loci in a distant East Asian population, we further supported the evidence of polygenic adaptation at height-associated loci among the Sardinians. In mainland Europeans, the adaptive signature was detected in haplotype-based analysis but not in polygenic score-based analysis.

Polygenic risk score and biochemical/environmental variables predict a low-risk profile of age-related macular degeneration in Sardinia.

PURPOSE: To ascertain the prevalence and clinical and genetic features of age-related macular degeneration (AMD) in subjects living in the Lanusei valley, Central Sardinia, Italy, involved in a study on ageing (SardiNIA project). METHODS: A total of 814 volunteers aged ≥ 50 years, randomly selected from the SardiNIA project dataset, were included. A color fundus (CF) photograph of the 30° central retina of each eye was obtained and graded according to the Age-Related Eye Disease Study system. Life-style choices were investigated using standardized questionnaires. The concentrations of several inflammatory biomarkers (i.e., complement component, fibrinogen, and C-reactive protein) were measured. Polygenic risk score (PRS) was calculated and compared with results obtained from a European cohort. RESULTS: A total of 756 subjects had gradable CF photographs for AMD detection. In 91.3%, no signs of AMD were observed. The prevalence rates of early and late AMDs were 6.9% and 0.6%, respectively. A total of 85% of subjects were physically active; only 13.5% were current smokers. Low concentrations of complement component, fibrinogen, and C-reactive protein were found. We calculated the polygenic risk scores (PRS) using 40 AMD markers distributed on several candidate genes in Europeans and Sardinians. The mean PRS value was significantly lower in Sardinians than in the Europeans (0.21 vs. 0.248, respectively, p = 1.18 × 10-77). CONCLUSIONS: In our cohort, most subjects showed no sign of any AMD type and late AMD was a condition rarely observed. Results of genetic, biochemical, and life-style investigation support the hypothesis that Sardinia population may present of a peculiar background with a protective effect against AMD development.

Effect of Genetic Factors, Age and Sex on Levels of Circulating Extracellular Vesicles and Platelets.

Extracellular vesicles (EVs) mediate cell interactions in biological processes, such as receptor activation or molecule transfer. Estimates of variation by age and sex have been limited by small sample size, and no report has assessed the contribution of genetic factors to levels of EVs. Here, we evaluated blood levels of 25 EV and 3 platelet traits in 974 individuals (933 genotyped) and reported the first genome-wide association study (GWAS) on levels of these traits. EV levels all decreased with age, whereas the trend for their surface markers was more heterogeneous. Platelets and CD31dim platelet EVs significantly increased in females compared to males, although CD31 expression on both platelets and platelet EVs decreased in females. Levels of the other EV subsets were similar between sexes. GWAS revealed three statistically significant genetic signals associated with EV levels in the F10 and GBP1 genes and in the intergenic region between LRIG1 and KBTBD8. These add to a signal in the 3'UTR of RHOF associated with CD31 expression on platelets that was previously found to be associated with other platelet traits. These findings suggest that EV formation is not a simple, constant adjunct of metabolism but is under both age-related and genetic control that can be independent of the regulation of the levels of the cells from which the EVs derive.

Mitochondrial genetic variation is enriched in G-quadruplex regions that stall DNA synthesis in vitro.

As the powerhouses of the eukaryotic cell, mitochondria must maintain their genomes which encode proteins essential for energy production. Mitochondria are characterized by guanine-rich DNA sequences that spontaneously form unusual three-dimensional structures known as G-quadruplexes (G4). G4 structures can be problematic for the essential processes of DNA replication and transcription because they deter normal progression of the enzymatic-driven processes. In this study, we addressed the hypothesis that mitochondrial G4 is a source of mutagenesis leading to base-pair substitutions. Our computational analysis of 2757 individual genomes from two Italian population cohorts (SardiNIA and InCHIANTI) revealed a statistically significant enrichment of mitochondrial mutations within sequences corresponding to stable G4 DNA structures. Guided by the computational analysis results, we designed biochemical reconstitution experiments and demonstrated that DNA synthesis by two known mitochondrial DNA polymerases (Pol γ, PrimPol) in vitro was strongly blocked by representative stable G4 mitochondrial DNA structures, which could be overcome in a specific manner by the ATP-dependent G4-resolving helicase Pif1. However, error-prone DNA synthesis by PrimPol using the G4 template sequence persisted even in the presence of Pif1. Altogether, our results suggest that genetic variation is enriched in G-quadruplex regions that impede mitochondrial DNA replication.

Genetic Variation among Pharmacogenes in the Sardinian Population.

Pharmacogenetics (PGx) aims to identify the genetic factors that determine inter-individual differences in response to drug treatment maximizing efficacy while decreasing the risk of adverse events. Estimating the prevalence of PGx variants involved in drug response, is a critical preparatory step for large-scale implementation of a personalized medicine program in a target population. Here, we profiled pharmacogenetic variation in fourteen clinically relevant genes in a representative sample set of 1577 unrelated sequenced Sardinians, an ancient island population that accounts for genetic variation in Europe as a whole, and, at the same time is enriched in genetic variants that are very rare elsewhere. To this end, we used PGxPOP, a PGx allele caller based on the guidelines created by the Clinical Pharmacogenetics Implementation Consortium (CPIC), to identify the main phenotypes associated with the PGx alleles most represented in Sardinians. We estimated that 99.43% of Sardinian individuals might potentially respond atypically to at least one drug, that on average each individual is expected to have an abnormal response to about 17 drugs, and that for 27 drugs the fraction of the population at risk of atypical responses to therapy is more than 40%. Finally, we identified 174 pharmacogenetic variants for which the minor allele frequency was at least 10% higher among Sardinians as compared to other European populations, a fact that may contribute to substantial interpopulation variability in drug response phenotypes. This study provides baseline information for further large-scale pharmacogenomic investigations in the Sardinian population and underlines the importance of PGx characterization of diverse European populations, such as Sardinians.

Genetic-Driven Druggable Target Identification and Validation.

Choosing the right biological target is the critical primary decision for the development of new drugs. Systematic genetic association testing of both human diseases and quantitative traits, along with resultant findings of coincident associations between them, is becoming a powerful approach to infer drug targetable candidates and generate in vitro tests to identify compounds that can modulate them therapeutically. Here, we discuss opportunities and challenges, and infer criteria for the optimal use of genetic findings in the drug discovery pipeline.

PRF1 mutation alters immune system activation, inflammation, and risk of autoimmunity.

BACKGROUND: Defective alleles within the PRF1 gene, encoding the pore-forming protein perforin, in combination with environmental factors, cause familial type 2 hemophagocytic lymphohistiocytosis (FHL2), a rare, severe autosomal recessive childhood disorder characterized by massive release of cytokines-cytokine storm. OBJECTIVE: The aim of this study was to determine the function of hypomorph PRF1:p.A91V g.72360387 G > A on multiple sclerosis (MS) and type 1 diabetes (T1D). METHODS: We cross-compare the association data for PRF1:p.A91V mutation derived from GWAS on adult MS and pediatric T1D in Sardinians. The novel association with T1D was replicated in metanalysis in 12,584 cases and 17,692 controls from Sardinia, the United Kingdom, and Scotland. To dissect this mutation function, we searched through the coincident association immunophenotypes in additional set of general population Sardinians. RESULTS: We report that PRF1:p.A91V, is associated with increase of lymphocyte levels, especially within the cytotoxic memory T-cells, at general population level with reduced interleukin 7 receptor expression on these cells. The minor allele increased risk of MS, in 2903 cases and 2880 controls from Sardinia p = 2.06 × 10-4, odds ratio OR = 1.29, replicating a previous finding, whereas it protects from T1D p = 1.04 × 10-5, OR = 0.82. CONCLUSION: Our results indicate opposing contributions of the cytotoxic T-cell compartment to MS and T1D pathogenesis.

Eda-activated RelB recruits an SWI/SNF (BAF) chromatin-remodeling complex and initiates gene transcription in skin appendage formation.

Ectodysplasin A (Eda) signaling activates NF-κB during skin appendage formation, but how Eda controls specific gene transcription remains unclear. Here, we find that Eda triggers the formation of an NF-κB-associated SWI/SNF (BAF) complex in which p50/RelB recruits a linker protein, Tfg, that interacts with BAF45d in the BAF complex. We further reveal that Tfg is initially induced by Eda-mediated RelB activation and then bridges RelB and BAF for subsequent gene regulation. The BAF component BAF250a is particularly up-regulated in skin appendages, and epidermal knockout of BAF250a impairs skin appendage development, resulting in phenotypes similar to those of Eda-deficient mouse models. Transcription profiling identifies several target genes regulated by Eda, RelB, and BAF. Notably, RelB and the BAF complex are indispensable for transcription of Eda target genes, and both BAF complex and Eda signaling are required to open chromatin of Eda targets. Our studies thus suggest that Eda initiates a signaling cascade and recruits a BAF complex to specific gene loci to facilitate transcription during organogenesis.

Relative impact of indels versus SNPs on complex disease.

It is unclear whether insertions and deletions (indels) are more likely to influence complex traits than abundant single-nucleotide polymorphisms (SNPs). We sought to understand which category of variation is more likely to impact health. Using the SardiNIA study as an exemplar, we characterized 478,876 common indels and 8,246,244 common SNPs in up to 5,949 well-phenotyped individuals from an isolated valley in Sardinia. We assessed association between 120 traits, resulting in 89 nonoverlapping-associated loci.We evaluated whether indels were enriched among credible sets of potential causal variants. These credible sets included 1,319 SNPs and 88 indels. We did not find indels to be significantly enriched. Indels were the most likely causal variant in seven loci, including one locus associated with monocyte count where an indel with causality and mechanism previously demonstrated (rs200748895:TGCTG/T) had a 0.999 posterior probability. Overall, our results show a very modest and nonsignificant enrichment for common indels in associated loci.

Meta-MultiSKAT: Multiple phenotype meta-analysis for region-based association test.

The power of genetic association analyses can be increased by jointly meta-analyzing multiple correlated phenotypes. Here, we develop a meta-analysis framework, Meta-MultiSKAT, that uses summary statistics to test for association between multiple continuous phenotypes and variants in a region of interest. Our approach models the heterogeneity of effects between studies through a kernel matrix and performs a variance component test for association. Using a genotype kernel, our approach can test for rare-variants and the combined effects of both common and rare-variants. To achieve robust power, within Meta-MultiSKAT, we developed fast and accurate omnibus tests combining different models of genetic effects, functional genomic annotations, multiple correlated phenotypes, and heterogeneity across studies. In addition, Meta-MultiSKAT accommodates situations where studies do not share exactly the same set of phenotypes or have differing correlation patterns among the phenotypes. Simulation studies confirm that Meta-MultiSKAT can maintain the type-I error rate at the exome-wide level of 2.5 × 10-6 . Further simulations under different models of association show that Meta-MultiSKAT can improve the power of detection from 23% to 38% on average over single phenotype-based meta-analysis approaches. We demonstrate the utility and improved power of Meta-MultiSKAT in the meta-analyses of four white blood cell subtype traits from the Michigan Genomics Initiative (MGI) and SardiNIA studies.

Kidney size in relation to ageing, gender, renal function, birthweight and chronic kidney disease risk factors in a general population.

BACKGROUND: The relationship of kidney size to ageing, kidney function and kidney disease risk factors is not fully understood. METHODS: Ultrasound length and parenchymal kidney volume were determined from a population-based sample of 3972 Sardinians (age range 18-100 years). We then identified the subset of 2256 'healthy' subjects to define age- and sex-specific reference ranges (2.5-97.5 percentile) of kidney volume. Logistic regression (accounting for family clustering) was used to identify the clinical characteristics associated with abnormally large kidneys or abnormally small kidneys. RESULTS: In the healthy subset, kidney volume and length increased up to the fourth to fifth decade of life followed by a progressive decrease in men, whereas there was a gradual kidney volume decrease throughout the lifespan of women. In the whole sample, independent predictors of lower kidney volume (<2.5 percentile for age and sex) were male sex, low body mass index, short height, low waist:hip ratio and high serum creatinine (SCr); the independent predictors of larger kidney volume (>97.5 percentile for age and sex) were younger age, female sex, diabetes, obesity, high height, high waist:hip ratio and lower SCr. Estimated heritability for kidney volume was 15%, and for length 27%; kidney volume correlated strongly with birthweight. CONCLUSIONS: Overall, in a general healthy population, kidney measures declined with age differently in men and women. The determinants of kidney parenchymal volume include genetic factors and modifiable clinical factors.

Variation in human chromosome 21 ribosomal RNA genes characterized by TAR cloning and long-read sequencing.

Despite the key role of the human ribosome in protein biosynthesis, little is known about the extent of sequence variation in ribosomal DNA (rDNA) or its pre-rRNA and rRNA products. We recovered ribosomal DNA segments from a single human chromosome 21 using transformation-associated recombination (TAR) cloning in yeast. Accurate long-read sequencing of 13 isolates covering ∼0.82 Mb of the chromosome 21 rDNA complement revealed substantial variation among tandem repeat rDNA copies, several palindromic structures and potential errors in the previous reference sequence. These clones revealed 101 variant positions in the 45S transcription unit and 235 in the intergenic spacer sequence. Approximately 60% of the 45S variants were confirmed in independent whole-genome or RNA-seq data, with 47 of these further observed in mature 18S/28S rRNA sequences. TAR cloning and long-read sequencing enabled the accurate reconstruction of multiple rDNA units and a new, high-quality 44 838 bp rDNA reference sequence, which we have annotated with variants detected from chromosome 21 of a single individual. The large number of variants observed reveal heterogeneity in human rDNA, opening up the possibility of corresponding variations in ribosome dynamics.

Arterial stiffness and multiple organ damage: a longitudinal study in population.

AIMS: Previous cross-sectional observation identified arterial aging, indexed as pulse-wave velocity (PWV), as a key determinant of the simultaneous multiple organ damage (heart, carotid artery, and kidney). The aim of the present cohort study is to investigate trajectories of repeated measures of PWV and traditional CV risk factors in subjects who eventually presented clinical evidence of multiple organ damage in the SardiNIA study. METHODS AND RESULTS: Organ damage was measured in the heart (left ventricular hypertrophy, LVH), the common carotid artery (intima-media thickness > 0.9 mm and/or plaque), and the kidney (eGFR < 60 ml/min/1.73 m2) of 2130 men and women of a broad age range participating the SardiNIA study. SHATS was defined as the simultaneous occurrence of all the three-organ damages. Trajectory in traditional CV risk factors and PWV was analyzed retrospectively (four observations over 9 years) according to the number of organ damage (from 0 to 3). Compared to subjects with no organ damage, after controlling for traditional CV risk factors, each 1 m/s increase in baseline PWV was accompanied by a 93% higher odds of developing SHATS; and each 1 cm/s (0.01 m/s) annual increase in PWV by a 31% greater odds of developing SHATS. CONCLUSIONS: Arterial stiffness, a proxy of arterial aging that can be measured clinically as PWV, is an integrated predictive marker of multiple age-associated organ damage recognized as clinical diseases.

Molecular dynamics of Dkk4 modulates Wnt action and regulates meibomian gland development.

Secreted Dickkopf (Dkk) proteins are major Wnt pathway modulators during organ development. Dkk1 has been widely studied and acts as a general Wnt inhibitor. However, the molecular function of other Dkks remains largely unknown. Here, we show that Dkk4 selectively inhibits a subset of Wnts, but is further inactivated by proteolytic cleavage. Meibomian gland (MG) formation is employed as a model where Dkk4 and its Wnt targets are expressed. Skin-specific expression of Dkk4 arrests MG growth at early germ phase, which is similar to that observed in Eda-ablated Tabby mice. Consistent with transient Dkk4 action, intact Dkk4 inhibits MG extension but the cleaved form progressively increases during MG development with a concomitant upswing in Wnt activity. Furthermore, both Dkk4 and its receptor (and Wnt co-receptor) Lrp6 are direct Eda targets during MG induction. In cell and organotypic cultures, Dkk4 inhibition is eliminated by elevation of Lrp6. Also, Lrp6 upregulation restores MG formation in Tabby mice. Thus, the dynamic state of Dkk4 itself and its interaction with Lrp6 modulates Wnt function during MG development, with a novel limitation of Dkk4 action by proteolytic cleavage.