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.

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.

Morphological and Transcriptional Changes in Human Bone Marrow During Natural Plasmodium vivax Malaria Infections.

BACKGROUND: The presence of Plasmodium vivax malaria parasites in the human bone marrow (BM) is still controversial. However, recent data from a clinical case and experimental infections in splenectomized nonhuman primates unequivocally demonstrated the presence of parasites in this tissue. METHODS: In the current study, we analyzed BM aspirates of 7 patients during the acute attack and 42 days after drug treatment. RNA extracted from CD71+ cell suspensions was used for sequencing and transcriptomic analysis. RESULTS: We demonstrated the presence of parasites in all patients during acute infections. To provide further insights, we purified CD71+ BM cells and demonstrated dyserythropoiesis and inefficient erythropoiesis in all patients. In addition, RNA sequencing from 3 patients showed that genes related to erythroid maturation were down-regulated during acute infections, whereas immune response genes were up-regulated. CONCLUSIONS: This study thus shows that during P. vivax infections, parasites are always present in the BM and that such infections induced dyserythropoiesis and ineffective erythropoiesis. Moreover, infections induce transcriptional changes associated with such altered erythropoietic response, thus highlighting the importance of this hidden niche during natural infections.

Characterization of Acute HCV Infection and Transmission Networks in People Who Currently Inject Drugs in Catalonia: Usefulness of Dried Blood Spots.

BACKGROUND AND AIMS: Accurate identification of recent HCV infections is critical for tracing the extent and mechanisms of ongoing transmission. We aimed to validate dried blood spot (DBS) samples for the assessment of Hepatitis C virus (HCV) genetic diversity and to determine epidemiological parameters including incidence, determinants of acute infection, and phylogenetic clustering in people who inject drugs (PWID). APPROACH AND RESULTS: HCV nonstructural protein 5B next-generation sequencing was performed from plasma and/or DBS in 220 viremic PWID from the HepCdetect II study. No significant differences were found in consensus sequences or Shannon entropy (SE) intrahost diversity estimate between paired plasma/DBS specimens. SE values were used to identify acute infections with 93.3% sensitivity (95% CI, 0.81-1.06) and 95.0% specificity (95% CI, 0.88-1.02) in a set of well-defined controls. An acute HCV infection (either primary infection or reinfection) was detected in 13.5% of viremic participants and was associated with age ≤30 years (OR, 8.09), injecting less than daily (OR, 4.35), ≤5 years of injected drug use (OR, 3.43), sharing cocaine snorting straws (OR, 2.89), and being unaware of their HCV status (OR, 3.62). Annualized HCV incidence was estimated between 31 and 59/100 person-years. On phylogenetic analysis, 46.8% of viremic cases were part of a transmission pair or cluster; age ≤30 years (OR, 6.16), acute infection (OR, 5.73), and infection with subtype 1a (OR, 4.78) were independently associated with this condition. CONCLUSIONS: The results obtained from plasma and DBS characterize PWID with acute infection and those involved in ongoing HCV transmission and allow estimating incidence from cross-sectional data. This information is critical for the design and assessment of targeted harm reduction programs and test-and-treat interventions and to facilitate monitoring of HCV elimination in this key population.

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.

A likelihood ratio approach for identifying three-quarter siblings in genetic databases.

The detection of family relationships in genetic databases is of interest in various scientific disciplines such as genetic epidemiology, population and conservation genetics, forensic science, and genealogical research. Nowadays, screening genetic databases for related individuals forms an important aspect of standard quality control procedures. Relatedness research is usually based on an allele sharing analysis of identity by state (IBS) or identity by descent (IBD) alleles. Existing IBS/IBD methods mainly aim to identify first-degree relationships (parent-offspring or full siblings) and second degree (half-siblings, avuncular, or grandparent-grandchild) pairs. Little attention has been paid to the detection of in-between first and second-degree relationships such as three-quarter siblings (3/4S) who share fewer alleles than first-degree relationships but more alleles than second-degree relationships. With the progressively increasing sample sizes used in genetic research, it becomes more likely that such relationships are present in the database under study. In this paper, we extend existing likelihood ratio (LR) methodology to accurately infer the existence of 3/4S, distinguishing them from full siblings and second-degree relatives. We use bootstrap confidence intervals to express uncertainty in the LRs. Our proposal accounts for linkage disequilibrium (LD) by using marker pruning, and we validate our methodology with a pedigree-based simulation study accounting for both LD and recombination. An empirical genome-wide array data set from the GCAT Genomes for Life cohort project is used to illustrate the method.

