Sequencing Depth Research Articles

An overlooked phenomenon: complex interactions of potential error sources on the quality of bacterial de novo genome assemblies

BackgroundParameters adversely affecting the contiguity and accuracy of the assemblies from Illumina next-generation sequencing (NGS) are well described. However, past studies generally focused on their additive effects, overlooking their potential interactions possibly exacerbating one another’s effects in a multiplicative manner. To investigate whether or not they act interactively on de novo genome assembly quality, we simulated sequencing data for 13 bacterial reference genomes, with varying levels of error rate, sequencing depth, PCR and optical duplicate ratios.ResultsWe assessed the quality of assemblies from the simulated sequencing data with a number of contiguity and accuracy metrics, which we used to quantify both additive and multiplicative effects of the four parameters. We found that the tested parameters are engaged in complex interactions, exerting multiplicative, rather than additive, effects on assembly quality. Also, the ratio of non-repeated regions and GC% of the original genomes can shape how the four parameters affect assembly quality.ConclusionsWe provide a framework for consideration in future studies using de novo genome assembly of bacterial genomes, e.g. in choosing the optimal sequencing depth, balancing between its positive effect on contiguity and negative effect on accuracy due to its interaction with error rate. Furthermore, the properties of the genomes to be sequenced also should be taken into account, as they might influence the effects of error sources themselves.

Large-scale single-virus genomics uncovers hidden diversity of river water viruses and diversified gene profiles.

Environmental viruses (primarily bacteriophages) are widely recognized as playing an important role in ecosystem homeostasis through the infection of host cells. However, the majority of environmental viruses are still unknown as their mosaic structure and frequent mutations in their sequences hinder genome construction in current metagenomics. To enable the large-scale acquisition of environmental viral genomes, we developed a new single-viral genome sequencing platform with microfluidic-generated gel beads. Amplification of individual DNA viral genomes in mass-produced gel beads allows high-throughput genome sequencing compared to conventional single-virus genomics. The sequencing analysis of river water samples yielded 1431 diverse viral single-amplified genomes, whereas viral metagenomics recovered 100 viral metagenome-assembled genomes at the comparable sequence depth. The 99.5% of viral single-amplified genomes were determined novel at the species level, most of which could not be recovered by a metagenomic assembly. The large-scale acquisition of diverse viral genomes identified protein clusters commonly detected in different viral strains, allowing the gene transfer to be tracked. Moreover, comparative genomics within the same viral species revealed that the profiles of various methyltransferase subtypes were diverse, suggesting an enhanced escape from host bacterial internal defense mechanisms. Our use of gel bead-based single-virus genomics will contribute to exploring the nature of viruses by accelerating the accumulation of draft genomes of environmental DNA viruses.

Genome stability assessment of PRRS vaccine strain with new ARTIC-style sequencing protocol.

A tiling amplicon sequencing protocol was developed to analyse the genome sequence stability of the modified live PRRSV vaccine strain, Porcilis MLV. The backbone of the ARTIC-style protocol was formed by 34 individual primer pairs, which were divided into two primer pools. Primer pairs were designed to amplify 532 to 588 bp fragments of the corresponding genomic region. The amplicons are suitable for sequencing on Illumina DNA sequencers with available 600-cycle sequencing kits. The concentration of primer pairs in the pools was optimized to obtain a balanced sequencing depth along the genome. Deep sequencing data of three vaccine batches were also analysed. All three vaccine batches were very similar to each other, although they also showed single nucleotide variations (SNVs) affecting less than 1 % of the genome. In the three vaccine strains, 113 to 122 SNV sites were identified; at these sites, the minority variants represented a frequency range of 1 to 48.7 percent. Additionally, the strains within the batches contained well-known length polymorphisms; the genomes of these minority deletion mutants were 135 to 222 bp shorter than the variant with the complete genome. Our results show the usefulness of ARTIC-style protocols in the evaluation of the genomic stability of PRRS MLV strains.

Whole-genome resequencing of Russian sturgeon (Acipenser gueldenstaedtii) reveals selection signatures associated with caviar color

