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.

Genome-wide characterization and expression profiling of the PDR gene family in tobacco (Nicotiana tabacum).

The pleiotropic drug resistance (PDR) proteins of the ATP-binding cassette (ABC) family play essential roles in physiological processes and have been characterized in many plant species. However, no comprehensive investigation of tobacco (Nicotiana tabacum), an important economic crop and a useful model plant for scientific research, has been presented. We identified 32 PDR genes in the tobacco genome and explored their domain organization, chromosomal distribution and evolution, promoter cis-elements, and expression profiles. A phylogenetic analysis revealed that tobacco has a significantly expanded number of PDR genes involved in plant defense. It also revealed that two tobacco PDR proteins may function as strigolactone transporters to regulate shoot branching, and several NtPDR genes may be involved in cadmium transport. Moreover, tissue expression profiles of NtPDR genes and their responses to several hormones and abiotic stresses were assessed using quantitative real-time PCR. Most of the NtPDR genes were regulated by jasmonate or salicylic acid, suggesting the important regulatory roles of NtPDRs in plant defense and secondary metabolism. They were also responsive to abiotic stresses, like drought and cold, and there was a strong correlation between the presence of promoter cis-elements and abiotic/biotic stress responses. These results provide useful clues for further in-depth studies on the functions of the tobacco PDR genes.

Karyotypic Evolution of Sauropsid Vertebrates Illuminated by Optical and Physical Mapping of the Painted Turtle and Slider Turtle Genomes.

Recent sequencing and software enhancements have advanced our understanding of the evolution of genomic structure and function, especially addressing novel evolutionary biology questions. Yet fragmentary turtle genome assemblies remain a challenge to fully decipher the genetic architecture of adaptive evolution. Here, we use optical mapping to improve the contiguity of the painted turtle (Chrysemys picta) genome assembly and use de novo fluorescent in situ hybridization (FISH) of bacterial artificial chromosome (BAC) clones, BAC-FISH, to physically map the genomes of the painted and slider turtles (Trachemys scripta elegans). Optical mapping increased C. picta's N50 by ~242% compared to the previous assembly. Physical mapping permitted anchoring ~45% of the genome assembly, spanning 5544 genes (including 20 genes related to the sex determination network of turtles and vertebrates). BAC-FISH data revealed assembly errors in C. picta and T. s. elegans assemblies, highlighting the importance of molecular cytogenetic data to complement bioinformatic approaches. We also compared C. picta's anchored scaffolds to the genomes of other chelonians, chicken, lizards, and snake. Results revealed a mostly one-to-one correspondence between chromosomes of painted and slider turtles, and high homology among large syntenic blocks shared with other turtles and sauropsids. Yet, numerous chromosomal rearrangements were also evident across chelonians, between turtles and squamates, and between avian and non-avian reptiles.

Gapless assembly of maize chromosomes using long-read technologies.

Creating gapless telomere-to-telomere assemblies of complex genomes is one of the ultimate challenges in genomics. We use two independent assemblies and an optical map-based merging pipeline to produce a maize genome (B73-Ab10) composed of 63 contigs and a contig N50 of 162 Mb. This genome includes gapless assemblies of chromosome 3 (236 Mb) and chromosome 9 (162 Mb), and 53 Mb of the Ab10 meiotic drive haplotype. The data also reveal the internal structure of seven centromeres and five heterochromatic knobs, showing that the major tandem repeat arrays (CentC, knob180, and TR-1) are discontinuous and frequently interspersed with retroelements.

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.

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.

Physical Map of FISH 5S rDNA and (AG3T3)3 Signals Displays Chimonanthus campanulatus R.H. Chang & C.S. Ding Chromosomes, Reproduces its Metaphase Dynamics and Distinguishes Its Chromosomes.

Chimonanthus campanulatus R.H. Chang & C.S. Ding is a good horticultural tree because of its beautiful yellow flowers and evergreen leaves. In this study, fluorescence in situ hybridization (FISH) was used to analyse mitotic metaphase chromosomes of Ch. campanulatus with 5S rDNA and (AG3T3)3 oligonucleotides. Twenty-two small chromosomes were observed. Weak 5S rDNA signals were observed only in proximal regions of two chromosomes, which were adjacent to the (AG3T3)3 proximal signals. Weak (AG3T3)3 signals were observed on both chromosome ends, which enabled accurate chromosome counts. A pair of satellite bodies was observed. (AG3T3)3 signals displayed quite high diversity, changing in intensity from weak to very strong as follows: far away from the chromosome ends (satellites), ends, subtelomeric regions, and proximal regions. Ten high-quality spreads revealed metaphase dynamics from the beginning to the end and the transition to anaphase. Chromosomes gradually grew larger and thicker into linked chromatids, which grew more significantly in width than in length. Based on the combination of 5S rDNA and (AG3T3)3 signal patterns, ten chromosomes were exclusively distinguished, and the remaining twelve chromosomes were divided into two distinct groups. Our physical map, which can reproduce dynamic metaphase progression and distinguish chromosomes, will powerfully guide cytogenetic research on Chimonanthus and other trees.

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.

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.

