Genetic architecture of microRNA expression and its link to complex diseases in the Japanese population

Common and rare genetic variants predisposing females to unexplained recurrent pregnancy loss

Recurrent pregnancy loss (RPL) is a major reproductive health issue with multifactorial causes, affecting 2.6% of all pregnancies worldwide. Nearly half of the RPL cases lack clinically identifiable causes (e.g., antiphospholipid syndrome, uterine anomalies, and parental chromosomal abnormalities), referred to as unexplained RPL (uRPL). Here, we perform a genome-wide association study focusing on uRPL in 1,728 cases and 24,315 female controls of Japanese ancestry. We detect significant associations in the major histocompatibility complex (MHC) region at 6p21 (lead variant=rs9263738; P = 1.4 × 10-10; odds ratio [OR] = 1.51 [95% CI: 1.33-1.72]; risk allele frequency = 0.871). The MHC associations are fine-mapped to the classical HLA alleles, HLA-C*12:02, HLA-B*52:01, and HLA-DRB1*15:02 (P = 1.1 × 10-10, 1.5 × 10-10, and 1.2 × 10-9, respectively), which constitute a population-specific common long-range haplotype with a protective effect (P = 2.8 × 10-10; OR = 0.65 [95% CI: 0.57-0.75]; haplotype frequency=0.108). Genome-wide copy-number variation (CNV) calling demonstrates rare predicted loss-of-function (pLoF) variants of the cadherin-11 gene (CDH11) conferring the risk of uRPL (P = 1.3 × 10-4; OR = 3.29 [95% CI: 1.78-5.76]). Our study highlights the importance of reproductive immunology and rare variants in the uRPL etiology.

microRNA-guided immunity against respiratory virus infection in human and mouse lung cells

Viral infectivity depends on multiple factors. Recent studies showed that the interaction between viral RNAs and endogenous microRNAs (miRNAs) regulates viral infectivity; viral RNAs function as a sponge of endogenous miRNAs and result in upregulation of its original target genes, while endogenous miRNAs target viral RNAs directly and result in repression of viral gene expression. In this study, we analyzed the possible interaction between parainfluenza virus RNA and endogenous miRNAs in human and mouse lungs. We showed that the parainfluenza virus can form base pairs with human miRNAs abundantly than mouse miRNAs. Furthermore, we analyzed that the sponge effect of endogenous miRNAs on viral RNAs may induce the upregulation of transcription regulatory factors. Then, we performed RNA-sequence analysis and observed the upregulation of transcription regulatory factors in the early stages of parainfluenza virus infection. Our studies showed how the differential expression of endogenous miRNAs in lungs could contribute to respiratory virus infection and species- or tissue-specific mechanisms and common mechanisms could be conserved in humans and mice and regulated by miRNAs during viral infection.

Body mass index stratification optimizes polygenic prediction of type 2 diabetes in cross-biobank analyses

Type 2 diabetes (T2D) shows heterogeneous body mass index (BMI) sensitivity. Here, we performed stratification based on BMI to optimize predictions for BMI-related diseases. We obtained BMI-stratified datasets using data from more than 195,000 individuals (nT2D = 55,284) from BioBank Japan (BBJ) and UK Biobank. T2D heritability in the low-BMI group was greater than that in the high-BMI group. Polygenic predictions of T2D toward low-BMI targets had pseudo-R2 values that were more than 22% higher than BMI-unstratified targets. Polygenic risk scores (PRSs) from low-BMI discovery outperformed PRSs from high BMI, while PRSs from BMI-unstratified discovery performed best. Pathway-specific PRSs demonstrated the biological contributions of pathogenic pathways. Low-BMI T2D cases showed higher rates of neuropathy and retinopathy. Combining BMI stratification and a method integrating cross-population effects, T2D predictions showed greater than 37% improvements over unstratified-matched-population prediction. We replicated findings in the Tohoku Medical Megabank (n = 26,000) and the second BBJ cohort (n = 33,096). Our findings suggest that target stratification based on existing traits can improve the polygenic prediction of heterogeneous diseases.

Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.

COLOCdb: a comprehensive resource for multi-model colocalization of complex traits

