Spectrum of PRSS1, SPINK1, CTRC, CFTR, and CPA1 Gene Variants in Chronic Pancreatitis Patients in Russia.

The aim of the study was to define the spectrum of genetic risk factors of chronic pancreatitis (CP) development in patients living in the European part of the Russian Federation. Materials and Methods: The study group included 105 patients with CP, with the age of the disease onset under 40 years old (the average age of onset was 26.9 years). The control group consisted of 76 persons without clinical signs of pancreatitis. The diagnosis of chronic pancreatitis in patients was made on the basis of clinical manifestations and the results of laboratory and instrumental investigations. Genetic examination of patients was conducted using the next-generation sequencing (NGS) technology and included targeted sequencing of all exons and exon-intron boundaries of the PRSS1, SPINK1, CTRC, CFTR, and CPA1 genes. The genotyping of the rs61734659 locus of the PRSS2 gene was also conducted. Results: Genetic risk factors of the CP development were found in 61% of patients. Pathogenic and likely-pathogenic variants associated with the risk of CP development were identified in the following genes: CTRC (37.1% of patients), CFTR (18.1%), SPINK1 (8.6%), PRSS1 (8.6%), and CPA1 (6.7%). The frequent gene variants in Russian patients with CP were as follows: CTRC gene - c.180C>T (rs497078), c.760C>T (rs121909293), c.738_761del24 (rs746224507); cumulative odds ratio (OR) for all risk alleles was 1.848 (95% CI: 1.054-3.243); CFTR gene - c.3485G>T (rs1800120), c.1521_1523delCTT (p.Phe508del, rs113993960), and c.650A>G (rs121909046); OR=2.432 (95% CI: 1.066-5.553). In the SPINK1, PRSS1, and CPA1 genes, pathogenic variants were found only in the group of patients with CP. The frequent variants of the SPINK1 gene include c.101A>G (p.Asn34Ser, rs17107315) and c.194+2T>C (rs148954387); of the PRSS1 gene - c.86A>T (p.Asn29Ile, rs111033566); of the CPA1 gene - c.586-30C>T (rs782335525) and c.696+23_696+24delGG. The OR for the CP development for the c.180TT genotype (rs497078) CTRC according to the recessive model (TT vs. CT+CC) was 7.05 (95% CI: 0.86-263, p=0.011). In the CTRC gene, the variant c.493+49G>C (rs6679763) appeared to be benign, the c.493+51C>A (rs10803384) variant was frequently detected among both the diseased and healthy persons and did not demonstrate a protective effect. The protective factor c.571G>A (p.Gly191Arg, rs61734659) of the PRSS2 gene was detected only in the group of healthy individuals and confirmed its protective role. 12.4% of the patients with CP had risk factors in 2 or 3 genes. Conclusion: Sequencing of the coding regions of the PRSS1, SPINK1, CTRC, CFTR, and CPA1 genes allowed to identify genetic risk factors of the CP development in 61% of cases. Determining the genetic cause of CP helps to predict the disease course, perform preventive measures in the proband's relatives, and facilitate a personalized treatment of the patient in future.

Multilocus Analysis of Phylogenetic Relationships in the Crocidura suaveolens Sensu Lato Species Complex: A Comparison with Mitochondrial Data.

Multilocus analysis was for the first time used to study the phylogeny of the Crocidura suaveolens s. l. species complex. Sequencing data for 16 nuclear genes indicated that several distinct forms exist within the species complex. The structure of the complex did generally not contradict its mitochondrial phylogeny. Siberian shrew showed certain specificity of the nuclear genome, but the degree of its genetic differentiation did not correspond to the species level. Relationships of Crocidura aff. suaveolens from South Gansu and Sichuan with other forms of the species complex were clarified. Shrews from Buryatia and Khentei also belong to this form, but their mtDNA apparently introgressed from C. shantungensis in the past. Hybridization of C. suaveolens s. str. with C. aff. suaveolens and C. güeldenstaedtii occurred recently. Due to multiple introgression events in the history of C. suaveolens s. l., a far larger set of loci is necessary for the analysis of the phylogenetic relationships between its forms.

