Variable chromatin structure revealed by in situ spatially correlated DNA cleavage mapping.

Chromatin structure at the length scale encompassing local nucleosome-nucleosome interactions is thought to play a crucial role in regulating transcription and access to DNA. However, this secondary structure of chromatin remains poorly understood compared with the primary structure of single nucleosomes or the tertiary structure of long-range looping interactions. Here we report the first genome-wide map of chromatin conformation in human cells at the 1-3 nucleosome (50-500 bp) scale, obtained using ionizing radiation-induced spatially correlated cleavage of DNA with sequencing (RICC-seq) to identify DNA-DNA contacts that are spatially proximal. Unbiased analysis of RICC-seq signal reveals regional enrichment of DNA fragments characteristic of alternating rather than adjacent nucleosome interactions in tri-nucleosome units, particularly in H3K9me3-marked heterochromatin. We infer differences in the likelihood of nucleosome-nucleosome contacts among open chromatin, H3K27me3-marked, and H3K9me3-marked repressed chromatin regions. After calibrating RICC-seq signal to three-dimensional distances, we show that compact two-start helical fibre structures with stacked alternating nucleosomes are consistent with RICC-seq fragmentation patterns from H3K9me3-marked chromatin, while non-compact structures and solenoid structures are consistent with open chromatin. Our data support a model of chromatin architecture in intact interphase nuclei consistent with variable longitudinal compaction of two-start helical fibres.

BAC-FISH Based Physical Map of Endangered Catfish Clarias magur for Chromosome Cataloguing and Gene Isolation through Positional Cloning.

Construction of a physical chromosome map of a species is important for positional cloning, targeted marker development, fine mapping of genes, selection of candidate genes for molecular breeding, as well as understanding the genome organization. The genomic libraries in the form of bacterial artificial chromosome (BAC) clones are also a very useful resource for physical mapping and identification and isolation of full-length genes and the related cis acting elements. Some BAC-FISH based studies reported in the past were gene based physical chromosome maps of Clarias magur (magur) to understand the genome organization of the species and to establish the relationships with other species in respect to genes' organization and evolution in the past. In the present study, we generated end sequences of the BAC clones and analyzed those end sequences within the scaffolds of the draft genome of magur to identify and map the genes bioinformatically for each clone. A total of 36 clones mostly possessing genes were identified and used in probe construction and their subsequent hybridization on the metaphase chromosomes of magur. This study successfully mapped all 36 specific clones on 16 chromosome pairs, out of 25 pairs of magur chromosomes. These clones are now recognized as chromosome-specific makers, which are an aid in individual chromosome identification and fine assembly of the genome sequence, and will ultimately help in developing anchored genes' map on the chromosomes of C. magur for understanding their organization, inheritance of important fishery traits and evolution of magur with respect to channel catfish, zebrafish and other species.

High-resolution mapping of a major and consensus quantitative trait locus for oil content to a ~ 0.8-Mb region on chromosome A08 in peanut (Arachis hypogaea L.).

KEY MESSAGE: ddRAD-seq-based high-density genetic map comprising 2595 loci identified a major and consensus QTL with a linked marker in a 0.8-Mb physical interval for oil content in peanut. Enhancing oil content is an important breeding objective in peanut. High-resolution mapping of quantitative trait loci (QTLs) with linked markers could facilitate marker-assisted selection in breeding for target traits. In the present study, a recombined inbred line population (Xuhua 13 × Zhonghua 6) was used to construct a genetic map based on double-digest restriction-site-associated DNA sequencing (ddRAD-seq). The resulting high-density genetic map contained 2595 loci, and spanned a length of 2465.62 cM, with an average distance of 0.95 cM/locus. Seven QTLs for oil content were identified on five linkage groups, including the major and stable QTL qOCA08.1 on chromosome A08 with 10.14-27.19% phenotypic variation explained. The physical interval of qOCA08.1 was further delimited to a ~ 0.8-Mb genomic region where two genes affecting oil synthesis had been annotated. The marker SNPOCA08 was developed targeting the SNP loci associated with oil content and validated in peanut cultivars with diverse oil contents. The major and stable QTL identified in the present study could be further dissected for gene discovery. Furthermore, the tightly linked marker for oil content would be useful in marker-assisted breeding in peanut.

