Mapping the pericentric heterochromatin by comparative genomic hybridization analysis and chromosome deletions in Drosophila melanogaster

Effects of Mutations in the Drosophila melanogaster Rif1 Gene on the Replication and Underreplication of Pericentromeric Heterochromatin in Salivary Gland Polytene Chromosomes

In Drosophila salivary gland polytene chromosomes, a substantial portion of heterochromatin is underreplicated. The combination of mutations SuUR^ES and Su(var)3-9⁰⁶ results in the polytenization of a substantial fraction of unique and moderately repeated sequences but has almost no effect on satellite DNA replication. The Rap1 interacting factor 1 (Rif) protein is a conserved regulator of replication timing, and in Drosophila, it affects underreplication in polytene chromosomes. We compared the morphology of pericentromeric regions and labeling patterns of in situ hybridization of heterochromatin-specific DNA probes between wild-type salivary gland polytene chromosomes and the chromosomes of Rif1 mutants and SuUR Su(var)3-9⁰⁶ double mutants. We show that, despite general similarities, heterochromatin zones exist that are polytenized only in the Rif1 mutants, and that there are zones that are under specific control of Su(var)3-9. In the Rif1 mutants, we found additional polytenization of the largest blocks of satellite DNA (in particular, satellite 1.688 of chromosome X and simple satellites in chromosomes X and 4) as well as partial polytenization of chromosome Y. Data on pulsed incorporation of 5-ethynyl-2′-deoxyuridine (EdU) into polytene chromosomes indicated that in the Rif1 mutants, just as in the wild type, most of the heterochromatin becomes replicated during the late S phase. Nevertheless, a significantly increased number of heterochromatin replicons was noted. These results suggest that Rif1 regulates the activation probability of heterochromatic origins in the satellite DNA region.

A Gene Expression Screen in Drosophila melanogaster Identifies Novel JAK/STAT and EGFR Targets During Oogenesis

The Janus Kinase/Signal Transducer and Activator of Transcription (JAK/STAT) and epidermal growth factor receptor (EGFR) signaling pathways are conserved regulators of tissue patterning, morphogenesis, and other cell biological processes. During Drosophila oogenesis, these pathways determine the fates of epithelial follicle cells (FCs). JAK/STAT and EGFR together specify a population of cells called the posterior follicle cells (PFCs), which signal to the oocyte to establish the embryonic axes. In this study, whole genome expression analysis was performed to identify genes activated by JAK/STAT and/or EGFR. We observed that 317 genes were transcriptionally upregulated in egg chambers with ectopic JAK/STAT and EGFR activity in the FCs. The list was enriched for genes encoding extracellular matrix (ECM) components and ECM-associated proteins. We tested 69 candidates for a role in axis establishment using RNAi knockdown in the FCs. We report that the signaling protein Semaphorin 1b becomes enriched in the PFCs in response to JAK/STAT and EGFR. We also identified ADAM metallopeptidase with thrombospondin type 1 motif A (AdamTS-A) as a novel target of JAK/STAT in the FCs that regulates egg chamber shape. AdamTS-A mRNA becomes enriched at the anterior and posterior poles of the egg chamber at stages 6 to 7 and is regulated by JAK/STAT. Altering AdamTS-A expression in the poles or middle of the egg chamber produces rounder egg chambers. We propose that AdamTS-A regulates egg shape by remodeling the basement membrane.

A Short History and Description of Drosophila melanogaster Classical Genetics: Chromosome Aberrations, Forward Genetic Screens, and the Nature of Mutations