Russian sturgeon, Acipenser gueldenstaedtii, a primitive freshwater fish, holds high economic value as the source of caviar. Remarkably, the economic significance of gold caviar surpasses that of its black counterpart. However, the genetic basis underlying the color variations in caviar remains unexplored. In this study, we conducted whole-genome resequencing of 62 individuals from gold caviar (G) and black caviar (B) populations, achieving an average sequencing depth of approximately 12.21×, and each individual produced approximately 32.47 million high-quality single nucleotide polymorphisms (SNPs). By combining FST and XPEHH methods, we identified 1138 candidate selection regions spanning 26.7 Mb, encompassing 771 genes. Furthermore, significantly enriched GO terms and KEGG pathways were related to pigmentation (e.g., Protein digestion and absorption and Retinol metabolism, etc.) and immunity (e.g., Purine metabolism, Phagosome and Vibrio cholerae infection). Based on the findings, 12 genes (ACTA1, ADCY3, ADCY5, ATP6V1C1, ATP6V1H, CHST11, COL6A3, GART, KCNJ13, RDH10, KCNN4 and XDH) were identified as candidate genes that may affect caviar color. Furthermore, four genes (ATP6V1H, GART, RDH10 and KCNN4) detected simultaneously by FST and XPEHH methods were validated using KASP. This study lays the groundwork for unraveling the molecular mechanisms underlying caviar color and provides insights for the selective breeding of sturgeons with the desired gold caviar trait.

Systematic comparison of variant calling pipelines of target genome sequencing cross multiple next-generation sequencers.

Targeted genomic sequencing (TS) greatly benefits precision oncology by rapidly detecting genetic variations with better accuracy and sensitivity owing to its high sequencing depth. Multiple sequencing platforms and variant calling tools are available for TS, making it excruciating for researchers to choose. Therefore, benchmarking study across different platforms and pipelines available for TS is imperative. In this study, we performed a TS of Reference OncoSpan FFPE (HD832) sample enriched by TSO500 panel using four commercially available sequencers, and analyzed the output 50 datasets using five commonly-used bioinformatics pipelines. We systematically investigated the sequencing quality and variant detection sensitivity, expecting to provide optimal recommendations for future research. Four sequencing platforms returned highly concordant results in terms of base quality (Q20 > 94%), sequencing coverage (>97%) and depth (>2000×). Benchmarking revealed good concordance of variant calling across different platforms and pipelines, among which, FASTASeq 300 platform showed the highest sensitivity (100%) and precision (100%) in high-confidence variants calling when analyzed by SNVer and VarScan 2 algorithms. Furthermore, this sequencer demonstrated the shortest sequencing time (∼21h) at the sequencing mode PE150. Through the intersection of 50 datasets generated in this study, we recommended a novel set of variant genes outside the truth set published by HD832, expecting to replenish HD832 for future research on tumor variant diagnosis. Besides, we applied these five tools to another panel (TargetSeq One) for Twist cfDNA Pan-cancer Reference Standard, comprehensive consideration of SNP and InDel sensitivity, SNVer and VarScan 2 performed best among them. Furthermore, SNVer and VarScan 2 also performed best for six cancer cell lines samples regarding SNP and InDel sensitivity. Considering the dissimilarity of variant calls across different pipelines for datasets from the same platform, we recommended an integration of multiple tools to improve variant calling sensitivity and accuracy for the cancer genome. Illumina and GeneMind technologies can be used independently or together by public health laboratories performing tumor TS. SNVer and VarScan 2 perform better regarding variant detection sensitivity for three typical tumor samples. Our study provides a standardized target sequencing resource to benchmark new bioinformatics protocols and sequencing platforms.

Exploring the cobia (Rachycentron canadum) genome: unveiling putative male heterogametic regions and identification of sex-specific markers.

