Genome-wide DNA methylome and transcriptome profiling reveals key genes involved in the dysregulation of adipose-stem cells in Crohn's disease.

BACKGROUND AND AIMS: Crohn's disease (CD) is characterised by the expansion of mesenteric adipose tissue (MAT), named creeping fat (CF), which seems to be directly related to disease activity. Adipose-stem cells (ASCs) isolated from the CF of patients with CD are extremely pro-inflammatory, which persists during disease remission. We hypothesised that the dysfunctional ASCs in CD accumulate epigenetic modifications triggered by the inflammatory environment that could serve as molecular markers. METHODS: Genome-wide DNA methylome and transcriptome profiling were performed in ASCs isolated from MAT adipose-tissue biopsies of patients with active and inactive disease and from non-Crohn's disease patients (non-CD). A validation cohort was used to test the main candidate genes via qPCR in other fat depots and immune cells. RESULTS: We found differences in DNA-methylation and gene expression between ASCs isolated from patients with CD and from non-CD subjects, but we found no differences related to disease activity. Pathway enrichment analysis revealed that oxidative stress and immune response were significantly enriched in active CD and integration analysis identified MAB21L2, a cell fate-determining gene, as the most affected gene in CD. Validation analysis confirmed the elevated gene expression of MAB21L2 in MAT and in adipose tissue macrophages in active CD. We also found a strong association between expression of the calcium channel subunit gene CACNA1H and disease remission, as CACNA1H expression was higher in ASCs and MAT from patients with inactive CD, and correlates negatively with C-reactive protein in peripheral blood mononuclear cells. CONCLUSION: We identified a potential gene signature of CD in ASCs obtained from MAT. Integration analysis highlighted two novel genes demonstrating a negative correlation between promoter DNA methylation and transcription: one linked to ASCs in CD (MAB21L2) and the other (CACNA1H) related to disease remission.

Altered Plasma microRNA Signature in Hospitalized COVID-19 Patients Requiring Oxygen Support.

To discover potential micro(mi)RNA biomarkers of SARS-CoV-2 infection and disease progression, large-scale deep-sequencing analysis of small RNA expression was performed on plasma samples from 40 patients hospitalized for SARS-CoV-2 infection (median 13.50 [IQR 9-24] days since symptoms initiation) and 21 healthy noninfected individuals. A total of 1218 different miRNAs were identified. When compared with healthy noninfected donors, SARS-CoV-2-infected patients showed significantly (fold change [FC] > 1.2 and adjusted p [padj] < 0.05) altered expression of 190 miRNAs. The top-10 differentially expressed (DE) miRNAs were miR-122-5p, let-7b-5p, miR-146a-5p, miR-342-3p, miR-146b-5p, miR-629-5p, miR-24-3p, miR-12136, let-7a-5p, and miR-191-5p, which displayed FC and padj values ranging from 153 to 5 and 2.51 × 10-32 to 2.21 × 10-21, respectively, which unequivocally diagnosed SARS-CoV-2 infection. No differences in blood cell counts and biochemical plasma parameters, including interleukin 6, ferritin, and D-dimer, were observed between COVID-19 patients on high-flow oxygen therapy, low-flow oxygen therapy, or not requiring oxygen therapy. Notably, 31 significantly deregulated miRNAs were found, when patients on high- and low-flow oxygen therapy were compared. SARS-CoV-2 infection generates a specific miRNA signature in hospitalized patients. Specific miRNA profiles are associated with COVID-19 prognosis in patients requiring oxygen flow.

Seminal extracellular vesicle sncRNA sequencing reveals altered miRNA/isomiR profiles as sperm retrieval biomarkers for azoospermia.

