Disruptions of topological chromatin domains cause pathogenic rewiring of gene-enhancer interactions.

Mammalian genomes are organized into megabase-scale topologically associated domains (TADs). We demonstrate that disruption of TADs can rewire long-range regulatory architecture and result in pathogenic phenotypes. We show that distinct human limb malformations are caused by deletions, inversions, or duplications altering the structure of the TAD-spanning WNT6/IHH/EPHA4/PAX3 locus. Using CRISPR/Cas genome editing, we generated mice with corresponding rearrangements. Both in mouse limb tissue and patient-derived fibroblasts, disease-relevant structural changes cause ectopic interactions between promoters and non-coding DNA, and a cluster of limb enhancers normally associated with Epha4 is misplaced relative to TAD boundaries and drives ectopic limb expression of another gene in the locus. This rewiring occurred only if the variant disrupted a CTCF-associated boundary domain. Our results demonstrate the functional importance of TADs for orchestrating gene expression via genome architecture and indicate criteria for predicting the pathogenicity of human structural variants, particularly in non-coding regions of the human genome.

Duplications disrupt chromatin architecture and rewire GPR101-enhancer communication in X-linked acrogigantism.

X-linked acrogigantism (X-LAG) is the most severe form of pituitary gigantism and is characterized by aggressive growth hormone (GH)-secreting pituitary tumors that occur in early childhood. X-LAG is associated with chromosome Xq26.3 duplications (the X-LAG locus typically includes VGLL1, CD40LG, ARHGEF6, RBMX, and GPR101) that lead to massive pituitary tumoral expression of GPR101, a novel regulator of GH secretion. The mechanism by which the duplications lead to marked pituitary misexpression of GPR101 alone was previously unclear. Using Hi-C and 4C-seq, we characterized the normal chromatin structure at the X-LAG locus. We showed that GPR101 is located within a topologically associating domain (TAD) delineated by a tissue-invariant border that separates it from centromeric genes and regulatory sequences. Next, using 4C-seq with GPR101, RBMX, and VGLL1 viewpoints, we showed that the duplications in multiple X-LAG-affected individuals led to ectopic interactions that crossed the invariant TAD border, indicating the existence of a similar and consistent mechanism of neo-TAD formation in X-LAG. We then identified several pituitary active cis-regulatory elements (CREs) within the neo-TAD and demonstrated in vitro that one of them significantly enhanced reporter gene expression. At the same time, we showed that the GPR101 promoter permits the incorporation of new regulatory information. Our results indicate that X-LAG is a TADopathy of the endocrine system in which Xq26.3 duplications disrupt the local chromatin architecture forming a neo-TAD. Rewiring GPR101-enhancer interaction within the new regulatory unit is likely to cause the high levels of aberrant expression of GPR101 in pituitary tumors caused by X-LAG.

Gain of gene regulatory network interconnectivity at the origin of vertebrates.

SignificanceIn this manuscript, we address an essential question in developmental and evolutionary biology: How have changes in gene regulatory networks contributed to the invertebrate-to-vertebrate transition? To address this issue, we perturbed four signaling pathways critical for body plan formation in the cephalochordate amphioxus and in zebrafish and compared the effects of such perturbations on gene expression and gene regulation in both species. Our data reveal that many developmental genes have gained response to these signaling pathways in the vertebrate lineage. Moreover, we show that the interconnectivity between these pathways is much higher in zebrafish than in amphioxus. We conclude that this increased signaling pathway complexity likely contributed to vertebrate morphological novelties during evolution.

Position effects at the FGF8 locus are associated with femoral hypoplasia.

