Mate-pair Libraries Research Articles

Optimization of the "in-silico" mate-pair method improves contiguity and accuracy of genome assembly.

A combination of short-insert paired-ended 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, was 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 developed a pipeline of the optimized in silico method and compared it with another reference-guided assembler, RagTag. We found that RagTag 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. This optimized in silico pipeline developed in this study should facilitate further studies on genomics, population genetics, and conservation of endangered species.

Identification and analysis of deletion breakpoints in four Mohr-Tranebjærg syndrome (MTS) patients

Mohr-Tranebjærg syndrome is an X-linked syndrome characterized by sensorineural hearing impairment in childhood, followed by progressive neurodegeneration leading to a broad phenotypic spectrum. Genetically MTS is caused by pathogenic variants in the TIMM8A gene, including gene deletions and larger contiguous gene deletions. Some of the latter involve the neighboring gene BTK, resulting in agammaglobulinemia. By next‐generation mate‐pair sequencing we have mapped the chromosomal deletion breakpoints of one MTS case and three XLA-MTS cases and used breakpoint-spanning PCR to fine map the breakpoints by Sanger sequencing. Two of the XLA-MTS cases presented with large deletions (63.5 and 27.2 kb), and the junctional regions were characterized by long stretches of microhomology, indicating that the events have emerged through homologous recombination. Conversely, the MTS case exhibited a small 2 bp region of microhomology, and the regions were not characterized by extensive microhomology. The third XLA-MTS case had a more complex breakpoint, including a 59 bp inverted insertion, thus at least four breakpoints were involved in this event. In conclusion, mate-pair library generation combined with next-generation sequencing is an efficient method for breakpoint identification, also in regions characterized by repetitive elements.

Evaluating Illumina-, Nanopore-, and PacBio-based genome assembly strategies with the bald notothen, Trematomus borchgrevinki.

For any genome-based research, a robust genome assembly is required. De novo assembly strategies have evolved with changes in DNA sequencing technologies and have been through at least 3 phases: (1) short-read only, (2) short- and long-read hybrid, and (3) long-read only assemblies. Each of the phases has its own error model. We hypothesized that hidden short-read scaffolding errors and erroneous long-read contigs degrade the quality of short- and long-read hybrid assemblies. We assembled the genome of Trematomus borchgrevinki from data generated during each of the 3 phases and assessed the quality problems we encountered. We developed strategies such as k-mer-assembled region replacement, parameter optimization, and long-read sampling to address the error models. We demonstrated that a k-mer-based strategy improved short-read assemblies as measured by Benchmarking Universal Single-Copy Ortholog while mate-pair libraries introduced hidden scaffolding errors and perturbed Benchmarking Universal Single-Copy Ortholog scores. Furthermore, we found that although hybrid assemblies can generate higher contiguity they tend to suffer from lower quality. In addition, we found long-read-only assemblies can be optimized for contiguity by subsampling length-restricted raw reads. Our results indicate that long-read contig assembly is the current best choice and that assemblies from phase I and phase II were of lower quality.

Genome sequencing of the NIES Cyanobacteria collection with a focus on the heterocyst-forming clade.

Cyanobacteria are a diverse group of Gram-negative prokaryotes that perform oxygenic photosynthesis. Cyanobacteria have been used for research on photosynthesis and have attracted attention as a platform for biomaterial/biofuel production. Cyanobacteria are also present in almost all habitats on Earth and have extensive impacts on global ecosystems. Given their biological, economical, and ecological importance, the number of high-quality genome sequences for Cyanobacteria strains is limited. Here, we performed genome sequencing of Cyanobacteria strains in the National Institute for Environmental Studies microbial culture collection in Japan. We sequenced 28 strains that can form a heterocyst, a morphologically distinct cell that is specialized for fixing nitrogen, and 3 non-heterocystous strains. Using Illumina sequencing of paired-end and mate-pair libraries with in silico finishing, we constructed highly contiguous assemblies. We determined the phylogenetic relationship of the sequenced genome assemblies and found potential difficulties in the classification of certain heterocystous clades based on morphological observation. We also revealed a bias on the sequenced strains by the phylogenetic analysis of the 16S rRNA gene including unsequenced strains. Genome sequencing of Cyanobacteria strains deposited in worldwide culture collections will contribute to understanding the enormous genetic and phenotypic diversity within the phylum Cyanobacteria.

