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.     

Single‐copy gene‐based chromosome painting in cucumber and its application for chromosome rearrangement analysis in Cucumis

Chromosome painting based on fluorescence in situ hybridization (FISH) has played an important role in chromosome identification and research into chromosome rearrangements, diagnosis of chromosome abnormalities and evolution in human and animal species. However, it has not been applied widely in plants due to the large amounts of dispersed repetitive sequences in chromosomes. In the present work, a chromosome painting method for single‐copy gene pools in Cucumis sativus was successfully developed. Gene probes with sizes above 2 kb were detected consistently. A cucumber karyotype was constructed based on FISH using a cocktail containing chromosome‐specific gene probes. This single‐copy gene‐based chromosome painting (ScgCP) technique was performed by PCR amplification, purification, pooling, labeling and hybridization onto chromosome spreads. Gene pools containing sequential genes with an interval less than 300 kb yielded painting patterns on pachytene chromosomes. Seven gene pools corresponding to individual chromosomes unambiguously painted each chromosome pair of C. sativus. Three mis‐aligned regions on chromosome 4 were identified by the painting patterns. A probe pool comprising 133 genes covering the 8 Mb distal end of chromosome 4 was used to evaluate the potential utility of the ScgCP technique for chromosome rearrangement research through cross‐species FISH in the Cucumis genus. Distinct painting patterns of this region were observed in C. sativus, C. melo and C. metuliferus species. A comparative chromosome map of this region was constructed between cucumber and melon. With increasing sequence resources, this ScgCP technique may be applied on any other sequenced species for chromosome painting research.     

Chromosome-Specific Painting in Cucumis Species Using Bulked Oligonucleotides

Chromosome-specific painting is a powerful technique in molecular cytogenetic and genome research. We developed an oligonucleotide (oligo)-based chromosome painting technique in cucumber (Cucumis sativus) that will be applicable in any plant species with a sequenced genome. Oligos specific to a single chromosome of cucumber were identified using a newly developed bioinformatic pipeline and then massively synthesized de novo in parallel. The synthesized oligos were amplified and labeled with biotin or digoxigenin for use in fluorescence in situ hybridization (FISH). We developed three different probes with each containing 23,000–27,000 oligos. These probes spanned 8.3–17 Mb of DNA on targeted cucumber chromosomes and had the densities of 1.5–3.2 oligos per kilobases. These probes produced FISH signals on a single cucumber chromosome and were used to paint homeologous chromosomes in other Cucumis species diverged from cucumber for up to 12 million years. The bulked oligo probes allowed us to track a single chromosome in early stages during meiosis. We were able to precisely map the pairing between cucumber chromosome 7 and chromosome 1 of Cucumis hystrix in a F₁ hybrid. These two homeologous chromosomes paired in 71% of prophase I cells but only 25% of metaphase I cells, which may provide an explanation of the higher recombination rates compared to the chiasma frequencies between homeologous chromosomes reported in plant hybrids.     

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.     

Phylogenetics of <Emphasis Type="Italic">Cucumis</Emphasis> (Cucurbitaceae): Cucumber (<Emphasis Type="Italic">C. sativus</Emphasis>) belongs in an Asian/Australian clade far from melon (<Emphasis Type="Italic">C. melo</Emphasis>)

Background  

Melon, Cucumis melo, and cucumber, C. sativus, are among the most widely cultivated crops worldwide. Cucumis, as traditionally conceived, is geographically centered in Africa, with C. sativus and C. hystrix thought to be the only Cucumis species in Asia. This taxonomy forms the basis for all ongoing Cucumis breeding and genomics efforts. We tested relationships among Cucumis and related genera based on DNA sequences from chloroplast gene, intron, and spacer regions (rbcL, matK, rpl20-rps12, trnL, and trnL-F), adding nuclear internal transcribed spacer sequences to resolve relationships within Cucumis.     

Chromosome rearrangements during domestication of cucumber as revealed by high-density genetic mapping and draft genome assembly

