Next-generation sequencing, FISH mapping and synteny-based modeling reveal mechanisms of decreasing dysploidy in Cucumis

Translocations and inversions: major chromosomal rearrangements during Vigna (Leguminosae) evolution

KEY MESSAGE

Inversions and translocations are the major chromosomal rearrangements involved in Vigna subgenera evolution, being Vigna vexillata the most divergent species. Centromeric repositioning seems to be frequent within the genus. Oligonucleotide-based fluorescence in situ hybridization (Oligo-FISH) provides a powerful chromosome identification system for inferring plant chromosomal evolution. Aiming to understand macrosynteny, chromosomal diversity, and the evolution of bean species from five Vigna subgenera, we constructed cytogenetic maps for eight taxa using oligo-FISH-based chromosome identification. We used oligopainting probes from chromosomes 2 and 3 of Phaseolus vulgaris L. and two barcode probes designed from V. unguiculata (L.) Walp. genome. Additionally, we analyzed genomic blocks among the Ancestral Phaseoleae Karyotype (APK), two V. unguiculata subspecies (V. subg. Vigna), and V. angularis (Willd.) Ohwi & Ohashi (V. subg. Ceratotropis). We observed macrosynteny for chromosomes 2, 3, 4, 6, 7, 8, 9, and 10 in all investigated taxa except for V. vexillata (L.) A. Rich (V. subg. Plectrotropis), in which only chromosomes 4, 7, and 9 were unambiguously identified. Collinearity breaks involved with chromosomes 2 and 3 were revealed. We identified minor differences in the painting pattern among the subgenera, in addition to multiple intra- and interblock inversions and intrachromosomal translocations. Other rearrangements included a pericentric inversion in chromosome 4 (V. subg. Vigna), a reciprocal translocation between chromosomes 1 and 5 (V. subg. Ceratotropis), a potential deletion in chromosome 11 of V. radiata (L.) Wilczek, as well as multiple intrablock inversions and centromere repositioning via genomic blocks. Our study allowed the visualization of karyotypic patterns in each subgenus, revealing important information for understanding intrageneric karyotypic evolution, and suggesting V. vexillata as the most karyotypically divergent species.

Genome size, chromosome number variation and its correlation with stomatal characters for assessment of ploidy levels in a core subset of mulberry (Morus spp.) germplasm

The large size of the germplasm collection along with scanty information on their cytological and genome constitution have hindered well-planned breeding schemes in mulberry. To address the issue, a study was undertaken to investigate the variability in DNA content and genome size, chromosome number, ploidy and its relation with important stomatal characteristics among 162 mulberry germplasm. These germplasm comprise a core subset of 150 collections along with a representative collection of different mulberry species including the wild. Among the germplasm belonging to 16 species, we identified 122 diploids (2n = 28), 4 aneuploids (2n = 30), 13 triploids (2n = 42), 15 tetraploids (2n = 56), 7 hexaploids (2n = 84) and 1 dodecosaploid (2n = 308) based on the chromosome count. Most of the cultivated mulberries are found to be diploids. The mean nuclear 2C DNA content estimated by Flow cytometry, varied from 0.723±0.006 pg (M. australis, 2n = 2x) to 7.732 pg (M. nigra, 2n = 22x). The 2C DNA content positively correlated with the ploidy status and stomatal length (r = 0.814, p<0.001). Based on the 1Cx value, the study also suggests that the majority of the polyploid species have experienced genome downsizing in relation to their diploid progenitors. This study provides the most essential information on chromosome number, ploidy and DNA content to facilitate the utilization of a core subset of germplasm in the mulberry breeding program.

Insight into the Organization of the B10v3 Cucumber Genome by Integration of Biological and Bioinformatic Data

The availability of a well-organized and annotated reference genome is essential for genome research and the analysis of re-sequencing approaches. The B10v3 cucumber (Cucumis sativus L.) reference genome has been sequenced and assembled into 8035 contigs, a small fraction of which have been mapped to individual chromosomes. Currently, bioinformatics methods based on comparative homology have made it possible to re-order the sequenced contigs by mapping them to the reference genomes. The B10v3 genome (North-European, Borszczagowski line) was rearranged against the genomes of cucumber 9930 ('Chinese Long' line) and Gy14 (North American line). Furthermore, a better insight into the organization of the B10v3 genome was obtained by integrating the data available in the literature on the assignment of contigs to chromosomes in the B10v3 genome with the results of the bioinformatic analysis. The combination of information on the markers used in the assembly of the B10v3 genome and the results of FISH and DArT-seq experiments confirmed the reliability of the in silico assignment. Approximately 98% of the protein-coding genes within the chromosomes were assigned and a significant proportion of the repetitive fragments in the sequenced B10v3 genome were identified using the RagTag programme. In addition, BLAST analyses provided comparative information between the B10v3 genome and the 9930 and Gy14 data sets. This revealed both similarities and differences in the functional proteins found between the coding sequences region in the genomes. This study contributes to better knowledge and understanding of cucumber genome line B10v3.

Comparative cytogenomics reveals genome reshuffling and centromere repositioning in the legume tribe Phaseoleae

The tribe Phaseoleae includes several legume crops with assembled genomes. Comparative genomic studies have evidenced the preservation of large genomic blocks among legumes, although chromosome dynamics during Phaseoleae evolution has not been investigated. We conducted a comparative genomic analysis to define an informative genomic block (GB) system and to reconstruct the ancestral Phaseoleae karyotype (APK). We identified GBs based on the orthologous genes between Phaseolus vulgaris and Vigna unguiculata and searched for GBs in different genomes of the Phaseolinae (P. lunatus) and Glycininae (Amphicarpaea edgeworthii) subtribes and Spatholobus suberectus (sister to Phaseolinae and Glycininae), using Medicago truncatula as the outgroup. We also used oligo-FISH probes of two P. vulgaris chromosomes to paint the orthologous chromosomes of two non-sequenced Phaseolinae species. We inferred the APK as having n = 11 and 19 GBs (A to S), hypothesizing five chromosome fusions that reduced the ancestral legume karyotype to n = 11. We identified the rearrangements among the APK and the subtribes and species, with extensive centromere repositioning in Phaseolus. We also reconstructed the chromosome number reduction in S. suberectus. The development of the GB system and the proposed APK provide useful approaches for future comparative genomic analyses of legume species.

The Divergence of Chromosome Structures and 45S Ribosomal DNA Organization in Cucumis debilis Inferred by Comparative Molecular Cytogenetic Mapping

Cucumis debilis W.J.de Wilde & Duyfjes is an annual and monoecious plant. This species is endemic to Southeast Asia, particularly Vietnam. However, C. debilis is rarely studied, and no detailed information is available regarding its basic chromosome number, 45S ribosomal DNA (rDNA) status, and divergence among other Cucumis species. In this study, we characterized the morphological characters and determined and investigated the basic chromosome number and chromosomal distribution of 45S rDNA of C. debilis using the fluorescent in situ hybridization (FISH) technique. A maximum likelihood tree was constructed by combining the chloroplast and internal transcribed spacer of 45S rDNAs to infer its relationship within Cucumis. C. debilis had an oval fruit shape, green fruit peel, and protrusion-like white spots during the immature fruit stage. FISH analysis using 45S rDNA probe showed three pairs of 45S rDNA loci located at the terminal region in C. debilis, similar to C. hystrix. Meanwhile, two, two, and five pairs of 45S rDNA loci were observed for C. melo, C. metuliferus, and C. sativus, respectively. One melon (P90) and cucumber accessions exhibited different chromosomal localizations compared with other members of Cucumis. The majority of Cucumis species showed the terminal location of 45S rDNA, but melon P90 and cucumber exhibited terminal–interstitial and all interstitial orientations of 45S rDNA loci. Based on molecular cytogenetics and phylogenetic evidence, C. debilis is more closely related to cucumber than melon. Therefore, C. debilis may serve as a potential parental accession for genetic improvement of cucumber through interspecific hybridization.

Celebrating Mendel, McClintock, and Darlington: On end-to-end chromosome fusions and nested chromosome fusions

The evolution of eukaryotic genomes is accompanied by fluctuations in chromosome number, reflecting cycles of chromosome number increase (polyploidy and centric fissions) and decrease (chromosome fusions). Although all chromosome fusions result from DNA recombination between two or more nonhomologous chromosomes, several mechanisms of descending dysploidy are exploited by eukaryotes to reduce their chromosome number. Genome sequencing and comparative genomics have accelerated the identification of inter-genome chromosome collinearity and gross chromosomal rearrangements and have shown that end-to-end chromosome fusions (EEFs) and nested chromosome fusions (NCFs) may have played a more important role in the evolution of eukaryotic karyotypes than previously thought. The present review aims to summarize the limited knowledge on the origin, frequency, and evolutionary implications of EEF and NCF events in eukaryotes and especially in land plants. The interactions between nonhomologous chromosomes in interphase nuclei and chromosome (mis)pairing during meiosis are examined for their potential importance in the origin of EEFs and NCFs. The remaining open questions that need to be addressed are discussed.

Identification and Characterization Roles of Phytoene Synthase (PSY) Genes in Watermelon Development