Genome report: a draft genome of Alliaria petiolata (garlic mustard) as a model system for invasion genetics.

The emerging field of invasion genetics examines the genetic causes and consequences of biological invasions, but few study systems are available that integrate deep ecological knowledge with genomic tools. Here, we report on the de novo assembly and annotation of a genome for the biennial herb Alliaria petiolata (M. Bieb.) Cavara and Grande (Brassicaceae), which is widespread in Eurasia and invasive across much of temperate North America. Our goal was to sequence and annotate a genome to complement resources available from hundreds of published ecological studies, a global field survey, and hundreds of genetic lines maintained in Germany and Canada. We sequenced a genotype (EFCC3-3-20) collected from the native range near Venice, Italy, and sequenced paired-end and mate pair libraries at ∼70 × coverage. A de novo assembly resulted in a highly continuous draft genome (N50 = 121 Mb; L50 = 2) with 99.7% of the 1.1 Gb genome mapping to scaffolds of at least 50 Kb in length. A total of 64,770 predicted genes in the annotated genome include 99% of plant BUSCO genes and 98% of transcriptome reads. Consistent with previous reports of (auto)hexaploidy in western Europe, we found that almost one-third of BUSCO genes (390/1440) mapped to two or more scaffolds despite <2% genome-wide average heterozygosity. The continuity and gene space quality of our draft assembly will enable molecular and functional genomic studies of A. petiolata to address questions relevant to invasion genetics and conservation strategies.

Whole-Genome Mate Pair Sequencing Transforms Cytogenetic Analysis of Hematologic Malignancies: Development of a Targeted Panel to Detect Diagnostic/Prognostic Chromosomal Rearrangements and Copy Number Changes in AML

Acute myeloid leukemia (AML) is the most common form of acute leukemia affecting both children and adults. Once recognized as a single disease, the World Health Organization now groups AML into ten subtypes based on genetic abnormality. This is based on differences among recurrent genetic abnormalities with respect to response to therapy, relapse risk and overall survival. Currently most testing of AML patients occurs by karyotype analysis, FISH or RT-PCR. Here we describe a next generation sequencing approach to better classify and accurately detect abnormalities associated with AML.For the validation of whole genome mate pair sequencing (MPseq) targeted panel to test patients with AML, twenty-five abnormal and ten normal samples extracted from each sample type (peripheral blood (PB) and bone marrow (BM)) were run with the MPseq assay. Samples were chosen based on reported abnormalities previously tested by conventional karyotyping and/or FISH. The mate pair panel is designed to detect rearrangements involving genes ABL1, BCR, CBFB, CREBBP, DEK, KAT6A(MYST3), MECOM(EVI1), MLF1, MLL(KMT2A), MYH11, NPM1, NUP214(CAN), NUP98, PML, RARA, RPN1, RUNX1(AML1), RUNX1T1(ETO) and copy number changes on chromosomes 5, 7, 8, 13, 17, and 20. DNA was processed using the Illumina Nextera Mate Pair library preparation kit, multiplexed at 2 samples/lane and sequenced on the Illumina HiSeq. Data was aligned to the reference genome (GRCh38) using BIMA and abnormalities were identified using SVAtools, both in-house developed tools. A targeted analysis approach was utilized using a well-defined list of recurrent chromosomal abnormalities associated with known diagnostic and prognostic significance as well as targeted therapies in AML. The karyotype, FISH and MPseq results were compared to determine concordance.To test the limit of detection based on percent tumor, samples previously tested with FISH were diluted. Tumor percentage was determined by the abnormal FISH result and was diluted with a normal sample to mimic various percentages (100%, 50%, 25%, 10%, 5%, 2%). For both PB and BM specimens, the rearrangements were detected down to the 10% dilution and the copy number changes were detected down to 25%. Based on these cutoffs, all 50 abnormal samples (25 PB, 25 BM) used for the validation were 100% concordant with abnormalities targeted by the panel that were previously detected by FISH and/or karyotype. In addition, MPseq identified additional abnormalities not detected by karyotype or FISH and further characterized several complex abnormalities. One example includes a case with MLLT10/MLL(KMT2A) fusion detected by MPseq which was missed by the AML FISH panel because the rearrangement involved an insertional translocation. Another example is the identification of BCR/FGFR1 fusion in a case when AML FISH identified three copies of BCR.In conclusion, we validated a whole genome mate pair sequencing targeted panel to detect diagnostic/prognostic chromosomal rearrangements and copy number changes in patients with AML. Due to the limitations in detecting very low level abnormalities, MPseq would not be recommended for follow-up post therapy or minimal residual disease testing unless the coverage is increased with deeper sequencing. However, our validation study demonstrate the clinical utility of MPseq as a potential replacement assay to conventional FISH on diagnostic AML samples and highlight the increased diagnostic yield. Thus, MPseq represents a next generation sequencing technology that has the potential to revolutionize the diagnosis of hematologic malignancies and provide an opportunity to advance precision medicine. DisclosuresNo relevant conflicts of interest to declare.

