Full-length 16S rRNA gene sequencing by PacBio improves taxonomic resolution in human microbiome samples.

BACKGROUND: Sequencing variable regions of the 16S rRNA gene (≃300 bp) with Illumina technology is commonly used to study the composition of human microbiota. Unfortunately, short reads are unable to differentiate between highly similar species. Considering that species from the same genus can be associated with health or disease it is important to identify them at the lowest possible taxonomic rank. Third-generation sequencing platforms such as PacBio SMRT, increase read lengths allowing to sequence the whole gene with the maximum taxonomic resolution. Despite its potential, full length 16S rRNA gene sequencing is not widely used yet. The aim of the current study was to compare the sequencing output and taxonomic annotation performance of the two approaches (Illumina short read sequencing and PacBio long read sequencing of 16S rRNA gene) in different human microbiome samples. DNA from saliva, oral biofilms (subgingival plaque) and faeces of 9 volunteers was isolated. Regions V3-V4 and V1-V9 were amplified and sequenced by Illumina Miseq and by PacBio Sequel II sequencers, respectively. RESULTS: With both platforms, a similar percentage of reads was assigned to the genus level (94.79% and 95.06% respectively) but with PacBio a higher proportion of reads were further assigned to the species level (55.23% vs 74.14%). Regarding overall bacterial composition, samples clustered by niche and not by sequencing platform. In addition, all genera with > 0.1% abundance were detected in both platforms for all types of samples. Although some genera such as Streptococcus tended to be observed at higher frequency in PacBio than in Illumina (20.14% vs 14.12% in saliva, 10.63% vs 6.59% in subgingival plaque biofilm samples) none of the differences were statistically significant when correcting for multiple testing. CONCLUSIONS: The results presented in the current manuscript suggest that samples sequenced using Illumina and PacBio are mostly comparable. Considering that PacBio reads were assigned at the species level with higher accuracy than Illumina, our data support the use of PacBio technology for future microbiome studies, although a higher cost is currently required to obtain an equivalent number of reads per sample.

Starkeya nomas sp. nov., a prosthecate and budding bacterium isolated from an immunocompromized patient.

Strain HF14-78462T is an environmental bacterium found in clinical samples from an immunocompromized patient in 2014 at Hospital Universitari i Politècnic La Fe (Valencia, Spain). Phenotypically, strain HF14-78462T cells were Gram-stain-negative, aerobic, non-spore forming and non-motile small rods which formed mucous and whitish-translucent colonies when incubated at 20-36 °C. Phylogenetic analyses based on the 16S rRNA genes and the whole genomes of closest sequenced relatives confirmed that strain HF14-78462T is affiliated with the genus Starkeya. The strain was oxidase, catalase and urease positive; but indole, lysine decarboxylase, ornithine decarboxylase and DNase negative, did not produce H2S and was able to utilize a wide variety of carbon sources including acetamide, adonitol, amygdalin, l-arabinose, citric acid, glucose, mannitol and melibiose. Unlike Starkeya novella and Starkeya koreensis, strain HF14-78462T failed to grow in thiosulphate-oxidizing media and had a narrower temperature growth range. Its genome was characterized by a size of 4.83 Mbp and a C+G content of 67.75 mol%. Major fatty acids were C18:1 ω7c, cyclo C19 : 0 and C16 : 0, its polar acids were diphosphatidylglycerol, phosphatidylcholine, phosphatidylethanolamine, phosphatidylglycerol and an aminophospholipid; while the ubiquinones were Q9 (1.8 %) and Q10 (98.2 %). Digital DNA-DNA hybridization values were 41 and 41.4 against S. novella and S. koreensis, respectively, while average nucleotide identity values were around 84 %. Phenotypic, average nucleotide identity and phylogenomic comparative studies suggest that strain HF14-78462T is a new representative of the genus Starkeya and the name Starkeya nomas sp. nov. is proposed. The type strain is HF14-78462T (=CECT 30124T=LMG 31874T).

MicroRNA-148b-3p and MicroRNA-25-3p Are Overexpressed in Fetuses with Late-Onset Fetal Growth Restriction.