BACKGROUND: Non-invasive molecular biomarkers for classifying azoospermia by origin into either obstructive or non-obstructive/secretory azoospermia, as well as for inferring the spermatogenic reserve of the testis of non-obstructive/secretory azoospermia patients, are of great interest for testicular sperm retrieval outcome prediction for assisted reproduction. Prior analyses of semen small non-coding RNA expression in azoospermia have focused on microRNAs, but there has been a lack of attention on other regulatory small RNA species. In this regard, studying more in-depth expression changes of small non-coding RNA subtypes in small extracellular vesicles from semen of azoospermic individuals could be useful to select additional non-invasive biomarkers with diagnostic/prognostic purposes. MATERIAL AND METHODS: A high-throughput small RNA profiling analysis to determine the expression pattern of seminal small extracellular vesicle microRNAs (analyzed at the isomiR level), PIWI-interacting RNAs, and transfer RNA-derived small RNAs in normozoospermic (n = 4) and azoospermic (obstructive azoospermia because of pathological occurring obstruction in the genital tract, n = 4; secretory azoospermic individuals with positive testicular sperm extraction value, n = 5; secretory azoospermic individuals with negative testicular sperm extraction value, n = 4) individuals was carried out. Reverse transcriptase-quantitative real-time polymerase chain reaction validation analysis of selected microRNAs was additionally performed in a larger number of individuals. RESULTS AND DISCUSSION: Clinically relevant quantitative changes in the small non-coding RNA levels contained in semen small extracellular vesicles can be used as biomarkers for the origin of azoospermia and for predicting the presence of residual spermatogenesis. In this regard, canonical isoform microRNAs (n = 185) but also other isomiR variants (n = 238) stand out in terms of numbers and fold-change differences in expression, underlining the need to consider isomiRs when investigating microRNA-based regulation. Conversely, although transfer RNA-derived small RNAs are shown in our study to represent a high proportion of small non-coding RNA sequences in seminal small extracellular vesicle samples, they are not able to discriminate the origin of azoospermia. PIWI-interacting RNA cluster profiles and individual PIWI-interacting RNAs with significant differential expression were also not able to discriminate. Our study demonstrated that expression values of individual and/or combined canonical isoform microRNAs (miR-10a-5p, miR-146a-5p, miR-31-5p, miR-181b-5p; area under the receiver operating characteristic curve >0.8) in small extracellular vesicles provide considerable clinical value in identifying samples with a high likelihood of sperm retrieval while discriminating azoospermia by origin. Although no individual microRNA showed sufficient discriminating power on its own to identify severe spermatogenic disorders with focal spermatogenesis, multivariate microRNA models in semen small extracellular vesicles have the potential to identify those individuals with residual spermatogenesis. Availability and adoption of such non-invasive molecular biomarkers would represent a great improvement in reproductive treatment decision protocols for azoospermia in clinical practice.

3'IsomiR Species Composition Affects Reliable Quantification of miRNA/isomiR Variants by Poly(A) RT-qPCR: Impact on Small RNA-Seq Profiling Validation.

Small RNA-sequencing (small RNA-seq) has revealed the presence of small RNA-naturally occurring variants such as microRNA (miRNA) isoforms or isomiRs. Due to their small size and the sequence similarity among miRNA isoforms, their validation by RT-qPCR is challenging. We previously identified two miR-31-5p isomiRs-the canonical and a 3'isomiR variant (3' G addition)-which were differentially expressed between individuals with azoospermia of different origin. Here, we sought to determine the discriminatory capacity between these two closely-related miRNA isoforms of three alternative poly(A) based-RT-qPCR strategies in both synthetic and real biological context. We found that these poly(A) RT-qPCR strategies exhibit a significant cross-reactivity between these miR-31-5p isomiRs which differ by a single nucleotide, compromising the reliable quantification of individual miRNA isoforms. Fortunately, in the biological context, given that the two miRNA variants show changes in the same direction, RT-qPCR results were consistent with the findings of small RNA-seq study. We suggest that miRNA selection for RT-qPCR validation should be performed with care, prioritizing those canonical miRNAs that, in small RNA-seq, show parallel/homogeneous expression behavior with their most prevalent isomiRs, to avoid confounding RT-qPCR-based results. This is suggested as the current best strategy for robust biomarker selection to develop clinically useful tests.

Rapid and accurate quantification of isomiRs by RT-qPCR.

Currently, microRNAs (miRs) are annotated as a single defined sequence (canonical), even though high-throughput small RNA sequencing has identified miR isoforms (isomiRs) that differ from their canonical counterparts in length, sequence, or both. Here we describe a simple reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR)-based assay for quantification of the miR-100-5p_iso_3p:-2 variant. We chose miR-100-5p because the canonical sequence was underrepresented in our evaluation of human plasma. The quantification of miR-100-5p_iso_3 p:-2 from 57 plasma samples demonstrated high concordance between high-throughput RNA sequencing and RT-qPCR results (r = 0.55, p < 0.0001). Of note, we could not detect or quantify miR-100-5p in our plasma samples using a commercial TaqMan canonical miR-100-5p RT-qPCR kit. With these 57 samples, we also adapted this assay to specifically quantify the canonical sequences of miR-122-5p and miR-192-5p. Similar to the results obtained with miR-100-5p_iso_3p:-2, RT-qPCR results for miR-122-5p and miR-192-5p highly correlated with high-throughput RNA sequencing data (miR-122-5p: r = 0.44, p = 0.0005; miR-192-5p: r = 0.72, p < 0.0001). The assay described here can be easily adapted to many different identified isomiRs. Because of the high specificity of isomiRs, their reliable RT-qPCR-based quantification could provide greater resolution and higher accuracy than using canonical sequences.

