Detection Of Single Nucleotide Polymorphisms Research Articles

Simultaneous selection of multiple important single nucleotide polymorphisms in familial genome wide association studies data

We propose a resampling-based fast variable selection technique for detecting relevant single nucleotide polymorphisms (SNP) in a multi-marker mixed effect model. Due to computational complexity, current practice primarily involves testing the effect of one SNP at a time, commonly termed as ‘single SNP association analysis’. Joint modeling of genetic variants within a gene or pathway may have better power to detect associated genetic variants, especially the ones with weak effects. In this paper, we propose a computationally efficient model selection approach—based on the e-values framework—for single SNP detection in families while utilizing information on multiple SNPs simultaneously. To overcome computational bottleneck of traditional model selection methods, our method trains one single model, and utilizes a fast and scalable bootstrap procedure. We illustrate through numerical studies that our proposed method is more effective in detecting SNPs associated with a trait than either single-marker analysis using family data or model selection methods that ignore the familial dependency structure. Further, we perform gene-level analysis in Minnesota Center for Twin and Family Research (MCTFR) dataset using our method to detect several SNPs using this that have been implicated to be associated with alcohol consumption.

Length-Dependent Gold Nanoparticle-Labeled Nucleic Acid Probe Based Dual Test Strip for Visualized Discrimination of Single Nucleotide Polymorphism

Single nucleotide polymorphism (SNP) is closely related to the susceptibility of various common diseases and adverse drug reactions. Conclusively, SNP detection is significant for early-stage disease diagnosis and personalized medicine. The gold nanoparticle nucleic acid strip has the advantages of low-cost and rapid visualization (within 5–10 min), and it was introduced into detecting cancer-related SNP sites (A/A, A/G, G/G) in the CYP1A1 gene. As an efficient visualized approach based on the observation of the test zone (TZ) on the strip, it is improper to develop color at the TZ toward the negative sample; otherwise, it may lead to invalid false positive results. Thus, it is crucial to avert such a situation. In this work, dual strips named A and G were prepared, which were completely complementary to the T/T site (complementary base of A/A site) and C/C site (complementary base of G/G site), respectively. Single-strand DNA (ssDNA) containing SNP sites was obtained by conventional two-step polymerase chain reaction (PCR). Then SNP sites could be directly visualized and accurately distinguished on the dual strips. Furthermore, the false positive result could be eliminated by optimizing the length of DNA modified on gold nanoparticle and the concentration of NaCl. The result for the detection of 21 human blood samples with this method was consistent with those of sequencing. The platform exhibits the benefits of rapidity, easily visual determination, is false-positive-free, and is low-cost and plays an important part in the initial screening of SNP-related diseases.

Genome-wide association screening and MassARRAY for detection of high-temperature resistance-related SNPs and genes in a hybrid abalone (Haliotis discus hannai ♀ × H. fulgens ♂) based on super genotyping-by-sequencing

High seawater temperatures in summer often lead to high mortality in farmed abalone, which greatly increases the cost of farming. The hybrid abalone (Haliotis discus hannai ♀ × H. fulgens ♂, DF), is one of the main varieties of abalone cultured in Fujian Province, China, to resist the high temperature. In this study, the high-temperature resistance of DF abalones was assessed under 28 °C - 32 °C. 31 °C was applied in the subsequent screening experiments of high-temperature resistant individuals. A total of 614,321 single nucleotide polymorphisms (SNPs) were identified via super genotyping-by-sequencing (S-GBS), and 69 SNPs involved in 42 candidate genes were identified by genome-wide association study (GWAS). In addition, according to gene functional annotation, the locus with the lowest P value was screened within the threshold range (−log10P = 4.98296). The genotyping was verified by Sequenom MassARRAY, and 7 SNPs were proved to be related to high temperature resistance in DF population of 100 individuals. 7 candidate genes (fabp, e3rglg4, rgs2, tanc2x6, osbp1, arsb, and uch1) were identified near these SNPs. Finally, these candidate genes were further proved to be associated with the high-temperature resistance trait of DF by the analysis of differential expression genes between sensitive and tolerant individuals. These candidate genes involved regulatory mechanisms such as cell stability, protein repair, lipid metabolism and immune regulation. Notably, the mechanism by which the e3rglg4 gene was involved in the degradation of abnormally folded proteins of abalone under heat stress, corroborating the reliability of the results of GWAS and Sequenom MassARRA. These results complement the shortage of molecular markers for heat resistance of hybrid abalone, provide a reference for marker-assisted selection (MAS) and accelerate the development of heat-resistant abalone varieties.