Cobia (Rachycentron canadum) is the only member of the Rachycentridae family and exhibits considerable sexual dimorphism in growth rate. Sex determination in teleosts has been a long-standing basic biological question, and the molecular mechanisms of sex determination/differentiation in cobia are completely unknown. Here, we reported 2 high-quality, chromosome-level annotated male and female cobia genomes with assembly sizes of 586.51 Mb (contig/scaffold N50: 86.0kb/24.3 Mb) and 583.88 Mb (79.9kb/22.5 Mb), respectively. Synteny inference among perciform genomes revealed that cobia and the remora Echeneis naucrates were sister groups. Further, whole-genome resequencing of 31 males and 60 females, genome-wide association study, and sequencing depth analysis identified 3 short male-specific regions within a 10.7-kb continuous genomic region on male chromosome 18, which hinted at an undifferentiated sex chromosome system with a putative XX/XY mode of sex determination in cobia. Importantly, the only 2 genes within/between the male-specific regions, epoxide hydrolase 1 (ephx1, renamed cephx1y) and transcription factor 24 (tcf24, renamed ctcf24y), showed testis-specific/biased gene expression, whereas their counterparts cephx1x and ctf24x, located in female chromosome 18, were similarly expressed in both sexes. In addition, male-specific PCR targeting the cephx1y gene revealed that this genomic feature is conserved in cobia populations from Panama, Brazil, Australia, and Japan. The first comprehensive genomic survey presented here is a valuable resource for future studies on cobia population structure and dynamics, conservation, and evolutionary history. Furthermore, it establishes evidence of putative male heterogametic regions with 2 genes playing a potential role in the sex determination of the species, and it provides further support for the rapid evolution of sex-determining mechanisms in teleost fish.

High-quality assembly of the T2T genome for Isodon rubescens f. lushanensis reveals genomic structure variations between 2 typical forms of Isodon rubescens.

Rabdosiae rubescentis herba (Isodon rubescens) is widely used as a folk medicine to treat esophageal cancer and sore throat in China. Its germplasm resources are abundant in China, with I. rubescens (Hemsl.) Hara and I. rubescens f. lushanensis as 2 typical forms. I. rubescens (Hemsl.) Hara is featured by biosynthesis of the diterpenoid oridonin with strong anticancer activity, while I. rubescens f. lushanensis produces another diterpenoid with anticancer activity, lushanrubescensin. However, the biosynthetic pathways of both still need to be fully understood. In particular, little is known about the genetic background of I. rubescens f. lushanensis. We used Pacific Biosciences (PacBio) single-molecule real-time and Nanopore Ultra-long sequencing platforms, respectively, and obtained 139.07 Gb of high-quality data, with a sequencing depth of about 328×. We also obtained a high-quality reference genome for I. rubescens f. lushanensis, with a genome size of 349 Mb and a contig N50 of 28.8 Mb. The heterozygosity of the genome is 1.7% and the repeatability is 83.43%. In total, 34,865 protein-coding genes were predicted. Moreover, we found that most of the variant or unique genes in the diterpenoid synthesis pathways of I. rubescens f. lushanensis and I. rubescens (Hemsl.) Hara were enriched in diterpene synthases. We provide the first genome sequence and gene annotation for the I. rubescens f. lushanensis, which provides molecular evidence for understanding the chemotypic differences of I. rubescens.

Spatial-temporal characterization of photoemission in a streak-mode dynamic transmission electron microscope.

A long-standing motivation driving high-speed electron microscopy development is to capture phase transformations and material dynamics in real time with high spatial and temporal resolution. Current dynamic transmission electron microscopes (DTEMs) are limited to nanosecond temporal resolution and the ability to capture only a few frames of a transient event. With the motivation to overcome these limitations, we present our progress in developing a streak-mode DTEM (SM-DTEM) and demonstrate the recovery of picosecond images with high frame sequence depth. We first demonstrate that a zero-dimensional (0D) SM-DTEM can provide temporal information on any local region of interest with a 0.37 μm diameter, a 20-GHz sampling rate, and 1200 data points in the recorded trace. We use this method to characterize the temporal profile of the photoemitted electron pulse, finding that it deviates from the incident ultraviolet laser pulse and contains an unexpected peak near its onset. Then, we demonstrate a two-dimensional (2D) SM-DTEM, which uses compressed-sensing-based tomographic imaging to recover a full spatiotemporal photoemission profile over a 1.85-μm-diameter field of view with nanoscale spatial resolution, 370-ps inter-frame interval, and 140-frame sequence depth in a 50-ns time window. Finally, a perspective is given on the instrumental modifications necessary to further develop this promising technique with the goal of decreasing the time to capture a 2D SM-DTEM dataset.

CD-Loop: a chromatin loop detection method based on the diffusion model.

In recent years, there have been significant advances in various chromatin conformation capture techniques, and annotating the topological structure from Hi-C contact maps has become crucial for studying the three-dimensional structure of chromosomes. However, the structure and function of chromatin loops are highly dynamic and diverse, influenced by multiple factors. Therefore, obtaining the three-dimensional structure of the genome remains a challenging task. Among many chromatin loop prediction methods, it is difficult to fully extract features from the contact map and make accurate predictions at low sequencing depths. In this study, we put forward a deep learning framework based on the diffusion model called CD-Loop for predicting accurate chromatin loops. First, by pre-training the input data, we obtain prior probabilities for predicting the classification of the Hi-C contact map. Then, by combining the denoising process based on the diffusion model and the prior probability obtained by pre-training, candidate loops were predicted from the input Hi-C contact map. Finally, CD-Loop uses a density-based clustering algorithm to cluster the candidate chromatin loops and predict the final chromatin loops. We compared CD-Loop with the currently popular methods, such as Peakachu, Chromosight, and Mustache, and found that in different cell types, species, and sequencing depths, CD-Loop outperforms other methods in loop annotation. We conclude that CD-Loop can accurately predict chromatin loops and reveal cell-type specificity. The code is available at https://github.com/wangyang199897/CD-Loop.

Analysis of the primary factors influencing donor derived cell-free DNA testing in kidney transplantation.

The donor-derived cell-free DNA (ddcfDNA) is found in the plasma and urine of kidney transplant recipients and displays notable potential in diagnosing rejection, specifically antibody-mediated rejection (ABMR). Nonetheless, the quantitative methods of ddcfDNA lacking standardization and diverse detection techniques can impact the test outcomes. Besides, both the fraction and absolute values of ddcfDNA have been reported as valuable markers for rejection diagnosis, but they carry distinct meanings and are special in various pathological conditions. Additionally, ddcfDNA is highly sensitive to kidney transplant injury. The various sampling times and combination with other diseases can indeed impact ddcfDNA detection values. This review comprehensively analyses the various factors affecting ddcfDNA detection in kidney transplantation, including the number of SNPs and sequencing depths. Furthermore, different pathological conditions, distinct sampling time points, and the presence of complex heterologous signals can influence ddcfDNA testing results in kidney transplantation. The review also provides insights into ddcfDNA testing on different platforms along with key considerations.

Synchrotron radiation X-ray diffraction computed tomography (XRDCT): a new tool in cultural heritage and stone conservation for 3D non-destructive probing and phase analysis of inorganic re-treatments.

The issue of preserving carbonatic stones of cultural heritage (CH) restored in the past that have undergone new decay phenomena is strongly emerging and conservation science has not yet found a reliable solution. In this paper, we propose the application of synchrotron radiation X-ray diffraction computed tomography (XRDCT) to explore the effects of using inorganic-mineral products (ammonium oxalate; ammonium phosphate) in sequence as a novel, compatible and effective re-treatment approach to consolidate decayed carbonatic stones already treated with inorganic-mineral treatments. High-quality XRDCT datasets were used to qualitatively/quantitatively investigate and 3D localize the complex mixture of crystalline phases formed after the conservation re-treatments within a porous carbonatic stone substrate. The XRDCT reconstruction images and the structural refinements of XRD patterns with the Rietveld methods showed that the phase composition of reaction products, their volume distribution, and weight fraction vary as a function of the treatment sequence and penetration depth. The high potential of XRDCT allows (i) assessment of peculiar trends of each treatment/treatment sequence; (ii) exploration of the reaction steps of the sequential treatments and (iii) demonstration of the consolidating effect of inorganic re-treatments, non-destructively and at the micron scale. Above all, our study (i) provides new analytical tools to support the conservation choices, (ii) showcases new analytical possibilities for XRDCT in conservation science, including in investigations of CH materials and decay processes, and (iii) opens up new perspectives in analytical chemistry and material characterisation for the non-destructive and non-invasive analysis of reactions within heterogeneous polycrystalline systems.

Variant calling and genotyping accuracy of ddRAD-seq: Comparison with 20X WGS in layers.