The purpose of this chapter in FlyBook is to acquaint the reader with the Drosophila genome and the ways in which it can be altered by mutation. Much of what follows will be familiar to the experienced Fly Pusher but hopefully will be useful to those just entering the field and are thus unfamiliar with the genome, the history of how it has been and can be altered, and the consequences of those alterations. I will begin with the structure, content, and organization of the genome, followed by the kinds of structural alterations (karyotypic aberrations), how they affect the behavior of chromosomes in meiotic cell division, and how that behavior can be used. Finally, screens for mutations as they have been performed will be discussed. There are several excellent sources of detailed information on Drosophila husbandry and screening that are recommended for those interested in further expanding their familiarity with Drosophila as a research tool and model organism. These are a book by Ralph Greenspan and a review article by John Roote and Andreas Prokop, which should be required reading for any new student entering a fly lab for the first time.

Single-molecule sequencing resolves the detailed structure of complex satellite DNA loci in Drosophila melanogaster

Highly repetitive satellite DNA (satDNA) repeats are found in most eukaryotic genomes. SatDNAs are rapidly evolving and have roles in genome stability and chromosome segregation. Their repetitive nature poses a challenge for genome assembly and makes progress on the detailed study of satDNA structure difficult. Here, we use single-molecule sequencing long reads from Pacific Biosciences (PacBio) to determine the detailed structure of all major autosomal complex satDNA loci in Drosophila melanogaster, with a particular focus on the 260-bp and Responder satellites. We determine the optimal de novo assembly methods and parameter combinations required to produce a high-quality assembly of these previously unassembled satDNA loci and validate this assembly using molecular and computational approaches. We determined that the computationally intensive PBcR-BLASR assembly pipeline yielded better assemblies than the faster and more efficient pipelines based on the MHAP hashing algorithm, and it is essential to validate assemblies of repetitive loci. The assemblies reveal that satDNA repeats are organized into large arrays interrupted by transposable elements. The repeats in the center of the array tend to be homogenized in sequence, suggesting that gene conversion and unequal crossovers lead to repeat homogenization through concerted evolution, although the degree of unequal crossing over may differ among complex satellite loci. We find evidence for higher-order structure within satDNA arrays that suggest recent structural rearrangements. These assemblies provide a platform for the evolutionary and functional genomics of satDNAs in pericentric heterochromatin.

Regulatory functions and chromatin loading dynamics of linker histone H1 during endoreplication in Drosophila

Eukaryotic DNA replicates asynchronously, with discrete genomic loci replicating during different stages of S phase. Drosophila larval tissues undergo endoreplication without cell division, and the latest replicating regions occasionally fail to complete endoreplication, resulting in underreplicated domains of polytene chromosomes. Here we show that linker histone H1 is required for the underreplication (UR) phenomenon in Drosophila salivary glands. H1 directly interacts with the Suppressor of UR (SUUR) protein and is required for SUUR binding to chromatin in vivo. These observations implicate H1 as a critical factor in the formation of underreplicated regions and an upstream effector of SUUR. We also demonstrate that the localization of H1 in chromatin changes profoundly during the endocycle. At the onset of endocycle S (endo-S) phase, H1 is heavily and specifically loaded into late replicating genomic regions and is then redistributed during the course of endoreplication. Our data suggest that cell cycle-dependent chromosome occupancy of H1 is governed by several independent processes. In addition to the ubiquitous replication-related disassembly and reassembly of chromatin, H1 is deposited into chromatin through a novel pathway that is replication-independent, rapid, and locus-specific. This cell cycle-directed dynamic localization of H1 in chromatin may play an important role in the regulation of DNA replication timing.

Comparative Genomic Analyses Provide New Insights into the Evolutionary Dynamics of Heterochromatin in Drosophila

