Genetic Alternatives for Experimental Adaptation to Colistin in Three Pseudomonas aeruginosa Lineages.

Pseudomonas aeruginosa is characterized by a high adaptive potential, developing resistance in response to antimicrobial pressure. We employed a spatiotemporal evolution model to disclose the pathways of adaptation to colistin, a last-resort polymyxin antimicrobial, among three unrelated P. aeruginosa lineages. The P. aeruginosa ATCC-27833 reference strain (Pa_ATCC), an environmental P. aeruginosa isolate (Pa_Environment), and a clinical isolate with multiple drug resistance (Pa_MDR) were grown over an increasing 5-step colistin concentration gradient from 0 to 400 mg/L. Pa_Environment demonstrated the highest growth pace, achieving the 400 mg/L band in 15 days, whereas it took 37 and 60 days for Pa_MDR and Pa_ATCC, respectively. To identify the genome changes that occurred during adaptation to colistin, the isolates selected during the growth of the bacteria (n = 185) were subjected to whole genome sequencing. In total, 17 mutation variants in eight lipopolysaccharide-synthesis-associated genes were detected. phoQ and lpxL/PA0011 were affected in all three lineages, whereas changes in pmrB were found in Pa_Environment and Pa_MDR but not in Pa_ATCC. In addition, mutations were detected in 34 general metabolism genes, and each lineage developed mutations in a unique set of such genes. Thus, the three examined distinct P. aeruginosa strains demonstrated different capabilities and genetic pathways of colistin adaptation.

Microbiome markers in HPV-positive and HPV-negative women of reproductive age with ASCUS and SIL determined by V4 region of 16S rRNA gene sequencing.

Introduction: Human papilloma virus (HPV) is the most common sexually transmitted infection worldwide. Cervicovaginal microbiota plays an important role in HPV infection and is associated with the development of squamous intraepithelial lesions (SIL). The natural history of cervical cancer involves reversible changes in the cervical tissue from a normal state, in which no neoplastic changes are detected in the squamous epithelium, to varying states of cellular abnormalities that ultimately lead to cervical cancer. Low-grade SIL (LSIL), like another cytological category - atypical squamous cells of undetermined significance (ASCUS), may progress to high-grade SIL (HSIL) and invasive cervical cancer or may regress to a normal state. Methods: In this work, we studied cervical canal microbiome in 165 HPV-positive and HPV-negative women of a reproductive age with ASCUS [HPV(+) n = 29; HPV(-) n = 11], LSIL [HPV(+) n = 32; HPV(-) n = 25], HSIL [HPV(+) n = 46], and the control group with negative for intraepithelial lesion malignancy (NILM) [HPV(-) n = 22]. Results and Discussion: HPV16 is the most prevalent HPV type. We have not found any differences between diversity in studied groups, but several genus [like Prevotella (p-value = 0.026), Gardnerella (p-value = 0.003), Fannyhessea (p-value = 0.024)] more often occurred in HSIL group compared by NILM or LSIL regardless of HPV. We have found statistically significant difference in occurrence or proportion of bacterial genus in studied groups. We also identified that increasing of the ratio of Lactobacillus iners or age of patient lead to higher chance to HSIL, while increasing of the ratio of Lactobacillus crispatus lead to higher chance to LSIL. Patients with a moderate dysbiosis equally often had either of three types of vaginal microbial communities (CST, Community State Type) with the prevalence of Lactobacillus crispatus (CST I), Lactobacillus gasseri (CST II), and Lactobacillus iners (CST III); whereas severe dysbiosis is linked with CST IV involving the microorganisms genera associated with bacterial vaginosis and aerobic vaginitis: Gardnerella, Fannyhessea, Dialister, Sneathia, Anaerococcus, Megasphaera, Prevotella, Finegoldia, Peptoniphilus, Porphyromonas, Parvimonas, and Streptococcus.

Clonal diversity, antimicrobial resistance, and genome features among nonfermenting gram-negative bacteria isolated from patients with cystic fibrosis in Russia.