GCAT|Panel, a comprehensive structural variant haplotype map of the Iberian population from high-coverage whole-genome sequencing.

The combined analysis of haplotype panels with phenotype clinical cohorts is a common approach to explore the genetic architecture of human diseases. However, genetic studies are mainly based on single nucleotide variants (SNVs) and small insertions and deletions (indels). Here, we contribute to fill this gap by generating a dense haplotype map focused on the identification, characterization, and phasing of structural variants (SVs). By integrating multiple variant identification methods and Logistic Regression Models (LRMs), we present a catalogue of 35 431 441 variants, including 89 178 SVs (≥50 bp), 30 325 064 SNVs and 5 017 199 indels, across 785 Illumina high coverage (30x) whole-genomes from the Iberian GCAT Cohort, containing a median of 3.52M SNVs, 606 336 indels and 6393 SVs per individual. The haplotype panel is able to impute up to 14 360 728 SNVs/indels and 23 179 SVs, showing a 2.7-fold increase for SVs compared with available genetic variation panels. The value of this panel for SVs analysis is shown through an imputed rare Alu element located in a new locus associated with Mononeuritis of lower limb, a rare neuromuscular disease. This study represents the first deep characterization of genetic variation within the Iberian population and the first operational haplotype panel to systematically include the SVs into genome-wide genetic studies.

Modulation of AggR levels reveals features of virulence regulation in enteroaggregative E. coli.

Enteroaggregative Escherichia coli (EAEC) strains are one of the diarrheagenic pathotypes. EAEC strains harbor a virulence plasmid (pAA2) that encodes, among other virulence determinants, the aggR gene. The expression of the AggR protein leads to the expression of several virulence determinants in both plasmids and chromosomes. In this work, we describe a novel mechanism that influences AggR expression. Because of the absence of a Rho-independent terminator in the 3'UTR, aggR transcripts extend far beyond the aggR ORF. These transcripts are prone to PNPase-mediated degradation. Structural alterations in the 3'UTR result in increased aggR transcript stability, leading to increased AggR levels. We therefore investigated the effect of increased AggR levels on EAEC virulence. Upon finding the previously described AggR-dependent virulence factors, we detected novel AggR-regulated genes that may play relevant roles in EAEC virulence. Mutants exhibiting high AggR levels because of structural alterations in the aggR 3'UTR show increased mobility and increased pAA2 conjugation frequency. Furthermore, among the genes exhibiting increased fold change values, we could identify those of metabolic pathways that promote increased degradation of arginine, fatty acids and gamma-aminobutyric acid (GABA), respectively. In this paper, we discuss how the AggR-dependent increase in specific metabolic pathways activity may contribute to EAEC virulence.

Paired-end small RNA sequencing reveals a possible overestimation in the isomiR sequence repertoire previously reported from conventional single read data analysis.

BACKGROUND: Next generation sequencing has allowed the discovery of miRNA isoforms, termed isomiRs. Some isomiRs are derived from imprecise processing of pre-miRNA precursors, leading to length variants. Additional variability is introduced by non-templated addition of bases at the ends or editing of internal bases, resulting in base differences relative to the template DNA sequence. We hypothesized that some component of the isomiR variation reported so far could be due to systematic technical noise and not real. RESULTS: We have developed the XICRA pipeline to analyze small RNA sequencing data at the isomiR level. We exploited its ability to use single or merged reads to compare isomiR results derived from paired-end (PE) reads with those from single reads (SR) to address whether detectable sequence differences relative to canonical miRNAs found in isomiRs are true biological variations or the result of errors in sequencing. We have detected non-negligible systematic differences between SR and PE data which primarily affect putative internally edited isomiRs, and at a much smaller frequency terminal length changing isomiRs. This is relevant for the identification of true isomiRs in small RNA sequencing datasets. CONCLUSIONS: We conclude that potential artifacts derived from sequencing errors and/or data processing could result in an overestimation of abundance and diversity of miRNA isoforms. Efforts in annotating the isomiRnome should take this into account.

Genotypic and Phenotypic Characterization of Staphylococcus aureus Isolates from the Respiratory Tract in Mechanically-Ventilated Patients.