Phytoene synthase (PSY) plays an essential role in carotenoid biosynthesis. In this study, three ClPSY genes were identified through the watermelon genome, and their full-length cDNA sequences were cloned. The deduced proteins of the three ClPSY genes were ranged from 355 to 421 amino acid residues. Phylogenetic analysis suggested that the ClPSYs are highly conserved with bottle gourd compared to other cucurbit crops PSY proteins. Variation in ClPSY1 expression in watermelon with different flesh colors was observed; ClPSY1 was most highly expressed in fruit flesh and associated with the flesh color formation. ClPSY1 expression was much lower in the white-fleshed variety than the colored fruits. Gene expression analysis of ClPSY genes in root, stem, leaf, flower, ovary and flesh of watermelon plants showed that the levels of ClPSY2 transcripts found in leaves was higher than other tissues; ClPSY3 was dominantly expressed in roots. Functional complementation assays of the three ClPSY genes suggested that all of them could encode functional enzymes to synthesize the phytoene from Geranylgeranyl Pyrophosphate (GGPP). Some of the homologous genes clustered together in the phylogenetic tree and located in the synteny chromosome region seemed to have similar expression profiles among different cucurbit crops. The findings provide a foundation for watermelon flesh color breeding with regard to carotenoid synthesis and also provide an insight for the further research of watermelon flesh color formation.

Functional copy number variation of CsSHINE1 is associated with fruit skin netting intensity in cucumber, Cucumis sativus

Fruit skin netting in cucumber (Cucumis sativus) is associated with important fruit quality attributes. Two simply inherited genes H (Heavy netting) and Rs (Russet skin) control skin netting, but their molecular basis is unknown. Here, we reported map-based cloning and functional characterization of the candidate gene for the Rs locus that encodes CsSHINE1 (CsSHN1), an AP2 domain containing ethylene-responsive transcription factor protein. Comparative phenotypic analysis in near-isogenic lines revealed that fruit with netted skin had different epidermal structures from that with smooth skin including thicker cuticles, smaller, palisade-shaped epidermal and sub-epidermal cells with heavily suberized and lignified cell walls, higher peroxidase activities, which suggests multiple functions of CsSHN1 in regulating fruit skin netting and epidermal cell patterning. Among three representative cucumber inbred lines, three haplotypes at three polymorphic sites were identified inside CsSHN1: a functional copy in Gy14 (wild type) with light fruit skin netting, a copy number variant with two tandemly arrayed functional copies in WI7120 with heavy skin netting, and a loss-of-function copy in 9930 with smooth skin. The expression level of CsSHN1 in fruit exocarp of three lines was positively correlated with the skin netting intensity. Comparative analysis between cucumber and melon revealed conserved and divergent genetic mechanisms underlying fruit skin netting/reticulation that may reflect the different selection histories in the two crops. A discussion was made on genetic basis of fruit skin netting in the context of natural and artificial selections of fruit quality-related epidermal features during cucumber breeding.

Chromosome restructuring and number change during the evolution of Morus notabilis and Morus alba

Mulberry (Morus spp.) is an economically important plant as the main food plant used for rearing domesticated silkworm and it has multiple uses in traditional Chinese medicine. Two basic chromosome numbers (Morus notabilis, n = 7, and Morus alba, n = 14) have been reported in the genus Morus, but the evolutionary history and relationship between them remain unclear. In the present study, a 335-Mb high-quality chromosome-scale genome was assembled for the wild mulberry species M. notabilis. Comparative genomic analyses indicated high chromosomal synteny between the 14 chromosomes of cultivated M. alba and the six chromosomes of wild M. notabilis. These results were successfully verified by fluorescence in situ hybridization. Chromosomal fission/fusion events played crucial roles in the chromosome restructuring process between M. notabilis and M. alba. The activity of the centromere was another key factor that ensured the stable inheritance of chromosomes. Our results also revealed that long terminal repeat retrotransposons were a major driver of the genome divergence and evolution of the mulberry genomes after they diverged from each other. This study provides important insights and a solid foundation for studying the evolution of mulberry, allowing the accelerated genetic improvement of cultivated mulberry species.

Genome-wide identification of small heat shock protein (HSP20) homologs in three cucurbit species and the expression profiles of CsHSP20s under several abiotic stresses

Small heat shock protein (HSP20) genes play important roles in biological processes of plants. In this study, a total of 47 CsHSP20 genes, 45 CmHSP20 genes, and 47 ClHSP20 genes were genome-wide identified by 'hmmsearch' and BLASTP using the latest versions of cucumber, melon, and watermelon genomes, respectively. According to the phylogenetic relationships and predicted subcellular localizations, HSP20s of these three cucurbit species were divided into 8 subfamilies (CI-CIV, CP, ER, M, and PX), in which some HSP20s were closely related with each other based on the collinearity analysis. Specific expression patterns of CsHSP20s were checked in 10 different tissues of cucumber plants. RNA-seq analysis of transcript levels, combined with cis-acting elements and GO enrichment analysis suggested that CsHSP20s were responsive to several different types of abiotic stresses, including chilling, temperature and photoperiod, high temperature and high humidity, and salinity. In conclusion, results of this work not only provided valuable information for exploring the regulating mechanisms of CsHSP20s in responding to abiotic stresses in cucumber, but also shed light on the potentially evolutional relations among cucumber, melon, and watermelon from a perspective of comparative genomics that specified on HSP20 gene families.

Oligo-FISH barcode in beans: a new chromosome identification system

An Oligo-FISH barcode system was developed for two model legumes, allowing the identification of all cowpea and common bean chromosomes in a single FISH experiment, and revealing new chromosome rearrangements. The FISH barcode system emerges as an effective tool to understand the chromosome evolution of economically important legumes and their related species. Current status on plant cytogenetic and cytogenomic research has allowed the selection and design of oligo-specific probes to individually identify each chromosome of the karyotype in a target species. Here, we developed the first chromosome identification system for legumes based on oligo-FISH barcode probes. We selected conserved genomic regions between Vigna unguiculata (Vu, cowpea) and Phaseolus vulgaris (Pv, common bean) (diverged ~ 9.7-15 Mya), using cowpea as a reference, to produce a unique barcode pattern for each species. We combined our oligo-FISH barcode pattern with a set of previously developed FISH probes based on BACs and ribosomal DNA sequences. In addition, we integrated our FISH maps with genome sequence data. Based on this integrated analysis, we confirmed two translocation events (involving chromosomes 1, 5, and 8; and chromosomes 2 and 3) between both species. The application of the oligo-based probes allowed us to demonstrate the participation of chromosome 5 in the translocation complex for the first time. Additionally, we detailed a pericentric inversion on chromosome 4 and identified a new paracentric inversion on chromosome 10. We also detected centromere repositioning associated with chromosomes 2, 3, 5, 7, and 9, confirming previous results for chromosomes 2 and 3. This first barcode system for legumes can be applied for karyotyping other Phaseolinae species, especially non-model, orphan crop species lacking genomic assemblies and cytogenetic maps, expanding our understanding of the chromosome evolution and genome organization of this economically important legume group.

Genetic Resources and Vulnerabilities of Major Cucurbit Crops

The Cucurbitaceae family provides numerous important crops including watermelons (Citrullus lanatus), melons (Cucumis melo), cucumbers (Cucumis sativus), and pumpkins and squashes (Cucurbita spp.). Centers of domestication in Africa, Asia, and the Americas were followed by distribution throughout the world and the evolution of secondary centers of diversity. Each of these crops is challenged by multiple fungal, oomycete, bacterial, and viral diseases and insects that vector disease and cause feeding damage. Cultivated varieties are constrained by market demands, the necessity for climatic adaptations, domestication bottlenecks, and in most cases, limited capacity for interspecific hybridization, creating narrow genetic bases for crop improvement. This analysis of crop vulnerabilities examines the four major cucurbit crops, their uses, challenges, and genetic resources. ex situ germplasm banks, the primary strategy to preserve genetic diversity, have been extensively utilized by cucurbit breeders, especially for resistances to biotic and abiotic stresses. Recent genomic efforts have documented genetic diversity, population structure, and genetic relationships among accessions within collections. Collection size and accessibility are impacted by historical collections, current ability to collect, and ability to store and maintain collections. The biology of cucurbits, with insect-pollinated, outcrossing plants, and large, spreading vines, pose additional challenges for regeneration and maintenance. Our ability to address ongoing and future cucurbit crop vulnerabilities will require a combination of investment, agricultural, and conservation policies, and technological advances to facilitate collection, preservation, and access to critical Cucurbitaceae diversity.

Reconstruction of ancestral karyotype illuminates chromosome evolution in the genus Cucumis

Karyotype dynamics driven by complex chromosome rearrangements constitute a fundamental issue in evolutionary genetics. The evolutionary events underlying karyotype diversity within plant genera, however, have rarely been reconstructed from a computed ancestral progenitor. Here, we developed a method to rapidly and accurately represent extant karyotypes with the genus, Cucumis, using highly customizable comparative oligo-painting (COP) allowing visualization of fine-scale genome structures of eight Cucumis species from both African-origin and Asian-origin clades. Based on COP data, an evolutionary framework containing a genus-level ancestral karyotype was reconstructed, allowing elucidation of the evolutionary events that account for the origin of these diverse genomes within Cucumis. Our results characterize the cryptic rearrangement hotspots on ancestral chromosomes, and demonstrate that the ancestral Cucumis karyotype (n = 12) evolved to extant Cucumis genomes by hybridizations and frequent lineage- and species-specific genome reshuffling. Relative to the African species, the Asian species, including melon (Cucumis melo, n = 12), Cucumis hystrix (n = 12) and cucumber (Cucumis sativus, n = 7), had highly shuffled genomes caused by large-scale inversions, centromere repositioning and chromothripsis-like rearrangement. The deduced reconstructed ancestral karyotype for the genus allowed us to propose evolutionary trajectories and specific events underlying the origin of these Cucumis species. Our findings highlight that the partitioned evolutionary plasticity of Cucumis karyotype is primarily located in the centromere-proximal regions marked by rearrangement hotspots, which can potentially serve as a reservoir for chromosome evolution due to their fragility.

