Familial recurrent hydatidiform mole (RHM) is a maternal-effect autosomal recessive disorder usually associated with mutations of the NLRP7 gene. It is characterized by HM with excessive trophoblastic proliferation, which mimics the appearance of androgenetic molar conceptuses despite their diploid biparental constitution. It has been proposed that the phenotypes of both types of mole are associated with aberrant genomic imprinting. However no systematic analyses for imprinting defects have been reported. Here, we present the genome-wide methylation profiles of both spontaneous androgenetic and biparental NLRP7 defective molar tissues. We observe total paternalization of all ubiquitous and placenta-specific differentially methylated regions (DMRs) in four androgenetic moles; namely gain of methylation at paternally methylated loci and absence of methylation at maternally methylated regions. The methylation defects observed in five RHM biopsies from NLRP7 defective patients are restricted to lack-of-methylation at maternal DMRs. Surprisingly RHMs from two sisters with the same missense mutations, as well as consecutive RHMs from one affected female show subtle allelic methylation differences, suggesting inter-RHM variation. These epigenotypes are consistent with NLRP7 being a maternal-effect gene and involved in imprint acquisition in the oocyte. In addition, bioinformatic screening of the resulting methylation datasets identified over sixty loci with methylation profiles consistent with imprinting in the placenta, of which we confirm 22 as novel maternally methylated loci. These observations strongly suggest that the molar phenotypes are due to defective placenta-specific imprinting and over-expression of paternally expressed transcripts, highlighting that maternal-effect mutations of NLRP7 are associated with the most severe form of multi-locus imprinting defects in humans.

Autism spectrum disorders are associated with defects in social response and communication that often occur in the context of intellectual disability. Rett syndrome is one example in which epilepsy, motor impairment, and motor disturbance may co-occur. Mutations in histone demethylases are known to occur in several of these syndromes. Herein, we aimed to identify whether mutations in the candidate histone demethylase JMJD1C (jumonji domain containing 1C) are implicated in these disorders.

Despite the increasing knowledge about DNA methylation, the understanding of human epigenome evolution is in its infancy. Using whole genome bisulfite sequencing we identified hundreds of differentially methylated regions (DMRs) in humans compared to non-human primates and estimated that ∼25% of these regions were detectable throughout several human tissues. Human DMRs were enriched for specific histone modifications and the majority were located distal to transcription start sites, highlighting the importance of regions outside the direct regulatory context. We also found a significant excess of endogenous retrovirus elements in human-specific hypomethylated.We reported for the first time a close interplay between inter-species genetic and epigenetic variation in regions of incomplete lineage sorting, transcription factor binding sites and human differentially hypermethylated regions. Specifically, we observed an excess of human-specific substitutions in transcription factor binding sites located within human DMRs, suggesting that alteration of regulatory motifs underlies some human-specific methylation patterns. We also found that the acquisition of DNA hypermethylation in the human lineage is frequently coupled with a rapid evolution at nucleotide level in the neighborhood of these CpG sites. Taken together, our results reveal new insights into the mechanistic basis of human-specific DNA methylation patterns and the interpretation of inter-species non-coding variation.

Here, we report DNA methylation and hydroxymethylation dynamics at nucleotide resolution using C/EBPα-enhanced reprogramming of B cells into induced pluripotent cells (iPSCs). We observed successive waves of hydroxymethylation at enhancers, concomitant with a decrease in DNA methylation, suggesting active demethylation. Consistent with this finding, ablation of the DNA demethylase Tet2 almost completely abolishes reprogramming. C/EBPα, Klf4, and Tfcp2l1 each interact with Tet2 and recruit the enzyme to specific DNA sites. During reprogramming, some of these sites maintain high levels of 5hmC, and enhancers and promoters of key pluripotency factors become demethylated as early as 1 day after Yamanaka factor induction. Surprisingly, methylation changes precede chromatin opening in distinct chromatin regions, including Klf4 bound sites, revealing a pioneer factor activity associated with alternation in DNA methylation. Rapid changes in hydroxymethylation similar to those in B cells were also observed during compound-accelerated reprogramming of fibroblasts into iPSCs, highlighting the generality of our observations.

Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ∼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

T47D_rep2 and b1913e6c1_51720e9cf were 2 Hi-C samples. They were born and processed at the same time, yet their fates were very different. The life of b1913e6c1_51720e9cf was simple and fruitful, while that of T47D_rep2 was full of accidents and sorrow. At the heart of these differences lies the fact that b1913e6c1_51720e9cf was born under a lab culture of Documentation, Automation, Traceability, and Autonomy and compliance with the FAIR Principles. Their lives are a lesson for those who wish to embark on the journey of managing high-throughput sequencing data.

Modern dog breeds display traits that are either breed-specific or shared by a few breeds as a result of genetic bottlenecks during the breed creation process and artificial selection for breed standards. Selective sweeps in the genome result from strong selection and can be detected as a reduction or elimination of polymorphism in a given region of the genome.

