guidedNOMe-seq quantifies chromatin states at single allele resolution for hundreds of custom regions in parallel.

Since the introduction of next generation sequencing technologies, the field of epigenomics has evolved rapidly. However, most commonly used assays are enrichment-based methods and thus only semi-quantitative. Nucleosome occupancy and methylome sequencing (NOMe-seq) allows for quantitative inference of chromatin states with single locus resolution, but this requires high sequencing depth and is therefore prohibitively expensive to routinely apply to organisms with large genomes. To overcome this limitation, we introduce guidedNOMe-seq, where we combine NOMe profiling with large scale sgRNA synthesis and Cas9-mediated region-of-interest (ROI) liberation. To facilitate quantitative comparisons between multiple samples, we additionally develop an R package to standardize differential analysis of any type of NOMe-seq data. We extensively benchmark guidedNOMe-seq in a proof-of-concept study, dissecting the interplay of ChAHP and CTCF on chromatin. In summary we present a cost-effective, scalable, and customizable target enrichment extension to the existing NOMe-seq protocol allowing genome-scale quantification of nucleosome occupancy and transcription factor binding at single allele resolution.

The ChAHP Complex Counteracts Chromatin Looping at CTCF Sites that Emerged from SINE Expansions in Mouse.

CCCTC-binding factor (CTCF) and cohesin are key players in three-dimensional chromatin organization. The topologically associating domains (TADs) demarcated by CTCF are remarkably well conserved between species, although genome-wide CTCF binding has diverged substantially following transposon-mediated motif expansions. Therefore, the CTCF consensus motif poorly predicts TADs, and additional factors must modulate CTCF binding and subsequent TAD formation. Here, we demonstrate that the ChAHP complex (CHD4, ADNP, HP1) competes with CTCF for a common set of binding motifs. In Adnp knockout cells, novel insulated regions are formed at sites normally bound by ChAHP, whereas proximal canonical boundaries are weakened. These data reveal that CTCF-mediated loop formation is modulated by a distinct zinc-finger protein complex. Strikingly, ChAHP-bound loci are mainly situated within less diverged SINE B2 transposable elements. This implicates ChAHP in maintenance of evolutionarily conserved spatial chromatin organization by buffering novel CTCF binding sites that emerged through SINE expansions.

Myocardial Fibrosis in Classical Low-Flow, Low-Gradient Aortic Stenosis.

Background Few data exist on the degree of interstitial myocardial fibrosis in patients with classical low-flow, low-gradient aortic stenosis (LFLG-AS) and its association with left ventricular flow reserve (FR) on dobutamine stress echocardiography. This study sought to evaluate the diffuse interstitial fibrosis measured by T1 mapping cardiac magnetic resonance technique in LFLG-AS patients with and without FR. Methods Prospective study including 65 consecutive patients (41 LFLG-AS [mean age, 67.1±8.4 years; 83% men] and 24 high-gradient aortic stenosis used as controls) undergoing dobutamine stress echocardiography to assess FR and cardiac magnetic resonance to determine the extracellular volume (ECV) fraction of the myocardium, indexed ECV (iECV) to body surface area and late gadolinium enhancement. Results Interstitial myocardial fibrosis measured by iECV was higher in patients with LFLG-AS with and without FR as compared with high-gradient aortic stenosis (35.25±9.75 versus 32.93±11.00 versus 21.19±6.47 mL/m2, respectively; P<0.001). However, both ECV and iECV levels were similar between LFLG-AS patients with and without FR ( P=0.950 and P=0.701, respectively). Also, FR did not correlate significantly with ECV (r=-0.16, P=0.31) or iECV (r=0.11, P=0.51). Late gadolinium enhancement mass was also similar in patients with versus without FR but lower in high-gradient aortic stenosis (13.3±10.2 versus 10.5±7.5 versus 4.8±5.9 g, respectively; P=0.018). Conclusions Patients with LFLG-AS have higher ECV, iECV, and late gadolinium enhancement mass compared with high-gradient aortic stenosis. Moreover, among patients with LFLG-AS, the degree of myocardial fibrosis was similar in patients with versus those without FR. These findings suggest that diffuse myocardial fibrosis may not be the main factor responsible for the absence of FR in LFLG-AS patients.

