Prospects of telomere-to-telomere assembly in barley: analysis of sequence gaps in the MorexV3 reference genome

The first gapless, telomere-to-telomere (T2T) sequence assemblies of plant chromosomes were reported recently. However, sequence assemblies of most plant genomes remain fragmented. Only recent breakthroughs in accurate long-read sequencing have made it possible to achieve highly contiguous sequence assemblies with a few tens of contigs per chromosome, i.e. a number small enough to allow for a systematic inquiry into the causes of the remaining sequence gaps and the approaches and resources needed to close them. Here, we analyze sequence gaps in the current reference genome sequence of barley cv. Morex (MorexV3). Optical map and sequence raw data, complemented by ChIP-seq data for centromeric histone variant CENH3, were used to estimate the abundance of centromeric, ribosomal DNA and subtelomeric repeats in the barley genome. These estimates were compared with copy numbers in the MorexV3 pseudomolecule sequence. We found that almost all centromeric sequences and 45S ribosomal DNA repeat arrays were absent from the MorexV3 pseudomolecules and that the majority of sequence gaps can be attributed to assembly breakdown in long stretches of satellite repeats. However, missing sequences cannot fully account for the difference between assembly size and flow cytometric genome size estimates. We discuss the prospects of gap closure with ultra-long sequence reads.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of living amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., centric) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding, surrounded by pericentromeric LINE/L1 elements. We explored chromosome structure across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (HiC) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible association of centromeric chromatin, and of telomeres, reflecting a Rabl configuration similar to the "bouquet" structure of meiotic cells. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

BRCA1 prevents R-loop-associated centromeric instability

Centromeres are defined by chromatin containing the histone H3 variant CENP-A assembled onto repetitive α-satellite sequences, which are actively transcribed throughout the cell cycle. Centromeres play an essential role in chromosome inheritance and genome stability through coordinating kinetochores assembly during mitosis. Structural and functional alterations of the centromeres cause aneuploidy and chromosome aberrations which can induce cell death. In human cells, the tumor suppressor BRCA1 associates with centromeric chromatin in the absence of exogenous damage. While we previously reported that BRCA1 contributes to proper centromere homeostasis, the mechanism underlying its centromeric function and recruitment was not fully understood. Here, we show that BRCA1 association with centromeric chromatin depends on the presence of R-loops, which are non-canonical three-stranded structures harboring a DNA:RNA hybrid and are frequently formed during transcription. Subsequently, BRCA1 counteracts the accumulation of R-loops at centromeric α-satellite repeats. Strikingly, BRCA1-deficient cells show impaired localization of CENP-A, higher transcription of centromeric RNA, increased breakage at centromeres and formation of acentric micronuclei, all these features being R-loop-dependent. Finally, BRCA1 depletion reveals a Rad52-dependent hyper-recombination process between centromeric satellite repeats, associated with centromere instability and missegregation. Altogether, our findings provide molecular insights into the key function of BRCA1 in maintaining centromere stability and identity.

Loss of centromeric RNA activates the spindle assembly checkpoint in mammalian female meiosis I

The repetitive sequences of DNA centromeric regions form the structural basis for kinetochore assembly. Recently they were found to be transcriptionally active in mitosis, with their RNAs providing noncoding functions. Here we explore the role, in mouse oocytes, of transcripts generated from within the minor satellite repeats. Depletion of minor satellite transcripts delayed progression through meiosis I by activation of the spindle assembly checkpoint. Arrested oocytes had poorly congressed chromosomes, and centromeres were frequently split by microtubules. Thus, we have demonstrated that the centromeric RNA plays a specific role in female meiosis I compared with mitosis and is required for maintaining the structural integrity of centromeres. This may contribute to the high aneuploidy rates observed in female meiosis.

Red queen's race: rapid evolutionary dynamics of an expanding family of meiotic drive factors and their hpRNA suppressors

Meiotic drivers are a class of selfish genetic elements that are widespread across eukaryotes. Their activities are often detrimental to organismal fitness and thus trigger drive suppression to ensure fair segregation during meiosis. Accordingly, their existence is frequently hidden in genomes, and their molecular functions are little known. Here, we trace evolutionary steps that generated the Dox meiotic drive system in Drosophila simulans (Dsim), which distorts male:female balance (sex-ratio) by depleting male progeny. We show that Dox emerged via stepwise mobilization and acquisition of portions of multiple D. melanogaster genes, including the sperm chromatin packaging gene protamine. Moreover, we reveal novel Dox homologs in Dsim and massive, recent, amplification of Dox superfamily genes specifically on X chromosomes of its closest sister species D. mauritiana (Dmau) and D. sechellia (Dsech). The emergence of Dox superfamily genes is tightly associated with 1.688 family satellite repeats that flank de novo genomic copies. In concert, we find coordinated emergence and diversification of autosomal hairpin RNA/siRNAs loci that target subsets of Dox superfamily genes across simulans clade species. Finally, an independent set of protamine amplifications the Y chromosome of D. melanogaster indicates that protamine genes are frequent and recurrent players in sex chromosome dynamics. Overall, we reveal fierce genetic arms races between meiotic drive factors and siRNA suppressors associated with recent speciation.

Genome Assembly of the Popular Korean Soybean Cultivar Hwangkeum

