Normalization of a chromosomal contact map

ABSTRACT: BACKGROUND: Chromatin organization has been increasingly studied in relation with its important influence on DNA-related metabolic processes such as replication or regulation of gene expression. Since its original design ten years ago, capture of chromosome conformation (3C) has become an essential tool to investigate the overall conformation of chromosomes. It relies on the capture of long-range trans and cis interactions of chromosomal segments whose relative proportions in the final bank reflect their frequencies of interactions, hence their average spatial proximity. The recent coupling of 3C with deep sequencing approaches now allows the generation of high resolution genome-wide chromosomal contact maps. Different protocols have been used to generate such maps in various organisms. This includes mammals, drosophila and yeast. The massive amount of raw data generated by the genomic 3C has to be carefully processed to alleviate the various biases and byproducts generated by the experiments. Our study aims at proposing a simple normalization procedure to take into account these unwanted but inevitable events. RESULTS: Careful analysis of the raw data generated previously for budding yeast Saccharomyces cerevisiae led to the identification of three main biases affecting the final datasets, including an original bias resulting from the circularization of DNA molecules exhibiting specific lengths in accordance with laws from polymer physics. We then developed a simple normalization procedure to process the data and allow the generation of a normalized, highly contrasted, chromosomal contact map for S. cerevisiae. The same method was then extended to the first human genome contact map. Using the normalized data, we revisited the preferential interactions originally described between subsets of discrete chromosomal features. Notably, the detection of preferential interactions between tRNA in yeast and CTCF, PolII binding sites in human can vary with the normalization procedure used. CONCLUSIONS: We quantitatively reanalyzed the genomic 3C data obtained for S. cerevisiae, identified some of the biases inherent to the technique and proposed a simple normalization procedure to analyze them. Such an approach can be easily generalized for genomic 3C experiments in other organisms. More experiments and analysis will be necessary to reach optimal resolution and accuracies of the maps generated through these approaches. Working with cell population presenting highest levels of homogeneity will prove useful in this regards.     

Environment,preferred habitats and potential refugia for Pleistocene Homo in Southeast Asia

Multiple lines of evidence suggest that Southeast Asia was dominated by a mix of savannah, open woodlands, and evergreen forests throughout much of the Pleistocene. These conditions are ideal for early hominin subsistence; however, they would have been rare for much of the rest of Asia during glacial periods. We explore the possibility that Southeast Asia would have served as a refugium for hominins during these periods. In particular, we draw parallels with the population source and sink model proposed for northern Europe (Dennell et al., 2010), with Southeast Asia acting as a population source and northern China acting as a population sink.     

Traffic instabilities in self-organized pedestrian crowds

In human crowds as well as in many animal societies, local interactions among individuals often give rise to self-organized collective organizations that offer functional benefits to the group. For instance, flows of pedestrians moving in opposite directions spontaneously segregate into lanes of uniform walking directions. This phenomenon is often referred to as a smart collective pattern, as it increases the traffic efficiency with no need of external control. However, the functional benefits of this emergent organization have never been experimentally measured, and the underlying behavioral mechanisms are poorly understood. In this work, we have studied this phenomenon under controlled laboratory conditions. We found that the traffic segregation exhibits structural instabilities characterized by the alternation of organized and disorganized states, where the lifetime of well-organized clusters of pedestrians follow a stretched exponential relaxation process. Further analysis show that the inter-pedestrian variability of comfortable walking speeds is a key variable at the origin of the observed traffic perturbations. We show that the collective benefit of the emerging pattern is maximized when all pedestrians walk at the average speed of the group. In practice, however, local interactions between slow- and fast-walking pedestrians trigger global breakdowns of organization, which reduce the collective and the individual payoff provided by the traffic segregation. This work is a step ahead toward the understanding of traffic self-organization in crowds, which turns out to be modulated by complex behavioral mechanisms that do not always maximize the group's benefits. The quantitative understanding of crowd behaviors opens the way for designing bottom-up management strategies bound to promote the emergence of efficient collective behaviors in crowds.     

Marburg Virus VP35 Can Both Fully Coat the Backbone and Cap the Ends of dsRNA for Interferon Antagonism

Filoviruses, including Marburg virus (MARV) and Ebola virus (EBOV), cause fatal hemorrhagic fever in humans and non-human primates. All filoviruses encode a unique multi-functional protein termed VP35. The C-terminal double-stranded (ds)RNA-binding domain (RBD) of VP35 has been implicated in interferon antagonism and immune evasion. Crystal structures of the VP35 RBD from two ebolaviruses have previously demonstrated that the viral protein caps the ends of dsRNA. However, it is not yet understood how the expanses of dsRNA backbone, between the ends, are masked from immune surveillance during filovirus infection. Here, we report the crystal structure of MARV VP35 RBD bound to dsRNA. In the crystal structure, molecules of dsRNA stack end-to-end to form a pseudo-continuous oligonucleotide. This oligonucleotide is continuously and completely coated along its sugar-phosphate backbone by the MARV VP35 RBD. Analysis of dsRNA binding by dot-blot and isothermal titration calorimetry reveals that multiple copies of MARV VP35 RBD can indeed bind the dsRNA sugar-phosphate backbone in a cooperative manner in solution. Further, MARV VP35 RBD can also cap the ends of the dsRNA in solution, although this arrangement was not captured in crystals. Together, these studies suggest that MARV VP35 can both coat the backbone and cap the ends, and that for MARV, coating of the dsRNA backbone may be an essential mechanism by which dsRNA is masked from backbone-sensing immune surveillance molecules.     

A meta-analysis and genome-wide association study of platelet count and mean platelet volume in african americans

Several genetic variants associated with platelet count and mean platelet volume (MPV) were recently reported in people of European ancestry. In this meta-analysis of 7 genome-wide association studies (GWAS) enrolling African Americans, our aim was to identify novel genetic variants associated with platelet count and MPV. For all cohorts, GWAS analysis was performed using additive models after adjusting for age, sex, and population stratification. For both platelet phenotypes, meta-analyses were conducted using inverse-variance weighted fixed-effect models. Platelet aggregation assays in whole blood were performed in the participants of the GeneSTAR cohort. Genetic variants in ten independent regions were associated with platelet count (N?=?16,388) with p<5×10(-8) of which 5 have not been associated with platelet count in previous GWAS. The novel genetic variants associated with platelet count were in the following regions (the most significant SNP, closest gene, and p-value): 6p22 (rs12526480, LRRC16A, p?=?9.1×10(-9)), 7q11 (rs13236689, CD36, p?=?2.8×10(-9)), 10q21 (rs7896518, JMJD1C, p?=?2.3×10(-12)), 11q13 (rs477895, BAD, p?=?4.9×10(-8)), and 20q13 (rs151361, SLMO2, p?=?9.4×10(-9)). Three of these loci (10q21, 11q13, and 20q13) were replicated in European Americans (N?=?14,909) and one (11q13) in Hispanic Americans (N?=?3,462). For MPV (N?=?4,531), genetic variants in 3 regions were significant at p<5×10(-8), two of which were also associated with platelet count. Previously reported regions that were also significant in this study were 6p21, 6q23, 7q22, 12q24, and 19p13 for platelet count and 7q22, 17q11, and 19p13 for MPV. The most significant SNP in 1 region was also associated with ADP-induced maximal platelet aggregation in whole blood (12q24). Thus through a meta-analysis of GWAS enrolling African Americans, we have identified 5 novel regions associated with platelet count of which 3 were replicated in other ethnic groups. In addition, we also found one region associated with platelet aggregation that may play a potential role in atherothrombosis.     

Genome Structure and Reproductive Behaviour Influence the Evolutionary Potential of a Fungal Phytopathogen

Modern agriculture favours the selection and spread of novel plant diseases. Furthermore, crop genetic resistance against pathogens is often rendered ineffective within a few years of its commercial deployment. Leptosphaeria maculans, the cause of phoma stem canker of oilseed rape, develops gene-for-gene interactions with its host plant, and has a high evolutionary potential to render ineffective novel sources of resistance in crops. Here, we established a four-year field experiment to monitor the evolution of populations confronted with the newly released Rlm7 resistance and to investigate the nature of the mutations responsible for virulence against Rlm7. A total of 2551 fungal isolates were collected from experimental crops of a Rlm7 cultivar or a cultivar without Rlm7. All isolates were phenotyped for virulence and a subset was genotyped with neutral genetic markers. Virulent isolates were investigated for molecular events at the AvrLm4-7 locus. Whilst virulent isolates were not found in neighbouring crops, their frequency had reached 36% in the experimental field after four years. An extreme diversity of independent molecular events leading to virulence was identified in populations, with large-scale Repeat Induced Point mutations or complete deletion of AvrLm4-7 being the most frequent. Our data suggest that increased mutability of fungal genes involved in the interactions with plants is directly related to their genomic environment and reproductive system. Thus, rapid allelic diversification of avirulence genes can be generated in L. maculans populations in a single field provided that large population sizes and sexual reproduction are favoured by agricultural practices.     

Genomic Variation and Its Impact on Gene Expression in Drosophila melanogaster

Understanding the relationship between genetic and phenotypic variation is one of the great outstanding challenges in biology. To meet this challenge, comprehensive genomic variation maps of human as well as of model organism populations are required. Here, we present a nucleotide resolution catalog of single-nucleotide, multi-nucleotide, and structural variants in 39 Drosophila melanogaster Genetic Reference Panel inbred lines. Using an integrative, local assembly-based approach for variant discovery, we identify more than 3.6 million distinct variants, among which were more than 800,000 unique insertions, deletions (indels), and complex variants (1 to 6,000 bp). While the SNP density is higher near other variants, we find that variants themselves are not mutagenic, nor are regions with high variant density particularly mutation-prone. Rather, our data suggest that the elevated SNP density around variants is mainly due to population-level processes. We also provide insights into the regulatory architecture of gene expression variation in adult flies by mapping cis-expression quantitative trait loci (cis-eQTLs) for more than 2,000 genes. Indels comprise around 10% of all cis-eQTLs and show larger effects than SNP cis-eQTLs. In addition, we identified two-fold more gene associations in males as compared to females and found that most cis-eQTLs are sex-specific, revealing a partial decoupling of the genomic architecture between the sexes as well as the importance of genetic factors in mediating sex-biased gene expression. Finally, we performed RNA-seq-based allelic expression imbalance analyses in the offspring of crosses between sequenced lines, which revealed that the majority of strong cis-eQTLs can be validated in heterozygous individuals.     

Smyd3 regulates cancer cell phenotypes and catalyzes histone H4 lysine 5 methylation

Smyd3 is a lysine methyltransferase implicated in chromatin and cancer regulation. Here we show that Smyd3 catalyzes histone H4 methylation at lysine 5 (H4K5me). This novel histone methylation mark is detected in diverse cell types and its formation is attenuated by depletion of Smyd3 protein. Further, Smyd3-driven cancer cell phenotypes require its enzymatic activity. Thus, Smyd3, via H4K5 methylation, provides a potential new link between chromatin dynamics and neoplastic disease.