Completeness Of Genome Assembly Research Articles

A chromosome-scale genome assembly of the nipa palm hispid beetle Octodonta nipae

Nipa palm hispid beetle (Octodonta nipae) is an insect species that is native to Malaysia but has spread to southern China and beyond, seriously threatening palm production. A lack of high-quality genome resources has hindered understanding of the insect’s invasive characteristics and ecological adaptations. Here, we combined Illumina short read, PacBio long-read, and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to generate a high-quality, chromosome-scale genome assembly of nipa palm hispid beetle. The genome assembly was 1.31 Gb in size, consisting of nine chromosomes. The contig and scaffold N50 values were 1.022 Mb and 148.6 Mb, respectively. The genome assembly completeness was estimated at 99.1%. Annotation revealed 16,305 protein-coding genes and 62.16% repeat sequences. This high-quality genome assembly is a valuable resource that will contribute to understanding of the genetic factors underlying the invasive characteristics of nipa palm hispid beetle, ultimately promoting development of efficient control policies.

The genome of Citrus australasica reveals disease resistance and other species specific genes

BackgroundThe finger lime (Citrus australasica), one of six Australian endemic citrus species shows a high natural phenotypic diversity and novel characteristics. The wide variation and unique horticultural features have made this lime an attractive candidate for domestication. Currently no haplotype resolved genome is available for this species. Here we present a high quality, haplotype-resolved reference genome for this species using PacBio HiFi and Hi-C sequencing.ResultsHifiasm assembly and SALSA scaffolding resulted in a collapsed genome size of 344.2 Mb and 321.1 Mb and 323.2 Mb size for the two haplotypes. The nine pseudochromosomes of the collapsed genome had an N50 of 35.2 Mb, 99.1% genome assembly completeness and 98.9% gene annotation completeness (BUSCO). A total of 41,304 genes were predicted in the nuclear genome. Comparison with C. australis revealed that 13,661 genes in pseudochromosomes were unique in C. australasica. These were mainly involved in plant-pathogen interactions, stress response, cellular metabolic and developmental processes, and signal transduction. The two genomes showed a syntenic arrangement at the chromosome level with large structural rearrangements in some chromosomes. Genetic variation among five C. australasica cultivars was analysed. Genes related to defense, synthesis of volatile compounds and red/yellow coloration were identified in the genome. A major expansion of genes encoding thylakoid curvature proteins was found in the C. australasica genome.ConclusionsThe genome of C. australasica present in this study is of high quality and contiguity. This genome helps deepen our understanding of citrus evolution and reveals disease resistance and quality related genes with potential to accelerate the genetic improvement of citrus.

Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery.

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

Compleasm: a faster and more accurate reimplementation of BUSCO.

Evaluating the gene completeness is critical to measuring the quality of a genome assembly. An incomplete assembly can lead to errors in gene predictions, annotation, and other downstream analyses. Benchmarking Universal Single-Copy Orthologs (BUSCO) is a widely used tool for assessing the completeness of genome assembly by testing the presence of a set of single-copy orthologs conserved across a wide range of taxa. However, BUSCO is slow particularly for large genome assemblies. It is cumbersome to apply BUSCO to a large number of assemblies. Here, we present compleasm, an efficient tool for assessing the completeness of genome assemblies. Compleasm utilizes the miniprot protein-to-genome aligner and the conserved orthologous genes from BUSCO. It is 14 times faster than BUSCO for human assemblies and reports a more accurate completeness of 99.6% than BUSCO's 95.7%, which is in close agreement with the annotation completeness of 99.5% for T2T-CHM13. https://github.com/huangnengCSU/compleasm.

Benchmarking multi-platform sequencing technologies for human genome assembly.

