The value of hybrid genomes: Building two highly contiguous reference genome assemblies to advance Canis genomic studies.

Previous studies of canid population and evolutionary genetics have relied on high-quality domestic dog reference genomes that have been produced primarily for biomedical and trait mapping studies in dog breeds. However, the absence of highly contiguous genomes from other Canis species like the gray wolf and coyote, that represent additional distinct demographic histories, may bias inferences regarding interspecific genetic diversity and phylogenetic relationships. Here, we present single haplotype de novo genome assemblies for the gray wolf and coyote, generated by applying the trio-binning approach to long sequence reads generated from the genome of a female first-generation hybrid produced from a gray wolf and coyote mating. The assemblies were highly contiguous, with contig N50 sizes of 44.6 and 42.0 Mb for the wolf and coyote, respectively. Genome scaffolding and alignments between the two Canis assemblies and published dog reference genomes showed near complete collinearity, with one exception: a coyote-specific chromosome fission of chromosome 13 and fusion of the proximal portion of that chromosome with chromosome 8, retaining the Canis-typical haploid chromosome number of 2n = 78. We evaluated mapping quality for previous RADseq data from 334 canids and found nearly identical mapping quality and patterns among canid species and regional populations regardless of the genome used for alignment (dog, coyote, or gray wolf). These novel wolf and coyote genome reference assemblies will be important resources for proper and accurate inference of Canis demography, taxonomic evaluation, and conservation genetics.

The role of structural variants in pest adaptation and genome evolution of the Colorado potato beetle, Leptinotarsa decemlineata (Say).

Structural variation has been associated with genetic diversity and adaptation. Despite these observations, it is not clear what their relative importance is for evolution, especially in rapidly adapting species. Here, we examine the significance of structural polymorphisms in pesticide resistance evolution of the agricultural super-pest, the Colorado potato beetle, Leptinotarsa decemlineata. By employing a parent offspring trio sequencing procedure, we develop highly contiguous reference genomes to characterize structural variation. These updated assemblies represent >100-fold improvement of contiguity and include derived pest and ancestral nonpest individuals. We identify >200,000 structural variations, which appear to be nonrandomly distributed across the genome as they co-occur with transposable elements and genes. Structural variations intersect with exons in a large proportion of gene annotations (~20%) that are associated with insecticide resistance (including cytochrome P450s), development, and transcription. To understand the role structural variations play in adaptation, we measure their allele frequencies among an additional 57 individuals using whole genome resequencing data, which represents pest and nonpest populations of North America. Incorporating multiple independent tests to detect the signature of natural selection using SNP data, we identify 14 genes that are probably under positive selection, include structural variations, and SNPs of elevated frequency within the pest lineages. Among these, three are associated with insecticide resistance based on previous research. One of these genes, CYP4g15, is coinduced during insecticide exposure with glycosyltransferase-13, which is a duplicated gene enclosed within a structural variant adjacent to the CYP4g15 genic region. These results demonstrate the significance of structural variations as a genomic feature to describe species history, genetic diversity, and adaptation.

Ultracontinuous Single Haplotype Genome Assemblies for the Domestic Cat (Felis catus) and Asian Leopard Cat (Prionailurus bengalensis).

In addition to including one of the most popular companion animals, species from the cat family Felidae serve as a powerful system for genetic analysis of inherited and infectious disease, as well as for the study of phenotypic evolution and speciation. Previous diploid-based genome assemblies for the domestic cat have served as the primary reference for genomic studies within the cat family. However, these versions suffered from poor resolution of complex and highly repetitive regions, with substantial amounts of unplaced sequence that is polymorphic or copy number variable. We sequenced the genome of a female F1 Bengal hybrid cat, the offspring of a domestic cat (Felis catus) x Asian leopard cat (Prionailurus bengalensis) cross, with PacBio long sequence reads and used Illumina sequence reads from the parents to phase >99.9% of the reads into the 2 species' haplotypes. De novo assembly of the phased reads produced highly continuous haploid genome assemblies for the domestic cat and Asian leopard cat, with contig N50 statistics exceeding 83 Mb for both genomes. Whole-genome alignments reveal the Felis and Prionailurus genomes are colinear, and the cytogenetic differences between the homologous F1 and E4 chromosomes represent a case of centromere repositioning in the absence of a chromosomal inversion. Both assemblies offer significant improvements over the previous domestic cat reference genome, with a 100% increase in contiguity and the capture of the vast majority of chromosome arms in 1 or 2 large contigs. We further demonstrated that comparably accurate F1 haplotype phasing can be achieved with members of the same species when one or both parents of the trio are not available. These novel genome resources will empower studies of feline precision medicine, adaptation, and speciation.

