Genome organization controls transcriptional dynamics during development.

Past studies offer contradictory claims for the role of genome organization in the regulation of gene activity. Here, we show through high-resolution chromosome conformation analysis that the Drosophila genome is organized by two independent classes of regulatory sequences, tethering elements and insulators. Quantitative live imaging and targeted genome editing demonstrate that this two-tiered organization is critical for the precise temporal dynamics of Hox gene transcription during development. Tethering elements mediate long-range enhancer-promoter interactions and foster fast activation kinetics. Conversely, the boundaries of topologically associating domains (TADs) prevent spurious interactions with enhancers and silencers located in neighboring TADs. These two levels of genome organization operate independently of one another to ensure precision of transcriptional dynamics and the reliability of complex patterning processes.

HP1 drives de novo 3D genome reorganization in early Drosophila embryos.

Fundamental features of 3D genome organization are established de novo in the early embryo, including clustering of pericentromeric regions, the folding of chromosome arms and the segregation of chromosomes into active (A-) and inactive (B-) compartments. However, the molecular mechanisms that drive de novo organization remain unknown1,2. Here, by combining chromosome conformation capture (Hi-C), chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq), 3D DNA fluorescence in situ hybridization (3D DNA FISH) and polymer simulations, we show that heterochromatin protein 1a (HP1a) is essential for de novo 3D genome organization during Drosophila early development. The binding of HP1a at pericentromeric heterochromatin is required to establish clustering of pericentromeric regions. Moreover, HP1a binding within chromosome arms is responsible for overall chromosome folding and has an important role in the formation of B-compartment regions. However, depletion of HP1a does not affect the A-compartment, which suggests that a different molecular mechanism segregates active chromosome regions. Our work identifies HP1a as an epigenetic regulator that is involved in establishing the global structure of the genome in the early embryo.

High-resolution single-cell 3D-models of chromatin ensembles during Drosophila embryogenesis.

Single-cell chromatin studies provide insights into how chromatin structure relates to functions of individual cells. However, balancing high-resolution and genome wide-coverage remains challenging. We describe a computational method for the reconstruction of large 3D-ensembles of single-cell (sc) chromatin conformations from population Hi-C that we apply to study embryogenesis in Drosophila. With minimal assumptions of physical properties and without adjustable parameters, our method generates large ensembles of chromatin conformations via deep-sampling. Our method identifies specific interactions, which constitute 5-6% of Hi-C frequencies, but surprisingly are sufficient to drive chromatin folding, giving rise to the observed Hi-C patterns. Modeled sc-chromatins quantify chromatin heterogeneity, revealing significant changes during embryogenesis. Furthermore, >50% of modeled sc-chromatin maintain topologically associating domains (TADs) in early embryos, when no population TADs are perceptible. Domain boundaries become fixated during development, with strong preference at binding-sites of insulator-complexes upon the midblastula transition. Overall, high-resolution 3D-ensembles of sc-chromatin conformations enable further in-depth interpretation of population Hi-C, improving understanding of the structure-function relationship of genome organization.

The Impact of Drosophila Awd/NME1/2 Levels on Notch and Wg Signaling Pathways.

Awd, the Drosophila homologue of NME1/2 metastasis suppressors, plays key roles in many signaling pathways. Mosaic analysis of the null awdJ2A4 allele showed that loss of awd gene function blocks Notch signaling and the expression of its target genes including the Wingless (Wg/Wnt1) morphogen. We also showed that RNA interference (RNAi)-mediated awd silencing (awdi) in larval wing disc leads to chromosomal instability (CIN) and to Jun amino-terminal kinases (JNK)-mediated cell death. Here we show that this cell death is independent of p53 activity. Based on our previous finding showing that forced survival of awdi-CIN cells leads to aneuploidy without the hyperproliferative effect, we investigated the Wg expression in awdi wing disc cells. Interestingly, the Wg protein is expressed in its correct dorso-ventral domain but shows an altered cellular distribution which impairs its signaling. Further, we show that RNAi-mediated knock down of awd in wing discs does not affect Notch signaling. Thus, our analysis of the hypomorphic phenotype arising from awd downregulation uncovers a dose-dependent effect of Awd in Notch and Wg signaling.

Enhanced nucleotide chemistry and toehold nanotechnology reveals lncRNA spreading on chromatin.

