Characterization of Seminal Microbiome of Infertile Idiopathic Patients Using Third-Generation Sequencing Platform.

Since the first description of a commensal seminal microbiome using sequencing, less than a decade ago, interest in the composition of this microbiome and its relationship with fertility has been growing. Articles using next-generation sequencing techniques agree on the identification of the most abundant bacterial phyla. However, at the genus level, there is still no consensus on which bacteria are most abundant in human seminal plasma. This discrepancy may be due to methodological variability such as sample collection, bacterial DNA extraction methodology, which hypervariable regions of 16S rRNA gene have been amplified, or bioinformatic analysis. In the present work, seminal microbiota of 14 control samples and 42 samples of idiopathic infertile patients were characterized based on full-length sequencing of the 16S rRNA gene using MinION platform from Oxford Nanopore. These same samples had been analyzed previously using Illumina's MiSeq sequencing platform. Comparison between the results obtained with the two platforms has been used to analyze the impact of sequencing method on the study of the seminal microbiome's composition. Seminal microbiota observed with MinION were mainly composed of the phyla Firmicutes, Proteobacteria, Bacteroidetes and Actinobacteria, with the most abundant genera being Peptoniphilus, Finegoldia, Staphylococcus, Anaerococcus, Campylobacter, Prevotella, Streptococcus, Lactobacillus, Ezakiella and Enterococcus. This composition was similar to that found by the Illumina platform, since these 10 most abundant genera were also among the most abundant genera detected by the Nanopore platform. In both cases, the top 10 genera represented more than 70% of the classified reads. However, relative abundance of each bacterium did not correlate between these two platforms, with intraindividual variations of up to 50 percentage points in some cases. Results suggest that the effect of the sequencing platform on the characterization of seminal microbiota is not very large at the phylum level, with slightly variances in Firmicutes and Actinobacteria, but presents differences at the genus level. These differences could alter the composition and diversity of bacterial profiles or posterior analyses. This indicates the importance of conducting multi-platform studies to better characterize seminal microbioma.

Key Genes of the Immune System and Predisposition to Acquired Hemophilia A: Evidence from a Spanish Cohort of 49 Patients Using Next-Generation Sequencing.

Acquired hemophilia A (AHA) is a rare bleeding disorder caused by the presence of autoantibodies against factor VIII (FVIII). As with other autoimmune diseases, its etiology is complex and its genetic basis is unknown. The aim of this study was to identify the immunogenetic background that predisposes individuals to AHA. HLA and KIR gene clusters, as well as KLRK1, were sequenced using next-generation sequencing in 49 AHA patients. Associations between candidate genes involved in innate and adaptive immune responses and AHA were addressed by comparing the alleles, genotypes, haplotypes, and gene frequencies in the AHA cohort with those in the donors' samples or Spanish population cohort. Two genes of the HLA cluster, as well as rs1049174 in KLRK1, which tags the natural killer (NK) cytotoxic activity haplotype, were found to be linked to AHA. Specifically, A*03:01 (p = 0.024; odds ratio (OR) = 0.26[0.06-0.85]) and DRB1*13:03 (p = 6.8 × 103, OR = 7.56[1.64-51.40]), as well as rs1049174 (p = 0.012), were significantly associated with AHA. In addition, two AHA patients were found to carry one copy each of the low-frequency allele DQB1*03:09 (nallele = 2, 2.04%), which was completely absent in the donors. To the best of our knowledge, this is the first time that the involvement of these specific alleles in the predisposition to AHA has been proposed. Further molecular and functional studies will be needed to unravel their specific contributions. We believe our findings expand the current knowledge on the genetic factors involved in susceptibility to AHA, which will contribute to improving the diagnosis and prognosis of AHA patients.

