Chromosome-level genome assembly of milk thistle (Silybum marianum (L.) Gaertn.).

Silybum marianum (L.) Gaertn., commonly known as milk thistle, is a medicinal plant belonging to the Asteraceae family. This plant has been recognized for its medicinal properties for over 2,000 years. However, the genome of this plant remains largely undiscovered, having no reference genome at a chromosomal level. Here, we assembled the chromosome-level genome of S. marianum, allowing for the annotation of 53,552 genes and the identification of transposable elements comprising 58% of the genome. The genome assembly from this study showed 99.1% completeness as determined by BUSCO assessment, while the previous assembly (ASM154182v1) showed 36.7%. Functional annotation of the predicted genes showed 50,329 genes (94% of total genes) with known protein functions in public databases. Comparative genome analysis among Asteraceae plants revealed a striking conservation of collinearity between S. marianum and C. cardunculus. The genomic information generated from this study will be a valuable resource for milk thistle breeding and for use by the larger research community.

Repeated co-option of HMG-box genes for sex determination in brown algae and animals.

In many eukaryotes, genetic sex determination is not governed by XX/XY or ZW/ZZ systems but by a specialized region on the poorly studied U (female) or V (male) sex chromosomes. Previous studies have hinted at the existence of a dominant male-sex factor on the V chromosome in brown algae, a group of multicellular eukaryotes distantly related to animals and plants. The nature of this factor has remained elusive. Here, we demonstrate that an HMG-box gene acts as the male-determining factor in brown algae, mirroring the role HMG-box genes play in sex determination in animals. Over a billion-year evolutionary timeline, these lineages have independently co-opted the HMG box for male determination, representing a paradigm for evolution's ability to recurrently use the same genetic "toolkit" to accomplish similar tasks.

Chromosome level genome assembly of endangered medicinal plant Anisodus tanguticus.

Anisodus tanguticus is a medicinal herb that belongs to the Anisodus genus of the Solanaceae family. This endangered herb is mainly distributed in Qinghai-Tibet Plateau. In this study, we combined the Illumina short-read, Nanopore long-read and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to de novo assemble the A. tanguticus genome. A high-quality chromosomal-level genome assembly was obtained with a genome size of 1.26 Gb and a contig N50 of 25.07 Mb. Of the draft genome sequences, 97.47% were anchored to 24 pseudochromosomes with a scaffold N50 of 51.28 Mb. In addition, 842.14 Mb of transposable elements occupying 66.70% of the genome assembly were identified and 44,252 protein-coding genes were predicted. The genome assembly of A. tanguticus will provide genetic repertoire to understand the adaptation strategy of Anisodus species in the plateau, which will further promote the conservation of endangered A. tanguticus resources.

Chromosome-level genome assembly of Prunella vulgaris L. provides insights into pentacyclic triterpenoid biosynthesis.

Prunella vulgaris is one of the bestselling and widely used medicinal herbs. It is recorded as an ace medicine for cleansing and protecting the liver in Chinese Pharmacopoeia and has been used as the main constitutions of many herbal tea formulas in China for centuries. It is also a traditional folk medicine in Europe and other countries of Asia. Pentacyclic triterpenoids are a major class of bioactive compounds produced in P. vulgaris. However, their biosynthetic mechanism remains to be elucidated. Here, we report a chromosome-level reference genome of P. vulgaris using an approach combining Illumina, ONT, and Hi-C technologies. It is 671.95 Mb in size with a scaffold N50 of 49.10 Mb and a complete BUSCO of 98.45%. About 98.31% of the sequence was anchored into 14 pseudochromosomes. Comparative genome analysis revealed a recent WGD in P. vulgaris. Genome-wide analysis identified 35 932 protein-coding genes (PCGs), of which 59 encode enzymes involved in 2,3-oxidosqualene biosynthesis. In addition, 10 PvOSC, 358 PvCYP, and 177 PvUGT genes were identified, of which five PvOSCs, 25 PvCYPs, and 9 PvUGTs were predicted to be involved in the biosynthesis of pentacyclic triterpenoids. Biochemical activity assay of PvOSC2, PvOSC4, and PvOSC6 recombinant proteins showed that they were mixed amyrin synthase (MAS), lupeol synthase (LUS), and ß-amyrin synthase (BAS), respectively. The results provide a solid foundation for further elucidating the biosynthetic mechanism of pentacyclic triterpenoids in P. vulgaris.