Understanding the targeting and spreading patterns of long non-coding RNAs (lncRNAs) on chromatin requires a technique that can detect both high-intensity binding sites and reveal genome-wide changes in spreading patterns with high precision and confidence. Here we determine lncRNA localization using biotinylated locked nucleic acid (LNA)-containing oligonucleotides with toehold architecture capable of hybridizing to target RNA through strand-exchange reaction. During hybridization, a protecting strand competitively displaces contaminating species, leading to highly specific RNA capture of individual RNAs. Analysis of Drosophila roX2 lncRNA using this approach revealed that heat shock, unlike the unfolded protein response, leads to reduced spreading of roX2 on the X chromosome, but surprisingly also to relocalization to sites on autosomes. Our results demonstrate that this improved hybridization capture approach can reveal previously uncharacterized changes in the targeting and spreading of lncRNAs on chromatin.

Cytogenetic markers reveal a reinforcement of variation in the tension zone between chromosome races in the brachypterous grasshopper Podisma sapporensis Shir. on Hokkaido Island.

The cytogenetic characteristics of the grasshopper Podisma sapporensis (two races 2n = 23♂ X0/XX and 2n = 22♂ neo-XY/neo-XX) were analysed through fluorescence in situ hybridization with rDNA and telomeric DNA probes, C-banding, fluorochrome and silver staining. For the first time, samples from the neighbourhood of a hybrid population (i.e., Mikuni Pass population) were studied. Our results indicated a significant degree of chromosomal differentiation between P. sapporensis races when comparing the number and position of the rDNA sites, as well as the heterochromatin composition and distribution obtained by C-banding and DAPI/CMA3 staining. Telomeric signals were usually detected at the distal and/or subdistal position of the autosomes; however, some chromosome ends lacked signals, probably due to a low number of telomeric repeats. On the other hand, telomeric DNA sequences were found as interstitial telomeric repeats in some autosomes, which can trigger a variety of genome instability. B chromosomes were found in specimens belonging to both main races from nine out of 22 localities. Four types of X chromosomes in the X0/XX race were identified. It was concluded that the physical mapping of rDNA sequences and heterochromatin are useful as additional markers for understanding the phylogeographic patterns of cytogenetic differentiation in P. sapporensis populations.

Large distances separate coregulated genes in living Drosophila embryos.

Transcriptional enhancers are short segments of DNA that switch genes on and off in response to a variety of cellular signals. Many enhancers map quite far from their target genes, on the order of tens or even hundreds of kilobases. There is extensive evidence that remote enhancers are brought into proximity with their target promoters via long-range looping interactions. However, the exact physical distances of these enhancer-promoter interactions remain uncertain. Here, we employ high-resolution imaging of living Drosophila embryos to visualize the distances separating linked genes that are coregulated by a shared enhancer. Cotransvection assays (linked genes on separate homologs) suggest a surprisingly large distance during transcriptional activity: at least 100-200 nm. Similar distances were observed when a shared enhancer was placed into close proximity with linked reporter genes in cis. These observations are consistent with the occurrence of "transcription hubs," whereby clusters (or condensates) of multiple RNA polymerase II complexes and associated cofactors are periodically recruited to active promoters. The dynamics of this process might be responsible for rapid fluctuations in the distances separating the transcription of coregulated reporter genes during transvection. We propose that enhancer-promoter communication depends on a combination of classical looping and linking models.

Unprecedented reorganization of holocentric chromosomes provides insights into the enigma of lepidopteran chromosome evolution.

Chromosome evolution presents an enigma in the mega-diverse Lepidoptera. Most species exhibit constrained chromosome evolution with nearly identical haploid chromosome counts and chromosome-level gene collinearity among species more than 140 million years divergent. However, a few species possess radically inflated chromosomal counts due to extensive fission and fusion events. To address this enigma of constraint in the face of an exceptional ability to change, we investigated an unprecedented reorganization of the standard lepidopteran chromosome structure in the green-veined white butterfly (Pieris napi). We find that gene content in P. napi has been extensively rearranged in large collinear blocks, which until now have been masked by a haploid chromosome number close to the lepidopteran average. We observe that ancient chromosome ends have been maintained and collinear blocks are enriched for functionally related genes suggesting both a mechanism and a possible role for selection in determining the boundaries of these genome-wide rearrangements.

Sister centromere fusion during meiosis I depends on maintaining cohesins and destabilizing microtubule attachments.

Sister centromere fusion is a process unique to meiosis that promotes co-orientation of the sister kinetochores, ensuring they attach to microtubules from the same pole during metaphase I. We have found that the kinetochore protein SPC105R/KNL1 and Protein Phosphatase 1 (PP1-87B) regulate sister centromere fusion in Drosophila oocytes. The analysis of these two proteins, however, has shown that two independent mechanisms maintain sister centromere fusion. Maintenance of sister centromere fusion by SPC105R depends on Separase, suggesting cohesin proteins must be maintained at the core centromeres. In contrast, maintenance of sister centromere fusion by PP1-87B does not depend on either Separase or WAPL. Instead, PP1-87B maintains sister centromeres fusion by regulating microtubule dynamics. We demonstrate that this regulation is through antagonizing Polo kinase and BubR1, two proteins known to promote stability of kinetochore-microtubule (KT-MT) attachments, suggesting that PP1-87B maintains sister centromere fusion by inhibiting stable KT-MT attachments. Surprisingly, C(3)G, the transverse element of the synaptonemal complex (SC), is also required for centromere separation in Pp1-87B RNAi oocytes. This is evidence for a functional role of centromeric SC in the meiotic divisions, that might involve regulating microtubule dynamics. Together, we propose two mechanisms maintain co-orientation in Drosophila oocytes: one involves SPC105R to protect cohesins at sister centromeres and another involves PP1-87B to regulate spindle forces at end-on attachments.

A determining factor for insect feeding preference in the silkworm, Bombyx mori.

Feeding preference is critical for insect adaptation and survival. However, little is known regarding the determination of insect feeding preference, and the genetic basis is poorly understood. As a model lepidopteran insect with economic importance, the domesticated silkworm, Bombyx mori, is a well-known monophagous insect that predominantly feeds on fresh mulberry leaves. This species-specific feeding preference provides an excellent model for investigation of host-plant selection of insects, although the molecular mechanism underlying this phenomenon remains unknown. Here, we describe the gene GR66, which encodes a putative bitter gustatory receptor (GR) that is responsible for the mulberry-specific feeding preference of B. mori. With the aid of a transposon-based, clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein-9 nuclease (Cas9) system, the GR66 locus was genetically mutated, and homozygous mutant silkworm strains with truncated gustatory receptor 66 (GR66) proteins were established. GR66 mutant larvae acquired new feeding activity, exhibiting the ability to feed on a number of plant species in addition to mulberry leaves, including fresh fruits and grain seeds that are not normally consumed by wild-type (WT) silkworms. Furthermore, a feeding choice assay revealed that the mutant larvae lost their specificity for mulberry. Overall, our findings provide the first genetic and phenotypic evidences that a single bitter GR is a major factor affecting the insect feeding preference.

Gene expression changes elicited by a parasitic B chromosome in the grasshopper Eyprepocnemis plorans are consistent with its phenotypic effects.

Parasitism evokes adaptive physiological changes in the host, many of which take place through gene expression changes. This response can be more or less local, depending on the organ or tissue affected by the parasite, or else systemic when the parasite affects the entire host body. The most extreme of the latter cases is intragenomic parasitism, where the parasite is present in all host nuclei as any other genomic element. Here, we show the molecular crosstalk between a parasitic chromosome (also named B chromosome) and the host genome, manifested through gene expression changes. The transcriptome analysis of 0B and 1B females of the grasshopper Eyprepocnemis plorans, validated by a microarray experiment performed on four B-lacking and five B-carrying females, revealed changes in gene expression for 188 unigenes being consistent in both experiments. Once discarded B-derived transcripts, there were 46 differentially expressed genes (30 up- and 16 downregulated) related with the adaptation of the host genome to the presence of the parasitic chromosome. Interestingly, the functions of these genes could explain some of the most important effects of B chromosomes, such as nucleotypic effects derived from the additional DNA they represent, chemical defense and detoxification, protein modification and response to stress, ovary function, and regulation of gene expression. Collectively, these changes uncover an intimate host-parasite interaction between A and B chromosomes during crucial steps of gene expression and protein function.

A new role for Drosophila Aurora-A in maintaining chromosome integrity.

Aurora-A is a conserved mitotic kinase overexpressed in many types of cancer. Growing evidence shows that Aurora-A plays a crucial role in DNA damage response (DDR) although this aspect has been less characterized. We isolated a new aur-A mutation, named aur-A949, in Drosophila, and we showed that it causes chromosome aberrations (CABs). In addition, aur-A949 mutants were sensitive to X-ray treatment and showed impaired γ-H2Av foci dissolution kinetics. To identify the pathway in which Aur-A works, we conducted an epistasis analysis by evaluating CAB frequencies in double mutants carrying aur-A949 mutation combined to mutations in genes related to DNA damage response (DDR). We found that mutations in tefu (ATM) and in the histone variant H2Av were epistatic over aur-A949 indicating that Aur-A works in DDR and that it is required for γ-H2Av foci dissolution. More interestingly, we found that a mutation in lig4, a gene belonging to the non-homologous end joining (NHEJ) repair pathway, was epistatic over aur-A949. Based on studies in other systems, which show that phosphorylation is important to target Lig4 for degradation, we hypothesized that in aur-A949 mutant cells, there is a persistence of Lig4 that could be, in the end, responsible for CABs. Finally, we observed a synergistic interaction between Aur-A and the homologous recombination (HR) repair system component Rad 51 in the process that converts chromatid deletions into isochromatid deletions. Altogether, these data indicate that Aur-A depletion can elicit chromosome damage. This conclusion should be taken into consideration, since some anticancer therapies are aimed at reducing Aurora-A expression.

Identification Key for the Chagas Disease Vectors of Five Brazilian States, Based on Cytogenetic Data.

Chagas disease is a public health problem caused by the protozoan Trypanosoma cruzi that affects about 8 million people worldwide. The main form of transmission of T. cruzi is vectorial, through triatomines feces contaminated with the parasite. All species are considered as potential vectors of T. cruzi. The main identification keys of these vectors are based only on morphological characters. However, there are very similar or even same species (cryptic species) that may lead to wrong classification of the vectors. Therefore, we developed an identification key using cytogenetic data, to aid and help the correct classification of triatomines. From the cytogenetic characters, identification keys were created for the five Brazilian states (Alagoas, Amapá, Ceará, Roraima, and Santa Catarina). These data are important because the correct classification of triatomines helps directly the activity of the vector control programs.

New Evidence of the Monophyletic Relationship of the Genus Psammolestes Bergroth, 1911 (Hemiptera, Reduviidae, Triatominae).

The genus Psammolestes within the subfamily Triatominae and tribe Rhodniini comprises the species Psammolestes arthuri, Psammolestes coreodes, and Psammolestes tertius, all potential vectors of Chagas disease. A feature of Psammolestes is their close association with birds, which makes them an interesting model for evolutionary studies. We analyzed cytogenetically Psammolestes spp., with the aim of contributing to the genetic and evolutionary knowledge of these vectors. All species of the Psammolestes showed the same chromosomal characteristics: chromocenter formed only by sex chromosomes X and Y, karyotype 2n = 22 and constitutive heterochromatin, and AT base pairs restricted to the sex chromosome Y. These results corroborate the monophyly of the genus and lead to the hypothesis that during the derivation of P. tertius, P. coreodes, and P. arthuri from their common ancestor, there was no reorganization in the number or structure of chromosomes.

Pervasive genetic interactions modulate neurodevelopmental defects of the autism-associated 16p11.2 deletion in Drosophila melanogaster.

As opposed to syndromic CNVs caused by single genes, extensive phenotypic heterogeneity in variably-expressive CNVs complicates disease gene discovery and functional evaluation. Here, we propose a complex interaction model for pathogenicity of the autism-associated 16p11.2 deletion, where CNV genes interact with each other in conserved pathways to modulate expression of the phenotype. Using multiple quantitative methods in Drosophila RNAi lines, we identify a range of neurodevelopmental phenotypes for knockdown of individual 16p11.2 homologs in different tissues. We test 565 pairwise knockdowns in the developing eye, and identify 24 interactions between pairs of 16p11.2 homologs and 46 interactions between 16p11.2 homologs and neurodevelopmental genes that suppress or enhance cell proliferation phenotypes compared to one-hit knockdowns. These interactions within cell proliferation pathways are also enriched in a human brain-specific network, providing translational relevance in humans. Our study indicates a role for pervasive genetic interactions within CNVs towards cellular and developmental phenotypes.

How epigenome drives chromatin folding and dynamics, insights from efficient coarse-grained models of chromosomes.

The 3D organization of chromosomes is crucial for regulating gene expression and cell function. Many experimental and polymer modeling efforts are dedicated to deciphering the mechanistic principles behind chromosome folding. Chromosomes are long and densely packed-topologically constrained-polymers. The main challenges are therefore to develop adequate models and simulation methods to investigate properly the multi spatio-temporal scales of such macromolecules. Here, we proposed a generic strategy to develop efficient coarse-grained models for self-avoiding polymers on a lattice. Accounting accurately for the polymer entanglement length and the volumic density, we show that our simulation scheme not only captures the steady-state structural and dynamical properties of the system but also tracks the same dynamics at different coarse-graining. This strategy allows a strong power-law gain in numerical efficiency and offers a systematic way to define reliable coarse-grained null models for chromosomes and to go beyond the current limitations by studying long chromosomes during an extended time period with good statistics. We use our formalism to investigate in details the time evolution of the 3D organization of chromosome 3R (20 Mbp) in drosophila during one cell cycle (20 hours). We show that a combination of our coarse-graining strategy with a one-parameter block copolymer model integrating epigenomic-driven interactions quantitatively reproduce experimental data at the chromosome-scale and predict that chromatin motion is very dynamic during the cell cycle.

Similarity in replication timing between polytene and diploid cells is associated with the organization of the Drosophila genome.

Morphologically, polytene chromosomes of Drosophila melanogaster consist of compact "black" bands alternating with less compact "grey" bands and interbands. We developed a comprehensive approach that combines cytological mapping data of FlyBase-annotated genes and novel tools for predicting cytogenetic features of chromosomes on the basis of their protein composition and determined the genomic coordinates for all black bands of polytene chromosome 2R. By a PCNA immunostaining assay, we obtained the replication timetable for all the bands mapped. The results allowed us to compare replication timing between polytene chromosomes in salivary glands and chromosomes from cultured diploid cell lines and to observe a substantial similarity in the global replication patterns at the band resolution level. In both kinds of chromosomes, the intervals between black bands correspond to early replication initiation zones. Black bands are depleted of replication initiation events and are characterized by a gradient of replication timing; therefore, the time of replication completion correlates with the band length. The bands are characterized by low gene density, contain predominantly tissue-specific genes, and are represented by silent chromatin types in various tissues. The borders of black bands correspond well to the borders of topological domains as well as to the borders of the zones showing H3K27me3, SUUR, and LAMIN enrichment. In conclusion, the characteristic pattern of polytene chromosomes reflects partitioning of the Drosophila genome into two global types of domains with contrasting properties. This partitioning is conserved in different tissues and determines replication timing in Drosophila.

TADs are 3D structural units of higher-order chromosome organization in Drosophila.

Deciphering the rules of genome folding in the cell nucleus is essential to understand its functions. Recent chromosome conformation capture (Hi-C) studies have revealed that the genome is partitioned into topologically associating domains (TADs), which demarcate functional epigenetic domains defined by combinations of specific chromatin marks. However, whether TADs are true physical units in each cell nucleus or whether they reflect statistical frequencies of measured interactions within cell populations is unclear. Using a combination of Hi-C, three-dimensional (3D) fluorescent in situ hybridization, super-resolution microscopy, and polymer modeling, we provide an integrative view of chromatin folding in Drosophila. We observed that repressed TADs form a succession of discrete nanocompartments, interspersed by less condensed active regions. Single-cell analysis revealed a consistent TAD-based physical compartmentalization of the chromatin fiber, with some degree of heterogeneity in intra-TAD conformations and in cis and trans inter-TAD contact events. These results indicate that TADs are fundamental 3D genome units that engage in dynamic higher-order inter-TAD connections. This domain-based architecture is likely to play a major role in regulatory transactions during DNA-dependent processes.

Chromosome-nuclear envelope attachments affect interphase chromosome territories and entanglement.

BACKGROUND: It is well recognized that the interphase chromatin of higher eukaryotes folds into non-random configurations forming territories within the nucleus. Chromosome territories have biologically significant properties, and understanding how these properties change with time during lifetime of the cell is important. Chromosome-nuclear envelope (Chr-NE) interactions play a role in epigenetic regulation of DNA replication, repair, and transcription. However, their role in maintaining chromosome territories remains unclear. RESULTS: We use coarse-grained molecular dynamics simulations to study the effects of Chr-NE interactions on the dynamics of chromosomes within a model of the Drosophila melanogaster regular (non-polytene) interphase nucleus, on timescales comparable to the duration of interphase. The model simulates the dynamics of chromosomes bounded by the NE. Initially, the chromosomes in the model are prearranged in fractal-like configurations with physical parameters such as nucleus size and chromosome persistence length taken directly from experiment. Time evolution of several key observables that characterize the chromosomes is quantified during each simulation: chromosome territories, chromosome entanglement, compactness, and presence of the Rabl (polarized) chromosome arrangement. We find that Chr-NE interactions help maintain chromosome territories by slowing down and limiting, but not eliminating, chromosome entanglement on biologically relevant timescales. At the same time, Chr-NE interactions have little effect on the Rabl chromosome arrangement as well as on how chromosome compactness changes with time. These results are rationalized by simple dimensionality arguments, robust to model details. All results are robust to the simulated activity of topoisomerase, which may be present in the interphase cell nucleus. CONCLUSIONS: Our study demonstrates that Chr-NE attachments may help maintain chromosome territories, while slowing down and limiting chromosome entanglement on biologically relevant timescales. However, Chr-NE attachments have little effect on chromosome compactness or the Rabl chromosome arrangement.