Tdrd6a Regulates the Aggregation of Buc into Functional Subcellular Compartments that Drive Germ Cell Specification.

Phase separation represents an important form of subcellular compartmentalization. However, relatively little is known about how the formation or disassembly of such compartments is regulated. In zebrafish, the Balbiani body (Bb) and the germ plasm (Gp) are intimately linked phase-separated structures essential for germ cell specification and home to many germ cell-specific mRNAs and proteins. Throughout development, these structures occur as a single large aggregate (Bb), which disperses throughout oogenesis and upon fertilization accumulates again into relatively large assemblies (Gp). Formation of the Bb requires Bucky ball (Buc), a protein with prion-like properties. We found that the multi-tudor domain-containing protein Tdrd6a interacts with Buc, affecting its mobility and aggregation properties. Importantly, lack of this regulatory interaction leads to significant defects in germ cell development. Our work presents insights into how prion-like protein aggregations can be regulated and highlights the biological relevance of such regulatory events.

Systemic Loss and Gain of Chromatin Architecture throughout Zebrafish Development.

The spatial organization of chromosomes is critical in establishing gene expression programs. We generated in situ Hi-C maps throughout zebrafish development to gain insight into higher-order chromatin organization and dynamics. Zebrafish chromosomes segregate in active and inactive chromatin (A/B compartments), which are further organized into topologically associating domains (TADs). Zebrafish A/B compartments and TADs have genomic features similar to those of their mammalian counterparts, including evolutionary conservation and enrichment of CTCF binding sites at TAD borders. At the earliest time point, when there is no zygotic transcription, the genome is highly structured. After zygotic genome activation (ZGA), the genome loses structural features, which are re-established throughout early development. Despite the absence of structural features, we see clustering of super-enhancers in the 3D genome. Our results provide insight into vertebrate genome organization and demonstrate that the developing zebrafish embryo is a powerful model system to study the dynamics of nuclear organization.

Enhancers reside in a unique epigenetic environment during early zebrafish development.

BACKGROUND: Enhancers, not promoters, are the most dynamic in their DNA methylation status throughout development and differentiation. Generally speaking, enhancers that are primed to or actually drive gene expression are characterized by relatively low levels of DNA methylation (hypo-methylation), while inactive enhancers display hyper-methylation of the underlying DNA. The direct functional significance of the DNA methylation state of enhancers is, however, unclear for most loci. RESULTS: In contrast to conventional epigenetic interactions at enhancers, we find that DNA methylation status and enhancer activity during early zebrafish development display very unusual correlation characteristics: hypo-methylation is a unique feature of primed enhancers whereas active enhancers are generally hyper-methylated. The hypo-methylated enhancers that we identify (hypo-enhancers) are enriched close to important transcription factors that act later in development. Interestingly, hypo-enhancers are de-methylated shortly before the midblastula transition and reside in a unique epigenetic environment. Finally, we demonstrate that hypo-enhancers do become active at later developmental stages and that they are physically associated with the transcriptional start site of target genes, irrespective of target gene activity. CONCLUSIONS: We demonstrate that early development in zebrafish embodies a time window characterized by non-canonical DNA methylation-enhancer relationships, including global DNA hypo-methylation of inactive enhancers and DNA hyper-methylation of active enhancers.

Retrospective analysis of the impact of HPV status and smoking on mucositis in patients with oropharyngeal squamous cell carcinoma treated with concurrent chemotherapy and radiotherapy.