Cucumber, Cucumis sativus L. is the only taxon with 2n = 2x = 14 chromosomes in the genus Cucumis. It consists of two cross‐compatible botanical varieties: the cultivated C. sativus var. sativus and the wild C. sativus var. hardwickii. There is no consensus on the evolutionary relationship between the two taxa. Whole‐genome sequencing of the cucumber genome provides a new opportunity to advance our understanding of chromosome evolution and the domestication history of cucumber. In this study, a high‐density genetic map for cultivated cucumber was developed that contained 735 marker loci in seven linkage groups spanning 707.8 cM. Integration of genetic and physical maps resulted in a chromosome‐level draft genome assembly comprising 193 Mbp, or 53% of the 367 Mbp cucumber genome. Strategically selected markers from the genetic map and draft genome assembly were employed to screen for fosmid clones for use as probes in comparative fluorescence in situ hybridization analysis of pachytene chromosomes to investigate genetic differentiation between wild and cultivated cucumbers. Significant differences in the amount and distribution of heterochromatins, as well as chromosomal rearrangements, were uncovered between the two taxa. In particular, six inversions, five paracentric and one pericentric, were revealed in chromosomes 4, 5 and 7. Comparison of the order of fosmid loci on chromosome 7 of cultivated and wild cucumbers, and the syntenic melon chromosome I suggested that the paracentric inversion in this chromosome occurred during domestication of cucumber. The results support the sub‐species status of these two cucumber taxa, and suggest that C. sativus var. hardwickii is the progenitor of cultivated cucumber.     

Molecular and cytogenetic analyses provide evidence of the introgression of chromosomal segments from the wild <Emphasis Type="Italic">Cucumis hystrix</Emphasis> into the cultivated cucumber through the bridge of a synthetic allotetraploid

Cucumis × hytivus (2n = 4× = 38) is a synthetic allotetraploid obtained from interspecific hybridization between the cucumber (2n = 2× = 14) and its wild relative C. hystrix (2n = 2× = 24). The synthesis of this species built a bridge for cucumber improvement through gene introgression. Allotriploid and introgression lines (ILs) have previously been produced and characterized with respect to morphology, cytology, and molecular markers. However, no clear evidence of how the chromosomal segments of C. hystrix were introgressed and inherited was found owing to the small size of chromosomes. In the present study, cucumber-C. hystrix introgression lines were developed by backcrossing the allotriploid to North China cucumber breeding line “P01” followed by self-pollination. The introgressed segments of C. hystrix in the ILs were revealed by meiotic pachytene chromosome analysis. Fluorescence in situ hybridization (FISH) was performed on pachytene chromosomes using fosmid clones from cucumber, which confirmed that introgression occurred in the long arm of chromosome 7. Molecular analysis using a set of 53 simple sequence repeats (SSRs) indicated that the chromosomal segments of C. hystrix were introduced into 4 cucumber chromosomes, the short arms of chromosomes 2 and 6, and long arms of chromosomes 3 and 7. The inheritance of alien sequences in the long arm of chromosome 7 was investigated with 21 SSRs in self-pollinated progenies. C. hystrix-specific bands of several SSRs were still present in some individuals, indicating that the introgressed segment was partially preserved. The first unambiguous identification of alien chromosome segments in cucumber ILs using combined molecular cytogenetics could facilitate the determination of effects of wild alleles and promote cucumber improvement.     

Identification and characterization of potential NBS-encoding resistance genes and induction kinetics of a putative candidate gene associated with downy mildew resistance in <Emphasis Type="Italic">Cucumis</Emphasis>

Background  

Due to the variation and mutation of the races of Pseudoperonospora cubensis, downy mildew has in recent years become the most devastating leaf disease of cucumber worldwide. Novel resistance to downy mildew has been identified in the wild Cucumis species, C. hystrix Chakr. After the successful hybridization between C. hystrix and cultivated cucumber (C. sativus L.), an introgression line (IL5211S) was identified as highly resistant to downy mildew. Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes are the largest class of disease resistance genes cloned from plant with highly conserved domains, which can be used to facilitate the isolation of candidate genes associated with downy mildew resistance in IL5211S.     

<Emphasis Type="Italic">Cucumis</Emphasis> monosomic alien addition lines: morphological,cytological, and genotypic analyses

