Integrated maps of the Drosophila genome: progress and prospects

A physical map of the Drosophila melanogaster genome is being assembled, consisting of ordered overlapping cosmid clones. The map is constructed in steps, separately for each chromosomal division. Gaps in this map are to be bridged with yeast artificial chromosome clones. Hybridization to previously cloned genes and extensive use of in situ hybridization to polytene chromosomes ensure that the cosmid map is firmly anchored to the wealth of available genetic and cytogenetic information. The intention is to make the physical map widely available as part of an overall, integrated genetic resource for the Drosophila research community.     

Polytene chromosomal maps of 11 Drosophila species: the order of genomic scaffolds inferred from genetic and physical maps

The sequencing of the 12 genomes of members of the genus Drosophila was taken as an opportunity to reevaluate the genetic and physical maps for 11 of the species, in part to aid in the mapping of assembled scaffolds. Here, we present an overview of the importance of cytogenetic maps to Drosophila biology and to the concepts of chromosomal evolution. Physical and genetic markers were used to anchor the genome assembly scaffolds to the polytene chromosomal maps for each species. In addition, a computational approach was used to anchor smaller scaffolds on the basis of the analysis of syntenic blocks. We present the chromosomal map data from each of the 11 sequenced non-Drosophila melanogaster species as a series of sections. Each section reviews the history of the polytene chromosome maps for each species, presents the new polytene chromosome maps, and anchors the genomic scaffolds to the cytological maps using genetic and physical markers. The mapping data agree with Muller's idea that the majority of Drosophila genes are syntenic. Despite the conservation of genes within homologous chromosome arms across species, the karyotypes of these species have changed through the fusion of chromosomal arms followed by subsequent rearrangement events.     

A Physical Map of the X Chromosome of Drosophila Melanogaster: Cosmid Contigs and Sequence Tagged Sites

A physical map of the euchromatic X chromosome of Drosophila melanogaster has been constructed by assembling contiguous arrays of cosmids that were selected by screening a library with DNA isolated from microamplified chromosomal divisions. This map, consisting of 893 cosmids, covers ~64% of the euchromatic part of the chromosome. In addition, 568 sequence tagged sites (STS), in aggregate representing 120 kb of sequenced DNA, were derived from selected cosmids. Most of these STSs, spaced at an average distance of ~35 kb along the euchromatic region of the chromosome, represent DNA tags that can be used as entry points to the fruitfly genome. Furthermore, 42 genes have been placed on the physical map, either through the hybridization of specific probes to the cosmids or through the fact that they were represented among the STSs. These provide a link between the physical and the genetic maps of D. melanogaster. Nine novel genes have been tentatively identified in Drosophila on the basis of matches between STS sequences and sequences from other species.     

Isolation of the Drosophila melanogaster dunce chromosomal region and recombinational mapping of dunce sequences with restriction site polymorphisms as genetic markers.

Using the method of chromosomal walking, we have isolated a contiguous region of the Drosophila melanogaster X chromosome which corresponds to salivary gland chromosome bands 3C12 to 3D4. This five-band region contains approximately 100 kilobases of DNA, including those sequences comprising dunce, a gene which functions in memory and cyclic nucleotide metabolism. Genome blots of DNA from flies carrying several different chromosomal aberrations with breakpoints in the region have been probed with the isolated clones to map the breakpoints on the cloned DNA and to delimit dunce sequences. This has localized dunce to a 50-kilobase region. In addition, we have searched this 50-kilobase region for restriction site polymorphisms between X chromosomes from different Drosophila strains by genome blotting experiments, and we have followed the segregation of detected polymorphisms and dunce alleles after meiotic recombination. The data map one dunce mutation between two polymorphisms located 10 to 12 kilobases apart.     

Functional and comparative genomic analysis of the piebald deletion region of mouse chromosome 14

Several developmentally important genomic regions map within the piebald deletion complex on distal mouse chromosome 14. We have combined computational gene prediction and comparative sequence analysis to characterize an approximately 4.3-Mb segment of the piebald region to identify candidate genes for the phenotypes presented by homozygous deletion mice. As a result we have ordered 13 deletion breakpoints, integrated the sequence with markers from a bacterial artificial chromosome (BAC) physical map, and identified 16 known or predicted genes and >1500 conserved sequence elements (CSEs) across the region. The candidate genes identified include Phr1 (formerly Pam) and Spry2, which are mouse homologs of genes required for development in Drosophila melanogaster. Gene content, order, and position are highly conserved between mouse chromosome 14 and the orthologous region of human chromosome 13. Our studies combining computational gene prediction with genetic and comparative genomic analyses provide insight regarding the functional composition and organization of this defined chromosomal region.     

Clustering of multiple specific genes and gene-rich R-bands around SC-35 domains: evidence for local euchromatic neighborhoods