Genome assembly is a computational technique that involves piecing together deoxyribonucleic acid (DNA) fragments generated by sequencing technologies to create a comprehensive and precise representation of the entire genome. Generating a high-quality human reference genome is a crucial prerequisite for comprehending human biology, and it is also vital for downstream genomic variation analysis. Many efforts have been made over the past few decades to create a complete and gapless reference genome for humans by using a diverse range of advanced sequencing technologies. Several available tools are aimed at enhancing the quality of haploid and diploid human genome assemblies, which include contig assembly, polishing of contig errors, scaffolding and variant phasing. Selecting the appropriate tools and technologies remains a daunting task despite several studies have investigated the pros and cons of different assembly strategies. The goal of this paper was to benchmark various strategies for human genome assembly by combining sequencing technologies and tools on two publicly available samples (NA12878 and NA24385) from Genome in a Bottle. We then compared their performances in terms of continuity, accuracy, completeness, variant calling and phasing. We observed that PacBio HiFi long-reads are the optimal choice for generating an assembly with low base errors. On the other hand, we were able to produce the most continuous contigs with Oxford Nanopore long-reads, but they may require further polishing to improve on quality. We recommend using short-reads rather than long-reads themselves to improve the base accuracy of contigs from Oxford Nanopore long-reads. Hi-C is the best choice for chromosome-level scaffolding because it can capture the longest-range DNA connectedness compared to 10× linked-reads and Bionano optical maps. However, a combination of multiple technologies can be used to further improve the quality and completeness of genome assembly. For diploid assembly, hifiasm is the best tool for human diploid genome assembly using PacBio HiFi and Hi-C data. Looking to the future, we expect that further advancements in human diploid assemblers will leverage the power of PacBio HiFi reads and other technologies with long-range DNA connectedness to enable the generation of high-quality, chromosome-level and haplotype-resolved human genome assemblies.

Accurate microRNA annotation of animal genomes using trained covariance models of curated microRNA complements in MirMachine

The annotation of microRNAs depends on the availability of transcriptomics data and expert knowledge. This has led to a gap between the availability of novel genomes and high-quality microRNA complements. Using >16,000 microRNAs from the manually curated microRNA gene database MirGeneDB, we generated trained covariance models for all conserved microRNA families. These models are available in our tool MirMachine, which annotates conserved microRNAs within genomes. We successfully applied MirMachine to a range of animal species, including those with large genomes and genome duplications and extinct species, where small RNA sequencing is hard to achieve. We further describe a microRNA score of expected microRNAs that can be used to assess the completeness of genome assemblies. MirMachine closes a long-persisting gap in the microRNA field by facilitating automated genome annotation pipelines and deeper studies into the evolution of genome regulation, even in extinct organisms.

Gapless provides combined scaffolding, gap filling, and assembly correction with long reads.

Continuity, correctness, and completeness of genome assemblies are important for many biological projects. Long reads represent a major driver towards delivering high-quality genomes, but not everybody can achieve the necessary coverage for good long read-only assemblies. Therefore, improving existing assemblies with low-coverage long reads is a promising alternative. The improvements include correction, scaffolding, and gap filling. However, most tools perform only one of these tasks and the useful information of reads that supported the scaffolding is lost when running separate programs successively. Therefore, we propose a new tool for combined execution of all three tasks using PacBio or Oxford Nanopore reads. gapless is available at: https://github.com/schmeing/gapless.

Haplotype resolved chromosome level genome assembly of Citrus australis reveals disease resistance and other citrus specific genes.

Recent advances in genome sequencing and assembly techniques have made it possible to achieve chromosome level reference genomes for citrus. Relatively few genomes have been anchored at the chromosome level and/or are haplotype phased, with the available genomes of varying accuracy and completeness. We now report a phased high-quality chromosome level genome assembly for an Australian native citrus species; Citrus australis (round lime) using highly accurate PacBio HiFi long reads, complemented with Hi-C scaffolding. Hifiasm with Hi-C integrated assembly resulted in a 331Mb genome of C. australis with two haplotypes of nine pseudochromosomes with an N50 of 36.3Mb and 98.8% genome assembly completeness (BUSCO). Repeat analysis showed that more than 50% of the genome contained interspersed repeats. Among them, LTR elements were the predominant type (21.0%), of which LTR Gypsy (9.8%) and LTR copia (7.7%) elements were the most abundant repeats. A total of 29 464 genes and 32 009 transcripts were identified in the genome. Of these, 28 222 CDS (25 753 genes) had BLAST hits and 21 401 CDS (75.8%) were annotated with at least one GO term. Citrus specific genes for antimicrobial peptides, defense, volatile compounds and acidity regulation were identified. The synteny analysis showed conserved regions between the two haplotypes with some structural variations in Chromosomes 2, 4, 7 and 8. This chromosome scale, and haplotype resolved C. australis genome will facilitate the study of important genes for citrus breeding and will also allow the enhanced definition of the evolutionary relationships between wild and domesticated citrus species.

JASPER: A fast genome polishing tool that improves accuracy of genome assemblies.