Chromosome-level genome assembly of the blue crab, Callinectes sapidus.

The blue crab, Callinectes sapidus (Rathbun, 1896) is an economically, culturally, and ecologically important species found across the temperate and tropical North and South American Atlantic coast. A reference genome will enable research for this high-value species. Initial assembly combined 200× coverage Illumina paired-end reads, a 60× 8 kb mate-paired library, and 50× PacBio data using the MaSuRCA assembler resulting in a 985 Mb assembly with a scaffold N50 of 153 kb. Dovetail Chicago and HiC sequencing with the 3d DNA assembler and Juicebox assembly tools were then used for chromosome scaffolding. The 50 largest scaffolds span 810 Mb are 1.5–37 Mb long and have a repeat content of 36%. The 190 Mb unplaced sequence is in 3921 sequences over 10 kb with a repeat content of 68%. The final assembly N50 is 18.9 Mb for scaffolds and 9317 bases for contigs. Of arthropod BUSCO, ∼88% (888/1013) were complete and single copies. Using 309 million RNAseq read pairs from 12 different tissues and developmental stages, 25,249 protein-coding genes were predicted. Between C. sapidus and Portunus trituberculatus genomes, 41 of 50 large scaffolds had high nucleotide identity and protein-coding synteny, but 9 scaffolds in both assemblies were not clear matches. The protein-coding genes included 9423 one-to-one putative orthologs, of which 7165 were syntenic between the two crab species. Overall, the two crab genome assemblies show strong similarities at the nucleotide, protein, and chromosome level and verify the blue crab genome as an excellent reference for this important seafood species.

Improved chromosome-level genome assembly and annotation of the seagrass, Zostera marina (eelgrass)

Background: Seagrasses (Alismatales) are the only fully marine angiosperms. Zostera marina (eelgrass) plays a crucial role in the functioning of coastal marine ecosystems and global carbon sequestration. It is the most widely studied seagrass and has become a marine model system for exploring adaptation under rapid climate change. The original draft genome (v.1.0) of the seagrass Z. marina (L.) was based on a combination of Illumina mate-pair libraries and fosmid-ends. A total of 25.55 Gb of Illumina and 0.14 Gb of Sanger sequence was obtained representing 47.7× genomic coverage. The assembly resulted in ~2000 unordered scaffolds (L50 of 486 Kb), a final genome assembly size of 203MB, 20,450 protein coding genes and 63% TE content. Here, we present an upgraded chromosome-scale genome assembly and compare v.1.0 and the new v.3.1, reconfirming previous results from Olsen et al. (2016), as well as pointing out new findings. Methods: The same high molecular weight DNA used in the original sequencing of the Finnish clone was used. A high-quality reference genome was assembled with the MECAT assembly pipeline combining PacBio long-read sequencing and Hi-C scaffolding. Results: In total, 75.97 Gb PacBio data was produced. The final assembly comprises six pseudo-chromosomes and 304 unanchored scaffolds with a total length of 260.5Mb and an N50 of 34.6 MB, showing high contiguity and few gaps (~0.5%). 21,483 protein-encoding genes are annotated in this assembly, of which 20,665 (96.2%) obtained at least one functional assignment based on similarity to known proteins. Conclusions: As an important marine angiosperm, the improved Z. marina genome assembly will further assist evolutionary, ecological, and comparative genomics at the chromosome level. The new genome assembly will further our understanding into the structural and physiological adaptations from land to marine life.

Improved chromosome-level genome assembly and annotation of the seagrass, Zostera marina (eelgrass).