Typically, eukaryotic nuclei contain 10-30 prominent domains (referred to here as SC-35 domains) that are concentrated in mRNA metabolic factors. Here, we show that multiple specific genes cluster around a common SC-35 domain, which contains multiple mRNAs. Nonsyntenic genes are capable of associating with a common domain, but domain "choice" appears random, even for two coordinately expressed genes. Active genes widely separated on different chromosome arms associate with the same domain frequently, assorting randomly into the 3-4 subregions of the chromosome periphery that contact a domain. Most importantly, visualization of six individual chromosome bands showed that large genomic segments ( approximately 5 Mb) have striking differences in organization relative to domains. Certain bands showed extensive contact, often aligning with or encircling an SC-35 domain, whereas others did not. All three gene-rich reverse bands showed this more than the gene-poor Giemsa dark bands, and morphometric analyses demonstrated statistically significant differences. Similarly, late-replicating DNA generally avoids SC-35 domains. These findings suggest a functional rationale for gene clustering in chromosomal bands, which relates to nuclear clustering of genes with SC-35 domains. Rather than random reservoirs of splicing factors, or factors accumulated on an individual highly active gene, we propose a model of SC-35 domains as functional centers for a multitude of clustered genes, forming local euchromatic "neighborhoods."     

Statistical analysis of domains in interacting protein pairs

MOTIVATION: Several methods have recently been developed to analyse large-scale sets of physical interactions between proteins in terms of physical contacts between the constituent domains, often with a view to predicting new pairwise interactions. Our aim is to combine genomic interaction data, in which domain-domain contacts are not explicitly reported, with the domain-level structure of individual proteins, in order to learn about the structure of interacting protein pairs. Our approach is driven by the need to assess the evidence for physical contacts between domains in a statistically rigorous way. RESULTS: We develop a statistical approach that assigns p-values to pairs of domain superfamilies, measuring the strength of evidence within a set of protein interactions that domains from these superfamilies form contacts. A set of p-values is calculated for SCOP superfamily pairs, based on a pooled data set of interactions from yeast. These p-values can be used to predict which domains come into contact in an interacting protein pair. This predictive scheme is tested against protein complexes in the Protein Quaternary Structure (PQS) database, and is used to predict domain-domain contacts within 705 interacting protein pairs taken from our pooled data set.     

Equilibrium folding-unfolding pathways of model proteins: effect of myoglobin-heme contacts

The equilibrium of protein folding–unfolding has been investigated with a simple two-state, native and random-coil, model; we have termed this the globule–coil model. Energies are calculated by favoring native long-range contact pairs. Most-probable native domains are obtained at all stages of the transition; plausible folding pathways are constructed by connecting these domains by assuming simple growth. Even though native heme–protein contacts represent less than 6% of the total number of native contact pairs, their inclusion appears to change the folding pathway of apomyoglobin from the growth and merging of two native domains to the growth of a single domain. This indicates that pathways derived with this method may be critically sensitive to the details of the contact map and physical constraints during the folding process.     

A plant-transformation-competent BIBAC/BAC-based map of rice for functional analysis and genetic engineering of its genomic sequence.

Sequencing of the rice genome has provided a platform for functional genomics research of rice and other cereal species. However, multiple approaches are needed to determine the functions of its genes and sequences and to use the genome sequencing results for genetic improvement of cereal crops. Here, we report a plant-transformation-competent, binary bacterial artificial chromosome (BIBAC) and bacterial artificial chromosome (BAC) based map of rice to facilitate these studies. The map was constructed from 20 835 BIBAC and BAC clones, and consisted of 579 overlapping BIBAC/BAC contigs. To facilitate functional analysis of chromosome 8 genomic sequence and cloning of the genes and QTLs mapped to the chromosome, we anchored the chromosomal contigs to the existing rice genetic maps. The chromosomal map consists of 11 contigs, 59 genetic markers, and 36 sequence tagged sites, spanning a total of ca. 38 Mb in physical length. Comparative analysis between the genetic and physical maps of chromosome 8 showed that there are 3 "hot" and 2 "cold" spots of genetic recombination along the chromosomal arms in addition to the "cold spot" in the centromeric region, suggesting that the sequence component contents of a chromosome may affect its local genetic recombination frequencies. Because of its plant transformability, the BIBAC/BAC map could provide a platform for functional analysis of the rice genome sequence and effective use of the sequencing results for gene and QTL cloning and molecular breeding.     

Physical and genetic mapping of the Rhodobacter sphaeroides 2.4.1 genome: presence of two unique circular chromosomes.