Cucumis hystrix Chakr. (HH, 2n=24), a wild relative of the cultivated cucumber, possesses several potentially valuable disease-resistance and abiotic stress-tolerance traits for cucumber (C. sativus L., CC, 2n=14) improvement. Numerous attempts have been made to transfer desirable traits since the successful interspecific hybridization between C. hystrix and C. sativus, one of which resulted in the production of an allotriploid (HCC, 2n=26: one genome of C. hystrix and two of C. sativus). When this genotype was treated with colchicine to induce polyploidy, two monosomic alien addition lines (MAALs) (plant nos. 87 and 517: 14 CC+1 H, 2n=15) were recovered among 252 viable plants. Each of these plants was morphologically distinct from allotriploids and cultivated cucumbers. Cytogenetic and molecular marker analyses were performed to confirm the genetic constitution and further characterize these two MAALs. Chromosome counts made from at least 30 meristematic cells from each plant confirmed 15 nuclear chromosomes. In pollen mother cells of plant nos. 87 and 517, seven bivalents and one univalent were observed at diakinesis and metaphase I; the frequency of trivalent formation was low (about 4–5%). At anaphase I and II, stochastic and asymmetric division led to the formation of two gamete classes: n=7 and n=8; however, pollen fertility was relatively high. Pollen stainability in plant no. 87 was 86.7% and in plant no. 517 was 93.2%. Random amplified polymorphic DNA analysis was performed using 100 random 10-base primers. Genotypes obtained with eight primers (A-9, A-11, AH-13, AI-19, AJ-18, AJ-20, E-19, and N-20) showed a band common to the two MAAL plants and C. hystrix that was absent in C. sativus, confirming that the alien chromosomes present in the MAALs were derived from C. hystrix. Morphological differences and differences in banding patterns were also observed between plant nos. 87 and 517 after amplification with primers AI-5, AJ-13, N-12, and N-20, suggesting that these plants may contain different C. hystrix chromosomes.Communicated by H. Nybom     

Sugarcane genome architecture decrypted with chromosome‐specific oligo probes

Sugarcane (Saccharum spp.) is probably the crop with the most complex genome. Modern cultivars (2n = 100–120) are highly polyploids and aneuploids derived from interspecific hybridization between Saccharum officinarum (2n = 80) and Saccharum spontaneum (2n = 40–128). Chromosome‐specific oligonucleotide probes were used in combination with genomic in situ hybridization to analyze the genome architecture of modern cultivars and representatives of their parental species. The results validated a basic chromosome number of x = 10 for S. officinarum. In S. spontaneum, rearrangements occurred from a basic chromosome of x = 10, probably in the Northern part of India, in two steps leading to x = 9 and then x = 8. Each step involved three chromosomes that were rearranged into two. Further polyploidization led to the wide geographical extension of clones with x = 8. We showed that the S. spontaneum contribution to modern cultivars originated from cytotypes with x = 8 and varied in proportion between cultivars (13–20%). Modern cultivars had mainly 12 copies for each of the first four basic chromosomes, and a more variable number for those basic chromosomes whose structure differs between the two parental species. One?four of these copies corresponded to entire S. spontaneum chromosomes or interspecific recombinant chromosomes. In addition, a few inter‐chromosome translocations were revealed. The new information and cytogenetic tools described in this study substantially improve our understanding of the extreme level of complexity of modern sugarcane cultivar genomes.     

Identification of all homoeologous chromosomes of newly synthetic allotetraploid <Emphasis Type="Italic">Cucumis</Emphasis> × <Emphasis Type="Italic">hytivus</Emphasis> 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.     

Chromosomal‐level genomes of three rice planthoppers provide new insights into sex chromosome evolution