Background:Seagrasses(Alismatales) are the only fully marine angiosperms. Zostera marina(eelgrass) plays a crucial role in the functioning of coastal marine ecosystems and global carbon sequestration. It is the most widely studied seagrass and has become a marine model system forexploringadaptation under rapid climate change.The original draft genome(v.1.0)of the seagrass Z. marina(L.)was based on a combination of Illumina mate-pair libraries andfosmid-ends. A total of 25.55 Gb of Illumina and 0.14 Gb of Sanger sequence was obtained representing 47.7×genomic coverage. The assembly resulted in~2000 unordered scaffolds (L50 of 486Kb), a final genome assembly size of 203MB, 20,450 protein coding genes and 63% TE content. Here,we present an upgraded chromosome-scale genome assemblyand comparev.1.0and the newv.3.1, reconfirming previous results from Olsen et al. (2016), as well as pointing out new findings. Methods:The same high molecular weight DNAused in the original sequencingof the Finnish clonewas used.Ahigh-qualityreferencegenome was assembledwiththeMECATassemblypipelinecombiningPacBiolong-readsequencingand Hi-C scaffolding. Results:In total,75.97 GbPacBiodata wasproduced.The final assembly comprisessixpseudo-chromosomesand304 unanchored scaffolds with a total length of 260.5Mb andanN50 of 34.6 MB, showing high contiguity and few gaps (~0.5%). 21,483 protein-encoding genes are annotated in this assembly, of which 20,665 (96.2%)obtained at least one functional assignment based on similarity to known proteins. Conclusions:Asanimportant marine angiosperm,the improved Z. marinagenomeassemblywillfurtherassistevolutionary, ecological,and comparative genomicsat the chromosome level.Thenew genome assemblywillfurther our understandinginto the structural and physiological adaptations from land to marinelife.

Refinement of Draft Genome Assemblies of Pigeonpea (Cajanus cajan).

Genome assembly of short reads from large plant genomes remains a challenge in computational biology despite major developments in next generation sequencing. Of late several draft assemblies have been reported in sequenced plant genomes. The reported draft genome assemblies of Cajanus cajan have different levels of genome completeness, a large number of repeats, gaps, and segmental duplications. Draft assemblies with portions of genome missing are shorter than the referenced original genome. These assemblies come with low map accuracy affecting further functional annotation and the prediction of gene components as desired by crop researchers. Genome coverage, i.e., the number of sequenced raw reads mapped onto a certain location of the genome is an important quality indicator of completeness and assembly quality in draft assemblies. The present work aimed to improve the coverage in reported de novo sequenced draft genomes (GCA_000340665.1 and GCA_000230855.2) of pigeonpea, a legume widely cultivated in India. The two recently sequenced assemblies, A1 and A2 comprised 72% and 75% of the estimated coverage of the genome, respectively. We employed an assembly reconciliation approach to compare the draft assemblies and merge them, filling the gaps by employing an algorithm size sorting mate-pair library to generate a high quality and near complete assembly with enhanced contiguity. The majority of gaps present within scaffolds were filled with right-sized mate-pair reads. The improved assembly reduced the number of gaps than those reported in draft assemblies resulting in an improved genome coverage of 82.4%. Map accuracy of the improved assembly was evaluated using various quality metrics and for the presence of specific trait-related functional genes. Employed pair-end and mate-pair local libraries helped us to reduce gaps, repeats, and other sequence errors resulting in lengthier scaffolds compared to the two draft assemblies. We reported the prediction of putative host resistance genes against Fusarium wilt disease by their performance and evaluated them both in wet laboratory and field phenotypic conditions.

Genome/transcriptome analysis of the chigger mite Leptotrombidium pallidum, a major vector for scrub typhus, with a special focus on genes more abundantly expressed in larval stage