Whole Genome Sequencing (WGS) remains a costly or unsuitable method for routine genotyping of laying hens. Until now, breeding companies have been using or developing SNP chips. Nevertheless, alternatives methods based on sequencing have been developed. Among these, reduced representation sequencing approaches can offer sequencing quality and cost-effectiveness by reducing the genomic regions covered by sequencing. The aim of this study was to evaluate the ability of double digested Restriction site Associated DNA sequencing (ddRAD-seq) to identify and genotype SNPs in laying hens, by comparison with a presumed reliable WGS approach. Firstly, the sensitivity and precision of variant calling and the genotyping reliability of ddRADseq were determined. Next, the SNP Call Rate (CRSNP) and mean depth of sequencing per SNP (DPSNP) were compared between both methods. Finally, the effect of multiple combinations of thresholds for these parameters on genotyping reliability and amount of remaining SNPs in ddRAD-seq was studied. In raw form, the ddRAD-seq identified 349,497 SNPs evenly distributed on the genome with a CRSNP of 0.55, a DPSNP of 11X and a mean genotyping reliability rate per SNP of 80%. Considering genomic regions covered by expected enzymatic fragments (EFs), the sensitivity of the ddRAD-seq was estimated at 32.4% and its precision at 96.4%. The low CRSNP and DPSNP values were explained by the detection of SNPs outside the EFs theoretically generated by the ddRAD-seq protocol. Indeed, SNPs outside the EFs had significantly lower CRSNP (0.25) and DPSNP (1X) values than SNPs within the EFs (0.7 and 17X, resp.). The study demonstrated the relationship between CRSNP, DPSNP, genotyping reliability and the number of SNPs retained, to provide a decision-support tool for defining filtration thresholds. Severe quality control over ddRAD-seq data allowed to retain a minimum of 40% of the SNPs with a CcR of 98%. Then, ddRAD-seq was defined as a suitable method for variant calling and genotyping in layers.

Deep RNA sequencing of muscle tissue reveals absence of viral signatures in dermatomyositis.

Objective: To explore a possible connection between active viral infections and manifestation of dermatomyositis (DM). Methods: Skeletal muscle biopsies were analyzed from patients diagnosed with juvenile (n=10) and adult (n=12) DM. Adult DM patients harbored autoantibodies against either TIF-1γ (n=7) or MDA5 (n=5). Additionally, we investigated skeletal muscle biopsies from non-diseased controls (NDC, n=5). We used an unbiased high-throughput RNA sequencing (HTS) approach to detect viral sequences. To further increase sequencing depth, a host depletion approach was applied. Results: In this observational study, no relevant viral sequences were detected either by native sequencing or after host depletion. The absence of detectable viral sequences makes an active viral infection of the muscle tissue unlikely to be the cause of DM in our cohorts. Discussion: Type I interferons (IFN) play a major role in the pathogenesis of both juvenile and adult DM. The IFN response is remarkably conserved between DM subtypes classified by specific autoantibodies. Certain acute viral infections are accompanied by a prominent type I IFN response involving similar downstream mechanisms as in DM. Aiming to elucidate the pathogenesis of DM in skeletal muscle tissue, we used deep RNA sequencing and a host depletion approach to detect possible causative viruses.

High-Throughput Mapping of Chromosomal Conformations in E. coli Under Physiological Conditions Using Massively Multiplexed Mu Transposition.

Many approaches for measuring three-dimensional chromosomal conformations rely upon formaldehyde crosslinking followed by subsequent proximity ligation, a family of methods exemplified by 3C, Hi-C, etc. Here we provide an alternative crosslinking-free procedure for high-throughput identification of long-range contacts in the chromosomes of enterobacteria, making use of contact-dependent transposition of phage Mu to identify distant loci in close contact. The procedure described here will suffice to provide a comprehensive map of transposition frequencies between tens of thousands of loci in a bacterial genome, with the resolution limited by the diversity of the insertion site library used and the sequencing depth applied.

Deciphering the population structure and genetic basis of growth traits from whole-genome resequencing of the leopard coral grouper ( Plectropomus leopardus).

