Analysis of blood and nasal epithelial transcriptomes to identify mechanisms associated with control of SARS-CoV-2 viral load in the upper respiratory tract.

OBJECTIVES: The amount of SARS-CoV-2 detected in the upper respiratory tract (URT viral load) is a key driver of transmission of infection. Current evidence suggests that mechanisms constraining URT viral load are different from those controlling lower respiratory tract viral load and disease severity. Understanding such mechanisms may help to develop treatments and vaccine strategies to reduce transmission. Combining mathematical modelling of URT viral load dynamics with transcriptome analyses we aimed to identify mechanisms controlling URT viral load. METHODS: COVID-19 patients were recruited in Spain during the first wave of the pandemic. RNA sequencing of peripheral blood and targeted NanoString nCounter transcriptome analysis of nasal epithelium were performed and gene expression analysed in relation to paired URT viral load samples collected within 15 days of symptom onset. Proportions of major immune cells in blood were estimated from transcriptional data using computational differential estimation. Weighted correlation network analysis (adjusted for cell proportions) and fixed transcriptional repertoire analysis were used to identify associations with URT viral load, quantified as standard deviations (z-scores) from an expected trajectory over time. RESULTS: Eighty-two subjects (50% female, median age 54 years (range 3-73)) with COVID-19 were recruited. Paired URT viral load samples were available for 16 blood transcriptome samples, and 17 respiratory epithelial transcriptome samples. Natural Killer (NK) cells were the only blood cell type significantly correlated with URT viral load z-scores (r = -0.62, P = 0.010). Twenty-four blood gene expression modules were significantly correlated with URT viral load z-score, the most significant being a module of genes connected around IFNA14 (Interferon Alpha-14) expression (r = -0.60, P = 1e-10). In fixed repertoire analysis, prostanoid-related gene expression was significantly associated with higher viral load. In nasal epithelium, only GNLY (granulysin) gene expression showed significant negative correlation with viral load. CONCLUSIONS: Correlations between the transcriptional host response and inter-individual variations in SARS-CoV-2 URT viral load, revealed many molecular mechanisms plausibly favouring or constraining viral replication. Existing evidence corroborates many of these mechanisms, including likely roles for NK cells, granulysin, prostanoids and interferon alpha-14. Inhibition of prostanoid production and administration of interferon alpha-14 may be attractive transmission-blocking interventions.

Unraveling The Effects of DICER1 Overexpression on Immune-Related Genes Expression in Mesenchymal Stromal/Stem Cells: Insights for Therapeutic Applications.

OBJECTIVE: The immunoregulatory properties of mesenchymal stromal/stem cells (MSCs) bring a promise for the treatment of inflammatory diseases. However, their ability to suppress the immune system is unstable. To enhance their effectiveness against immune responses, it may be necessary to manipulate MSCs. Although some dsRNA transcripts come from invading viruses, the majority of dsRNA has an endogenous origin and is known as endo-siRNA. DICER1 is a ribonuclease protein that can generate small RNAs to modulate gene expression at the post-transcriptional level. We aimed to evaluate the expression of several immune-related genes at mRNA and protein levels in MSCs overexpressing DICER1 exogenously. MATERIALS AND METHODS: In this comparative transcriptomic experimental study, the adipose-derived MSCs (Ad-MSCs) were transfected using the pCAGGS-Flag-hsDicer vector for the DICER1 overexpression. Following the RNA extraction, mRNA expression level of DICER1 and several inflammatory cytokines were examined. We performed a relative real-time polymerase chain reaction (PCR) assay and transcriptome analysis between two groups including DICER1- transfected MSCs and control MSCs. Moreover, media from the transfected MSCs were evaluated for various interferon response factors by ELISA. RESULTS: The overexpression of DICER1 is associated with a significant increase in the mRNA expression level of COX-2, DDX-58, IFIH1, MYD88, RNase L, TLR3/4, and TDO2 genes and a downregulation of the TSG-6 gene in MSCs. Moreover, the expression levels of IL-1, 6, 8, 17, 18, CCL2, INF-γ, TGF-ß, and TNF-α were higher in the DICER1-transfected MSCs group. CONCLUSION: It seems that the ectopic expression of DICER1 in Ad-MSCs is linked to alterations in the expression level of immune-related genes. It is suggested that the manipulation of immune-related pathways in MSCs via the Dicer1 overexpression could facilitate the development of MSCs with distinct immunoregulatory phenotypes.