Combining GS-assisted GWAS and transcriptome analysis to mine candidate genes for nitrogen utilization efficiency in Populus cathayana

BackgroundForest trees such as poplar, shrub willow, et al. are essential natural resources for sustainable and renewable energy production, and their wood can reduce dependence on fossil fuels and reduce environmental pollution. However, the productivity of forest trees is often limited by the availability of nitrogen (N), improving nitrogen use efficiency (NUE) is an important way to address it. Currently, NUE genetic resources are scarce in forest tree research, and more genetic resources are urgently needed.ResultsHere, we performed genome-wide association studies (GWAS) using the mixed linear model (MLM) to identify genetic loci regulating growth traits in Populus cathayana at two N levels, and attempted to enhance the signal strength of single nucleotide polymorphism (SNP) detection by performing genome selection (GS) assistance GWAS. The results of the two GWAS analyses identified 55 and 40 SNPs that were respectively associated with plant height (PH) and ground diameter (GD), and 92 and 69 candidate genes, including 30 overlapping genes. The prediction accuracy of the GS model (rrBLUP) for phenotype exceeds 0.9. Transcriptome analysis of 13 genotypes under two N levels showed that genes related to carbon and N metabolism, amino acid metabolism, energy metabolism, and signal transduction were differentially expressed in the xylem of P. cathayana under N treatment. Furthermore, we observed strong regional patterns in gene expression levels of P. cathayana, with significant differences between different regions. Among them, P. cathayana in Longquan region exhibited the highest response to N. Finally, through weighted gene co-expression network analysis (WGCNA), we identified a module closely related to the N metabolic process and eight hub genes.ConclusionsIntegrating the GWAS, RNA-seq and WGCNA data, we ultimately identified four key regulatory genes (PtrNAC123, PtrNAC025, Potri.002G233100, and Potri.006G236200) involved in the wood formation process, and they may affect P. cathayana growth and wood formation by regulating nitrogen metabolism. This study will provide strong evidence for N regulation mechanisms, and reliable genetic resources for growth and NUE genetic improvement in poplar.

Comprehensive analysis of clinical indications and viral strain variants among patients infected with SARS-CoV-2 in Inner Mongolia, China

Since the first appearance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in 2019, the virus is still evolving and mutating until now. In this study, we collected 6 throat swabs from patients who diagnosed with COVID-19 in Inner Mongolia, China, to understand the entry of multiple SARS-CoV-2 variants into Inner Mongolia and analyze the relationships between variants and clinical features observed in infected patients. In addition, we performed a combined analysis of clinical parameters associated with SARS-CoV-2 variants of interest, pedigree analysis, and detection of single-nucleotide polymorphisms. Our results showed that the clinical symptoms were generally mild although some patients demonstrated some degree of liver function abnormalities, and the SARS-CoV-2 strain was related to the Delta variant (B.1.617.2), AY.122 lineage. The epidemiological investigations and clinical manifestations confirmed that the variant exhibits strong transmission, a high viral load, and moderate clinical symptoms. SARS-CoV-2 has undergone extensive mutations in various hosts and countries. Timely monitoring of virus mutation can help to monitor the spread of infection and characterize the diversity of genomic variants, thus limiting future waves of SARS-CoV-2 infection.

Association of the rs2111234, rs3135499, rs8057341 polymorphisms in the NOD2 gene with leprosy: A case-control study in the Norte de Santander, Colombia population.