The term heterochromatin has been long considered synonymous with gene silencing, but it is now clear that the presence of transcribed genes embedded in pericentromeric heterochromatin is a conserved feature in the evolution of eukaryotic genomes. Several studies have addressed the epigenetic changes that enable the expression of genes in pericentric heterochromatin, yet little is known about the evolutionary processes through which this has occurred. By combining genome annotation analysis and high-resolution cytology, we have identified and mapped 53 orthologs of D. melanogaster heterochromatic genes in the genomes of two evolutionarily distant species, D. pseudoobscura and D. virilis. Our results show that the orthologs of the D. melanogaster heterochromatic genes are clustered at three main genomic regions in D. virilis and D. pseudoobscura. In D. virilis, the clusters lie in the middle of euchromatin, while those in D. pseudoobscura are located in the proximal portion of the chromosome arms. Some orthologs map to the corresponding Muller C element in D. pseudoobscura and D. virilis, while others localize on the Muller B element, suggesting that chromosomal rearrangements that have been instrumental in the fusion of two separate elements involved the progenitors of genes currently located in D. melanogaster heterochromatin. These results demonstrate an evolutionary repositioning of gene clusters from ancestral locations in euchromatin to the pericentromeric heterochromatin of descendent D. melanogaster chromosomes. Remarkably, in both D. virilis and D. pseudoobscura the gene clusters show a conserved association with the HP1a protein, one of the most highly evolutionarily conserved epigenetic marks. In light of these results, we suggest a new scenario whereby ancestral HP1-like proteins (and possibly other epigenetic marks) may have contributed to the evolutionary repositioning of gene clusters into heterochromatin.

Flow-dependent myosin recruitment during Drosophila cellularization requires zygotic dunk activity

Actomyosin contractility underlies force generation in morphogenesis ranging from cytokinesis to epithelial extension or invagination. In Drosophila, the cleavage of the syncytial blastoderm is initiated by an actomyosin network at the base of membrane furrows that invaginate from the surface of the embryo. It remains unclear how this network forms and how it affects tissue mechanics. Here, we show that during Drosophila cleavage, myosin recruitment to the cleavage furrows proceeds in temporally distinct phases of tension-driven cortical flow and direct recruitment, regulated by different zygotic genes. We identify the gene dunk, which we show is transiently transcribed when cellularization starts and functions to maintain cortical myosin during the flow phase. The subsequent direct myosin recruitment, however, is Dunk-independent but requires Slam. The Slam-dependent direct recruitment of myosin is sufficient to drive cleavage in the dunk mutant, and the subsequent development of the mutant is normal. In the dunk mutant, cortical myosin loss triggers misdirected flow and disrupts the hexagonal packing of the ingressing furrows. Computer simulation coupled with laser ablation suggests that Dunk-dependent maintenance of cortical myosin enables mechanical tension build-up, thereby providing a mechanism to guide myosin flow and define the hexagonal symmetry of the furrows.

Summary: During Drosophila cellularisation, myosin recruitment to the cleavage furrows proceeds in temporally and mechanistically distinct phases separately regulated by dunk and slam.

Identification of Bari Transposons in 23 Sequenced Drosophila Genomes Reveals Novel Structural Variants, MITEs and Horizontal Transfer

Bari elements are members of the Tc1-mariner superfamily of DNA transposons, originally discovered in Drosophila melanogaster, and subsequently identified in silico in 11 sequenced Drosophila genomes and as experimentally isolated in four non-sequenced Drosophila species. Bari-like elements have been also studied for their mobility both in vivo and in vitro. We analyzed 23 Drosophila genomes and carried out a detailed characterization of the Bari elements identified, including those from the heterochromatic Bari1 cluster in D. melanogaster. We have annotated 401 copies of Bari elements classified either as putatively autonomous or inactive according to the structure of the terminal sequences and the presence of a complete transposase-coding region. Analyses of the integration sites revealed that Bari transposase prefers AT-rich sequences in which the TA target is cleaved and duplicated. Furthermore evaluation of transposon’s co-occurrence near the integration sites of Bari elements showed a non-random distribution of other transposable elements. We also unveil the existence of a putatively autonomous Bari1 variant characterized by two identical long Terminal Inverted Repeats, in D. rhopaloa. In addition, we detected MITEs related to Bari transposons in 9 species. Phylogenetic analyses based on transposase gene and the terminal sequences confirmed that Bari-like elements are distributed into three subfamilies. A few inconsistencies in Bari phylogenetic tree with respect to the Drosophila species tree could be explained by the occurrence of horizontal transfer events as also suggested by the results of dS analyses. This study further clarifies the Bari transposon’s evolutionary dynamics and increases our understanding on the Tc1-mariner elements’ biology.