Quantitative trait and transcriptome analysis of genetic complexity underpinning cardiac interatrial septation in mice using an advanced intercross line.

Unlike single-gene mutations leading to Mendelian conditions, common human diseases are likely to be emergent phenomena arising from multilayer, multiscale, and highly interconnected interactions. Atrial and ventricular septal defects are the most common forms of cardiac congenital anomalies in humans. Atrial septal defects (ASD) show an open communication between the left and right atria postnatally, potentially resulting in serious hemodynamic consequences if untreated. A milder form of atrial septal defect, patent foramen ovale (PFO), exists in about one-quarter of the human population, strongly associated with ischaemic stroke and migraine. The anatomic liabilities and genetic and molecular basis of atrial septal defects remain unclear. Here, we advance our previous analysis of atrial septal variation through quantitative trait locus (QTL) mapping of an advanced intercross line (AIL) established between the inbred QSi5 and 129T2/SvEms mouse strains, that show extremes of septal phenotypes. Analysis resolved 37 unique septal QTL with high overlap between QTL for distinct septal traits and PFO as a binary trait. Whole genome sequencing of parental strains and filtering identified predicted functional variants, including in known human congenital heart disease genes. Transcriptome analysis of developing septa revealed downregulation of networks involving ribosome, nucleosome, mitochondrial, and extracellular matrix biosynthesis in the 129T2/SvEms strain, potentially reflecting an essential role for growth and cellular maturation in septal development. Analysis of variant architecture across different gene features, including enhancers and promoters, provided evidence for the involvement of non-coding as well as protein-coding variants. Our study provides the first high-resolution picture of genetic complexity and network liability underlying common congenital heart disease, with relevance to human ASD and PFO.

Genomic investigations of unexplained acute hepatitis in children.

Since its first identification in Scotland, over 1,000 cases of unexplained paediatric hepatitis in children have been reported worldwide, including 278 cases in the UK1. Here we report an investigation of 38 cases, 66 age-matched immunocompetent controls and 21 immunocompromised comparator participants, using a combination of genomic, transcriptomic, proteomic and immunohistochemical methods. We detected high levels of adeno-associated virus 2 (AAV2) DNA in the liver, blood, plasma or stool from 27 of 28 cases. We found low levels of adenovirus (HAdV) and human herpesvirus 6B (HHV-6B) in 23 of 31 and 16 of 23, respectively, of the cases tested. By contrast, AAV2 was infrequently detected and at low titre in the blood or the liver from control children with HAdV, even when profoundly immunosuppressed. AAV2, HAdV and HHV-6 phylogeny excluded the emergence of novel strains in cases. Histological analyses of explanted livers showed enrichment for T cells and B lineage cells. Proteomic comparison of liver tissue from cases and healthy controls identified increased expression of HLA class 2, immunoglobulin variable regions and complement proteins. HAdV and AAV2 proteins were not detected in the livers. Instead, we identified AAV2 DNA complexes reflecting both HAdV-mediated and HHV-6B-mediated replication. We hypothesize that high levels of abnormal AAV2 replication products aided by HAdV and, in severe cases, HHV-6B may have triggered immune-mediated hepatic disease in genetically and immunologically predisposed children.

eQTL Colocalization Analyses Identify NTN4 as a Candidate Breast Cancer Risk Gene.