Haplotype-resolved chromosomal-level genome assembly of Buzhaye (Microcos paniculata).

Microcos paniculata is a shrub used traditionally as folk medicine and to make herbal teas. Previous research into this species has mainly focused on its chemical composition and medicinal value. However, the lack of a reference genome limits the study of the molecular mechanisms of active compounds in this species. Here, we assembled a haplotype-resolved chromosome-level genome of M. paniculata based on PacBio HiFi and Hi-C data. The assembly contains two haploid genomes with sizes 399.43 Mb and 393.10 Mb, with contig N50 lengths of 43.44 Mb and 30.17 Mb, respectively. About 99.93% of the assembled sequences could be anchored to 18 pseudo-chromosomes. Additionally, a total of 482 Mb repeat sequences were identified, accounting for 60.76% of the genome. A total of 49,439 protein-coding genes were identified, of which 48,979 (99%) were functionally annotated. This haplotype-resolved chromosome-level assembly and annotation of M. paniculata will serve as a valuable resource for investigating the biosynthesis and genetic basis of active compounds in this species, as well as advancing evolutionary phylogenomic studies in Malvales.

Chromosome-level genome assembly of bunching onion illuminates genome evolution and flavor formation in Allium crops.

The Allium genus is cultivated globally as vegetables, condiments, or medicinal plants and is characterized by large genomes and strong pungency. However, the genome evolution and genomic basis underlying their unique flavor formation remain poorly understood. Herein, we report an 11.27-Gb chromosome-scale genome assembly for bunching onion (A. fistulosum). The uneven bursts of long-terminal repeats contribute to diversity in genome constituents, and dispersed duplication events largely account for gene expansion in Allium genomes. The extensive duplication and differentiation of alliinase and lachrymatory factor synthase manifest as important evolutionary events during flavor formation in Allium crops. Furthermore, differential selective preference for flavor-related genes likely lead to the variations in isoalliin content in bunching onions. Moreover, we reveal that China is the origin and domestication center for bunching onions. Our findings provide insights into Allium genome evolution, flavor formation and domestication history and enable future genome-assisted breeding of important traits in these crops.

Integrating Genetic and Chromosome Maps of Allium cepa: From Markers Visualization to Genome Assembly Verification.

The ability to directly look into genome sequences has opened great opportunities in plant breeding. Yet, the assembly of full-length chromosomes remains one of the most difficult problems in modern genomics. Genetic maps are commonly used in de novo genome assembly and are constructed on the basis of a statistical analysis of the number of recombinations. This may affect the accuracy of the ordering and orientation of scaffolds within the chromosome, especially in the region of recombination suppression. Moreover, it is impossible to assign contigs lacking DNA markers. Here, we report the use of Tyr-FISH to determine the position of the short DNA sequence of markers and non-mapped unique copy sequence on the physical chromosomes of a large-genome onion (Allium cepa L.). In order to minimize potential background masking of the target signal, we improved our earlier developed pipeline for probe design. A total of 23 markers were located on physical chromosomes 2 and 6. The order of markers was corrected by the integration of genetic, pseudochromosome maps and cytogenetic maps. Additionally, the position of the mlh1 gene, which was not on the genetic map, was defined on physical chromosome 2. Tyr-FISH mapping showed that the order of 23.1% (chromosome 2) and 27.3% (chromosome 6) of the tested genes differed between physical chromosomes and pseudochromosomes. The results can be used for the improvement of pseudochromosome 2 and 6 assembly. The present study aims to demonstrate the value of the in situ visualization of DNA sequences in chromosome-scaffold genome assembly.

Comparative analysis reveals chromosome number reductions in the evolution of African bermudagrass (Cynodon transvaalensis Burtt-Davy).

African bermudagrass (Cynodon transvaalensis Burtt-Davy) (2n = 2x = 18) belongs to the genus Cynodon, tribe Cynodonteae, subfamily Chloridoideae in the grass family Poaceae. The species is frequently crossed with common bermudagrass (Cynodon dactylon Pers.) in developing high-quality hybrid turf cultivars. Molecular resources for C. transvaalensis are scarce; thus, its genomic evolution is unknown. Recently, a linkage map consisting of 1278 markers provided a powerful tool for African bermudagrass genomic research. The objective of this study was to investigate chromosome number reduction events that resulted in the nine haploid chromosomes in this species. Tag sequences of mapped single nucleotide polymorphism markers in C. transvaalensis were compared against genome sequences of Oropetium thomaeum (L.f.) Trin. (2n = 2x = 20), a genomic model in the Cynodonteae tribe. The comparative genomic analyses revealed broad collinearity between the genomes of these two species. The analyses further revealed that two major interchromosomal rearrangements of the paleochromosome ρ12 (ρ1-ρ12-ρ1 and ρ6-ρ12-ρ6) resulted in nine chromosomes in the genome of C. transvaalensis. The findings provide novel information regarding the formation of the initial diploid species in the Cynodon genus.