Copy-number variations (CNVs) are a common cause of congenital limb malformations and are interpreted primarily on the basis of their effect on gene dosage. However, recent studies show that CNVs also influence the 3D genome chromatin organization. The functional interpretation of whether a phenotype is the result of gene dosage or a regulatory position effect remains challenging. Here, we report on two unrelated families with individuals affected by bilateral hypoplasia of the femoral bones, both harboring de novo duplications on chromosome 10q24.32. The ∼0.5 Mb duplications include FGF8, a key regulator of limb development and several limb enhancer elements. To functionally characterize these variants, we analyzed the local chromatin architecture in the affected individuals' cells and re-engineered the duplications in mice by using CRISPR-Cas9 genome editing. We found that the duplications were associated with ectopic chromatin contacts and increased FGF8 expression. Transgenic mice carrying the heterozygous tandem duplication including Fgf8 exhibited proximal shortening of the limbs, resembling the human phenotype. To evaluate whether the phenotype was a result of gene dosage, we generated another transgenic mice line, carrying the duplication on one allele and a concurrent Fgf8 deletion on the other allele, as a control. Surprisingly, the same malformations were observed. Capture Hi-C experiments revealed ectopic interaction with the duplicated region and Fgf8, indicating a position effect. In summary, we show that duplications at the FGF8 locus are associated with femoral hypoplasia and that the phenotype is most likely the result of position effects altering FGF8 expression rather than gene dosage effects.

Formation of new chromatin domains determines pathogenicity of genomic duplications.

Chromosome conformation capture methods have identified subchromosomal structures of higher-order chromatin interactions called topologically associated domains (TADs) that are separated from each other by boundary regions. By subdividing the genome into discrete regulatory units, TADs restrict the contacts that enhancers establish with their target genes. However, the mechanisms that underlie partitioning of the genome into TADs remain poorly understood. Here we show by chromosome conformation capture (capture Hi-C and 4C-seq methods) that genomic duplications in patient cells and genetically modified mice can result in the formation of new chromatin domains (neo-TADs) and that this process determines their molecular pathology. Duplications of non-coding DNA within the mouse Sox9 TAD (intra-TAD) that cause female to male sex reversal in humans, showed increased contact of the duplicated regions within the TAD, but no change in the overall TAD structure. In contrast, overlapping duplications that extended over the next boundary into the neighbouring TAD (inter-TAD), resulted in the formation of a new chromatin domain (neo-TAD) that was isolated from the rest of the genome. As a consequence of this insulation, inter-TAD duplications had no phenotypic effect. However, incorporation of the next flanking gene, Kcnj2, in the neo-TAD resulted in ectopic contacts of Kcnj2 with the duplicated part of the Sox9 regulatory region, consecutive misexpression of Kcnj2, and a limb malformation phenotype. Our findings provide evidence that TADs are genomic regulatory units with a high degree of internal stability that can be sculptured by structural genomic variations. This process is important for the interpretation of copy number variations, as these variations are routinely detected in diagnostic tests for genetic disease and cancer. This finding also has relevance in an evolutionary setting because copy-number differences are thought to have a crucial role in the evolution of genome complexity.

Optimization of the Mainzelliste software for fast privacy-preserving record linkage.

BACKGROUND: Data analysis for biomedical research often requires a record linkage step to identify records from multiple data sources referring to the same person. Due to the lack of unique personal identifiers across these sources, record linkage relies on the similarity of personal data such as first and last names or birth dates. However, the exchange of such identifying data with a third party, as is the case in record linkage, is generally subject to strict privacy requirements. This problem is addressed by privacy-preserving record linkage (PPRL) and pseudonymization services. Mainzelliste is an open-source record linkage and pseudonymization service used to carry out PPRL processes in real-world use cases. METHODS: We evaluate the linkage quality and performance of the linkage process using several real and near-real datasets with different properties w.r.t. size and error-rate of matching records. We conduct a comparison between (plaintext) record linkage and PPRL based on encoded records (Bloom filters). Furthermore, since the Mainzelliste software offers no blocking mechanism, we extend it by phonetic blocking as well as novel blocking schemes based on locality-sensitive hashing (LSH) to improve runtime for both standard and privacy-preserving record linkage. RESULTS: The Mainzelliste achieves high linkage quality for PPRL using field-level Bloom filters due to the use of an error-tolerant matching algorithm that can handle variances in names, in particular missing or transposed name compounds. However, due to the absence of blocking, the runtimes are unacceptable for real use cases with larger datasets. The newly implemented blocking approaches improve runtimes by orders of magnitude while retaining high linkage quality. CONCLUSION: We conduct the first comprehensive evaluation of the record linkage facilities of the Mainzelliste software and extend it with blocking methods to improve its runtime. We observed a very high linkage quality for both plaintext as well as encoded data even in the presence of errors. The provided blocking methods provide order of magnitude improvements regarding runtime performance thus facilitating the use in research projects with large datasets and many participants.