Epigenetic footprint enables molecular risk stratification of hepatoblastoma with clinical implications.

BACKGROUND & AIMS: Hepatoblastoma (HB) is a rare disease. Nevertheless, it is the predominant pediatric liver cancer, with limited therapeutic options for patients with aggressive tumors. Herein, we aimed to uncover the mechanisms of HB pathobiology and to identify new biomarkers and therapeutic targets in a move towards precision medicine for patients with advanced HB. METHODS: We performed a comprehensive genomic, transcriptomic and epigenomic characterization of 159 clinically annotated samples from 113 patients with HB, using high-throughput technologies. RESULTS: We discovered a widespread epigenetic footprint of HB that includes hyperediting of the tumor suppressor BLCAP concomitant with a genome-wide dysregulation of RNA editing and the overexpression of mainly non-coding genes of the oncogenic 14q32 DLK1-DIO3 locus. By unsupervised analysis, we identified 2 epigenomic clusters (Epi-CA, Epi-CB) with distinct degrees of DNA hypomethylation and CpG island hypermethylation that are associated with the C1/C2/C2B transcriptomic subtypes. Based on these findings, we defined the first molecular risk stratification of HB (MRS-HB), which encompasses 3 main prognostic categories and improves the current clinical risk stratification approach. The MRS-3 category (28%), defined by strong 14q32 locus expression and Epi-CB methylation features, was characterized by CTNNB1 and NFE2L2 mutations, a progenitor-like phenotype and clinical aggressiveness. Finally, we identified choline kinase alpha as a promising therapeutic target for intermediate and high-risk HBs, as its inhibition in HB cell lines and patient-derived xenografts strongly abrogated tumor growth. CONCLUSIONS: These findings provide a detailed insight into the molecular features of HB and could be used to improve current clinical stratification approaches and to develop treatments for patients with HB. LAY SUMMARY: Hepatoblastoma is a rare childhood liver cancer that has been understudied. We have used cutting-edge technologies to expand our molecular knowledge of this cancer. Our biological findings can be used to improve clinical management and pave the way for the development of novel therapies for this cancer.

Multitrait genome association analysis identifies new susceptibility genes for human anthropometric variation in the GCAT cohort.

BACKGROUND: Heritability estimates have revealed an important contribution of SNP variants for most common traits; however, SNP analysis by single-trait genome-wide association studies (GWAS) has failed to uncover their impact. In this study, we applied a multitrait GWAS approach to discover additional factor of the missing heritability of human anthropometric variation. METHODS: We analysed 205 traits, including diseases identified at baseline in the GCAT cohort (Genomes For Life- Cohort study of the Genomes of Catalonia) (n=4988), a Mediterranean adult population-based cohort study from the south of Europe. We estimated SNP heritability contribution and single-trait GWAS for all traits from 15 million SNP variants. Then, we applied a multitrait-related approach to study genome-wide association to anthropometric measures in a two-stage meta-analysis with the UK Biobank cohort (n=336 107). RESULTS: Heritability estimates (eg, skin colour, alcohol consumption, smoking habit, body mass index, educational level or height) revealed an important contribution of SNP variants, ranging from 18% to 77%. Single-trait analysis identified 1785 SNPs with genome-wide significance threshold. From these, several previously reported single-trait hits were confirmed in our sample with LINC01432 (p=1.9×10-9) variants associated with male baldness, LDLR variants with hyperlipidaemia (ICD-9:272) (p=9.4×10-10) and variants in IRF4 (p=2.8×10-57), SLC45A2 (p=2.2×10-130), HERC2 (p=2.8×10-176), OCA2 (p=2.4×10-121) and MC1R (p=7.7×10-22) associated with hair, eye and skin colour, freckling, tanning capacity and sun burning sensitivity and the Fitzpatrick phototype score, all highly correlated cross-phenotypes. Multitrait meta-analysis of anthropometric variation validated 27 loci in a two-stage meta-analysis with a large British ancestry cohort, six of which are newly reported here (p value threshold <5×10-9) at ZRANB2-AS2, PIK3R1, EPHA7, MAD1L1, CACUL1 and MAP3K9. CONCLUSION: Considering multiple-related genetic phenotypes improve associated genome signal detection. These results indicate the potential value of data-driven multivariate phenotyping for genetic studies in large population-based cohorts to contribute to knowledge of complex traits.