Natural killer cell receptors and ligand variants modulate response to tyrosine kinase inhibitors in patients with chronic myeloid leukemia.

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm treated with tyrosine kinase inhibitors (TKIs). Although survival rates have improved, response to these treatments is highly heterogeneous. Variations in response rates may be due to different causes such as, treatment adherence, mutations in the BCR-ABL1 gene, clonal evolution and amplification of the BCR-ABL1 gene, but innate immune response is also considered to play a very important role and, specifically, NK cell activity through their receptors and ligands, could be determinant. The aim of this retrospective study was to explore the role of different activating and inhibiting KIR genes as well as the activating NKG2D receptor, present in NK cells, and also their respective ligands, HLA-A, -B, -C, -G, -F, MICA and MICB, in the progression of 190 patients with CML and treated at two hospitals from Barcelona between 2000 and 2019. Early molecular response (EMR), major molecular response (MMR) or MR3.0 and deep molecular response (DMR) or MR4.0 were correlated. As control samples, healthy donors from the Barcelona Blood Bank were analyzed. The presence of KIR2DL2/KIR2DS2 was associated with the achievement of EMR, MR3.0, and MR4.0. Carriers of the higher expression NKG2D variant and MICA*009:01 were also likely to achieve molecular response (MR). The most remarkable difference between CML patients and controls was a higher frequency of the lower expression NKG2D variant in CML patients. In summary, our results showed that activating NK receptor phenotypes might help to achieve MR and DMR in CML patients treated with TKIs although confirmatory studies are necessary.

Seminal Microbiota of Idiopathic Infertile Patients and Its Relationship With Sperm DNA Integrity.

The development of new biomarkers for human male infertility is crucial to improve the diagnosis and the prognosis of this disease. Recently, seminal microbiota was shown to be related to sperm quality parameters, suggesting an effect in human fertility and postulating it as a biomarker candidate. However, its relationship to sperm DNA integrity has not been studied yet. The aim of the present study is to characterize the seminal microbiota of a western Mediterranean population and to evaluate its relationship to sperm chromatin integrity parameters, and oxidative stress. For that purpose, 14 samples from sperm donors and 42 samples from infertile idiopathic patients were obtained and were analyzed to assess the composition of the microbiota through full-length 16S rRNA gene sequencing (Illumina MiSeq platform). Microbial diversity and relative abundances were compared to classic sperm quality parameters (macroscopic semen parameters, motility, morphology and concentration), chromatin integrity (global DNA damage, double-stranded DNA breaks and DNA protamination status) and oxidative stress levels (oxidation-reduction potential). The seminal microbiota observed of these samples belonged to the phyla Firmicutes, Proteobacteria, Actinobacteria and Bacteroidetes. The most abundant genera were Finegoldia, Peptoniphilus, Anaerococcus, Campylobacter, Streptococcus, Staphylococcus, Moraxella, Prevotella, Ezakiella, Corynebacterium and Lactobacillus. To our knowledge, this is the first detection of Ezakiella genus in seminal samples. Two clusters of microbial profiles were built based on a clustering analysis, and specific genera were found with different frequencies in relation to seminal quality defects. The abundances of several bacteria negatively correlate with the sperm global DNA fragmentation, most notably Moraxella, Brevundimonas and Flavobacterium. The latter two were also associated with higher sperm motility and Brevundimonas additionally with lower oxidative-reduction potential. Actinomycetaceae, Ralstonia and Paenibacillus correlated with reduced chromatin protamination status and increased double-stranded DNA fragmentation. These effects on DNA integrity coincide in many cases with the metabolism or enzymatic activities of these genera. Significant differences between fertile and infertile men were found in the relative presence of the Propionibacteriaceae family and the Cutibacterium, Rhodopseudomonas and Oligotropha genera, which supports its possible involvement in male fertility. Our findings sustain the hypothesis that the seminal microbiome has an effect on male fertility.

High-throughput genotyping of HLA-G, HLA-F, MICA, and MICB and analysis of frequency distributions in healthy blood donors from Catalonia.