The brown planthopper Nilaparvata lugens, white‐backed planthopper Sogatella furcifera, and small brown planthopper Laodelphax striatellus are three major insect pests of rice. They are genetically close; however, they differ in several ecological traits such as host range, migration capacity, and in their sex chromosomes. Though the draft genome of these three planthoppers have been previously released, the quality of genome assemblies need to be improved. The absence of chromosome‐level genome resources has hindered in‐depth research of these three species. Here, we performed a de novo genome assembly for N. lugens to increase its genome assembly quality with PacBio and Illumina platforms, increasing the contig N50 to 589.46 Kb. Then, with the new N. lugens genome and previously reported S. furcifera and L. striatellus genome assemblies, we generated chromosome‐level scaffold assemblies of these three planthopper species using HiC scaffolding technique. The scaffold N50s significantly increased to 77.63 Mb, 43.36 Mb and 29.24 Mb for N. lugens, S. furcifera and L. striatellus, respectively. To identify sex chromosomes of these three planthopper species, we carried out genome re‐sequencing of males and females and successfully determined the X and Y chromosomes for N. lugens, and X chromosome for S. furcifera and L. striatellus. The gene content of the sex chromosomes showed high diversity among these three planthoppers suggesting the rapid evolution of sex‐linked genes, and all chromosomes showed high synteny. The chromosome‐level genome assemblies of three planthoppers would provide a valuable resource for a broad range of future research in molecular ecology, and subsequently benefits development of modern pest control strategies.     

XY<Subscript>1</Subscript>Y<Subscript>2</Subscript> chromosome system in <Emphasis Type="Italic">Salinomys delicatus</Emphasis> (Rodentia,Cricetidae)

Salinomys delicatus is considered a rare species due to its restricted and patchy distribution, poor records and low abundances. It is also the phyllotine with the lowest known diploid chromosome number (2n = 18), however its sex chromosome system has never been described. Here, we studied the chromosomes of six females and three males with bands G, C, DAPI/CMA₃ and meiosis. In males, the chromosome number was 2n = 19, with one large metacentric X-chromosome and two medium-sized acrocentrics absent in females. The karyotype of females was the same as previously described (2n = 18, FN = 32), with X-chromosomes being metacentric and the largest elements of the complement. In males, the two acrocentrics and the large metacentric form a trivalent in meiotic prophase. This indicates that S. delicatus has XY₁Y₂ sex chromosomes, which is confirmed by G and DAPI bands. Constitutive heterochromatin (CH) is restricted to small pericentromeric blocks in all chromosomes. The X-chromosome shows the largest block of centromeric CH, which could favor the establishment of this X-autosome translocation. This sex chromosome system is rare in mammals and, compared with other phyllotine rodents, S. delicatus seems to have undergone a major chromosome restructuring during its karyotypic evolution.     

Production and characterization of intergeneric somatic hybrids between <Emphasis Type="Italic">Brassica napus</Emphasis> and <Emphasis Type="Italic">Orychophragmus violaceus</Emphasis> and their backcrossing progenies

Alien chromosome addition lines have been widely used for identifying gene linkage groups, assigning species-specific characters to a particular chromosome and comparing gene synteny between related species. In plant breeding, their utilization lies in introgressing characters of agronomic value. The present investigation reports the production of intergeneric somatic hybrids Brassica napus (2n = 38) + Orychophragmus violaceus (2n = 24) through asymmetric fusions of mesophyll protoplasts and subsequent development of B. napus-O. violaceous chromosome addition lines. Somatic hybrids showed variations in morphology and fertility and were mixoploids (2n = 51–67) with a range of 19–28 O. violaceus chromosomes identified by genomic in situ hybridization (GISH). After pollinated with B. napus parent and following embryo rescue, 20 BC₁ plants were obtained from one hybrid. These exhibited typical serrated leaves of O. violaceus or B. napus-type leaves. All BC₁ plants were partially male fertile but female sterile because of abnormal ovules. These were mixoploids (2n = 41–54) with 9–16 chromosomes from O. violaceus. BC₂ plants showed segregations for female fertility, leaf shape and still some chromosome variation (2n = 39–43) with 2–5 O. violaceus chromosomes, but mainly containing the whole complement from B. napus. Among the selfed progenies of BC₂ plants, monosomic addition lines (2n = 39, AACC + 1O) with or without the serrated leaves of O. violaceus or female sterility were established. The complete set of additions is expected from this investigation. In addition, O. violaceus plants at diploid and tetraploid levels with some variations in morphology and chromosome numbers were regenerated from the pretreated protoplasts by iodoacetate and UV-irradiation. Z. Zhao and T. Hu make equal contributions to this work.     

Fissions,fusions, and translocations shaped the karyotype and multiple sex chromosome constitution of the northeast‐Asian wood white butterfly,Leptidea amurensis

