The genetic and epigenetic landscape of the Arabidopsis centromeres

Telomere-to-telomere carrot (Daucus carota) genome assembly reveals carotenoid characteristics

Carrot (Daucus carota) is an Apiaceae plant with multi-colored fleshy roots that provides a model system for carotenoid research. In this study, we assembled a 430.40 Mb high-quality gapless genome to the telomere-to-telomere (T2T) level of "Kurodagosun" carrot. In total, 36 268 genes were identified and 34 961 of them were functionally annotated. The proportion of repeat sequences in the genome was 55.3%, mainly long terminal repeats. Depending on the coverage of the repeats, 14 telomeres and 9 centromeric regions on the chromosomes were predicted. A phylogenetic analysis showed that carrots evolved early in the family Apiaceae. Based on the T2T genome, we reconstructed the carotenoid metabolic pathway and identified the structural genes that regulate carotenoid biosynthesis. Among the 65 genes that were screened, 9 were newly identified. Additionally, some gene sequences overlapped with transposons, suggesting replication and functional differentiation of carotenoid-related genes during carrot evolution. Given that some gene copies were barely expressed during development, they might be functionally redundant. Comparison of 24 cytochrome P450 genes associated with carotenoid biosynthesis revealed the tandem or proximal duplication resulting in expansion of CYP gene family. These results provided molecular information for carrot carotenoid accumulation and contributed to a new genetic resource.

The gap-free genome of mulberry elucidates the architecture and evolution of polycentric chromosomes

Mulberry is a fundamental component of the global sericulture industry, and its positive impact on our health and the environment cannot be overstated. However, the mulberry reference genomes reported previously remained unassembled or unplaced sequences. Here, we report the assembly and analysis of the telomere-to-telomere gap-free reference genome of the mulberry species, Morus notabilis, which has emerged as an important reference in mulberry gene function research and genetic improvement. The mulberry gap-free reference genome produced here provides an unprecedented opportunity for us to study the structure and function of centromeres. Our results revealed that all mulberry centromeric regions share conserved centromeric satellite repeats with different copies. Strikingly, we found that M. notabilis is a species with polycentric chromosomes and the only reported polycentric chromosome species up to now. We propose a compelling model that explains the formation mechanism of new centromeres and addresses the unsolved scientific question of the chromosome fusion-fission cycle in mulberry species. Our study sheds light on the functional genomics, chromosome evolution, and genetic improvement of mulberry species.

A complete telomere-to-telomere assembly of the maize genome

A complete telomere-to-telomere (T2T) finished genome has been the long pursuit of genomic research. Through generating deep coverage ultralong Oxford Nanopore Technology (ONT) and PacBio HiFi reads, we report here a complete genome assembly of maize with each chromosome entirely traversed in a single contig. The 2,178.6 Mb T2T Mo17 genome with a base accuracy of over 99.99% unveiled the structural features of all repetitive regions of the genome. There were several super-long simple-sequence-repeat arrays having consecutive thymine-adenine-guanine (TAG) tri-nucleotide repeats up to 235 kb. The assembly of the entire nucleolar organizer region of the 26.8 Mb array with 2,974 45S rDNA copies revealed the enormously complex patterns of rDNA duplications and transposon insertions. Additionally, complete assemblies of all ten centromeres enabled us to precisely dissect the repeat compositions of both CentC-rich and CentC-poor centromeres. The complete Mo17 genome represents a major step forward in understanding the complexity of the highly recalcitrant repetitive regions of higher plant genomes.

Holocentromeres can consist of merely a few megabase-sized satellite arrays

The centromere is the chromosome region where microtubules attach during cell division. In contrast to monocentric chromosomes with one centromere, holocentric species usually distribute hundreds of centromere units along the entire chromatid. We assembled the chromosome-scale reference genome and analyzed the holocentromere and (epi)genome organization of the lilioid Chionographis japonica. Remarkably, each of its holocentric chromatids consists of only 7 to 11 evenly spaced megabase-sized centromere-specific histone H3-positive units. These units contain satellite arrays of 23 and 28 bp-long monomers capable of forming palindromic structures. Like monocentric species, C. japonica forms clustered centromeres in chromocenters at interphase. In addition, the large-scale eu- and heterochromatin arrangement differs between C. japonica and other known holocentric species. Finally, using polymer simulations, we model the formation of prometaphase line-like holocentromeres from interphase centromere clusters. Our findings broaden the knowledge about centromere diversity, showing that holocentricity is not restricted to species with numerous and small centromere units.

Gap-free genome assembly of anadromous Coilia nasus

The Chinese tapertail anchovy, Coilia nasus, is a socioeconomically important anadromous fish that migrates from near ocean waters to freshwater to spawn every spring. The analysis of genomic architecture and information of C. nasus were hindered by the previously released versions of reference genomes with gaps. Here, we report the assembly of a chromosome-level gap-free genome of C. nasus by incorporating high-coverage and accurate long-read sequence data with multiple assembly strategies. All 24 chromosomes were assembled without gaps, representing the highest completeness and assembly quality. We assembled the genome with a size of 851.67 Mb and used BUSCO to estimate the completeness of the assembly as 92.5%. Using a combination of de novo prediction, protein homology and RNA-seq annotation, 21,900 genes were functionally annotated, representing 99.68% of the total predicted protein-coding genes. The availability of gap-free reference genomes for C. nasus will provide the opportunity for understanding genome structure and function, and will also lay a solid foundation for further management and conservation of this important species.

Scalable telomere-to-telomere assembly for diploid and polyploid genomes with double graph

Despite recent advances in the length and the accuracy of long-read data, building haplotype-resolved genome assemblies from telomere to telomere still requires considerable computational resources. In this study, we present an efficient de novo assembly algorithm that combines multiple sequencing technologies to scale up population-wide telomere-to-telomere assemblies. By utilizing twenty-two human and two plant genomes, we demonstrate that our algorithm is around an order of magnitude cheaper than existing methods, while producing better diploid and haploid assemblies. Notably, our algorithm is the only feasible solution to the haplotype-resolved assembly of polyploid genomes.

Recent Advances in Assembly of Complex Plant Genomes