Staphylococcus aureus is a commensal and frequent colonizer of the upper respiratory tract. When mechanical ventilation disrupts natural defenses, S. aureus is frequently isolated from the lower airways, but distinguishing between colonization and infection is difficult. The objectives of this study were (1) to investigate the bacterial genome sequence in consecutive isolates in order to identify changes related to the pathological adaptation to the lower respiratory tract and (2) to explore the relationship between specific phenotypic and genotypic features with the patient's study group, persistence of the clinical isolate and clinical outcome. A set of 94 clinical isolates were selected and corresponded to 34 patients that were classified as having pneumonia (10), tracheobronchitis (11) and bronchial colonization (13). Clinical strains were phenotypically characterized by conventional identification and susceptibility testing methods. Isolates underwent whole genome sequencing using Illumina HiSeq4000. Genotypic characterization was performed with an in-house pipeline (BacterialTyper). Genomic variation arising within-host was determined by comparing mapped sequences and de novo assemblies. Virulence factors important in staphylococcal colonization and infection were characterized using previously established functional assays. (1) Toxin production was assessed using a THP-1 cytotoxicity assay, which reports on the gross cytotoxicity of individual isolates. In addition, we investigated the expression of the major virulence factor, alpha-toxin (Hla) by Western blot. (2) Adhesion to the important extracellular matrix molecule, fibronectin, was determined using a standardized microtitre plate assay. Finally, invasion experiments using THP-1 and A539 cell lines and selected clinical strains were also performed. Repeated isolation of S. aureus from endotracheal aspirate usually reflects persistence of the same strain. Within-host variation is detectable in this setting, but it shows no evidence of pathological adaptation related to virulence, resistance or niche adaptations. Cytotoxicity was variable among isolates with 14 strains showing no cytotoxicity, with these latter presenting an unaltered Fn binding capacity. No changes on cytotoxicity were reported when comparing study groups. Fn binding capacity was reported for almost all strains, with the exception of two strains that presented the lowest values. Strains isolated from patients with pneumonia presented a lower capacity of adhesion in comparison to those isolated during tracheobronchitis (p = 0.002). Hla was detected in 71 strains (75.5%), with most of the producer strains in pneumonia and bronchial colonization group (p = 0.06). In our cohort, Hla expression (presence or absence) in sequential isolates was usually preserved (70%) although in seven cases the expression varied over time. No relationship was found between low cytotoxicity and intracellular persistence in invasion experiments. In our study population, persistent S. aureus isolation from airways in ventilated patients does not reflect pathological adaptation. There is an important diversity of sequence types. Cytotoxicity is variable among strains, but no association with study groups was found, whereas isolates from patients with pneumonia had lower adhesion capability. Favorable clinical outcome correlated with increased bacterial adhesion in vitro. Most of the strains isolated from the lower airways were Hla producers and no correlation with an adverse outcome was reported. The identification of microbial factors that contribute to virulence is relevant to optimize patient management during lower respiratory tract infections.

Gene Duplications in the Genomes of Staphylococci and Enterococci.

Gene duplications are a feature of bacterial genomes. In the present work we analyze the extent of gene duplications in the genomes of three microorganisms that belong to the Firmicutes phylum and that are etiologic agents of several nosocomial infections: Staphylococcus aureus, Enterococcus faecium, and Enterococcus faecalis. In all three groups, there is an irregular distribution of duplications in the genomes of the strains analyzed. Whereas in some of the strains duplications are scarce, hundreds of duplications are present in others. In all three species, mobile DNA accounts for a large percentage of the duplicated genes: phage DNA in S. aureus, and plasmid DNA in the enterococci. Duplicates also include core genes. In all three species, a reduced group of genes is duplicated in all strains analyzed. Duplication of the deoC and rpmG genes is a hallmark of S. aureus genomes. Duplication of the gene encoding the PTS IIB subunit is detected in all enterococci genomes. In E. faecalis it is remarkable that the genomes of some strains encode duplicates of the prgB and prgU genes. They belong to the prgABCU cluster, which responds to the presence of the peptide pheromone cCF10 by expressing the surface adhesins PrgA, PrgB, and PrgC.

The draft genome sequence of the spider Dysdera silvatica (Araneae, Dysderidae): A valuable resource for functional and evolutionary genomic studies in chelicerates.