A Reference Genome Assembly of Simmental Cattle, Bos taurus taurus.

Genomics research has relied principally on the establishment and curation of a reference genome for the species. However, it is increasingly recognized that a single reference genome cannot fully describe the extent of genetic variation within many widely distributed species. Pangenome representations are based on high-quality genome assemblies of multiple individuals and intended to represent the broadest possible diversity within a species. A Bovine Pangenome Consortium (BPC) has recently been established to begin assembling genomes from more than 600 recognized breeds of cattle, together with other related species to provide information on ancestral alleles and haplotypes. Previously reported de novo genome assemblies for Angus, Brahman, Hereford, and Highland breeds of cattle are part of the initial BPC effort. The present report describes a complete single haplotype assembly at chromosome-scale for a fullblood Simmental cow from an F1 bison-cattle hybrid fetus by trio binning. Simmental cattle, also known as Fleckvieh due to their red and white spots, originated in central Europe in the 1830s as a triple-purpose breed selected for draught, meat, and dairy production. There are over 50 million Simmental cattle in the world, known today for their fast growth and beef yields. This assembly (ARS_Simm1.0) is similar in length to the other bovine assemblies at 2.86 Gb, with a scaffold N50 of 102 Mb (max scaffold 156.8 Mb) and meets or exceeds the continuity of the best Bos taurus reference assemblies to date.

A Reference Genome Assembly of American Bison, Bison bison bison.

Bison are an icon of the American West and an ecologically, commercially, and culturally important species. Despite numbering in the hundreds of thousands today, conservation concerns remain for the species, including the impact on genetic diversity of a severe bottleneck around the turn of the 20th century and genetic introgression from domestic cattle. Genetic diversity and admixture are best evaluated at genome-wide scale, for which a high-quality reference is necessary. Here, we use trio binning of long reads from a bison-Simmental cattle (Bos taurus taurus) male F1 hybrid to sequence and assemble the genome of the American plains bison (Bison bison bison). The male haplotype genome is chromosome-scale, with a total length of 2.65 Gb across 775 scaffolds (839 contigs) and a scaffold N50 of 87.8 Mb. Our bison genome is ~13× more contiguous overall and ~3400× more contiguous at the contig level than the current bison reference genome. The bison genome sequence presented here (ARS-UCSC_bison1.0) will enable new research into the evolutionary history of this iconic megafauna species and provide a new tool for the management of bison populations in federal and commercial herds.

Hybrid de novo and haplotype-resolved genome assembly of Vechur cattle - elucidating genetic variation.

Cattle contribute to the nutritional needs and economy of a place. The performance and fitness of cattle depend on the response and adaptation to local climatic conditions. Genomic and genetic studies are important for advancing cattle breeding, and availability of relevant reference genomes is essential. In the present study, the genome of a Vechur calf was sequenced on both short-read Illumina and long-read Nanopore sequencing platforms. The hybrid de novo assembly approach was deployed to obtain an average contig length of 1.97 Mbp and an N50 of 4.94 Mbp. By using a short-read genome sequence of the corresponding sire and dam, a haplotype-resolved genome was also assembled. In comparison to the taurine reference genome, we found 28,982 autosomal structural variants and 16,926,990 SNVs, with 883,544 SNVs homozygous in the trio samples. Many of these SNPs have been reported to be associated with various QTLs including growth, milk yield, and milk fat content, which are crucial determinants of cattle production. Furthermore, population genotype data analysis indicated that the present sample belongs to an Indian cattle breed forming a unique cluster of Bos indicus. Subsequent FST analysis revealed differentiation of the Vechur cattle genome at multiple loci, especially those regions related to whole body growth and cell division, especially IGF1, HMGA2, RRM2, and CD68 loci, suggesting a possible role of these genes in its small stature and better disease resistance capabilities in comparison with the local crossbreeds. This provides an opportunity to select and engineer cattle breeds optimized for local conditions.

Representing true plant genomes: haplotype-resolved hybrid pepper genome with trio-binning.

