The development of chromosome microdissection and microcloning technique and its applications in genomic research

Autoploid origin and rapid diploidization of the tetraploid Thinopyrum elongatum revealed by genome differentiation and chromosome pairing in meiosis

Polyploidy is a common mode of evolution in flowering plants. Both the natural tetraploid Thinopyrum elongatum and the diploid one from the same population show a diploid-like pairing in meiosis. However, debate on the chromosome composition and origin of the tetraploid Th. elongatum is ongoing. In the present study, we obtained the induced tetraploid Th. elongatum and found that the induced and natural tetraploids are morphologically close, except for slower development and lower seed setting. Using probes developed from single chromosome microdissection and a Fosmid library, obvious differentiations were discovered between two chromosome sets (E₁ and E₂ ) of the natural tetraploid Th. elongatum but not the induced one. Interestingly, hybrid F₁ derived from the two different wheat-tetraploid Th. elongatum amphiploids 8802 and 8803 produced seeds well. More importantly, analysis of meiosis in F₂ individuals revealed that chromosomes from E₁ and E₂ could pair well on the durum wheat background with the presence of Ph1. No chromosome set differentiation on the FISH level was discovered from the S1 to S4 generations in the induced one. In metaphase of the meiosis first division in the natural tetraploid, more pairings were bivalents and fewer quadrivalents with ratio of 13.94 II + 0.03 IV (n = 31). Chromosome pairing configuration in the induced tetraploid is 13.05 II + 0.47 IV (n = 19), with the quadrivalent ratio being only slightly higher than the ratio in the natural tetraploid. Therefore, the natural tetraploid Th. elongatum is of autoploid origin and the induced tetraploid Th. elongatum evolutionarily underwent rapid diploidization in the low generation.

Towards accurate, contiguous and complete alignment-based polyploid phasing algorithms

Phasing, and in particular polyploid phasing, have been challenging problems held back by the limited read length of high-throughput short read sequencing methods which can't overcome the distance between heterozygous sites and labor high cost of alternative methods such as the physical separation of chromosomes for example. Recently developed single molecule long-read sequencing methods provide much longer reads which overcome this previous limitation. Here we review the alignment-based methods of polyploid phasing that rely on four main strategies: population inference methods, which leverage the genetic information of several individuals to phase a sample; objective function minimization methods, which minimize a function such as the Minimum Error Correction (MEC); graph partitioning methods, which represent the read data as a graph and split it into k haplotype subgraphs; cluster building methods, which iteratively grow clusters of similar reads into a final set of clusters that represent the haplotypes. We discuss the advantages and limitations of these methods and the metrics used to assess their performance, proposing that accuracy and contiguity are the most meaningful metrics. Finally, we propose the field of alignment-based polyploid phasing would greatly benefit from the use of a well-designed benchmarking dataset with appropriate evaluation metrics. We consider that there are still significant improvements which can be achieved to obtain more accurate and contiguous polyploid phasing results which reflect the complexity of polyploid genome architectures.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

nPhase: an accurate and contiguous phasing method for polyploids