OBJECTIVE: It was the aim of this study to describe a micro-RNA (miRNA) profile characteristic of late-onset fetal growth restriction (FGR) and to investigate the pathways involved in their biochemical action. METHODS: In this prospective study, 25 fetuses (16 normal and 9 with FGR [estimated fetal weight <10th centile plus cerebroplacental ratio <0.6765 multiples of the median]) were evaluated with Doppler ultrasound after 36 weeks. Afterwards, for every fetus, plasma from umbilical vein blood was collected at birth, miRNA was extracted, and full miRNA sequencing was performed. Subsequently, comparisons were done in order to obtain those miRNAs that were differentially expressed. RESULTS: The FGR fetuses expressed upregulation of two miRNAs: miR-25-3p and, especially, miR-148b-3p, a miRNA directly involved in Schwann cell migration, neuronal plasticity, and energy metabolism (p = 0.0072, p = 0.0013). CONCLUSIONS: FGR fetuses express a different miRNA profile, which includes overexpression of miR-25-3p and miR-148b-3p. This information might improve our understanding of the pathophysiological processes involved in late-onset FGR. Future validation and feasibility studies will be required to propose miRNAs as a valid tool in the diagnosis and management of FGR.

Tobacco rattle virus 16K silencing suppressor binds ARGONAUTE 4 and inhibits formation of RNA silencing complexes.

The cysteine-rich 16K protein of tobacco rattle virus (TRV), the type member of the genus Tobravirus, is known to suppress RNA silencing. However, the mechanism of action of the 16K suppressor is not well understood. In this study, we used a GFP-based sensor strategy and an Agrobacterium-mediated transient assay in Nicotiana benthamiana to show that 16K was unable to inhibit the activity of existing small interfering RNA (siRNA)- and microRNA (miRNA)-programmed RNA-induced silencing effector complexes (RISCs). In contrast, 16K efficiently interfered with de novo formation of miRNA- and siRNA-guided RISCs, thus preventing cleavage of target RNA. Interestingly, we found that transiently expressed endogenous miR399 and miR172 directed sequence-specific silencing of complementary sequences of viral origin. 16K failed to bind small RNAs, although it interacted with ARGONAUTE 4, as revealed by bimolecular fluorescence complementation and immunoprecipitation assays. Site-directed mutagenesis demonstrated that highly conserved cysteine residues within the N-terminal and central regions of the 16K protein are required for protein stability and/or RNA silencing suppression.

Urgent Implementation in a Hospital Setting of a Strategy To Rule Out Secondary Cases Caused by Imported Extensively Drug-Resistant Mycobacterium tuberculosis Strains at Diagnosis.

Current migratory movements require new strategies for rapidly tracking the transmission of high-risk imported Mycobacterium tuberculosis strains. Whole-genome sequencing (WGS) enables us to identify single-nucleotide polymorphisms (SNPs) and therefore design PCRs to track specific relevant strains. However, fast implementation of these strategies in the hospital setting is difficult because professionals working in diagnostics, molecular epidemiology, and genomics are generally at separate institutions. In this study, we describe the urgent implementation of a system that integrates genomics and molecular tools in a genuine high-risk epidemiological alert involving 2 independent importations of extensively drug resistant (XDR) and pre-XDR Beijing M. tuberculosis strains from Russia into Spain. Both cases involved commercial sex workers with long-standing tuberculosis (TB). The system was based on strain-specific PCRs tailored from WGS data that were transferred to the local node that was managing the epidemiological alert. The optimized tests were available for prospective implementation in the local node 33 working days after receiving the primary cultures of the XDR strains and were applied to all 42 new incident cases. An interpretable result was obtained in each case (directly from sputum for 27 stain-positive cases) and corresponded to the amplification profiles for strains other than the targeted pre-XDR and XDR strains, which made it possible to prospectively rule out transmission of these high-risk strains at diagnosis.

Two strawberry miR159 family members display developmental-specific expression patterns in the fruit receptacle and cooperatively regulate Fa-GAMYB.

• We have reported previously that the gibberellin (GA) content in strawberry receptacle is high, peaking at specific stages, pointing to a role of this hormone in fruit development. In Arabidopsis, miR159 levels are dependent on GA concentration. This prompted us to investigate the role of two members of the miR159 family and their putative strawberry target gene, GAMYB, in relation to changes in GA content during the course of fruit development. • The highest expression level of the two Fa-MIR159 genes was in the fruit's receptacle tissue, with dramatic changes observed throughout development. The lowest levels of total mature miR159 (a and b) were observed during the white stage of receptacle development, which was concurrent with the highest expression of Fa-GAMYB. A functional interaction between miR159 and Fa-GAMYB has been demonstrated in receptacle tissue. • The application of bioactive GA (i.e. GA(3) ) to strawberry plants caused the down-regulated expression of Fa-MIR159a, but the expression of Fa-MIR159b was not affected significantly. Clear discrepancies between Fa-MIR159b and mature Fa-miR159b levels were indicative of post-transcriptional regulation of Fa-MIR159b gene expression. • We propose that Fa-miR159a and Fa-miR159b interact with Fa-GAMYB during the course of strawberry receptacle development, and that they act in a cooperative fashion to respond, in part, to changes in GA endogenous levels.