Nonfermenting gram-negative (NFGN) bacteria were isolated from cystic fibrosis (CF) patients and subjected to susceptibility testing and whole-genome sequencing. Among 170 enrolled CF patients, 112 (65.9%) were colonized with at least 1 key NFGN species. The species-specific infection rate was highest for Pseudomonas aeruginosa (40.6%) followed by Stenotrophomonas maltophilia (14.1%), Achromobacter spp. (9.4%), and Burkholderia cepacia complex (Bcc, 8.2%) demonstrating a significant age-dependent increase for P. aeruginosa and Achromobacter spp., but not for S. maltophilia or Bcc. P. aeruginosa sequence types (STs) related to high-risk epidemic and global CF clones were carried by 12 (7.1%) and 13 (7.6%) patients, respectively. In total, 47% NFGN isolates, predominantly P. aeruginosa, harbored at least 1 plasmid-borne resistance gene; 5 ST235 isolates carried blaVIM2. Pathogenicity island-borne virulence genes were harbored by 9% NFGN isolates. These findings in conjunction with frequent early colonization by Bcc raised serious concerns regarding infection control in Russian CF centers.

Benchmarking DNA isolation methods for marine metagenomics.

Metagenomics is a powerful tool to study marine microbial communities. However, obtaining high-quality environmental DNA suitable for downstream sequencing applications is a challenging task. The quality and quantity of isolated DNA heavily depend on the choice of purification procedure and the type of sample. Selection of an appropriate DNA isolation method for a new type of material often entails a lengthy trial and error process. Further, each DNA purification approach introduces biases and thus affects the composition of the studied community. To account for these problems and biases, we systematically investigated efficiency of DNA purification from three types of samples (water, sea sediment, and digestive tract of a model invertebrate Magallana gigas) with eight commercially available DNA isolation kits. For each kit-sample combination we measured the quantity of purified DNA, extent of DNA fragmentation, the presence of PCR-inhibiting contaminants, admixture of eukaryotic DNA, alpha-diversity, and reproducibility of the resulting community composition based on 16S rRNA amplicons sequencing. Additionally, we determined a "kitome", e.g., a set of contaminating taxa inherent for each type of purification kit used. The resulting matrix of evaluated parameters allows one to select the best DNA purification procedure for a given type of sample.

Candidate Genes for IgA Nephropathy in Pediatric Patients: Exome-Wide Association Study.

IgA nephropathy (IgAN) is an autoimmune disorder which is believed to be non-monogenic. We performed an exome-wide association study of 70 children with IgAN and 637 healthy donors. The HLA allele frequencies were compared between the patients and healthy donors from the bone marrow registry of the Pirogov University. We tested 78,020 gene markers for association and performed functional enrichment analysis and transcription factor binding preference detection. We identified 333 genetic variants, employing three inheritance models. The most significant association with the disorder was observed for rs143409664 (PRAG1) in the case of the additive and dominant models (PBONF = 1.808 × 10-15 and PBONF = 1.654 × 10-15, respectively), and for rs13028230 (UBR3) in the case of the recessive model (PBONF = 1.545 × 10-9). Enrichment analysis indicated the strongly overrepresented "immune system" and "kidney development" terms. The HLA-DQA1*01:01:01G allele (p = 0.0076; OR, 2.021 [95% CI, 1.322-3.048]) was significantly the most frequent among IgAN patients. Here, we characterized, for the first time, the genetic background of Russian IgAN patients, identifying the risk alleles typical of the population. The most important signals were detected in previously undescribed loci.

The Search of Association of HLA Class I and Class II Alleles with COVID-19 Mortality in the Russian Cohort.

HLA genes play a pivotal role in the immune response via presenting the pathogen peptides on the cell surface in a host organism. Here, we studied the association of HLA allele variants of class I (loci A, B, C) and class II (loci DRB1, DQB1, DPB1) genes with the outcome of COVID-19 infection. We performed high-resolution sequencing of class HLA I and class II genes based on the sample population of 157 patients who died from COVID-19 and 76 patients who survived despite severe symptoms. The results were further compared with HLA genotype frequencies in the control population represented by 475 people from the Russian population. Although the obtained data revealed no significant differences between the samples at a locus level, they allowed one to uncover a set of notable alleles potentially contributing to the COVID-19 outcome. Our results did not only confirm the previously discovered fatal role of age or association of DRB1*01:01:01G and DRB1*01:02:01G alleles with severe symptoms and survival, but also allowed us to single out the DQB1*05:03:01G allele and B*14:02:01G~C*08:02:01G haplotype, which were associated with survival. Our findings showed that not only separate allele, but also their haplotype, could serve as potential markers of COVID-19 outcome and be used during triage for hospital admission.