Improved protocols for BAC insert DNA isolation, BAC end sequencing and FISH for construction of BAC based physical map of genes on the chromosomes.

Bacterial artificial chromosome (BAC) library is an important genomic resource useful in targeted marker development, positional cloning, physical mapping and a substrate for genome sequencing for better understanding the genome organization of a species. The present manuscript elucidates the improvement in protocols for economical and efficient BAC insert DNA isolation, BAC end sequencing and FISH for physical localization on the metaphase chromosome complements. BAC clones of Clarias magur, maintained in 384-well plate format in our laboratory, were used in this study. The protocols gave consistent and efficient results. We use routinely these protocols for BAC insert DNA extraction, generating end sequence data of the clone and constructing DNA probes to hybridize on the metaphase spreads of C. magur using FISH for physical their localization.

DDmap: a MATLAB package for the double digest problem using multiple genetic operators.

BACKGROUND: In computational biology, the physical mapping of DNA is a key problem. We know that the double digest problem (DDP) is NP-complete. Many algorithms have been proposed for solving the DDP, although it is still far from being resolved. RESULTS: We present DDmap, an open-source MATLAB package for solving the DDP, based on a newly designed genetic algorithm that combines six genetic operators in searching for optimal solutions. We test the performance of DDmap by using a typical DDP dataset, and we depict exact solutions to these DDP instances in an explicit manner. In addition, we propose an approximate method for solving some hard DDP scenarios via a scaling-rounding-adjusting process. CONCLUSIONS: For typical DDP test instances, DDmap finds exact solutions within approximately 1 s. Based on our simulations on 1000 random DDP instances by using DDmap, we find that the maximum length of the combining fragments has observable effects towards genetic algorithms for solving the DDP problem. In addition, a Maple source code for illustrating DDP solutions as nested pie charts is also included.

A high-density BAC physical map covering the entire MHC region of addax antelope genome.

BACKGROUND: The mammalian major histocompatibility complex (MHC) harbours clusters of genes associated with the immunological defence of animals against infectious pathogens. At present, no complete MHC physical map is available for any of the wild ruminant species in the world. RESULTS: The high-density physical map is composed of two contigs of 47 overlapping bacterial artificial chromosome (BAC) clones, with an average of 115 Kb for each BAC, covering the entire addax MHC genome. The first contig has 40 overlapping BAC clones covering an approximately 2.9 Mb region of MHC class I, class III, and class IIa, and the second contig has 7 BAC clones covering an approximately 500 Kb genomic region that harbours MHC class IIb. The relative position of each BAC corresponding to the MHC sequence was determined by comparative mapping using PCR screening of the BAC library of 192,000 clones, and the order of BACs was determined by DNA fingerprinting. The overlaps of neighboring BACs were cross-verified by both BAC-end sequencing and co-amplification of identical PCR fragments within the overlapped region, with their identities further confirmed by DNA sequencing. CONCLUSIONS: We report here the successful construction of a high-quality physical map for the addax MHC region using BACs and comparative mapping. The addax MHC physical map we constructed showed one gap of approximately 18 Mb formed by an ancient autosomal inversion that divided the MHC class II into IIa and IIb. The autosomal inversion provides compelling evidence that the MHC organizations in all of the ruminant species are relatively conserved.

Fine mapping of qGL5H, a major grain length locus in barley (Hordeum vulgare L.).

KEY MESSAGE: A major grain length QTL on chromosome 5H was fine mapped from 180.5 to 1.7 Mb. Quantitative trait loci (QTLs) mapping has been used extensively in barley to detect QTLs that underlie complex traits such as grain size. In the present study, we utilised 312 double haploid lines derived from a cross between two Australian malting varieties, Vlamingh and Buloke, to dissect the genetic control of a number of grain size characteristics. Digital image analysis was used to measure grain size characteristics including length, width, thickness and plumpness which are important traits influencing barley yield and grain physical quality. Using data from four independent environments and molecular marker genotype data, we identified 23 significant QTLs for these four traits, ten of which were consensus QTLs and identified in two or more environments. A QTL region on chromosome 5H designated qGL5H that was associated with grain size was fine mapped to a 1.7 Mb interval. qGL5H was able to explain 21.6% of phenotypic variation for grain length within the population. This major QTL is an appropriate candidate for further genetic dissection.

