Benchmarking datasets for assembly-based variant calling using high-fidelity long reads

Highly accurate Korean draft genomes reveal structural variation highlighting human telomere evolution

Given the presence of highly repetitive genomic regions such as subtelomeric regions, understanding human genomic evolution remains challenging. Recently, long-read sequencing technology has facilitated the identification of complex genetic variants, including structural variants (SVs), at the single-nucleotide level. Here, we resolved SVs and their underlying DNA damage-repair mechanisms in subtelomeric regions, which are among the most uncharted genomic regions. We generated ∼20 × high-fidelity long-read sequencing data from three Korean individuals and their partially phased high-quality de novo genome assemblies (contig N50: 6.3-58.2 Mb). We identified 131 138 deletion and 121 461 insertion SVs, 41.6% of which were prevalent in the East Asian population. The commonality of the SVs identified among the Korean population was examined by short-read sequencing data from 103 Korean individuals, providing the first comprehensive SV set representing the population based on the long-read assemblies. Manual investigation of 19 large subtelomeric SVs (≥5 kb) and their associated repair signatures revealed the potential repair mechanisms leading to the formation of these SVs. Our study provides mechanistic insight into human telomere evolution and can facilitate our understanding of human SV formation.

HiFi long-read genomes for difficult-to-detect, clinically relevant variants

Clinical short-read exome and genome sequencing approaches have positively impacted diagnostic testing for rare diseases. Yet, technical limitations associated with short reads challenge their use for the detection of disease-associated variation in complex regions of the genome. Long-read sequencing (LRS) technologies may overcome these challenges, potentially qualifying as a first-tier test for all rare diseases. To test this hypothesis, we performed LRS (30× high-fidelity [HiFi] genomes) for 100 samples with 145 known clinically relevant germline variants that are challenging to detect using short-read sequencing and necessitate a broad range of complementary test modalities in diagnostic laboratories. We show that relevant variant callers readily re-identified the majority of variants (120/145, 83%), including ∼90% of structural variants, SNVs/insertions or deletions (indels) in homologous sequences, and expansions of short tandem repeats. Another 10% (n = 14) was visually apparent in the data but not automatically detected. Our analyses also identified systematic challenges for the remaining 7% (n = 11) of variants, such as the detection of AG-rich repeat expansions. Titration analysis showed that 90% of all automatically called variants could also be identified using 15-fold coverage. Long-read genomes thus identified 93% of challenging pathogenic variants from our dataset. Even with reduced coverage, the vast majority of variants remained detectable, possibly enhancing cost-effective diagnostic implementation. Most importantly, we show the potential to use a single technology to accurately identify all types of clinically relevant variants.

Units containing telomeric repeats are prevalent in subtelomeric regions of a Mesorhabditis isolate collected from the Republic of Korea

BACKGROUND

Mesorhabditis is known for its somatic genome being only a small portion of the germline genome due to programmed DNA elimination. This phenotype may be associated with the maintenance of telomeres at the ends of fragmented somatic chromosomes.

OBJECTIVE

To comprehensively investigate the telomeric regions of Mesorhabditis nematodes at the sequence level, we endeavored to collect a Mesorhabditis nematode in the Republic of Korea and acquire its highly contiguous genome sequences.

METHODS

We isolated a Mesorhabditis nematode and assembled its 108-Mb draft genome using both 6.3 Gb (53 ×) of short-read and 3.0 Gb (25 × , N50 = 5.7 kb) of nanopore-based long-read sequencing data. Our genome assembly exhibits comparable quality to the public genome of Mesorhabditis belari in terms of contiguity and evolutionary conserved genes.

RESULTS

Unexpectedly, our Mesorhabditis genome has many more interstitial telomeric sequences (ITSs), specifically subtelomeric ones, compared to the genomes of Caenorhabditis elegans and M. belari. Moreover, several subtelomeric sequences containing ITSs had 4-26 homologous sequences, implying they are highly repetitive. Based on this highly repetitive nature, we hypothesize that subtelomeric ITSs might have accumulated through the action of transposable elements containing ITSs.

CONCLUSIONS

It still remains elusive whether these ITS-containing units are associated with programmed DNA elimination, but they may facilitate new telomere formation after DNA elimination. Our genomic resources for Mesorhabditis can aid in understanding how its distinct phenotypes have evolved.

Fully phased genome assemblies and graph-based genetic variants of the olive flounder, Paralichthys olivaceus

The olive flounder, Paralichthys olivaceus, also known as the Korean halibut, is an economically important flatfish in East Asian countries. Here, we provided four fully phased genome assemblies of two different olive flounder individuals using high-fidelity long-read sequencing and their parental short-read sequencing data. We obtained 42-44 Gb of ~15-kb and ~Q30 high-fidelity long reads, and their assembly quality values were ~53. We annotated ~30 K genes, ~170-Mb repetitive sequences, and ~3 M 5-methylcytosine positions for each genome assembly, and established a graph-based draft pan-genome of the olive flounder. We identified 5 M single-nucleotide variants and 100 K structural variants with their genotype information, where ~13% of the variants were possibly fixed in the two Korean individuals. Based on our chromosome-level genome assembly, we also explored chromosome evolution in the Pleuronectiformes family, as reported earlier. Our high-quality genomic resources will contribute to future genomic selection for accelerating the breeding process of the olive flounder.