MiR-21 is required for efficient kidney regeneration in fish.

BACKGROUND: Acute kidney injury in mammals, which is caused by cardiovascular diseases or the administration of antibiotics with nephrotoxic side-effects is a life-threatening disease, since loss of nephrons is irreversible in mammals. In contrast, fish are able to generate new nephrons even in adulthood and thus provide a good model to study renal tubular regeneration. RESULTS: Here, we investigated the early response after gentamicin-induced renal injury, using the short-lived killifish Nothobranchius furzeri. A set of microRNAs was differentially expressed after renal damage, among them miR-21, which was up-regulated. A locked nucleic acid-modified antimiR-21 efficiently knocked down miR-21 activity and caused a lag in the proliferative response, enhanced apoptosis and an overall delay in regeneration. Transcriptome profiling identified apoptosis as a process that was significantly affected upon antimiR-21 administration. Together with functional data this suggests that miR-21 acts as a pro-proliferative and anti-apoptotic factor in the context of kidney regeneration in fish. Possible downstream candidate genes that mediate its effect on proliferation and apoptosis include igfbp3 and fosl1, among other genes. CONCLUSION: In summary, our findings extend the role of miR-21 in the kidney. For the first time we show its functional involvement in regeneration indicating that fast proliferation and reduced apoptosis are important for efficient renal tubular regeneration.

Homeotic arm-to-leg transformation associated with genomic rearrangements at the PITX1 locus.

The study of homeotic-transformation mutants in model organisms such as Drosophila revolutionized the field of developmental biology, but how these mutants relate to human developmental defects remains to be elucidated. Here, we show that Liebenberg syndrome, an autosomal-dominant upper-limb malformation, shows features of a homeotic limb transformation in which the arms have acquired morphological characteristics of a leg. Using high-resolution array comparative genomic hybridization and paired-end whole-genome sequencing, we identified two deletions and a translocation 5' of PITX1. The structural changes are likely to remove active PITX1 forelimb suppressor and/or insulator elements and thereby move active enhancer elements in the vicinity of the PITX1 regulatory landscape. We generated transgenic mice in which PITX1 was misexpressed under the control of a nearby enhancer and were able to recapitulate the Liebenberg phenotype.

Case report: Management of pediatric gigantism caused by the TADopathy, X-linked acrogigantism.

X-linked acrogigantism (X-LAG) is a rare form of pituitary gigantism that is associated with growth hormone (GH) and prolactin-secreting pituitary adenomas/pituitary neuroendocrine tumors (PitNETs) that develop in infancy. It is caused by a duplication on chromosome Xq26.3 that leads to the misexpression of the gene GPR101, a constitutively active stimulator of pituitary GH and prolactin secretion. GPR101 normally exists within its own topologically associating domain (TAD) and is insulated from surrounding regulatory elements. X-LAG is a TADopathy in which the duplication disrupts a conserved TAD border, leading to a neo-TAD in which ectopic enhancers drive GPR101 over-expression, thus causing gigantism. Here we trace the full diagnostic and therapeutic pathway of a female patient with X-LAG from 4C-seq studies demonstrating the neo-TAD through medical and surgical interventions and detailed tumor histopathology. The complex nature of treating young children with X-LAG is illustrated, including the achievement of hormonal control using a combination of neurosurgery and adult doses of first-generation somatostatin analogs.

Chromatin conformation capture in the clinic: 4C-seq/HiC distinguishes pathogenic from neutral duplications at the GPR101 locus.

