scRNA-sequencing in chick suggests a probabilistic model for cell fate allocation at the neural plate border

The vertebrate 'neural plate border' is a transient territory located at the edge of the neural plate containing precursors for all ectodermal derivatives: the neural plate, neural crest, placodes and epidermis. Elegant functional experiments in a range of vertebrate models have provided an in-depth understanding of gene regulatory interactions within the ectoderm. However, these experiments conducted at tissue level raise seemingly contradictory models for fate allocation of individual cells. Here, we carry out single cell RNA sequencing of chick ectoderm from primitive streak to neurulation stage, to explore cell state diversity and heterogeneity. We characterise the dynamics of gene modules, allowing us to model the order of molecular events which take place as ectodermal fates segregate. Furthermore, we find that genes previously classified as neural plate border 'specifiers' typically exhibit dynamic expression patterns and are enriched in either neural, neural crest or placodal fates, revealing that the neural plate border should be seen as a heterogeneous ectodermal territory and not a discrete transitional transcriptional state. Analysis of neural, neural crest and placodal markers reveals that individual NPB cells co-express competing transcriptional programmes suggesting that their ultimate identify is not yet fixed. This population of 'border located undecided progenitors' (BLUPs) gradually diminishes as cell fate decisions take place. Considering our findings, we propose a probabilistic model for cell fate choice at the neural plate border. Our data suggest that the probability of a progenitor's daughters to contribute to a given ectodermal derivative is related to the balance of competing transcriptional programmes, which in turn are regulated by the spatiotemporal position of a progenitor.

nf-core/clipseq - a robust Nextflow pipeline for comprehensive CLIP data analysis

Crosslinking and immunoprecipitation (CLIP) technologies have become a central component of the molecular biologists' toolkit to study protein-RNA interactions and thus to uncover core principles of RNA biology. There has been a proliferation of CLIP-based experimental protocols, as well as computational tools, especially for peak-calling. Consequently, there is an urgent need for a well-documented bioinformatic pipeline that enshrines the principles of robustness, reproducibility, scalability, portability and flexibility while embracing the diversity of experimental and computational CLIP tools. To address this, we present nf-core/clipseq - a robust Nextflow pipeline for quality control and analysis of CLIP sequencing data. It is part of the international nf-core community effort to develop and curate a best-practice, gold-standard set of pipelines for data analysis. The standards enabled by Nextflow and nf-core, including workflow management, version control, continuous integration and containerisation ensure that these key needs are met. Furthermore, multiple tools are implemented ( e.g. for peak-calling), alongside visualisation of quality control metrics to empower the user to make their own informed decisions based on their data. nf-core/clipseq remains under active development, with plans to incorporate newly released tools to ensure that pipeline remains up-to-date and relevant for the community. Engagement with users and developers is encouraged through the nf-core GitHub repository and Slack channel to promote collaboration. It is available at https://nf-co.re/clipseq.

A Genetic Risk Score Distinguishes Different Types of Autoantibody-Mediated Membranous Nephropathy

Introduction

Membranous nephropathy (MN) is the leading cause of nephrotic syndrome in adults and is characterized by detectable autoantibodies against glomerular antigens, most commonly phospholipase A2 receptor 1 (PLA2R1) and thrombospondin type-1 domain containing 7A (THSD7A). In Europeans, genetic variation in at least five loci, PLA2R1, HLA-DRB1, HLA-DQA1, IRF4, and NFKB1, affects the risk of disease. Here, we investigated the genetic risk differences between different autoantibody states.

Methods

1,409 MN individuals were genotyped genome-wide with a dense SNV array. The genetic risk score (GRS) was calculated utilizing the previously identified European MN loci, and results were compared with 4,929 healthy controls and 422 individuals with steroid-sensitive nephrotic syndrome.

Results

GRS was calculated in the 759 MN individuals in whom antibody status was known. The GRS for MN was elevated in the anti-PLA2R1 antibody-positive (N = 372) compared with both the unaffected control (N = 4,929) and anti-THSD7A-positive (N = 31) groups (p < 0.0001 for both comparisons), suggesting that this GRS reflects anti-PLA2R1 MN. Among PLA2R1-positive patients, GRS was inversely correlated with age of disease onset (p = 0.009). Further, the GRS in the dual antibody-negative group (N = 355) was intermediate between controls and the PLA2R1-positive group (p < 0.0001).

Conclusion

We demonstrate that the genetic risk factors for PLA2R1- and THSD7A-antibody-associated MN are different. A higher GRS is associated with younger age of onset of disease. Further, a proportion of antibody-negative MN cases have an elevated GRS similar to PLA2R1-positive disease. This suggests that in some individuals with negative serology the disease is driven by autoimmunity against PLA2R1.

Cell-intrinsic Aryl Hydrocarbon Receptor signalling is required for the resolution of injury-induced colonic stem cells

The aryl hydrocarbon receptor (AHR) is an environmental sensor that integrates microbial and dietary cues to influence physiological processes within the intestinal microenvironment, protecting against colitis and colitis-associated colorectal cancer development. Rapid tissue regeneration upon injury is important for the reinstatement of barrier integrity and its dysregulation promotes malignant transformation. Here we show that AHR is important for the termination of the regenerative response and the reacquisition of mature epithelial cell identity post injury in vivo and in organoid cultures in vitro. Using an integrative multi-omics approach in colon organoids, we show that AHR is required for timely termination of the regenerative response through direct regulation of transcription factors involved in epithelial cell differentiation as well as restriction of chromatin accessibility to regeneration-associated Yap/Tead transcriptional targets. Safeguarding a regulated regenerative response places AHR at a pivotal position in the delicate balance between controlled regeneration and malignant transformation.