Leprosy is a chronic infectious disease caused by Mycobacterium leprae. The development of leprosy involves several factors, including the causative agent, the individual host's immune response, environmental factors, and the genetic background of the host. Specifically, the host's innate immune response, encoded by genes, determines their susceptibility to developing leprosy post-infection. Polymorphic variants in the nucleotide-binding oligomerization domain 2 (NOD2) gene are associated with leprosy among populations in a variety of endemic areas around the world. Colombia, a country located in the tropical zone, has several leprosy-endemic regions, including Norte de Santander. The aim of this study was to analyze the rs7194886, rs2111234, rs3135499, and rs8057341 single nucleotide polymorphisms (SNPs) in the NOD2 gene using a case-control study to determine whether they confer greater or lesser susceptibility to the development of leprosy. The TaqMan qPCR amplification system was used for SNPs detection. An association between the A-rs8057341 SNP (p = 0,006286) and resistance to leprosy was found. However, the rs3135499 (p = 0,9063) and rs2111234 (p = 0.1492) were not found to be associated with leprosy susceptibility. In addition, the rs7194886 SNP was not found to be in Hardy-Weinberg equilibrium (HWE) in the study population. The GAG haplotype, consisting of SNPs rs2111234-G, rs3135499-A, and rs8057341G, acts as a susceptibility factor for the development of leprosy in women. SNPs rs3135499 and rs8057341 are functionally related to decreased NOD2 expression according to an in-silico analysis. The SNPs rs8057341-A was related with resistance to leprosy and the haplotype rs2111234-G, rs3135499-A and rs8057341-G SNPs was related with susceptibility in the Norte de Santander Colombia, studied population.

Single-Nucleotide Polymorphisms, c.415C > T (Arg139Cys) and c.416G > A (Arg139His), in the NUDT15 Gene Are Associated with Thiopurine-Induced Leukopenia.

While nucleoside diphosphate-linked moiety X-type motif 15 (NUDT15) gene polymorphism Arg139Cys (rs116855232) is known to be a risk factor for thiopurine-induced severe leukopenia, association with the NUDT15 gene polymorphism Arg139His (rs147390019) has not yet been clarified. In addition, the accuracy of TaqMan PCR to assess these two polymorphisms has not been investigated. In this study, we evaluated TaqMan PCR for detection of the NUDT15 single-nucleotide polymorphisms (SNPs) and examined the clinical impact of Arg139His on thiopurine-induced leukopenia. First, we demonstrated that a TaqMan PCR assay successfully detected the Arg139His polymorphism of NUDT15 in clinical samples. Next, the NUDT15 gene polymorphisms (Arg139Cys and Arg139His) were separately analyzed by TaqMan Real-Time PCR in 189 patients from August 2018 to July 2019. The incidences of leukopenia within 2 years were 16.2, 57.9, and 100% for arginine (Arg)/Arg, Arg/cysteine (Cys), and Arg/histidine (His), respectively. The leukopenia was significantly increased in Arg/Cys and Arg/His compared with Arg/Arg. This retrospective clinical study indicated that, in addition to Arg139Cys, Arg139His may be clinically associated with a high risk of leukopenia. Pharmacogenomics will help in selecting drugs and determining the individualized dosage of thiopurine drugs.

Detection of single nucleotide polymorphism in HTR4 and its relationship with growth traits in sheep

In this study, a single nucleotide polymorphism of HTR4 (hydroxytryptamine receptor 4) was detected using DNA sequencing and KASPar (Kompetitive Allele-Specific PCR) technique with the aim of analyzing its effect on growth traits in 1102 sheep. A synonymous mutation g.101220 C > T located on the fifth intron of the ovis HTR4 gene was detected, and association analysis showed that this mutation was significantly associated with growth traits in sheep (p ＜.05), with TT genotypes having significantly lower body weight, height, length and chest circumference than TC and CC genotypes. It showed that the polymorphism of this locus was significantly associated with growth traits in sheep. In addition, qRT-PCR results showed that HTR4 was expressed in different tissues of sheep. It is highly expressed in the liver, spleen and duodenum. As important metabolic, immune and digestive absorption organs in animals, the above tissues can regulate the excitability of intestinal smooth muscle by participating in the body metabolism and nutrient metabolism of sheep, so that sheep can show better growth characteristics. In conclusion, the polymorphic locus identified in HTR4 gene can be used as candidate molecular marker in sheep breeding.