• The genome size of L. pallidum was determined to be 193 Mb and 15,842 genes were predicted. • A total of 10,885 protein-coding genes were constitutively expressed in all four stages. • A total of 4,078 genes were differentially expressed in all four stages. • Whole genome amplification is a reliable option for genome sequencing of a single tiny individual. Leptotrombidium pallidum is the major vector mite for Orientia tsutsugamushi , the causative agent of scrub typhus, in Asian countries, including Korea. Despite its medical importance, L. pallidum has little genetic information available to date. To analyze the L. pallidum genome, we extracted genomic DNA (gDNA) from a single female of a 7-generation inbred L. pallidum colony and amplified the gDNA by whole genome amplification (WGA). The resulting amplified gDNA was used to construct paired-end and mate-pair libraries that were sequenced using Illumina platforms (HiSeq2000 and MiSeq). More than 45 Gb of sequence reads from both paired-end and mate-pair libraries of the WGA gDNA were trimmed and then de novo assembled using CLC Assembly Cell v.4.0 for contig assembly and SSPACE for scaffolding. The assembly generated approximately 6,545 scaffolds with an N50 value of 92,945 and total size of ~ 193 Mb. For gene predictions, the PASA and GeneWise models were used, and ab initio gene predictions were performed independently, resulting in the prediction of 15,842 genes. RNA-Seq expression profiles revealed constitutive expression of 11,572 unique protein-coding genes in larva, 12,364 in protonymph, 12,872 in male adult, and 12,617 in female adult stages. Of the 15,842 predicted genes, 10,885 were commonly expressed through all L. pallidum stages. Genes selectively over-transcribed in the larval stage, which is when host parasitization and disease transmission occur, were further annotated, and their putative roles were discussed.

De Novo Assembly of a High-Quality Reference Genome for the Horned Lark (Eremophila alpestris).

The Horned Lark (Eremophila alpestris) is a small songbird that exhibits remarkable geographic variation in appearance and habitat across an expansive distribution. While E. alpestris has been the focus of many ecological and evolutionary studies, we still lack a highly contiguous genome assembly for the Horned Lark and related taxa (Alaudidae). Here, we present CLO_EAlp_1.0, a highly contiguous assembly for E. alpestris generated from a blood sample of a wild, male bird captured in the Altiplano Cundiboyacense of Colombia. By combining short-insert and mate-pair libraries with the ALLPATHS-LG genome assembly pipeline, we generated a 1.04 Gb assembly comprised of 2713 scaffolds, with a largest scaffold size of 31.81 Mb, a scaffold N50 of 9.42 Mb, and a scaffold L50 of 30. These scaffolds were assembled from 23685 contigs, with a largest contig size of 1.69 Mb, a contig N50 of 193.81 kb, and a contig L50 of 1429. Our assembly pipeline also produced a single mitochondrial DNA contig of 14.00 kb. After polishing the genome, we identified 94.5% of single-copy gene orthologs from an Aves data set and 97.7% of single-copy gene orthologs from a vertebrata data set, which further demonstrates the high quality of our assembly. We anticipate that this genomic resource will be useful to the broader ornithological community and those interested in studying the evolutionary history and ecological interactions of larks, which comprise a widespread, yet understudied lineage of songbirds.

The Genome Sequence of the Eastern Woodchuck (Marmota monax) – A Preclinical Animal Model for Chronic Hepatitis B

The Eastern woodchuck (Marmota monax) has been extensively used in research of chronic hepatitis B and liver cancer because its infection with the woodchuck hepatitis virus closely resembles a human hepatitis B virus infection. Development of novel immunotherapeutic approaches requires genetic information on immune pathway genes in this animal model. The woodchuck genome was assembled with a combination of high-coverage whole-genome shotgun sequencing of Illumina paired-end, mate-pair libraries and fosmid pool sequencing. The result is a 2.63 Gigabase (Gb) assembly with a contig N50 of 74.5 kilobases (kb), scaffold N50 of 892 kb, and genome completeness of 99.2%. RNA sequencing (RNA-seq) from seven different tissues aided in the annotation of 30,873 protein-coding genes, which in turn encode 41,826 unique protein products. More than 90% of the genes have been functionally annotated, with 82% of them containing open reading frames. This genome sequence and its annotation will enable further research in chronic hepatitis B and hepatocellular carcinoma and contribute to the understanding of immunological responses in the woodchuck.

Genomic analyses reveal an absence of contemporary introgressive admixture between fin whales and blue whales, despite known hybrids.

Fin whales (Balaenoptera physalus) and blue whales (B. musculus) are the two largest species on Earth and are widely distributed across the world’s oceans. Hybrids between these species appear to be relatively widespread and have been reported in both the North Atlantic and North Pacific; they are also relatively common, and have been proposed to occur once in every thousand fin whales. However, despite known hybridization, fin and blue whales are not sibling species. Rather, the closest living relative of fin whales are humpback whales (Megaptera novaeangliae). To improve the quality of fin whale data available for analysis, we assembled and annotated a fin whale nuclear genome using in-silico mate pair libraries and previously published short-read data. Using this assembly and genomic data from a humpback, blue, and bowhead whale, we investigated whether signatures of introgression between the fin and blue whale could be found. We find no signatures of contemporary admixture in the fin and blue whale genomes, although our analyses support ancestral gene flow between the species until 2.4–1.3 Ma. We propose the following explanations for our findings; i) fin/blue whale hybridization does not occur in the populations our samples originate from, ii) contemporary hybrids are a recent phenomenon and the genetic consequences have yet to become widespread across populations, or iii) fin/blue whale hybrids are under large negative selection, preventing them from backcrossing and contributing to the parental gene pools.