Rapid intestinal regeneration after injury is critical to maintain barrier integrity and homeostasis, but must be tightly controlled to prevent tumorigenesis. Here they show that the aryl hydrocarbon receptor is required to terminate the regenerative response after wound healing.

Identification of a novel genetic locus associated with immune mediated thrombotic thrombocytopenic purpura

Immune thrombotic thrombocytopenic purpura (iTTP) is an ultra-rare, life-threatening disorder, mediated through severe ADAMTS13 deficiency causing multi-system micro-thrombi formation, and has specific human leukocyte antigen associations. We undertook a large genome-wide association study to investigate additional genetically distinct associations in iTTP. We compared two iTTP patient cohorts with controls, following standardized genome-wide quality control procedures for single-nucleotide polymorphisms and imputed HLA types. Associations were functionally investigated using expression quantitative trait loci (eQTL), and motif binding prediction software. Independent associations consistent with previous findings in iTTP were detected at the HLA locus and in addition a novel association was detected on chromosome 3 (rs9884090, P=5.22x10^-10, odds ratio 0.40) in the UK discovery cohort. Meta-analysis, including the French replication cohort, strengthened the associations. The haploblock containing rs9884090 is associated with reduced protein O-glycosyltransferase 1 (POGLUT1) expression (eQTL P<0.05), and functional annotation suggested a potential causative variant (rs71767581). This work implicates POGLUT1 in iTTP pathophysiology and suggests altered post-translational modification of its targets may influence disease susceptibility.

Genetic Identification of Two Novel Loci Associated with Steroid-Sensitive Nephrotic Syndrome

BACKGROUND

Steroid-sensitive nephrotic syndrome (SSNS), the most common form of nephrotic syndrome in childhood, is considered an autoimmune disease with an established classic HLA association. However, the precise etiology of the disease is unclear. In other autoimmune diseases, the identification of loci outside the classic HLA region by genome-wide association studies (GWAS) has provided critical insights into disease pathogenesis. Previously conducted GWAS of SSNS have not identified non-HLA loci achieving genome-wide significance.

METHODS

In an attempt to identify additional loci associated with SSNS, we conducted a GWAS of a large cohort of European ancestry comprising 422 ethnically homogeneous pediatric patients and 5642 ethnically matched controls.

RESULTS

The GWAS found three loci that achieved genome-wide significance, which explain approximately 14% of the genetic risk for SSNS. It confirmed the previously reported association with the HLA-DR/DQ region (lead single-nucleotide polymorphism [SNP] rs9273542, P=1.59×10^-43; odds ratio [OR], 3.39; 95% confidence interval [95% CI], 2.86 to 4.03) and identified two additional loci outside the HLA region on chromosomes 4q13.3 and 6q22.1. The latter contains the calcium homeostasis modulator family member 6 gene CALHM6 (previously called FAM26F). CALHM6 is implicated in immune response modulation; the lead SNP (rs2637678, P=1.27×10^-17; OR, 0.51; 95% CI, 0.44 to 0.60) exhibits strong expression quantitative trait loci effects, the risk allele being associated with lower lymphocytic expression of CALHM6.

CONCLUSIONS

Because CALHM6 is implicated in regulating the immune response to infection, this may provide an explanation for the typical triggering of SSNS onset by infections. Our results suggest that a genetically conferred risk of immune dysregulation may be a key component in the pathogenesis of SSNS.

Polycystic Kidney Disease with Hyperinsulinemic Hypoglycemia Caused by a Promoter Mutation in Phosphomannomutase 2

Hyperinsulinemic hypoglycemia (HI) and congenital polycystic kidney disease (PKD) are rare, genetically heterogeneous disorders. The co-occurrence of these disorders (HIPKD) in 17 children from 11 unrelated families suggested an unrecognized genetic disorder. Whole-genome linkage analysis in five informative families identified a single significant locus on chromosome 16p13.2 (logarithm of odds score 6.5). Sequencing of the coding regions of all linked genes failed to identify biallelic mutations. Instead, we found in all patients a promoter mutation (c.-167G>T) in the phosphomannomutase 2 gene (PMM2), either homozygous or in trans with PMM2 coding mutations. PMM2 encodes a key enzyme in N-glycosylation. Abnormal glycosylation has been associated with PKD, and we found that deglycosylation in cultured pancreatic β cells altered insulin secretion. Recessive coding mutations in PMM2 cause congenital disorder of glycosylation type 1a (CDG1A), a devastating multisystem disorder with prominent neurologic involvement. Yet our patients did not exhibit the typical clinical or diagnostic features of CDG1A. In vitro, the PMM2 promoter mutation associated with decreased transcriptional activity in patient kidney cells and impaired binding of the transcription factor ZNF143. In silico analysis suggested an important role of ZNF143 for the formation of a chromatin loop including PMM2 We propose that the PMM2 promoter mutation alters tissue-specific chromatin loop formation, with consequent organ-specific deficiency of PMM2 leading to the restricted phenotype of HIPKD. Our findings extend the spectrum of genetic causes for both HI and PKD and provide insights into gene regulation and PMM2 pleiotropy.