BACKGROUND: We present the draft genome sequence of Dysdera silvatica, a nocturnal ground-dwelling spider from a genus that has undergone a remarkable adaptive radiation in the Canary Islands. RESULTS: The draft assembly was obtained using short (Illumina) and long (PaciBio and Nanopore) sequencing reads. Our de novo assembly (1.36 Gb), which represents 80% of the genome size estimated by flow cytometry (1.7 Gb), is constituted by a high fraction of interspersed repetitive elements (53.8%). The assembly completeness, using BUSCO and core eukaryotic genes, ranges from 90% to 96%. Functional annotations based on both ab initio and evidence-based information (including D. silvatica RNA sequencing) yielded a total of 48,619 protein-coding sequences, of which 36,398 (74.9%) have the molecular hallmark of known protein domains, or sequence similarity with Swiss-Prot sequences. The D. silvatica assembly is the first representative of the superfamily Dysderoidea, and just the second available genome of Synspermiata, one of the major evolutionary lineages of the "true spiders" (Araneomorphae). CONCLUSIONS: Dysderoids, which are known for their numerous instances of adaptation to underground environments, include some of the few examples of trophic specialization within spiders and are excellent models for the study of cryptic female choice. This resource will be therefore useful as a starting point to study fundamental evolutionary and functional questions, including the molecular bases of the adaptation to extreme environments and ecological shifts, as well of the origin and evolution of relevant spider traits, such as the venom and silk.

Gene duplications in the E. coli genome: common themes among pathotypes.

BACKGROUND: Gene duplication underlies a significant proportion of gene functional diversity and genome complexity in both eukaryotes and prokaryotes. Although several reports in the literature described the duplication of specific genes in E. coli, a detailed analysis of the extent of gene duplications in this microorganism is needed. RESULTS: The genomes of the E. coli enteroaggregative strain 042 and other pathogenic strains contain duplications of the gene that codes for the global regulator Hha. To determine whether the presence of additional copies of the hha gene correlates with the presence of other genes, we performed a comparative genomic analysis between E. coli strains with and without hha duplications. The results showed that strains harboring additional copies of the hha gene also encode the yeeR irmA (aec69) gene cluster, which, in turn, is also duplicated in strain 042 and several other strains. The identification of these duplications prompted us to obtain a global map of gene duplications, first in strain 042 and later in other E. coli genomes. Duplications in the genomes of the enteroaggregative strain 042, the uropathogenic strain CFT073 and the enterohemorrhagic strain O145:H28 have been identified by a BLASTp protein similarity search. This algorithm was also used to evaluate the distribution of the identified duplicates among the genomes of a set of 28 representative E. coli strains. Despite the high genomic diversity of E. coli strains, we identified several duplicates in the genomes of almost all studied pathogenic strains. Most duplicated genes have no known function. Transcriptomic analysis also showed that most of these duplications are regulated by the H-NS/Hha proteins. CONCLUSIONS: Several duplicated genes are widely distributed among pathogenic E. coli strains. In addition, some duplicated genes are present only in specific pathotypes, and others are strain specific. This gene duplication analysis shows novel relationships between E. coli pathotypes and suggests that newly identified genes that are duplicated in a high percentage of pathogenic E. coli isolates may play a role in virulence. Our study also shows a relationship between the duplication of genes encoding regulators and genes encoding their targets.

Streptococcus gallolyticus subsp. gallolyticus from human and animal origins: genetic diversity, antimicrobial susceptibility, and characterization of a vancomycin-resistant calf isolate carrying a vanA-Tn1546-like element.

The aim of this work was to characterize the antibiotic susceptibility and genetic diversity of 41 Streptococcus gallolyticus subsp. gallolyticus isolates: 18 isolates obtained from animals and 23 human clinical isolates. Antibiotic susceptibility was determined by the semiautomatic Wider system and genetic diversity by pulsed-field gel electrophoresis (PFGE) with SmaI. Animal isolates grouped separately in the PFGE analysis, but no statistical differences in antimicrobial resistance were found between the two groups. The LMG 17956 sequence type 28 (ST28) strain recovered from the feces of a calf exhibited high levels of resistance to vancomycin and teicoplanin (MIC, ≥256 mg/liter). Its glycopeptide resistance mechanism was characterized by Southern blot hybridization and a primer-walking strategy, and finally its genome, determined by whole-genome sequencing, was compared with four closely related S. gallolyticus subsp. gallolyticus genomes. Hybridization experiments demonstrated that a Tn1546-like element was integrated into the bacterial chromosome. In agreement with this finding, whole-genome sequencing confirmed a partial deletion of the vanY-vanZ region and partial duplication of the vanH gene. The comparative genomic analyses revealed that the LMG 17956 ST28 strain had acquired an unusually high number of transposable elements and had experienced extensive chromosomal rearrangements, as well as gene gain and loss events. In conclusion, S. gallolyticus subsp. gallolyticus isolates from animals seem to belong to lineages separate from those infecting humans. In addition, we report a glycopeptide-resistant isolate from a calf carrying a Tn1546-like element integrated into its chromosome.