As sequencing costs decrease and availability of high fidelity long-read sequencing increases, generating experiment specific de novo genome assemblies becomes feasible. In many crop species, obtaining the genome of a hybrid or heterozygous individual is necessary for systems that do not tolerate inbreeding or for investigating important biological questions, such as hybrid vigor. However, most genome assembly methods that have been used in plants result in a merged single sequence representation that is not a true biologically accurate representation of either haplotype within a diploid individual. The resulting genome assembly is often fragmented and exhibits a mosaic of the two haplotypes, referred to as haplotype-switching. Important haplotype level information, such as causal mutations and structural variation is therefore lost causing difficulties in interpreting downstream analyses. To overcome this challenge, we have applied a method developed for animal genome assembly called trio-binning to an intra-specific hybrid of chili pepper (Capsicum annuum L. cv. HDA149 x Capsicum annuum L. cv. HDA330). We tested all currently available softwares for performing trio-binning, combined with multiple scaffolding technologies including Bionano to determine the optimal method of producing the best haplotype-resolved assembly. Ultimately, we produced highly contiguous biologically true haplotype-resolved genome assemblies for each parent, with scaffold N50s of 266.0 Mb and 281.3 Mb, with 99.6% and 99.8% positioned into chromosomes respectively. The assemblies captured 3.10 Gb and 3.12 Gb of the estimated 3.5 Gb chili pepper genome size. These assemblies represent the complete genome structure of the intraspecific hybrid, as well as the two parental genomes, and show measurable improvements over the currently available reference genomes. Our manuscript provides a valuable guide on how to apply trio-binning to other plant genomes.

Phased grapevine genome sequence of an Rpv12 carrier for biotechnological exploration of resistance to Plasmopara viticola.

The downy mildew disease caused by the oomycete Plasmopara viticola is a serious threat for grapevine and can cause enormous yield losses in viticulture. The quantitative trait locus Rpv12, mediating resistance against P. viticola, was originally found in Asian Vitis amurensis. This locus and its genes were analyzed here in detail. A haplotype-separated genome sequence of the diploid Rpv12-carrier Gf.99-03 was created and annotated. The defense response against P. viticola was investigated in an infection time-course RNA-seq experiment, revealing approximately 600 upregulated Vitis genes during host-pathogen interaction. The Rpv12 regions of the resistance and the sensitivity encoding Gf.99-03 haplotype were structurally and functionally compared with each other. Two different clusters of resistance-related genes were identified within the Rpv12 locus. One cluster carries a set of four differentially expressed genes with three ACCELERATED CELL DEATH 6-like genes. The other cluster carries a set of six resistance gene analogs related to qualitative pathogen resistance. The Rpv12 locus and its candidate genes for P. viticola resistance provide a precious genetic resource for P. viticola resistance breeding. Newly developed co-segregating simple sequence repeat markers in close proximity to the R-genes enable its improved applicability in marker-assisted grapevine breeding.

A new haplotype-resolved turkey genome to enable turkey genetics and genomics research.

BACKGROUND: The domesticated turkey (Meleagris gallopavo) is a species of significant agricultural importance and is the second largest contributor, behind broiler chickens, to world poultry meat production. The previous genome is of draft quality and partly based on the chicken (Gallus gallus) genome. A high-quality reference genome of M. gallopavo is essential for turkey genomics and genetics research and the breeding industry. RESULTS: By adopting the trio-binning approach, we were able to assemble a high-quality chromosome-level F1 assembly and 2 parental haplotype assemblies, leveraging long-read technologies and genome-wide chromatin interaction data (Hi-C). From a total of 40 chromosomes (2n = 80), we captured 35 chromosomes in a single scaffold, showing much improved genome completeness and continuity compared to the old assembly build. The 3 assemblies are of higher quality than the previous draft quality assembly and comparable to the chicken assemblies (GRCg7) shown by the largest contig N50 (26.6 Mb) and comparable BUSCO gene set completeness scores (96-97%). Comparative analyses confirm a previously identified large inversion of around 19 Mbp on the Z chromosome not found in other Galliformes. Structural variation between the parent haplotypes was identified, which poses potential new target genes for breeding. CONCLUSIONS: We contribute a new high-quality turkey genome at the chromosome level, benefiting turkey genetics and other avian genomics research as well as the turkey breeding industry.