High-quality chromosome-level genomes of Cucumis metuliferus and Cucumis melo provide insight into Cucumis genome evolution

Cucumis metuliferus (African horned cucumber), a wild relative of Cucumis sativus (cucumber) and Cucumis melo (melon), displays high-level resistance to several important plant pathogens (e.g., root-knot nematodes and several viruses). Here, we report a chromosome-level genome assembly for C. metuliferus, with a 316 Mb genome sequence comprising 29 039 genes. Phylogenetic analysis of related species in family Cucurbitaceae indicated that the divergence time between C. metuliferus and melon was 17.8 million years ago. Comparisons between the C. metuliferus and melon genomes revealed large structural variations (inversions and translocations >1 Mb) in eight chromosomes of these two species. Gene family comparison showed that C. metuliferus has the largest number of resistance-related nucleotide-binding site leucine-rich repeat (NBS-LRR) genes in Cucurbitaceae. The loss of NBS-LRR loci caused by large insertions or deletions (indels) and pseudogenization caused by small indels explained the loss of NBS-LRR genes in Cucurbitaceae. Population structure analysis suggested that C. metuliferus originated in Zimbabwe, then spread to other southern African regions where it likely underwent similar domestic selection as melon. This C. metuliferus reference sequence will accelerate the understanding of the molecular evolution of resistance-related genes and enhance cucurbit crop improvement efforts.

Genome Wide Characterization, Comparative and Genetic Diversity Analysis of Simple Sequence Repeats in Cucurbita Species

Simple sequence repeats (SSRs) are widely used in mapping constructions and comparative and genetic diversity analyses. Here, 103,056 SSR loci were found in Cucurbita species by in silico PCR. In general, the frequency of these SSRs decreased with the increase in the motif length, and di-nucleotide motifs were the most common type. For the same repeat types, the SSR frequency decreased sharply with the increase in the repeat number. The majority of the SSR loci were suitable for marker development (84.75% in Cucurbita moschata, 94.53% in Cucurbita maxima, and 95.09% in Cucurbita pepo). Using these markers, the cross-species transferable SSR markers between C. pepo and other Cucurbitaceae species were developed, and the complicated mosaic relationships among them were analyzed. Especially, the main syntenic relationships between C. pepo and C. moschata or C. maxima indicated that the chromosomes in the Cucurbita genomes were highly conserved during evolution. Furthermore, 66 core SSR markers were selected to measure the genetic diversity in 61 C. pepo germplasms, and they were divided into two groups by structure and unweighted pair group method with arithmetic analysis. These results will promote the utilization of SSRs in basic and applied research of Cucurbita species.

A Dual-Color Tyr-FISH Method for Visualizing Genes/Markers on Plant Chromosomes to Create Integrated Genetic and Cytogenetic Maps

In situ imaging of molecular markers on a physical chromosome is an indispensable tool for refining genetic maps and validation genome assembly at the chromosomal level. Despite the tremendous progress in genome sequencing, the plant genome assembly at the chromosome level remains a challenge. Recently developed optical and Hi-C mapping are aimed at assistance in genome assembly. For high confidence in the genome assembly at chromosome level, more independent approaches are required. The present study is aimed at refining an ultrasensitive Tyr-FISH technique and developing a reliable and simple method of in situ mapping of a short unique DNA sequences on plant chromosomes. We have carefully analyzed the critical steps of the Tyr-FISH to find out the reasons behind the flaws of this technique. The accurate visualization of markers/genes appeared to be significantly dependent on the means of chromosome slide preparation, probe design and labeling, and high stringency washing. Appropriate adjustment of these steps allowed us to detect a short DNA sequence of 1.6 Kb with a frequency of 51.6%. Based on our results, we developed a more reliable and simple protocol for dual-color Tyr-FISH visualization of unique short DNA sequences on plant chromosomes. This new protocol can allow for more accurate determination of the physical distance between markers and can be applied for faster integration of genetic and cytogenetic maps.

Chromosome-scale genome assembly of Cucumis hystrix-a wild species interspecifically cross-compatible with cultivated cucumber

Cucumis hystrix Chakr. (2n = 2x = 24) is a wild species that can hybridize with cultivated cucumber (C. sativus L., 2n = 2x = 14), a globally important vegetable crop. However, cucumber breeding is hindered by its narrow genetic base. Therefore, introgression from C. hystrix has been anticipated to bring a breakthrough in cucumber improvement. Here, we report the chromosome-scale assembly of C. hystrix genome (289 Mb). Scaffold N50 reached 14.1 Mb. Over 90% of the sequences were anchored onto 12 chromosomes. A total of 23,864 genes were annotated using a hybrid method. Further, we conducted a comprehensive comparative genomic analysis of cucumber, C. hystrix, and melon (C. melo L., 2n = 2x = 24). Whole-genome comparisons revealed that C. hystrix is phylogenetically closer to cucumber than to melon, providing a molecular basis for the success of its hybridization with cucumber. Moreover, expanded gene families of C. hystrix were significantly enriched in "defense response," and C. hystrix harbored 104 nucleotide-binding site-encoding disease resistance gene analogs. Furthermore, 121 genes were positively selected, and 12 (9.9%) of these were involved in responses to biotic stimuli, which might explain the high disease resistance of C. hystrix. The alignment of whole C. hystrix genome with cucumber genome and self-alignment revealed 45,417 chromosome-specific sequences evenly distributed on C. hystrix chromosomes. Finally, we developed four cucumber-C. hystrix alien addition lines and identified the exact introgressed chromosome using molecular and cytological methods. The assembled C. hystrix genome can serve as a valuable resource for studies on Cucumis evolution and interspecific introgression breeding of cucumber.

Agarose-resolvable InDel markers based on whole genome re-sequencing in cucumber

Insertion and Deletion (InDel) are common features in genomes and are associated with genetic variation. The whole-genome re-sequencing data from two parents (X1 and X2) of the elite cucumber (Cucumis sativus) hybrid variety Lvmei No.1 was used for genome-wide InDel polymorphisms analysis. Obtained sequence reads were mapped to the genome reference sequence of Chinese fresh market type inbred line ‘9930’ and gaps conforming to InDel were pinpointed. Further, the level of cross-parents polymorphism among five pairs of cucumber breeding parents and their corresponding hybrid varieties were used for evaluating hybrid seeds purity test efficiency of InDel markers. A panel of 48 cucumber breeding lines was utilized for PCR amplification versatility and phylogenetic analysis of these markers. In total, 10,470 candidate InDel markers were identified for X1 and X2. Among these, 385 markers with more than 30 nucleotide difference were arbitrary chosen. These markers were selected for experimental resolvability through electrophoresis on an Agarose gel. Two hundred and eleven (211) accounting for 54.81% of markers could be validated as single and clear polymorphic pattern while 174 (45.19%) showed unclear or monomorphic genetic bands between X1 and X2. Cross-parents polymorphism evaluation recorded 68 (32.23%) of these markers, which were designated as cross-parents transferable (CPT) InDel markers. Interestingly, the marker InDel114 presented experimental transferability between cucumber and melon. A panel of 48 cucumber breeding lines including parents of Lvmei No. 1 subjected to PCR amplification versatility using CPT InDel markers successfully clustered them into fruit and common cucumber varieties based on phylogenetic analysis. It is worth noting that 16 of these markers were predominately associated to enzymatic activities in cucumber. These agarose-based InDel markers could constitute a valuable resource for hybrid seeds purity testing, germplasm classification and marker-assisted breeding in cucumber.

An efficient Oligo-FISH painting system for revealing chromosome rearrangements and polyploidization in Triticeae

A chromosome-specific painting technique has been developed which combines the most recent approaches of the companion disciplines of molecular cytogenetics and genome research. We developed seven oligonucleotide (oligo) pools derivd from single-copy sequences on chromosomes 1 to 7 of barley (Hordeum vulgare L.) and corresponding collinear regions of wheat (Triticum aestivum L.). The seven groups of pooled oligos comprised between 10 986 and 12 496 45-bp monomers, and these then produced stable fluorescence in situ hybridization (FISH) signals on chromosomes of each linkage group of wheat and barley. The pooled oligo probes were applied to high-throughput karyotyping of the chromosomes of other Triticeae species in the genera Secale, Aegilops, Thinopyrum, and Dasypyrum, and the study also extended to some wheat-alien amphiploids and derived lines. We demonstrated that a complete set of whole-chromosome oligo painting probes facilitated the study of inter-species chromosome homologous relationships and visualized non-homologous chromosomal rearrangements in Triticeae species and some wheat-alien species derivatives. When combined with other non-denaturing FISH procedures using tandem-repeat oligos, the newly developed oligo painting techniques provide an efficient tool for the study of chromosome structure, organization, and evolution among any wild Triticeae species with non-sequenced genomes.

Quantitative trait loci for horticulturally important traits defining the Sikkim cucumber, Cucumis sativus var. sikkimensis