Detecting Melanocortin 1 Receptor Gene's SNPs by CRISPR/enAsCas12a.

Beyond its powerful genome-editing capabilities, the CRISPR/Cas system has opened up a new era of molecular diagnostics due to its highly specific base recognition and trans-cleavage activity. However, most CRISPR/Cas detection systems are mainly used to detect nucleic acids of bacteria or viruses, while the application of single nucleotide polymorphism (SNP) detection is limited. The MC1R SNPs were investigated by CRISPR/enAsCas12a and are not limited to the protospacer adjacent motif (PAM) sequence in vitro. Specifically, we optimized the reaction conditions, which proved that the enAsCas12a has a preference for divalent magnesium ion (Mg2+) and can effectively distinguish the genes with a single base difference in the presence of Mg2+, and the Melanocortin l receptor (MC1R) gene with three kinds of SNP sites (T305C, T363C, and G727A) was quantitatively detected. Since the enAsCas12a is not limited by PAM sequence in vitro, the method shown here can extend this extraordinary CRISPR/enAsCas12a detection system to other SNP targets, thus providing a general SNP detection toolbox.

Polymorphisms and gene expression of metalloproteinases and their inhibitors associated with cerebral ischemic stroke in young patients with sickle cell anemia.

Sickle cell anemia (SCA) is a genetic disease with great clinical heterogeneity and few viable strategies for treatment; hydroxyurea (HU) is the only widely used drug. Thus, the study of single nucleotide polymorphisms (SNPs) and the gene expression of MMPs 1, 2, 9, 7 and TIMPs 1 and 2, which are involved in the regulation of extracellular matrix, inflammation, and neuropathies, may provide further insights into the pathophysiology of the disease and elucidate biomarkers and molecules as potential therapeutic targets for patients with SCA. We evaluated 251 young individuals with SCA from northeastern Brazil. The groups were divided according to vaso-occlusive crisis (VOC) and cerebrovascular disease (CVD), compared to control individuals. SNP detection and gene expression assays were performed by real-time PCR, TaqMan system®. Both the expression levels of MMP1 gene, and the SNP MMP1-1607 1G/2G were associated with the risk of cerebral ischemic stroke (IS), and the expression of MMP1 was also associated with a higher frequency of VOC/year. Expression levels of MMP7, TIMP1, and TIMP2 were increased in patients conditioned to IS. The SNP 372T>C (rs4898) TIMP1 T alleles were more frequent in patients with > 5 VOC events/year. The SNP rs17576 of MMP9 showed differences in gene expression levels; it was increased in the genotypes AG, and AG+GG. The findings of this study, the SNPs, and expression provide initial support for understanding the role of MMPs-TIMPs in the pathophysiology of SCA in young patients.

SARS-CoV-2 genotyping and sequencing following a simple and economical RNA extraction and storage protocol.

Since the beginning of the SARS-CoV-2 pandemic, supply chain shortages have caused major disruptions in sourcing the materials needed for laboratory-based molecular assays. With increasing demand for molecular testing, these disruptions have limited testing capacity and hindered efforts to mitigate spread of the virus and new variants. Here we evaluate an economical and reliable protocol for the extraction and short-term ambient temperature storage of SARS-CoV-2 RNA. Additional objectives of the study were to evaluate RNA from this protocol for 1) detection of single nucleotide polymorphisms (SNPs) in the spike gene and 2) whole genome sequencing of SARS-CoV-2. The RNAES protocol was evaluated with residual nasopharyngeal (NP) samples collected from Emory Healthcare and Emory Student Health services. All RNAES extractions were performed in duplicate and once with a commercial extraction robot for comparison. Following extraction, eluates were immediately tested by rRT-PCR. SARS-CoV-2 RNA was successfully detected in 56/60 (93.3%) RNAES replicates, and Ct values corresponded with comparator results. Upon testing in spike SNP assays, three genotypes were identified, and all variant calls were consistent with those previously obtained after commercial extraction. Additionally, the SARS-RNAES protocol yield eluate pure enough for downstream whole genome sequencing, and results were consistent with SARS-CoV-2 whole genome sequencing of eluates matched for Ct value. With reproducible results across a range of virus concentrations, the SARS-RNAES protocol could help increase SARS-CoV-2 diagnostic testing and monitoring for emerging variants in resource-constrained communities.