Polymer coil-globule phase transition is a universal folding principle of Drosophila epigenetic domains.

BACKGROUND: Localized functional domains within chromosomes, known as topologically associating domains (TADs), have been recently highlighted. In Drosophila, TADs are biochemically defined by epigenetic marks, this suggesting that the 3D arrangement may be the "missing link" between epigenetics and gene activity. Recent observations (Boettiger et al. in Nature 529(7586):418-422, 2016) provide access to structural features of these domains with unprecedented resolution thanks to super-resolution experiments. In particular, they give access to the distribution of the radii of gyration for domains of different linear length and associated with different transcriptional activity states: active, inactive or repressed. Intriguingly, the observed scaling laws lack consistent interpretation in polymer physics. RESULTS: We develop a new methodology conceived to extract the best information from such super-resolution data by exploiting the whole distribution of gyration radii, and to place these experimental results on a theoretical framework. We show that the experimental data are compatible with the finite-size behavior of a self-attracting polymer. The same generic polymer model leads to quantitative differences between active, inactive and repressed domains. Active domains behave as pure polymer coils, while inactive and repressed domains both lie at the coil-globule crossover. For the first time, the "color-specificity" of both the persistence length and the mean interaction energy are estimated, leading to important differences between epigenetic states. CONCLUSION: These results point toward a crucial role of criticality to enhance the system responsivity, resulting in both energy transitions and structural rearrangements. We get strong indications that epigenetically induced changes in nucleosome-nucleosome interaction can cause chromatin to shift between different activity states.

Integrated physical map of bread wheat chromosome arm 7DS to facilitate gene cloning and comparative studies.

Bread wheat (Triticum aestivum L.) is a staple food for a significant part of the world's population. The growing demand on its production can be satisfied by improving yield and resistance to biotic and abiotic stress. Knowledge of the genome sequence would aid in discovering genes and QTLs underlying these traits and provide a basis for genomics-assisted breeding. Physical maps and BAC clones associated with them have been valuable resources from which to generate a reference genome of bread wheat and to assist map-based gene cloning. As a part of a joint effort coordinated by the International Wheat Genome Sequencing Consortium, we have constructed a BAC-based physical map of bread wheat chromosome arm 7DS consisting of 895 contigs and covering 94% of its estimated length. By anchoring BAC contigs to one radiation hybrid map and three high resolution genetic maps, we assigned 73% of the assembly to a distinct genomic position. This map integration, interconnecting a total of 1713 markers with ordered and sequenced BAC clones from a minimal tiling path, provides a tool to speed up gene cloning in wheat. The process of physical map assembly included the integration of the 7DS physical map with a whole-genome physical map of Aegilops tauschii and a 7DS Bionano genome map, which together enabled efficient scaffolding of physical-map contigs, even in the non-recombining region of the genetic centromere. Moreover, this approach facilitated a comparison of bread wheat and its ancestor at BAC-contig level and revealed a reconstructed region in the 7DS pericentromere.

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.

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.

Physical Genome Mapping Using Fluorescence In Situ Hybridization with Mosquito Chromosomes.

The development of genomic resources and tools is an important step in designing novel approaches to genetic control of mosquitoes. Physical genome maps enhance the quality of the genome assemblies, improve gene annotation, and provide a better framework for comparative and population genomics studies in mosquitoes. In this chapter, we describe protocols for an important procedure in physical genome mapping-fluorescence in situ hybridization (FISH). We provide details on (1) dissection of salivary glands, ovaries, and imaginal discs for obtaining high-quality polytene or mitotic chromosome preparations; (2) DNA-labeling procedures and extraction of repetitive DNA fractions; and (3) approaches to FISH on polytene and mitotic chromosomes.

A chromosome-scale assembly of the sorghum genome using nanopore sequencing and optical mapping.

Long-read sequencing technologies have greatly facilitated assemblies of large eukaryotic genomes. In this paper, Oxford Nanopore sequences generated on a MinION sequencer are combined with Bionano Genomics Direct Label and Stain (DLS) optical maps to generate a chromosome-scale de novo assembly of the repeat-rich Sorghum bicolor Tx430 genome. The final assembly consists of 29 scaffolds, encompassing in most cases entire chromosome arms. It has a scaffold N50 of 33.28 Mbps and covers 90% of the expected genome length. A sequence accuracy of 99.85% is obtained after aligning the assembly against Illumina Tx430 data and 99.6% of the 34,211 public gene models align to the assembly. Comparisons of Tx430 and BTx623 DLS maps against the public BTx623 v3.0.1 genome assembly suggest substantial discrepancies whose origin remains to be determined. In summary, this study demonstrates that informative assemblies of complex plant genomes can be generated by combining nanopore sequencing with DLS optical maps.

Challenges and guidelines toward 4D nucleome data and model standards.