Rapid Oxford Nanopore Technologies MinION Sequencing Workflow for Campylobacter jejuni Identification in Broilers on Site­A Proof-of-Concept Study

Campylobacter is recognised as one of the most important foodborne bacteria, with a worldwide health and socioeconomic impact. This bacterium is one of the most important zoonotic players in poultry, where efficient and fast detection methods are required. Current official culture methods for Campylobacter enumeration in poultry usually include >44 h of culture and >72 h for identification, thus requiring at least five working shifts (ISO/TS 10272-2:2017). Here, we have assembled a portable sequencing kit composed of the Bento Lab and the MinION and developed a workflow for on-site farm use that is able to detect and report the presence of Campylobacter from caecal samples in less than five hours from sampling time, as well as the relationship of Campylobacter with other caecal microbes. Beyond that, our workflow may offer a cost-effective and practical method of microbiologically monitoring poultry at the farm. These results would demonstrate the possibility of carrying out rapid on-site screening to monitor the health status of the poultry farm/flock during the production chain.

Profiling the Bladder Microbiota in Patients With Bladder Cancer.

Evidence suggests that microbiota may contribute to the pathogenesis of several diseases, including cancer. In the case of bladder cancer, preliminary studies have found alterations in the urinary microbiota of patients with urothelial carcinoma compared with healthy individuals. Conversely, the urinary microbiota differ between men and women, and it has been hypothesized that these differences are associated with the lower incidence of bladder cancers in women. The objective of this study was to characterize the bladder microbiota in paired samples of tumor and non-tumor mucosa of patients with malignant bladder neoplasia using next-generation sequencing. In addition, we aimed to study potential differences in microbial composition in tumor samples according to clinical and pathological variables, and to determine possible microbial profiles. We found significant differences in microbial richness at the genus level, with a higher richness observed in the non-tumor compared with the tumor mucosa. It was also shown that Actinobacteria were significantly more enriched in the non-tumor compared with the tumor mucosa (P = 0.014). In the multivariate analysis, we found significant differences in microbial composition according to tumor grade (P = 0.03 and 0.04 at the phylum and genus levels, respectively). Moreover, we detected a higher microbial richness in non-tumor vs. tumor tissues which agrees with the global assumption that microbial richness is an indicator of health. The greater abundance of members of the Actinobacteria phylum in the non-neoplastic bladder mucosa samples supports the hypothesis that a higher abundance of Actinomycetes is associated with a lower rate of bladder cancer in women and suggests a protective role for these microbiota. We detected a microbial profile that was enriched for Enterococcus in low-grade tumors. Finally, we identified the presence of two clusters in the microbial composition of the tumor mucosa samples, significantly enriched for the genera Barnesiella, Parabacteroides, Prevotella, Alistipes, and Lachnospiracea_incertae_sedis (Cluster 1), or Staphylococcus (Cluster 2). Further longitudinal studies are needed to assess the role of the bladder microbiota in carcinogenesis.

Evolutionary and Phenotypic Characterization of Two Spike Mutations in European Lineage 20E of SARS-CoV-2.

We have detected two mutations in the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at amino acid positions 1163 and 1167 that appeared independently in multiple transmission clusters and different genetic backgrounds. Furthermore, both mutations appeared together in a cluster of 1,627 sequences belonging to clade 20E. This cluster is characterized by 12 additional single nucleotide polymorphisms but no deletions. The available structural information on the S protein in the pre- and postfusion conformations predicts that both mutations confer rigidity, which could potentially decrease viral fitness. Accordingly, we observed reduced infectivity of this spike genotype relative to the ancestral 20E sequence in vitro, and the levels of viral RNA in nasopharyngeal swabs were not significantly higher. Furthermore, the mutations did not impact thermal stability or antibody neutralization by sera from vaccinated individuals but moderately reduce neutralization by convalescent-phase sera from the early stages of the pandemic. Despite multiple successful appearances of the two spike mutations during the first year of SARS-CoV-2 evolution, the genotype with both mutations was displaced upon the expansion of the 20I (Alpha) variant. The midterm fate of the genotype investigated was consistent with the lack of advantage observed in the clinical and experimental data. IMPORTANCE We observed repeated, independent emergence of mutations in the SARS-CoV-2 spike involving amino acids 1163 and 1167, within the HR2 functional motif. Conclusions derived from evolutionary and genomic diversity analysis suggest that the co-occurrence of both mutations might pose an advantage for the virus and therefore a threat to effective control of the epidemic. However, biological characterization, including in vitro experiments and analysis of clinical data, indicated no clear benefit in terms of stability or infectivity. In agreement with this, continuous epidemiological surveillance conducted months after the first observations revealed that both mutations did not successfully outcompete other variants and stopped circulating 9 months after their initial detection. Additionally, we evaluated the potential of both mutations to escape neutralizing antibodies, finding that the presence of these two mutations on their own is not likely to confer antibody escape. Our results provide an example of how newly emerged spike mutations can be assessed to better understand the risk posed by new variants and indicate that some spike mutations confer no clear advantage to the virus despite independently emerging multiple times and are eventually displaced by fitter variants.