Advances in long-read sequencing technologies have dramatically improved the contiguity and completeness of genome assemblies. Using the latest nanopore-based sequencers, we can generate enough data for the assembly of a human genome from a single flow cell. With the long-read data from these sequences, we can now routinely produce de novo genome assemblies in which half or more of a genome is contained in megabase-scale contigs. Assemblies produced from nanopore data alone, though, have relatively high error rates and can benefit from a process called polishing, in which more-accurate reads are used to correct errors in the consensus sequence. In this manuscript, we present a novel tool for genome polishing called JASPER (Jellyfish-based Assembly Sequence Polisher for Error Reduction). In contrast to many other polishing methods, JASPER gains efficiency by avoiding the alignment of reads to the assembly. Instead, JASPER uses a database of k-mer counts that it creates from the reads to detect and correct errors in the consensus. Our experiments demonstrate that JASPER is faster than alignment-based polishers, and both faster and more accurate than other k-mer based polishing methods. We also introduce the idea of using a polishing tool to create population-specific reference genomes, and illustrate this idea using sequence data from multiple individuals from Tokyo, Japan.

Chromosome-Level Genome Assembly of the Speckled Blue Grouper (Epinephelus cyanopodus) Provides Insight into Its Adaptive Evolution.

Epinephelus cyanopodus is a coral reef-dwelling grouper with important economic and ecological value and is widely distributed in the western Pacific Ocean. The lack of genomic resources for E. cyanopodus hinders its adaptive evolution and phylogeny research. We constructed the first high-quality genome of E. cyanopodus based on DNBSEQ, PacBio, and Hic sequencing technologies, with a genome size of 998.82 Mb, contig N50 of 5.855 Mb, and scaffold N50 of 41.98 Mb. More than 99.7% of contigs were anchored to 24 pseudochromosomes, and 94.2% of BUSCO genes were found in the E. cyanopodus genome, indicating a high genome assembly completeness. A total of 26,337 protein-coding genes were predicted, of which 98.77% were functionally annotated. Phylogenetic analysis showed that E. cyanopodus separated from its closely related species Epinephelus akaara about 11.5-26.5 million years ago, and the uplift of the Indo-Australian archipelago may have provided an opportunity for its rapid radiation. Moreover, several gene families associated with innate and adaptive immunity were significantly expanded in speckled blue grouper compared to other teleost genomes. Additionally, we identified several genes associated with immunity, growth and reproduction that are under positive selection in E. cyanopodus compared to other groupers, suggesting that E. cyanopodus has evolved broad adaptability in response to complex survival environment, which may provide the genetic basis for its rapid radiation. In brief, the high-quality reference genome of the speckled blue grouper provides a foundation for research on its biological traits and adaptive evolution and will be an important genetic tool to guide aquaculture and resolve its taxonomic controversies in future studies.

High-Quality Genome Assembly and Annotation Resource of Three Botryosphaeria Pathogens Causing Chinese Hickory Canker.

Tree trunk canker is the most critical disease threatening the Chinese hickory (Carya cathayensis). Except for the well-known pathogen Botryosphaeria dothidea, our group has previously reported three new pathogenic species, including B. corticis, B. fabicerciana, and B. qingyuanensis. Here, we reported almost chromosome-scale genome assemblies of the representative strains of these Botryosphaeria species using the Oxford Nanopore Technologies (ONT) sequencing platform. The assemblies were 44 to 46 Mb in size, and contained less than 14 contigs (13 - 14) with the contig N50 of around 4 Mb. The BUSCO completeness of genome assemblies was around 99% at the Ascomycota level. These genome assemblies contained no more than 8% repeats and around 13,000 (12,876 - 13,677) protein-coding genes. Gene functional annotations, including pathogen-host interaction, carbohydrate active enzymes (CAZymes), cytochrome P450 enzymes (CYPs), putative secreted proteins, and secondary metabolite biosynthesis gene clusters (SMBGCs) were performed. Comparative genomics analysis shows that these stains shared 11,489 core orthologous gene clusters with our previously reported B. dothidea strain BDLA16-7. The almost chromosomal-level and well-annotated genome assemblies will provide valuable comparative genomic resources for Botryosphaeria species.

High-Quality Genome Assembly and Annotation Resource of Botryosphaeria dothidea Strain BDLA16-7, Causing Trunk Canker Disease on Chinese Hickory.

Botryosphaeria dothidea is a latent pathogen with global importance to woody plant health, which causes serious tree trunk cankers on Chinese hickory. To date, only one Illumina short-read-based genome assembly of strain CK16 is available for host Chinese hickory. To address this problem, we reported a near telomere-to-telomere genome assembly of strain BDLA16-7 (46.05 Mb, N50 3.87 Mb) using Oxford Nanopore Sequencing Technology. Our genome assembly was consisted of 15 contigs, of which, 3 were assembled into chromosomal level and the maximum contig length was 6.19 Mb. The assembly contained 7.96% repeats and 12,815 protein-coding genes (10,274 genes were functional annotated). We also identified 3,642 pathogen-host interaction (PHI) genes, 250 carbohydrateactive enzymes (CAZymes), 252 cytochrome P450 enzymes (CYPs), 752 putative secreted proteins and 63 secondary metabolite biosynthesis gene clusters (SMBGCs). The BUSCO completeness of genome assembly and predicted genes was 99.34% and 97.50%, respectively, at fungal level (n=758). The almost chromosomal-level and well-annotated genome assembly will provide a valuable genetic resource for understanding of the infection mechanisms of B. dothidea in future.