Genome Sequencing of Pleurozium schreberi: The Assembled and Annotated Draft Genome of a Pleurocarpous Feather Moss

The pleurocarpous feather moss Pleurozium schreberi is a ubiquitous moss species which plays a fundamental role in many terrestrial ecosystems, for instance within the boreal forest, the Earth’s largest terrestrial biome, this species plays a significant role in driving ecosystem nitrogen and carbon inputs and fluxes. By hosting dinitrogen (N2)-fixing cyanobacteria, the moss-cyanobacteria symbiosis constitutes the main nitrogen input into the ecosystem and by the high productivity and the low decomposability of the moss litter, P. schreberi contributes significantly to build-up soil organic matter, and therefore long-term C sequestration. Knowledge on P. schreberi genome will facilitate the development of ‘omics’ and system’s biology approaches to gain a more complete understanding of the physiology and ecological adaptation of the moss and the mechanisms underpinning the establishment of the symbiosis. Here we present the de novo assembly and annotation of P. schreberi genome that will help investigating these questions. The sequencing was performed using the HiSeq X platform with Illumina paired-end and mate-pair libraries prepared with CTAB extracted DNA. In total, the assembled genome was approximately 318 Mb, while repetitive elements account for 28.42% of the genome and 15,992 protein-coding genes were predicted from the genome, of which 84.23% have been functionally annotated. We anticipate that the genomic data generated will constitute a significant resource to study ecological and evolutionary genomics of P. schreberi, and will be valuable for evo-devo investigations as well as our understanding of the evolution of land plants by providing the genome of a pleurocarpous moss.

A Reference Genome Sequence for the European Silver Fir (Abies alba Mill.): A Community-Generated Genomic Resource.

Silver fir (Abies alba Mill.) is a keystone conifer of European montane forest ecosystems that has experienced large fluctuations in population size during during the Quaternary and, more recently, due to land-use change. To forecast the species’ future distribution and survival, it is important to investigate the genetic basis of adaptation to environmental change, notably to extreme events. For this purpose, we here provide a first draft genome assembly and annotation of the silver fir genome, established through a community-based initiative. DNA obtained from haploid megagametophyte and diploid needle tissue was used to construct and sequence Illumina paired-end and mate-pair libraries, respectively, to high depth. The assembled A. alba genome sequence accounted for over 37 million scaffolds corresponding to 18.16 Gb, with a scaffold N50 of 14,051 bp. Despite the fragmented nature of the assembly, a total of 50,757 full-length genes were functionally annotated in the nuclear genome. The chloroplast genome was also assembled into a single scaffold (120,908 bp) that shows a high collinearity with both the A. koreana and A. sibirica complete chloroplast genomes. This first genome assembly of silver fir is an important genomic resource that is now publicly available in support of a new generation of research. By genome-enabling this important conifer, this resource will open the gate for new research and more precise genetic monitoring of European silver fir forests.

A draft nuclear-genome assembly of the acoel flatworm Praesagittifera naikaiensis.

BackgroundAcoels are primitive bilaterians with very simple soft bodies, in which many organs, including the gut, are not developed. They provide platforms for studying molecular and developmental mechanisms involved in the formation of the basic bilaterian body plan, whole-body regeneration, and symbiosis with photosynthetic microalgae. Because genomic information is essential for future research on acoel biology, we sequenced and assembled the nuclear genome of an acoel, Praesagittifera naikaiensis.FindingsTo avoid sequence contamination derived from symbiotic microalgae, DNA was extracted from embryos that were free of algae. More than 290x sequencing coverage was achieved using a combination of Illumina (paired-end and mate-pair libraries) and PacBio sequencing. RNA sequencing and Iso-Seq data from embryos, larvae, and adults were also obtained. First, a preliminary ∼17–kilobase pair (kb) mitochondrial genome was assembled, which was deleted from the nuclear sequence assembly. As a result, a draft nuclear genome assembly was ∼656 Mb in length, with a scaffold N50 of 117 kb and a contig N50 of 57 kb. Although ∼70% of the assembled sequences were likely composed of repetitive sequences that include DNA transposons and retrotransposons, the draft genome was estimated to contain 22,143 protein-coding genes, ∼99% of which were substantiated by corresponding transcripts. We could not find horizontally transferred microalgal genes in the acoel genome. Benchmarking Universal Single-Copy Orthologs analyses indicated that 77% of the conserved single-copy genes were complete. Pfam domain analyses provided a basic set of gene families for transcription factors and signaling molecules.ConclusionsOur present sequencing and assembly of the P. naikaiensis nuclear genome are comparable to those of other metazoan genomes, providing basic information for future studies of genic and genomic attributes of this animal group. Such studies may shed light on the origins and evolution of simple bilaterians.

Extensive and deep sequencing of the Venter/HuRef genome for developing and benchmarking genome analysis tools

We produced an extensive collection of deep re-sequencing datasets for the Venter/HuRef genome using the Illumina massively-parallel DNA sequencing platform. The original Venter genome sequence is a very-high quality phased assembly based on Sanger sequencing. Therefore, researchers developing novel computational tools for the analysis of human genome sequence variation for the dominant Illumina sequencing technology can test and hone their algorithms by making variant calls from these Venter/HuRef datasets and then immediately confirm the detected variants in the Sanger assembly, freeing them of the need for further experimental validation. This process also applies to implementing and benchmarking existing genome analysis pipelines. We prepared and sequenced 200 bp and 350 bp short-insert whole-genome sequencing libraries (sequenced to 100x and 40x genomic coverages respectively) as well as 2 kb, 5 kb, and 12 kb mate-pair libraries (49x, 122x, and 145x physical coverages respectively). Lastly, we produced a linked-read library (128x physical coverage) from which we also performed haplotype phasing.

SGTK: a toolkit for visualization and assessment of scaffold graphs.

Scaffolding is an important step in every genome assembly pipeline, which allows to order contigs into longer sequences using various types of linkage information, such as mate-pair libraries and long reads. In this work, we operate with a notion of a scaffold graph-a graph, vertices of which correspond to the assembled contigs and edges represent connections between them. We present a software package called Scaffold Graph ToolKit that allows to construct and visualize scaffold graphs using different kinds of sequencing data. We show that the scaffold graph appears to be useful for analyzing and assessing genome assemblies, and demonstrate several use cases that can be helpful for both assembly software developers and their users. SGTK is implemented in C++, Python and JavaScript and is freely available at https://github.com/olga24912/SGTK. Supplementary data are available at Bioinformatics online.

Genomes of trombidid mites reveal novel predicted allergens and laterally transferred genes associated with secondary metabolism.

BackgroundTrombidid mites have a unique life cycle in which only the larval stage is ectoparasitic. In the superfamily Trombiculoidea (“chiggers”), the larvae feed preferentially on vertebrates, including humans. Species in the genus Leptotrombidium are vectors of a potentially fatal bacterial infection, scrub typhus, that affects 1 million people annually. Moreover, chiggers can cause pruritic dermatitis (trombiculiasis) in humans and domesticated animals. In the Trombidioidea (velvet mites), the larvae feed on other arthropods and are potential biological control agents for agricultural pests. Here, we present the first trombidid mites genomes, obtained both for a chigger, Leptotrombidium deliense, and for a velvet mite, Dinothrombium tinctorium.ResultsSequencing was performed using Illumina technology. A 180 Mb draft assembly for D. tinctorium was generated from two paired-end and one mate-pair library using a single adult specimen. For L. deliense, a lower-coverage draft assembly (117 Mb) was obtained using pooled, engorged larvae with a single paired-end library. Remarkably, both genomes exhibited evidence of ancient lateral gene transfer from soil-derived bacteria or fungi. The transferred genes confer functions that are rare in animals, including terpene and carotenoid synthesis. Thirty-seven allergenic protein families were predicted in the L. deliense genome, of which nine were unique. Preliminary proteomic analyses identified several of these putative allergens in larvae.ConclusionsTrombidid mite genomes appear to be more dynamic than those of other acariform mites. A priority for future research is to determine the biological function of terpene synthesis in this taxon and its potential for exploitation in disease control.