Com probe implemented STexS II greatly enhances specificity in SARS-CoV-2 variant detection

The initial introduction of utilizing double helix structural oligonucleotides known as SNP typing with excellent specificity (STexS) in a standard PCR greatly improved the detection of single nucleotide polymorphisms (SNP) by enhancing amplification rates of primer-matching strands and interrupting mismatched strands by constant instability of kinetics regarding alignment attaching and detaching. The model was beneficial overall in detecting SNP variants consisting of large amounts of wildtype strands such as EGFR mutation genotyping for early detection of non-small cell lung cancer. While the STexS PCR is advantageous in detecting SNPs and biomarkers, limitations were yet observed. Despite the ability to detect variants 10 times more effective than a typical amplification-refractory mutation system PCR, it could only perform optimally in DNA concentrations around 101 ~ 105. To further enhance STexS specificity to perform detecting viral-RNA variants such as the infamous SARS-CoV-2, a novel improvement of the regular TaqMan Probe using Com-probes to inhibit high copy wild targets and amplify low copy mutant targets. By introducing the novel STexS II, omicron variants of SARS-CoV-2 were able to be successfully detected in high concentrations of normal genes.

A protein enzyme-free strategy for fluorescence detection of single nucleotide polymorphisms using asymmetric MNAzymes

To establish protein enzyme-free and simple approach for sensitive detection of single nucleotide polymorphisms (SNPs), the nucleic acid amplification reactions were developed to reduce the dependence on protein enzymes (polymerase, endonuclease, ligase). These methods, while enabling highly amplified analysis for the short sequences, cannot be generalized to long genomic sequences. Herein, we develop a protein enzyme-free and general SNPs assay based on asymmetric MNAzyme probes. The multi-arm probe (MNAzyme-9M-13) with two asymmetric recognition arms, containing a short (9 nt) and a long (13 nt) arm, is designed to detect EGFR T790 M mutation (MT). Owing to the excellent selectivity of short recognition arm, MNAzyme-9M-13 probe can efficiently avoid interferences from wild-type target (WT) and various single-base mutations. Through a one-pot mixing, MNAzyme-9M-13 probe enables the sensitive detection of MT, without protein enzyme or multi-step operation. The calculated detection limit for MT is 0.59 nM and 0.83%. Moreover, this asymmetric MNAzyme strategy can be applied for SNPs detection in long genomic sequences as well as short microRNAs (miRNAs) only by changing the low-cost unlabeled recognition arms. Therefore, along with simple operation, low-cost, protein enzyme-free and strong versatility, our asymmetric MNAzyme strategy provides a novel solution for SNPs detection and genes analysis.

Identification of single nucleotide polymorphisms by a peptide nucleic acid-based sandwich hybridization assay coupled with toehold-mediated strand displacement reactions

In this work, we developed a simple and accurate peptide nucleic acid (PNA)-based sandwich hybridization assay for single nucleotide polymorphisms (SNPs) in the p53 gene. Our approach combines the enzyme-free toehold-mediated strand displacement reaction (SDR) with real-time enzyme-linked immunosorbent assay (ELISA) to detect SNPs with high sensitivity and specificity. A PNA-DNA heteroduplex with an external toehold is designed and fixed on well surface of a 96-well plate. The strand displacement from PNA-DNA heteroduplexes is initiated by the hybridization of target sequence with the toehold domain and ends with the fully displacing of the incumbent DNA. Finally, the as formed PNA-target DNA duplex with overhang at its 5′-end hybridizes with a biotin-labeled reporter PNA to form a sandwich structure on surface for signal amplification. The proposed PNA-based sandwich biosensor displays high sensitivity and greatly enhanced discriminability to target p53 gene segments against single-base mutant sequences compared to its all-DNA counterpart. Furthermore, the probe design is elegantly simple and the sensing procedure is easy to operate. We believe that this strategy may provide a simple and universal strategy for SNPs detection through easily altering the sequences of probes according to the sequences around target SNPs.

