ABSTRACT
According to various estimates, only a small percentage of existing viruses have been discovered, naturally much less being represented in the genomic databases. High-throughput sequencing technologies develop rapidly, empowering large-scale screening of various biological samples for the presence of pathogen-associated nucleotide sequences, but many organisms are yet to be attributed specific loci for identification. This problem particularly impedes viral screening, due to vast heterogeneity in viral genomes. In this paper, we present a new bioinformatic pipeline, VirIdAl, for detecting and identifying viral pathogens in sequencing data. We also demonstrate the utility of the new software by applying it to viral screening of the feces of bats collected in the Moscow region, which revealed a significant variety of viruses associated with bats, insects, plants, and protozoa. The presence of alpha and beta coronavirus reads, including the MERS-like bat virus, deserves a special mention, as it once again indicates that bats are indeed reservoirs for many viral pathogens. In addition, it was shown that alignment-based methods were unable to identify the taxon for a large proportion of reads, and we additionally applied other approaches, showing that they can further reveal the presence of viral agents in sequencing data. However, the incompleteness of viral databases remains a significant problem in the studies of viral diversity, and therefore necessitates the use of combined approaches, including those based on machine learning methods.
Subject(s)
Alphacoronavirus/isolation & purification , Betacoronavirus/isolation & purification , Chiroptera/virology , Genome, Viral/genetics , Metagenome/genetics , Alphacoronavirus/classification , Alphacoronavirus/genetics , Animals , Betacoronavirus/classification , Betacoronavirus/genetics , Chiroptera/genetics , Computational Biology/methods , Feces/virology , High-Throughput Nucleotide Sequencing , Metagenomics/methods , Moscow , Phycodnaviridae/classification , Phycodnaviridae/genetics , Phycodnaviridae/isolation & purification , Sequence Analysis, DNAABSTRACT
Kidney stone disease is an urgent medical and social problem. Genetic factors play an important role in the disease development. This study aims to establish an association between polymorphisms in genes coding for proteins involved in calcium metabolism and the development of calcium urolithiasis in Russian population. In this case-control study, we investigated 50 patients with calcium urolithiasis (experimental group) and 50 persons lacking signs of kidney stone disease (control group). For molecular genetic analysis we used a previously developed gene panel consisting of 33 polymorphisms in 15 genes involved in calcium metabolism: VDR, CASR, CALCR, OPN, MGP, PLAU, AQP1, DGKH, SLC34A1, CLDN14, TRPV6, KLOTHO, ORAI1, ALPL, and RGS14. High-throughput target sequencing was utilized to study the loci of interest. Odds ratios and 95% confidence intervals were used to estimate the association between each SNP and risk of urolithiasis development. Multifactor dimensionality reduction analysis was also carried out to analyze the gene-gene interaction. We found statistically significant (unadjusted p-value < 0.05) associations between calcium urolithiasis and the polymorphisms in the following genes: CASR rs1042636 (OR = 3.18 for allele A), CALCR rs1801197 (OR = 6.84 for allele A), and ORAI1 rs6486795 (OR = 2.25 for allele C). The maximum OR was shown for AA genotypes in loci rs1042636 (CASR) and rs1801197 (CALCR) (OR = 4.71, OR = 11.8, respectively). After adjustment by Benjamini-Hochberg FDR we found only CALCR (rs1801197) was significantly associated with the risk of calcium urolithiasis development. There was no relationship between recurrent course of the disease and family history of urolithiasis in investigated patients. Thus we found a statistically significant association of polymorphism rs1801197 (gene CALCR) with calcium urolithiasis in Russian population.