Gene expression changes in blood RNA after swimming in a chlorinated pool.

Exposure to disinfection by-products (DBP) such as trihalomethanes (THM) in swimming pools has been linked to adverse health effects in humans, but their biological mechanisms are unclear. We evaluated short-term changes in blood gene expression of adult recreational swimmers after swimming in a chlorinated pool. Volunteers swam 40min in an indoor chlorinated pool. Blood samples were drawn and four THM (chloroform, bromodichloromethane, dibromochloromethane and bromoform) were measured in exhaled breath before and after swimming. Intensity of physical activity was measured as metabolic equivalents (METs). Gene expression in whole blood mRNA was evaluated using IlluminaHumanHT-12v3 Expression-BeadChip. Linear mixed models were used to evaluate the relationship between gene expression changes and THM exposure. Thirty-seven before-after pairs were analyzed. The median increase from baseline to after swimming were: 0.7 to 2.3 for MET, and 1.4 to 7.1µg/m3 for exhaled total THM (sum of the four THM). Exhaled THM increased on average 0.94µg/m3 per 1 MET. While 1643 probes were differentially expressed post-exposure. Of them, 189 were also associated with exhaled levels of individual/total THM or MET after False Discovery Rate. The observed associations with the exhaled THM were low to moderate (Log-fold change range: -0.17 to 0.15). In conclusion, we identified short-term gene expression changes associated with swimming in a pool that were minor in magnitude and their biological meaning was unspecific. The high collinearity between exhaled THM levels and intensity of physical activity precluded mutually adjusted models with both covariates. These exploratory results should be validated in future studies.

Cyclic AMP signaling restricts activation and promotes maturation and antioxidant defenses in astrocytes.

BACKGROUND: cAMP signaling produces dramatic changes in astrocyte morphology and physiology. However, its involvement in phenotype acquisition and the transcriptionally mediated mechanisms of action are largely unknown. RESULTS: Here we analyzed the global transcriptome of cultured astroglial cells incubated with activators of cAMP pathways. A bulk of astroglial transcripts, 6221 annotated genes, were differentially regulated by cAMP signaling. cAMP analogs strongly upregulated genes involved in typical functions of mature astrocytes, such as homeostatic control, metabolic and structural support to neurons, antioxidant defense and communication, whereas they downregulated a considerable number of proliferating and immaturity-related transcripts. Moreover, numerous genes typically activated in reactive cells, such as scar components and immunological mediators, were repressed by cAMP. GSEA analysis contrasting gene expression profiles with transcriptome signatures of acutely isolated astrocytes and in situ evaluation of protein levels in these cells showed that cAMP signaling conferred mature and in vivo-like transcriptional features to cultured astrocytes. CONCLUSIONS: These results indicate that cAMP signaling is a key pathway promoting astrocyte maturation and restricting their developmental and activation features. Therefore, a positive modulation of cAMP signaling may promote the normal state of differentiated astrocytes and favor the protection and function of neuronal networks.

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.

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.

Interleukin-16 is increased in obesity and alters adipogenesis and inflammation in vitro.