Fine mapping and molecular marker development of the Sm gene conferring resistance to gray leaf spot (Stemphylium spp.) in tomato.

KEY MESSAGE: The tomato gray leaf spot resistance gene Sm was fine-mapped in a 185-kb region through a map-based cloning strategy and genome-wide association study; a candidate gene was proved to be involved in Sm-mediated resistance through transient gene silencing. Gray leaf spot, caused by Stemphylium spp., is a warm weather foliar disease in tomato (Solanum lycopersicum L). Resistance against gray leaf spot is conferred by a single incompletely dominant gene (Sm) located on chromosome 11. This study aimed to map and identify molecular marker tightly linked to the Sm gene for the use of marker-assisted selection in breeding. Using an F2 population derived from a cross between the resistant line '9706' and the susceptible line 'Heinz 1706', the Sm gene was mapped to a 185-kb interval between two markers, InDel343 and InDel-FT-32 on chromosome 11, which was consistent with the result of a genome-wide association study using 289 diverse accessions. An ORF predicted in this region was proved to be involved in Sm-mediated resistance through transient gene silencing and seems to be a good candidate of the Sm locus. To clone the Sm gene, a bacterial artificial chromosome (BAC) library was screened and one BAC clone B80B15 containing the predicted ORF was identified. The analysis of sequence and structure characteristics demonstrated that the candidate gene was not a typical type resistance gene. Additionally, a co-dominant marker Sm-InDel, which produced a 122-bp or 140-bp fragment for resistant or susceptible alleles, respectively, was developed. This marker was validated in 289 germplasm and could be used in marker-assisted selection for gray leaf spot resistance.

pBACode: a random-barcode-based high-throughput approach for BAC paired-end sequencing and physical clone mapping.

Applications that use Bacterial Artificial Chromosome (BAC) libraries often require paired-end sequences and knowledge of the physical location of each clone in plates. To facilitate obtaining this information in high-throughput, we generated pBACode vectors: a pool of BAC cloning vectors, each with a pair of random barcodes flanking its cloning site. In a pBACode BAC library, the BAC ends and their linked barcodes can be sequenced in bulk. Barcode pairs are determined by sequencing the empty pBACode vectors, which allows BAC ends to be paired according to their barcodes. For physical clone mapping, the barcodes are used as unique markers for their linked genomic sequence. After multi-dimensional pooling of BAC clones, the barcodes are sequenced and deconvoluted to locate each clone. We generated a pBACode library of 94,464 clones for the flounder Paralichthys olivaceus and obtained paired-end sequence from 95.4% of the clones. Incorporating BAC paired-ends into the genome preassembly improved its continuity by over 10-fold. Furthermore, we were able to use the barcodes to map the physical locations of each clone in just 50 pools, with up to 11 808 clones per pool. Our physical clone mapping located 90.2% of BAC clones, enabling targeted characterization of chromosomal rearrangements.

Reconstruction of an SSR-based Magnaporthe oryzae physical map to locate avirulence gene AvrPi12.

BACKGROUND: Pathogen avirulence (Avr) genes can evolve rapidly when challenged by the widespread deployment of host genes for resistance. They can be effectively isolated by positional cloning provided a robust and well-populated genetic map is available. RESULTS: An updated, SSR-based physical map of the rice blast pathogen Magnaporthe oryzae (Mo) has been constructed based on 116 of the 120 SSRs used to assemble the last map, along with 18 newly developed ones. A comparison between the two versions of the map has revealed an altered marker content and order within most of the Mo chromosomes. The avirulence gene AvrPi12 was mapped in a population of 219 progeny derived from a cross between the two Mo isolates CHL42 and CHL357. A bulked segregant analysis indicated that the gene was located on chromosome 6, a conclusion borne out by an analysis of the pattern of segregation shown by individual isolates. Six additional PCR-based markers were developed to improve the map resolution in the key region. AvrPi12 was finally located within the sub-telomeric region of chromosome 6, distal to the SSR locus LSM6-5. CONCLUSIONS: The improved SSR-based linkage map should be useful as a platform for gene mapping and isolation in Mo. It was used to establish the location of AvrPi12, thereby providing a starting point for its positional cloning.