Chromosome-scale genome assembly of an important medicinal plant honeysuckle.

Lonicera japonica (honeysuckle) is one of the most important medicinal plants and widely utilized in traditional Chinese medicine. At present, there are many varieties of honeysuckle used in cultivation, among which Sijihua variety are widely cultivated due to its wide adaptability, stress resistance, early flowering and high yield. In this study, we assembled the genome of Sijihua, which was approximately 886.04 Mb in size with a scaffold N50 of 79.5 Mb. 93.28% of the total assembled sequences were anchored to 9 pseudo-chromosomes by using PacBio long reads and Hi-C sequencing data. We predicted 39,320 protein-coding genes and 92.87% of them could be annotated in NR, GO, KOG, KEGG and other databases. In addition, we identified 644 tRNAs, 2,156 rRNAs, 109 miRNAs and 5,502 pseudogenes from the genome. The chromosome-scale genome of Sijihua will be a significant resource for understanding the genetic basis of high stress-resistance, which will facilitate further study of the genetic diversity and accelerate the genetic improvement and breeding of L. japonica.

The DcPS1 cooperates with OSDLa during pollen development and 2n gamete production in carnation meiosis.

BACKGROUND: Deciphering the mechanisms of meiosis has important implications for potential applications in plant breeding programmes and species evolution. However, the process of meiosis is poorly understood in carnation, which is famous for its cut flowers. RESULTS: We report that Dianthus caryophyllus parallel spindle 1 (DcPS1) regulates omission of second division like a (OSDLa) during pollen development and 2n gamete production in carnation meiosis. In DcPS1 and OSDLa RNAi lines, an absence of the second meiotic division and the abnormal orientation of spindles at meiosis II might be the main reason for dyad/triad formation, resulting in unreduced gametes. We also found that carnation OSDLa interacted with DcPS1 and DcRAD51D. In the DcPS1 RNAi lines, a decrease in OSDLa and DcRAD51D expression was observed. In the OSDLa RNAi lines, a decrease in DcPS1 and DcRAD51D expression was also observed. We propose that DcPS1 regulates OSDLa expression, allowing entry into meiosis II and the proper orientation of the metaphase II spindle in meiosis II. We also propose that OSDLa regulates DcRAD51D expression, allowing for homologous recombination. CONCLUSIONS: These results suggest a critical role for DcPS1 and OSDLa in diplogamete production during meiosis and open a new pathway for meiosis-related studies.

A Solanum lycopersicoides reference genome facilitates insights into tomato specialized metabolism and immunity.

Wild relatives of tomato are a valuable source of natural variation in tomato breeding, as many can be hybridized to the cultivated species (Solanum lycopersicum). Several, including Solanum lycopersicoides, have been crossed to S. lycopersicum for the development of ordered introgression lines (ILs), facilitating breeding for desirable traits. Despite the utility of these wild relatives and their associated ILs, few finished genome sequences have been produced to aid genetic and genomic studies. Here we report a chromosome-scale genome assembly for S. lycopersicoides LA2951, which contains 37 938 predicted protein-coding genes. With the aid of this genome assembly, we have precisely delimited the boundaries of the S. lycopersicoides introgressions in a set of S. lycopersicum cv. VF36 × LA2951 ILs. We demonstrate the usefulness of the LA2951 genome by identifying several quantitative trait loci for phenolics and carotenoids, including underlying candidate genes, and by investigating the genome organization and immunity-associated function of the clustered Pto gene family. In addition, syntenic analysis of R2R3MYB genes sheds light on the identity of the Aubergine locus underlying anthocyanin production. The genome sequence and IL map provide valuable resources for studying fruit nutrient/quality traits, pathogen resistance, and environmental stress tolerance. We present a new genome resource for the wild species S. lycopersicoides, which we use to shed light on the Aubergine locus responsible for anthocyanin production. We also provide IL boundary mappings, which facilitated identifying novel carotenoid quantitative trait loci of which one was likely driven by an uncharacterized lycopene ß-cyclase whose function we demonstrate.