The leopard coral grouper ( Plectropomus leopardus) is a species of significant economic importance. Although artificial cultivation of P. leopardus has thrived in recent decades, the advancement of selective breeding has been hindered by the lack of comprehensive population genomic data. In this study, we identified over 8.73 million single nucleotide polymorphisms (SNPs) through whole-genome resequencing of 326 individuals spanning six distinct groups. Furthermore, we categorized 226 individuals with high-coverage sequencing depth (≥14×) into eight clusters based on their genetic profiles and phylogenetic relationships. Notably, four of these clusters exhibited pronounced genetic differentiation compared with the other populations. To identify potentially advantageous loci for P. leopardus, we examined genomic regions exhibiting selective sweeps by analyzing the nucleotide diversity ( θπ) and fixation index ( F ST) in these four clusters. Using these high-coverage resequencing data, we successfully constructed the first haplotype reference panel specific to P. leopardus. This achievement holds promise for enabling high-quality, cost-effective imputation methods. Additionally, we combined low-coverage sequencing data with imputation techniques for a genome-wide association study, aiming to identify candidate SNP loci and genes associated with growth traits. A significant concentration of these genes was observed on chromosome 17, which is primarily involved in skeletal muscle and embryonic development and cell proliferation. Notably, our detailed investigation of growth-related SNPs across the eight clusters revealed that cluster 5 harbored the most promising candidate SNPs, showing potential for genetic selective breeding efforts. These findings provide a robust toolkit and valuable insights into the management of germplasm resources and genome-driven breeding initiatives targeting P. leopardus.

CHAMP delivers accurate taxonomic profiles of the prokaryotes, eukaryotes, and bacteriophages in the human microbiome.

Accurate taxonomic profiling of the human microbiome composition is crucial for linking microbial species to health outcomes. Therefore, we created the Clinical Microbiomics Human Microbiome Profiler (CHAMP), a comprehensive tool designed for the profiling of prokaryotes, eukaryotes, and viruses across all body sites. CHAMP uses a reference database derived from 30,382 human microbiome samples, covering 6,567 prokaryotic and 244 eukaryotic species, as well as 64,003 viruses. We benchmarked CHAMP against established profiling tools (MetaPhlAn 4, Bracken 2, mOTUs 3, and Phanta) using a diverse set of in silico metagenomes and DNA mock communities. CHAMP demonstrated unparalleled species recall, F1 score, and significantly reduced false positives compared to all other tools benchmarked. The false positive relative abundance (FPRA) for CHAMP was, on average, 50-fold lower than the second-best performing profiler. CHAMP also proved to be more robust than other tools at low sequencing depths, highlighting its application for low biomass samples. Taken together, this establishes CHAMP as a best-in-class human microbiome profiler of prokaryotes, eukaryotes, and viruses in diverse and complex communities across low and high biomass samples. CHAMP profiling is offered as a service by Clinical Microbiomics A/S and is available for a fee at https://cosmosidhub.com.

Retrieving complete plastid genomes of endangered Guibourtia timber using hybridization capture for forensic identification and phylogenetic analysis

The high economic value and increased demand for timber have led to illegal logging and overexploitation, threatening wild populations. In this context, there is an urgent need to develop effective and accurate forensic tools for identifying endangered Guibourtia timber species to protect forest ecosystem resources and regulate their trade. In this study, a hybridization capture method was developed and applied to explore the feasibility of retrieving complete plastid genomes from Guibourtia sapwood and heartwood specimens stored in a xylarium (wood collection). We then carried out forensic identification and phylogenetic analyses of Guibourtia within the subfamily Detarioideae. This study is the first to successfully retrieve high-quality plastid genomes from xylarium specimens, with 76.95–99.97% coverage. The enrichment efficiency, sequence depth, and coverage of plastid genomes from sapwood were 16.73 times, 70.47 times and 1.14 times higher, respectively, than those from heartwood. Although the DNA capture efficiency of heartwood was lower than that of sapwood, the hybridization capture method used in this study is still suitable for heartwood DNA analysis. Based on the complete plastid genome, we identified six endangered or commonly traded Guibourtia woods at the species level. This technique also has the potential for geographical traceability, especially for Guibourtia demeusei and Guibourtia ehie. Meanwhile, Bayesian phylogenetic analysis suggested that these six Guibourtia species diverged from closely related species within the subfamily Detarioideae ca. 18 Ma during the Miocene. The DNA reference database established based on the xylarium specimens provides admissible evidence for diversity conservation and evolutionary analyses of endangered Guibourtia species.