Dense gene physical maps of the non-model species Drosophila subobscura.

The comparative analysis of genetic and physical maps as well as of whole genome sequences had revealed that in the Drosophila genus, most structural rearrangements occurred within chromosomal elements as a result of paracentric inversions. Genome sequence comparison would seem the best method to estimate rates of chromosomal evolution, but the high-quality reference genomes required for this endeavor are still scanty. Here, we have obtained dense physical maps for Muller elements A, C, and E of Drosophila subobscura, a species with an extensively studied rich and adaptive chromosomal polymorphism. These maps are based on 462 markers: 115, 236, and 111 markers for elements A, C, and E, respectively. The availability of these dense maps will facilitate genome assembly and will thus greatly contribute to obtaining a good reference genome, which is a required step for D. subobscura to attain the model species status. The comparative analysis of these physical maps and those obtained from the D. pseudoobscura and D. melanogaster genomes allowed us to infer the number of fixed inversions and chromosomal evolutionary rates for each pairwise comparison. For all three elements, rates inferred from the more closely related species were higher than those inferred from the more distantly related species, which together with results of relative-rate tests point to an acceleration in the D. subobscura lineage at least for elements A and E.

High-resolution in situ hybridization analysis on the chromosomal interval 61C7-61C8 of Drosophila melanogaster reveals interbands as open chromatin domains.

Eukaryotic chromatin is organized in contiguous domains that differ in protein binding, histone modifications, transcriptional activity, and in their degree of compaction. Genome-wide comparisons suggest that, overall, the chromatin organization is similar in different cells within an organism. Here, we compare the structure and activity of the 61C7-61C8 interval in polytene and diploid cells of Drosophila. By in situ hybridization on polytene chromosomes combined with high-resolution microscopy, we mapped the boundaries of the 61C7-8 interband and of the 61C7 and C8 band regions, respectively. Our results demonstrate that the 61C7-8 interband is significantly larger than estimated previously. This interband extends over 20 kbp and is in the range of the flanking band domains. It contains several active genes and therefore can be considered as an open chromatin domain. Comparing the 61C7-8 structure of Drosophila S2 cells and polytene salivary gland cells by ChIP for chromatin protein binding and histone modifications, we observe a highly consistent domain structure for the proximal 13 kbp of the domain in both cell types. However, the distal 7 kbp of the open domain differs in protein binding and histone modification between both tissues. The domain contains four protein-coding genes in the proximal part and two noncoding transcripts in the distal part. The differential transcriptional activity of one of the noncoding transcripts correlates with the observed differences in the chromatin structure between both tissues. The significance of our findings for the organization and structure of open chromatin domains will be discussed.

Mapping centromeres of microchromosomes in the zebra finch (Taeniopygia guttata) using half-tetrad analysis.

Centromeres usually consist of hundreds of kilobases of repetitive sequence which renders them difficult to assemble. As a consequence, centromeres are often missing from assembled genomes and their locations on physical chromosome maps have to be inferred from flanking sequences via fluorescence in situ hybridization (FISH). Alternatively, centromere positions can be mapped using linkage analyses in accidentally triploid individuals formed by half-tetrads (resulting from the inheritance of two chromatids from a single meiosis). The current genome assembly of the zebra finch (Taeniopygia guttata) comprises 32 chromosomes, but only for the ten largest chromosomes centromere positions have been mapped using FISH. We here map the positions of most of the remaining centromeres using half-tetrad analyses. For this purpose, we genotyped 37 zebra finches that were triploid or tetraploid due to inheritance errors (and mostly died as embryos) together with their parents at 64 microsatellite markers (at least two per chromosome). Using the information on centromere positions on the ten largest chromosomes, we were able to identify 12 cases of non-disjunction in maternal meiosis I and 10 cases of non-disjunction in maternal meiosis II. These 22 informative cases allowed us to infer centromere positions on additional 19 microchromosomes in reference to the current genome assembly. This knowledge will be valuable for studies of chromosome evolution, meiotic drive and species divergence in the avian lineage.