QTL mapping identified simply inherited genes and quantitative trait loci underlying morphologically characteristic traits of the Sikkim cucumber, which reveals their genetic basis during crop evolution. The data suggest the Sikkim cucumber as an ecotype of cultivated cucumber not worthy of formal taxonomic recognition. The Sikkim cucumber, Cucumis sativus var. sikkimensis, is featured with some morphological traits like black spine, brown fruit with fine and heavy netting, as well as large hollow in mature fruit. Despite its establishment as a botanical variety ~ 150 years ago, and its wide use as an important source of disease resistances in cucumber breeding, little is known about its taxonomic status and genetic basis of those characteristic traits. Here we reported QTL mapping with segregating populations derived from two Sikkim-type inbred lines, WI7088D and WI7120, and identification of 48 QTL underlying phenotypic variation for 18 horticulturally important traits. We found that the fruit spine and skin colors in the two populations were controlled by the previously cloned pleiotropic B (black spine) locus. The fruit netting in WI7088D and WI7120 was controlled by the well-known H (Heavy netting) and a novel Rs (Russet skin) locus, which was delimited to a 271-kb region on Chr5 and ~ 736-kb region on Chr1, respectively. A single major-effect QTL was detected for flowering time in each population (ft1.1 for WI7088D and ft6.2 for WI7120). Fifteen, six and five QTL were identified for fruit size, hollow size and flesh thickness variation in the two populations, respectively. No major structural changes were found between the Sikkim and cultivated cucumbers. Except for the rare allele at the Rs locus, there seem no private QTL/alleles identified from this study supporting the Sikkim cucumber as an ecotype of C. sativus, not worthy of formal taxonomic recognition.

Resequencing of 297 melon accessions reveals the genomic history of improvement and loci related to fruit traits in melon

Domestication and improvement are two important stages in crop evolution. Melon (Cucumis melo L.) is an important vegetable crop with wide phenotypic diversity in many horticultural traits, especially fruit size, flesh thickness and aroma, which are likely the results of long-term extensive selection during its evolution. However, selective signals in domestication and improvement stages for these remarkable variations remain unclear. We resequenced 297 wild, landrace and improved melon accessions and obtained 2 045 412 high-quality SNPs. Population structure and genetic diversity analyses revealed independent and two-step selections in two subspecies of melon: ssp. melo and ssp. agrestis during melon breeding. We detected 233 (~18.35 Mbp) and 159 (~17.71 Mbp) novel potential selective signals during the improvement stage in ssp. agrestis and spp. melo, respectively. Two alcohol acyltransferase genes (CmAATs) unique to the melon genome compared with other cucurbit crops may have undergone stronger selection in ssp. agrestis for the characteristic aroma as compared with other cucurbits. Genome-wide association analysis identified eight fruit size and seven flesh thickness signals overlapping with selective sweeps. Compared with thin-skinned ssp. agrestis, thick-skinned ssp. melo has undergone a stronger selection for thicker flesh. In most melon accessions, CmCLV3 has pleiotropic effects on carpel number and fruit shape. Findings from this study provide novel insights into melon crop evolution, and new tools to advance melon breeding.

Development of alien addition lines from Cucumis hystrix in Cucumis sativus: cytological and molecular marker analyses

Transferring desired genes from wild species to cultivars through alien addition lines (AALs) has been shown to be an effective method for genetic improvement. Cucumis hystrix Chakr. (HH, 2n = 24) is a wild species of Cucumis that possesses many resistant genes. A synthetic allotetraploid species, C. hytivus (HHCC, 2n = 38), was obtained from the cross between cultivated cucumber, C. sativus (CC, 2n = 14), and C. hystrix followed by chromosome doubling. Cucumis sativus - C. hystrix AALs were developed by continuous backcrossing to the cultivated cucumbers. In this study, 10 different types of AALs (CC-H01, CC-H06, CC-H08, CC-H10, CC-H12, CC-H06+H09, CC-H06+H10, CC-H06+H12, CC-H08+H10, CC-H01+H06+H10) were identified based on the analysis of fluorescence in situ hybridization (FISH) and molecular markers specific to C. hystrix chromosomes. And the behavior of the alien chromosomes in three AALs (CC-H01, CC-H06+H10, CC-H01+H06+H10) at meiosis was investigated. The results showed that alien chromosomes paired with C. sativus chromosome in few pollen mother cells (PMCs). Further, disomic alien addition lines (DAALs) carrying a pair of C. hystrix chromosome H10 were screened from the selfed progenies of CC-H10. Chromosome pairing between genomes provides cytological evidence for the possible introgression of alien chromosome segments. The development of AALs could serve as a key step for exploiting and utilizing valuable genes from C. hystrix.

QTL for horticulturally important traits associated with pleiotropic andromonoecy and carpel number loci, and a paracentric inversion in cucumber

The legendary cucumber inbred line WI2757 possesses a rare combination of resistances against nine pathogens, which is an important germplasm for cucumber breeding. However, WI2757 flowers late and does not perform well under field conditions. The genetic basis for horticulturally important traits other than disease resistances in WI2757 is largely unknown. In this study, we conducted QTL mapping using F₂ and recombinant inbred line (RIL) populations from the WI2757 × True Lemon cross that were segregating for multiple traits. Phenotypic data were collected in replicated field trials across multiple years for seven traits including fruit carpel number (CN) and sex expression. A high-density SNP-based genetic map was developed with genotyping by sequencing of the RIL population, which revealed a region on chromosome 1 with strong recombination suppression. The reduced recombination in this region was due to a ~ 10-Mbp paracentric inversion in WI2757 that was confirmed with additional segregation and cytological (FISH) analyses. Thirty-six QTL were detected for flowering time, fruit length (FL), fruit diameter (FD), fruit shape (LD), fruit number (FN), CN, and powdery mildew resistance. Five moderate- or major-effect QTL for FL, FD, LD, and FN inside the inversion are likely the pleiotropic effects of the andromonoecy (m), or the cn locus. The major-effect flowering time QTL ft1.1 was also mapped inside the inversion, which seems to be different from the previously assigned delayed flowering in WI2757. Implications of these findings on the use of WI2757 in cucumber breeding are discussed.

Flexible chromosome painting based on multiplex PCR of oligonucleotides and its application for comparative chromosome analyses in Cucumis

Chromosome painting is a powerful technique for chromosome and genome studies. We developed a flexible chromosome painting technique based on multiplex PCR of a synthetic oligonucleotide (oligo) library in cucumber (Cucumis sativus L., 2n = 14). Each oligo in the library was associated with a universal as well as nested specific primers for amplification, which allow the generation of different probes from the same oligo library. We were also able to generate double-stranded labelled oligos, which produced much stronger signals than single-stranded labelled oligos, by amplification using fluorophore-conjugated primer pairs. Oligos covering cucumber chromosome 1 (Chr1) and chromosome 4 (Chr4) consisting of eight segments were synthesized in one library. Different oligo probes generated from the library painted the corresponding chromosomes/segments unambiguously, especially on pachytene chromosomes. This technique was then applied to study the homoeologous relationships among cucumber, C. hystrix and C. melo chromosomes based on cross-species chromosome painting using Chr4 probes. We demonstrated that the probe was feasible to detect interspecies chromosome homoeologous relationships and chromosomal rearrangement events. Based on its advantages and great convenience, we anticipate that this flexible oligo-painting technique has great potential for the studies of the structure, organization, and evolution of chromosomes in any species with a sequenced genome.

Quantitative Trait Loci for Seed Size Variation in Cucurbits - A Review

Cucurbits (Cucurbitaceae family) include many economically important fruit vegetable crops such as watermelon, pumpkin/squash, cucumber, and melon. Seed size (SS) is an important trait in cucurbits breeding, which is controlled by quantitative trait loci (QTL). Recent advances have deciphered several signaling pathways underlying seed size variation in model plants such as Arabidopsis and rice, but little is known on the genetic basis of SS variation in cucurbits. Here we conducted literature review on seed size QTL identified in watermelon, pumpkin/squash, cucumber and melon, and inferred 14, 9 and 13 consensus SS QTL based on their physical positions in respective draft genomes. Among them, four from watermelon (ClSS2.2, ClSS6.1, ClSS6.2, and ClSS8.2), two from cucumber (CsSS4.1 and CsSS5.1), and one from melon (CmSS11.1) were major-effect, stable QTL for seed size and weight. Whole genome sequence alignment revealed that these major-effect QTL were located in syntenic regions across different genomes suggesting possible structural and functional conservation of some important genes for seed size control in cucurbit crops. Annotation of genes in the four watermelon consensus SS QTL regions identified genes that are known to play important roles in seed size control including members of the zinc finger protein and the E3 ubiquitin-protein ligase families. The present work highlights the utility of comparative analysis in understanding the genetic basis of seed size variation, which may help future mapping and cloning of seed size QTL in cucurbits.

A Major QTL Located in Chromosome 8 of Cucurbita moschata Is Responsible for Resistance to Tomato Leaf Curl New Delhi Virus

Tomato leaf curl New Delhi virus (ToLCNDV) is a bipartite whitefly transmitted begomovirus, responsible since 2013 of severe damages in cucurbit crops in Southeastern Spain. Zucchini (Cucurbita pepo) is the most affected species, but melon (Cucumis melo) and cucumber (Cucumis sativus) are also highly damaged by the infection. The virus has spread across Mediterranean basin and European countries, and integrated control measures are not being enough to reduce economic losses. The identification of resistance genes is required to develop resistant cultivars. In this assay, we studied the inheritance of the resistance to ToLCNDV previously identified in two Cucurbita moschata accessions. We generated segregating populations crossing both resistant pumpkins, an American improved cultivar Large Cheese (PI 604506) and an Indian landrace (PI 381814), with a susceptible C. moschata genotype (PI 419083). The analysis of symptoms and viral titers of all populations established the same monogenic recessive genetic control in both resistant accessions, and the allelism tests suggest the occurrence of alleles of the same locus. By genotyping with a single nucleotide polymorphism (SNP) collection evenly distributed along the C. moschata genome, a major quantitative trait locus (QTL) was identified in chromosome 8 controlling resistance to ToLCNDV. This major QTL was also confirmed in the interspecific C. moschata × C. pepo segregating populations, although C. pepo genetic background affected the resistance level. Molecular markers here identified, linked to the ToLCNDV resistance locus, are highly valuable for zucchini breeding programs, allowing the selection of improved commercial materials. The duplication of the candidate region within the C. moschata genome was studied, and genes with paralogs or single-copy genes were identified. Its synteny with the region of chromosome 17 of the susceptible C. pepo revealed an INDEL including interesting candidate genes. The chromosome 8 candidate region of C. moschata was also syntenic to the region in chromosome 11 of melon, previously described as responsible of ToLCNDV resistance. Common genes in the candidate regions of both cucurbits, with high- or moderate-impact polymorphic SNPs between resistant and susceptible C. moschata accessions, are interesting to study the mechanisms involved in this recessive resistance.

Genetic architecture of fruit size and shape variation in cucurbits: a comparative perspective

The Cucurbitaceae family hosts many economically important fruit vegetables (cucurbits) such as cucumber, melon, watermelon, pumpkin/squash, and various gourds. The cucurbits are probably best known for the diverse fruit sizes and shapes, but little is known about their genetic basis and molecular regulation. Here, we reviewed the literature on fruit size (FS), shape (FSI), and fruit weight (FW) QTL identified in cucumber, melon, and watermelon, from which 150 consensus QTL for these traits were inferred. Genome-wide survey of the three cucurbit genomes identified 253 homologs of eight classes of fruit or grain size/weight-related genes cloned in Arabidopsis, tomato, and rice that encode proteins containing the characteristic CNR (cell number regulator), CSR (cell size regulator), CYP78A (cytochrome P450), SUN, OVATE, TRM (TONNEAU1 Recruiting Motif), YABBY, and WOX domains. Alignment of the consensus QTL with candidate gene homologs revealed widespread structure and function conservation of fruit size/shape gene homologs in cucurbits, which was exemplified with the fruit size/shape candidate genes CsSUN25-26-27a and CsTRM5 in cucumber, CmOFP1a in melon, and ClSUN25-26-27a in watermelon. In cucurbits, the andromonoecy (for 1-aminocyclopropane-1-carboxylate synthase) and the carpel number (for CLAVATA3) loci are known to have pleiotropic effects on fruit shape, which may complicate identification of fruit size/shape candidate genes in these regions. The present work illustrates the power of comparative analysis in understanding the genetic architecture of fruit size/shape variation, which may facilitate QTL mapping and cloning for fruit size-related traits in cucurbits. The limitations and perspectives of this approach are also discussed.

J, A. S, Bhat KV, Gaikwad AB. Genome wide identification and characterization of microsatellite markers in black pepper (Piper nigrum): A valuable resource for boosting genomics applications

Black pepper is one of the most valued and widely used spices in the world and dominates multi-billion dollar global spices trade. India is amongst the major producers, consumers and exporters of black pepper. In spite of its commercial and cultural importance, black pepper has received meagre attention in terms of generation of genomic resources. Availability of markers distributed throughout the genome would facilitate and accelerate genetic studies, QTL identification, genetic enhancement and crop improvement in black pepper. In this perspective, the sequence information from the recently sequenced black pepper (Piper nigrum) genome has been used for identification and characterisation of Simple Sequence Repeats (SSRs). Total 69,126 SSRs were identified from assembled genomic sequence of P. nigrum. The SSR frequency was 158 per MB making it, one SSR for every 6.3 kb in the assembled genome. Among the different types of microsatellite repeat motifs, dinucleotides were the most abundant (48.6%), followed by trinucleotide (23.7%) and compound repeats (20.62%). A set of 85 SSRs were used for validation, of which 74 produced amplification products of expected size. Genetic diversity of 30 black pepper accessions using 50 SSRs revealed four distinct clusters. Further, the cross species transferability of the SSRs was checked in nine other Piper species. Out of 50 SSRs used, 19 and 31 SSRs were amplified in nine and seven species, respectively. Thus the identified SSRs may have application in other species of the genus Piper where genome sequence is not available yet. Present study reports the first NGS based genomic SSRs in black pepper and thus constitute a valuable resource for a whole fleet of applications in genetics and plant breeding studies such as genetic map construction, QTL identification, map-based gene cloning, marker-assisted selection and evolutionary studies in Piper nigrum and related species.

Chromosomal dynamics in space and time: evolutionary history of Mycetophylax ants across past climatic changes in the Brazilian Atlantic coast

Fungus-farming ants of the genus Mycetophylax exhibit intra and interspecific chromosome variability, which makes them suitable for testing hypotheses about possible chromosomal rearrangements that endure lineage diversification. We combined cytogenetic and molecular data from Mycetophylax populations from coastal environments to trace the evolutionary history of the clade in light of chromosomal changes under a historical and geographic context. Our cytogenetic analyses revealed chromosomal differences within and among species. M. morschi exhibited three distinct karyotypes and considerable variability in the localization of 45S rDNA clusters. The molecular phylogeny was congruent with our cytogenetic findings. Biogeographical and divergence time dating analyses estimated that the most recent common ancestor of Mycetophylax would have originated at about 30 Ma in an area including the Amazon and Southern Grasslands, and several dispersion and vicariance events may have occurred before the colonization of the Brazilian Atlantic coast. Diversification of the psammophilous Mycetophylax first took place in the Middle Miocene (ca. 18-10 Ma) in the South Atlantic coast, while "M. morschi" lineages diversified during the Pliocene-Pleistocene transition (ca. 3-2 Ma) through founder-event dispersal for the Northern coastal regions. Psammophilous Mycetophylax diversification fits into the major global climatic events that have had a direct impact on the changes in sea level as well as deep ecological impact throughout South America. We assume therefore that putative chromosomal rearrangements correlated with increased ecological stress during the past climatic transitions could have intensified and/or accompanied the divergence of the psammophilous Mycetophylax. We further reiterate that "M. morschi" comprises a complex of at least three well-defined lineages, and we emphasize the role of this integrative approach for the identification and delimitation of evolutionary lineages.

Genetic and epigenetic effects on centromere establishment

Endogenous chromosomes contain centromeres to direct equal chromosomal segregation in mitosis and meiosis. The location and function of existing centromeres is usually maintained through cell cycles and generations. Recent studies have investigated how the centromere-specific histone H3 variant CENP-A is assembled and replenished after DNA replication to epigenetically propagate the centromere identity. However, existing centromeres occasionally become inactivated, with or without change in underlying DNA sequences, or lost after chromosomal rearrangements, resulting in acentric chromosomes. New centromeres, known as neocentromeres, may form on ectopic, non-centromeric chromosomal regions to rescue acentric chromosomes from being lost, or form dicentric chromosomes if the original centromere is still active. In addition, de novo centromeres can form after chromatinization of purified DNA that is exogenously introduced into cells. Here, we review the phenomena of naturally occurring and experimentally induced new centromeres and summarize the genetic (DNA sequence) and epigenetic features of these new centromeres. We compare the characteristics of new and native centromeres to understand whether there are different requirements for centromere establishment and propagation. Based on our understanding of the mechanisms of new centromere formation, we discuss the perspectives of developing more stably segregating human artificial chromosomes to facilitate gene delivery in therapeutics and research.

Comparatively Barcoded Chromosomes of Brachypodium Perennials Tell the Story of Their Karyotype Structure and Evolution

The Brachypodium genus is an informative model system for studying grass karyotype organization. Previous studies of a limited number of species and reference chromosomes have not provided a comprehensive picture of the enigmatic phylogenetic relationships in the genus. Comparative chromosome barcoding, which enables the reconstruction of the evolutionary history of individual chromosomes and their segments, allowed us to infer the relationships between putative ancestral karyotypes of extinct species and extant karyotypes of current species. We used over 80 chromosome-specific BAC (bacterial artificial chromosome) clones derived from five reference chromosomes of B. distachyon as probes against the karyotypes of twelve accessions representing five diploid and polyploid Brachypodium perennials. The results showed that descending dysploidy is common in Brachypodium and occurs primarily via nested chromosome fusions. Brachypodiumdistachyon was rejected as a putative ancestor for allotetraploid perennials and B. stacei for B. mexicanum. We propose two alternative models of perennial polyploid evolution involving either the incorporation of a putative x = 5 ancestral karyotype with different descending dysploidy patterns compared to B. distachyon chromosomes or hybridization of two x = 9 ancestors followed by genome doubling and descending dysploidy. Details of the karyotype structure and evolution in several Brachypodium perennials are revealed for the first time.

Oligo-painting and GISH reveal meiotic chromosome biases and increased meiotic stability in synthetic allotetraploid Cucumis ×hytivus with dysploid parental karyotypes

BACKGROUND

Meiosis of newly formed allopolyploids frequently encounter perturbations induced by the merging of divergent and hybridizable genomes. However, to date, the meiotic properties of allopolyploids with dysploid parental karyotypes have not been studied in detail. The allotetraploid Cucumis ×hytivus (HHCC, 2n = 38) was obtained from interspecific hybridization between C. sativus (CC, 2n = 14) and C. hystrix (HH, 2n = 24) followed by chromosome doubling. The results of this study thus offer an excellent opportunity to explore the meiotic properties of allopolyploids with dysploid parental karyotypes.

RESULTS

In this report, we describe the meiotic properties of five chromosomes (C5, C7, H1, H9 and H10) and two genomes in interspecific hybrids and C. ×hytivus (the 4th and 14th inbred family) through oligo-painting and genomic in situ hybridization (GISH). We show that 1) only two translocations carrying C5-oligo signals were detected on the chromosomes C2 and C4 of one 14th individual by the karyotyping of eight 4th and 36 14th plants based on C5- and C7-oligo painting, and possible cytological evidence was observed in meiosis of the 4th generation; 2) individual chromosome have biases for homoeologous pairing and univalent formation in F₁ hybrids and allotetraploids; 3) extensive H-chromosome autosyndetic pairings (e.g., H-H, 25.5% PMCs) were observed in interspecific F₁ hybrid, whereas no C-chromosome autosyndetic pairings were observed (e.g. C-C); 4) the meiotic properties of two subgenomes have significant biases in allotetraploids: H-subgenome exhibits higher univalent and chromosome lagging frequencies than C-subgenome; and 5) increased meiotic stability in the S₁₄ generation compared with the S₄ generation, including synchronous meiosis behavior, reduced incidents of univalent and chromosome lagging.

CONCLUSIONS

These results suggest that the meiotic behavior of two subgenomes has dramatic biases in response to interspecific hybridization and allopolyploidization, and the meiotic behavior harmony of subgenomes is a key subject of meiosis evolution in C. ×hytivus. This study helps to elucidate the meiotic properties and evolution of nascent allopolyploids with the dysploid parental karyotypes.

Hybridization-facilitated genome merger and repeated chromosome fusion after 8 million years

The small genus Ricotia (nine species, Brassicaceae) is confined to the eastern Mediterranean. By comparative chromosome painting and a dated multi-gene chloroplast phylogeny, we reconstructed the origin and subsequent evolution of Ricotia. The ancestral Ricotia genome originated through hybridization between two older genomes with n = 7 and n = 8 chromosomes, respectively, on the Turkish mainland during the Early Miocene (c. 17.8 million years ago, Ma). Since then, the allotetraploid (n = 15) genome has been altered by two independent descending dysploidies (DD) to n = 14 in Ricotia aucheri and the Tenuifolia clade (2 spp.). By the Late Miocene (c. 10 Ma), the latter clade started to evolve in the most diverse Ricotia core clade (6 spp.), the process preceded by a DD event to n = 13. It is noteworthy that this dysploidy was mediated by a unique chromosomal rearrangement, merging together the same two chromosomes as were merged during the origin of a fusion chromosome within the paternal n = 7 genome c. 20 Ma. This shows that within a time period of c. 8 Myr genome evolution can repeat itself and that structurally very similar chromosomes may originate repeatedly from the same ancestral chromosomes by different pathways (end-to-end translocation versus nested chromosome insertion).

The USDA cucumber (Cucumis sativus L.) collection: genetic diversity, population structure, genome-wide association studies, and core collection development

Germplasm collections are a crucial resource to conserve natural genetic diversity and provide a source of novel traits essential for sustained crop improvement. Optimal collection, preservation and utilization of these materials depends upon knowledge of the genetic variation present within the collection. Here we use the high-throughput genotyping-by-sequencing (GBS) technology to characterize the United States National Plant Germplasm System (NPGS) collection of cucumber (Cucumis sativus L.). The GBS data, derived from 1234 cucumber accessions, provided more than 23 K high-quality single-nucleotide polymorphisms (SNPs) that are well distributed at high density in the genome (~1 SNP/10.6 kb). The SNP markers were used to characterize genetic diversity, population structure, phylogenetic relationships, linkage disequilibrium, and population differentiation of the NPGS cucumber collection. These results, providing detailed genetic analysis of the U.S. cucumber collection, complement NPGS descriptive information regarding geographic origin and phenotypic characterization. We also identified genome regions significantly associated with 13 horticulturally important traits through genome-wide association studies (GWAS). Finally, we developed a molecularly informed, publicly accessible core collection of 395 accessions that represents at least 96% of the genetic variation present in the NPGS. Collectively, the information obtained from the GBS data enabled deep insight into the diversity present and genetic relationships among accessions within the collection, and will provide a valuable resource for genetic analyses, gene discovery, crop improvement, and germplasm preservation.

The Evolutionary Consequences of Transposon-Related Pericentromer Expansion in Melon

Transposable elements (TEs) are a major driver of plant genome evolution. A part from being a rich source of new genes and regulatory sequences, TEs can also affect plant genome evolution by modifying genome size and shaping chromosome structure. TEs tend to concentrate in heterochromatic pericentromeric regions and their proliferation may expand these regions. Here, we show that after the split of melon and cucumber, TEs have expanded the pericentromeric regions of melon chromosomes that, probably as a consequence, show a very low recombination frequency. In contrast, TEs have not proliferated to a high extent in cucumber, which has small TE-dense pericentromeric regions and shows a relatively constant recombination rate along chromosomes. These differences in chromosome structure also translate in differences in gene nucleotide diversity. Although gene nucleotide diversity is essentially constant along cucumber chromosomes, melon chromosomes show a bimodal pattern of genetic variability, with a gene-poor region where variability is negatively correlated with gene density. Interestingly, genes are not homogeneously distributed in melon, and the high variable low-recombining pericentromeric regions show a higher concentration of melon-specific genes whereas genes shared with cucumber and other plants are essentially found in gene-rich chromosomal arms. The results presented here suggest that melon pericentromeric regions may allow gene sequences to evolve more freely than in other chromosomal compartments which may allow new ORFs to arise and eventually be selected. These results show that TEs can drastically change the structure of chromosomes creating different chromosomal compartments imposing different constraints for gene evolution.

Chromosome identification in Cucumis anguria revealed by cross-species single-copy gene FISH

Cucumis anguria is a potential genetic resource for improving crops of the genus Cucumis, owing to its broad-spectrum resistance. However, few cytogenetic studies on C. anguria have been reported because of its small metaphase chromosomes and the scarcity of distinguished chromosomal landmarks. In this study, 14 single-copy genes from cucumber and rDNAs were used as probes for FISH to identify the individual chromosomes of C. anguria. The distinctive signal distribution patterns of the probes allowed us to distinguish each chromosome of C. anguria (A01–A12). Further, detailed chromosome characteristics were obtained through pachytene chromosome FISH. The lengths of pachytene chromosomes varied from 54.80 to 143.41 μm. The proportion of heterochromatin regions varied from 13.56% to 63.86%. Finally, the chromosomal homeologous relationship between C. anguria and cucumber (C1–C7) was analyzed. The results showed that A06 + A09, A03 + A12, A02 + A04, and A01 + A11 were homeologs of C1, C2, C3, and C6, respectively. Furthemore, chromosomes A08, A10, and A05 were homeologs of C4, C5, and C7, respectively. Chromosome identification and homeologous relationship analysis between C. anguria and cucumber lay the foundation for further research of genome structure evolution in species of Cucumis.

GLABROUS (CmGL) encodes a HD-ZIP IV transcription factor playing roles in multicellular trichome initiation in melon

Map-based cloning identified CmGL that encodes a HD-ZIP type IV transcription factor that controls multicellular trichome initiation in melon. Trichomes are small hairs covering the aerial parts of plants that originate from the epidermal cells, which can protect plants against the damage by insects and pathogens. The regulatory pathway of unicellular trichomes has been well studied in the model plant Arabidopsis. Little is known about the genetic control and regulation of trichome development in melon (Cucumis melo L.) which has multicellular trichomes. In this study, we identified a melon mutant, cmgl, which showed completely glabrous on all aerial organs. A bulked segregant analysis was conducted to identify polymorphic markers for linkage analysis in a population with 256 F₂ plants, which allowed to locate the cmgl locus in melon chromosome VIII. Next-generation sequencing-aided marker discovery and fine mapping in a large population with 1536 F₂ plants narrowed the candidate gene region to 12 kb that harbored only one candidate gene for cmgl, which encoded a class IV homeodomain-associated leucine zipper transcription factor. Four SNPs in the coding region of the CmGL gene were identified between the two parental lines; a single base substitution from C to A resulted in a premature termination codon and a truncated protein in the cmgl. The SNP was converted into a dCAPS marker, which showed co-segregation in the F₂ population and 564 melon accessions. Result of this study will be helpful for better understanding of genetic control of trichome development in melon and marker-assisted selection in developing new cultivars.

Organization and evolution of four differentially amplified tandem repeats in the Cucumis hystrix genome