Chromosome-level assembly of the Neolamarckia cadamba genome provides insights into the evolution of cadambine biosynthesis.

Neolamarckia cadamba (Roxb.), a close relative of Coffea canephora and Ophiorrhiza pumila, is an important traditional medicine in Southeast Asia. Three major glycosidic monoterpenoid indole alkaloids (MIAs), cadambine and its derivatives 3ß-isodihydrocadambine and 3ß-dihydrocadambine, accumulate in the bark and leaves, and exhibit antimalarial, antiproliferative, antioxidant, anticancer and anti-inflammatory activities. Here, we report a chromosome-scale N. cadamba genome, with 744.5 Mb assembled into 22 pseudochromosomes with contig N50 and scaffold N50 of 824.14 Kb and 29.20 Mb, respectively. Comparative genomic analysis of N. cadamba with Co. canephora revealed that N. cadamba underwent a relatively recent whole-genome duplication (WGD) event after diverging from Co. canephora, which contributed to the evolution of the MIA biosynthetic pathway. We determined the key intermediates of the cadambine biosynthetic pathway and further showed that NcSTR1 catalyzed the synthesis of strictosidine in N. cadamba. A new component, epoxystrictosidine (C27H34N2O10, m/z 547.2285), was identified in the cadambine biosynthetic pathway. Combining genome-wide association study (GWAS), population analysis, multi-omics analysis and metabolic gene cluster prediction, this study will shed light on the evolution of MIA biosynthetic pathway genes. This N. cadamba reference sequence will accelerate the understanding of the evolutionary history of specific metabolic pathways and facilitate the development of tools for enhancing bioactive productivity by metabolic engineering in microbes or by molecular breeding in plants.

Identification of starch candidate genes using SLAF-seq and BSA strategies and development of related SNP-CAPS markers in tetraploid potato.

Potato starch is an essential nutrient for humans and is widely used worldwide. Locating relevant genomic regions, mining stable genes and developing candidate gene markers can promote the breeding of new high-starch potato varieties. A total of 106 F1 individuals and their parents (YSP-4 × MIN-021) were used as test materials, from which 20 plants with high starch content and 20 with low starch content were selected to construct DNA pools for site-specific amplified fragment sequencing (SLAF-seq) and bulked segregation analysis (BSA). A genomic region related to the starch traits was first identified in the 0-5.62 Mb of chromosome 2 in tetraploid potato. In this section, a total of 41 non-synonymous genes, which were considered as candidate genes related to the starch trait, were annotated through a basic local alignment search tool (BLAST) search of multiple databases. Six candidate genes for starch (PGSC0003DMG400017793, PGSC0003DMG400035245, PGSC0003DMG400036713, PGSC0003DMG400040452, PGSC0003DMG400006636 and PGSC0003DMG400044547) were further explored. In addition, cleaved amplified polymorphic sequence (CAPS) markers were developed based on single nucleotide polymorphism (SNP) sites associated with the starch candidate genes. SNP-CAPS markers chr2-CAPS6 and chr2-CAPS21 were successfully developed and validated with the F2 population and 24 tetraploid potato varieties (lines). Functional analysis and cloning of the candidate genes associated with potato starch will be performed in further research, and the SNP-CAPS markers chr2-CAPS6 and chr2-CAPS21 can be further used in marker-assisted selection breeding of tetraploid potato varieties with high starch content.

QTL mapping in outbred tetraploid (and diploid) diallel populations.

Over the last decade, multiparental populations have become a mainstay of genetics research in diploid species. Our goal was to extend this paradigm to autotetraploids by developing software for quantitative trait locus (QTL) mapping in connected F1 populations derived from a set of shared parents. For QTL discovery, phenotypes are regressed on the dosage of parental haplotypes to estimate additive effects. Statistical properties of the model were explored by simulating half-diallel diploid and tetraploid populations with different population sizes and numbers of parents. Across scenarios, the number of progeny per parental haplotype (pph) largely determined the statistical power for QTL detection and accuracy of the estimated haplotype effects. Multiallelic QTL with heritability 0.2 were detected with 90% probability at 25 pph and genome-wide significance level 0.05, and the additive haplotype effects were estimated with over 90% accuracy. Following QTL discovery, the software enables a comparison of models with multiple QTL and nonadditive effects. To illustrate, we analyzed potato tuber shape in a half-diallel population with three tetraploid parents. A well-known QTL on chromosome 10 was detected, for which the inclusion of digenic dominance lowered the Deviance Information Criterion (DIC) by 17 points compared to the additive model. The final model also contained a minor QTL on chromosome 1, but higher-order dominance and epistatic effects were excluded based on the DIC. In terms of practical impacts, the software is already being used to select offspring based on the effect and dosage of particular haplotypes in breeding programs.