[Monitoring the spread of the SARS-CoV-2 (Coronaviridae: Coronavirinae: Betacoronavirus; Sarbecovirus) variants in the Moscow region using targeted high-throughput sequencing].

INTRODUCTION: Since the outbreak of the COVID-19 pandemic caused by SARS-CoV-2 novel coronavirus, the international community has been concerned about the emergence of mutations altering some biological properties of the pathogen like increasing its infectivity or virulence. Particularly, since the end of 2020, several variants of concern have been identified around the world, including Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), and Delta (B.1.617.2). However, the existing mechanism of detecting important mutations are not always effective enough, since only a relatively small part of all pathogen samples can be examined by whole genome sequencing due to its high cost. MATERIAL AND METHODS: In this study, we have designed special primer panel and used it for targeted highthroughput sequencing of several significant S-gene (spike) regions of SARS-CoV-2. The Illumina platform averaged approximately 50,000 paired-end reads with a length of ≥150 bp per sample. This method was used to examine 579 random samples obtained from COVID-19 patients in Moscow and the Moscow region from February to June 2021. RESULTS: This study demonstrated the dynamics of distribution of several SARS-CoV-2 strains and its some single mutations. It was found that the Delta strain appeared in the region in May 2021, and became prevalent in June, partially displacing other strains. DISCUSSION: The obtained results provide an opportunity to assign the viral samples to one of the strains, including the previously mentioned in time- and cost-effective manner. The approach can be used for standardization of the procedure of searching for mutations in individual regions of the SARS-CoV-2 genome. It allows to get a more detailed data about the epidemiological situation in a region.

Evaluating the accuracy and sensitivity of detecting minority HIV-1 populations by Illumina next-generation sequencing.

The accuracy and sensitivity of deep sequencing were assessed using viral standards (pNL4-3 and pLAI.2) of both DNA and RNA. The sequencing accuracy did not depend on the type of nucleic acid, but critically depended on the number of reads and threshold of sensitivity to minor viral populations. With coverage of more than 236 reads, the accuracy of viral RNA sequencing was equal to or exceeded 99.9%, with a sensitivity threshold to minor nucleotides of 20%. When the sensitivity threshold was below 1%, reduced accuracy dynamics were clearly visible even when the coverage was massive (more than 9.000 reads). It was found that the floating sensitivity threshold allowed the sequencing accuracy to be maintained at an acceptable level in cases of low coverage (less than 1.500-2.000) of reads. These results indicate the quality that can be expected with a specific number of reads and sensitivity threshold. Deep sequencing is a very powerful tool that can significantly improve the value of study results, but despite its superior performance, it should be used with caution regarding its sensitivity to minor populations below 1%.

[GENETIC CHARACTERISTICS OF INFLUENZA A/H3N2 AND B VIRUSES THAT HAD CIRCULATED IN RUSSIA IN 2013 - 2015].