Introduction: Obesity is a chronic condition associated with low-grade inflammation mainly due to immune cell infiltration of white adipose tissue (WAT). WAT is distributed into two main depots: subcutaneous WAT (sWAT) and visceral WAT (vWAT), each with different biochemical features and metabolic roles. Proinflammatory cytokines including interleukin (IL)-16 are secreted by both adipocytes and infiltrated immune cells to upregulate inflammation. IL-16 has been widely studied in the peripheral proinflammatory immune response; however, little is known about its role in adipocytes in the context of obesity. Aim & Methods: We aimed to study the levels of IL-16 in WAT derived from sWAT and vWAT depots of humans with obesity and the role of this cytokine in palmitate-exposed 3T3-L1 adipocytes. Results: The results demonstrated that IL-16 expression was higher in vWAT compared with sWAT in individuals with obesity. In addition, IL-16 serum levels were higher in patients with obesity compared with normal-weight individuals, increased at 6 months after bariatric surgery, and at 12 months after surgery decreased to levels similar to before the intervention. Our in vitro models showed that IL-16 could modulate markers of adipogenesis (Pref1), lipid metabolism (Plin1, Cd36, and Glut4), fibrosis (Hif1a, Col4a, Col6a, and Vegf), and inflammatory signaling (IL6) during adipogenesis and in mature adipocytes. In addition, lipid accumulation and glycerol release assays suggested lipolysis alteration. Discussion: Our results suggest a potential role of IL-16 in adipogenesis, lipid and glucose homeostasis, fibrosis, and inflammation in an obesity context.

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.

Claudin-1 as a novel target gene induced in obesity and associated to inflammation, fibrosis, and cell differentiation.

OBJECTIVE: T lymphocytes from visceral and subcutaneous white adipose tissues (vWAT and sWAT, respectively) can have opposing roles in the systemic metabolic changes associated with obesity. However, few studies have focused on this subject. Claudin-1 (CLDN1) is a protein involved canonically in tight junctions and tissue paracellular permeability. We evaluated T-lymphocyte gene expression in vWAT and sWAT and in the whole adipose depots in human samples. METHODS: A Clariom D-based transcriptomic analysis was performed on T lymphocytes magnetically separated from vWAT and sWAT from patients with obesity (Cohort 1; N = 11). Expression of candidate genes resulting from that analysis was determined in whole WAT from individuals with and without obesity (Cohort 2; patients with obesity: N = 13; patients without obesity: N = 14). RESULTS: We observed transcriptional differences between T lymphocytes from sWAT compared with vWAT. Specifically, CLDN1 expression was found to be dramatically induced in vWAT T cells relative to those isolated from sWAT in patients with obesity. CLDN1 was also induced in obesity in vWAT and its expression correlates with genes involved in inflammation, fibrosis, and adipogenesis. CONCLUSION: These results suggest that CLDN1 is a novel marker induced in obesity and differentially expressed in T lymphocytes infiltrated in human vWAT as compared with sWAT. This protein may have a crucial role in the crosstalk between T lymphocytes and other adipose tissue cells and may contribute to inflammation, fibrosis, and alter homeostasis and promote metabolic disease in obesity.

Microarray and deep sequencing cross-platform analysis of the mirRNome and isomiR variation in response to epidermal growth factor.

BACKGROUND: Epidermal Growth Factor (EGF) plays an important function in the regulation of cell growth, proliferation, and differentiation by binding to its receptor (EGFR) and providing cancer cells with increased survival responsiveness. Signal transduction carried out by EGF has been extensively studied at both transcriptional and post-transcriptional levels. Little is known about the involvement of microRNAs (miRNAs) in the EGF signaling pathway. miRNAs have emerged as major players in the complex networks of gene regulation, and cancer miRNA expression studies have evidenced a direct involvement of miRNAs in cancer progression. RESULTS: In this study, we have used an integrative high content analysis approach to identify the specific miRNAs implicated in EGF signaling in HeLa cells as potential mediators of cancer mediated functions. We have used microarray and deep-sequencing technologies in order to obtain a global view of the EGF miRNA transcriptome with a robust experimental cross-validation. By applying a procedure based on Rankprod tests, we have delimited a solid set of EGF-regulated miRNAs. After validating regulated miRNAs by reverse transcription quantitative PCR, we have derived protein networks and biological functions from the predicted targets of the regulated miRNAs to gain insight into the potential role of miRNAs in EGF-treated cells. In addition, we have analyzed sequence heterogeneity due to editing relative to the reference sequence (isomiRs) among regulated miRNAs. CONCLUSIONS: We propose that the use of global genomic miRNA cross-validation derived from high throughput technologies can be used to generate more reliable datasets inferring more robust networks of co-regulated predicted miRNA target genes.

A myriad of miRNA variants in control and Huntington's disease brain regions detected by massively parallel sequencing.