OBJECTIVES: The standard concurrent radiotherapy and chemotherapy regimens for patients with oropharyngeal cancer are highly toxic. Human papillomavirus (HPV)-associated oropharyngeal squamous cell carcinoma (OPSCC) has recently emerged as a distinct biological and clinical entity with improved response to treatment and prognosis. A tailored therapeutic approach is needed to optimize patient care. The aim of our study was to investigate the impact of HPV and smoking status on early toxicities (primarily mucositis) associated with concurrent chemotherapy and radiotherapy in patients with OPSCC. MATERIALS AND METHODS: We retrospectively evaluated 72 consecutive patients with OPSCC and known HPV status treated with concurrent radiotherapy and chemotherapy at our institution. Treatment-related toxicities were stratified by smoking and HPV status and compared using univariate and multivariate logistic regression. RESULTS: HPV-positive patients had a 6.86-fold increase in the risk of having severe, grade 3-4 mucositis. This effect was preserved after adjusting for patient smoking status, nodal stage, radiotherapy technique and radiotherapy maximum dose. Additionally, HPV status had significant effect on the objective weight loss during treatment and at three months after treatment. Consistently, non-smokers had a significant 2.70-fold increase in the risk of developing severe mucositis. CONCLUSION: Risk factors for OPSCC modify the incidence of treatment-related early toxicities, with HPV-positive and non-smoking status correlating with increased risk of high grade mucositis and associated outcomes. Retrospective single-institution studies need to be interpreted cautiously. However, this finding is important to consider when designing therapeutic strategies for HPV-positive patients and merits further investigation in prospective clinical trials.

PID-1 is a novel factor that operates during 21U-RNA biogenesis in Caenorhabditis elegans.

The Piwi-piRNA pathway represents a small RNA-based mechanism responsible for the recognition and silencing of invading DNA. Biogenesis of piRNAs (21U-RNAs) is poorly understood. In Caenorhabditis elegans, the piRNA-binding Argonaute protein PRG-1 is the only known player acting downstream from precursor transcription. From a screen aimed at the isolation of piRNA-induced silencing-defective (Pid) mutations, we identified, among known Piwi pathway components, PID-1 as a novel player. PID-1 is a mostly cytoplasmic, germline-specific factor essential for 21U-RNA biogenesis, affecting an early step in the processing or transport of 21U precursor transcripts. We also show that maternal 21U-RNAs are essential to initiate silencing.

piRNA dynamics in divergent zebrafish strains reveal long-lasting maternal influence on zygotic piRNA profiles.

Transposable elements (TEs) are mobile genetic elements that can have many deleterious effects on the fitness of their host. The germline-specific PIWI pathway guards the genome against TEs, deriving its specificity from sequence complementarity between PIWI-bound small RNAs (piRNAs) and the TEs. The piRNAs are derived from so-called piRNA clusters. Recent studies have demonstrated that the piRNA repertoire can be adjusted to accommodate recent TE invasions by capturing invading TEs in piRNA loci. Thus far, no information concerning piRNA divergence is available from vertebrates. We present piRNA analyses of two relatively divergent zebrafish strains. We find that significant differences in the piRNA populations have accumulated, most notably among active class I TEs. This divergence can be split into differences in piRNA abundance per element and differences in sense/antisense polarity ratios. In crosses between animals of the different strains, many of these differences are resolved in the progeny. However, some differences remain, often leaning to the maternally contributed piRNA population. These differences can be detected at least two generations later. Our data illustrate, for the first time, the fluidity of piRNA populations in vertebrates and how the established diversity is transmitted to future generations.

Extremely stable Piwi-induced gene silencing in Caenorhabditis elegans.

In recent years, the Piwi pathway has been shown to regulate the silencing of mobile genetic elements. However, we know little about how Piwi pathways impose silencing and even less about trans-generational stability of Piwi-induced silencing. We demonstrate that the Caenorhabditis elegans Piwi protein PRG-1 can initiate an extremely stable form of gene silencing on a transgenic, single-copy target. This type of silencing is faithfully maintained over tens of generations in the absence of a functional Piwi pathway. Interestingly, RNAi can also trigger permanent gene silencing of a single-copy transgene and the phenomenon will be collectively referred to as RNA-induced epigenetic silencing (RNAe). RNAe can act in trans and is dependent on endogenous RNAi factors. The involvement of factors known to act in nuclear RNAi and the fact that RNAe is accompanied by repressive chromatin marks indicate that RNAe includes a transcriptional silencing component. Our results demonstrate that, at least in C. elegans, the Piwi pathway can impose a state of gene silencing that borders on 'permanently silent'. Such a property may be more widely conserved among Piwi pathways in different animals.