The first wave of the COVID-19 epidemic in Spain was associated with early introductions and fast spread of a dominating genetic variant.

The coronavirus disease 2019 (COVID-19) pandemic has affected the world radically since 2020. Spain was one of the European countries with the highest incidence during the first wave. As a part of a consortium to monitor and study the evolution of the epidemic, we sequenced 2,170 samples, diagnosed mostly before lockdown measures. Here, we identified at least 500 introductions from multiple international sources and documented the early rise of two dominant Spanish epidemic clades (SECs), probably amplified by superspreading events. Both SECs were related closely to the initial Asian variants of SARS-CoV-2 and spread widely across Spain. We inferred a substantial reduction in the effective reproductive number of both SECs due to public-health interventions (Re < 1), also reflected in the replacement of SECs by a new variant over the summer of 2020. In summary, we reveal a notable difference in the initial genetic makeup of SARS-CoV-2 in Spain compared with other European countries and show evidence to support the effectiveness of lockdown measures in controlling virus spread, even for the most successful genetic variants.

Structural and genetic requirements for the biogenesis of tobacco rattle virus-derived small interfering RNAs.

In plants, small RNA-guided processes referred to as RNA silencing control gene expression and serve as an efficient antiviral mechanism. Plant viruses are inducers and targets of RNA silencing as infection involves the production of functional virus-derived small interfering RNAs (siRNAs). Here we investigate the structural and genetic components influencing the formation of Tobacco rattle virus (TRV)-derived siRNAs. TRV siRNAs are mostly 21 nucleotides in length and derive from positive and negative viral RNA strands, although TRV siRNAs of positive polarity are significantly more abundant. This asymmetry appears not to correlate with the presence of highly structured regions of single-stranded viral RNA. The Dicer-like enzyme DCL4, DCL3, or DCL2 targets, alone or in combination, viral templates to promote synthesis of siRNAs of both polarities from all regions of the viral genome. The heterogeneous distribution profile of TRV siRNAs reveals differential contributions throughout the TRV genome to siRNA formation. Indirect evidence suggests that DCL2 is responsible for production of a subset of siRNAs derived from the 3' end region of TRV. TRV siRNA biogenesis and antiviral silencing are strongly dependent on the combined activity of the host-encoded RNA-dependent RNA polymerases RDR1, RDR2, and RDR6, thus providing evidence that perfectly complementary double-stranded RNA serves as a substrate for siRNA production. We conclude that the overall composition of viral siRNAs in TRV-infected plants reflects the combined action of several interconnected pathways involving different DCL and RDR activities.

Draft Genome Sequences of 12 Monophasic Salmonella enterica subsp. enterica Serotype Typhimurium 1,4,[5],12:i:- Strains Isolated from Wild Griffon Vultures in Eastern Spain.

Monophasic Salmonella enterica subsp. enterica serovar Typhimurium is one of the most common zoonotic pathogens. Salmonella species reside in a wide variety of hosts, including wild animals. Thus, we report here the genome sequences of 12 monophasic S. Typhimurium strains isolated from healthy wild vultures to gain better insight into their epidemiology and host-pathogen interactions.

Genetic Structure of the Population of Pepino mosaic virus Infecting Tomato Crops in Spain.

ABSTRACT The population structure of Pepino mosaic virus (PepMV), which has caused severe epidemics in tomato in Spain since 2000, was analyzed. Isolates were characterized by the nucleotide sequence of the triple gene block and coat protein gene and, for a subset of isolates, a part of the RNA-dependent RNA polymerase gene. The full-length sequence of the genomic RNA of a Solanum muricatum isolate from Peru also was determined. In spite of high symptom diversity, the Spanish population of PepMV mostly comprised highly similar isolates belonging to the strain reported in Europe (European tomato strain), which has been the most prevalent genotype in Spain. The Spanish PepMV population was not structured spatially or temporally. Also, isolates highly similar to those from nontomato hosts from Peru (Peruvian strain) or to isolate US2 from the United States (US2 strain) were detected at lower frequency relative to the European strain. These two strains were detected in peninsular Spain only in 2004, but the Peruvian strain has been detected in the Canary Islands since 2000. These results suggest that PepMV was introduced into Spain more than once. Isolates from the Peruvian and US2 strains always were found in mixed infections with the European tomato strain, and interstrain recombinants were detected. The presence of different strains of the virus, and of recombinant isolates, should be considered for the development of control strategies based on genetic resistance.