Breast cancer genome-wide association studies (GWASs) have identified 150 genomic risk regions containing more than 13,000 credible causal variants (CCVs). The CCVs are predominantly noncoding and enriched in regulatory elements. However, the genes underlying breast cancer risk associations are largely unknown. Here, we used genetic colocalization analysis to identify loci at which gene expression could potentially explain breast cancer risk phenotypes. Using data from the Breast Cancer Association Consortium (BCAC) and quantitative trait loci (QTL) from the Genotype-Tissue Expression (GTEx) project and The Cancer Genome Project (TCGA), we identify shared genetic relationships and reveal novel associations between cancer phenotypes and effector genes. Seventeen genes, including NTN4, were identified as potential mediators of breast cancer risk. For NTN4, we showed the rs61938093 CCV at this region was located within an enhancer element that physically interacts with the NTN4 promoter, and the risk allele reduced NTN4 promoter activity. Furthermore, knockdown of NTN4 in breast cells increased cell proliferation in vitro and tumor growth in vivo. These data provide evidence linking risk-associated variation to genes that may contribute to breast cancer predisposition.

Chromatin interactome mapping at 139 independent breast cancer risk signals.

BACKGROUND: Genome-wide association studies have identified 196 high confidence independent signals associated with breast cancer susceptibility. Variants within these signals frequently fall in distal regulatory DNA elements that control gene expression. RESULTS: We designed a Capture Hi-C array to enrich for chromatin interactions between the credible causal variants and target genes in six human mammary epithelial and breast cancer cell lines. We show that interacting regions are enriched for open chromatin, histone marks for active enhancers, and transcription factors relevant to breast biology. We exploit this comprehensive resource to identify candidate target genes at 139 independent breast cancer risk signals and explore the functional mechanism underlying altered risk at the 12q24 risk region. CONCLUSIONS: Our results demonstrate the power of combining genetics, computational genomics, and molecular studies to rationalize the identification of key variants and candidate target genes at breast cancer GWAS signals.

Non-coding RNAs underlie genetic predisposition to breast cancer.

BACKGROUND: Genetic variants identified through genome-wide association studies (GWAS) are predominantly non-coding and typically attributed to altered regulatory elements such as enhancers and promoters. However, the contribution of non-coding RNAs to complex traits is not clear. RESULTS: Using targeted RNA sequencing, we systematically annotated multi-exonic non-coding RNA (mencRNA) genes transcribed from 1.5-Mb intervals surrounding 139 breast cancer GWAS signals and assessed their contribution to breast cancer risk. We identify more than 4000 mencRNA genes and show their expression distinguishes normal breast tissue from tumors and different breast cancer subtypes. Importantly, breast cancer risk variants, identified through genetic fine-mapping, are significantly enriched in mencRNA exons, but not the promoters or introns. eQTL analyses identify mencRNAs whose expression is associated with risk variants. Furthermore, chromatin interaction data identify hundreds of mencRNA promoters that loop to regions that contain breast cancer risk variants. CONCLUSIONS: We have compiled the largest catalog of breast cancer-associated mencRNAs to date and provide evidence that modulation of mencRNAs by GWAS variants may provide an alternative mechanism underlying complex traits.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes.

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

Hsp70 inhibitors: Implications for the treatment of colorectal cancer.

Colorectal cancer (CRC) is one of the most common malignancies in the world. Despite intensive advances in diagnosis and treatment of CRC, it is yet one of the leading cause of cancer related morbidity and mortality. Therefore, there is an urgent medical need for alternative therapeutic approaches to treat CRC. The 70 kDa heat shock proteins (Hsp70s) are a family of evolutionary conserved heat shock proteins, which play an important role in cell homeostasis and survival. They overexpress in various types of malignancy including CRC and are typically accompanied with poor prognosis. Hence, inhibition of Hsp70 may be considered as a striking chemotherapeutic avenue. This review summarizes the current knowledge on the progress made so far to discover compounds, which target the Hsp70 family, with particular emphasis on their efficacy in treatment of CRC. We also briefly explain the induction of Hsp70 as a strategy to prevent CRC.