SARS-CoV-2 evolution in a patient with secondary B-cell immunodeficiency: A clinical case.

The article highlights the course of long-term SARS-CoV-2 infection in a patient with a secondary immunodeficiency developed with B-cell-depleting therapy of the underlying disease. Analysis of the intrapatient virus evolution revealed an inpatient S:G75A mutation that alters the 72GTNGTKR78 motif of the S-protein, with a possible role in binding to alternative cellular receptors. Therapy with a ready-made COVID-19-globulin preparation (native human immunoglobulin G (IgG) derived from the plasma of convalescent COVID-19-patients) resulted in rapid improvement of the patient's condition, fast, and stable elimination of the virus, and passive immunization of the patient for at least 30 days. The results suggest the use of products containing neutralizing antibodies opens new prospects for treatment algorithms for patients with persistent coronavirus infection, as well as for passive immunization schemes for patients with a presumably reduced specific response to vaccination.

Performance comparison of Agilent new SureSelect All Exon v8 probes with v7 probes for exome sequencing.

Exome sequencing is becoming a routine in health care, because it increases the chance of pinpointing the genetic cause of an individual patient's condition and thus making an accurate diagnosis. It is important for facilities providing genetic services to keep track of changes in the technology of exome capture in order to maximize throughput while reducing cost per sample. In this study, we focused on comparing the newly released exome probe set Agilent SureSelect Human All Exon v8 and the previous probe set v7. In preparation for higher throughput of exome sequencing using the DNBSEQ-G400, we evaluated target design, coverage statistics, and variants across these two different exome capture products. Although the target size of the v8 design has not changed much compared to the v7 design (35.24 Mb vs 35.8 Mb), the v8 probe design allows you to call more of SNVs (+ 3.06%) and indels (+ 8.49%) with the same number of raw reads per sample on the common target regions (34.84 Mb). Our results suggest that the new Agilent v8 probe set for exome sequencing yields better data quality than the current Agilent v7 set.

Diplocloster agilis gen. nov., sp. nov. and Diplocloster modestus sp. nov., two novel anaerobic fermentative members of Lachnospiraceae isolated from human faeces.

Three novel strains of Gram-stain-negative, obligately anaerobic, spore-forming straight or slightly curved rods with pointed ends occurring singly or in pairs were isolated from the faeces of healthy human children. The strains were characterized by mesophilic fermentative metabolism and production of acetate, ethanol and H2 as the end metabolic products. Strains ASD3451 and ASD5720T were motile, fermented lactose and raffinose, and weakly fermented maltose. Strain ASD4241T was non-motile and did not ferment the carbohydrates listed above but fermented starch. Strains ASD3451 and ASD5720T shared average nucleotide identity higher than 98.5 % with each other, while ASD4241T had only 88.5-89 % identity to them. Based on phylogenetic and chemotaxonomic analyses, we propose Diplocloster agilis gen. nov., sp. nov. (ASD5720T=JCM 34353T=VKM B-3497T) and Diplocloster modestus sp. nov. (ASD4241T=JCM 34351T=VKM B-3498T) within the family Lachnospiraceae.

System analysis of the sequencing quality of human whole exome samples on BGI NGS platform.

Human exome sequencing is a classical method used in most medical genetic applications. The leaders in the field are the manufacturers of enrichment kits based on hybridization of cRNA or cDNA biotinylated probes specific for a genomic region of interest. Recently, the platforms manufactured by the Chinese company MGI Tech have become widespread in Europe and Asia. The reliability and quality of the obtained data are already beyond any doubt. However, only a few kits compatible with these sequencers can be used for such specific tasks as exome sequencing. We developed our own solution for library pre-capture pooling and exome enrichment with Agilent probes. In this work, using a set of the standard benchmark samples from the Platinum Genome collection, we demonstrate that the qualitative and quantitative parameters of our protocol which we called "RSMU_exome" exceed those of the MGI Tech kit. Our protocol allows for identifying more SNV and indels, generates fewer PCR duplicates, enables pooling of more samples in a single enrichment procedure, and requires less raw data to obtain results comparable with the MGI Tech's protocol. The cost of our protocol is also lower than that of MGI Tech's solution.