Identification and characterization of structural variants related to meat quality in pigs using chromosome-level genome assemblies

BACKGROUND

Many studies have been performed to identify various genomic loci and genes associated with the meat quality in pigs. However, the full genetic architecture of the trait still remains unclear in part because of the lack of accurate identification of related structural variations (SVs) which resulted from the shortage of target breeds, the limitations of sequencing data, and the incompleteness of genome assemblies. The recent generation of a new pig breed with superior meat quality, called Nanchukmacdon, and its chromosome-level genome assembly (the NCMD assembly) has provided new opportunities.

RESULTS

By applying assembly-based SV calling approaches to various genome assemblies of pigs including Nanchukmacdon, the impact of SVs on meat quality was investigated. Especially, by checking the commonality of SVs with other pig breeds, a total of 13,819 Nanchukmacdon-specific SVs (NSVs) were identified, which have a potential effect on the unique meat quality of Nanchukmacdon. The regulatory potentials of NSVs for the expression of nearby genes were further examined using transcriptome- and epigenome-based analyses in different tissues.

CONCLUSIONS

Whole-genome comparisons based on chromosome-level genome assemblies have led to the discovery of SVs affecting meat quality in pigs, and their regulatory potentials were analyzed. The identified NSVs will provide new insights regarding genetic architectures underlying the meat quality in pigs. Finally, this study confirms the utility of chromosome-level genome assemblies and multi-omics analysis to enhance the understanding of unique phenotypes.

Ribotin: automated assembly and phasing of rDNA morphs

MOTIVATION

The ribosomal DNA (rDNA) arrays are highly repetitive and homogenous regions which exist in all life. Due to their repetitiveness, current assembly methods do not fully assemble the rDNA arrays in humans and many other eukaryotes, and so variation within the rDNA arrays cannot be effectively studied.

RESULTS

Here, we present the tool ribotin to assemble full length rDNA copies, or morphs. Ribotin uses a combination of highly accurate long reads and extremely long nanopore reads to resolve the variation between rDNA morphs. We show that ribotin successfully recovers the most abundant morphs in human and nonhuman genomes. We also find that genome wide consensus sequences of the rDNA arrays frequently produce a mosaic sequence that does not exist in the genome.

AVAILABILITY AND IMPLEMENTATION

Ribotin is available on https://github.com/maickrau/ribotin and as a package on bioconda.

A reference genome for ecological restoration of the sunflower sea star, Pycnopodia helianthoides

Wildlife diseases, such as the sea star wasting (SSW) epizootic that outbroke in the mid-2010s, appear to be associated with acute and/or chronic abiotic environmental change; dissociating the effects of different drivers can be difficult. The sunflower sea star, Pycnopodia helianthoides, was the species most severely impacted during the SSW outbreak, which overlapped with periods of anomalous atmospheric and oceanographic conditions, and there is not yet a consensus on the cause(s). Genomic data may reveal underlying molecular signatures that implicate a subset of factors and, thus, clarify past events while also setting the scene for effective restoration efforts. To advance this goal, we used Pacific Biosciences HiFi long sequencing reads and Dovetail Omni-C proximity reads to generate a highly contiguous genome assembly that was then annotated using RNA-seq-informed gene prediction. The genome assembly is 484 Mb long, with contig N50 of 1.9 Mb, scaffold N50 of 21.8 Mb, BUSCO completeness score of 96.1%, and 22 major scaffolds consistent with prior evidence that sea star genomes comprise 22 autosomes. These statistics generally fall between those of other recently assembled chromosome-scale assemblies for two species in the distantly related asteroid genus Pisaster. These novel genomic resources for P. helianthoides will underwrite population genomic, comparative genomic, and phylogenomic analyses-as well as their integration across scales-of SSW and environmental stressors.

Genomic Analysis of Antifungal Drug Resistance Induced in Trichophyton rubrum After Prolonged Culture with Terbinafine

In this study, we induced terbinafine (TRF) resistance in a T. rubrum strain in vitro for 18 months then compared the genomes of the TRF-resistant (N42-3) and TRF-susceptible wild-type (N42WT) strains to identify mutations. In the SQLE gene, N42WT had no mutation while N42-3 had a F397L mutation. We sequenced approximately 22.53 Mb of the genomes of the N43WT and N42-3 strains. Other than the F397L mutation in SQLE, there were three other genetic mutations in three different genes that were found in N42-3, but not in N43WT; however, these three mutations were not detected in other TRF-resistant T. rubrum strains. From this genome sequencing analysis, the only variation that was confirmed to be associated with drug resistance in the genome of the TRF-resistant T. rubrum was a hotspot mutation in SQLE.