AIM: Establish genetic characteristics, carry out phylogenetic analysis and determination of molecular markers of resistance to etiotropic preparations against influenza A/H3N2 and B viruses that had circulated in Russia in 2013 - 2015. MATERIALS AND METHODS: 80 biological samples containing influenza A/H3N2 virus RNA and 31 samples containing influenza B virus RNA were studied. Sequencing of PCR fragments was carried out inABI-3 100 PRIZMTM GeneticAnalyzer (AppliedBiosystems, USA) and using MiSeq (Illumina, USA). Data treatment and analysis was carried out using CLC v.3.6.5., DNASTAR and BioNumerics v.6.5. programs. RESULTS: In 2013 -2014 A/Texas/50/2012-like-clade 3C.3 influenza A/H3N2 viruses dominated, 10% belonged to subclade 3C.2a and 10% - to 3C.3b. Most of the viruses (8 1%) of 2014 - 2015 were of 3C.2a clade, the portion of viruses belonging to 3C.3b and 3C.3a was 9 and 10%. Yamagata-like viruses predominated among the studied influenza B viruses, only 1 virus of 2014 - 2015 belonged to Victoria lineage, 1 reassortant of Yamagata and Victoria lineages was detected. Rimantadine- resistance mutationS3 lN(M2 protein) was detected in all the influenza A/H3N2 viruses. Mutations determining resistance to oseltamivir (NA gene) were not detected in influenza A/H3N2 and B viruses. CONCLUSION: Increase of influenza morbidity in 2014 - 2015 was determined by the emergence of influenza A/H3N2 and B viruses, antigenically distinct from those that had circulated previously and those included into the vaccine, thus resulting in the WHO decision to change A/ H3N2 and B components of the 2015 - 2016 vaccine: Simultaneous circulation of 2 lineages of influenza B virus and emergence of their reassortants gives evidence on the necessity of use of quadrivalent vaccines, containing both lineages.

A Rapid and Cost-Effective Method for DNA Extraction from Archival Herbarium Specimens.

Here we report a rapid and cost-effective method for the extraction of total DNA from herbarium specimens up to 50-90-year-old. The method takes about 2 h, uses AMPure XP magnetic beads diluted by PEG-8000- containing buffer, and does not require use of traditional volatile components like chloroform, phenol, and liquid nitrogen. It yields up to 4 µg of total nucleic acid with high purity from about 30 mg of dry material. The quality of the extracted DNA was tested by PCR amplification of 5S rRNA and rbcL genes (nuclear and chloroplast DNA markers) and compared against the traditional chloroform/isoamyl alcohol method. Our results demonstrate that the use of the magnetic beads is crucial for extraction of DNA suitable for subsequent PCR from herbarium samples due to the decreasing inhibitor concentrations, reducing short fragments of degraded DNA, and increasing median DNA fragment sizes.

A new reliable reference gene UBA52 for quantitative real-time polymerase chain reaction studies in pyloric cecal tissues of the starfish Asterias rubens.

The starfish Asterias rubens is one of the most abundant echinoderm species in the White, Barents, North, and Baltic Seas. This species is an important component of marine ecosystems and a model object for certain biological studies, in particular those requiring quantitative estimation of gene expression. As a rule, expression at the transcriptional level is estimated by real-time qPCR using the ΔΔCt method, which allows the comparison of the copy number of target gene transcripts in samples with unknown mRNA/cDNA concentration. Application of this method requires normalization of the results relative to genes with stable expression levels (reference genes). The identification of reference genes is still a challenging task since data of this kind are missing for certain taxa, whereas the use of "standard" endogenous control genes without additional tests might lead to erroneous conclusions. We performed a preliminary analysis of the expression of many housekeeping genes in the pyloric ceca of A. rubens by high-throughput sequencing under normal and heat shock conditions. For one of them, the ubiquitin gene UBA52, low variation of expression (not greater than 2-fold) was shown using real-time qPCR. Tissues of pyloric ceca of normal adults and underyearlings and of adults after heat shock were used. The data obtained suggest that the UBA52 gene may be used as reference for normalization of gene expression at the mRNA level in the starfish A. rubens and probably in closely related species.

[Differential expression of genes that encode glycolysis enzymes in kidney and lung cancer in humans].