LD-CNV: rapid and simple discovery of chromosomal translocations using linkage disequilibrium between copy number variable loci.

Large-scale structural variations, such as chromosomal translocations, can have profound effects on fitness and phenotype, but are difficult to identify and characterize. Here, we describe a simple and effective method aimed at identifying translocations using only the dosage of sequence reads mapped on the reference genome. We binned reads on genomic segments sized according to sequencing coverage and identified instances when copy number segregated in populations. For each dosage-polymorphic 1 Mb bin, we tested independence, effectively an apparent linkage disequilibrium (LD), with other variable bins. In nine potato (Solanum tuberosum) dihaploid families translocations affecting pericentromeric regions were common and in two cases were due to genomic misassembly. In two populations, we found evidence for translocation affecting euchromatic arms. In cv. PI 310467, a nonreciprocal translocation between chromosomes (chr.) 7 and 8 resulted in a 5-3 copy number change affecting several Mb at the respective chromosome tips. In cv. "Alca Tarma," the terminal arm of chr. 4 translocated to the tip of chr. 1. Using oligonucleotide-based fluorescent in situ hybridization painting probes (oligo-FISH), we tested and confirmed the predicted arrangement in PI 310467. In 192 natural accessions of Arabidopsis thaliana, dosage haplotypes tended to vary continuously and resulted in higher noise, while apparent LD between pericentromeric regions suggested the effect of repeats. This method, LD-CNV, should be useful in species where translocations are suspected because it tests linkage without the need for genotyping.

In-silico genome wide analysis of Mitogen activated protein kinase kinase kinase gene family in C. sinensis.

Mitogen activated protein kinase kinase kinase (MAPKKK) form the upstream component of MAPK cascade. It is well characterized in several plants such as Arabidopsis and rice however the knowledge about MAPKKKs in tea plant is largely unknown. In the present study, MAPKKK genes of tea were obtained through a genome wide search using Arabidopsis thaliana as the reference genome. Among 59 candidate MAPKKK genes in tea, 17 genes were MEKK-like, 31 genes were Raf-like and 11 genes were ZIK- like. Additionally, phylogenetic relationships were established along with structural analysis, which includes gene structure, its location as well as conserved motifs, cis-acting regulatory elements and functional domain signatures that were systematically examined. Also, on the basis of one orthologous gene found between tea and Arabidopsis, functional interaction was carried out in C. sinensis based on an Arabidopsis association model. The expressional profiles indicated major involvement of MAPKKK genes from tea in response to various abiotic stress factors. Taken together, this study provides the targets for additional inclusive identification, functional study, and provides comprehensive knowledge for a better understanding of the MAPKKK cascade regulatory network in C. sinensis.

Linkage and QTL mapping for tuber shape and specific gravity in a tetraploid mapping population of potato representing the russet market class.

BACKGROUND: Tuber shape and specific gravity (dry matter) are important agronomic traits in potato processing and impact production costs, quality, and consistency of the final processed food products such as French fries and potato chips. In this study, linkage and QTL mapping were performed for these two traits to allow for the implementation of marker-assisted selection to facilitate breeding efforts in the russet market class. Two parents, Rio Grande Russet (female) and Premier Russet (male) and their 205 F1 progenies were initially phenotyped for tuber shape and specific gravity in field trials conducted in Idaho and North Carolina in 2010 and 2011, with specific gravity also being measured in Minnesota in 2011. Progenies and parents were previously genotyped using the Illumina SolCAP Infinium 8303 Potato SNP array, with ClusterCall and MAPpoly (R-packages) subsequently used for autotetraploid SNP calling and linkage mapping in this study. The 12 complete linkage groups and phenotypic data were then imported into QTLpoly, an R-package designed for polyploid QTL analyses. RESULTS: Significant QTL for tuber shape were detected on chromosomes 4, 7, and 10, with heritability estimates ranging from 0.09 to 0.36. Significant tuber shape QTL on chromosomes 4 and 7 were specific to Idaho and North Carolina environments, respectively, whereas the QTL on chromosome 10 was significant regardless of growing environment. Single marker analyses identified alleles in the parents associated with QTL on chromosomes 4, 7, and 10 that contributed to significant differences in tuber shape among progenies. Significant QTL were also identified for specific gravity on chromosomes 1 and 5 with heritability ranging from 0.12 to 0.21 and were reflected across environments. CONCLUSION: Fully automated linkage mapping and QTL analysis were conducted to identify significant QTL for tuber shape and dry matter in a tetraploid mapping population representing the russet market class. The findings are important for the development of molecular markers useful to potato breeders for marker-assisted selection for the long tuber shape and acceptable dry matter required by the potato industry within this important market class.