Bioinformatic method for determining single nucleotide polymorphisms on the example of gene <i>WIN</i> in <i>Glycine max</i>

In this paper, a hypothetical method for locating SNPs (single nucleotide polymorphisms) on the example of the ribonuclease gene WIN was proposed. Ribonuclease comprises an enzyme that participates in defence reactions against fungal infections in soybeans, as well as other protective responses to biotic stress. Its belonging to the RNA-ases group determines the specific properties, namely the ability to degrade foreign nucleic acids. This ability provides for a general nonspecific immune response of the plant to the invasion of antigenic structures. Modern biotechnology calls for the development of molecular methods and approaches that will increase the resistance of a culture or accelerate the processes of its adaptation in the field. This problem can be solved by using technologies of SNP artificial induction in those parts of the genome that encode proteins capable of acting in protective reactions against biotic stress. In the study, 5 single-nucleotide polymorphisms were proposed using bioinformatic analysis. Since the localisation and detection of SNPs comprise a challenging task due to the presence of a single nucleotide change, in the biotechnological practice, predictive analysis is carried out in order to localise the potential sequence of occurring single-nucleotide polymorphism. Following the identification of the hypothetical SNP location, they can be further detected using complex molecular methods, such as real-time PCR or local sequencing. This technology can become a powerful tool for breeding soybean varieties having predetermined properties. Such theoretical and predictive models will allow for a quicker response to the dynamic environment under manmade load on plants.

Universal probe-based SNP genotyping with visual readout: a robust and versatile method.

Detection of single nucleotide polymorphisms (SNPs) is critical for personalized clinical diagnosis, treatment, and medication. Current clinical detection methods suffer from primer dimerization and require the redesigning of reaction systems for different targets, resulting in a time-consuming and laborious process. Here, we present a robust and versatile method for SNP typing by using tailed primers and universal small molecule probes in combination with a visualized lateral flow assay (LFA). This approach enables not only rapid typing of different targets, but also eliminates the interference of primer dimers and enhances the accuracy and reliability of the results. Our proposed universal assay has been successfully applied to the typing of four SNP loci of clinical samples to verify the accuracy and universality, and the results are consistent with those obtained by Sanger sequencing. Therefore, our study establishes a new universal "typing formula" using nucleic acid tags and small molecule probes that provides a powerful genotyping platform for genetic analysis and molecular diagnostics.

High-Throughput Association Mapping in Brassica napus L.: Methods and Applications.

Oil seed rape (Braasica napus L.) is ranked second among oil seed crops cultivated globally for edible oil for human, and seed cake for animal consumption. Recent genetic and genomics advancements highlighted the diversity that exists within B. napus, which is largely discovered using the most promising genetic markers called single nucleotide polymorphism (SNP). Their calling rate is also enhanced to ~100 folds after the continuous advancements in the next generation sequencing (NGS) technologies. As the high throughput of NGS resulted in multi-Giga bases data, the detailed quality control (QC) prior to downstream analyses is a pre-requisite. It mainly involved the removal of false positives, missing proportions, filtering of low-quality SNPs, and adjustments of minor-allele frequency and heterozygosity. After marker-trait association, for conformation of target SNPs, validations of SNPs can be performed using various methods, especially allele-specific PCR assay-based methods have been utilized for SNP genotyping of genes targeting agronomic traits and somaclonal variations occurred during transgenic studies. In the present study, the authors mainly argue on the genotypic progress, and pipelines/methods that are being used for detection, calling, filtering, and validation of SNPs. Also, insight is provided into the application of SNPs in linkage and association mapping, including QTL mapping and genome-wide association studies targeting mainly developmental traits related to the root system and plant architecture, flowering time, silique, and oil quality. Briefly, the present study provides the recent information and recommendations on the SNP genotyping methods and its applications, which can be useful for marker-assisted breeding in B. napus and other crops.