Effect of Bifidobacterium breve on the Intestinal Microbiota of Coeliac Children on a Gluten Free Diet: A Pilot Study.

Coeliac disease (CD) is associated with alterations of the intestinal microbiota. Although several Bifidobacterium strains showed anti-inflammatory activity and prevention of toxic gliadin peptides generation in vitro, few data are available on their efficacy when administered to CD subjects. This study evaluated the effect of administration for three months of a food supplement based on two Bifidobacterium breve strains (B632 and BR03) to restore the gut microbial balance in coeliac children on a gluten free diet (GFD). Microbial DNA was extracted from faeces of 40 coeliac children before and after probiotic or placebo administration and 16 healthy children (Control group). Sequencing of the amplified V3-V4 hypervariable region of 16S rRNA gene as well as qPCR of Bidobacterium spp., Lactobacillus spp., Bacteroides fragilis group Clostridiumsensu stricto and enterobacteria were performed. The comparison between CD subjects and Control group revealed an alteration in the intestinal microbial composition of coeliacs mainly characterized by a reduction of the Firmicutes/Bacteroidetes ratio, of Actinobacteria and Euryarchaeota. Regarding the effects of the probiotic, an increase of Actinobacteria was found as well as a re-establishment of the physiological Firmicutes/Bacteroidetes ratio. Therefore, a three-month administration of B. breve strains helps in restoring the healthy percentage of main microbial components.

Direct squencing from the minimal number of DNA molecules needed to fill a 454 picotiterplate.

The large amount of DNA needed to prepare a library in next generation sequencing protocols hinders direct sequencing of small DNA samples. This limitation is usually overcome by the enrichment of such samples with whole genome amplification (WGA), mostly by multiple displacement amplification (MDA) based on φ29 polymerase. However, this technique can be biased by the GC content of the sample and is prone to the development of chimeras as well as contamination during enrichment, which contributes to undesired noise during sequence data analysis, and also hampers the proper functional and/or taxonomic assignments. An alternative to MDA is direct DNA sequencing (DS), which represents the theoretical gold standard in genome sequencing. In this work, we explore the possibility of sequencing the genome of Escherichia coli fs 24 from the minimum number of DNA molecules required for pyrosequencing, according to the notion of one-bead-one-molecule. Using an optimized protocol for DS, we constructed a shotgun library containing the minimum number of DNA molecules needed to fill a selected region of a picotiterplate. We gathered most of the reference genome extension with uniform coverage. We compared the DS method with MDA applied to the same amount of starting DNA. As expected, MDA yielded a sparse and biased read distribution, with a very high amount of unassigned and unspecific DNA amplifications. The optimized DS protocol allows unbiased sequencing to be performed from samples with a very small amount of DNA.

Silencing suppressor activity of the Tobacco rattle virus-encoded 16-kDa protein and interference with endogenous small RNA-guided regulatory pathways.

Higher plants use RNA silencing as a defense mechanism against viral infections, but viruses may encode suppressor proteins that counteract these defenses. Several virus-encoded suppressors also exert an inhibitory effect on endogenous small RNA regulatory pathways. Here we characterized the Tobacco rattle virus-encoded 16-kDa (TRV-16K) protein as a suppressor that blocked local RNA silencing induced by single (s)- and double-stranded (ds) RNA, indicating that TRV-16K interfered with a step in the silencing pathway downstream of dsRNA formation. The suppressor activity of TRV-16K was severely compromised by moderate to high dosages of dsRNA inducer. When silencing was locally triggered by ssRNA or low levels of dsRNA, silencing suppression by TRV-16K was associated with reduced accumulation of silencing-related siRNAs. TRV-16K also prevented partially cell-to-cell movement and systemic propagation of silencing but not transitive amplification of RNA silencing. We showed that neither TRV nor TRV-16K caused a global deregulation of the microRNA-regulatory pathway in Arabidopsis, suggesting that interference with microRNA biology was not a prerequisite for TRV, and probably many other plant viruses, to develop systemic infections in plants.