Genetic Dissection of Phosphorus Use Efficiency in a Maize Association Population under Two P Levels in the Field.

Phosphorus (P) deficiency is an important challenge the world faces while having to increase crop yields. It is therefore necessary to select maize (Zea may L.) genotypes with high phosphorus use efficiency (PUE). Here, we extensively analyzed the biomass, grain yield, and PUE-related traits of 359 maize inbred lines grown under both low-P and normal-P conditions. A significant decrease in grain yield per plant and biomass, an increase in PUE under low-P condition, as well as significant correlations between the two treatments were observed. In a genome-wide association study, 49, 53, and 48 candidate genes were identified for eleven traits under low-P, normal-P conditions, and in low-P tolerance index (phenotype under low-P divided by phenotype under normal-P condition) datasets, respectively. Several gene ontology pathways were enriched for the genes identified under low-P condition. In addition, seven key genes related to phosphate transporter or stress response were molecularly characterized. Further analyses uncovered the favorable haplotype for several core genes, which is less prevalent in modern lines but often enriched in a specific subpopulation. Collectively, our research provides progress in the genetic dissection and molecular characterization of PUE in maize.

Genome-wide identification, evolution, and transcript profiling of Aldehyde dehydrogenase superfamily in potato during development stages and stress conditions.

The Aldehyde dehydrogenase (ALDH) superfamily comprises a group of enzymes involved in the scavenging of toxic aldehyde molecules by converting them into their corresponding non-toxic carboxylic acids. A genome-wide study in potato identified a total of 22 ALDH genes grouped into ten families that are presented unevenly throughout all the 12 chromosomes. Based on the evolutionary analysis of ALDH proteins from different plant species, ALDH2 and ALDH3 were found to be the most abundant families in the plant, while ALDH18 was found to be the most distantly related one. Gene expression analysis revealed that the expression of StALDH genes is highly tissue-specific and divergent in various abiotic, biotic, and hormonal treatments. Structural modelling and functional analysis of selected StALDH members revealed conservancy in their secondary structures and cofactor binding sites. Taken together, our findings provide comprehensive information on the ALDH gene family in potato that will help in developing a framework for further functional studies.

Non-homologous chromosome pairing during meiosis in haploid Brassica rapa.

KEY MESSAGE: Cytological observations of chromosome pairing showed that evolutionarily genome duplication might reshape non-homologous pairing during meiosis in haploid B. rapa. A vast number of flowering plants have evolutionarily undergone whole genome duplication (WGD) event. Typically, Brassica rapa is currently considered as an evolutionary mesohexaploid, which has more complicated genomic constitution among flowering plants. In this study, we demonstrated chromosome behaviors in haploid B. rapa to understand how meiosis proceeds in presence of a single homolog. The findings showed that a diploid-like chromosome pairing was generally adapted during meiosis in haploid B. rapa. Non-homologous chromosomes in haploid cells paired at a high-frequency at metaphase I, over 50% of examined meiocytes showed at least three pairs of bivalents then equally segregated at anaphase I during meiosis. The fluorescence immunostaining showed that the cytoskeletal configurations were mostly well-organized during meiosis. Moreover, the expressed genes identified at meiosis in floral development was rather similar between haploid and diploid B. rapa, especially the expression of known hallmark genes pivotal to chromosome synapsis and homologous recombination were mostly in haploid B. rapa. Whole-genome duplication evolutionarily homology of genomic segments might be an important reason for this phenomenon, which would reshape the first division course of meiosis and influence pollen development in plants.