Hawaiian Drosophila genomes: size variation and evolutionary expansions

This paper reports genome sizes of one Hawaiian Scaptomyza and 16 endemic Hawaiian Drosophila species that include five members of the antopocerus species group, one member of the modified mouthpart group, and ten members of the picture wing clade. Genome size expansions have occurred independently multiple times among Hawaiian Drosophila lineages, and have resulted in an over 2.3-fold range of genome sizes among species, with the largest observed in Drosophila cyrtoloma (1C = 0.41 pg). We find evidence that these repeated genome size expansions were likely driven by the addition of significant amounts of heterochromatin and satellite DNA. For example, our data reveal that the addition of seven heterochromatic chromosome arms to the ancestral haploid karyotype, and a remarkable proportion of ~70 % satellite DNA, account for the greatly expanded size of the D. cyrtoloma genome. Moreover, the genomes of 13/17 Hawaiian picture wing species are composed of substantial proportions (22-70 %) of detectable satellites (all but one of which are AT-rich). Our results suggest that in this tightly knit group of recently evolved species, genomes have expanded, in large part, via evolutionary amplifications of satellite DNA sequences in centric and pericentric domains (especially of the X and dot chromosomes), which have resulted in longer acrocentric chromosomes or metacentrics with an added heterochromatic chromosome arm. We discuss possible evolutionary mechanisms that may have shaped these patterns, including rapid fixation of novel expanded genomes during founder-effect speciation.

The Release 6 reference sequence of the Drosophila melanogaster genome

Drosophila melanogaster plays an important role in molecular, genetic, and genomic studies of heredity, development, metabolism, behavior, and human disease. The initial reference genome sequence reported more than a decade ago had a profound impact on progress in Drosophila research, and improving the accuracy and completeness of this sequence continues to be important to further progress. We previously described improvement of the 117-Mb sequence in the euchromatic portion of the genome and 21 Mb in the heterochromatic portion, using a whole-genome shotgun assembly, BAC physical mapping, and clone-based finishing. Here, we report an improved reference sequence of the single-copy and middle-repetitive regions of the genome, produced using cytogenetic mapping to mitotic and polytene chromosomes, clone-based finishing and BAC fingerprint verification, ordering of scaffolds by alignment to cDNA sequences, incorporation of other map and sequence data, and validation by whole-genome optical restriction mapping. These data substantially improve the accuracy and completeness of the reference sequence and the order and orientation of sequence scaffolds into chromosome arm assemblies. Representation of the Y chromosome and other heterochromatic regions is particularly improved. The new 143.9-Mb reference sequence, designated Release 6, effectively exhausts clone-based technologies for mapping and sequencing. Highly repeat-rich regions, including large satellite blocks and functional elements such as the ribosomal RNA genes and the centromeres, are largely inaccessible to current sequencing and assembly methods and remain poorly represented. Further significant improvements will require sequencing technologies that do not depend on molecular cloning and that produce very long reads.

Simple method for fluorescence DNA in situ hybridization to squashed chromosomes

DNA in situ hybridization (DNA ISH) is a commonly used method for mapping sequences to specific chromosome regions. This approach is particularly effective at mapping highly repetitive sequences to heterochromatic regions, where computational approaches face prohibitive challenges. Here we describe a streamlined protocol for DNA ISH that circumvents formamide washes that are standard steps in other DNA ISH protocols. Our protocol is optimized for hybridization with short single strand DNA probes that carry fluorescent dyes, which effectively mark repetitive DNA sequences within heterochromatic chromosomal regions across a number of different insect tissue types. However, applications may be extended to use with larger probes and visualization of single copy (non-repetitive) DNA sequences. We demonstrate this method by mapping several different repetitive sequences to squashed chromosomes from Drosophila melanogaster neural cells and Nasonia vitripennis spermatocytes. We show hybridization patterns for both small, commercially synthesized probes and for a larger probe for comparison. This procedure uses simple laboratory supplies and reagents, and is ideal for investigators who have little experience with performing DNA ISH.