Three subtelomeric satellites and one interstitial 5S rDNA were characterized in Cucumis hystrix, and the pericentromeric signals of two C. hystrix subtelomeric satellites along C. sativus chromosomes supported the hypothesis of chromosome fusion in Cucumis. Tandem repeats are chromosome structural fractions consisting of highly repetitive sequences organized in large tandem arrays in most eukaryotes. Differentiation of tandem repeats directly affects the chromosome structure, which contributes to species formation and evolution. Cucumis hystrix (2n = 2x = 24) is the only wild Cucumis species grouped into the same subgenus with C. sativus (2n = 2x = 14), hence its phylogenetic position confers a vital role for C. hystrix to understand the chromosome evolution in Cucumis. However, our knowledge of C. hystrix tandem repeats is insufficient for a detailed understanding of the chromosome evolution in Cucumis. Based on de novo tandem repeat characterization using bioinformatics and in situ hybridization (ISH), we identified and characterized four differentially amplified tandem repeats, Cucumis hystrix satellite 1-3 (CuhySat1-CuhySat3) located at the subtelomeric regions of all chromosomes, and Cucumis hystrix 5S (Cuhy5S) located at the interstitial regions of one single chromosome pair. Comparative ISH mapping using CuhySat1-3 and Cuhy5S revealed high homology of tandem repeats between C. hystrix and C. sativus. Intriguingly, we found signal distribution variations of CuhySat2 and CuhySat3 on C. sativus chromosomes. In comparison to their subtelomeric signal distribution on C. hystrix chromosomes, CuhySat3 showed a pericentromeric signal distribution and CuhySat2 showed both subtelomeric and pericentromeric signal distributions on C. sativus chromosomes. This detailed characterization of four C. hystrix tandem repeats significantly widens our knowledge of the C. hystrix chromosome structure, and the observed signal distribution variations will be helpful for understanding the chromosome evolution of Cucumis.

FISH-based mitotic and meiotic diakinesis karyotypes of Morus notabilis reveal a chromosomal fusion-fission cycle between mitotic and meiotic phases

Mulberry (Morus spp.), in family Moraceae, is a plant with important economic value. Many polyploid levels of mulberry have been determined. In the present study, the fluorescence in situ hybridization (FISH) technique was applied in Morus notabilis, using four single-copy sequences, telomere repeats, and 5S and 25S rDNAs as probes. All the mitotic chromosomes were clearly identified and grouped into seven pairs of homologous chromosomes. Three dot chromosome pairs were distinguished by the FISH patterns of the 25S rDNA probe and a simple sequence repeat (SSR2524). According to the FISH signals, chromosome length and morphology, detailed meiotic diakinesis karyotype was constructed. Interestingly, only six bivalent chromosomes were observed in diakinesis cells. The 25S rDNA probe was used to illustrate chromosome alterations. The results indicated that mitotic chromosomes 5 and 7 fused into diakinesis chromosome 5 during the meiotic phase. In mitotic cells, the fused chromosome 5 broke into chromosomes 5 and 7. A chromosomal fusion-fission cycle between the meiotic and mitotic phases in the same individual is reported here for the first time. This finding will contribute to the understanding of karyotype evolution in plants.

Identification of all homoeologous chromosomes of newly synthetic allotetraploid Cucumis × hytivus and its wild parent reveals stable subgenome structure

Allopolyploidy and homoeologous recombination are two important processes in reshaping genomes and generating evolutionary novelties. Newly formed allopolyploids usually display chromosomal perturbations as a result of pairing errors at meiosis. To understand mechanisms of stabilization of allopolyploid species derived from distant chromosome bases, we investigated mitotic stability of a synthetic Cucumis allotetraploid species in relation to meiosis chromosome behavior. The Cucumis × hytivus is an allotetraploid synthesized from interspecific hybridization between cucumber (Cucumis sativus, 2n = 14) and its wild relative Cucumis hystrix (2n = 24) followed by spontaneous chromosome doubling. In the present study, we analyzed the wild parent C. hystrix and the latest generation of C. hytivus using GISH (genomic in situ hybridization) and cross-species FISH (fluorescence in situ hybridization). The karyotype of C. hystrix was constructed with two methods using cucumber fosmid clones and repetitive sequences. Using repeat-element probe mix in two successive hybridizations allowed for routine identification of all 19 homoeologous chromosomes of allotetraploid C. hytivus. No aneuploids were identified in any C. hytivus individuals that were characterized, and no large-scale chromosomal rearrangements were identified in this synthetic allotetraploid. Meiotic irregularities, such as homoeologous pairing, were frequently observed, resulting in univalent and intergenomic multivalent formation. The relatively stable chromosome structure of the synthetic Cucumis allotetraploid may be explained by more deleterious chromosomal viable gametes compared with other allopolyploids. The knowledge of genetic and genomic information of Cucumis allotetraploid species could provide novel insights into the establishment of allopolyploids with different chromosome bases.

Round fruit shape in WI7239 cucumber is controlled by two interacting quantitative trait loci with one putatively encoding a tomato SUN homolog

QTL analysis revealed two interacting loci, FS1.2 and FS2.1, underlying round fruit shape in WI7239 cucumber; CsSUN , a homolog of tomato fruit shape gene SUN , was a candidate for FS1.2. Fruit size is an important quality and yield trait in cucumber, but its genetic basis remains poorly understood. Here we reported QTL mapping results on fruit size with segregating populations derived from the cross between WI7238 (long fruit) and WI7239 (round fruit) inbred cucumber lines. Phenotypic data of fruit length and diameter were collected at anthesis, immature and mature fruit stages in four environments. Ten major-effect QTL were detected for six traits; synthesis of information from these QTL supported two genes, FS1.2 and FS2.1, underlying fruit size variation in the examined populations. Under the two-gene model, deviation from expected segregation ratio in fruit length and diameter among segregating populations was observed, which could be explained mainly by the interactions between FS1.2 and FS2.1, and segregation distortion in the FS2.1 region. Genome-wide candidate gene search identified CsSUN, a homolog of the tomato fruit shape gene SUN, as the candidate for FS1.2. The round-fruited WI7239 had a 161-bp deletion in the first exon of CsSUN, and its expression in WI7239 was significantly lower than that in WI7238. A marker derived from this deletion was mapped at the peak location of FS1.2 in QTL analysis. Comparative analysis suggested the melon gene CmSUN-14, a homolog of CsSUN as a candidate of the fl2/fd2/fw2 QTL in melon. This study revealed the unique genetic architecture of round fruit shape in WI7239 cucumber. It also highlights the power of QTL analysis for traits with a simple genetic basis but their expression is complicated by other factors.

Comparative High-Density Linkage Mapping Reveals Conserved Genome Structure but Variation in Levels of Heterochiasmy and Location of Recombination Cold Spots in the Common Frog

By combining 7077 SNPs and 61 microsatellites, we present the first linkage map for some of the early diverged lineages of the common frog, Rana temporaria, and the densest linkage map to date for this species. We found high homology with the published linkage maps of the Eastern and Western lineages but with differences in the order of some markers. Homology was also strong with the genome of the Tibetan frog Nanorana parkeri and we found high synteny with the clawed frog Xenopus tropicalis. We confirmed marked heterochiasmy between sexes and detected nonrecombining regions in several groups of the male linkage map. Contrary to the expectations set by the male heterogamety of the common frog, we did not find male heterozygosity excess in the chromosome previously shown to be linked to sex determination. Finally, we found blocks of loci showing strong transmission ratio distortion. These distorted genomic regions might be related to genetic incompatibilities between the parental populations, and are promising candidates for further investigation into the genetic basis of speciation and adaptation in the common frog.

A valid name for the Xishuangbanna gourd, a cucumber with carotene-rich fruits

Herbarium specimens deposited in publicly accessible collections are the basis for all scientific names because only permanent specimens can be re-studied by independent researchers, the very essence of science. Re-investigations may be done with morphological, chemical, genomic, computer-tomographic, or other methods. Based on new herbarium material, I here provide a name for the Xishuangbanna gourd, a plant long cultivated in Yunnan because of its large non-bitter fruits, rich in β-carotene. Genome re-sequencing of numerous accessions has shown that this cucumber mutant is closer to Cucumis sativus var. sativus than is the wild bitter-fruited progenitor C. sativus var. hardwickii, and two dozen studies have further clarified the genetics of key traits, including pulp color, fruit shape, and flowering times. Morphological and molecular diagnoses of the new variety are provided and museum-quality specimens have been distributed to the World's major herbaria.

Genome wide characterization of simple sequence repeats in watermelon genome and their application in comparative mapping and genetic diversity analysis

BACKGROUND

Microsatellite markers are one of the most informative and versatile DNA-based markers used in plant genetic research, but their development has traditionally been difficult and costly. The whole genome sequencing with next-generation sequencing (NGS) technologies provides large amounts of sequence data to develop numerous microsatellite markers at whole genome scale. SSR markers have great advantage in cross-species comparisons and allow investigation of karyotype and genome evolution through highly efficient computation approaches such as in silico PCR. Here we described genome wide development and characterization of SSR markers in the watermelon (Citrullus lanatus) genome, which were then use in comparative analysis with two other important crop species in the Cucurbitaceae family: cucumber (Cucumis sativus L.) and melon (Cucumis melo L.). We further applied these markers in evaluating the genetic diversity and population structure in watermelon germplasm collections.

RESULTS

A total of 39,523 microsatellite loci were identified from the watermelon draft genome with an overall density of 111 SSRs/Mbp, and 32,869 SSR primers were designed with suitable flanking sequences. The dinucleotide SSRs were the most common type representing 34.09 % of the total SSR loci and the AT-rich motifs were the most abundant in all nucleotide repeat types. In silico PCR analysis identified 832 and 925 SSR markers with each having a single amplicon in the cucumber and melon draft genome, respectively. Comparative analysis with these cross-species SSR markers revealed complicated mosaic patterns of syntenic blocks among the genomes of three species. In addition, genetic diversity analysis of 134 watermelon accessions with 32 highly informative SSR loci placed these lines into two groups with all accessions of C.lanatus var. citorides and three accessions of C. colocynthis clustered in one group and all accessions of C. lanatus var. lanatus and the remaining accessions of C. colocynthis clustered in another group. Furthermore, structure analysis was consistent with the dendrogram indicating the 134 watermelon accessions were classified into two populations.