Over the past 20 years, tremendous advances in sequencing technologies and computational algorithms have spurred plant genomic research into a thriving era with hundreds of genomes decoded already, ranging from those of nonvascular plants to those of flowering plants. However, complex plant genome assembly is still challenging and remains difficult to fully resolve with conventional sequencing and assembly methods due to high heterozygosity, highly repetitive sequences, or high ploidy characteristics of complex genomes. Herein, we summarize the challenges of and advances in complex plant genome assembly, including feasible experimental strategies, upgrades to sequencing technology, existing assembly methods, and different phasing algorithms. Moreover, we list actual cases of complex genome projects for readers to refer to and draw upon to solve future problems related to complex genomes. Finally, we expect that the accurate, gapless, telomere-to-telomere, and fully phased assembly of complex plant genomes could soon become routine.

Cycles of satellite and transposon evolution in Arabidopsis centromeres

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation^1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox^3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

A near-complete genome assembly of Brassica rapa provides new insights into the evolution of centromeres

Brassica rapa comprises many important cultivated vegetables and oil crops. However, Chiifu v3.0, the current B. rapa reference genome, still contains hundreds of gaps. Here, we presented a near-complete genome assembly of B. rapa Chiifu v4.0, which was 424.59 Mb with only two gaps, using Oxford Nanopore Technology (ONT) ultralong-read sequencing and Hi-C technologies. The new assembly contains 12 contigs, with a contig N50 of 38.26 Mb. Eight of the ten chromosomes were entirely reconstructed in a single contig from telomere to telomere. We found that the centromeres were mainly invaded by ALE and CRM long terminal repeats (LTRs). Moreover, there is a high divergence of centromere length and sequence among B. rapa genomes. We further found that centromeres are enriched for Copia invaded at 0.14 MYA on average, while pericentromeres are enriched for Gypsy LTRs invaded at 0.51 MYA on average. These results indicated the different invasion mechanisms of LTRs between the two structures. In addition, a novel repetitive sequence PCR630 was identified in the pericentromeres of B. rapa. Overall, the near-complete genome assembly, B. rapa Chiifu v4.0, offers valuable tools for genomic and genetic studies of Brassica species and provides new insights into the evolution of centromeres.

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1^In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish.

The complete reference genome for grapevine (Vitis vinifera L.) genetics and breeding

Grapevine is one of the most economically important crops worldwide. However, the previous versions of the grapevine reference genome tipically consist of thousands of fragments with missing centromeres and telomeres, limiting the accessibility of the repetitive sequences, the centromeric and telomeric regions, and the study of inheritance of important agronomic traits in these regions. Here, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar PN40024 using PacBio HiFi long reads. The T2T reference genome (PN_T2T) is 69 Mb longer with 9018 more genes identified than the 12X.v0 version. We annotated 67% repetitive sequences, 19 centromeres and 36 telomeres, and incorporated gene annotations of previous versions into the PN_T2T assembly. We detected a total of 377 gene clusters, which showed associations with complex traits, such as aroma and disease resistance. Even though PN40024 derives from nine generations of selfing, we still found nine genomic hotspots of heterozygous sites associated with biological processes, such as the oxidation-reduction process and protein phosphorylation. The fully annotated complete reference genome therefore constitutes an important resource for grapevine genetic studies and breeding programs.

A chromosome-scale genome assembly of Castanopsis hystrix provides new insights into the evolution and adaptation of Fagaceae species

Fagaceae species dominate forests and shrublands throughout the Northern Hemisphere, and have been used as models to investigate the processes and mechanisms of adaptation and speciation. Compared with the well-studied genus Quercus, genomic data is limited for the tropical-subtropical genus Castanopsis. Castanopsis hystrix is an ecologically and economically valuable species with a wide distribution in the evergreen broad-leaved forests of tropical-subtropical Asia. Here, we present a high-quality chromosome-scale reference genome of C. hystrix, obtained using a combination of Illumina and PacBio HiFi reads with Hi-C technology. The assembled genome size is 882.6 Mb with a contig N50 of 40.9 Mb and a BUSCO estimate of 99.5%, which are higher than those of recently published Fagaceae species. Genome annotation identified 37,750 protein-coding genes, of which 97.91% were functionally annotated. Repeat sequences constituted 50.95% of the genome and LTRs were the most abundant repetitive elements. Comparative genomic analysis revealed high genome synteny between C. hystrix and other Fagaceae species, despite the long divergence time between them. Considerable gene family expansion and contraction were detected in Castanopsis species. These expanded genes were involved in multiple important biological processes and molecular functions, which may have contributed to the adaptation of the genus to a tropical-subtropical climate. In summary, the genome assembly of C. hystrix provides important genomic resources for Fagaceae genomic research communities, and improves understanding of the adaptation and evolution of forest trees.

Structural variation among assembled genomes facilitates development of rapid and low-cost NOR-linked markers and NOR-telomere junction mapping in Arabidopsis

KEY MESSAGE

Genome-wide structural variants we identified and new NOR-linked markers we developed would be useful for future genome-wide association studies (GWAS), and for new gene/trait mapping purposes. Bioinformatic alignment of the assembled genomes of Col-0 and Sha ecotypes of Arabidopsis thaliana revealed ~ 13,000 genome-wide structural variants involving simple insertions or deletions and repeat contractions or expansions. Using some of these structural variants, we developed new, rapid, and low-cost PCR-based molecular markers that are genetically linked to the nucleolus organizer regions (NORs). A. thaliana has two NORs, one each on chromosome 2 (NOR2) and chromosome 4 (NOR4). Both NORs are ~ 4 Mb each, and hundreds of 45S ribosomal RNA (rRNA) genes are tandemly arrayed at these loci. Using previously characterized recombinant inbred lines (RILs) derived from Sha x Col-0  crosses, we validated the utility of the newly developed NOR-linked markers in genetically mapping rRNA genes and the associated telomeres to either NOR2 or NOR4. Lastly, we sequenced Sha genome using Oxford Nanopore Technology (ONT) and used the data to obtain sequences of NOR-telomere junctions, and with the help of RILs, we mapped them as new genetic markers to their respective NORs (NOR2-TEL2N and NOR4-TEL4N). The structural variants obtained from this study would serve as valuable data for genome-wide association studies (GWAS), and to rapidly design more genome-wide genetic (molecular) markers for new gene/trait mapping purposes.

De novo reconstruction of satellite repeat units from sequence data