Previous studies have shown a dynamic karyotype evolution and the presence of complex sex chromosome systems in three cryptic Leptidea species from the Western Palearctic. To further explore the chromosomal particularities of Leptidea butterflies, we examined the karyotype of an Eastern Palearctic species, Leptidea amurensis. We found a high number of chromosomes that differed between the sexes and slightly varied in females (i.e. 2n = 118–119 in females and 2n = 122 in males). The analysis of female meiotic chromosomes revealed multiple sex chromosomes with three W and six Z chromosomes. The curious sex chromosome constitution [i.e. W_1–3/Z_1–6 (females) and Z_1–6/Z_1–6 (males)] and the observed heterozygotes for a chromosomal fusion are together responsible for the sex‐specific and intraspecific variability in chromosome numbers. However, in contrast to the Western Palearctic Leptidea species, the single chromosomal fusion and static distribution of cytogenetic markers (18S rDNA and H3 histone genes) suggest that the karyotype of L. amurensis is stable. The data obtained for four Leptidea species suggest that the multiple sex chromosome system, although different among species, is a common feature of the genus Leptidea. Furthermore, inter‐ and intraspecific variations in chromosome numbers and the complex meiotic pairing of these multiple sex chromosomes indicate the role of chromosomal fissions, fusions, and translocations in the karyotype evolution of Leptidea butterflies.     

Improving Illumina assemblies with Hi‐C and long reads: An example with the North African dromedary

Researchers have assembled thousands of eukaryotic genomes using Illumina reads, but traditional mate‐pair libraries cannot span all repetitive elements, resulting in highly fragmented assemblies. However, both chromosome conformation capture techniques, such as Hi‐C and Dovetail Genomics Chicago libraries and long‐read sequencing, such as Pacific Biosciences and Oxford Nanopore, help span and resolve repetitive regions and therefore improve genome assemblies. One important livestock species of arid regions that does not have a high‐quality contiguous reference genome is the dromedary (Camelus dromedarius). Draft genomes exist but are highly fragmented, and a high‐quality reference genome is needed to understand adaptation to desert environments and artificial selection during domestication. Dromedaries are among the last livestock species to have been domesticated, and together with wild and domestic Bactrian camels, they are the only representatives of the Camelini tribe, which highlights their evolutionary significance. Here we describe our efforts to improve the North African dromedary genome. We used Chicago and Hi‐C sequencing libraries from Dovetail Genomics to resolve the order of previously assembled contigs, producing almost chromosome‐level scaffolds. Remaining gaps were filled with Pacific Biosciences long reads, and then scaffolds were comparatively mapped to chromosomes. Long reads added 99.32 Mbp to the total length of the new assembly. Dovetail Chicago and Hi‐C libraries increased the longest scaffold over 12‐fold, from 9.71 Mbp to 124.99 Mbp and the scaffold N50 over 50‐fold, from 1.48 Mbp to 75.02 Mbp. We demonstrate that Illumina de novo assemblies can be substantially upgraded by combining chromosome conformation capture and long‐read sequencing.     

Genome‐wide association study for pigmentation traits in Chinese Holstein population

With the Illumina BovineSNP50K BeadChip, we performed a genome‐wide association study (GWAS) for two pigmentation traits in a Chinese Holstein population: proportion of black (PB) and teat colour (TC). A case–control design was used. Cases were the cows with PB <0.30 (n = 129) and TC <2 points (n = 140); controls were those with PB >0.90 (n = 58) and TC >4 points (n = 281). The RM test of roadtrips (version 1.2) was applied to detect SNPs for the two traits with 42 883 and 42 741 SNPs respectively. A total of nine and 12 genome‐wide significant (P < 0.05) SNPs associated with PB and TC respectively were identified. Of these, two SNPs for PB were located within the KIT and IGFBP7 genes, and the other four SNPs were 23~212 kb away from the PDGFRA gene on BTA6; nine SNPs associated with TC were located within or 21~78.8 kb away from known genes on chromosomes 4, 11, 22, 23 and 24. By combing through our GWAS results and the biological functions of the genes, we suggest that the KIT, IGFBP7, PDGFRA, MITF, ING3 and WNT16 genes are promising candidates for PB and TC in Holstein cattle, providing a basis for further investigation on the genetic mechanism of pigmentation formation.