Long Noncoding RNAs CUPID1 and CUPID2 Mediate Breast Cancer Risk at 11q13 by Modulating the Response to DNA Damage.

Breast cancer risk is strongly associated with an intergenic region on 11q13. We have previously shown that the strongest risk-associated SNPs fall within a distal enhancer that regulates CCND1. Here, we report that, in addition to regulating CCND1, this enhancer regulates two estrogen-regulated long noncoding RNAs, CUPID1 and CUPID2. We provide evidence that the risk-associated SNPs are associated with reduced chromatin looping between the enhancer and the CUPID1 and CUPID2 bidirectional promoter. We further show that CUPID1 and CUPID2 are predominantly expressed in hormone-receptor-positive breast tumors and play a role in modulating pathway choice for the repair of double-strand breaks. These data reveal a mechanism for the involvement of this region in breast cancer.

Evidence that the 5p12 Variant rs10941679 Confers Susceptibility to Estrogen-Receptor-Positive Breast Cancer through FGF10 and MRPS30 Regulation.

Genome-wide association studies (GWASs) have revealed increased breast cancer risk associated with multiple genetic variants at 5p12. Here, we report the fine mapping of this locus using data from 104,660 subjects from 50 case-control studies in the Breast Cancer Association Consortium (BCAC). With data for 3,365 genotyped and imputed SNPs across a 1 Mb region (positions 44,394,495-45,364,167; NCBI build 37), we found evidence for at least three independent signals: the strongest signal, consisting of a single SNP rs10941679, was associated with risk of estrogen-receptor-positive (ER+) breast cancer (per-g allele OR ER+ = 1.15; 95% CI 1.13-1.18; p = 8.35 × 10-30). After adjustment for rs10941679, we detected signal 2, consisting of 38 SNPs more strongly associated with ER-negative (ER-) breast cancer (lead SNP rs6864776: per-a allele OR ER- = 1.10; 95% CI 1.05-1.14; p conditional = 1.44 × 10-12), and a single signal 3 SNP (rs200229088: per-t allele OR ER+ = 1.12; 95% CI 1.09-1.15; p conditional = 1.12 × 10-05). Expression quantitative trait locus analysis in normal breast tissues and breast tumors showed that the g (risk) allele of rs10941679 was associated with increased expression of FGF10 and MRPS30. Functional assays demonstrated that SNP rs10941679 maps to an enhancer element that physically interacts with the FGF10 and MRPS30 promoter regions in breast cancer cell lines. FGF10 is an oncogene that binds to FGFR2 and is overexpressed in ∼10% of human breast cancers, whereas MRPS30 plays a key role in apoptosis. These data suggest that the strongest signal of association at 5p12 is mediated through coordinated activation of FGF10 and MRPS30, two candidate genes for breast cancer pathogenesis.

Fine scale mapping of the 17q22 breast cancer locus using dense SNPs, genotyped within the Collaborative Oncological Gene-Environment Study (COGs).

Genome-wide association studies have found SNPs at 17q22 to be associated with breast cancer risk. To identify potential causal variants related to breast cancer risk, we performed a high resolution fine-mapping analysis that involved genotyping 517 SNPs using a custom Illumina iSelect array (iCOGS) followed by imputation of genotypes for 3,134 SNPs in more than 89,000 participants of European ancestry from the Breast Cancer Association Consortium (BCAC). We identified 28 highly correlated common variants, in a 53 Kb region spanning two introns of the STXBP4 gene, that are strong candidates for driving breast cancer risk (lead SNP rs2787486 (OR = 0.92; CI 0.90-0.94; P = 8.96 × 10(-15))) and are correlated with two previously reported risk-associated variants at this locus, SNPs rs6504950 (OR = 0.94, P = 2.04 × 10(-09), r(2) = 0.73 with lead SNP) and rs1156287 (OR = 0.93, P = 3.41 × 10(-11), r(2) = 0.83 with lead SNP). Analyses indicate only one causal SNP in the region and several enhancer elements targeting STXBP4 are located within the 53 kb association signal. Expression studies in breast tumor tissues found SNP rs2787486 to be associated with increased STXBP4 expression, suggesting this may be a target gene of this locus.