Satellite DNA are long tandemly repeating sequences in a genome and may be organized as high-order repeats (HORs). They are enriched in centromeres and are challenging to assemble. Existing algorithms for identifying satellite repeats either require the complete assembly of satellites or only work for simple repeat structures without HORs. Here we describe Satellite Repeat Finder (SRF), a new algorithm for reconstructing satellite repeat units and HORs from accurate reads or assemblies without prior knowledge on repeat structures. Applying SRF to real sequence data, we showed that SRF could reconstruct known satellites in human and well-studied model organisms. We also found satellite repeats are pervasive in various other species, accounting for up to 12% of their genome contents but are often underrepresented in assemblies. With the rapid progress on genome sequencing, SRF will help the annotation of new genomes and the study of satellite DNA evolution even if such repeats are not fully assembled.

A complete gap-free diploid genome in Saccharum complex and the genomic footprints of evolution in the highly polyploid Saccharum genus

A diploid genome in the Saccharum complex facilitates our understanding of evolution in the highly polyploid Saccharum genus. Here we have generated a complete, gap-free genome assembly of Erianthus rufipilus, a diploid species within the Saccharum complex. The complete assembly revealed that centromere satellite homogenization was accompanied by the insertions of Gypsy retrotransposons, which drove centromere diversification. An overall low rate of gene transcription was observed in the palaeo-duplicated chromosome EruChr05 similar to other grasses, which might be regulated by methylation patterns mediated by homologous 24 nt small RNAs, and potentially mediating the functions of many nucleotide-binding site genes. Sequencing data for 211 accessions in the Saccharum complex indicated that Saccharum probably originated in the trans-Himalayan region from a diploid ancestor (x = 10) around 1.9-2.5 million years ago. Our study provides new insights into the origin and evolution of Saccharum and accelerates translational research in cereal genetics and genomics.

HiCAT: a tool for automatic annotation of centromere structure

Significant improvements in long-read sequencing technologies have unlocked complex genomic areas, such as centromeres, in the genome and introduced the centromere annotation problem. Currently, centromeres are annotated in a semi-manual way. Here, we propose HiCAT, a generalizable automatic centromere annotation tool, based on hierarchical tandem repeat mining to facilitate decoding of centromere architecture. We apply HiCAT to simulated datasets, human CHM13-T2T and gapless Arabidopsis thaliana genomes. Our results are generally consistent with previous inferences but also greatly improve annotation continuity and reveal additional fine structures, demonstrating HiCAT's performance and general applicability.

Combination of bacterial N-acyl homoserine lactones primes Arabidopsis defenses via jasmonate metabolism

N-acyl homoserine lactones (AHLs) are important players in plant-bacteria interactions. Different AHL-producing bacteria can improve plant growth and resistance against plant pathogens. In nature, plants may host a variety of AHL-producing bacteria and frequently experience numerous AHLs at the same time. Therefore, a coordinated response to combined AHL molecules is necessary. The purpose of this study was to explore the mechanism of AHL-priming using combined AHL molecules including N-(3-oxo-hexanoyl)-L-homoserine lactone, N-3-oxo-octanoyl-L-homoserine lactone, N-3-oxo-dodecanoyl-L-homoserine lactone, and N-3-oxo-tetradecanoyl-L-homoserine lactone and AHL-producing bacteria including Serratia plymuthica HRO-C48, Rhizobium etli CFN42, Burkholderia graminis DSM17151, and Ensifer meliloti (Sinorhizobium meliloti) Rm2011. We used transcriptome analysis, phytohormone measurements, as well as genetic and microbiological approaches to assess how the combination of structurally diverse AHL molecules influence Arabidopsis (Arabidopsis thaliana). Our findings revealed a particular response to a mixture of AHL molecules (AHL mix). Different expression patterns indicated that the reaction of plants exposed to AHL mix differs from that of plants exposed to single AHL molecules. In addition, different content of jasmonic acid (JA) and derivatives revealed that jasmonates play an important role in AHL mix-induced priming. The fast and stable decreased concentration of COOH-JA-Ile after challenge with the flagellin-derived peptide flg22 indicated that AHL mix modifies the metabolism of jasmonates. Study of various JA- and salicylic acid-related Arabidopsis mutants strengthened the notion that JA homeostasis is involved in AHL-priming. Understanding how the combination of AHLs primes plants for enhanced resistance has the potential to broaden our approaches in sustainable agriculture and will help to effectively protect plants against pathogens.

A gap-free genome assembly of Chlamydomonas reinhardtii and detection of translocations induced by CRISPR-mediated mutagenesis

Genomic assemblies of the unicellular green alga Chlamydomonas reinhardtii have provided important resources for researchers. However, assembly errors, large gaps, and unplaced scaffolds as well as strain-specific variants currently impede many types of analysis. By combining PacBio HiFi and Oxford Nanopore long-read technologies, we generated a de novo genome assembly for strain CC-5816, derived from crosses of strains CC-125 and CC-124. Multiple methods of evaluating genome completeness and base-pair error rate suggest that the final telomere-to-telomere assembly is highly accurate. The CC-5816 assembly enabled previously difficult analyses that include characterization of the 17 centromeres, rDNA arrays on three chromosomes, and 56 insertions of organellar DNA into the nuclear genome. Using Nanopore sequencing, we identified sites of cytosine (CpG) methylation, which are enriched at centromeres. We analyzed CRISPR-Cas9 insertional mutants in the PF23 gene. Two of the three alleles produced progeny that displayed patterns of meiotic inviability that suggested the presence of a chromosomal aberration. Mapping Nanopore reads from pf23-2 and pf23-3 onto the CC-5816 genome showed that these two strains each carry a translocation that was initiated at the PF23 gene locus on chromosome 11 and joined with chromosomes 5 or 3, respectively. The translocations were verified by demonstrating linkage between loci on the two translocated chromosomes in meiotic progeny. The three pf23 alleles display the expected short-cilia phenotype, and immunoblotting showed that pf23-2 lacks the PF23 protein. Our CC-5816 genome assembly will undoubtedly provide an important tool for the Chlamydomonas research community.

Technology-driven approaches for meiosis research in tomato and wild relatives