Haplogenome assembly reveals structural variation in Eucalyptus interspecific hybrids.

BACKGROUND: De novo phased (haplo)genome assembly using long-read DNA sequencing data has improved the detection and characterization of structural variants (SVs) in plant and animal genomes. Able to span across haplotypes, long reads allow phased, haplogenome assembly in highly outbred organisms such as forest trees. Eucalyptus tree species and interspecific hybrids are the most widely planted hardwood trees with F1 hybrids of Eucalyptus grandis and E. urophylla forming the bulk of fast-growing pulpwood plantations in subtropical regions. The extent of structural variation and its effect on interspecific hybridization is unknown in these trees. As a first step towards elucidating the extent of structural variation between the genomes of E. grandis and E. urophylla, we sequenced and assembled the haplogenomes contained in an F1 hybrid of the two species. FINDINGS: Using Nanopore sequencing and a trio-binning approach, we assembled the separate haplogenomes (566.7 Mb and 544.5 Mb) to 98.0% BUSCO completion. High-density SNP genetic linkage maps of both parents allowed scaffolding of 88.0% of the haplogenome contigs into 11 pseudo-chromosomes (scaffold N50 of 43.8 Mb and 42.5 Mb for the E. grandis and E. urophylla haplogenomes, respectively). We identify 48,729 SVs between the two haplogenomes providing the first detailed insight into genome structural rearrangement in these species. The two haplogenomes have similar gene content, 35,572 and 33,915 functionally annotated genes, of which 34.7% are contained in genome rearrangements. CONCLUSIONS: Knowledge of SV and haplotype diversity in the two species will form the basis for understanding the genetic basis of hybrid superiority in these trees.

Construction and analysis of the chromosome-level haplotype-resolved genomes of two Crassostrea oyster congeners: Crassostrea angulata and Crassostrea gigas.

BACKGROUND: The Portuguese oyster Crassostrea angulata and the Pacific oyster C. gigas are two major Crassostrea species that are naturally distributed along the Northwest Pacific coast and possess great ecological and economic value. Here, we report the construction and comparative analysis of the chromosome-level haplotype-resolved genomes of the two oyster congeners. FINDINGS: Based on a trio-binning strategy, the PacBio high-fidelity and Illumina Hi-C reads of the offspring of the hybrid cross C. angulata (♂) × C. gigas (♀) were partitioned and independently assembled to construct two chromosome-level fully phased genomes. The assembly size (contig N50 size, BUSCO completeness) of the two genomes were 582.4 M (12.8 M, 99.1%) and 606.4 M (5.46 M, 98.9%) for C. angulata and C. gigas, respectively, ranking at the top of mollusk genomes with high contiguity and integrity. The general features of the two genomes were highly similar, and 15,475 highly conserved ortholog gene pairs shared identical gene structures and similar genomic locations. Highly similar sequences can be primarily identified in the coding regions, whereas most noncoding regions and introns of genes in the same ortholog group contain substantial small genomic and/or structural variations. Based on population resequencing analysis, a total of 2,756 species-specific single-nucleotide polymorphisms and 1,088 genes possibly under selection were identified. CONCLUSIONS: This is the first report of trio-binned fully phased chromosome-level genomes in marine invertebrates. The study provides fundamental resources for the research on mollusk genetics, comparative genomics, and molecular evolution.

A haploid diamondback moth (Plutella xylostella L.) genome assembly resolves 31 chromosomes and identifies a diamide resistance mutation.

The diamondback moth, Plutella xylostella (L.), is a highly mobile brassica crop pest with worldwide distribution and can rapidly evolve resistance to insecticides, including group 28 diamides. Reference genomes assembled using Illumina sequencing technology have provided valuable resources to advance our knowledge regarding the biology, origin and movement of diamondback moth, and more recently with its sister species, Plutella australiana. Here we apply a trio binning approach to sequence and annotate a chromosome level reference genome of P. xylostella using PacBio Sequel and Dovetail Hi-C sequencing technology and identify a point mutation that causes resistance to commercial diamides. A P. xylostella population collected from brassica crops in the Lockyer Valley, Australia (LV-R), was reselected for chlorantraniliprole resistance then a single male was crossed to a P. australiana female and a hybrid pupa sequenced. A chromosome level 328 Mb P. xylostella genome was assembled with 98.1% assigned to 30 autosomes and the Z chromosome. The genome was highly complete with 98.4% of BUSCO Insecta genes identified and RNAseq informed protein prediction annotated 19,002 coding genes. The LV-R strain survived recommended field application doses of chlorantraniliprole, flubendiamide and cyclaniliprole. Some hybrids also survived these doses, indicating significant departure from recessivity, which has not been previously documented for diamides. Diamide chemicals modulate insect Ryanodine Receptors (RyR), disrupting calcium homeostasis, and we identified an amino acid substitution (I4790K) recently reported to cause diamide resistance in a strain from Japan. This chromosome level assembly provides a new resource for insect comparative genomics and highlights the emergence of diamide resistance in Australia. Resistance management plans need to account for the fact that resistance is not completely recessive.