Large-scale genome-wide association studies (GWAS) have provided profound insights into complex traits and diseases. Yet, deciphering the fine-scale molecular mechanisms of how genetic variants manifest to cause the phenotypes remains a daunting task. Here, we present COLOCdb (https://ngdc.cncb.ac.cn/colocdb), a comprehensive genetic colocalization database by integrating more than 3000 GWAS summary statistics and 13 types of xQTL to date. By employing two representative approaches for the colocalization analysis, COLOCdb deposits results from three key components: (i) GWAS-xQTL, pair-wise colocalization between GWAS loci and different types of xQTL, (ii) GWAS-GWAS, pair-wise colocalization between the trait-associated genetic loci from GWASs and (iii) xQTL-xQTL, pair-wise colocalization between the genetic loci associated with molecular phenotypes in xQTLs. These results together represent the most comprehensive colocalization analysis, which also greatly expands the list of shared variants with genetic pleiotropy. We expect that COLOCdb can serve as a unique and useful resource in advancing the discovery of new biological mechanisms and benefit future functional studies.

Analysis of gut microbiome, host genetics, and plasma metabolites reveals gut microbiome-host interactions in the Japanese population

Interaction between the gut microbiome and host plays a key role in human health. Here, we perform a metagenome shotgun-sequencing-based analysis of Japanese participants to reveal associations between the gut microbiome, host genetics, and plasma metabolome. A genome-wide association study (GWAS) for microbial species (n = 524) identifies associations between the PDE1C gene locus and Bacteroides intestinalis and between TGIF2 and TGIF2-RAB5IF gene loci and Bacteroides acidifiaciens. In a microbial gene ortholog GWAS, agaE and agaS, which are related to the metabolism of carbohydrates forming the blood group A antigen, are associated with blood group A in a manner depending on the secretor status determined by the East Asian-specific FUT2 variant. A microbiome-metabolome association analysis (n = 261) identifies associations between bile acids and microbial features such as bile acid metabolism gene orthologs including bai and 7β-hydroxysteroid dehydrogenase. Our publicly available data will be a useful resource for understanding gut microbiome-host interactions in an underrepresented population.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

Systems genetics identifies miRNA-mediated regulation of host response in COVID-19

BACKGROUND

Individuals infected with SARS-CoV-2 vary greatly in their disease severity, ranging from asymptomatic infection to severe disease. The regulation of gene expression is an important mechanism in the host immune response and can modulate the outcome of the disease. miRNAs play important roles in post-transcriptional regulation with consequences on downstream molecular and cellular host immune response processes. The nature and magnitude of miRNA perturbations associated with blood phenotypes and intensive care unit (ICU) admission in COVID-19 are poorly understood.

RESULTS

We combined multi-omics profiling-genotyping, miRNA and RNA expression, measured at the time of hospital admission soon after the onset of COVID-19 symptoms-with phenotypes from electronic health records to understand how miRNA expression contributes to variation in disease severity in a diverse cohort of 259 unvaccinated patients in Abu Dhabi, United Arab Emirates. We analyzed 62 clinical variables and expression levels of 632 miRNAs measured at admission and identified 97 miRNAs associated with 8 blood phenotypes significantly associated with later ICU admission. Integrative miRNA-mRNA cross-correlation analysis identified multiple miRNA-mRNA-blood endophenotype associations and revealed the effect of miR-143-3p on neutrophil count mediated by the expression of its target gene BCL2. We report 168 significant cis-miRNA expression quantitative trait loci, 57 of which implicate miRNAs associated with either ICU admission or a blood endophenotype.

CONCLUSIONS

This systems genetics study has given rise to a genomic picture of the architecture of whole blood miRNAs in unvaccinated COVID-19 patients and pinpoints post-transcriptional regulation as a potential mechanism that impacts blood traits underlying COVID-19 severity. The results also highlight the impact of host genetic regulatory control of miRNA expression in early stages of COVID-19 disease.

Reconstruction of the personal information from human genome reads in gut metagenome sequencing data

Human DNA present in faecal samples can result in a small number of human reads in gut shotgun metagenomic sequencing data. However, it is presently unclear how much personal information can be reconstructed from such reads, and this has not been quantitatively evaluated. Such a quantitative evaluation is necessary to clarify the ethical concerns related to data sharing and to enable efficient use of human genetic information in stool samples, such as for research and forensics. Here we used genomic approaches to reconstruct personal information from the faecal metagenomes of 343 Japanese individuals with associated human genotype data. Genetic sex could be accurately predicted based on the sequencing depth of sex chromosomes for 97.3% of the samples. Individuals could be re-identified from the matched genotype data based on human reads recovered from the faecal metagenomic data with 93.3% sensitivity using a likelihood score-based method. This method also enabled us to predict the ancestries of 98.3% of the samples. Finally, we performed ultra-deep shotgun metagenomic sequencing of five faecal samples as well as whole-genome sequencing of blood samples. Using genotype-calling approaches, we demonstrated that the genotypes of both common and rare variants could be reconstructed from faecal samples. This included clinically relevant variants. Our approach can be used to quantify personal information contained within gut metagenome data.

The 8q24 region hosts miRNAs altered in biospecimens of colorectal and bladder cancer patients

BACKGROUND

The 8q24 locus is enriched in cancer-associated polymorphisms and, despite containing relatively few protein-coding genes, it hosts the MYC oncogene and other genetic elements connected to tumorigenesis, including microRNAs (miRNAs). Research on miRNAs may provide insights into the transcriptomic regulation of this multiple cancer-associated region.

MATERIAL AND METHODS

We profiled all miRNAs located in the 8q24 region in 120 colorectal cancer (CRC) patients and 80 controls. miRNA profiling was performed on cancer/non-malignant adjacent mucosa, stool, and plasma extracellular vesicles (EVs), and the results validated with The Cancer Genome Atlas (TCGA) data. To verify if the 8q24-annotated miRNAs altered in CRC were dysregulated in other cancers and biofluids, we evaluated their levels in bladder cancer (BC) cases from the TCGA dataset and in urine and plasma EVs from a set of BC cases and healthy controls.

RESULTS

Among the detected mature miRNAs in the region, 12 were altered between CRC and adjacent mucosa (adj. p < 0.05). Five and four miRNAs were confirmed as dysregulated in the CRC and BC TCGA dataset, respectively. A co-expression analysis of tumor/adjacent tissue data from the CRC group revealed a correlation between the dysregulated miRNAs and CRC-related genes (PVT1 and MYC) annotated in 8q24 region. miR-30d-5p and miR-151a-3p, altered in CRC tissue, were also dysregulated in stool of CRC patients and urine of BC cases, respectively. Functional enrichment of dysregulated miRNA target genes highlighted terms related to TP53-mediated cell cycle control.

CONCLUSIONS

Altered expression of 8q24-annotated miRNAs may be relevant for the initiation and/or progression of cancer.

An atlas of expression quantitative trait loci of microRNA in longissimus muscle of eight-way crossbred pigs

MicroRNAs (miRNAs) are key regulators of myocyte development and traits, yet insight into the genetic basis of variation in miRNA expression is still limited. Here, we present a systematic analysis of expression quantitative trait loci (eQTL) for miRNA profiling in longissimus muscle of pigs from an eight-breed crossed heterogeneous population. By integrating whole-genome sequencing and miRNAomics data, we map 54 cis- and 292 trans-eQTLs at high resolution that are associated with the expression of 54 and 92 miRNAs, respectively. Twenty-three trans-acting loci are identified to affect the expression of nine myomiRs (known muscle-specific miRNAs). MiRNAs in mammalian conserved miRNA clusters are found to be subjected to regulation by shared cis-eQTLs, while the expression of mature miRNA-5p/-3p counterparts is more likely to be regulated by different cis-eQTLs. Fine mapping and bioinformatics analyses pinpoint the peak cis-eSNP of miR-4331-5p, rs344650810, which is located in its seed region, as a causal variant for the changes in expression and function of this miRNA. Additionally, rs344650810 is significantly (P < 0.01) correlated with the density and percentage of type I muscle fibers. Altogether, this study provides a comprehensive atlas of miRNA-eQTLs in porcine skeletal muscle and new insights into regulatory mechanisms of miRNA expression.

From 'Omics to Multi-omics Technologies: the Discovery of Novel Causal Mediators

PURPOSE OF REVIEW

'Omics studies provide a comprehensive characterisation of a biological entity, such as the genome, epigenome, transcriptome, proteome, metabolome, or microbiome. This review covers the unique properties of these types of 'omics and their roles as causal mediators in cardiovascular disease. Moreover, applications and challenges of integrating multiple types of 'omics data to increase predictive power, improve causal inference, and elucidate biological mechanisms are discussed.

RECENT FINDINGS

Multi-omics approaches are growing in adoption as they provide orthogonal evidence and overcome the limitations of individual types of 'omics data. Studies with multiple types of 'omics data have improved the diagnosis and prediction of disease states and afforded a deeper understanding of underlying pathophysiological mechanisms, beyond any single type of 'omics data. For instance, disease-associated loci in the genome can be supplemented with other 'omics to prioritise causal genes and understand the function of non-coding variants. Alternatively, techniques, such as Mendelian randomisation, can leverage genetics to provide evidence supporting a causal role for disease-associated molecules, and elucidate their role in disease pathogenesis. As technologies improve, costs for 'omics studies will continue to fall and datasets will become increasingly accessible to researchers. The intrinsically unbiased nature of 'omics data is well-suited to exploratory analyses that discover causal mediators of disease, and multi-omics is an emerging discipline that leverages the strengths of each type of 'omics data to provide insights greater than the sum of its parts.