Correlating the Genetic and Physical Map of Barley Chromosome 3H Revealed Limitations of the FISH-Based Mapping of Nearby Single-Copy Probes Caused by the Dynamic Structure of Metaphase Chromosomes.

Genetic maps are based on the recombination frequency of molecular markers which often show different positions in comparison to the corresponding physical maps. To decipher the position and order of DNA sequences genetically mapped to terminal and interstitial regions of barley (Hordeum vulgare) chromosome 3H, fluorescence in situ hybridization (FISH) on mitotic metaphase chromosomes was performed with 16 genomic single-copy probes derived from fingerprinted BAC contigs. Long genetic distances at subterminal regions translated into short physical distances, confirming that recombination events occur more often at distal regions of chromosome 3H. Nonoverlapping FISH signals were frequently obtained for probes with a physical distance of at least 30-60 kb. Only 8% of the analyzed chromosomes showed a symmetric order of FISH signals on both sister chromatids. Due to the dynamic packing of metaphase chromatin, the order of 2 adjacent single-copy signals along the chromosome arms outside the (peri)centromeric region can only reliably be determined if the cytological distance is approximately 3%, corresponding to 21.6 Mb.

OMWare: a tool for efficient assembly of genome-wide physical maps.

BACKGROUND: Physical mapping of DNA with restriction enzymes allows for the characterization and assembly of much longer molecules than is feasible with sequencing. However, assemblies of physical map data are sensitive to input parameters, which describe noise inherent in the data collection process. One possible way to determine the parameter values that best describe a dataset is by trial and error. RESULTS: Here we present OMWare, a tool that efficiently generated 405 de novo map assemblies of a single datasets collected from the cotton species Gossypium raimondii. The assemblies were generated using various input parameter values, and were completed more efficiently by re-using compatible intermediate results. These assemblies were assayed for contiguity, internal consistency, and accuracy. CONCLUSIONS: Resulting assemblies had variable qualities. Although highly accurate assemblies were found, contiguity and internal consistency metrics were poor predictors of accuracy.

The Arabidopsis TAC Position Viewer: a high-resolution map of transformation-competent artificial chromosome (TAC) clones aligned with the Arabidopsis thaliana Columbia-0 genome.

We present a high-resolution map of genomic transformation-competent artificial chromosome (TAC) clones extending over all Arabidopsis thaliana (Arabidopsis) chromosomes. The Arabidopsis genomic TAC clones have been valuable genetic tools. Previously, we constructed an Arabidopsis genomic TAC library consisting of more than 10,000 TAC clones harboring large genomic DNA fragments extending over the whole Arabidopsis genome. Here, we determined 13,577 end sequences from 6987 Arabidopsis TAC clones and mapped 5937 TAC clones to precise locations, covering approximately 90% of the Arabidopsis chromosomes. We present the large-scale data set of TAC clones with high-resolution mapping information as a Java application tool, the Arabidopsis TAC Position Viewer, which provides ready-to-go transformable genomic DNA clones corresponding to certain loci on Arabidopsis chromosomes. The TAC clone resources will accelerate genomic DNA cloning, positional walking, complementation of mutants and DNA transformation for heterologous gene expression.

Multi-ethnic fine-mapping of 14 central adiposity loci.