Massive resequencing efforts have been undertaken to catalog allelic variants in major crop species including soybean, but the scope of the information for genetic variation often depends on short sequence reads mapped to the extant reference genome. Additional de novo assembled genome sequences provide a unique opportunity to explore a dispensable genome fraction in the pan-genome of a species. Here, we report the de novo assembly and annotation of Hwangkeum, a popular soybean cultivar in Korea. The assembly was constructed using PromethION nanopore sequencing data and two genetic maps, and was then error-corrected using Illumina short-reads and PacBio SMRT reads. The 933.12 Mb assembly was annotated as containing 79,870 transcripts for 58,550 genes using RNA-Seq data and the public soybean annotation set. Comparison of the Hwangkeum assembly with the Williams 82 soybean reference genome sequence (Wm82.a2.v1) revealed 1.8 million single-nucleotide polymorphisms, 0.5 million indels, and 25 thousand putative structural variants. However, there was no natural megabase-scale chromosomal rearrangement. Incidentally, by adding two novel subfamilies, we found that soybean contains four clearly separated subfamilies of centromeric satellite repeats. Analyses of satellite repeats and gene content suggested that the Hwangkeum assembly is a high-quality assembly. This was further supported by comparison of the marker arrangement of anthocyanin biosynthesis genes and of gene arrangement at the Rsv3 locus. Therefore, the results indicate that the de novo assembly of Hwangkeum is a valuable additional reference genome resource for characterizing traits for the improvement of this important crop species.

Impact of Repetitive DNA Elements on Snake Genome Biology and Evolution

The distinctive biology and unique evolutionary features of snakes make them fascinating model systems to elucidate how genomes evolve and how variation at the genomic level is interlinked with phenotypic-level evolution. Similar to other eukaryotic genomes, large proportions of snake genomes contain repetitive DNA, including transposable elements (TEs) and satellite repeats. The importance of repetitive DNA and its structural and functional role in the snake genome, remain unclear. This review highlights the major types of repeats and their proportions in snake genomes, reflecting the high diversity and composition of snake repeats. We present snakes as an emerging and important model system for the study of repetitive DNA under the impact of sex and microchromosome evolution. We assemble evidence to show that certain repetitive elements in snakes are transcriptionally active and demonstrate highly dynamic lineage-specific patterns as repeat sequences. We hypothesize that particular TEs can trigger different genomic mechanisms that might contribute to driving adaptive evolution in snakes. Finally, we review emerging approaches that may be used to study the expression of repetitive elements in complex genomes, such as snakes. The specific aspects presented here will stimulate further discussion on the role of genomic repeats in shaping snake evolution.

Nde1 is Required for Heterochromatin Compaction and Stability in Neocortical Neurons

The NDE1 gene encodes a scaffold protein essential for brain development. While biallelic NDE1 loss of function (LOF) causes microcephaly with profound mental retardation, NDE1 missense mutations and copy number variations are associated with multiple neuropsychiatric disorders. However, the etiology of the diverse phenotypes resulting from NDE1 aberrations remains elusive. Here we show Nde1 controls neurogenesis through heterochromatin compaction via histone H4K20 trimethylation. This mechanism patterns diverse chromatin landscapes and stabilizes constitutive heterochromatin of neocortical neurons. We show NDE1 undergoes dynamic liquid-liquid phase separation, partitioning to the nucleus and interacting with pericentromeric and centromeric satellite repeats. Nde1 LOF results in nuclear architecture aberrations, DNA double strand breaks, as well as instability and derepression of pericentromeric satellite repeats in neocortical neurons. These findings uncover a pivotal role of NDE1/Nde1 in establishing and maintaining neuronal heterochromatin. They suggest that heterochromatin impairments underlie a wide range of brain dysfunction.

High-quality Arabidopsis thaliana genome assembly with Nanopore and HiFi long reads

Here, we report a high-quality (HQ) and almost complete genome assembly with a single gap and quality value (QV) larger than 60 of the model plant Arabidopsis thaliana ecotype Columbia (Col-0), generated using combination of Oxford Nanopore Technology (ONT) ultra-long reads, high fidelity (HiFi) reads and Hi-C data. The total genome assembly size is 133,877,291 bp (chr1: 32,659,241 bp, chr2: 22,712,559 bp, chr3: 26,161,332 bp, chr4: 22,250,686 bp and chr5: 30,093,473 bp), and introduces 14.73 Mb (96% belong to centromere) novel sequences compared to TAIR10.1 reference genome. All five chromosomes of our HQ assembly are highly accurate with QV larger than 60, ranging from QV62 to QV68, which is significantly higher than TAIR10.1 referecne (44-51) and a recent published genome (41-43). We have completely resolved chr3 and chr5 from telomere-to-telomere. For chr2 and chr4, we have completely resolved apart from the nucleolar organizing regions, which are composed of highly long-repetitive DNA fragments. It has been reported that the length of centromere 1 is about 9 Mb and it is hard to assembly since tens of thousands of CEN180 satellite repeats. Based on the cutting-edge sequencing data, we assembled about 4Mb continuous sequence of centromere 1. We found different identity patterns across five centromeres, and all centromeres were significantly enriched with CENH3 ChIP-seq signals, confirming the accuracy of the assembly. We obtained four clusters of CEN180 repeats, and found CENH3 presented a strong preference for a cluster 3. Moreover, we observed hypomethylation patterns in CENH3 enriched regions. This high-quality assembly genome will be a valuable reference to assist us in the understanding of global pattern of centromeric polymorphism, genetic and epigenetic in naturally inbred lines of Arabidopsis thaliana.