Meiosis is a specialized cell division during reproduction where one round of chromosomal replication is followed by genetic recombination and two rounds of segregation to generate recombined, ploidy-reduced spores. Meiosis is crucial to the generation of new allelic combinations in natural populations and artificial breeding programs. Several plant species are used in meiosis research including the cultivated tomato (Solanum lycopersicum) which is a globally important crop species. Here we outline the unique combination of attributes that make tomato a powerful model system for meiosis research. These include the well-characterized behavior of chromosomes during tomato meiosis, readily available genomics resources, capacity for genome editing, clonal propagation techniques, lack of recent polyploidy and the possibility to generate hybrids with twelve related wild species. We propose that further exploitation of genome bioinformatics, genome editing and artificial intelligence in tomato will help advance the field of plant meiosis research. Ultimately this will help address emerging themes including the evolution of meiosis, how recombination landscapes are determined, and the effect of temperature on meiosis.

Satellite DNA evolution in Corvoidea inferred from short and long reads

Satellite DNA (satDNA) is a fast-evolving portion of eukaryotic genomes. The homogeneous and repetitive nature of such satDNA causes problems during the assembly of genomes, and therefore it is still difficult to study it in detail in nonmodel organisms as well as across broad evolutionary timescales. Here, we combined the use of short- and long-read data to explore the diversity and evolution of satDNA between individuals of the same species and between genera of birds spanning ~40 millions of years of bird evolution using birds-of-paradise (Paradisaeidae) and crow (Corvus) species. These avian species highlighted the presence of a GC-rich Corvoidea satellitome composed of 61 satellite families and provided a set of candidate satDNA monomers for being centromeric on the basis of length, abundance, homogeneity and transcription. Surprisingly, we found that the satDNA of crow species rapidly diverged between closely related species while the satDNA appeared more similar between birds-of-paradise species belonging to different genera.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Gap-free genome assembly and comparative analysis reveal the evolution and anthocyanin accumulation mechanism of Rhodomyrtus tomentosa

Rhodomyrtus tomentosa is an important fleshy-fruited tree and a well-known medicinal plant of the Myrtaceae family that is widely cultivated in tropical and subtropical areas of the world. However, studies on the evolution and genomic breeding of R. tomentosa were hindered by the lack of a reference genome. Here, we presented a chromosome-level gap-free T2T genome assembly of R. tomentosa using PacBio and ONT long read sequencing. We assembled the genome with size of 470.35 Mb and contig N50 of ~43.80 Mb with 11 pseudochromosomes. A total of 33 382 genes and 239.31 Mb of repetitive sequences were annotated in this genome. Phylogenetic analysis elucidated the independent evolution of R. tomentosa starting from 14.37MYA and shared a recent WGD event with other Myrtaceae species. We identified four major compounds of anthocyanins and their synthetic pathways in R. tomentosa. Comparative genomic and gene expression analysis suggested the coloring and high anthocyanin accumulation in R. tomentosa tends to be determined by the activation of anthocyanin synthesis pathway. The positive selection and up-regulation of MYB transcription factors were the implicit factors in this process. The copy number increase of downstream anthocyanin transport-related OMT and GST gene were also detected in R. tomentosa. Expression analysis and pathway identification enriched the importance of starch degradation, response to stimuli, effect of hormones, and cell wall metabolism during the fleshy fruit development in Myrtaceae. Our genome assembly provided a foundation for investigating the origins and differentiation of Myrtaceae species and accelerated the genetic improvement of R. tomentosa.

A proposed metric set for evaluation of genome assembly quality

Quality control is essential for genome assemblies; however, a consensus has yet to be reached on what metrics should be adopted for the evaluation of assembly quality. N50 is widely used for contiguity measurement, but its effectiveness is constantly in question. Prevailing metrics for the completeness evaluation focus on gene space, yet challenging areas such as tandem repeats are commonly overlooked. Achieving correctness has become an indispensable dimension for quality control, while prevailing assembly releases lack scores reflecting this aspect. We propose a metric set with a set of statistic indexes for effective, comprehensive evaluation of assemblies and provide a score of a finished assembly for each metric, which can be utilized as a benchmark for achieving high-quality genome assemblies.

The roles of recombination and selection in shaping genomic divergence in an incipient ecological species complex

Speciation genomic studies have revealed that genomes of diverging lineages are shaped jointly by the actions of gene flow and selection. These evolutionary forces acting in concert with processes such as recombination and genome features such as gene density shape a mosaic landscape of divergence. We investigated the roles of recombination and gene density in shaping the patterns of differentiation and divergence between the cyclically parthenogenetic ecological sister-taxa, Daphnia pulicaria and Daphnia pulex. First, we assembled a phased chromosome-scale genome assembly using trio-binning for D. pulicaria and constructed a genetic map using an F2-intercross panel to understand sex-specific recombination rate heterogeneity. Finally, we used a ddRADseq data set with broad geographic sampling of D. pulicaria, D. pulex, and their hybrids to understand the patterns of genome-scale divergence and demographic parameters. Our study provides the first sex-specific estimates of recombination rates for a cyclical parthenogen, and unlike other eukaryotic species, we observed male-biased heterochiasmy in D. pulicaria, which may be related to this somewhat unique breeding mode. Additionally, regions of high gene density and recombination are generally more divergent than regions of suppressed recombination. Outlier analysis indicated that divergent genomic regions are probably driven by selection on D. pulicaria, the derived lineage colonizing a novel lake habitat. Together, our study supports a scenario of selection acting on genes related to local adaptation shaping genome-wide patterns of differentiation despite high local recombination rates in this species complex. Finally, we discuss the limitations of our data in light of demographic uncertainty.

Telomere-to-telomere and haplotype-resolved genome of the kiwifruit Actinidia eriantha

Actinidia eriantha is a characteristic fruit tree featuring with great potential for its abundant vitamin C and strong disease resistance. It has been used in a wide range of breeding programs and functional genomics studies. Previously published genome assemblies of A. eriantha are quite fragmented and not highly contiguous. Using multiple sequencing strategies, we get the haplotype-resolved and gap-free genomes of an elite breeding line "Midao 31" (MD), termed MDHAPA and MDHAPB. The new assemblies anchored to 29 pseudochromosome pairs with a length of 619.3 Mb and 611.7 Mb, as well as resolved 27 and 28 gap-close chromosomes in a telomere-to-telomere (T2T) manner. Based on the haplotype-resolved genome, we found that most alleles experienced purifying selection and coordinately expressed. Owing to the high continuity of assemblies, we defined the centromeric regions of A. eriantha, and identified the major repeating monomer, which is designated as Ae-CEN153. This resource lays a solid foundation for further functional genomics study and horticultural traits improvement in kiwifruit.