The organization and evolution of the Responder satellite in species of the Drosophila melanogaster group: dynamic evolution of a target of meiotic drive

Background

Satellite DNA can make up a substantial fraction of eukaryotic genomes and has roles in genome structure and chromosome segregation. The rapid evolution of satellite DNA can contribute to genomic instability and genetic incompatibilities between species. Despite its ubiquity and its contribution to genome evolution, we currently know little about the dynamics of satellite DNA evolution. The Responder (Rsp) satellite DNA family is found in the pericentric heterochromatin of chromosome 2 of Drosophila melanogaster. Rsp is well-known for being the target of Segregation Distorter (SD)— an autosomal meiotic drive system in D. melanogaster. I present an evolutionary genetic analysis of the Rsp family of repeats in D. melanogaster and its closely-related species in the melanogaster group (D. simulans, D. sechellia, D. mauritiana, D. erecta, and D. yakuba) using a combination of available BAC sequences, whole genome shotgun Sanger reads, Illumina short read deep sequencing, and fluorescence in situ hybridization.

Results

I show that Rsp repeats have euchromatic locations throughout the D. melanogaster genome, that Rsp arrays show evidence for concerted evolution, and that Rsp repeats exist outside of D. melanogaster, in the melanogaster group. The repeats in these species are considerably diverged at the sequence level compared to D. melanogaster, and have a strikingly different genomic distribution, even between closely-related sister taxa.

Conclusions

The genomic organization of the Rsp repeat in the D. melanogaster genome is complex—it exists of large blocks of tandem repeats in the heterochromatin and small blocks of tandem repeats in the euchromatin. My discovery of heterochromatic Rsp-like sequences outside of D. melanogaster suggests that SD evolved after its target satellite and that the evolution of the Rsp satellite family is highly dynamic over a short evolutionary time scale (<240,000 years).

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0233-9) contains supplementary material, which is available to authorized users.

On the evolution of Yeti, a Drosophila melanogaster heterochromatin gene

Constitutive heterochromatin is a ubiquitous and still unveiled component of eukaryotic genomes, within which it comprises large portions. Although constitutive heterochromatin is generally considered to be transcriptionally silent, it contains a significant variety of sequences that are expressed, among which about 300 single-copy coding genes have been identified by genetic and genomic analyses in the last decades. Here, we report the results of the evolutionary analysis of Yeti, an essential gene of Drosophila melanogaster located in the deep pericentromeric region of chromosome 2R. By FISH, we showed that Yeti maintains a heterochromatin location in both D. simulans and D. sechellia species, closely related to D. melanogaster, while in the more distant species e.g., D. pseudoobscura and D. virilis, it is found within euchromatin, in the syntenic chromosome Muller C, that corresponds to the 2R arm of D. melanogaster chromosome 2. Thus, over evolutionary time, Yeti has been resident on the same chromosomal element, but it progressively moved closer to the pericentric regions. Moreover, in silico reconstruction of the Yeti gene structure in 19 Drosophila species and in 5 non-drosophilid dipterans shows a rather stable organization during evolution. Accordingly, by PCR analysis and sequencing, we found that the single intron of Yeti does not undergo major intraspecies or interspecies size changes, unlike the introns of other essential Drosophila heterochromatin genes, such as light and Dbp80. This implicates diverse evolutionary forces in shaping the structural organization of genes found within heterochromatin. Finally, the results of d_S - d_N tests show that Yeti is under negative selection both in heterochromatin and euchromatin, and indicate that the change in genomic location did not affected significantly the molecular evolution of the gene. Together, the results of this work contribute to our understanding of the evolutionary dynamics of constitutive heterochromatin in the genomes of higher eukaryotes.