Huntington disease (HD) is a neurodegenerative disorder that predominantly affects neurons of the forebrain. We have applied the Illumina massively parallel sequencing to deeply analyze the small RNA populations of two different forebrain areas, the frontal cortex (FC) and the striatum (ST) of healthy individuals and individuals with HD. More than 80% of the small-RNAs were annotated as microRNAs (miRNAs) in all samples. Deep sequencing revealed length and sequence heterogeneity (IsomiRs) for the vast majority of miRNAs. Around 80-90% of the miRNAs presented modifications in the 3'-terminus mainly in the form of trimming and/or as nucleotide addition variants, while the 5'-terminus of the miRNAs was specially protected from changes. Expression profiling showed strong miRNA and isomiR expression deregulation in HD, most being common to both FC and ST. The analysis of the upstream regulatory regions in co-regulated miRNAs suggests a role for RE1-Silencing Transcription Factor (REST) and P53 in miRNAs downregulation in HD. The putative targets of deregulated miRNAs and seed-region IsomiRs strongly suggest that their altered expression contributes to the aberrant gene expression in HD. Our results show that miRNA variability is a ubiquitous phenomenon in the adult human brain, which may influence gene expression in physiological and pathological conditions.

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.

Monitoring SARS-CoV-2 variant transitions using differences in diagnostic cycle threshold values of target genes.

Monitoring the emergence of new SARS-CoV-2 variants is important to detect potential risks of increased transmission or disease severity. We investigated the identification of SARS-CoV-2 variants from real-time reverse transcriptase polymerase chain reaction (RT-PCR) routine diagnostics data. Cycle threshold (Ct) values of positive samples were collected from April 2021 to January 2022 in the Northern Metropolitan Area of Barcelona (n = 15,254). Viral lineage identification from whole genome sequencing (WGS) was available for 4618 (30.3%) of these samples. Pairwise differences in the Ct values between gene targets (ΔCt) were analyzed for variants of concern or interest circulating in our area. A specific delay in the Ct of the N-gene compared to the RdRp-gene (ΔCtNR) was observed for Alpha, Delta, Eta and Omicron. Temporal differences in ΔCtNR correlated with the dynamics of viral replacement of Alpha by Delta and of Delta by Omicron according to WGS results. Using ΔCtNR, prediction of new variants of concern at early stages of circulation was achieved with high sensitivity and specificity (91.1% and 97.8% for Delta; 98.5% and 90.8% for Omicron). Thus, tracking population-wide trends in ΔCt values obtained from routine diagnostics testing in combination with WGS could be useful for real-time management and response to local epidemics.

A Polygenic Risk Score Based on a Cardioembolic Stroke Multitrait Analysis Improves a Clinical Prediction Model for This Stroke Subtype.

Background: Occult atrial fibrillation (AF) is one of the major causes of embolic stroke of undetermined source (ESUS). Knowing the underlying etiology of an ESUS will reduce stroke recurrence and/or unnecessary use of anticoagulants. Understanding cardioembolic strokes (CES), whose main cause is AF, will provide tools to select patients who would benefit from anticoagulants among those with ESUS or AF. We aimed to discover novel loci associated with CES and create a polygenetic risk score (PRS) for a more efficient CES risk stratification. Methods: Multitrait analysis of GWAS (MTAG) was performed with MEGASTROKE-CES cohort (n = 362,661) and AF cohort (n = 1,030,836). We considered significant variants and replicated those variants with MTAG p-value < 5 × 10-8 influencing both traits (GWAS-pairwise) with a p-value < 0.05 in the original GWAS and in an independent cohort (n = 9,105). The PRS was created with PRSice-2 and evaluated in the independent cohort. Results: We found and replicated eleven loci associated with CES. Eight were novel loci. Seven of them had been previously associated with AF, namely, CAV1, ESR2, GORAB, IGF1R, NEURL1, WIPF1, and ZEB2. KIAA1755 locus had never been associated with CES/AF, leading its index variant to a missense change (R1045W). The PRS generated has been significantly associated with CES improving discrimination and patient reclassification of a model with age, sex, and hypertension. Conclusion: The loci found significantly associated with CES in the MTAG, together with the creation of a PRS that improves the predictive clinical models of CES, might help guide future clinical trials of anticoagulant therapy in patients with ESUS or AF.