CONCLUSION

The large number of genome wide SSR markers developed herein from the watermelon genome provides a valuable resource for genetic map construction, QTL exploration, map-based gene cloning and marker-assisted selection in watermelon which has a very narrow genetic base and extremely low polymorphism among cultivated lines. Furthermore, the cross-species transferable SSR markers identified herein should also have practical uses in many applications in species of Cucurbitaceae family whose whole genome sequences are not yet available.

Chromosome number reduction in the sister clade of Carica papaya with concomitant genome size doubling

PREMISE OF THE STUDY

Caricaceae include six genera and 34 species, among them papaya, a model species in plant sex chromosome research. The family was held to have a conserved karyotype with 2n = 18 chromosomes, an assumption based on few counts. We examined the karyotypes and genome size of species from all genera to test for possible cytogenetic variation.

METHODS

We used fluorescent in situ hybridization using standard telomere, 5S, and 45S rDNA probes. New and published data were combined with a phylogeny, molecular clock dating, and C values (available for ∼50% of the species) to reconstruct genome evolution.

KEY RESULTS

The African genus Cylicomorpha, which is sister to the remaining Caricaceae (all neotropical), has 2n = 18, as do the species in two other genera. A Mexican clade of five species that includes papaya, however, has 2n = 18 (papaya), 2n = 16 (Horovitzia cnidoscoloides), and 2n = 14 (Jarilla caudata and J. heterophylla; third Jarilla not counted), with the phylogeny indicating that the dysploidy events occurred ∼16.6 and ∼5.5 million years ago and that Jarilla underwent genome size doubling (∼450 to 830-920 Mbp/haploid genome). Pericentromeric interstitial telomere repeats occur in both Jarilla adjacent to 5S rDNA sites, and the variability of 5S rDNA sites across all genera is high.

CONCLUSIONS

On the basis of outgroup comparison, 2n = 18 is the ancestral number, and repeated chromosomal fusions with simultaneous genome size increase as a result of repetitive elements accumulating near centromeres characterize the papaya clade. These results have implications for ongoing genome assemblies in Caricaceae.

Chromosome identification for the carnivorous plant Genlisea margaretae

Genlisea margaretae, subgenus Genlisea, section Recurvatae (184 Mbp/1C), belongs to a plant genus with a 25-fold genome size difference and an extreme genome plasticity. Its 19 chromosome pairs could be distinguished individually by an approach combining optimized probe pooling and consecutive rounds of multicolor fluorescence in situ hybridization (mcFISH) with bacterial artificial chromosomes (BACs) selected for repeat-free inserts. Fifty-one BACs were assigned to 18 chromosome pairs. They provide a tool for future assignment of genomic sequence contigs to distinct chromosomes as well as for identification of homeologous chromosome regions in other species of the carnivorous Lentibulariaceae family, and potentially of chromosome rearrangements, in cases where more than one BAC per chromosome pair was identified.

Integrated analysis in bi-parental and natural populations reveals CsCLAVATA3 (CsCLV3) underlying carpel number variations in cucumber

Carpel number variation in cucumber was controlled by a single gene, Cn . Linkage and association analysis revealed CsCLV3 as the candidate gene of the Cn locus. Carpel number (CN) is an important fruit quality trait of cucumber, but the genetic basis of CN variations is largely unknown. In the present study, segregating analysis in multiple bi-parental mapping populations (F2, F3, and RILs) derived from WI2757 (CN = 3) × True Lemon (CN = 5) suggested that CN is controlled by a simply inherited gene, Cn, with CN = 3 being incompletely dominant to CN = 5. Initial linkage mapping located Cn in a 1.9-Mb region of cucumber chromosome 1. Exploration of DNA sequence variations in this region with in silico bulked segregant analysis among eight re-sequenced lines allowed delimiting the Cn locus to a 16-kb region with five predicted genes including CsCLV3, a homolog of the Arabidopsis gene CLAVATA3. Fine genetic mapping in F2 and RIL populations and association analysis in natural populations confirmed CsCLV3 as the candidate gene for Cn, which was further evidenced from gene expression analysis and microscopic examination of floral meristem size in the two parent lines. This study highlights the importance of integrated use of linkage and association analysis as well as next-gen high-throughput sequencing in mapping and cloning genes that are difficult in accurate genotyping. The results provide new insights into the genetic control of CN variations in cucumber, which were discussed in the context of the well-characterized CLAVATA pathway for stem cell homeostasis and regulation of meristem sizes in plants. The associations of carpel number with fruit shape, size, and weight in cucumber and melon are also discussed.

Advancing Eucalyptus Genomics: Cytogenomics Reveals Conservation of Eucalyptus Genomes

The genus Eucalyptus encloses several species with high ecological and economic value, being the subgenus Symphyomyrtus one of the most important. Species such as E. grandis and E. globulus are well characterized at the molecular level but knowledge regarding genome and chromosome organization is very scarce. Here we characterized and compared the karyotypes of three economically important species, E. grandis, E. globulus, and E. calmadulensis, and three with ecological relevance, E. pulverulenta, E. cornuta, and E. occidentalis, through an integrative approach including genome size estimation, fluorochrome banding, rDNA FISH, and BAC landing comprising genes involved in lignin biosynthesis. All karyotypes show a high degree of conservation with pericentromeric 35S and 5S rDNA loci in the first and third pairs, respectively. GC-rich heterochromatin was restricted to the 35S rDNA locus while the AT-rich heterochromatin pattern was species-specific. The slight differences in karyotype formulas and distribution of AT-rich heterochromatin, along with genome sizes estimations, support the idea of Eucalyptus genome evolution by local expansions of heterochromatin clusters. The unusual co-localization of both rDNA with AT-rich heterochromatin was attributed mainly to the presence of silent transposable elements in those loci. The cinnamoyl CoA reductase gene (CCR1) previously assessed to linkage group 10 (LG10) was clearly localized distally at the long arm of chromosome 9 establishing an unexpected correlation between the cytogenetic chromosome 9 and the LG10. Our work is novel and contributes to the understanding of Eucalyptus genome organization which is essential to develop successful advanced breeding strategies for this genus.

Next generation sequencing and omics in cucumber (Cucumis sativus L.) breeding directed research

In the post-genomic era the availability of genomic tools and resources is leading us to novel generation methods in plant breeding, as they facilitate the study of the genotype and its relationship with the phenotype, in particular for complex traits. In this study we have mainly concentrated on the Cucumis sativus and (but much less) Cucurbitaceae family several important vegetable crops. There are many reports on research conducted in Cucurbitaceae plant breeding programs on the ripening process, phloem transport, disease resistance, cold tolerance and fruit quality traits. This paper presents the role played by new omic technologies in the creation of knowledge on the mechanisms of the formation of the breeding features. The analysis of NGS (NGS-next generation sequencing) data allows the discovery of new genes and regulatory sequences, their positions, and makes available large collections of molecular markers. Genome-wide expression studies provide breeders with an understanding of the molecular basis of complex traits. Firstly a high density map should be created for the reference genome, then each re-sequencing data could be mapped and new markers brought out into breeding populations. The paper also presents methods that could be used in the future for the creation of variability and genomic modification of the species in question. It has been shown also the state and usefulness in breeding the chloroplastomic and mitochondriomic study.

Chromosomal structures and repetitive sequences divergence in Cucumis species revealed by comparative cytogenetic mapping

Background

Differentiation and copy number of repetitive sequences affect directly chromosome structure which contributes to reproductive isolation and speciation. Comparative cytogenetic mapping has been verified an efficient tool to elucidate the differentiation and distribution of repetitive sequences in genome. In present study, the distinct chromosomal structures of five Cucumis species were revealed through genomic in situ hybridization (GISH) technique and comparative cytogenetic mapping of major satellite repeats.

Results

Chromosome structures of five Cucumis species were investigated using GISH and comparative mapping of specific satellites. Southern hybridization was employed to study the proliferation of satellites, whose structural characteristics were helpful for analyzing chromosome evolution. Preferential distribution of repetitive DNAs at the subtelomeric regions was found in C. sativus, C hystrix and C. metuliferus, while majority was positioned at the pericentromeric heterochromatin regions in C. melo and C. anguria. Further, comparative GISH (cGISH) through using genomic DNA of other species as probes revealed high homology of repeats between C. sativus and C. hystrix. Specific satellites including 45S rDNA, Type I/II, Type III, Type IV, CentM and telomeric repeat were then comparatively mapped in these species. Type I/II and Type IV produced bright signals at the subtelomeric regions of C. sativus and C. hystrix simultaneously, which might explain the significance of their amplification in the divergence of Cucumis subgenus from the ancient ancestor. Unique positioning of Type III and CentM only at the centromeric domains of C. sativus and C. melo, respectively, combining with unique southern bands, revealed rapid evolutionary patterns of centromeric DNA in Cucumis. Obvious interstitial telomeric repeats were observed in chromosomes 1 and 2 of C. sativus, which might provide evidence of the fusion hypothesis of chromosome evolution from x = 12 to x = 7 in Cucumis species. Besides, the significant correlation was found between gene density along chromosome and GISH band intensity in C. sativus and C. melo.

Conclusions

In summary, comparative cytogenetic mapping of major satellites and GISH revealed the distinct differentiation of chromosome structure during species formation. The evolution of repetitive sequences was the main force for the divergence of Cucumis species from common ancestor.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1877-6) contains supplementary material, which is available to authorized users.