BACKGROUND: X-linked acrogigantism (X-LAG; MIM: 300942) is a severe form of pituitary gigantism caused by chromosome Xq26.3 duplications involving GPR101. X-LAG-associated duplications disrupt the integrity of the topologically associating domain (TAD) containing GPR101 and lead to the formation of a neo-TAD that drives pituitary GPR101 misexpression and gigantism. As X-LAG is fully penetrant and heritable, duplications involving GPR101 identified on prenatal screening studies, like amniocentesis, can pose an interpretation challenge for medical geneticists and raise important concerns for patients and families. Therefore, providing robust information on the functional genomic impact of such duplications has important research and clinical value with respect to gene regulation and triplosensitivity traits. METHODS: We employed 4C/HiC-seq as a clinical tool to determine the functional impact of incidentally discovered GPR101 duplications on TAD integrity in three families. After defining duplications and breakpoints around GPR101 by clinical-grade and high-density aCGH, we constructed 4C/HiC chromatin contact maps for our study population and compared them with normal and active (X-LAG) controls. RESULTS: We showed that duplications involving GPR101 that preserved the centromeric invariant TAD boundary did not generate a pathogenic neo-TAD and that ectopic enhancers were not adopted. This allowed us to discount presumptive/suspected X-LAG diagnoses and GPR101 misexpression, obviating the need for intensive clinical follow-up. CONCLUSIONS: This study highlights the importance of TAD boundaries and chromatin interactions in determining the functional impact of copy number variants and provides proof-of-concept for using 4C/HiC-seq as a clinical tool to acquire crucial information for genetic counseling and to support clinical decision-making in cases of suspected TADopathies.

In Situ Inclusion Detection and Material Characterization in an Electron Beam Powder Bed Fusion Process Using Electron Optical Imaging.

Electron Beam Powder Bed Fusion (PBF-EB) is an Additive Manufacturing (AM) method that utilizes an electron beam to melt and consolidate metal powder. The beam, combined with a backscattered electron detector, enables advanced process monitoring, a method termed Electron Optical Imaging (ELO). ELO is already known to provide great topographical information, but its capabilities regarding material contrast are less studied. In this article the extents of material contrast using ELO are investigated, focusing mainly on identifying powder contamination. It will be shown that an ELO detector is capable of distinguishing a single 100 µm foreign powder particle, during an PBF-EB process, if the backscattering coefficient of the inclusion is sufficiently higher than its surroundings. Additionally, it is investigated how the material contrast can be used for material characterization. A mathematical framework is provided to describe the relationship between the signal intensity in the detector and the effective atomic number Zeff of the imaged alloy. The approach is verified with empirical data from twelve different materials, demonstrating that the effective atomic number of an alloy can be predicted to within one atomic number from its ELO intensity.

Combinatorial effects on gene expression at the Lbx1/Fgf8 locus resolve split-hand/foot malformation type 3.

Split-Hand/Foot Malformation type 3 (SHFM3) is a congenital limb malformation associated with tandem duplications at the LBX1/FGF8 locus. Yet, the disease patho-mechanism remains unsolved. Here we investigate the functional consequences of SHFM3-associated rearrangements on chromatin conformation and gene expression in vivo in transgenic mice. We show that the Lbx1/Fgf8 locus consists of two separate, but interacting, regulatory domains. Re-engineering of a SHFM3-associated duplication and a newly reported inversion in mice results in restructuring of the chromatin architecture. This leads to ectopic activation of the Lbx1 and Btrc genes in the apical ectodermal ridge (AER) in an Fgf8-like pattern induced by AER-specific enhancers of Fgf8. We provide evidence that the SHFM3 phenotype is the result of a combinatorial effect on gene misexpression in the developing limb. Our results reveal insights into the molecular mechanism underlying SHFM3 and provide conceptual framework for how genomic rearrangements can cause gene misexpression and disease.

Evaluation of the Multiple Assessment Subjective Shoulder Score (MASSS). A pilot study.

