Optimization of the “in-silico” mate-pair method improved contiguity and accuracy of genome assembly

A combination of next-generation sequencing technologies and mate-pair libraries of large insert sizes is used as a standard method to generate genome assemblies with high contiguity. The third-generation sequencing techniques also are used to improve the quality of assembled genomes. However, both mate-pair libraries and the third-generation libraries require high-molecular-weight DNA, making the use of these libraries inappropriate for samples with only degraded DNA. An in silico method that generates mate-pair libraries using a reference genome was devised for the task of assembling target genomes. Although the contiguity and completeness of assembled genomes were significantly improved by this method, a high level of errors manifested in the assembly, further to which the methods for using reference genomes were not optimized. Here, we tested different strategies for using reference genomes to generate in silico mate-pairs. The results showed that using a closely related reference genome from the same genus was more effective than using divergent references. Conservation of in silico mate-pairs by comparing two references and using those to guide genome assembly reduced the number of misassemblies (18.6% – 46.1%) and increased the contiguity of assembled genomes (9.7% – 70.7%), while maintaining gene completeness at a level that was either similar or marginally lower than that obtained via the current method. Finally, we compared the optimized method with another reference-guided assembler, RaGOO. We found that RaGOO produced longer scaffolds (17.8 Mbp vs 3.0 Mbp), but resulted in a much higher misassembly rate (85.68%) than our optimized in silico mate-pair method.

Recombination Marks the Evolutionary Dynamics of a Recently Endogenized Retrovirus

All vertebrate genomes have been colonized by retroviruses along their evolutionary trajectory. Although endogenous retroviruses (ERVs) can contribute important physiological functions to contemporary hosts, such benefits are attributed to long-term coevolution of ERV and host because germline infections are rare and expansion is slow, and because the host effectively silences them. The genomes of several outbred species including mule deer (Odocoileus hemionus) are currently being colonized by ERVs, which provides an opportunity to study ERV dynamics at a time when few are fixed. We previously established the locus-specific distribution of cervid ERV (CrERV) in populations of mule deer. In this study, we determine the molecular evolutionary processes acting on CrERV at each locus in the context of phylogenetic origin, genome location, and population prevalence. A mule deer genome was de novo assembled from short- and long-insert mate pair reads and CrERV sequence generated at each locus. We report that CrERV composition and diversity have recently measurably increased by horizontal acquisition of a new retrovirus lineage. This new lineage has further expanded CrERV burden and CrERV genomic diversity by activating and recombining with existing CrERV. Resulting interlineage recombinants then endogenize and subsequently expand. CrERV loci are significantly closer to genes than expected if integration were random and gene proximity might explain the recent expansion of one recombinant CrERV lineage. Thus, in mule deer, retroviral colonization is a dynamic period in the molecular evolution of CrERV that also provides a burst of genomic diversity to the host population.

A de novo 10q11.23q22.1 deletion detected by whole genome mate-pair sequencing: a case report

Background Interstitial deletions of chromosome band 10q11-q22 was a genomic disorder distinguished by developmental delay, congenital cleft palate and muscular hypotonia. The phenotypes involved were heterogeneous, hinge on the variable breakpoints and size. Case presentation Here, we presented a patient with soft palate cleft, growth and development delay. The patient was a 2 years and 5 months girl who was not able to walk unless using a children’s crutches to support herself. Whole-exome sequencing (WES) and whole-genome mate-pair sequencing (WGMS) were both performed by next generation sequencing (NGS). A 20.76 Mb deletion at 10q11.23q22.1 (seq[GRCh37/hg19]del(10)(50,319,387-71,083,899) × 1) was revealed by the WGMS, which was verified as de novo by quantitative polymerase chain reaction (QPCR). Conclusion Children with 10q11-q22 deletions greater than 20 MB have never been reported before, and we are the first to report and provide a detailed clinical phenotype, which brings further knowledge of 10q11-q22 deletions.

A chromosome-scale strawberry genome assembly of a Japanese variety, Reikou

Cultivated strawberry (Fragaria × ananassa) is an octoploid species (2n = 8x= 56) that is widely consumed around the world as both fresh and processed fruit. In this study, we report a chromosome-scale strawberry genome assembly of a Japanese variety, Reikou. The Illumina short reads derived from paired-end, mate-pair, and 10X Genomics libraries were assembled using Denovo MAGIC 3.0. The generated phased scaffolds consisted of 32,715 sequences with a total length of 1.4 Gb and an N50 length of 3.9 Mb. A total of 63 pseudomolecules including chr0 were created by aligning the scaffolds onto the Reikou S1 linkage maps with the IStraw90 Axiom SNP array and ddRAD-Seq. Meanwhile, genomes of diploid Fragaria species were resequenced and compared with the most similar chromosome-scale scaffolds to investigate the possible progenitor of each subgenome. Clustering analysis suggested that the most likely progenitors were F. vesca and F. iinumae. The phased pseudomolecules were assigned the scaffolds names with Av, Bi, and X, representing sequence similarity with F. vesca (Av), F. iinumae (Bi), and others (X), respectively. The result of a comparison with the Camerosa genome suggested the possibility of subgenome structure differences between the two varieties.