Due to recent advances in experimental and theoretical approaches, the dynamic three-dimensional organization (3D) of the nucleus has become a very active area of research in life sciences. We now understand that the linear genome is folded in ways that may modulate how genes are expressed during the basic functioning of cells. Importantly, it is now possible to build 3D models of how the genome folds within the nucleus and changes over time (4D). Because genome folding influences its function, this opens exciting new possibilities to broaden our understanding of the mechanisms that determine cell fate. However, the rapid evolution of methods and the increasing complexity of data can result in ambiguity and reproducibility challenges, which may hamper the progress of this field. Here, we describe such challenges ahead and provide guidelines to think about strategies for shared standardized validation of experimental 4D nucleome data sets and models.

Optical and physical mapping with local finishing enables megabase-scale resolution of agronomically important regions in the wheat genome.

BACKGROUND: Numerous scaffold-level sequences for wheat are now being released and, in this context, we report on a strategy for improving the overall assembly to a level comparable to that of the human genome. RESULTS: Using chromosome 7A of wheat as a model, sequence-finished megabase-scale sections of this chromosome were established by combining a new independent assembly using a bacterial artificial chromosome (BAC)-based physical map, BAC pool paired-end sequencing, chromosome-arm-specific mate-pair sequencing and Bionano optical mapping with the International Wheat Genome Sequencing Consortium RefSeq v1.0 sequence and its underlying raw data. The combined assembly results in 18 super-scaffolds across the chromosome. The value of finished genome regions is demonstrated for two approximately 2.5 Mb regions associated with yield and the grain quality phenotype of fructan carbohydrate grain levels. In addition, the 50 Mb centromere region analysis incorporates cytological data highlighting the importance of non-sequence data in the assembly of this complex genome region. CONCLUSIONS: Sufficient genome sequence information is shown to now be available for the wheat community to produce sequence-finished releases of each chromosome of the reference genome. The high-level completion identified that an array of seven fructosyl transferase genes underpins grain quality and that yield attributes are affected by five F-box-only-protein-ubiquitin ligase domain and four root-specific lipid transfer domain genes. The completed sequence also includes the centromere.

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.

Development of a quantitative pachytene chromosome map and its unification with somatic chromosome and linkage maps of rice (Oryza sativa L.).

A quantitative pachytene chromosome map of rice (Oryza sativa L.) was developed using imaging methods. The map depicts not only distribution patterns of chromomeres specific to pachytene chromosomes, but also the higher order information of chromosomal structures, such as heterochromatin (condensed regions), euchromatin (decondensed regions), the primary constrictions (centromeres), and the secondary constriction (nucleolar organizing regions, NOR). These features were image analyzed and quantitatively mapped onto the map by Chromosome Image Analyzing System ver. 4.0 (CHIAS IV). Correlation between H3K9me2, an epigenetic marker and formation and/or maintenance of heterochromatin, thus was, clearly visualized. Then the pachytene chromosome map was unified with the existing somatic chromosome and linkage maps by physically mapping common DNA markers among them, such as a rice A genome specific tandem repeat sequence (TrsA), 5S and 45S ribosomal RNA genes, five bacterial artificial chromosome (BAC) clones, four P1 bacteriophage artificial chromosome (PAC) clones using multicolor fluorescence in situ hybridization (FISH). Detailed comparison between the locations of the DNA probes on the pachytene chromosomes using multicolor FISH, and the linkage map enabled determination of the chromosome number and short/long arms of individual pachytene chromosomes using the chromosome number and arm assignment designated for the linkage map. As a result, the quantitative pachytene chromosome map was unified with two other major rice chromosome maps representing somatic prometaphase chromosomes and genetic linkages. In conclusion, the unification of the three rice maps serves as an indispensable basic information, not only for an in-depth comparison between genetic and chromosomal data, but also for practical breeding programs.

Investigation of the spatial structure and interactions of the genome at sub-kilobase-pair resolution using T2C.

Chromosome conformation capture (3C) and its derivatives (e.g., 4C, 5C and Hi-C) are used to analyze the 3D organization of genomes. We recently developed targeted chromatin capture (T2C), an inexpensive method for studying the 3D organization of genomes, interactomes and structural changes associated with gene regulation, the cell cycle, and cell survival and development. Here, we present the protocol for T2C based on capture, describing all experimental steps and bio-informatic tools in full detail. T2C offers high resolution, a large dynamic interaction frequency range and a high signal-to-noise ratio. Its resolution is determined by the resulting fragment size of the chosen restriction enzyme, which can lead to sub-kilobase-pair resolution. T2C's high coverage allows the identification of the interactome of each individual DNA fragment, which makes binning of reads (often used in other methods) basically unnecessary. Notably, T2C requires low sequencing efforts. T2C also allows multiplexing of samples for the direct comparison of multiple samples. It can be used to study topologically associating domains (TADs), determining their position, shape, boundaries, and intra- and inter-domain interactions, as well as the composition of aggregated loops, interactions between nucleosomes, individual transcription factor binding sites, and promoters and enhancers. T2C can be performed by any investigator with basic skills in molecular biology techniques in ∼7-8 d. Data analysis requires basic expertise in bioinformatics and in Linux and Python environments.