BC-store: A program for MGISEQ barcode sets analysis.

Here we present the devised BC-store-a program for analyzing and selecting sets of barcodes for sequencing on platforms manufactured by MGI Tech (China). The app is available as an open source in Python3 and as a desktop version. The application allows analyzing the compatibility of barcodes on a single lane of a flow cell in a set in the case of equal and arbitrary fractions. In addition, with the help of this tool barcodes can be added to an existing set with custom share options. In this paper we describe how BC-store works for different tasks and consider the effectiveness of using BC-store in sequence lab routine tasks.

The novel HLA-DQB1*05:02:24 allele, identified in a Russian bone marrow donor.

The new allele HLA-DQB1*05:02:24 showed one synonymous nucleotide difference with HLA-DQB1*05:02:01:01 in codon 140.

HLA-A*11:382N, a novel HLA-A null allele identified by next-generation sequencing.

The new allele HLA-A*11:382N showed one nucleotide difference with HLA-A*11:01:01:01 at codon 254 (nonsense mutation).

NIPT Technique Based on the Use of Long Chimeric DNA Reads.

Non-invasive prenatal testing (NIPT) for aneuploidy on Chromosomes 21 (T21), 18 (T18) and 13 (T13) is actively used in clinical practice around the world. One of the limitations of the wider implementation of this test is the high cost of the analysis itself, as high-throughput sequencing is still relatively expensive. At the same time, there is an increasing trend in the length of reads yielded by sequencers. Since extracellular DNA is short, in the order of 140-160 bp, it is not possible to effectively use long reads. The authors used high-performance sequencing of cell-free DNA (cfDNA) libraries that went through additional stages of enzymatic fragmentation and random ligation of the resulting products to create long chimeric reads. The authors used a controlled set of samples to analyze a set of cfDNA samples from pregnant women with a high risk of fetus aneuploidy according to the results of the first trimester screening and confirmed by invasive karyotyping of the fetus using laboratory and analytical approaches developed by the authors. They evaluated the sensitivity, specificity, PPV (positive predictive value), and NPV (negative predictive value) of the results. The authors developed a technique for constructing long chimeric reads from short cfDNA fragments and validated the test using a control set of extracellular DNA samples obtained from pregnant women. The obtained sensitivity and specificity parameters of the NIPT developed by the authors corresponded to the approaches proposed earlier (99.93% and 99.14% for T21; 100% and 98.34% for T18; 100% and 99.17% for T13, respectively).

Comparative analysis of novel MGISEQ-2000 sequencing platform vs Illumina HiSeq 2500 for whole-genome sequencing.

The MGISEQ-2000 developed by MGI Tech Co. Ltd. (a subsidiary of the BGI Group) is a new competitor of such next-generation sequencing platforms as NovaSeq and HiSeq (Illumina). Its sequencing principle is based on the DNB and the cPAS technologies, which were also used in the previous version of the BGISEQ-500 device. However, the reagents for MGISEQ-2000 have been refined and the platform utilizes updated software. The cPAS method is an advanced technology based on the cPAL previously created by Complete Genomics. In this paper, the authors compare the results of the whole-genome sequencing of a DNA sample from a Russian female donor performed on MGISEQ-2000 and Illumina HiSeq 2500 (both PE150). Two platforms were compared in terms of sequencing quality, number of errors and performance. Additionally, we performed variant calling using four different software packages: Samtools mpileaup, Strelka2, Sentieon, and GATK. The accuracy of SNP detection was similar in the data generated by MGISEQ-2000 and HiSeq 2500, which was used as a reference. At the same time, a separate indel analysis of the overall error rate revealed similar FPR values and lower sensitivity. It may be concluded with confidence that the data generated by the analyzed sequencing systems is characterized by comparable magnitudes of error and that MGISEQ-2000 and HiSeq 2500 can be used interchangeably for similar tasks like whole genome sequencing.