FlyBase: introduction of the Drosophila melanogaster Release 6 reference genome assembly and large-scale migration of genome annotations

Release 6, the latest reference genome assembly of the fruit fly Drosophila melanogaster, was released by the Berkeley Drosophila Genome Project in 2014; it replaces their previous Release 5 genome assembly, which had been the reference genome assembly for over 7 years. With the enormous amount of information now attached to the D. melanogaster genome in public repositories and individual laboratories, the replacement of the previous assembly by the new one is a major event requiring careful migration of annotations and genome-anchored data to the new, improved assembly. In this report, we describe the attributes of the new Release 6 reference genome assembly, the migration of FlyBase genome annotations to this new assembly, how genome features on this new assembly can be viewed in FlyBase (http://flybase.org) and how users can convert coordinates for their own data to the corresponding Release 6 coordinates.

Repetitive centromeric satellite RNA is essential for kinetochore formation and cell division

SAT III RNA binds to the kinetochore component CENP-C and is required for correct assembly and function of the kinetochore at centromeres.

Chromosome segregation requires centromeres on every sister chromatid to correctly form and attach the microtubule spindle during cell division. Even though centromeres are essential for genome stability, the underlying centromeric DNA is highly variable in sequence and evolves quickly. Epigenetic mechanisms are therefore thought to regulate centromeres. Here, we show that the 359-bp repeat satellite III (SAT III), which spans megabases on the X chromosome of Drosophila melanogaster, produces a long noncoding RNA that localizes to centromeric regions of all major chromosomes. Depletion of SAT III RNA causes mitotic defects, not only of the sex chromosome but also in trans of all autosomes. We furthermore find that SAT III RNA binds to the kinetochore component CENP-C, and is required for correct localization of the centromere-defining proteins CENP-A and CENP-C, as well as outer kinetochore proteins. In conclusion, our data reveal that SAT III RNA is an integral part of centromere identity, adding RNA to the complex epigenetic mark at centromeres in flies.

Absence/presence calling in microarray-based CGH experiments with non-model organisms

Structural variations in genomes are commonly studied by (micro)array-based comparative genomic hybridization. The data analysis methods to infer copy number variation in model organisms (human, mouse) are established. In principle, the procedures are based on signal ratios between test and reference samples and the order of the probe targets in the genome. These procedures are less applicable to experiments with non-model organisms, which frequently comprise non-sequenced genomes with an unknown order of probe targets. We therefore present an additional analysis approach, which does not depend on the structural information of a reference genome, and quantifies the presence or absence of a probe target in an unknown genome. The principle is that intensity values of target probes are compared with the intensities of negative-control probes and positive-control probes from a control hybridization, to determine if a probe target is absent or present. In a test, analyzing the genome content of a known bacterial strain: Staphylococcus aureus MRSA252, this approach proved to be successful, demonstrated by receiver operating characteristic area under the curve values larger than 0.9995. We show its usability in various applications, such as comparing genome content and validating next-generation sequencing reads from eukaryotic non-model organisms.

The Hmr and Lhr hybrid incompatibility genes suppress a broad range of heterochromatic repeats

Hybrid incompatibilities (HIs) cause reproductive isolation between species and thus contribute to speciation. Several HI genes encode adaptively evolving proteins that localize to or interact with heterochromatin, suggesting that HIs may result from co-evolution with rapidly evolving heterochromatic DNA. Little is known, however, about the intraspecific function of these HI genes, the specific sequences they interact with, or the evolutionary forces that drive their divergence. The genes Hmr and Lhr genetically interact to cause hybrid lethality between Drosophila melanogaster and D. simulans, yet mutations in both genes are viable. Here, we report that Hmr and Lhr encode proteins that form a heterochromatic complex with Heterochromatin Protein 1 (HP1a). Using RNA-Seq analyses we discovered that Hmr and Lhr are required to repress transcripts from satellite DNAs and many families of transposable elements (TEs). By comparing Hmr and Lhr function between D. melanogaster and D. simulans we identify several satellite DNAs and TEs that are differentially regulated between the species. Hmr and Lhr mutations also cause massive overexpression of telomeric TEs and significant telomere lengthening. Hmr and Lhr therefore regulate three types of heterochromatic sequences that are responsible for the significant differences in genome size and structure between D. melanogaster and D. simulans and have high potential to cause genetic conflicts with host fitness. We further find that many TEs are overexpressed in hybrids but that those specifically mis-expressed in lethal hybrids do not closely correlate with Hmr function. Our results therefore argue that adaptive divergence of heterochromatin proteins in response to repetitive DNAs is an important underlying force driving the evolution of hybrid incompatibility genes, but that hybrid lethality likely results from novel epistatic genetic interactions that are distinct to the hybrid background.

Sister species capable of mating often produce hybrids that are sterile or die during development. This reproductive isolation is caused by incompatibilities between the two sister species' genomes. Some hybrid incompatibilities involve genes that encode rapidly evolving proteins that localize to heterochromatin. Heterochromatin is largely made up of highly repetitive transposable elements and satellite DNAs. It has been hypothesized that rapid changes in heterochromatic DNA drives the changes in these HI genes and thus the evolution of reproductive isolation. In support of this model, we show that two rapidly evolving HI proteins, Lhr and Hmr, which reproductively isolate the fruit fly sister species D. melanogaster and D. simulans, repress transposable elements and satellite DNAs. These proteins also help regulate the length of the atypical Drosophila telomeres, which are themselves made of domesticated transposable elements. Our data suggest that these proteins are part of the adaptive machinery that allows the host to respond to changes and increases in heterochromatin and to maintain the activity of genes located within or adjacent to heterochromatin.

Efficient identification of Y chromosome sequences in the human and Drosophila genomes

Notwithstanding their biological importance, Y chromosomes remain poorly known in most species. A major obstacle to their study is the identification of Y chromosome sequences; due to its high content of repetitive DNA, in most genome projects, the Y chromosome sequence is fragmented into a large number of small, unmapped scaffolds. Identification of Y-linked genes among these fragments has yielded important insights about the origin and evolution of Y chromosomes, but the process is labor intensive, restricting studies to a small number of species. Apart from these fragmentary assemblies, in a few mammalian species, the euchromatic sequence of the Y is essentially complete, owing to painstaking BAC mapping and sequencing. Here we use female short-read sequencing and k-mer comparison to identify Y-linked sequences in two very different genomes, Drosophila virilis and human. Using this method, essentially all D. virilis scaffolds were unambiguously classified as Y-linked or not Y-linked. We found 800 new scaffolds (totaling 8.5 Mbp), and four new genes in the Y chromosome of D. virilis, including JYalpha, a gene involved in hybrid male sterility. Our results also strongly support the preponderance of gene gains over gene losses in the evolution of the Drosophila Y. In the intensively studied human genome, used here as a positive control, we recovered all previously known genes or gene families, plus a small amount (283 kb) of new, unfinished sequence. Hence, this method works in large and complex genomes and can be applied to any species with sex chromosomes.

Something silent this way forms: the functional organization of the repressive nuclear compartment

The repressive compartment of the nucleus is comprised primarily of telomeric and centromeric regions, the silent portion of ribosomal RNA genes, the majority of transposable element repeats, and facultatively repressed genes specific to different cell types. This compartment localizes into three main regions: the peripheral heterochromatin, perinucleolar heterochromatin, and pericentromeric heterochromatin. Both chromatin remodeling proteins and transcription of noncoding RNAs are involved in maintenance of repression in these compartments. Global reorganization of the repressive compartment occurs at each cell division, during early development, and during terminal differentiation. Differential action of chromatin remodeling complexes and boundary element looping activities are involved in mediating these organizational changes. We discuss the evidence that heterochromatin formation and compartmentalization may drive nuclear organization.