INTRODUCTION: This study compares an entirely patient reported modification of the Constant-Murley Score (CS) (the Multiple Assessment Subjective Shoulder Score [MASSS]), with the CS in symptomatic patients with non-instability shoulder pathology. HYPOTHESIS: The hypothesis is that the sub scores and overall score of the MASSS will be well correlated and in high agreement to the corresponding scores of the CS. MATERIALS AND METHODS: The MASSS, the CS and adjusted-CS (a-CS) were administered to 74 patients. The MASSS replaced the strength and ROM domains of the CS with SANE scores. Correlations, differences, agreement, test-retest reliability, and internal consistency were assessed using Pearson's R test, paired t-test, Bland-Altman method, intraclass correlations (ICC) and Cronbach's alpha. RESULTS: Correlation between MASSS and CS (0.834) and a-CS (0.824) were excellent. Total MASSS (54.0) was higher than CS (41.8) (p<0.001) and a-CS (47.8) (p<0.05). MASSS strength and ROM domains were higher than those of the CS (Strength: 13.4 v 5.3 [p<0.001], ROM: 24.8 v 20.8 [p<0.001]. There was a floor effect for the strength domain of the CS (46.7% minimum score), but not the MASSS (0% minimum score). MASSS internal consistency (0.626-0.734) was better than the CS (0.401-0.643). Test re-test reliability of the total MASSS (ICC 0.93) and the strength (ICC 0.90) and ROM (ICC 0.86) domains was excellent. DISCUSSION: The MASSS has several advantages over the CS. Although the correlation of the total MASSS with the CS and a-CS was high, there was poor agreement in the strength and the ROM domains and therefore the MASSS value is not interchangeable with the CS, but rather should be considered as an alternative. LEVEL OF EVIDENCE: II; Prospective cohort study.

Cannabis Use for Endometriosis: Clinical and Legal Challenges in Australia and New Zealand.

Introduction: Endometriosis is a difficult to manage condition associated with a significant disease burden. High levels of illicit cannabis use for therapeutic purposes have been previously reported by endometriosis patients in Australia and New Zealand (NZ). Although access to legal medicinal cannabis (MC) is available through medical prescription via multiple federal schemes, significant barriers to patient access remain. Methods: An anonymous cross-sectional online survey was developed and distributed through social media via endometriosis advocacy groups worldwide. Respondents were asked about legal versus illicit cannabis usage, their understanding of access pathways and legal status, and their interactions with health care professionals. Results: Of 237 respondents who reported cannabis use with a medical diagnosis of endometriosis, 186 (72.0%) of Australian and 51 (88.2%) NZ respondents reported self-administering cannabis illicitly. Only 23.1% of Australian and 5.9% of NZ respondents accessed cannabis through a doctor's prescription, with 4.8% of Australian and no NZ respondents reporting to legally self-administer cannabis. Substantial substitution effects (>50% reduction) were observed in users of nonopioid analgesia (63.1%), opioid analgesia (66.1%), hormonal therapies (27.5%), antineuropathics (61.7%), antidepressants (28.2%) and antianxiety medications (47.9%). Of Australian respondents, 18.8% and of NZ respondents, 23.5% reported not disclosing their cannabis use to their medical doctor, citing concern over legal repercussions, societal judgment, or their doctors' reaction and presumed unwillingness to prescribe legal MC. Conclusions: Respondents self-reported positive outcomes when using cannabis for management of endometriosis, demonstrating a therapeutic potential for MC. Despite this, many are using cannabis without medical supervision. While evidence for a substantial substitution effect by cannabis was demonstrated in these data, of particular concern are the clinical consequences of using cannabis without medical supervision, particularly with regard to drug interactions and the tapering or cessation of certain medications without that supervision. Improving doctor and patient communication about MC use may improve levels of medical oversight, the preference for legal MC adoption over acquisition via illicit supply and reducing cannabis-associated stigma.

Regulatory Architecture of the RCA Gene Cluster Captures an Intragenic TAD Boundary, CTCF-Mediated Chromatin Looping and a Long-Range Intergenic Enhancer.

The Regulators of Complement Activation (RCA) gene cluster comprises several tandemly arranged genes with shared functions within the immune system. RCA members, such as complement receptor 2 (CR2), are well-established susceptibility genes in complex autoimmune diseases. Altered expression of RCA genes has been demonstrated at both the functional and genetic level, but the mechanisms underlying their regulation are not fully characterised. We aimed to investigate the structural organisation of the RCA gene cluster to identify key regulatory elements that influence the expression of CR2 and other genes in this immunomodulatory region. Using 4C, we captured extensive CTCF-mediated chromatin looping across the RCA gene cluster in B cells and showed these were organised into two topologically associated domains (TADs). Interestingly, an inter-TAD boundary was located within the CR1 gene at a well-characterised segmental duplication. Additionally, we mapped numerous gene-gene and gene-enhancer interactions across the region, revealing extensive co-regulation. Importantly, we identified an intergenic enhancer and functionally demonstrated this element upregulates two RCA members (CR2 and CD55) in B cells. We have uncovered novel, long-range mechanisms whereby autoimmune disease susceptibility may be influenced by genetic variants, thus highlighting the important contribution of chromatin topology to gene regulation and complex genetic disease.

Transduction of peptides and proteins into live cells by cell penetrating peptides.

Internalization of peptides and proteins into live cells is an essential prerequisite for studies on intracellular signal pathways, for treatment of certain microbial diseases and for signal transduction therapy, especially for cancer treatment. Cell penetrating peptides (CPPs) facilitate the transport of cargo-proteins through the cell membrane into live cells. CPPs which allow formation of non-covalent complexes with the cargo are used primarily in this study due to the relatively easy handling procedure. Efficiency of the protein uptake is estimated qualitatively by fluorescence microscopy and quantitatively by SDS-PAGE. Using the CPP cocktail JBS-Proteoducin, the intracellular concentrations of a secondary antibody and bovine serum albumin can reach the micromolar range. Internalization of antibodies allows mediation of intracellular pathways including knock down of signal transduction. The high specificity and affinity of antibodies makes them potentially more powerful than siRNA. Thus, CPPs represent a significant new possibility to study signal transduction processes in competition or in comparison to the commonly used other techniques. To estimate the highest attainable intracellular concentrations of cargo proteins, the CPPs are tested for cytotoxicity. Cell viability and membrane integrity relative to concentration of CPPs are investigated. Viability as estimated by the reductive activity of mitochondria (MTT-test) is more sensitive to higher concentrations of CPPs versus membrane integrity, as measured by the release of dead cell protease. Distinct differences in uptake efficiency and cytotoxic effects are found using six different CPPs and six different adhesion and suspension cell lines.

CTCF knockout in zebrafish induces alterations in regulatory landscapes and developmental gene expression.

Coordinated chromatin interactions between enhancers and promoters are critical for gene regulation. The architectural protein CTCF mediates chromatin looping and is enriched at the boundaries of topologically associating domains (TADs), which are sub-megabase chromatin structures. In vitro CTCF depletion leads to a loss of TADs but has only limited effects over gene expression, challenging the concept that CTCF-mediated chromatin structures are a fundamental requirement for gene regulation. However, how CTCF and a perturbed chromatin structure impacts gene expression during development remains poorly understood. Here we link the loss of CTCF and gene regulation during patterning and organogenesis in a ctcf knockout zebrafish model. CTCF absence leads to loss of chromatin structure and affects the expression of thousands of genes, including many developmental regulators. Our results demonstrate the essential role of CTCF in providing the structural context for enhancer-promoter interactions, thus regulating developmental genes.

Functional dissection of the Sox9-Kcnj2 locus identifies nonessential and instructive roles of TAD architecture.

The genome is organized in three-dimensional units called topologically associating domains (TADs), through a process dependent on the cooperative action of cohesin and the DNA-binding factor CTCF. Genomic rearrangements of TADs have been shown to cause gene misexpression and disease, but genome-wide depletion of CTCF has no drastic effects on transcription. Here, we investigate TAD function in vivo in mouse limb buds at the Sox9-Kcnj2 locus. We show that the removal of all major CTCF sites at the boundary and within the TAD resulted in a fusion of neighboring TADs, without major effects on gene expression. Gene misexpression and disease phenotypes, however, were achieved by redirecting regulatory activity through inversions and/or the repositioning of boundaries. Thus, TAD structures provide robustness and precision but are not essential for developmental gene regulation. Aberrant disease-related gene activation is not induced by a mere loss of insulation but requires CTCF-dependent redirection of enhancer-promoter contacts.

An evolutionary perspective of regulatory landscape dynamics in development and disease.

The organization of animal genomes into topologically associating domains (TADs) provides a structural scaffold in which cis-regulatory elements (CREs) operate on their target genes. Determining the position of CREs and genes relative to TADs has become instrumental to trace gene expression changes during evolution and in diseases. Here we will review recent studies and discuss TADs as structural units with respect to their conservation and stability during genome reorganization. Furthermore, we describe how TAD restructuring contributed to morphological novelties during evolution but also their deleterious effects associated with disease. Despite considering TADs as structural units, the nested and dynamic scaffold within TADs contributes to tissue-specific gene expression, implying that such changes can also account for gene expression differences during evolution.