Differential impact of the HEN1 homolog HENN-1 on 21U and 26G RNAs in the germline of Caenorhabditis elegans.

RNA interference (RNAi)-related pathways affect gene activity by sequence-specific recruitment of Ago proteins to mRNA target molecules. The sequence specificity of this process stems from small RNA (sRNA) co-factors bound by the Ago protein. Stability of sRNA molecules in some pathways is in part regulated by Hen1-mediated methylation of their 3' ends. Here we describe the effects of the Caenorhabditis elegans HEN1 RNA-methyl-transferase homolog, HENN-1, on the different RNAi pathways in this nematode. We reveal differential effects of HENN-1 on the two pathways that are known to employ methylated sRNA molecules: the 26G and 21U pathways. Surprisingly, in the germline, stability of 21U RNAs, the C. elegans piRNAs, is only mildly affected by loss of methylation; and introduction of artificial 21U target RNA does not further destabilize non-methylated 21U RNAs. In contrast, most 26G RNAs display reduced stability and respond to loss of HENN-1 by displaying increased 3'-uridylation frequencies. Within the 26G RNA class, we find that specifically ERGO-1-bound 26G RNAs are modified by HENN-1, while ALG-3/ALG-4-bound 26G RNAs are not. Global gene expression analysis of henn-1 mutants reveals mild effects, including down-regulation of many germline-expressed genes. Our data suggest that, apart from direct effects of reduced 26G RNA levels of henn-1 on gene expression, most effects on global gene expression are indirect. These studies further refine our understanding of endogenous RNAi in C. elegans and the roles for Hen1 like enzymes in these pathways.

Tdrd1 acts as a molecular scaffold for Piwi proteins and piRNA targets in zebrafish.

Piwi proteins function in an RNAi-like pathway that silences transposons. Piwi-associated RNAs, also known as piRNAs, act as a guide to identify Piwi targets. The tudor domain-containing protein Tdrd1 has been linked to this pathway but its function has thus far remained unclear. We show that zebrafish Tdrd1 is required for efficient Piwi-pathway activity and proper nuage formation. Furthermore, we find that Tdrd1 binds both zebrafish Piwi proteins, Ziwi and Zili, and reveals sequence specificity in the interaction between Tdrd1 tudor domains and symmetrically dimethylated arginines (sDMAs) in Zili. Finally, we show that Tdrd1 complexes contain piRNAs and RNA molecules that are longer than piRNAs. We name these longer transcripts Tdrd1-associated transcripts (TATs). TATs likely represent cleaved Piwi pathway targets and may serve as piRNA biogenesis intermediates. Altogether, our data suggest that Tdrd1 acts as a molecular scaffold for Piwi proteins, bound through specific tudor domain-sDMA interactions, piRNAs and piRNA targets.

The inactive X chromosome adopts a unique three-dimensional conformation that is dependent on Xist RNA.

Three-dimensional topology of DNA in the cell nucleus provides a level of transcription regulation beyond the sequence of the linear DNA. To study the relationship between the transcriptional activity and the spatial environment of a gene, we used allele-specific chromosome conformation capture-on-chip (4C) technology to produce high-resolution topology maps of the active and inactive X chromosomes in female cells. We found that loci on the active X form multiple long-range interactions, with spatial segregation of active and inactive chromatin. On the inactive X, silenced loci lack preferred interactions, suggesting a unique random organization inside the inactive territory. However, escapees, among which is Xist, are engaged in long-range contacts with each other, enabling identification of novel escapees. Deletion of Xist results in partial refolding of the inactive X into a conformation resembling the active X without affecting gene silencing or DNA methylation. Our data point to a role for Xist RNA in shaping the conformation of the inactive X chromosome at least partially independent of transcription.

Interfacial friction and adhesion of polymer brushes.