Two haplotype-resolved, gap-free genome assemblies for Actinidia latifolia and Actinidia chinensis shed light on the regulatory mechanisms of vitamin C and sucrose metabolism in kiwifruit

Kiwifruit is a recently domesticated horticultural fruit crop with substantial economic and nutritional value, especially because of the high content of vitamin C in its fruit. In this study, we de novo assembled two telomere-to-telomere kiwifruit genomes from Actinidia chinensis var. 'Donghong' (DH) and Actinidia latifolia 'Kuoye' (KY), with total lengths of 608 327 852 and 640 561 626 bp for 29 chromosomes, respectively. With a burst of structural variants involving inversion, translocations, and duplications within 8.39 million years, the metabolite content of DH and KY exhibited differences in saccharides, lignans, and vitamins. A regulatory ERF098 transcription factor family has expanded in KY and Actinidia eriantha, both of which have ultra-high vitamin C content. With each assembly phased into two complete haplotypes, we identified allelic variations between two sets of haplotypes, leading to protein sequence variations in 26 494 and 27 773 gene loci and allele-specific expression of 4687 and 12 238 homozygous gene pairs. Synchronized metabolome and transcriptome changes during DH fruit development revealed the same dynamic patterns in expression levels and metabolite contents; free fatty acids and flavonols accumulated in the early stages, but sugar substances and amino acids accumulated in the late stages. The AcSWEET9b gene that exhibits allelic dominance was further identified to positively correlate with high sucrose content in fruit. Compared with wild varieties and other Actinidia species, AcSWEET9b promoters were selected in red-flesh kiwifruits that have increased fruit sucrose content, providing a possible explanation on why red-flesh kiwifruits are sweeter. Collectively, these two gap-free kiwifruit genomes provide a valuable genetic resource for investigating domestication mechanisms and genome-based breeding of kiwifruit.

Centromere repositioning and shifts in wheat evolution

The centromere is the region of a chromosome that directs its separation and plays an important role in cell division and reproduction of organisms. Elucidating the dynamics of centromeres is an alternative strategy for exploring the evolution of wheat. Here, we comprehensively analyzed centromeres from the de novo-assembled common wheat cultivar Aikang58 (AK58), Chinese Spring (CS), and all sequenced diploid and tetraploid ancestors by chromatin immunoprecipitation sequencing, whole-genome bisulfite sequencing, RNA sequencing, assay for transposase-accessible chromatin using sequencing, and comparative genomics. We found that centromere-associated sequences were concentrated during tetraploidization and hexaploidization. Centromeric repeats of wheat (CRWs) have undergone expansion during wheat evolution, with strong interweaving between the A and B subgenomes post tetraploidization. We found that CENH3 prefers to bind with younger CRWs, as directly supported by immunocolocalization on two chromosomes (1A and 2A) of wild emmer wheat with dicentromeric regions, only one of which bound with CENH3. In a comparison of AK58 with CS, obvious centromere repositioning was detected on chromosomes 1B, 3D, and 4D. The active centromeres showed a unique combination of lower CG but higher CHH and CHG methylation levels. We also found that centromeric chromatin was more open than pericentromeric chromatin, with higher levels of gene expression but lower gene density. Frequent introgression between tetraploid and hexaploid wheat also had a strong influence on centromere position on the same chromosome. This study also showed that active wheat centromeres were genetically and epigenetically determined.

Telomere-to-telomere and gap-free reference genome assembly of the kiwifruit Actinidia chinensis

Kiwifruit is an economically and nutritionally important fruit crop with extremely high contents of vitamin C. However, the previously released versions of kiwifruit genomes all have a mass of unanchored or missing regions. Here, we report a highly continuous and completely gap-free reference genome of Actinidia chinensis cv. 'Hongyang', named Hongyang v4.0, which is the first to achieve two de novo haploid-resolved haplotypes, HY4P and HY4A. HY4P and HY4A have a total length of 606.1 and 599.6 Mb, respectively, with almost the entire telomeres and centromeres assembled in each haplotype. In comparison with Hongyang v3.0, the integrity and contiguity of Hongyang v4.0 is markedly improved by filling all unclosed gaps and correcting some misoriented regions, resulting in ~38.6-39.5 Mb extra sequences, which might affect 4263 and 4244 protein-coding genes in HY4P and HY4A, respectively. Furthermore, our gap-free genome assembly provides the first clue for inspecting the structure and function of centromeres. Globally, centromeric regions are characterized by higher-order repeats that mainly consist of a 153-bp conserved centromere-specific monomer (Ach-CEN153) with different copy numbers among chromosomes. Functional enrichment analysis of the genes located within centromeric regions demonstrates that chromosome centromeres may not only play physical roles for linking a pair of sister chromatids, but also have genetic features for participation in the regulation of cell division. The availability of the telomere-to-telomere and gap-free Hongyang v4.0 reference genome lays a solid foundation not only for illustrating genome structure and functional genomics studies but also for facilitating kiwifruit breeding and improvement.

Transgenerational epigenetic control of constitutive heterochromatin, transposons, and centromeres

Epigenetic mechanisms are important not only for development but also for genome stability and chromosome dynamics. The latter types of epigenetic controls can often be transgenerational. Here, we review recent progress in two examples of transgenerational epigenetic control: i) the control of constitutive heterochromatin and transposable elements and ii) epigenetic mechanisms that regulate centromere specification and functions. We also discuss the biological significance of enigmatic associations among centromeres, transposons, and constitutive heterochromatin.

Going broad and deep: sequencing-driven insights into plant physiology, evolution, and crop domestication

Deep sequencing is a term that has become embedded in the plant genomic literature in recent years and with good reason. A torrent of (largely) high-quality genomic and transcriptomic data has been collected and most of this has been publicly released. Indeed, almost 1000 plant genomes have been reported (www.plabipd.de) and the 2000 Plant Transcriptomes Project has long been completed. The EarthBioGenome project will dwarf even these milestones. That said, massive progress in understanding plant physiology, evolution, and crop domestication has been made by sequencing broadly (across a species) as well as deeply (within a single individual). We will outline the current state of the art in genome and transcriptome sequencing before we briefly review the most visible of these broad approaches, namely genome-wide association and transcriptome-wide association studies, as well as the compilation of pangenomes. This will include both (i) the most commonly used methods reliant on single nucleotide polymorphisms and short InDels and (ii) more recent examples which consider structural variants. We will subsequently present case studies exemplifying how their application has brought insight into either plant physiology or evolution and crop domestication. Finally, we will provide conclusions and an outlook as to the perspective for the extension of such approaches to different species, tissues, and biological processes.