Cells have the ability to sense, respond and adapt to environmental fluctuations. Stress causes a massive reorganization of the transcriptional program. Many examples of histone post-translational modifications (PTMs) have been associated with transcriptional activation or repression under steady-state growth conditions. Comparatively less is known about the role of histone PTMs in the cellular adaptive response to stress. Here, we performed high-throughput genetic screenings that provide a novel global map of the histone residues required for transcriptional reprogramming in response to heat and osmotic stress. Of note, we observed that the histone residues needed depend on the type of gene and/or stress, thereby suggesting a 'personalized', rather than general, subset of histone requirements for each chromatin context. In addition, we identified a number of new residues that unexpectedly serve to regulate transcription. As a proof of concept, we characterized the function of the histone residues H4-S47 and H4-T30 in response to osmotic and heat stress, respectively. Our results uncover novel roles for the kinases Cla4 and Ste20, yeast homologs of the mammalian PAK2 family, and the Ste11 MAPK as regulators of H4-S47 and H4-T30, respectively. This study provides new insights into the role of histone residues in transcriptional regulation under stress conditions.

Living organisms are continuously challenged by changes in their environment that can propagate to stresses at the cellular level, such as rapid changes in osmolarity or oxygen tension. To survive these sudden changes, cells have developed stress-responsive mechanisms that tune cellular processes. The response of Saccharomyces cerevisiae to osmostress includes a massive reprogramming of gene expression. Identifying the inherent features of stress-responsive genes is of significant interest for understanding the basic principles underlying the rewiring of gene expression upon stress. Here, we generated a comprehensive catalog of osmostress-responsive genes from 5 independent RNA-seq experiments. We explored 30 features of yeast genes and found that 25 (83%) were distinct in osmostress-responsive genes. We then identified 13 non-redundant minimal osmostress gene traits and used statistical modeling to rank the most stress-predictive features. Intriguingly, the most relevant features of osmostress-responsive genes are the number of transcription factors targeting them and gene conservation. Using data on HeLa samples, we showed that the same features that define yeast osmostress-responsive genes can predict osmostress-responsive genes in humans, but with changes in the rank-ordering of feature-importance. Our study provides a holistic understanding of the basic principles of the regulation of stress-responsive gene expression across eukaryotes.

Regulation of cell volume is essential for tissue homeostasis and cell viability. In response to hypertonic stress, cells need rapid electrolyte influx to compensate water loss and to prevent cell death in a process known as regulatory volume increase (RVI). However, the molecular component able to trigger such a process was unknown to date. Using a genome-wide CRISPR/Cas9 screen, we identified LRRC8A, which encodes a chloride channel subunit, as the gene most associated with cell survival under hypertonic conditions. Hypertonicity activates the p38 stress-activated protein kinase pathway and its downstream MSK1 kinase, which phosphorylates and activates LRRC8A. LRRC8A-mediated Cl- efflux facilitates activation of the with-no-lysine (WNK) kinase pathway, which in turn, promotes electrolyte influx via Na+/K+/2Cl- cotransporter (NKCC) and RVI under hypertonic stress. LRRC8A-S217A mutation impairs channel activation by MSK1, resulting in reduced RVI and cell survival. In summary, LRRC8A is key to bidirectional osmotic stress responses and cell survival under hypertonic conditions.

Here, we describe how the speed of C/EBPα-induced B cell to macrophage transdifferentiation (BMT) can be regulated, using both mouse and human models. The identification of a mutant of C/EBPα (C/EBPαR35A) that greatly accelerates BMT helped to illuminate the mechanism. Thus, incoming C/EBPα binds to PU.1, an obligate partner expressed in B cells, leading to the release of PU.1 from B cell enhancers, chromatin closing and silencing of the B cell program. Released PU.1 redistributes to macrophage enhancers newly occupied by C/EBPα, causing chromatin opening and activation of macrophage genes. All these steps are accelerated by C/EBPαR35A, initiated by its increased affinity for PU.1. Wild-type C/EBPα is methylated by Carm1 at arginine 35 and the enzyme's perturbations modulate BMT velocity as predicted from the observations with the mutant. Increasing the proportion of unmethylated C/EBPα in granulocyte/macrophage progenitors by inhibiting Carm1 biases the cell's differentiation toward macrophages, suggesting that cell fate decision velocity and lineage directionality are closely linked processes.

Although a variety of genetic changes have been implicated in causing phenotypic differences among dogs, the role of copy number variants (CNVs) and their impact on phenotypic variation is still poorly understood. Further, very limited knowledge exists on structural variation in the gray wolf, the ancestor of the dog, or other closely related wild canids. Documenting CNVs variation in wild canids is essential to identify ancestral states and variation that may have appeared after domestication.