Compound heterozygous POMGNT1 mutations leading to muscular dystrophy-dystroglycanopathy type A3: a case report.

BACKGROUND: Dystroglycanopathies, which are caused by reduced glycosylation of alpha-dystroglycan, are a heterogeneous group of neurodegenerative disorders characterized by variable brain and skeletal muscle involvement. Muscle-eye-brain disease (or muscular dystrophy-dystroglycanopathy type 3 A) is an autosomal recessive disorder characterized by congenital muscular dystrophy, ocular abnormalities, and lissencephaly. CASE PRESENTATION: We report clinical and genetic characteristics of a 6-year-old boy affected by muscular dystrophy-dystroglycanopathy. He has severe a delay in psychomotor and speech development, muscle hypotony, congenital myopia, partial atrophy of the optic nerve disc, increased level of creatine kinase, primary-muscle lesion, polymicrogyria, ventriculomegaly, hypoplasia of the corpus callosum, cysts of the cerebellum. Exome sequencing revealed compound heterozygous mutations in POMGNT1 gene (transcript NM_001243766.1): c.1539 + 1G > A and c.385C > T. CONCLUSIONS: The present case report shows diagnostic algorithm step by step and helps better understand the clinical and genetic features of congenital muscular dystrophy.

Contribution of functional KIR3DL1 to ankylosing spondylitis.

Increasing evidence points to a role for killer immunoglobulin-like receptors (KIRs) in the development of autoimmune diseases. In particular, a positive association of KIR3DS1 (activating receptor) and a negative association of KIR3DL1 (inhibitory receptor) alleles with ankylosing spondylitis (AS) have been reported by several groups. However, none of the studies analyzed these associations in the context of functionality of polymorphic KIR3DL1. To better understand how the KIR3DL1/3DS1 genes determine susceptibility to AS, we analyzed the frequencies of alleles and genotypes encoding functional (KIR3DL1*F) and non-functional (KIR3DL1*004) receptors. We genotyped 83 AS patients and 107 human leukocyte antigen (HLA)-B27-positive healthy controls from the Russian Caucasian population using a two-stage sequence-specific primer PCR, which distinguishes KIR3DS1, KIR3DL1*F and KIR3DL1*004 alleles. For the patients carrying two functional KIR3DL1 alleles, those alleles were additionally genotyped to identify KIR3DL1*005 and KIR3DL1*007 alleles, which are functional but are expressed at low levels. KIR3DL1 was negatively associated with AS at the expense of KIR3DL1*F but not of KIR3DL1*004. This finding indicates that the inhibitory KIR3DL1 receptor protects against the development of AS and is not simply a passive counterpart of the segregating KIR3DS1 allele encoding the activating receptor. However, analysis of genotype frequencies indicates that the presence of KIR3DS1 is a more important factor for AS susceptibility than the absence of KIR3DL1*F. The activation of either natural killer (NK) or T cells via the KIR3DS1 receptor can be one of the critical events in AS development, while the presence of the functional KIR3DL1 receptor has a protective effect. Nevertheless, even individuals with a genotype that carried two inhibitory KIR3DL1 alleles expressed at high levels could develop AS.

Torque Teno Virus (TTV) distribution in healthy Russian population.

BACKGROUND: Torque teno virus (TTV) is a circular, single-stranded DNA virus that chronically infects healthy individuals of all ages worldwide. There is a lot of data on the prevalence and genetic heterogeneity of TTV in healthy populations and in patients with various diseases now available. However, little is known about TTV load among healthy human population. In this study we analyzed TTV load in the group of 512 Russian elite athletes, who are supposed to be, by some standards, the healthiest part of the human population. RESULTS: The prevalence rate of TTV among the Russian Olympic Reserve members was 94% (for test sensitivity about 1000 genome equivalents per 1 ml of blood). Quantities varied from 103 (which corresponded to detection limit) to 1010 copies per 1 ml of blood, with median at 2.7 x 106 copies. CONCLUSION: About 94% of healthy individuals in Russian population have more than 1000 TTV genome copies per 1 ml of blood. This result exceeds the previously published data, and can be explained by either more sensitive PCR test system or by higher TTV distribution in Russian population or both. TTV viral load neither depends on gender, nor age.