A chromosome-level genome assembly of Cairina moschata and comparative genomic analyses

BackgroundThe Muscovy duck (Cairina moschata) is an economically important duck species, with favourable growth and carcass composition parameters in comparison to other ducks. However, limited genomic resources for Muscovy duck hinder our understanding of its evolution and genetic diversity.ResultsWe combined linked-reads sequencing technology and reference-guided methods for de novo genome assembly. The final draft assembly was 1.12 Gbp with 29 autosomes, one sex chromosome and 4,583 unlocalized scaffolds with an N50 size of 77.35 Mb. Based on universal single-copy orthologues (BUSCO), the draft genome assembly completeness was estimated to be 93.30 %. Genome annotation identified 15,580 genes, with 15,537 (99.72 %) genes annotated in public databases. We conducted comparative genomic analyses and found that species-specific and rapidly expanding gene families (compared to other birds) in Muscovy duck are mainly involved in Calcium signaling, Adrenergic signaling in cardiomyocytes, and GnRH signaling pathways. In comparison to the common domestic duck (Anas platyrhynchos), we identified 104 genes exhibiting strong signals of adaptive evolution (Ka/Ks > 1). Most of these genes were associated with immune defence pathways (e.g. IFNAR1 and TLR5). This is indicative of the existence of differences in the immune responses between the two species. Additionally, we combined divergence and polymorphism data to demonstrate the “faster-Z effect” of chromosome evolution.ConclusionsThe chromosome-level genome assembly of Muscovy duck and comparative genomic analyses provide valuable resources for future molecular ecology studies, as well as the evolutionary arms race between the host and influenza viruses.

Nanopore Sequencing and Hi-C Based De Novo Assembly of Trachidermus fasciatus Genome.

Trachidermus fasciatus is a roughskin sculpin fish widespread across the coastal areas of East Asia. Due to environmental destruction and overfishing, the population of this species is under threat. In order to protect this endangered species, it is important to have the genome sequenced. Reference genomes are essential for studying population genetics, domestic farming, and genetic resource protection. However, currently, no reference genome is available for Trachidermus fasciatus, and this has greatly hindered the research on this species. In this study, we integrated nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods to thoroughly assemble the Trachidermus fasciatus genome. Our results provided a chromosome-level high-quality genome assembly with a predicted genome size of 542.6 Mbp (2n = 40) and a scaffold N50 of 24.9 Mbp. The BUSCO value for genome assembly completeness was higher than 96%, and the single-base accuracy was 99.997%. Based on EVM-StringTie genome annotation, a total of 19,147 protein-coding genes were identified, including 35,093 mRNA transcripts. In addition, a novel gene-finding strategy named RNR was introduced, and in total, 51 (82) novel genes (transcripts) were identified. Lastly, we present here the first reference genome for Trachidermus fasciatus; this sequence is expected to greatly facilitate future research on this species.

Genome Sequence Resource of Phytophthora vignae, the Causal Agent of Stem and Root Rot of Cowpea.