The Genetic Investigation of Anthropometric Traits (GIANT) consortium identified 14 loci in European Ancestry (EA) individuals associated with waist-to-hip ratio (WHR) adjusted for body mass index. These loci are wide and narrowing the signals remains necessary. Twelve of 14 loci identified in GIANT EA samples retained strong associations with WHR in our joint EA/individuals of African Ancestry (AA) analysis (log-Bayes factor >6.1). Trans-ethnic analyses at five loci (TBX15-WARS2, LYPLAL1, ADAMTS9, LY86 and ITPR2-SSPN) substantially narrowed the signals to smaller sets of variants, some of which are in regions that have evidence of regulatory activity. By leveraging varying linkage disequilibrium structures across different populations, single-nucleotide polymorphisms (SNPs) with strong signals and narrower credible sets from trans-ethnic meta-analysis of central obesity provide more precise localizations of potential functional variants and suggest a possible regulatory role. Meta-analysis results for WHR were obtained from 77 167 EA participants from GIANT and 23 564 AA participants from the African Ancestry Anthropometry Genetics Consortium. For fine mapping we interrogated SNPs within ± 250 kb flanking regions of 14 previously reported index SNPs from loci discovered in EA populations by performing trans-ethnic meta-analysis of results from the EA and AA meta-analyses. We applied a Bayesian approach that leverages allelic heterogeneity across populations to combine meta-analysis results and aids in fine-mapping shared variants at these locations. We annotated variants using information from the ENCODE Consortium and Roadmap Epigenomics Project to prioritize variants for possible functionality.

SOAPindel: efficient identification of indels from short paired reads.

We present a new approach to indel calling that explicitly exploits that indel differences between a reference and a sequenced sample make the mapping of reads less efficient. We assign all unmapped reads with a mapped partner to their expected genomic positions and then perform extensive de novo assembly on the regions with many unmapped reads to resolve homozygous, heterozygous, and complex indels by exhaustive traversal of the de Bruijn graph. The method is implemented in the software SOAPindel and provides a list of candidate indels with quality scores. We compare SOAPindel to Dindel, Pindel, and GATK on simulated data and find similar or better performance for short indels (<10 bp) and higher sensitivity and specificity for long indels. A validation experiment suggests that SOAPindel has a false-positive rate of ∼10% for long indels (>5 bp), while still providing many more candidate indels than other approaches.

Nucleotide-resolution DNA double-strand break mapping by next-generation sequencing.

We present a genome-wide approach to map DNA double-strand breaks (DSBs) at nucleotide resolution by a method we termed BLESS (direct in situ breaks labeling, enrichment on streptavidin and next-generation sequencing). We validated and tested BLESS using human and mouse cells and different DSBs-inducing agents and sequencing platforms. BLESS was able to detect telomere ends, Sce endonuclease-induced DSBs and complex genome-wide DSB landscapes. As a proof of principle, we characterized the genomic landscape of sensitivity to replication stress in human cells, and we identified >2,000 nonuniformly distributed aphidicolin-sensitive regions (ASRs) overrepresented in genes and enriched in satellite repeats. ASRs were also enriched in regions rearranged in human cancers, with many cancer-associated genes exhibiting high sensitivity to replication stress. Our method is suitable for genome-wide mapping of DSBs in various cells and experimental conditions, with a specificity and resolution unachievable by current techniques.

High-resolution tyramide-FISH mapping of markers tightly linked to the male-fertility restoration (Ms) locus of onion.

KEY MESSAGE: Tyramide FISH was used to locate relatively small genomic amplicons from molecular markers linked to Ms locus onto onion chromosome 2 near the centromere, a region of relatively low recombination. Fluorescence in situ hybridization (FISH) has not been readily exploited for physical mapping of molecular markers in plants due to the technical challenge of visualizing small single-copy probes. Signal amplification using tyramide (tyr) FISH can increase sensitivity up to 100-fold. We used tyr-FISH to physically locate molecular markers tightly linked to the nuclear male-fertility (Ms) restoration locus of onion onto mitotic metaphase, pachytene, and super-stretched pachytene chromosomes. Relatively short genomic amplicons (846-2251 bp) and a cDNA clone (666 bp) were visualized in 9-42 % of observed cells. The markers were assigned to proximal locations close to the centromere on the long arm of chromosome 2, a region of lower recombination, revealing that tightly linked markers may be physically distant from Ms. This result explains why several labs have identified molecular markers tightly linked to the Ms locus after screening relatively few DNA clones or primers and segregating progenies. Although these markers are still useful for marker-aided selection, our results indicate that map-based cloning of Ms will likely be difficult due to reduced recombination near this gene.