A bead-probe lateral force microscopy (LFM) technique is used to characterize the interfacial friction and adhesion properties of polymer brushes. Our measurements attempt to relate the physical structure and chemical characteristics of the brush to their properties as thin-film, tethered lubricants. Brushes are synthesized at several chain lengths and surface coverages from polymer chains of polydimethylsiloxane (PDMS), polystyrene (PS), and a poly(propylene glycol)-poly(ethylene glycol) block copolymer (PPG/PEG). At high surface coverage, PDMS brushes manifest friction coefficients (COFs) that are among the lowest recorded for a dry lubricant film (µ ≈ 0.0024) and close to 1 order of magnitude lower than the COF of a bare silicon surface. Brushes synthesized from higher molar mass chains exhibit higher friction forces than those created using lower molar mass polymers. Increased grafting density of chains in the brush significantly reduces the COF by creating a uniform surface of stretched chains with a decreased surface viscosity. Brushes with lower surface tension and interfacial shear stresses manifest the lowest COF. In particular, PDMS chains exhibit COFs lower than PS by a factor of 3.7 and lower than PPG/PEG by a factor of 4.7. A scaling analysis conducted on the surface coverage (σ) in relation to the fraction (ε) of the friction force developing from adhesion predicts a universal relation ε ~ σ(4/3), which is supported by our experimental data.

Effect of pendent chains on the interfacial properties of thin polydimethylsiloxane (PDMS) networks.

The interfacial properties of end-linked polydimethylsiloxane (PDMS) films on silicon are examined. Thin cross-linked PDMS films (∼10 µm thick) were synthesized over a self-assembled monolayer supported on a silicon wafer. By systematically varying the concentration of monofunctional PDMS in a mixture with telechelic precursor molecules, structures ranging from near-ideal elastic networks to poorly cross-linked networks composed of a preponderance of dangling/pendent chains were synthesized. Lateral force microscopy (LFM) employing bead probes was used to quantify the effect of network structure on the interfacial friction coefficient and residual force. Indentation measurements employing an AFM in force mode were used to characterize the elastic modulus and the pull-off force for the films as a function of pendent chain content. These measurements were complemented with conventional mechanical rheometry measurements on similar thick network films to determine their bulk rheological properties. All networks studied manifested interfacial friction coefficients substantially lower than that of bare silicon. PDMS networks with the lowest pendent chain content displayed friction coefficients close to 1 order of magnitude lower than that of bare silicon, whereas networks with the highest pendent chain content manifested friction coefficients about 3 times lower than that of bare silicon. At intermediate sliding velocities, a crossover in the interfacial friction coefficient was observed, wherein cross-linked PDMS films with the least amount of pendent chains exhibit the highest friction coefficient. These observations are discussed in terms of the structure of the films and relaxation dynamics of elastic strands and dangling chains in tethered network films.

Hen1 is required for oocyte development and piRNA stability in zebrafish.

Piwi-interacting RNAs (piRNAs) are germ line-specific small RNA molecules that have a function in genome defence and germ cell development. They associate with a specific class of Argonaute proteins, named Piwi, and function through an RNA interference-like mechanism. piRNAs carry a 2'-O-methyl modification at their 3' end, which is added by the Hen1 enzyme. We show that zebrafish hen1 is specifically expressed in germ cells and is essential for maintaining a female germ line, whereas it is dispensable in the testis. Hen1 protein localizes to nuage through its C-terminal domain, but is not required for nuage formation. In hen1 mutant testes, piRNAs become uridylated and adenylated. Uridylation frequency is highest on retro-transposon-derived piRNAs and is accompanied by decreased piRNA levels and mild derepression of transposon transcripts. Altogether, our data suggest the existence of a uridylation-mediated 3'-5' exonuclease activity acting on piRNAs in zebrafish germ cells, which is counteracted by nuage-bound Hen1 protein. This system discriminates between piRNA targets and is required for ovary development and fully efficient transposon silencing.

Regulation of invasive behavior by vascular endothelial growth factor is HEF1-dependent.