The causal agent of stem and root rot of cowpea, Phytophthora vignae, is a widely distributed species of the Phytophthora genus. Here, we generate a high-quality complete genome assembly of P. vignae PSY2020 (89.39 Mb, N50 2.99 Mb) from China, using Oxford Nanopore Technologies (ONT) sequencing. The genome assembly completeness as evaluated by benchmarking universal single-copy orthologs was 94.51% at the eukaryote level. We identified 42.54% as repeat sequences and a total of 20,536 protein-encoding genes, of which 15,184 genes could be annotated. And we also identified 924 candidate RXLR effectors in the genome assembly. The described genome sequence will provide a valuable resource for better understanding of pathogenicity mechanisms of P. vignae and help in uncovering phylogenetical classification of Phytophthora species.[Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2021.

A first draft genome of the Sugarcane borer, Diatraea saccharalis.

Background: The sugarcane borer (Diatraea saccharalis), a widely distributed moth throughout the Americas, is a pest that affects economically important crops such as sugarcane, sorghum, wheat, maize and rice. Given its significant impact on yield reduction, whole-genome information of the species is needed. Here, we report the first draft assembly of the D. saccharalis genome. Methods: The genomic sequences were obtained using the Illumina HiSeq 2500 whole-genome sequencing of a single adult male specimen. We assembled the short-reads using the SPAdes software and predicted protein-coding genes using MAKER. Genome assembly completeness was assessed through BUSCO and the repetitive content by RepeatMasker. Results: The 453 Mb assembled sequences contain 1,445 BUSCO gene orthologs and 1,161 predicted gene models identified based on homology evidence to the domestic silk moth, Bombyx mori. The repeat content composes 41.18% of the genomic sequences which is in the range of other lepidopteran species. Conclusions: Functional annotation reveals that predicted gene models are involved in important cellular mechanisms such as metabolic pathways and protein synthesis. Thus, the data generated in this study expands our knowledge on the genomic characteristics of this devastating pest and provides essential resources for future genetic studies of the species.

An evaluation of methodology to determine algal genome completeness

The advancement of sequencing technologies has resulted in a rapid expansion of genome sequencing. One challenge in processing and analyzing this abundance of genomics data is development and application of tools to accurately assess the quality of novel genomes. The quality of a given genome assembly is commonly assessed using the presence of conserved orthologs. However, the applicability and accuracy of identifying conserved orthologs applied to algal genome assemblies, which are uniquely diverse, is unclear. Here, we analyze the utility of a genome content assessment tool, Benchmarking Universal Single-Copy Orthologs (BUSCO), to evaluate algal genome assembly completeness. We find that Chlorophyta and Stramenopile BUSCO databases are effective tools to analyze genome completeness of these lineages. For other lineages, the Eukaryota database should be utilized until additional lineage-specific databases are developed. Based on these results, suggested best practices for current usage and necessary future improvements are provided.

De novo Assembly of the Brugia malayi Genome Using Long Reads from a Single MinION Flowcell

Filarial nematode infections cause a substantial global disease burden. Genomic studies of filarial worms can improve our understanding of their biology and epidemiology. However, genomic information from field isolates is limited and available reference genomes are often discontinuous. Single molecule sequencing technologies can reduce the cost of genome sequencing and long reads produced from these devices can improve the contiguity and completeness of genome assemblies. In addition, these new technologies can make generation and analysis of large numbers of field isolates feasible. In this study, we assessed the performance of the Oxford Nanopore Technologies MinION for sequencing and assembling the genome of Brugia malayi, a human parasite widely used in filariasis research. Using data from a single MinION flowcell, a 90.3 Mb nuclear genome was assembled into 202 contigs with an N50 of 2.4 Mb. This assembly covered 96.9% of the well-defined B. malayi reference genome with 99.2% identity. The complete mitochondrial genome was obtained with individual reads and the nearly complete genome of the endosymbiotic bacteria Wolbachia was assembled alongside the nuclear genome. Long-read data from the MinION produced an assembly that approached the quality of a well-established reference genome using comparably fewer resources.

SLR: a scaffolding algorithm based on long reads and contig classification

BackgroundScaffolding is an important step in genome assembly that orders and orients the contigs produced by assemblers. However, repetitive regions in contigs usually prevent scaffolding from producing accurate results. How to solve the problem of repetitive regions has received a great deal of attention. In the past few years, long reads sequenced by third-generation sequencing technologies (Pacific Biosciences and Oxford Nanopore) have been demonstrated to be useful for sequencing repetitive regions in genomes. Although some stand-alone scaffolding algorithms based on long reads have been presented, scaffolding still requires a new strategy to take full advantage of the characteristics of long reads.ResultsHere, we present a new scaffolding algorithm based on long reads and contig classification (SLR). Through the alignment information of long reads and contigs, SLR classifies the contigs into unique contigs and ambiguous contigs for addressing the problem of repetitive regions. Next, SLR uses only unique contigs to produce draft scaffolds. Then, SLR inserts the ambiguous contigs into the draft scaffolds and produces the final scaffolds. We compare SLR to three popular scaffolding tools by using long read datasets sequenced with Pacific Biosciences and Oxford Nanopore technologies. The experimental results show that SLR can produce better results in terms of accuracy and completeness. The open-source code of SLR is available at https://github.com/luojunwei/SLR.ConclusionIn this paper, we describes SLR, which is designed to scaffold contigs using long reads. We conclude that SLR can improve the completeness of genome assembly.

A hybrid de novo genome assembly of the honeybee, Apis mellifera, with chromosome-length scaffolds

BackgroundThe ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybee Apis mellifera. We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map.ResultsEach of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor > 98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds of AvaI and AluI repeats associated with these features.ConclusionsThe improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics.