Low-coverage Genome Research Articles

Multiomics analysis of pediatric solid tumors within the INFORM precision oncology study: From functional drug profiling to biomarker identification.

10036 Background: Within the INFORM (INdividualized Therapy FOr Relapsed Malignancies in Childhood) registry over 1200 childhood malignancies were molecularly profiled using next generation sequencing to identify therapeutic targets. Nevertheless, high evidence targets were only detected in 5% of cases and only 50% of the patients were identified with druggable pathways, while the remaining cases lacked druggable alterations. Thus, an ex-vivo functional drug response profiling platform for pediatric solid tumors has been established within the INFORM program to identify novel biomarkers and unravel molecular mechanisms associated with drug response profiles for clinical translation. Methods: Solid tumors from 97 paediatric patients were screened against a library of 76 drugs. A quality control classification decision tree was designed to select the samples for further analysis. Whole exome sequencing, low-coverage whole genome sequencing, and RNA-seq were used to analyze the molecular profiles of the samples. Molecular features were filtered to retain oncogenic and druggable events. Five data feature views (mutations, mRNA, gene fusions, CNVs, and drug responses) were used to train a Multi-Omics Factor Analysis (MOFA) model to identify prominent latent factors. Results: 81 samples passed the quality control inclusion criteria of which 76 samples had available omics data profiles. Quantitative drug profiling measurements were reported using the selective asymmetric drug sensitivity score.The multi-omics analysis captured five entity-specific omics signatures of Ewing sarcoma, Wilms tumor, BCOR sarcoma, neuroblastoma and ependymoma. Moreover, the analysis revealed sensitivity to navitoclax in neuroblastoma samples with PHOX2B-GATA3 overexpression, and an association between MEK inhibitor sensitivity and an expression signature in Wilms tumors. Conclusions: The combination of functional drug and multi-omics profiling enabled the identification of novel biomarkers for drug sensitivities in pediatric solid tumors. As we continue to expand the number of patient samples evaluated with our drug sensitivity platform, this dataset will provide insights for novel drug targets, and could unravel key molecular events and mechanisms acting towards personalized therapies.

Testing the systematic status of Homalictus and Rostrohalictus with weakened cross‐vein groups within Halictini (Hymenoptera: Halictidae) using low‐coverage whole‐genome sequencing

The halictid genus Lasioglossum, as one of the most species-rich bee groups with persistently contentious subgeneric boundaries, is one of the most challenging bee groups from a systematic standpoint. An enduring question is the relationship of Lasioglossum and Homalictus, whether all halictine bees with weakened distal wing venation comprise one or multiple genera. Here, we analyzed the phylogenetic relationships among the subgroups within Lasioglossum s.l. based on thousands of single-copy orthologs and ultraconserved elements, which were extracted from 23 newly sequenced low-coverage whole genomes alongside a published genome (22 ingroups plus 2 outgroups). Both marker sets provided consistent results across maximum likelihood and coalescent-based species tree approaches. The phylogenetic and topology test results show that the Lasioglossum and Hemihalictus series are reciprocally monophyletic and Homalictus and Rostrohalictus are valid subgenera of Lasioglossum. Consequently, we lower Homalictus to subgenus status within Lasioglossum again, and we also raise Rostrohalictus to subgenus status from its prior synonymy with subgenus Hemihalictus. Lasioglossum przewalskyi is also transferred to the subgenus Hemihalictus. Ultimately, we redefine Lasioglossum to include all halictine bees with weakened distal wing venation.

Using ultraconserved elements to reconstruct the termite tree of life

The phylogenetic history of termites has been investigated using mitochondrial genomes and transcriptomes. However, both sets of markers have specific limitations. Mitochondrial genomes represent a single genetic marker likely to yield phylogenetic trees presenting incongruences with species trees, and transcriptomes can only be obtained from well-preserved samples. In contrast, ultraconserved elements (UCEs) include a great many independent markers that can be retrieved from poorly preserved samples. Here, we designed termite-specific baits targeting 50,616 UCE loci. We tested our UCE bait set on 42 samples of termites and three samples of Cryptocercus, for which we generated low-coverage highly-fragmented genome assemblies and successfully extracted in silico between 3,426 to 42,860 non-duplicated UCEs per sample. Our maximum likelihood phylogenetic tree, reconstructed using the 5,934 UCE loci retrieved from upward of 75% of samples, was congruent with transcriptome-based phylogenies, demonstrating that our UCE bait set is reliable and phylogenetically informative. Combined with non-destructive DNA extraction protocols, our UCE bait set provides the tool needed to carry out a global taxonomic revision of termites based on poorly preserved specimens such as old museum samples. The Termite UCE database is maintained at: https://github.com/oist/TER-UCE-DB/.

Copy number variation of urine exfoliated cells by low-coverage whole genome sequencing for diagnosis of prostate adenocarcinoma: a prospective cohort study.

BackgroundNon-invasive, especially the urine-based diagnosis of prostate cancer (PCa) remains challenging. Although prostate cancer antigen (PSA) is widely used in prostate cancer screening, the false positives may result in unnecessary invasive procedures. PSA elevated patients are triaged to further evaluation of free/total PSA ratio (f/t PSA), to find out potential clinically significant PCa before undergoing invasive procedures. Genomic instability, especially chromosomal copy number variations (CNVs) were proved much more tumor specific. Here we performed a prospective study to evaluate the diagnostic value of CNV via urine-exfoliated cell DNA analysis in PCa.MethodsWe enrolled 28 PSA elevated patients (≥ 4 ng/ml), including 16 PCa, 9 benign prostate hypertrophy (BPH) and 3 prostatic intraepithelial neoplasia (PIN). Fresh initial portion urine was collected after hospital admission. Urine exfoliated cell DNA was analyzed by low coverage Whole Genome Sequencing, followed by CNV genotyping by the prostate cancer chromosomal aneuploidy detector (ProCAD). CNVs were quantified in absolute z-score (|Z|). Serum free/total PSA ratio (f/t PSA) was reported altogether.ResultsIn patients with PCa, the most frequent CNV events were chr3q gain (n = 2), chr8q gain (n = 2), chr2q loss (n = 4), and chr18q loss (n = 3). CNVs were found in 81.2% (95% Confidence Interval (CI) 53.7–95.0%) PCa. No CNV was identified in BPH patients. A diagnosis model was established by incorporating all CNVs. At the optimal cutoff of |Z|≥ 2.50, the model reached an AUC of 0.91 (95% CI 0.83–0.99), a sensitivity of 81.2% and a specificity of 100%. The CNV approach significantly outperformed f/t PSA (AUC = 0.62, P = 0.012). Further analyses showed that the CNV positive rate was significantly correlated with tumor grade. CNVs were found in 90.9% (95% CI 57.1–99.5%) high grade tumors and 60.0% (95% CI 17.0–92.7%) low grade tumors. No statistical significance was found for patient age, BMI, disease history and family history.ConclusionsUrine exfoliated cells harbor enriched CNV features in PCa patients. Urine detection of CNV might be a biomarker for PCa diagnosis, especially in terms of the clinically significant high-grade tumors.Supplementary InformationThe online version contains supplementary material available at 10.1186/s12920-022-01253-5.

Comparison of SARS-CoV-2 sequencing using the ONT GridION and the Illumina MiSeq

BackgroundOver 4 million SARS-CoV-2 genomes have been sequenced globally in the past 2 years. This has been crucial in elucidating transmission chains within communities, the development of new diagnostic methods, vaccines, and antivirals. Although several sequencing technologies have been employed, Illumina and Oxford Nanopore remain the two most commonly used platforms. The sequence quality between these two platforms warrants a comparison of the genomes produced by the two technologies. Here, we compared the SARS-CoV-2 consensus genomes obtained from the Oxford Nanopore Technology GridION and the Illumina MiSeq for 28 sequencing runs.ResultsOur results show that the MiSeq had a significantly higher number of consensus genomes classified by Nextclade as good and mediocre compared to the GridION. The MiSeq also had a significantly higher genome coverage and mutation counts than the GridION.ConclusionDue to the low genome coverage, high number of indels, and sensitivity to SARS-CoV-2 viral load noted with the GridION when compared to MiSeq, we can conclude that the MiSeq is more favourable for SARS-CoV-2 genomic surveillance, as successful genomic surveillance is dependent on high quality, near-whole consensus genomes.

Genomic evaluation of hybridization in historic and modern North American Bison (Bison bison)

During the late nineteenth century North American bison underwent a significant population bottleneck resulting in a reduction in population size of over 99% and a species-level near-extinction event. Factors responsible for this destruction included indiscriminate killing, loss of access to suitable habitat, and diseases. At the nadir of this population crash, very few wild plains bison survived and were restricted to Yellowstone National Park, USA and a small number of wild wood bison remained in Wood Buffalo National Park, Canada. However, most surviving bison in the late 1800’s were maintained by cattle ranchers in private herds where hybridization between bison with various breeds of domestic cattle was often encouraged. Over the last 20 years, the legacy of this introgression has been identified using mitochondrial DNA and limited nuclear microsatellite analyses. However, no genome-wide assessment has been performed, and some herds were believed to be free of introgression based on current genetic testing strategies. Herein, we report detailed analyses using whole genome sequencing from nineteen modern and six historical bison, chosen to represent the major lineages of bison, to identify and quantitate signatures of nuclear introgression in their recent (within 200 years) history. Both low and high coverage genomes provided evidence for recent introgression, including animals from Yellowstone, Wind Cave, and Elk Island National Parks which were previously thought to be free from hybridization with domestic cattle. We employed multiple approaches, including one developed for this work, to identify putative cattle haplotypes in each bison genome. These regions vary greatly in size and frequency by sample and herd, though we detected domestic cattle introgression in all bison genomes tested. Since our sampling strategy spanned across the diversity of modern bison populations, these finding are best explained by multiple historical hybridization events between these two species with significant genetic recombination over the last 200 years. Our results demonstrate that whole genome sequencing approaches are required to accurately quantitate cattle introgression in bison.

A setback into a success: What can batch effects tell us about best practices in genomics?

The increasing access to high-throughput sequencing is certainly one of the major changes that molecular ecology has gone through over the last decade. With the positive trend towards more open science, most sequencing data sets are now available on public databases, which holds amazing potential, but also risks of introducing batch effects in studies combining data sets. In this issue of Molecular Ecology Resources, Lou and Therkildsen (2022) offer a timely discussion on the matter by analyzing an imperfect low-coverage Whole Genome Sequencing data set, in which they test the effects of differences in sequencing choices, DNA degradation, and read depth on routine population genomics analyses. Through a series of diagnostic tools, they uncover multiple factors producing technical artefacts that can bias estimates of genetic diversity, inference of population structure, and selection scans. For each confounding factor, they demonstrate the effectiveness of mitigation approaches and suggest other avenues to deal with the issue. In this perspective, we highlight considerations regarding (1) effects that arise from differences between batches of sequencing; (2) unavoidable heterogeneity within data sets; and (3) more general concerns around the use of next-generation sequencing in population genomics. Altogether, by exploring what may have appeared at first glimpse as a "failed" sequencing project, Lou and Therkildsen (2022) end up setting a standard of best practices to make the most of heterogeneous whole-genome sequences, opening a promising avenue towards efficient reuse of published data sets.

Prevalence and circulation patterns of SARS-CoV-2 variants in European sewage mirror clinical data of 54 European cities

For community-level monitoring, the European Commission under the EU Sewage Sentinel System recommends wastewater-based SARS-CoV-2 surveillance. Tracking SARS-CoV-2 variants in a community is pivotal for appropriate public health response. Genome sequencing of SARS-CoV-2 in wastewater samples for tracking variants is challenging, often resulting in low coverage genome sequences, thereby impeding the detection of the SARS-CoV-2 mutations. Therefore, we aimed at high-coverage SARS-CoV-2 genome sequences from sewage samples which we successfully accomplished. This first pan-European surveillance compared the mutation profiles associated with the variants of concerns: B.1.1.7, P.1, B.1.351 and B.1.617.2 across 20 European countries, including 54 municipalities. The results highlight that SARS-CoV-2 variants detected in the wastewater samples mirror the variants profiles reported in clinical data. This study demonstrated that >98% coverage of SARS-CoV-2 genomic sequences is possible and can be used to track SARS-CoV-2 mutations in wastewater to support identifying variants circulating in a city at the community level.

Biogeography and ecological drivers of evolution in the Andes: resolving the phylogenetic backbone forCalceolaria(Calceolariaceae)

AbstractCalceolaria (Calceolariaceae) is an emblematic and diverse genus in the Americas. Despite being one of the most easily recognized genera in the region and a system with great potential to improve our understanding of different drivers of species diversification in the Andes, its intrageneric evolutionary relationships are still poorly understood. Responding to the need for additional molecular markers to resolve the phylogenetic relationships of the group, we perform plastome analyses and resolve the backbone of the genus. Specifically, using low-coverage genomes for 14 species, we assembled plastomes, estimated and dated phylogenetic hypotheses and evaluated evolutionary trends in the group. Our approach allowed us to resolve the backbone of the genus, identify two main clades and estimate a timing of diversification contemporaneous to major climatic and orogenic events. Our biogeographic reconstructions suggest an independent colonization of the whole range of the genus by both clades. Finally, our evaluations of floral morphology reveal future avenues for investigating the relationship between the pollination biology and diversification of the group.

ORTHOSKIM: In silico sequence capture from genomic and transcriptomic libraries for phylogenomic and barcoding applications.

Low-coverage whole genome shotgun sequencing (or genome skimming) has emerged as a cost-effective method for acquiring genomic data in nonmodel organisms. This method provides sequence information on chloroplast genome (cpDNA), mitochondrial genome (mtDNA) and nuclear ribosomal regions (rDNA), which are over-represented within cells. However, numerous bioinformatic challenges remain to accurately and rapidly obtain such data in organisms with complex genomic structures and rearrangements, in particular for mtDNA in plants or for cpDNA in some plant families. Here we introduce the pipeline ORTHOSKIM, which performs in silico capture of targeted sequences from genomic and transcriptomic libraries without assembling whole organelle genomes. ORTHOSKIM proceeds in three steps: (i) global sequence assembly, (ii) mapping against reference sequences and (iii) target sequence extraction; importantly it also includes a range of quality control tests. Different modes are implemented to capture both coding and noncoding regions of cpDNA, mtDNA and rDNA sequences, along with predefined nuclear sequences (e.g., ultraconserved elements) or collections of single-copy orthologue genes. Moreover, aligned DNA matrices are produced for phylogenetic reconstructions, by performing multiple alignments of the captured sequences. While ORTHOSKIM is suitable for any eukaryote, a case study is presented here, using 114 genome-skimming libraries and four RNA sequencing libraries obtained for two plant families, Primulaceae and Ericaceae, the latter being a well-known problematic family for cpDNA assemblies. ORTHOSKIM recovered with high success rates cpDNA, mtDNA and rDNA sequences, well suited to accurately infer evolutionary relationships within these families. ORTHOSKIM is released under a GPL-3 licence and is available at: https://github.com/cpouchon/ORTHOSKIM.

Oceanic Crustal Fluid Single Cell Genomics Complements Metagenomic and Metatranscriptomic Surveys With Orders of Magnitude Less Sample Volume.

Fluids circulating through oceanic crust play important roles in global biogeochemical cycling mediated by their microbial inhabitants, but studying these sites is challenged by sampling logistics and low biomass. Borehole observatories installed at the North Pond study site on the western flank of the Mid-Atlantic Ridge have enabled investigation of the microbial biosphere in cold, oxygenated basaltic oceanic crust. Here we test a methodology that applies redox-sensitive fluorescent molecules for flow cytometric sorting of cells for single cell genomic sequencing from small volumes of low biomass (approximately 103 cells ml–1) crustal fluid. We compare the resulting genomic data to a recently published paired metagenomic and metatranscriptomic analysis from the same site. Even with low coverage genome sequencing, sorting cells from less than one milliliter of crustal fluid results in similar interpretation of dominant taxa and functional profiles as compared to ‘omics analysis that typically filter orders of magnitude more fluid volume. The diverse community dominated by Gammaproteobacteria, Bacteroidetes, Desulfobacterota, Alphaproteobacteria, and Zetaproteobacteria, had evidence of autotrophy and heterotrophy, a variety of nitrogen and sulfur cycling metabolisms, and motility. Together, results indicate fluorescence activated cell sorting methodology is a powerful addition to the toolbox for the study of low biomass systems or at sites where only small sample volumes are available for analysis.

Whole genome deep sequencing analysis of cell-free DNA in samples with low tumour content

BackgroundCirculating cell-free DNA (cfDNA) in the plasma of cancer patients contains cell-free tumour DNA (ctDNA) derived from tumour cells and it has been widely recognized as a non-invasive source of tumour DNA for diagnosis and prognosis of cancer. Molecular profiling of ctDNA is often performed using targeted sequencing or low-coverage whole genome sequencing (WGS) to identify tumour specific somatic mutations or somatic copy number aberrations (sCNAs). However, these approaches cannot efficiently detect all tumour-derived genomic changes in ctDNA.MethodsWe performed WGS analysis of cfDNA from 4 breast cancer patients and 2 patients with benign tumours. We sequenced matched germline DNA for all 6 patients and tumour samples from the breast cancer patients. All samples were sequenced on Illumina HiSeqXTen sequencing platform and achieved approximately 30x, 60x and 100x coverage on germline, tumour and plasma DNA samples, respectively.ResultsThe mutational burden of the plasma samples (1.44 somatic mutations/Mb of genome) was higher than the matched tumour samples. However, 90% of high confidence somatic cfDNA variants were not detected in matched tumour samples and were found to comprise two background plasma mutational signatures. In contrast, cfDNA from the di-nucleosome fraction (300 bp–350 bp) had much higher proportion (30%) of variants shared with tumour. Despite high coverage sequencing we were unable to detect sCNAs in plasma samples.ConclusionsDeep sequencing analysis of plasma samples revealed higher fraction of unique somatic mutations in plasma samples, which were not detected in matched tumour samples. Sequencing of di-nucleosome bound cfDNA fragments may increase recovery of tumour mutations from plasma.

Batch effects in population genomic studies with low-coverage whole genome sequencing data: Causes, detection and mitigation.

Over the past few decades, there has been an explosion in the amount of publicly available sequencing data. This opens new opportunities for combining data sets to achieve unprecedented sample sizes, spatial coverage or temporal replication in population genomic studies. However, a common concern is that nonbiological differences between data sets may generate patterns of variation in the data that can confound real biological patterns, a problem known as batch effects. In this paper, we compare two batches of low-coverage whole genome sequencing (lcWGS) data generated from the same populations of Atlantic cod (Gadus morhua). First, we show that with a "batch-effect-naive" bioinformatic pipeline, batch effects systematically biased our genetic diversity estimates, population structure inference and selection scans. We then demonstrate that these batch effects resulted from multiple technical differences between our data sets, including the sequencing chemistry (four-channel vs. two-channel), sequencing run, read type (single-end vs. paired-end), read length (125 vs. 150bp), DNA degradation level (degraded vs. well preserved) and sequencing depth (0.8× vs. 0.3× on average). Lastly, we illustrate that a set of simple bioinformatic strategies (such as different read trimming and single nucleotide polymorphism filtering) can be used to detect batch effects in our data and substantially mitigate their impact. We conclude that combining data sets remains a powerful approach as long as batch effects are explicitly accounted for. We focus on lcWGS data in this paper, which may be particularly vulnerable to certain causes of batch effects, but many of our conclusions also apply to other sequencing strategies.

Quantified CIN Score From Cell-free DNA as a Novel Noninvasive Predictor of Survival in Patients With Spinal Metastasis.

Purpose: Most currently available scores for survival prediction of patients with bone metastasis lack accuracy. In this study, we present a novel quantified CIN (Chromosome Instability) score modeled from cfDNA copy number variation (CNV) for survival prediction. Experimental Design: Plasma samples collected from 67 patients with bone metastases from 11 different cancer types between November 2015 and May 2016 were sent through low-coverage whole genome sequencing followed by CIN computation to make a correlation analysis between the CIN score and survival prognosis. The results were validated in an independent cohort of 213 patients. Results: During the median follow-up period of 598 (95% CI 364–832) days until December 25, 2018, 124 (44.3%) of the total 280 patients died. Analysis of the discovery dataset showed that CIN score = 12 was the optimal CIN cutoff. Validation dataset showed that CIN was elevated (score ≥12) in 87 (40.8%) patients, including 5 (5.75%) with head and neck cancer, 11 (12.6%) with liver and gallbladder cancer, 11 (12.6%) with cancer from unidentified sites, 21 (24.1%) with lung cancer, 7 (8.05%) with breast cancer, 4 (4.60%) with thyroid cancer, 6 (6.90%) with colorectal cancer, 4 (4.60%) with kidney cancer, 2 (2.30%) with prostate cancer, and 16 (18.4%) with other types of cancer. Further analysis showed that patients with elevated CIN were associated with worse survival (p < 0.001). For patients with low Tokuhashi score (≤8) who had predictive survival of less than 6 months, the CIN score was able to distinguish patients with a median overall survival (OS) of 443 days (95% CI 301–585) from those with a median OS of 258 days (95% CI 184–332). Conclusion: CNV examination in bone metastatic cancer from cfDNA is superior to the traditional predictive model in that it provides a noninvasive and objective method of monitoring the survival of patients with spine metastasis.

Comprehensive and accurate genetic variant identification from contaminated and low-coverage Mycobacterium tuberculosis whole genome sequencing data.

Improved understanding of the genomic variants that allow Mycobacterium tuberculosis (Mtb) to acquire drug resistance, or tolerance, and increase its virulence are important factors in controlling the current tuberculosis epidemic. Current approaches to Mtb sequencing, however, cannot reveal Mtb’s full genomic diversity due to the strict requirements of low contamination levels, high Mtb sequence coverage and elimination of complex regions. We have developed the XBS (compleX Bacterial Samples) bioinformatics pipeline, which implements joint calling and machine-learning-based variant filtering tools to specifically improve variant detection in the important Mtb samples that do not meet these criteria, such as those from unbiased sputum samples. Using novel simulated datasets, which permit exact accuracy verification, XBS was compared to the UVP and MTBseq pipelines. Accuracy statistics showed that all three pipelines performed equally well for sequence data that resemble those obtained from culture isolates of high depth of coverage and low-level contamination. In the complex genomic regions, however, XBS accurately identified 9.0 % more SNPs and 8.1 % more single nucleotide insertions and deletions than the WHO-endorsed unified analysis variant pipeline. XBS also had superior accuracy for sequence data that resemble those obtained directly from sputum samples, where depth of coverage is typically very low and contamination levels are high. XBS was the only pipeline not affected by low depth of coverage (5–10×), type of contamination and excessive contamination levels (>50 %). Simulation results were confirmed using whole genome sequencing (WGS) data from clinical samples, confirming the superior performance of XBS with a higher sensitivity (98.8%) when analysing culture isolates and identification of 13.9 % more variable sites in WGS data from sputum samples as compared to MTBseq, without evidence for false positive variants when rRNA regions were excluded. The XBS pipeline facilitates sequencing of less-than-perfect Mtb samples. These advances will benefit future clinical applications of Mtb sequencing, especially WGS directly from clinical specimens, thereby avoiding in vitro biases and making many more samples available for drug resistance and other genomic analyses. The additional genetic resolution and increased sample success rate will improve genome-wide association studies and sequence-based transmission studies.

Estimating repeat spectra and genome length from low-coverage genome skims with RESPECT.

The cost of sequencing the genome is dropping at a much faster rate compared to assembling and finishing the genome. The use of lightly sampled genomes (genome-skims) could be transformative for genomic ecology, and results using k-mers have shown the advantage of this approach in identification and phylogenetic placement of eukaryotic species. Here, we revisit the basic question of estimating genomic parameters such as genome length, coverage, and repeat structure, focusing specifically on estimating the k-mer repeat spectrum. We show using a mix of theoretical and empirical analysis that there are fundamental limitations to estimating the k-mer spectra due to ill-conditioned systems, and that has implications for other genomic parameters. We get around this problem using a novel constrained optimization approach (Spline Linear Programming), where the constraints are learned empirically. On reads simulated at 1X coverage from 66 genomes, our method, REPeat SPECTra Estimation (RESPECT), had 2.2% error in length estimation compared to 27% error previously achieved. In shotgun sequenced read samples with contaminants, RESPECT length estimates had median error 4%, in contrast to other methods that had median error 80%. Together, the results suggest that low-pass genomic sequencing can yield reliable estimates of the length and repeat content of the genome. The RESPECT software will be publicly available at https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_shahab-2Dsarmashghi_RESPECT.git&d=DwIGAw&c=-35OiAkTchMrZOngvJPOeA&r=ZozViWvD1E8PorCkfwYKYQMVKFoEcqLFm4Tg49XnPcA&m=f-xS8GMHKckknkc7Xpp8FJYw_ltUwz5frOw1a5pJ81EpdTOK8xhbYmrN4ZxniM96&s=717o8hLR1JmHFpRPSWG6xdUQTikyUjicjkipjFsKG4w&e=.

Functional DNA annotation from a preliminary de novo genome assembly of Brycon orbignyanus, an endangered Neotropical migratory fish

The predicted sequence for thousands of genes revealed by a preliminary low-coverage genome assembly is presented for Brycon orbignyanus, an endangered migratory fish. Neotropical migratory fish stocks have been drastically reduced due to accumulated environmental pressure. Brycon orbignyanus, once one of the main fisheries species in the Platine Basin, is now very rare in nature and relies on spawning programs and a few well preserved or still untouched sites. The use of high-throughput DNA sequencing is still untapped regarding the functional genome information from B. orbignyanus. In order to help bridging this gap, we present a dataset resulting from the first functional annotation from a de novo genome assembly for B. orbignyanus, from short reads (90 bp), obtained by the HiSeq 2000 platform (Illumina). The annotation was performed for scaffolds over 10 kb using the Maker pipeline, with reference sequences taken from the NCBI for the Characiformes order. This annotation resulted in the prediction of 12,734 genes, classified with the aid of PANTHER. The data presented here can facilitate the development of basic research in this threatened species, along with practical biotechnological tools for different areas, such as commercial and environmental fish spawning operations (e.g. hormonal induction, growth) and human health.

Low-Coverage Whole Genome Sequencing Using Laser Capture Microscopy with Combined Digital Droplet PCR: An Effective Tool to Study Copy Number and Kras Mutations in Early Lung Adenocarcinoma Development.