The intestinal stem cell niche provides cues that actively maintain gut homeostasis. Dysregulation of these cues may compromise intestinal regeneration upon tissue insult and/or promote tumor growth. Here, we identify secreted phospholipases A2 (sPLA2s) as stem cell niche factors with context-dependent functions in the digestive tract. We show that group IIA sPLA2, a known genetic modifier of mouse intestinal tumorigenesis, is expressed by Paneth cells in the small intestine, while group X sPLA2 is expressed by Paneth/goblet-like cells in the colon. During homeostasis, group IIA/X sPLA2s inhibit Wnt signaling through intracellular activation of Yap1. However, upon inflammation they are secreted into the intestinal lumen, where they promote prostaglandin synthesis and Wnt signaling. Genetic ablation of both sPLA2s improves recovery from inflammation but increases colon cancer susceptibility due to release of their homeostatic Wnt-inhibitory role. This "trade-off" effect suggests sPLA2s have important functions as genetic modifiers of inflammation and colon cancer.

Genomic experiments (e.g. differential gene expression, single-nucleotide polymorphism association) typically produce ranked list of genes. We present a simple but powerful approach which uses protein-protein interaction data to detect sub-networks within such ranked lists of genes or proteins. We performed an exhaustive study of network parameters that allowed us concluding that the average number of components and the average number of nodes per component are the parameters that best discriminate between real and random networks. A novel aspect that increases the efficiency of this strategy in finding sub-networks is that, in addition to direct connections, also connections mediated by intermediate nodes are considered to build up the sub-networks. The possibility of using of such intermediate nodes makes this approach more robust to noise. It also overcomes some limitations intrinsic to experimental designs based on differential expression, in which some nodes are invariant across conditions. The proposed approach can also be used for candidate disease-gene prioritization. Here, we demonstrate the usefulness of the approach by means of several case examples that include a differential expression analysis in Fanconi Anemia, a genome-wide association study of bipolar disorder and a genome-scale study of essentiality in cancer genes. An efficient and easy-to-use web interface (available at http://www.babelomics.org) based on HTML5 technologies is also provided to run the algorithm and represent the network.

Understanding the aspects of the cell functionality that account for disease or drug action mechanisms is one of the main challenges in the analysis of genomic data and is on the basis of the future implementation of precision medicine.

In conformational disorders, it is not evident which amyloid aggregates affect specific molecular mechanisms or cellular pathways, which cause disease because of their quantity and mechanical features and which states in aggregate formation are pathogenic. Due to the increasing consensus that prefibrillar oligomers play a major role in conformational diseases, there is a growing interest in understanding the characteristics of metastable polypeptide associations.

Forced transcription factor expression can transdifferentiate somatic cells into other specialised cell types or reprogram them into induced pluripotent stem cells (iPSCs) with variable efficiency. To better understand the heterogeneity of these processes, we used single-cell RNA sequencing to follow the transdifferentation of murine pre-B cells into macrophages as well as their reprogramming into iPSCs. Even in these highly efficient systems, there was substantial variation in the speed and path of fate conversion. We predicted and validated that these differences are inversely coupled and arise in the starting cell population, with Mychigh large pre-BII cells transdifferentiating slowly but reprogramming efficiently and Myclow small pre-BII cells transdifferentiating rapidly but failing to reprogram. Strikingly, differences in Myc activity predict the efficiency of reprogramming across a wide range of somatic cell types. These results illustrate how single cell expression and computational analyses can identify the origins of heterogeneity in cell fate conversion processes.

A hallmark of chromosome organization is the partition into transcriptionally active A and repressed B compartments, and into topologically associating domains (TADs). Both structures were regarded to be absent from the inactive mouse X chromosome, but to be re-established with transcriptional reactivation and chromatin opening during X-reactivation. Here, we combine a tailor-made mouse iPSC reprogramming system and high-resolution Hi-C to produce a time course combining gene reactivation, chromatin opening and chromosome topology during X-reactivation. Contrary to previous observations, we observe A/B-like compartments on the inactive X harbouring multiple subcompartments. While partial X-reactivation initiates within a compartment rich in X-inactivation escapees, it then occurs rapidly along the chromosome, concomitant with downregulation of Xist. Importantly, we find that TAD formation precedes transcription and initiates from Xist-poor compartments. Here, we show that TAD formation and transcriptional reactivation are causally independent during X-reactivation while establishing Xist as a common denominator.

The COVID-19 pandemic is raging. It revealed the importance of rapid scientific advancement towards understanding and treating new diseases. To address this challenge, we adapt an explainable artificial intelligence algorithm for data fusion and utilize it on new omics data on viral-host interactions, human protein interactions, and drugs to better understand SARS-CoV-2 infection mechanisms and predict new drug-target interactions for COVID-19. We discover that in the human interactome, the human proteins targeted by SARS-CoV-2 proteins and the genes that are differentially expressed after the infection have common neighbors central in the interactome that may be key to the disease mechanisms. We uncover 185 new drug-target interactions targeting 49 of these key genes and suggest re-purposing of 149 FDA-approved drugs, including drugs targeting VEGF and nitric oxide signaling, whose pathways coincide with the observed COVID-19 symptoms. Our integrative methodology is universal and can enable insight into this and other serious diseases.

Fusion genes are typically identified by RNA sequencing (RNA-seq) without elucidating the causal genomic breakpoints. However, non-poly(A)-enriched RNA-seq contains large proportions of intronic reads that also span genomic breakpoints.