While genome sequencing and assembly are now routine, we do not have a full, precise picture of polyploid genomes. No existing polyploid phasing method provides accurate and contiguous haplotype predictions. We developed nPhase, a ploidy agnostic tool that leverages long reads and accurate short reads to solve alignment-based phasing for samples of unspecified ploidy (https://github.com/OmarOakheart/nPhase). nPhase is validated by tests on simulated and real polyploids. nPhase obtains on average over 95% accuracy and a contiguous 1.25 haplotigs per haplotype to cover more than 90% of each chromosome (heterozygosity rate ≥ 0.5%). nPhase allows population genomics and hybrid studies of polyploids.

Cloning and physical localization of male-biased repetitive DNA sequences in Spinacia oleracea (Amaranthaceae)

Spinach (Spinacia oleracea Linnaeus, 1753) is an ideal material for studying molecular mechanisms of early-stage sex chromosome evolution in dioecious plants. Degenerate oligonucleotide-primed polymerase chain reaction (DOP-PCR) technique facilitates the retrotransposon-relevant studies by enriching specific repetitive DNA sequences from a micro-dissected single chromosome. We conducted genomic subtractive hybridization to screen sex-biased DNA sequences by using the DOP-PCR amplification products of micro-dissected spinach Y chromosome. The screening yielded 55 male-biased DNA sequences with 30 576 bp in length, of which, 32 DNA sequences (12 049 bp) contained repeat DNA sequences, including LTR/Copia, LTR/Gypsy, simple repeats, and DNA/CMC-EnSpm. Among these repetitive DNA sequences, four DNA sequences that contained a fragment of Ty3-gypsy retrotransposons (SP73, SP75, SP76, and SP77) were selected as fluorescence probes to hybridization on male and female spinach karyotypes. Fluorescence in situ hybridization (FISH) signals of SP73 and SP75 were captured mostly on the centromeres and their surrounding area for each homolog. Hybridization signals primarily appeared near the putative centromeres for each homologous chromosome pair by using SP76 and SP77 probes for FISH, and sporadic signals existed on the long arms. Results can be served as a basis to study the function of repetitive DNA sequences in sex chromosome evolution in spinach.

Bridging the Gap between Vertebrate Cytogenetics and Genomics with Single-Chromosome Sequencing (ChromSeq)

The study of vertebrate genome evolution is currently facing a revolution, brought about by next generation sequencing technologies that allow researchers to produce nearly complete and error-free genome assemblies. Novel approaches however do not always provide a direct link with information on vertebrate genome evolution gained from cytogenetic approaches. It is useful to preserve and link cytogenetic data with novel genomic discoveries. Sequencing of DNA from single isolated chromosomes (ChromSeq) is an elegant approach to determine the chromosome content and assign genome assemblies to chromosomes, thus bridging the gap between cytogenetics and genomics. The aim of this paper is to describe how ChromSeq can support the study of vertebrate genome evolution and how it can help link cytogenetic and genomic data. We show key examples of ChromSeq application in the refinement of vertebrate genome assemblies and in the study of vertebrate chromosome and karyotype evolution. We also provide a general overview of the approach and a concrete example of genome refinement using this method in the species Anolis carolinensis.

Microdissection and whole chromosome painting confirm karyotype transformation in cryptic species of the Lariophagus distinguendus (Förster, 1841) complex (Hymenoptera: Pteromalidae)

Karyotypes of two cryptic species of parasitoid Hymenoptera with n = 5 and 6 belonging to the Lariophagus distinguendus (Förster, 1841) complex, which includes cosmopolitan parasitoids of coleopteran stored-product pests, were studied using glass-needle based microdissection, reverse and cross-species fluorescence in situ hybridisation (FISH). This experiment strongly indicates that the largest metacentric chromosome in the karyotype with n = 5 originated from a particular fusion between the only acrocentric and a smaller metacentric chromosome of the set with n = 6, therefore confirming our previous hypothesis based on the karyotypic analysis using chromosome morphometrics. This study represents the first successful application of both microdissection and whole chromosome painting for the reconstruction of karyotypic rearrangements in closely related species of parasitoids, as well as in the order Hymenoptera in general.

RecoverY: k-mer-based read classification for Y-chromosome-specific sequencing and assembly

Motivation

The haploid mammalian Y chromosome is usually under-represented in genome assemblies due to high repeat content and low depth due to its haploid nature. One strategy to ameliorate the low coverage of Y sequences is to experimentally enrich Y-specific material before assembly. As the enrichment process is imperfect, algorithms are needed to identify putative Y-specific reads prior to downstream assembly. A strategy that uses k-mer abundances to identify such reads was used to assemble the gorilla Y. However, the strategy required the manual setting of key parameters, a time-consuming process leading to sub-optimal assemblies.

Results

We develop a method, RecoverY, that selects Y-specific reads by automatically choosing the abundance level at which a k-mer is deemed to originate from the Y. This algorithm uses prior knowledge about the Y chromosome of a related species or known Y transcript sequences. We evaluate RecoverY on both simulated and real data, for human and gorilla, and investigate its robustness to important parameters. We show that RecoverY leads to a vastly superior assembly compared to alternate strategies of filtering the reads or contigs. Compared to the preliminary strategy used by Tomaszkiewicz et al., we achieve a 33% improvement in assembly size and a 20% improvement in the NG50, demonstrating the power of automatic parameter selection.

Availability and implementation

Our tool RecoverY is freely available at https://github.com/makovalab-psu/RecoverY.

Contact

kmakova@bx.psu.edu or pashadag@cse.psu.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Protocol for chromosome-specific probe construction using PRINS, micromanipulation and DOP-PCR techniques

Chromosome-specific probes have been widely used in molecular cytogenetics, being obtained with different methods. In this study, a reproducible protocol for construction of chromosome-specific probes is proposed which associates in situ amplification (PRINS), micromanipulation and degenerate oligonucleotide-primed PCR (DOP-PCR). Human lymphocyte cultures were used to obtain metaphases from male and female individuals. The chromosomes were amplified via PRINS, and subcentromeric fragments of the X chromosome were microdissected using microneedles coupled to a phase contrast microscope. The fragments were amplified by DOP-PCR and labeled with tetramethyl-rhodamine-5-dUTP. The probes were used in fluorescent in situ hybridization (FISH) procedure to highlight these specific regions in the metaphases. The results show one fluorescent red spot in male and two in female X chromosomes and interphase nuclei.

Sequence analysis of cultivated strawberry (Fragaria × ananassa Duch.) using microdissected single somatic chromosomes

BACKGROUND

Cultivated strawberry (Fragaria × ananassa Duch.) has homoeologous chromosomes because of allo-octoploidy. For example, two homoeologous chromosomes that belong to different sub-genome of allopolyploids have similar base sequences. Thus, when conducting de novo assembly of DNA sequences, it is difficult to determine whether these sequences are derived from the same chromosome. To avoid the difficulties associated with homoeologous chromosomes and demonstrate the possibility of sequencing allopolyploids using single chromosomes, we conducted sequence analysis using microdissected single somatic chromosomes of cultivated strawberry.

RESULTS

Three hundred and ten somatic chromosomes of the Japanese octoploid strawberry 'Reiko' were individually selected under a light microscope using a microdissection system. DNA from 288 of the dissected chromosomes was successfully amplified using a DNA amplification kit. Using next-generation sequencing, we decoded the base sequences of the amplified DNA segments, and on the basis of mapping, we identified DNA sequences from 144 samples that were best matched to the reference genomes of the octoploid strawberry, F. × ananassa, and the diploid strawberry, F. vesca. The 144 samples were classified into seven pseudo-molecules of F. vesca. The coverage rates of the DNA sequences from the single chromosome onto all pseudo-molecular sequences varied from 3 to 29.9%.

CONCLUSION

We demonstrated an efficient method for sequence analysis of allopolyploid plants using microdissected single chromosomes. On the basis of our results, we believe that whole-genome analysis of allopolyploid plants can be enhanced using methodology that employs microdissected single chromosomes.

How Next-Generation Sequencing Has Aided Our Understanding of the Sequence Composition and Origin of B Chromosomes

Accessory, supernumerary, or-most simply-B chromosomes, are found in many eukaryotic karyotypes. These small chromosomes do not follow the usual pattern of segregation, but rather are transmitted in a higher than expected frequency. As increasingly being demonstrated by next-generation sequencing (NGS), their structure comprises fragments of standard (A) chromosomes, although in some plant species, their sequence also includes contributions from organellar genomes. Transcriptomic analyses of various animal and plant species have revealed that, contrary to what used to be the common belief, some of the B chromosome DNA is protein-encoding. This review summarizes the progress in understanding B chromosome biology enabled by the application of next-generation sequencing technology and state-of-the-art bioinformatics. In particular, a contrast is drawn between a direct sequencing approach and a strategy based on a comparative genomics as alternative routes that can be taken towards the identification of B chromosome sequences.

A time- and cost-effective strategy to sequence mammalian Y Chromosomes: an application to the de novo assembly of gorilla Y

The mammalian Y Chromosome sequence, critical for studying male fertility and dispersal, is enriched in repeats and palindromes, and thus, is the most difficult component of the genome to assemble. Previously, expensive and labor-intensive BAC-based techniques were used to sequence the Y for a handful of mammalian species. Here, we present a much faster and more affordable strategy for sequencing and assembling mammalian Y Chromosomes of sufficient quality for most comparative genomics analyses and for conservation genetics applications. The strategy combines flow sorting, short- and long-read genome and transcriptome sequencing, and droplet digital PCR with novel and existing computational methods. It can be used to reconstruct sex chromosomes in a heterogametic sex of any species. We applied our strategy to produce a draft of the gorilla Y sequence. The resulting assembly allowed us to refine gene content, evaluate copy number of ampliconic gene families, locate species-specific palindromes, examine the repetitive element content, and produce sequence alignments with human and chimpanzee Y Chromosomes. Our results inform the evolution of the hominine (human, chimpanzee, and gorilla) Y Chromosomes. Surprisingly, we found the gorilla Y Chromosome to be similar to the human Y Chromosome, but not to the chimpanzee Y Chromosome. Moreover, we have utilized the assembled gorilla Y Chromosome sequence to design genetic markers for studying the male-specific dispersal of this endangered species.

The Chromosome Microdissection and Microcloning Technique

Chromosome microdissection followed by microcloning is an efficient tool combining cytogenetics and molecular genetics that can be used for the construction of the high density molecular marker linkage map and fine physical map, the generation of probes for chromosome painting, and the localization and cloning of important genes. Here, we describe a modified technique to microdissect a single chromosome, paint individual chromosomes, and construct single-chromosome DNA libraries.

Laser microdissection-based analysis of the Y sex chromosome of the Antarctic fish Chionodracohamatus (Notothenioidei, Channichthyidae)

Microdissection, DOP-PCR amplification and microcloning were used to study the large Y chromosome of Chionodracohamatus, an Antarctic fish belonging to the Notothenioidei, the dominant component of the Southern Ocean fauna. The species has evolved a multiple sex chromosome system with digametic males showing an X1YX2 karyotype and females an X1X1X2X2 karyotype. Fluorescence in situ hybridization, performed with a painting probe made from microdissected Y chromosomes, allowed a deeper insight on the chromosomal rearrangement, which underpinned the fusion event that generated the Y. Then, we used a DNA library established by microdissection and microcloning of the whole Y chromosome of Chionodracohamatus for searching sex-linked sequences. One clone provided preliminary information on the presence on the Y chromosome of the CHD1 gene homologue, which is sex-linked in birds but in no other vertebrates. Several clones from the Y-chromosome mini-library contained microsatellites and transposable elements, one of which mapped to the q arm putative fusion region of the Y chromosome. The findings confirm that interspersed repetitive sequences might have fostered chromosome rearrangements and the emergence of the Y chromosome in Chionodracohamatus. Detection of the CHD1 gene in the Y sex-determining region could be a classical example of convergent evolution in action.

Intragenomic distribution of RTE retroelements suggests intrachromosomal movement

Much is known about the abundance of transposable elements (TEs) in eukaryotic genomes, but much is still unknown on their behaviour within cells. We employ here a combination of cytological, molecular and genomic approaches providing information on the intragenomic distribution and behaviour of non-long terminal repeat (LTR) retrotransposon-like elements (RTE). We microdissected every chromosome in a single first meiotic metaphase cell of the grasshopper Eyprepocnemis plorans and polymerase chain reaction (PCR) amplified a fragment of the RTE reverse transcriptase gene with specific primers. PCR products were cloned and 139 clones were sequenced. Analysis of molecular variance (AMOVA) showed significant intragenomic structure for these elements, with 4.6 % of molecular variance being found between chromosomes. A maximum likelihood tree built with the RTE sequences revealed the frequent presence of two or more elements showing very high similarity and being located on the same chromosome, thus suggesting intrachromosome movement. The 454 pyrosequencing of genomic DNA gave strong support to the microdissection results and provided evidence for the existence of 5' truncated elements. Our results thus indicate a tendency of RTE elements to reinsert into the same chromosome from where they were transcribed, which could be achieved if retrotranscription and insertion takes place immediately after transcription.

Flow sorting and sequencing meadow fescue chromosome 4F

The analysis of large genomes is hampered by a high proportion of repetitive DNA, which makes the assembly of short sequence reads difficult. This is also the case in meadow fescue (Festuca pratensis), which is known for good abiotic stress resistance and has been used in intergeneric hybridization with ryegrasses (Lolium spp.) to produce Festulolium cultivars. In this work, we describe a new approach to analyze the large genome of meadow fescue, which involves the reduction of sample complexity without compromising information content. This is achieved by dissecting the genome to smaller parts: individual chromosomes and groups of chromosomes. As the first step, we flow sorted chromosome 4F and sequenced it by Illumina with approximately 50× coverage. This provided, to our knowledge, the first insight into the composition of the fescue genome, enabled the construction of the virtual gene order of the chromosome, and facilitated detailed comparative analysis with the sequenced genomes of rice (Oryza sativa), Brachypodium distachyon, sorghum (Sorghum bicolor), and barley (Hordeum vulgare). Using GenomeZipper, we were able to confirm the collinearity of chromosome 4F with barley chromosome 4H and the long arm of chromosome 5H. Several new tandem repeats were identified and physically mapped using fluorescence in situ hybridization. They were found as robust cytogenetic markers for karyotyping of meadow fescue and ryegrass species and their hybrids. The ability to purify chromosome 4F opens the way for more efficient analysis of genomic loci on this chromosome underlying important traits, including freezing tolerance. Our results confirm that next-generation sequencing of flow-sorted chromosomes enables an overview of chromosome structure and evolution at a resolution never achieved before.

Common descent of B chromosomes in two species of the fish genus Prochilodus (Characiformes, Prochilodontidae)

To ascertain the origin of B chromosomes in 2 fish species of the genus Prochilodus, i.e. P. lineatus and P. nigricans, we microdissected them and generated B-specific DNA probes. These probes were used to perform chromosome painting in both species and in 3 further ones belonging to the same genus (P. argenteus, P. brevis and P. costatus). Both probes hybridized with the B chromosomes in P. lineatus and P. nigricans, but with none of the chromosomes in the 5 species. This indicates that the B chromosomes have low similarity with DNAs located in the A chromosomes and suggests the possibility that the B chromosomes in the 2 species have a common origin. The most parsimonious explanation would imply intergeneric hybridization in an ancestor of P. lineatus and P. nigricans yielding the B chromosome as a byproduct, which remained in these 2 species after their phylogenetic origin, but was perhaps lost in other Prochilodus species. This hypothesis predicts that B chromosomes are old genomic elements in this genus, and this could be tested once a species from a relative genus would be found showing homology of its A chromosomes with the B-probes employed here, through a comparison of B chromosome DNA sequences with those in the A chromosomes of this other species.

Microdissection and painting of the Y chromosome in spinach (Spinacia oleracea)

Spinach has long been used as a model for genetic and physiological studies of sex determination and expression. Although trisomic analysis from a cross between diploid and triploid plants identified the XY chromosome as the largest chromosome, no direct evidence has been provided to support this at the molecular level. In this study, the largest chromosomes of spinach from mitotic metaphase spreads were microdissected using glass needles. Degenerate oligonucleotide primed polymerase chain reaction was used to amplify the dissected chromosomes. The amplified products from the Y chromosome were identified using the male-specific marker T11A. For the first time, the largest spinach chromosome was confirmed to be a sex chromosome at the molecular level. PCR products from the isolated chromosomes were used in an in situ probe mixture for painting the Y chromosome. The fluorescence signals were mainly distributed on all chromosomes and four pair of weaker punctate fluorescence signal sites were observed on the terminal region of two pair of autosomes. These findings provide a foundation for the study of sex chromosome evolution in spinach.

Combined fluorescent-chromogenic in situ hybridization for identification and laser microdissection of interphase chromosomes

Chromosome territories constitute the most conspicuous feature of nuclear architecture, and they exhibit non-random distribution patterns in the interphase nucleus. We observed that in cell nuclei from humans with Down Syndrome two chromosomes 21 frequently localize proximal to one another and distant from the third chromosome. To systematically investigate whether the proximally positioned chromosomes were always the same in all cells, we developed an approach consisting of sequential FISH and CISH combined with laser-microdissection of chromosomes from the interphase nucleus and followed by subsequent chromosome identification by microsatellite allele genotyping. This approach identified proximally positioned chromosomes from cultured cells, and the analysis showed that the identity of the chromosomes proximally positioned varies. However, the data suggest that there may be a tendency of the same chromosomes to be positioned close to each other in the interphase nucleus of trisomic cells. The protocol described here represents a powerful new method for genome analysis.

Origin and molecular organization of supernumerary chromosomes of Prochilodus lineatus (characiformes, prochilodontidae) obtained by DNA probes

In Prochilodus lineatus B-chromosomes are visualized as reduced size extra elements identified as microchromosomes and are variable in morphology and number. We describe the specific total probe (B-chromosome probe) in P. lineatus obtained by chromosome microdissection and a whole genomic probe (genomic probe) from an individual without B-chromosome. The specific B-chromosome was scraped and processed to obtain DNA with amplification by DOP-PCR, and so did the genomic probe DNA. Fluorescence in situ hybridization using the B-chromosome probe labeled with dUTP-Tetramethyl-rhodamine and the genomic probe labeled with digoxigenin-FITC permitted to establish that in this species supernumerary chromosomes with varying number and morphology had different structure of chromatin when compared to that of the regular chromosomes or A complement, since only these extra elements were labeled in the metaphases. The present findings suggest that modifications in the chromatin structure of B-chromosomes to differentiate them from the A chromosomes could occur along their dispersion in the individuals of the population.

Isolation of expressed sequences from a specific chromosome of Thinopyrum intermedium infected by BYDV

To map important ESTs to specific chromosomes and (or) chromosomal regions is difficult in hexaploid wheat because of its large genome size and serious interference of homoeologous sequences. Large-scale EST sequencing and subsequent chromosome localization are both laborious and time-consuming. The wheat alien addition line TAi-27 contains a pair of chromosomes of Thinopyrum intermedium (Host) Barkworth & D.R. Dewey that carry the resistance gene against barley yellow dwarf virus. In this research, we developed a modified technique based on chromosome microdissection and hybridization-specific amplification to isolate expressed sequences from the alien chromosome of TAi-27 by hybridization between the DNA of the microdissected alien chromosome and cDNA of Th. intermedium infected by barley yellow dwarf virus. Twelve clones were selected, sequenced, and analyzed. Three of them were unknown genes without any hit in the GenBank database and the other nine were highly homologous with ESTs of wheat, barley, and (or) other plants in Gramineae induced by abiotic or biotic stress. The method used in this research to isolate expressed sequences from a specific chromosome has the following advantages: (i) the obtained expressed sequences are larger in size and have 3' end information and (ii) the operation is less complicated. It would be an efficient improved method for genomics and functional genomics research of polyploid plants, especially for EST development and mapping. The obtained expressed sequence data are also informative in understanding the resistance genes on the alien chromosome of TAi-27.