We previously reported a vascular endothelial growth factor (VEGF) autocrine loop in head and neck squamous cell carcinoma (HNSCC) cell lines, supporting a role for VEGF in HNSCC tumorigenesis. Using a phosphotyrosine proteomics approach, we screened the HNSCC cell line, squamous cell carcinoma-9 for effectors of VEGFR2 signaling. A cluster of proteins involved in cell migration and invasion, including the p130Cas paralog, human enhancer of filamentation 1 (HEF1/Cas-L/Nedd9) was identified. HEF1 silencing and overexpression studies revealed a role for VEGF in regulating cell migration, invasion and matrix metalloproteinase (MMP) expression in a HEF1-dependent manner. Moreover, cells plated on extracellular matrix-coated coverslips showed enhanced invadopodia formation in response to VEGF that was HEF1-dependent. Immunolocalization revealed that HEF1 colocalized to invadopodia with MT1-MMP. Analysis of HNSCC tissue microarrays for HEF1 immunoreactivity revealed a 6.5-fold increase in the odds of having a metastasis with a high HEF1 score compared with a low HEF1 score. These findings suggest that HEF1 may be prognostic for advanced stage HNSCC. They also show for the first time that HEF1 is required for VEGF-mediated HNSCC cell migration and invasion, consistent with HEF1's recent identification as a metastatic regulator. These results support a strategy targeting VEGF:VEGFR2 in HNSCC therapeutics.

Dcr1 tracked down.

RNAi is essential for pericentromeric heterochromatic formation in S. pombe, and although Dcr1, the initiator protein of this process, has been biochemically well described, its subcellular localization has remained elusive. In this issue of Developmental Cell, Emmerth et al. now show that Dcr1 is dynamically shuttling between nucleus and cytoplasm, adding new insight into the subcellular mechanics of RNAi.

Cardiac transplantation: first-year experience in a community hospital.

The issue of decentralizing heart transplant services, formerly restricted to a few large medical centers, is currently under review by federal and state governments. We present the results of the first year of cardiac transplantation at a 385-bed community hospital. Twelve patients were selected according to generally accepted criteria from a pool of 24 referrals, all from within 75 miles of our institution. All patients were in New York Heart Association Class IV preoperatively. The one-year survival rate was found to be 82%, which is equivalent to that reported by established centers. All surviving patients were fully rehabilitated. Rates of infection and rejection were lower than expected, and costs were about half the national average. This series, in all likelihood, tests the limits to which the decentralization of cardiac transplant services can be taken. We conclude that cardiac transplantation can be accomplished at a community hospital with results, even for the first patients undergoing transplantation, comparable to those obtained by established programs at major medical centers.

Recovery from spinal cord injury.

There are many classifications of spinal cord injury. The anatomic definitions help in understanding the pathophysiology but are less useful in calculating recovery. Basic classifications, such as those by Guttman and Frankel, are useful for epidemiological studies, especially if one wishes to compare patients internationally or to compare data bases but this information, except in complete lesions, is too limited to allow more refined calculations as to the effectiveness of treatment modalities, especially in partial lesions. An accurate motor index, such as we have described, or such as the Yale index (2), is needed for more specific calculations of recovery if any prediction about the efficacy of various treatments are to be made. We have described our data from the mid-1970s in detail to explain our current treatment protocols. While nothing can be done for the completely paralyzed patient, our treatment of the acute (transiently concussed) near complete lesion and partial lesion must be tailored to maximize recovery. The fact that patients with initially complete motor lesion show some motor recovery in 15% indicates that some lesions are reversible and that we should try to learn to identify these patients and treat them to maximize their recovery. Certain severe partial lesions show poor recovery; we must improve therapy for these patients. Patients with partial lesions and 25-30% motor power below the injury have an 80% (+/- 20%) recovery; we should strive to attain and improve on these good results by appropriate medical and surgical management. In sharing with you our methods and data, we have provided these preliminary conclusions which influence our behavior as neurosurgeons, but it is equally clear that much work still needs to be done. This chapter should represent a reference source with which we can continue to work.