Assembly of the 81.6 Mb centromere of pea chromosome 6 elucidates the structure and evolution of metapolycentric chromosomes

Centromeres in the legume genera Pisum and Lathyrus exhibit unique morphological characteristics, including extended primary constrictions and multiple separate domains of centromeric chromatin. These so-called metapolycentromeres resemble an intermediate form between monocentric and holocentric types, and therefore provide a great opportunity for studying the transitions between different types of centromere organizations. However, because of the exceedingly large and highly repetitive nature of metapolycentromeres, highly contiguous assemblies needed for these studies are lacking. Here, we report on the assembly and analysis of a 177.6 Mb region of pea (Pisum sativum) chromosome 6, including the 81.6 Mb centromere region (CEN6) and adjacent chromosome arms. Genes, DNA methylation profiles, and most of the repeats were uniformly distributed within the centromere, and their densities in CEN6 and chromosome arms were similar. The exception was an accumulation of satellite DNA in CEN6, where it formed multiple arrays up to 2 Mb in length. Centromeric chromatin, characterized by the presence of the CENH3 protein, was predominantly associated with arrays of three different satellite repeats; however, five other satellites present in CEN6 lacked CENH3. The presence of CENH3 chromatin was found to determine the spatial distribution of the respective satellites during the cell cycle. Finally, oligo-FISH painting experiments, performed using probes specifically designed to label the genomic regions corresponding to CEN6 in Pisum, Lathyrus, and Vicia species, revealed that metapolycentromeres evolved via the expansion of centromeric chromatin into neighboring chromosomal regions and the accumulation of novel satellite repeats. However, in some of these species, centromere evolution also involved chromosomal translocations and centromere repositioning.

Molecular cytogenetic map visualizes the heterozygotic genome and identifies translocation chromosomes in Citrus sinensis

Citrus sinensis is the most cultivated and economically valuable Citrus species in the world, whose genome has been assembled by three generation sequencings. However, chromosome recognition remains a problem due to the small size of chromosomes, and difficulty in differentiating between pseudo and real chromosomes because of a highly heterozygous genome. Here, we employ fluorescence in situ hybridization (FISH) with 9 chromosome painting probes, 30 oligo pools, and 8 repetitive sequences to visualize 18 chromosomes. Then, we develop an approach to identify each chromosome in one cell through single experiment of oligo-FISH and Chromoycin A3 (CMA) staining. By this approach, we construct a high-resolution molecular cytogenetic map containing the physical positions of CMA banding and 38 sequences of FISH including centromere regions, which enable us to visualize significant differences between homologous chromosomes. Based on the map, we locate several highly repetitive sequences on chromosomes and estimate sizes and copy numbers of each site. In particular, we discover the translocation regions of chromosomes 4 and 9 in C. sinensis "Valencia." The high-resolution molecular cytogenetic map will help improve understanding of sweet orange genome assembly and also provide a fundamental reference for investigating chromosome evolution and chromosome engineering for genetic improvement in Citrus.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

Graph-based pan-genomes: increased opportunities in plant genomics

Due to the development of sequencing technology and the great reduction in sequencing costs, an increasing number of plant genomes have been assembled, and numerous genomes have revealed large amounts of variations. However, a single reference genome does not allow the exploration of species diversity, and therefore the concept of pan-genome was developed. A pan-genome is a collection of all sequences available for a species, including a large number of consensus sequences, large structural variations, and small variations including single nucleotide polymorphisms and insertions/deletions. A simple linear pan-genome does not allow these structural variations to be intuitively characterized, so graph-based pan-genomes have been developed. These pan-genomes store sequence and structural variation information in the form of nodes and paths to store and display species variation information in a more intuitive manner. The key role of graph-based pan-genomes is to expand the coordinate system of the linear reference genome to accommodate more regions of genetic diversity. Here, we review the origin and development of graph-based pan-genomes, explore their application in plant research, and further highlight the application of graph-based pan-genomes for future plant breeding.

Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.

Approaching complete genomes, transcriptomes and epi-omes with accurate long-read sequencing

The year 2022 will be remembered as the turning point for accurate long-read sequencing, which now establishes the gold standard for speed and accuracy at competitive costs. We discuss the key bioinformatics techniques needed to power long reads across application areas and close with our vision for long-read sequencing over the coming years.

Telomere-to-telomere genome assembly of bitter melon (Momordica charantia L. var. abbreviata Ser.) reveals fruit development, composition and ripening genetic characteristics

Momordica charantia L. var. abbreviata Ser. (Mca), known as bitter gourd or bitter melon, is a Momordica variety with medicinal value and belongs to the Cucurbitaceae family. In view of the lack of genomic information on bitter gourd and other Momordica species and to promote Mca genomic research, we assembled a 295.6-Mb telomere-to-telomere (T2T) high-quality Mca genome with six gap-free chromosomes after Hi-C correction. This genome is anchored to 11 chromosomes, which is consistent with the karyotype information, and comprises 98 contigs (N50 of 25.4 Mb) and 95 scaffolds (N50 of 25.4 Mb). The Mca genome harbors 19 895 protein-coding genes, of which 45.59% constitute predicted repeat sequences. Synteny analysis revealed variations involved in fruit quality during the divergence of bitter gourd. In addition, assay for transposase-accessible chromatin by high-throughput sequencing and metabolic analysis showed that momordicosides and other substances are characteristic of Mca fruit pulp. A combined transcriptomic and metabolomic analysis revealed the mechanisms of pigment accumulation and cucurbitacin biosynthesis in Mca fruit peels, providing fundamental molecular information for further research on Mca fruit ripening. This report provides a new genetic resource for Momordica genomic studies and contributes additional insights into Cucurbitaceae phylogeny.

Chromosome Painting Using Chromosome-Specific BAC Clones

Chromosome painting (CP) refers to visualization of large chromosome regions, chromosome arms or entire chromosomes via fluorescence in situ hybridization (FISH) of chromosome-specific DNA sequences. For CP in crucifers (Brassicaceae), typically contigs of chromosome-specific bacterial artificial chromosomes (BAC) from Arabidopsis thaliana are applied as painting probes on chromosomes of A. thaliana or other species (comparative chromosome painting, CCP). CP/CCP enables to identify and trace particular chromosome regions and/or chromosomes throughout all mitotic and meiotic stages as well as corresponding interphase chromosome territories. However, extended pachytene chromosomes provide the highest resolution of CP/CCP. Fine-scale chromosome structure, structural chromosome rearrangements (such as inversions, translocations, centromere repositioning), and chromosome breakpoints can be investigated by CP/CCP. BAC DNA probes can be accompanied by other types of DNA probes, such as repetitive DNA, genomic DNA, or synthetic oligonucleotide probes. Here, we describe a robust step-by-step protocol of CP and CCP which proved to be efficient across the family Brassicaceae, but which is also applicable to other angiosperm families.

Automated assembly scaffolding using RagTag elevates a new tomato system for high-throughput genome editing

Advancing crop genomics requires efficient genetic systems enabled by high-quality personalized genome assemblies. Here, we introduce RagTag, a toolset for automating assembly scaffolding and patching, and we establish chromosome-scale reference genomes for the widely used tomato genotype M82 along with Sweet-100, a new rapid-cycling genotype that we developed to accelerate functional genomics and genome editing in tomato. This work outlines strategies to rapidly expand genetic systems and genomic resources in other plant species.

Centromere drive: chromatin conflict in meiosis

Centromeres are essential loci in eukaryotes that are necessary for the faithful segregation of chromosomes in mitosis and meiosis. Centromeres organize the kinetochore, the protein machine that attaches sister chromatids or homologous chromosomes to spindle microtubules and regulates their disjunction. Centromeres have both genetic and epigenetic determinants, which can come into conflict in asymmetric female meiosis in seed plants and animals. The centromere drive model was proposed to describe this conflict and explain how it leads to the rapid evolution of both centromeres and kinetochores. Recent studies confirm key aspects of the centromere drive model, clarify its mechanisms, and implicate rapid centromere/kinetochore evolution in hybrid inviability between species.

Pushing the limits of HiFi assemblies reveals centromere diversity between two Arabidopsis thaliana genomes

Although long-read sequencing can often enable chromosome-level reconstruction of genomes, it is still unclear how one can routinely obtain gapless assemblies. In the model plant Arabidopsis thaliana, other than the reference accession Col-0, all other accessions de novo assembled with long-reads until now have used PacBio continuous long reads (CLR). Although these assemblies sometimes achieved chromosome-arm level contigs, they inevitably broke near the centromeres, excluding megabases of DNA from analysis in pan-genome projects. Since PacBio high-fidelity (HiFi) reads circumvent the high error rate of CLR technologies, albeit at the expense of read length, we compared a CLR assembly of accession Eyach15-2 to HiFi assemblies of the same sample. The use of five different assemblers starting from subsampled data allowed us to evaluate the impact of coverage and read length. We found that centromeres and rDNA clusters are responsible for 71% of contig breaks in the CLR scaffolds, while relatively short stretches of GA/TC repeats are at the core of >85% of the unfilled gaps in our best HiFi assemblies. Since the HiFi technology consistently enabled us to reconstruct gapless centromeres and 5S rDNA clusters, we demonstrate the value of the approach by comparing these previously inaccessible regions of the genome between the Eyach15-2 accession and the reference accession Col-0.

Oat chromosome and genome evolution defined by widespread terminal intergenomic translocations in polyploids

Structural chromosome rearrangements involving translocations, fusions and fissions lead to evolutionary variation between species and potentially reproductive isolation and variation in gene expression. While the wheats (Triticeae, Poaceae) and oats (Aveneae) all maintain a basic chromosome number of x=7, genomes of oats show frequent intergenomic translocations, in contrast to wheats where these translocations are relatively rare. We aimed to show genome structural diversity and genome relationships in tetraploid, hexaploid and octoploid Avena species and amphiploids, establishing patterns of intergenomic translocations across different oat taxa using fluorescence in situ hybridization (FISH) with four well-characterized repetitive DNA sequences: pAs120, AF226603, Ast-R171 and Ast-T116. In A. agadiriana (2n=4x=28), the selected probes hybridized to all chromosomes indicating that this species originated from one (autotetraploid) or closely related ancestors with the same genomes. Hexaploid amphiploids were confirmed as having the genomic composition AACCDD, while octoploid amphiploids showed three different genome compositions: AACCCCDD, AAAACCDD or AABBCCDD. The A, B, C, and D genomes of oats differ significantly in their involvement in non-centromeric, intercalary translocations. There was a predominance of distal intergenomic translocations from the C- into the D-genome chromosomes. Translocations from A- to C-, or D- to C-genome chromosomes were less frequent, proving that at least some of the translocations in oat polyploids are non-reciprocal. Rare translocations from A- to D-, D- to A- and C- to B-genome chromosomes were also visualized. The fundamental research has implications for exploiting genomic biodiversity in oat breeding through introgression from wild species potentially with contrasting chromosomal structures and hence deleterious segmental duplications or large deletions in amphiploid parental lines.

Frequent spontaneous structural rearrangements promote rapid genome diversification in a Brassica napus F1 generation

In a cross between two homozygous Brassica napus plants of synthetic and natural origin, we demonstrate that novel structural genome variants from the synthetic parent cause immediate genome diversification among F1 offspring. Long read sequencing in twelve F1 sister plants revealed five large-scale structural rearrangements where both parents carried different homozygous alleles but the heterozygous F1 genomes were not identical heterozygotes as expected. Such spontaneous rearrangements were part of homoeologous exchanges or segmental deletions and were identified in different, individual F1 plants. The variants caused deletions, gene copy-number variations, diverging methylation patterns and other structural changes in large numbers of genes and may have been causal for unexpected phenotypic variation between individual F1 sister plants, for example strong divergence of plant height and leaf area. This example supports the hypothesis that spontaneous de novo structural rearrangements after de novo polyploidization can rapidly overcome intense allopolyploidization bottlenecks to re-expand crops genetic diversity for ecogeographical expansion and human selection. The findings imply that natural genome restructuring in allopolyploid plants from interspecific hybridization, a common approach in plant breeding, can have a considerably more drastic impact on genetic diversity in agricultural ecosystems than extremely precise, biotechnological genome modifications.

Centromeres: From chromosome biology to biotechnology applications and synthetic genomes in plants

Centromeres are the genomic regions that organize and regulate chromosome behaviours during cell cycle, and their variations are associated with genome instability, karyotype evolution and speciation in eukaryotes. The highly repetitive and epigenetic nature of centromeres were documented during the past half century. With the aid of rapid expansion in genomic biotechnology tools, the complete sequence and structural organization of several plant and human centromeres were revealed recently. Here, we systematically summarize the current knowledge of centromere biology with regard to the DNA compositions and the histone H3 variant (CENH3)-dependent centromere establishment and identity. We discuss the roles of centromere to ensure cell division and to maintain the three-dimensional (3D) genomic architecture in different species. We further highlight the potential applications of manipulating centromeres to generate haploids or to induce polyploids offspring in plant for breeding programs, and of targeting centromeres with CRISPR/Cas for chromosome engineering and speciation. Finally, we also assess the challenges and strategies for de novo design and synthesis of centromeres in plant artificial chromosomes. The biotechnology applications of plant centromeres will be of great potential for the genetic improvement of crops and precise synthetic breeding in the future.

Long-read sequencing characterizes mitochondrial and plastid genome variants in Arabidopsis msh1 mutants

The abundant repeats in plant mitochondrial genomes can cause rapid genome rearrangements and are also a major obstacle in short-read sequencing studies. Nuclear-encoded proteins such as MSH1 are known to suppress the generation of repeat-associated mitochondrial genome variants, but our understanding of these mechanisms has been constrained by the limitations of short-read technologies. Here, we used highly accurate long-read sequencing (PacBio HiFi) to characterize mitochondrial and plastid genome variants in Arabidopsis thaliana msh1 mutant individuals. The HiFi reads provided a global view of recombination dynamics with detailed quantification of parental and crossover recombination products for both large and small repeats. We found that recombination breakpoints were distributed relatively evenly across the length of repeated sequences and detected widespread internal exchanges of sequence variants between pairs of imperfect repeats in the mitochondrial genome of msh1 mutants. Long-read assemblies of mitochondrial genomes from seven other A. thaliana wild-type accessions differed by repeat-mediated structural rearrangements similar to those observed in msh1 mutants, but they were all in a simple low-heteroplasmy state. The Arabidopsis plastid genome generally lacks small repeats and exhibited a very different pattern of variant accumulation in msh1 mutants compared with the mitochondrial genome. Our data illustrate the power of HiFi technology in studying repeat-mediated recombination in plant organellar genomes and improved the sequence resolution for recombinational processes suppressed by MSH1. Plant organellar genomes can undergo rapid rearrangements. Long-read sequencing provides a detailed and quantitative view of mitochondrial and plastid genome variants normally suppressed by MSH1, advancing our understanding of plant organellar genome dynamics.

Haplotype-resolved Chinese male genome assembly based on high-fidelity sequencing

The advantages of both the length and accuracy of high-fidelity (HiFi) reads enable chromosome-scale haplotype-resolved genome assembly. In this study, we sequenced a cell line named HJ, established from a Chinese Han male individual by using HiFi and Hi-C. We assembled two high-quality haplotypes of the HJ genome (haplotype 1 (H1): 3.1 Gb, haplotype 2 (H2): 2.9 Gb). The continuity (H1: contig N50 = 28.2 Mb, H2: contig N50 = 25.9 Mb) and completeness (BUSCO: H1 = 94.9%, H2 = 93.5%) are substantially better than those of other Chinese genomes, for example, HX1, NH1.0, and YH2.0. By comparing HJ genome with GRCh38, we reported the mutation landscape of HJ and found that 176 and 213 N-gaps were filled in H1 and H2, respectively. In addition, we detected 12.9 Mb and 13.4 Mb novel sequences containing 246 and 135 protein-coding genes in H1 and H2, respectively. Our results demonstrate the advantages of HiFi reads in haplotype-resolved genome assembly and provide two high-quality haplotypes of a potential Chinese genome as a reference for the Chinese Han population.

Using the Arabidopsis Information Resource (TAIR) to Find Information About Arabidopsis Genes

The Arabidopsis Information Resource (TAIR; http://arabidopsis.org) is a comprehensive web resource of Arabidopsis biology for plant scientists. TAIR curates and integrates information about genes, proteins, gene function, orthologs, gene expression, mutant phenotypes, biological materials such as clones and seed stocks, genetic markers, genetic and physical maps, genome organization, images of mutant plants, protein sub-cellular localizations, publications, and the research community. The various data types are extensively interconnected and can be accessed through a variety of web-based search and display tools. This article primarily focuses on some basic methods for searching, browsing, visualizing, and analyzing information about Arabidopsis genes and genomes. Additionally, we describe how members of the community can share data via JBrowse and the Generic Online Annotation Submission Tool (GOAT) in order to make their published research more accessible and visible. © 2022 Wiley Periodicals LLC. Basic Protocol 1: TAIR homepage, sitemap, and navigation Basic Protocol 2: Finding comprehensive information about Arabidopsis genes Basic Protocol 3: Using the Arabidopsis genome browser: JBrowse Basic Protocol 4: Using the Gene Ontology annotations for gene discovery and gene function analysis Basic Protocol 5: Using gene lists to download bulk datasets Basic Protocol 6: Using TAIR's analysis tools to find short sequences and motifs Basic Protocol 7: Using the TAIR generic online annotation tool (GOAT) to submit functional annotations for Arabidopsis (or any other species) genes Basic Protocol 8: Using PhyloGenes to visualize gene families and predict functions Basic Protocol 9: Using TAIR to browse Arabidopsis literature Basic Protocol 10: Using the synteny viewer to find and display syntenic regions from Arabidopsis and other plant species.

Epigenetic regulation in the shoot apical meristem

Epigenetic mechanisms form the basis of cellular memory, developmental decisions, and the cellular immune system that defends against transposons and viruses. Organs develop from the shoot apical meristem (SAM) to shape the plant's areal phenotype, and stem cells in the SAM serve as a functional germline. While many details on the regulation of stem cell pool size, organ initiation, and patterning at the meristem periphery are known, we know surprisingly little about the molecular characteristics of SAM cells, including their epigenome and how it changes during development. Here, we summarize information on epigenetic regulation of selected genes necessary for stem cell maintenance. As recent evidence suggests that SAM stem cells might be a hotspot of transposon activation, we discuss this aspect of epigenetic control in the meristem and speculate on mechanisms that maintain the flexibility of SAM stem cells in response to developmental or environmental cues.