QTL Mapping and Genetic Map for the Ornamental Sunflower in China

Quantitative trait locus (QTL) mapping and genetic map are of great significance for ornamental sunflowers in China. In this study, a total of 956.50 Mbp data were obtained, the average Q30 was 93.76%, the average GC content was 42.43%, and the GC distribution of the parents and F2 population of the ornamental sunflowers was normal. At the same time, the double-end comparison efficiency of control data was 90.28%, and the enzyme digestion efficiency was 92.01%. The Specific-Locus Amplified Fragment (SLAF) library construction was normal. Furthermore, a total of 734,893 SLAF markers were obtained, among which 127,855 were polymorphic SLAF markers and 38,908 could be used for genetic map construction, and the effective polymorphism of the parents was 5.29%. Moreover, we constructed a total of 17 linkage groups, with 6181 markers in the QTL mapping, the total map distance was 2608.66 cM, the marker integrity in the figure above was 99%, the proportion of double exchange was 0.05, the sequencing depth of the parents was 42.455 × , and the progeny was 9.24 × . The relationship of traits (plant height, stem diameter, disk diameter, number of petals, leaf number, stigma color, petal color, petiole color) and QTL mapping was closely related to show the best of ornamental effect.

ScDMV: a zero-one inflated beta mixture model for DNA methylation variability with scBS-seq data.

The utilization of single-cell bisulfite sequencing (scBS-seq) methods allows for precise analysis of DNA methylation patterns at the individual cell level, enabling the identification of rare populations, revealing cell-specific epigenetic changes, and improving differential methylation analysis. Nonetheless, the presence of sparse data and an overabundance of zeros and ones, attributed to limited sequencing depth and coverage, frequently results in reduced precision accuracy during the process of differential methylation detection using scBS-seq. Consequently, there is a pressing demand for an innovative differential methylation analysis approach that effectively tackles these data characteristics and enhances recognition accuracy. We propose a novel beta mixture approach called scDMV for analyzing methylation differences in single-cell bisulfite sequencing data, which effectively handles excess zeros and ones and accommodates low-input sequencing. Our extensive simulation studies demonstrate that the scDMV approach outperforms several alternative methods in terms of sensitivity, precision, and controlling the false positive rate. Moreover, in real data applications, we observe that scDMV exhibits higher precision and sensitivity in identifying differentially methylated regions, even with low-input samples. In addition, scDMV reveals important information for GO enrichment analysis with single-cell whole-genome sequencing data that are often overlooked by other methods. The scDMV method, along with a comprehensive tutorial, can be accessed as an R package on the following GitHub repository: https://github.com/PLX-m/scDMV.

Development and evaluation of a haplotype reference panel of Zhikong scallop (Chlamys farreri) for genotype imputation

Single nucleotide polymorphisms (SNPs) across the genome are remarkably useful resources for investigating the genome architecture, studying the genetic basis of complex traits and implementation of genomic selection. However, information on high-density SNPs is unavailable for most aquaculture species owing to the high costs of whole genome sequencing. Genotype imputation is an attractive approach for obtaining large number of genomic variants at a low cost. In this study, the first haplotype reference panel was constructed for Chlamys farreri, an economically important mollusk species and a major contributor to aquaculture in China. A total of 174 scallops from six populations in the main scallop culture area were selected for whole-genome resequencing. More than 87.47 million SNPs were identified and 3,968,417 SNPs genome-wide were selected for constructing the haplotype reference panel for C. farreri. The number of SNPs per chromosome ranged from 125,000 to 279,981 with an average spacing of 242 bases. Several factors were evaluated based on the imputation performance of the reference panel. The imputation accuracies of 19 chromosomes ranged from 0.91 to 0.95, and there were no significant differences among them (p > 0.05). The imputation accuracies of populations could reach 0.95, 0.93, and 0.90 for ZY, PR, and CS respectively, which confirmed a positive relationship between the kinship within the population and the imputation accuracy. The imputation accuracies increased with sequencing depth and plateaued at 0.5×, reaching accuracies of 0.91–0.93 for artificially breeding populations. Taken together, we constructed a haplotype reference panel for C. farreri, which can provide a cost-effective way to extend the limited number of genotypes for boosting the power of downstream analyses such as genome-wide association studies and genomic selection.