A macrorestriction map representing the complete physical map of the Rhodobacter sphaeroides 2.4.1 chromosomes has been constructed by ordering the chromosomal DNA fragments from total genomic DNA digested with the restriction endonucleases AseI, SpeI, DraI, and SnaBI. Junction fragments and multiple restriction endonuclease digestions of the chromosomal DNAs derived from wild-type and various mutant strains, in conjunction with Southern hybridization analysis, have been used to order all of the chromosomal DNA fragments. Our results indicate that R. sphaeroides 2.4.1 carries two different circular chromosomes of 3,046 +/- 95 and 914 +/- 17 kilobases (kb). Both chromosome I (3,046 kb) and chromosome II (914 kb) contain rRNA cistrons. It appears that only a single copy of the rRNA genes is contained on chromosome I (rrnA) and that two copies are present on chromosome II (rrnB, rrnC). Additionally, genes for glyceraldehyde 3-phosphate dehydrogenase (gapB) and delta-aminolevulinic acid synthase (hemT) are found on chromosome II. In each instance, there appears to be a second copy of each of these genes on chromosome I, but the extent of the DNA homology is very low. Genes giving rise to enzymes involved in CO2 fixation and linked to the gene encoding the form I enzyme (i.e., the form I region) are on chromosome I, whereas those genes representing the form II region are on chromosome II. The complete physical and partial genetic maps for each chromosome are presented.     

Docking protein domains in contact space

Background  

Many biological processes involve the physical interaction between protein domains. Understanding these functional associations requires knowledge of the molecular structure. Experimental investigations though present considerable difficulties and there is therefore a need for accurate and reliable computational methods. In this paper we present a novel method that seeks to dock protein domains using a contact map representation. Rather than providing a full three dimensional model of the complex, the method predicts contacting residues across the interface. We use a scoring function that combines structural, physicochemical and evolutionary information, where each potential residue contact is assigned a value according to the scoring function and the hypothesis is that the real configuration of contacts is the one that maximizes the score. The search is performed with a simulated annealing algorithm directly in contact space.     

Models that include supercoiling of topological domains reproduce several known features of interphase chromosomes

Understanding the structure of interphase chromosomes is essential to elucidate regulatory mechanisms of gene expression. During recent years, high-throughput DNA sequencing expanded the power of chromosome conformation capture (3C) methods that provide information about reciprocal spatial proximity of chromosomal loci. Since 2012, it is known that entire chromatin in interphase chromosomes is organized into regions with strongly increased frequency of internal contacts. These regions, with the average size of ∼1 Mb, were named topological domains. More recent studies demonstrated presence of unconstrained supercoiling in interphase chromosomes. Using Brownian dynamics simulations, we show here that by including supercoiling into models of topological domains one can reproduce and thus provide possible explanations of several experimentally observed characteristics of interphase chromosomes, such as their complex contact maps.     

Physical Origin of the Contact Frequency in Chromosome Conformation Capture Data

Physical proximity between each pair of genomic loci in a nucleus is measured as a form of contact frequency in chromosome conformation capture-based methods. Complexity of chromosome structure in interphase can be characterized by measuring a statistical property of physical distance between genomic loci according to genomic separation along single chromatids. To find a relationship between the physical distance and the contact frequency, we propose a polymer model derived from the Langevin equation. The model is derived by considering a structure of a chromosome as a trajectory of a particle, where each consecutive segment in the chromosome corresponds to a transient position in the trajectory over time. Using chromosome conformation capture data, we demonstrate the functional relationship between the two quantities. The physical distances derived from the mean contact frequencies by the model show a good correlation with those from experimental data. From the model, we present that the mean contact frequency curve can be divided into three components that arise from different physical origins and show that the contact frequency is proportional to the contact surface area, not to the volume of segments suggested by the fractal globule model. The model explains both a decaying pattern of the contact frequency and the biphasic relationship between the physical distance and the genomic length.     

Physical Origin of the Contact Frequency in Chromosome Conformation Capture Data

Physical proximity between each pair of genomic loci in a nucleus is measured as a form of contact frequency in chromosome conformation capture-based methods. Complexity of chromosome structure in interphase can be characterized by measuring a statistical property of physical distance between genomic loci according to genomic separation along single chromatids. To find a relationship between the physical distance and the contact frequency, we propose a polymer model derived from the Langevin equation. The model is derived by considering a structure of a chromosome as a trajectory of a particle, where each consecutive segment in the chromosome corresponds to a transient position in the trajectory over time. Using chromosome conformation capture data, we demonstrate the functional relationship between the two quantities. The physical distances derived from the mean contact frequencies by the model show a good correlation with those from experimental data. From the model, we present that the mean contact frequency curve can be divided into three components that arise from different physical origins and show that the contact frequency is proportional to the contact surface area, not to the volume of segments suggested by the fractal globule model. The model explains both a decaying pattern of the contact frequency and the biphasic relationship between the physical distance and the genomic length.