Similarly to HLA class I molecules, certain non-classical HLA class I genes and MHC class I polypeptide-related sequences A and B (MICA and MICB) act as ligands for KIR and NKG2D natural killer receptors. Although these genes are less polymorphic than HLA class I, few studies have analyzed their association with diseases. Information on allele frequencies in healthy donors is needed to map their distribution worldwide. This study is the first to analyze high-resolution HLA-G, HLA-F, MICA, and MICB allele frequencies using a novel high-throughput next generation-sequencing method. We analyzed DNA samples from 96 unrelated blood donors resident in Catalonia, Spain, and registered in the Barcelona Blood and Tissue Bank. Using the first two fields of the HLA nomenclature, we detected six HLA-G and two HLA-F alleles. The most frequent alleles were HLA-G*01:01 (77.08%) and HLA-F*01:01(84.90%). When the four fields were analyzed, we detected 16 and 10 alleles, respectively. Nineteen alleles were detected for MICA and 10 for MICB. The most frequent alleles in these cases were MICA*008:01 (16.15%) and MICB*005:02 (46.84%). All frequencies were in Hardy Weinberg equilibrium except MICA. We also estimated maximum-likelihood haplotype frequencies and calculated corresponding linkage disequilibrium (LD) values and found that few allele pairs were in disequilibrium. Strong LD between MICA and HLA-B (using data from a previous study) was observed. Our findings will be useful for guiding further research evaluating the functional role of these genes in different diseases and populations.

Distribution of human killer cell immunoglobulin-like receptors and ligands among blood donors of Catalonia.

Killer cell immunoglobulin-like receptors (KIR) belong to a diverse family of receptors present in NK cells and certain subsets of T cells. They bind to HLA class I molecules and, such as these, are highly polymorphic. This study is the first to determine KIR gene content diversity and perform high-resolution genotyping of corresponding HLA class I ligands in 445 unrelated blood donors resident in Catalonia (northeast Spanish region) and registered in the Blood and Tissue Bank of Barcelona between 2017 and 2018. The study was performed using a new next-generation sequencing method developed and validated by our group, with an important hands-on-time reduction and very cost-effectiveness. Thirty-eight KIR genotypes were observed. Haplotype combination frequencies were 24.49% for AA, 75.51% for BX. The most frequent genotype, found in 109 individuals (24.49%), was genotype ID 1. Seventeen genotypes occurred only once. Thirty individuals carried all KIR genes (genotype ID 6). We detected 46 HLA-A, 63 HLA-A-B, and 40 HLA-A-C alleles. A majority of individuals were heterozygous for C1/C2 (42.24%), 38.43% were homozygous for C1, and 19.33% were homozygous for C2. The most common KIR-HLA ligand pair was KIR3DL1-Bw4+, present in 75.73% of cases. No differences were found in KIR gene frequencies between the Catalan cohort and other Iberian Peninsula populations. Our findings will be useful for guiding further research evaluating the functional significance of KIR-ligand associations in specific diseases.

Design and Validation of a Multiplex KIR and HLA Class I Genotyping Method Using Next Generation Sequencing.

Killer cell immunoglobulin-like receptors (KIR), considered the most polymorphic natural killer (NK) cell regulators, bind HLA class-I molecules or still unknown ligands. Interest in KIR genotyping is increasing because of the importance of these receptors for identifying the best possible donor in hematopoietic stem cell transplantation to obtain a graft-versus-leukemia effect. Currently, routine protocols to determine the gene content of the KIR cluster are exclusively performed by PCR-SSO and PCR-SSP. To improve the study of these genes, we developed a multiplex, long-range PCR strategy suitable for simultaneous, high-resolution HLA class I and KIR genotyping by next generation sequencing (NGS). This protocol allows amplification of the 14 KIR genes, 2 KIR pseudogenes, and HLA class I genes, with subsequent sequencing on an Illumina platform. The bioinformatics analysis for KIR genotyping was performed by virtual hybridization of gene-specific probes, and HLA genotyping was done by GenDx NGSengine software. To validate the method reliability, 192 genomic DNA samples previously characterized by PCR-SSO were used. When a specific KIR gene was present, a large number of gene-specific virtual probes were detected, whereas when it was absent, very few or none were found, enabling cutoff establishment. Concordance for both the KIR and HLA assignments as compared with the previous characterization was 100%. In conclusion, the multiplex PCR NGS-based strategy presented could provide an efficient, less costly method for KIR-ligand genotyping by gene presence/absence. Furthermore, allele resolution will be possible when KIR-specific software becomes available.