Glycolysis is a main catabolic pathway of glucose metabolism, accompanied by ATP synthesis. More than 30 enzymes are involved in glycolysis, and genes that encode them can be considered housekeeping genes due to the high conservatism and evolutionary antiquity of the process. We studied the expression of these genes in kidney papillary cancer and planocellular lung cancer via the bioinformatic analysis of transcriptome database and method of quantitative real time PCR. Quantitative analysis of mRNA level demonstrated that only a part ofgenes that encode glycolysis enzymes maintain relatively stable mRNA level, including the HK1, ADPGK, GPI, PGK1, and PKM2 genes in kidney papillary cancer and the ADPGK, ALDOA, GAPDH, PGK1, BPGM, ENO1, and PKM2 genes in planocellular lung cancer. The frequent increase in the mRNA expression of PFKP, ALDOA, and GAPDH genes in kidney cancer, as well as the GPI gene in lung cancer, were detected for the first time by real time PCR. For other genes, their differential expression was demonstrated; the cases of both a decrease and increase in the mRNA level were detected. Thus, several genes that can be used as control genes in transcriptome analysis by real time PCR in kidney and lung cancer, as well as a number of differentially expressed genes that can be potential oncomarkers, were identified.

Heterologous expression, purification, and properties of a potato protein inhibitor of serine proteinases.

The gene PKPI-B10 [AF536175] encoding in potato (Solanum tuberosum L., cv. Istrinskii) a Kunitz-type protein inhibitor of proteinases (PKPI) has been cloned into the pET23a vector and then expressed in Escherichia coli. The recombinant protein PKPI-B10 obtained as inclusion bodies was denatured, separated from admixtures by ion-exchange fast protein liquid chromatography (FPLC) on MonoQ under denaturing conditions, and renatured. The native protein was additionally purified by ion-exchange FPLC on DEAE-Toyopearl. The PKPI-B10 protein effectively inhibits the activity of trypsin, significantly weaker suppresses the activity of chymotrypsin, and has no effect on other serine proteinases: human leukocyte elastase, subtilisin Carlsberg, and proteinase K, and also the plant cysteine proteinase papain.

Effects of misoprostol and taurine on monochloramine ulcerogenesis in rats.

We compared gastroprotective characteristics of synthetic prostaglandin E1 misoprostol and amino acid taurine on rat model of monochloramine injury to the gastric mucosa. Both substances exhibited a pronounced gastroprotective effect.

Molecular cloning of Kunitz-type proteinase inhibitor group B genes from potato.

Eighteen clones representing copies of four Kunitz-type proteinase inhibitor group B genes (PKPI-B) obtained by polymerase chain reaction cloning of potato (Solanum tuberosum L. cv. Istrinskii) genomic DNA were sequenced and analyzed. Three new genes were found and named PKPI-B1, PKPI-B2, and PKPI-B10: these were represented by five, two, and seven clones, respectively. The remaining four clones corresponded to the formerly characterized PKPI-B9 gene. These data show that at least four PKPI-B encoding genes are harbored in the genome of potato cv. Istrinskii. Their analysis suggests that variability of PKPI-B encoding genes in potato is limited and could be explained by cross-hybridization events in the ancestor forms rather than by random mutagenesis.

Subtilisin protein inhibitor from potato tubers.

A protein with molecular weight of 21 kD denoted as PKSI has been isolated from potato tubers (Solanum tuberosum L., cv. Istrinskii). The isolation procedure includes precipitation with (NH4)2SO4, gel chromatography on Sephadex G-75, and ion-exchange chromatography on CM-Sepharose CL-6B. The protein effectively inhibits the activity of subtilisin Carlsberg (Ki = 1.67 +/- 0.2 nM) by stoichiometric complexing with the enzyme at the molar ratio of 1 : 1. The inhibitor has no effect on trypsin, chymotrypsin, and the cysteine proteinase papain. The N-terminal sequence of the protein consists of 19 amino acid residues and is highly homologous to sequences of the known inhibitors from group C of the subfamily of potato Kunitz-type proteinase inhibitors (PKPIs-C). By cloning PCR products from the genomic DNA of potato, a gene denoted as PKPI-C2 was isolated and sequenced. The N-terminal sequence (residues from 15 to 33) of the protein encoded by the PKPI-C2 gene is identical to the N-terminal sequence (residues from 1 to 19) of the isolated protein PKSI. Thus, the inhibitor PKSI is very likely encoded by this gene.