Nucleic acid amplification test (NAAT) conducted in a microfluidic chip to differentiate between various ginseng species.

Panax ginseng and Panax quinquefolius have different medicinal properties and market values; however, they can be difficult to distinguish from one another based on physical appearances alone. Therefore, a molecular test that can be performed in commercial settings is needed to overcome this difficulty. A locus that contains a single nucleotide polymorphism (SNP) site to differentiate between P. ginseng and P. quinquefolius has been selected. An isothermal nucleic acid amplification test (NAAT) has been developed for use in a microfluidic chip; this NAAT method, which is based on lesion-induced DNA amplification (LIDA), amplifies the extracted plant genomic samples and enhances the detection of specific SNPs. This NAAT method was used to authenticate five ginseng root samples which indicated that two of the five samples appear to be mislabeled. These authentication results were consistent with those obtained from next generation sequencing (NGS) although this molecular test is more affordable and faster than NGS.

Modified Allele-Specific qPCR (ASQ) Genotyping.

The allele-specific qPCR (ASQ) method for SNP (single nucleotide polymorphism) detection is based on the FRET (fluorescence resonance energy transfer) system, a system using position-dependent fluorescent dyes and quenches. The modified ASQ method requires two separate components: (1) the allele-specific part, two AS primers targeting the SNP with identity in the penultimate positions atthe 3'-end and specific tags in the 5'-end, and (2) the universal part, two universal probes (UPs) with corresponding tags and different fluorescent dyes in the 5'-end and a single common universal probe with a quencher in the 3'-ends (Uni-Q), complementary to all UP tags. There are two major variations of the ASQ method, with either short 4-bp tags (variant A) or longer 6-bp tags (variant B), both of which have been successfully used for SNP genotyping in plants. The modified ASQ method is much cheaper compared to other similar FRET-based methods because the most expensive parts, the universal probes, have a short and linear structure, where fluorophores and quenchers are located in the ends but not incorporated inside of the sequences.

NanoSNP: a progressive and haplotype-aware SNP caller on low-coverage nanopore sequencing data.

Oxford Nanopore sequencing has great potential and advantages in population-scale studies. Due to the cost of sequencing, the depth of whole-genome sequencing for per individual sample must be small. However, the existing single nucleotide polymorphism (SNP) callers are aimed at high-coverage Nanopore sequencing reads. Detecting the SNP variants on low-coverage Nanopore sequencing data is still a challenging problem. We developed a novel deep learning-based SNP calling method, NanoSNP, to identify the SNP sites (excluding short indels) based on low-coverage Nanopore sequencing reads. In this method, we design a multi-step, multi-scale and haplotype-aware SNP detection pipeline. First, the pileup model in NanoSNP utilizes the naive pileup feature to predict a subset of SNP sites with a Bi-long short-term memory (LSTM) network. These SNP sites are phased and used to divide the low-coverage Nanopore reads into different haplotypes. Finally, the long-range haplotype feature and short-range pileup feature are extracted from each haplotype. The haplotype model combines two features and predicts the genotype for the candidate site using a Bi-LSTM network. To evaluate the performance of NanoSNP, we compared NanoSNP with Clair, Clair3, Pepper-DeepVariant and NanoCaller on the low-coverage (∼16×) Nanopore sequencing reads. We also performed cross-genome testing on six human genomes HG002-HG007, respectively. Comprehensive experiments demonstrate that NanoSNP outperforms Clair, Pepper-DeepVariant and NanoCaller in identifying SNPs on low-coverage Nanopore sequencing data, including the difficult-to-map regions and major histocompatibility complex regions in the human genome. NanoSNP is comparable to Clair3 when the coverage exceeds 16×. https://github.com/huangnengCSU/NanoSNP.git. Supplementary data are available at Bioinformatics online.