Breast cancer risk variants at 6q25 display different phenotype associations and regulate ESR1, RMND1 and CCDC170.

We analyzed 3,872 common genetic variants across the ESR1 locus (encoding estrogen receptor α) in 118,816 subjects from three international consortia. We found evidence for at least five independent causal variants, each associated with different phenotype sets, including estrogen receptor (ER(+) or ER(-)) and human ERBB2 (HER2(+) or HER2(-)) tumor subtypes, mammographic density and tumor grade. The best candidate causal variants for ER(-) tumors lie in four separate enhancer elements, and their risk alleles reduce expression of ESR1, RMND1 and CCDC170, whereas the risk alleles of the strongest candidates for the remaining independent causal variant disrupt a silencer element and putatively increase ESR1 and RMND1 expression.

Investigation of association between PFO complicated by cryptogenic stroke and a common variant of the cardiac transcription factor GATA4.

Patent foramen ovale (PFO) is associated with clinical conditions including cryptogenic stroke, migraine and varicose veins. Data from studies in humans and mouse suggest that PFO and the secundum form of atrial septal defect (ASDII) exist in an anatomical continuum of septal dysmorphogenesis with a common genetic basis. Mutations in multiple members of the evolutionarily conserved cardiac transcription factor network, including GATA4, cause or predispose to ASDII and PFO. Here, we assessed whether the most prevalent variant of the GATA4 gene, S377G, was significantly associated with PFO or ASD. Our analysis of world indigenous populations showed that GATA4 S377G was largely Caucasian-specific, and so subjects were restricted to those of Caucasian descent. To select for patients with larger PFO, we limited our analysis to those with cryptogenic stroke in which PFO was a subsequent finding. In an initial study of Australian subjects, we observed a weak association between GATA4 S377G and PFO/Stroke relative to Caucasian controls in whom ASD and PFO had been excluded (OR = 2.16; p = 0.02). However, in a follow up study of German Caucasians no association was found with either PFO or ASD. Analysis of combined Australian and German data confirmed the lack of a significant association. Thus, the common GATA4 variant S377G is likely to be relatively benign in terms of its participation in CHD and PFO/Stroke.

Lack of beneficial effects of granulocyte colony-stimulating factor in patients with subacute myocardial infarction undergoing late revascularization: a double-blind, randomized, placebo-controlled clinical trial.

BACKGROUND: Early clinical studies have suggested that administration of granulocyte-colony stimulating factor (G-CSF) may improve the clinical condition of patients suffering from myocardial infarction (MI).This prospective, randomized, double-blind, placebo-controlled single-centre trial aims to assess the safety and clinical efficacy of G-CSF administration in patients with subacute MI and impaired LV function undergoing delayed primary percutaneous coronary intervention (PCI). METHODS: A total of 16 patients (13 men, mean age 51 years) with subacute ST-segment elevation MI and a left ventricular (LV) ejection fraction (EF) of less than 45% at baseline who underwent late revascularization, were included in the study. Patients were randomized in a double-blind fashion to receive either G-CSF (at a dose of 10 microg/kg body weight) or placebo for five consecutive days. End points consisted of assessment of safety parameters as well as changes of global and regional myocardial function from baseline until six months following PCI. RESULTS: G-CSF administration resulted in a significant mobilization of different cell populations (four-fold increase in WBC count and a six-fold increase in CD34+ cells). G-CSF treatment was well tolerated in most patients and no major adverse cardiac events or severe G-CSF-related side effects were identified during hospitalization and at follow-up. No significant differences were observed between the G-CSF and placebo groups regarding global and regional myocardial function parameters. CONCLUSION: G-CSF administration is safe, but not effective, in improving impaired LV functional parameters in patients with subacute MI who had an impaired baseline EF of less than 45%.