Defining detailed genomic characterization of early tumor progression is critical to identifying key regulators and pathways in carcinogenesis as potentially druggable targets. In human lung cancer, work to characterize early cancer development has mainly focused on squamous cancer, as the earliest lesions are more proximal in the airways and often accessible by repeated bronchoscopy. Adenocarcinomas are typically located distally in the lung, limiting accessibility for biopsy of pre-malignant and early stages. Mouse lung cancer models recapitulate many human genomic features and provide a model for tumorigenesis with pre-malignant atypical adenomatous hyperplasia and in situ adenocarcinomas often developing contemporaneously within the same animal. Here, we combined tissue characterization and collection by laser capture microscopy (LCM) with digital droplet PCR (ddPCR) and low-coverage whole genome sequencing (LC-WGS). ddPCR can be used to identify specific missense mutations in Kras (Kirsten rat sarcoma viral oncogene homolog, here focused on Kras Q61) and estimate the percentage of mutation predominance. LC-WGS is a cost-effective method to infer localized copy number alterations (CNAs) across the genome using low-input DNA. Combining these methods, the histological stage of lung cancer can be correlated with appearance of Kras mutations and CNAs. The utility of this approach is adaptable to other mouse models of human cancer.

Towards a Cost-Effective Implementation of Genomic Prediction Based on Low Coverage Whole Genome Sequencing in Dezhou Donkey.

Low-coverage whole genome sequencing is a low-cost genotyping technology. Combined with genotype imputation approaches, it is likely to become a critical component of cost-effective genomic selection programs in agricultural livestock. Here, we used the low-coverage sequence data of 617 Dezhou donkeys to investigate the performance of genotype imputation for low-coverage whole genome sequence data and genomic prediction based on the imputed genotype data. The specific aims were as follows: 1) to measure the accuracy of genotype imputation under different sequencing depths, sample sizes, minor allele frequency (MAF), and imputation pipelines and 2) to assess the accuracy of genomic prediction under different marker densities derived from the imputed sequence data, different strategies for constructing the genomic relationship matrixes, and single-vs. multi-trait models. We found that a high imputation accuracy (>0.95) can be achieved for sequence data with a sequencing depth as low as 1x and the number of sequenced individuals ≥400. For genomic prediction, the best performance was obtained by using a marker density of 410K and a G matrix constructed using expected marker dosages. Multi-trait genomic best linear unbiased prediction (GBLUP) performed better than single-trait GBLUP. Our study demonstrates that low-coverage whole genome sequencing would be a cost-effective approach for genomic prediction in Dezhou donkey.

Comparative genomic and transcriptomic analyses of transposable elements in polychaetous annelids highlight LTR retrotransposon diversity and evolution

BackgroundWith the expansion of high throughput sequencing, we now have access to a larger number of genome-wide studies analyzing the Transposable elements (TEs) composition in a wide variety of organisms. However, genomic analyses often remain too limited in number and diversity of species investigated to study in depth the dynamics and evolutionary success of the different types of TEs among metazoans. Therefore, we chose to investigate the use of transcriptomes to describe the diversity of TEs in phylogenetically related species by conducting the first comparative analysis of TEs in two groups of polychaetes and evaluate the diversity of TEs that might impact genomic evolution as a result of their mobility.ResultsWe present a detailed analysis of TEs distribution in transcriptomes extracted from 15 polychaetes depending on the number of reads used during assembly, and also compare these results with additional TE scans on associated low-coverage genomes. We then characterized the clades defined by 1021 LTR-retrotransposon families identified in 26 species. Clade richness was highly dependent on the considered superfamily. Copia elements appear rare and are equally distributed in only three clades, GalEa, Hydra and CoMol. Among the eight BEL/Pao clades identified in annelids, two small clades within the Sailor lineage are new for science. We characterized 17 Gypsy clades of which only 4 are new; the C-clade largely dominates with a quarter of the families. Finally, all species also expressed for the majority two distinct transcripts encoding PIWI proteins, known to be involved in control of TEs mobilities.ConclusionsThis study shows that the use of transcriptomes assembled from 40 million reads was sufficient to access to the diversity and proportion of the transposable elements compared to those obtained by low coverage sequencing. Among LTR-retrotransposons Gypsy elements were unequivocally dominant but results suggest that the number of Gypsy clades, although high, may be more limited than previously thought in metazoans. For BEL/Pao elements, the organization of clades within the Sailor lineage appears more difficult to establish clearly. The Copia elements remain rare and result from the evolutionary consistent success of the same three clades.