Merqury: reference-free quality, completeness, and phasing assessment for genome assemblies.

Recent long-read assemblies often exceed the quality and completeness of available reference genomes, making validation challenging. Here we present Merqury, a novel tool for reference-free assembly evaluation based on efficient k-mer set operations. By comparing k-mers in a de novo assembly to those found in unassembled high-accuracy reads, Merqury estimates base-level accuracy and completeness. For trios, Merqury can also evaluate haplotype-specific accuracy, completeness, phase block continuity, and switch errors. Multiple visualizations, such as k-mer spectrum plots, can be generated for evaluation. We demonstrate on both human and plant genomes that Merqury is a fast and robust method for assembly validation.

Draft Genomes of Amaranthus tuberculatus, Amaranthus hybridus, and Amaranthus palmeri.

Amaranthus tuberculatus, Amaranthus hybridus, and Amaranthus palmeri are agronomically important weed species. Here, we present the most contiguous draft assemblies of these three species to date. We utilized a combination of Pacific Biosciences long-read sequencing and chromatin contact mapping information to assemble and order sequences of A. palmeri to near-chromosome-level resolution, with scaffold N50 of 20.1 Mb. To resolve the issues of heterozygosity and coassembly of alleles in diploid species, we adapted the trio binning approach to produce haplotype assemblies of A. tuberculatus and A. hybridus. This approach resulted in an improved assembly of A. tuberculatus, and the first genome assembly for A. hybridus, with contig N50s of 2.58 and 2.26 Mb, respectively. Species-specific transcriptomes and information from related species were used to predict transcripts within each assembly. Syntenic comparisons of these species and Amaranthus hypochondriacus identified sites of genomic rearrangement, including duplication and translocation, whereas genetic map construction within A. tuberculatus highlighted the need for further ordering of the A. hybridus and A. tuberculatus contigs. These multiple reference genomes will accelerate genomic studies in these species to further our understanding of weedy evolution within Amaranthus.

A haplotype-resolved, de novo genome assembly for the wood tiger moth (Arctia plantaginis) through trio binning.

BACKGROUND: Diploid genome assembly is typically impeded by heterozygosity because it introduces errors when haplotypes are collapsed into a consensus sequence. Trio binning offers an innovative solution that exploits heterozygosity for assembly. Short, parental reads are used to assign parental origin to long reads from their F1 offspring before assembly, enabling complete haplotype resolution. Trio binning could therefore provide an effective strategy for assembling highly heterozygous genomes, which are traditionally problematic, such as insect genomes. This includes the wood tiger moth (Arctia plantaginis), which is an evolutionary study system for warning colour polymorphism. FINDINGS: We produced a high-quality, haplotype-resolved assembly for Arctia plantaginis through trio binning. We sequenced a same-species family (F1 heterozygosity ∼1.9%) and used parental Illumina reads to bin 99.98% of offspring Pacific Biosciences reads by parental origin, before assembling each haplotype separately and scaffolding with 10X linked reads. Both assemblies are contiguous (mean scaffold N50: 8.2 Mb) and complete (mean BUSCO completeness: 97.3%), with annotations and 31 chromosomes identified through karyotyping. We used the assembly to analyse genome-wide population structure and relationships between 40 wild resequenced individuals from 5 populations across Europe, revealing the Georgian population as the most genetically differentiated with the lowest genetic diversity. CONCLUSIONS: We present the first invertebrate genome to be assembled via trio binning. This assembly is one of the highest quality genomes available for Lepidoptera, supporting trio binning as a potent strategy for assembling heterozygous genomes. Using our assembly, we provide genomic insights into the geographic population structure of A. plantaginis.