DEGAP: Dynamic elongation of a genome assembly path.

Genome assembly remains to be a major task in genomic research. Despite the development over the past decades of different assembly software programs and algorithms, it is still a great challenge to assemble a complete genome without any gaps. With the latest DNA circular consensus sequencing (CCS) technology, several assembly programs can now build a genome from raw sequencing data to contigs; however, some complex sequence regions remain as unresolved gaps. Here, we present a novel gap-filling software, DEGAP (Dynamic Elongation of a Genome Assembly Path), that resolves gap regions by utilizing the dual advantages of accuracy and length of high-fidelity (HiFi) reads. DEGAP identifies differences between reads and provides 'GapFiller' or 'CtgLinker' modes to eliminate or shorten gaps in genomes. DEGAP adopts an iterative elongation strategy that automatically and dynamically adjusts parameters according to three complexity factors affecting the genome to determine the optimal extension path. DEGAP has already been successfully applied to decipher complex genomic regions in several projects and may be widely employed to generate more gap-free genomes.

Three near-complete genome assemblies reveal substantial centromere dynamics from diploid to tetraploid in Brachypodium genus.

BACKGROUND: Centromeres are critical for maintaining genomic stability in eukaryotes, and their turnover shapes genome architectures and drives karyotype evolution. However, the co-evolution of centromeres from different species in allopolyploids over millions of years remains largely unknown. RESULTS: Here, we generate three near-complete genome assemblies, a tetraploid Brachypodium hybridum and its two diploid ancestors, Brachypodium distachyon and Brachypodium stacei. We detect high degrees of sequence, structural, and epigenetic variations of centromeres at base-pair resolution between closely related Brachypodium genomes, indicating the appearance and accumulation of species-specific centromere repeats from a common origin during evolution. We also find that centromere homogenization is accompanied by local satellite repeats bursting and retrotransposon purging, and the frequency of retrotransposon invasions drives the degree of interspecies centromere diversification. We further investigate the dynamics of centromeres during alloploidization process, and find that dramatic genetics and epigenetics architecture variations are associated with the turnover of centromeres between homologous chromosomal pairs from diploid to tetraploid. Additionally, our pangenomes analysis reveals the ongoing variations of satellite repeats and stable evolutionary homeostasis within centromeres among individuals of each Brachypodium genome with different polyploidy levels. CONCLUSIONS: Our results provide unprecedented information on the genomic, epigenomic, and functional diversity of highly repetitive DNA between closely related species and their allopolyploid genomes at both coarse and fine scale.

Dissecting the molecular basis of spike traits by integrating gene regulatory networks and genetic variation in wheat.

Spike architecture influences both grain weight and grain number per spike, which are the two major components of grain yield in bread wheat (Triticum aestivum L.). However, the complex wheat genome and the influence of various environmental factors pose challenges in mapping the causal genes that affect spike traits. Here, we systematically identified genes involved in spike trait formation by integrating information on genomic variation and gene regulatory networks controlling young spike development in wheat. We identified 170 loci that are responsible for variations in spike length, spikelet number per spike, and grain number per spike through genome-wide association study and meta-QTL analyses. We constructed gene regulatory networks for young inflorescences at the double ridge stage and the floret primordium stage, in which the spikelet meristem and the floret meristem are predominant, respectively, by integrating transcriptome, histone modification, chromatin accessibility, eQTL, and protein-protein interactome data. From these networks, we identified 169 hub genes located in 76 of the 170 QTL regions whose polymorphisms are significantly associated with variation in spike traits. The functions of TaZF-B1, VRT-B2, and TaSPL15-A/D in establishment of wheat spike architecture were verified. This study provides valuable molecular resources for understanding spike traits and demonstrates that combining genetic analysis and developmental regulatory networks is a robust approach for dissection of complex traits.

Chemical-tag-based semi-annotated metabolomics facilitates gene identification and specialized metabolic pathway elucidation in wheat.

The importance of metabolite modification and species-specific metabolic pathways has long been recognized. However, linking the chemical structure of metabolites to gene function in order to explore the genetic and biochemical basis of metabolism has not yet been reported in wheat (Triticum aestivum). Here, we profiled metabolic fragment enrichment in wheat leaves and consequently applied chemical-tag-based semi-annotated metabolomics in a genome-wide association study in accessions of wheat. The studies revealed that all 1,483 quantified metabolites have at least one known functional group whose modification is tailored in an enzyme-catalyzed manner and eventually allows efficient candidate gene mining. A Triticeae crop-specific flavonoid pathway and its underlying metabolic gene cluster were elucidated in further functional studies. Additionally, upon overexpressing the major effect gene of the cluster TraesCS2B01G460000 (TaOMT24), the pathway was reconstructed in rice (Oryza sativa), which lacks this pathway. The reported workflow represents an efficient and unbiased approach for gene mining using forward genetics in hexaploid wheat. The resultant candidate gene list contains vast molecular resources for decoding the genetic architecture of complex traits and identifying valuable breeding targets and will ultimately aid in achieving wheat crop improvement.

Centromere Plasticity With Evolutionary Conservation and Divergence Uncovered by Wheat 10+ Genomes.

Centromeres (CEN) are the chromosomal regions that play a crucial role in maintaining genomic stability. The underlying highly repetitive DNA sequences can evolve quickly in most eukaryotes, and promote karyotype evolution. Despite their variability, it is not fully understood how these widely variable sequences ensure the homeostasis of centromere function. In this study, we investigated the genetics and epigenetics of CEN in a population of wheat lines from global breeding programs. We captured a high degree of sequences, positioning, and epigenetic variations in the large and complex wheat CEN. We found that most CENH3-associated repeats are Cereba element of retrotransposons and exhibit phylogenetic homogenization across different wheat lines, but the less-associated repeat sequences diverge on their own way in each wheat line, implying specific mechanisms for selecting certain repeat types as functional core CEN. Furthermore, we observed that CENH3 nucleosome structures display looser wrapping of DNA termini on complex centromeric repeats, including the repositioned CEN. We also found that strict CENH3 nucleosome positioning and intrinsic DNA features play a role in determining centromere identity among different lines. Specific non-B form DNAs were substantially associated with CENH3 nucleosomes for the repositioned centromeres. These findings suggest that multiple mechanisms were involved in the adaptation of CENH3 nucleosomes that can stabilize CEN. Ultimately, we proposed a remarkable epigenetic plasticity of centromere chromatin within the diverse genomic context, and the high robustness is crucial for maintaining centromere function and genome stability in wheat 10+ lines as a result of past breeding selections.

Deciphering genetic basis of developmental and agronomic traits by integrating high-throughput optical phenotyping and genome-wide association studies in wheat.

Dissecting the genetic basis of complex traits such as dynamic growth and yield potential is a major challenge in crops. Monitoring the growth throughout growing season in a large wheat population to uncover the temporal genetic controls for plant growth and yield-related traits has so far not been explored. In this study, a diverse wheat panel composed of 288 lines was monitored by a non-invasive and high-throughput phenotyping platform to collect growth traits from seedling to grain filling stage and their relationship with yield-related traits was further explored. Whole genome re-sequencing of the panel provided 12.64 million markers for a high-resolution genome-wide association analysis using 190 image-based traits and 17 agronomic traits. A total of 8327 marker-trait associations were detected and clustered into 1605 quantitative trait loci (QTLs) including a number of known genes or QTLs. We identified 277 pleiotropic QTLs controlling multiple traits at different growth stages which revealed temporal dynamics of QTLs action on plant development and yield production in wheat. A candidate gene related to plant growth that was detected by image traits was further validated. Particularly, our study demonstrated that the yield-related traits are largely predictable using models developed based on i-traits and provide possibility for high-throughput early selection, thus to accelerate breeding process. Our study explored the genetic architecture of growth and yield-related traits by combining high-throughput phenotyping and genotyping, which further unravelled the complex and stage-specific contributions of genetic loci to optimize growth and yield in wheat.

New insights on the evolution of nucleolar dominance in newly resynthesized hexaploid wheat Triticum zhukovskyi.

Nucleolar dominance (ND) is a widespread epigenetic phenomenon in hybridizations where nucleolus transcription fails at the nucleolus organizer region (NOR). However, the dynamics of NORs during the formation of Triticum zhukovskyi (GGAu Au Am Am ), another evolutionary branch of allohexaploid wheat, remains poorly understood. Here, we elucidated genetic and epigenetic changes occurring at the NOR loci within the Am , G, and D subgenomes during allopolyploidization by synthesizing hexaploid wheat GGAu Au Am Am and GGAu Au DD. In T. zhukovskyi, Au genome NORs from T. timopheevii (GGAu Au ) were lost, while the second incoming NORs from T. monococcum (Am Am ) were retained. Analysis of the synthesized T. zhukovskyi revealed that rRNA genes from the Am genome were silenced in F1 hybrids (GAu Am ) and remained inactive after genome doubling and subsequent self-pollinations. We observed increased DNA methylation accompanying the inactivation of NORs in the Am genome and found that silencing of NORs in the S1 generation could be reversed by a cytidine methylase inhibitor. Our findings provide insights into the ND process during the evolutionary period of T. zhukovskyi and highlight that inactive rDNA units may serve as a 'first reserve' in the form of R-loops, contributing to the successful evolution of T. zhukovskyi.

Asymmetric gene expression and cell-type-specific regulatory networks in the root of bread wheat revealed by single-cell multiomics analysis.

BACKGROUND: Homoeologs are defined as homologous genes resulting from allopolyploidy. Bread wheat, Triticum aestivum, is an allohexaploid species with many homoeologs. Homoeolog expression bias, referring to the relative contribution of homoeologs to the transcriptome, is critical for determining the traits that influence wheat growth and development. Asymmetric transcription of homoeologs has been so far investigated in a tissue or organ-specific manner, which could be misleading due to a mixture of cell types. RESULTS: Here, we perform single nuclei RNA sequencing and ATAC sequencing of wheat root to study the asymmetric gene transcription, reconstruct cell differentiation trajectories and cell-type-specific gene regulatory networks. We identify 22 cell types. We then reconstruct cell differentiation trajectories that suggest different origins between epidermis/cortex and endodermis, distinguishing bread wheat from Arabidopsis. We show that the ratio of asymmetrically transcribed triads varies greatly when analyzing at the single-cell level. Hub transcription factors determining cell type identity are also identified. In particular, we demonstrate that TaSPL14 participates in vasculature development by regulating the expression of BAM1. Combining single-cell transcription and chromatin accessibility data, we construct the pseudo-time regulatory network driving root hair differentiation. We find MYB3R4, REF6, HDG1, and GATAs as key regulators in this process. CONCLUSIONS: Our findings reveal the transcriptional landscape of root organization and asymmetric gene transcription at single-cell resolution in polyploid wheat.

Systemic development of wheat-Thinopyrum elongatum translocation lines and their deployment in wheat breeding for Fusarium head blight resistance.

Fusarium head blight (FHB), mainly caused by Fusarium graminearum, is one of the most destructive diseases of wheat (Triticum aestivum) around the world. FHB causes significant yield losses and reduces grain quality. The lack of resistance resources is a major bottleneck for wheat FHB resistance breeding. As a wheat relative, Thinopyrum elongatum contains many genes that can be used for wheat improvement. Although the novel gene Fhb-7EL was mapped on chromosome 7EL of Th. elongatum, successful transfer of the FHB resistance gene into commercial wheat varieties has not been reported. In this study, we developed 836 wheat-Th. elongatum translocation lines of various types by irradiating the pollen of the wheat-Th. elongatum addition line CS-7EL at the flowering stage, among which 81 were identified as resistant to FHB. By backcrossing the FHB-resistant lines with the main cultivar Jimai 22, three wheat-Th. elongatum translocation lines, Zhongke 1878, Zhongke 166, and Zhongke 545, were successfully applied in wheat breeding without yield penalty. Combining karyotype and phenotype analyses, we mapped the Fhb-7EL gene to the distal end of chromosome 7EL. Five molecular markers linked with the FHB resistance interval were developed, which facilitates molecular marker-assisted breeding. Altogether, we successfully applied alien chromatin with FHB resistance from Th. elongatum in wheat breeding without yield penalty. These newly developed FHB-resistant wheat-Th. elongatum translocation lines, Zhongke 1878, Zhongke 166, and Zhongke 545, can be used as novel resistance resources for wheat breeding.

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

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

Non-B-form DNA tends to form in centromeric regions and has undergone changes in polyploid oat subgenomes.

Centromeres are the specialized regions of the chromosomes that direct faithful chromosome segregation during cell division. Despite their functional conservation, centromeres display features of rapidly evolving DNA and wide evolutionary diversity in size and organization. Previous work found that the noncanonical B-form DNA structures are abundant in the centromeres of several eukaryotic species with a possible implication for centromere specification. Thus far, systematic studies into the organization and function of non-B-form DNA in plants remain scarce. Here, we applied the oat system to investigate the role of non-B-form DNA in centromeres. We conducted chromatin immunoprecipitation sequencing using an antibody to the centromere-specific histone H3 variant (CENH3); this accurately positioned oat centromeres with different ploidy levels and identified a series of centromere-specific sequences including minisatellites and retrotransposons. To define genetic characteristics of oat centromeres, we surveyed the repeat sequences and found that dyad symmetries were abundant in oat centromeres and were predicted to form non-B-DNA structures in vivo. These structures including bent DNA, slipped DNA, Z-DNA, G-quadruplexes, and R-loops were prone to form within CENH3-binding regions. Dynamic conformational changes of predicted non-B-DNA occurred during the evolution from diploid to tetraploid to hexaploid oat. Furthermore, we applied the single-molecule technique of AFM and DNA:RNA immunoprecipitation with deep sequencing to validate R-loop enrichment in oat centromeres. Centromeric retrotransposons exhibited strong associations with R-loop formation. Taken together, our study elucidates the fundamental character of non-B-form DNA in the oat genome and reveals its potential role in centromeres.

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.

Genome-wide mapping reveals R-loops associated with centromeric repeats in maize.

R-loops are stable chromatin structures comprising a DNA:RNA hybrid and a displaced single-stranded DNA. R-loops have been implicated in gene expression and chromatin structure, as well as in replication blocks and genome instability. Here, we conducted a genome-wide identification of R-loops and identified more than 700,000 R-loop peaks in the maize (Zea mays) genome. We found that sense R-loops were mainly enriched in promoters and transcription termination sites and relatively less enriched in gene bodies, which is different from the main gene-body localization of sense R-loops in Arabidopsis and Oryza sativa At the chromosome scale, maize R-loops were enriched in pericentromeric heterochromatin regions, and a significant portion of R-loops were derived from transposable elements. In centromeres, R-loops preferentially formed within the binding regions of the centromere-specific histone CENH3, and centromeric retrotransposons were strongly associated with R-loop formation. Furthermore, centromeric retrotransposon R-loops were observed by applying the single-molecule imaging technique of atomic force microscopy. These findings elucidate the fundamental character of R-loops in the maize genome and reveal the potential role of R-loops in centromeres.

Sequence of the supernumerary B chromosome of maize provides insight into its drive mechanism and evolution.

B chromosomes are enigmatic elements in thousands of plant and animal genomes that persist in populations despite being nonessential. They circumvent the laws of Mendelian inheritance but the molecular mechanisms underlying this behavior remain unknown. Here we present the sequence, annotation, and analysis of the maize B chromosome providing insight into its drive mechanism. The sequence assembly reveals detailed locations of the elements involved with the cis and trans functions of its drive mechanism, consisting of nondisjunction at the second pollen mitosis and preferential fertilization of the egg by the B-containing sperm. We identified 758 protein-coding genes in 125.9 Mb of B chromosome sequence, of which at least 88 are expressed. Our results demonstrate that transposable elements in the B chromosome are shared with the standard A chromosome set but multiple lines of evidence fail to detect a syntenic genic region in the A chromosomes, suggesting a distant origin. The current gene content is a result of continuous transfer from the A chromosomal complement over an extended evolutionary time with subsequent degradation but with selection for maintenance of this nonvital chromosome.

Knl1 participates in spindle assembly checkpoint signaling in maize.

The Knl1-Mis12-Ndc80 (KMN) network is an essential component of the kinetochore-microtubule attachment interface, which is required for genomic stability in eukaryotes. However, little is known about plant Knl1 proteins because of their complex evolutionary history. Here, we cloned the Knl1 homolog from maize (Zea mays) and confirmed it as a constitutive central kinetochore component. Functional assays demonstrated their conserved role in chromosomal congression and segregation during nuclear division, thus causing defective cell division during kernel development when Knl1 transcript was depleted. A 145 aa region in the middle of maize Knl1, that did not involve the MELT repeats, was associated with the interaction of spindle assembly checkpoint (SAC) components Bub1/Mad3 family proteins 1 and 2 (Bmf1/2) but not with the Bmf3 protein. They may form a helical conformation with a hydrophobic interface with the TPR domain of Bmf1/2, which is similar to that of vertebrates. However, this region detected in monocots shows extensive divergence in eudicots, suggesting that distinct modes of the SAC to kinetochore connection are present within plant lineages. These findings elucidate the conserved role of the KMN network in cell division and a striking dynamic of evolutionary patterns in the SAC signaling and kinetochore network.

Emerging roles of centromeric RNAs in centromere formation and function.

BACKGROUND: Centromeres are specialized chromosomal domains involved in kinetochore formation and faithful chromosome segregation. Despite a high level of functional conservation, centromeres are not identified by DNA sequences, but by epigenetic means. Universally, centromeres are typically formed on highly repetitive DNA, which were previously considered to be silent. However, recent studies have shown that transcription occurs in this region, known as centromeric-derived RNAs (cenRNAs). CenRNAs that contribute to fundamental aspects of centromere function have been recently investigated in detail. However, the distribution, behavior and contributions of centromeric transcripts are still poorly understood. OBJECTIVE: The aim of this article is to provide an overview of the roles of cenRNAs in centromere formation and function. METHODS: We describe the structure and DNA sequence of centromere from yeast to human. In addition, we briefly introduce the roles of cenRNAs in centromere formation and function, kinetochore structure, accurate chromosome segregation, and pericentromeric heterochromatin assembly. Centromeric circular RNAs (circRNAs) and R-loops are rising stars in centromere function. CircRNAs have been successfully identified in various species with the assistance of high-throughput sequencing and novel computational approaches for non-polyadenylated RNA transcripts. Centromeric R-loops can be identified by the single-strand DNA ligation-based library preparation technique. But the molecular features and function of these centromeric R-loops and circRNAs are still being investigated. CONCLUSION: In this review, we summarize recent findings on the epigenetic regulation of cenRNAs across species, which would provide useful information about cenRNAs and interesting hints for further studies.

Phosphorylation of histone H3 by Haspin regulates chromosome alignment and segregation during mitosis in maize.

In human cells, Haspin-mediated histone H3 threonine 3 (H3T3) phosphorylation promotes centromeric localization of the chromosomal passenger complex, thereby ensuring proper kinetochore-microtubule attachment. Haspin also binds to PDS5 cohesin-associated factor B (Pds5B), antagonizing the Wings apart-like protein homolog (Wapl)-Pds5B interaction and thus preventing Wapl from releasing centromeric cohesion during mitosis. However, the role of Haspin in plant chromosome segregation is not well understood. Here, we show that in maize (Zea mays) mitotic cells, ZmHaspin localized to the centromere during metaphase and anaphase, whereas it localized to the telomeres during meiosis. These results suggest that ZmHaspin plays different roles during mitosis and meiosis. Knockout of ZmHaspin led to decreased H3T3 phosphorylation and histone H3 serine 10 phosphorylation, and defects in chromosome alignment and segregation in mitosis. These lines of evidence suggest that Haspin regulates chromosome segregation in plants via the mechanism described for humans, namely, H3T3 phosphorylation. Plant Haspin proteins lack the RTYGA and PxVxL motifs needed to bind Pds5B and heterochromatin protein 1, and no obvious cohesion defects were detected in ZmHaspin knockout plants. Taken together, these results highlight the conserved but slightly different roles of Haspin proteins in cell division in plants and in animals.