Genetic regulation and heritability of miRNA and mRNA expression link to phosphorus utilization and gut microbiome

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.

Genetic regulation and variation of expression of miRNA and mRNA transcripts in fetal muscle tissue in the context of sex, dam and variable fetal weight

BACKGROUND

Impaired skeletal muscle growth in utero can result in reduced birth weight and pathogenesis of intrauterine growth restriction. Fetal and placental growth is influenced by many factors including genetic, epigenetic and environmental factors. In fact, the sex and genotype of the fetus itself, as well as the mother providing it with a suitable environment, influence the growth of the fetus. Hence, our goal was to decipher and elucidate the molecular pathways of developmental processes mediated by miRNAs and mRNAs in fetal muscle tissue in the context of sex, dam, and fetal weight. Therefore, we analyse the variation of miRNA and mRNA expression in relation to these factors. In addition, the coincidence of genetic regulation of these mRNAs and miRNAs, as revealed by expression quantitative trait loci (eQTL) analyses, with sex-, mother- and weight-associated expression was investigated.

METHODS

A three-generation pig F2 population (n = 118) based on reciprocal crossing of German Landrace (DL) and Pietrain (Pi) was used. Genotype information and transcriptomic data (mRNA and miRNA) from longissimus dorsi muscle (LDM) of pig fetuses sampled at 63 days post-conception (dpc) were used for eQTL analyses.

RESULTS

The transcript abundances of 13, 853, and 275 probe-sets were influenced by sex, dam and fetal weight at 63 dpc, respectively (FDR < 5%). Most of significant transcripts affected by sex were located on the sex chromosomes including KDM6A and ANOS1 or autosomes including ANKS1B, LOC100155138 and miR-153. The fetal muscle transcripts associated with fetal weight indicated clearer metabolic directions than maternally influenced fetal muscle transcripts. Moreover, coincidence of genetic regulation (eQTL) and variation in transcript abundance due to sex, dam and fetal weight were identified.

CONCLUSIONS

Integrating information on eQTL, sex-, dam- and weight-associated differential expression and QTL for fetal weight allowed us to identify molecular pathways and shed light on the basic biological processes associated with differential muscle development in males and females, with implications for adaptive fetal programming.

Genetic Polymorphisms in miR-604A>G, miR-938G>A, miR-1302-3C>T and the Risk of Idiopathic Recurrent Pregnancy Loss

The purpose of this study was to investigate whether polymorphisms in five microRNAs (miRNAs), miR-604A>G, miR-608C>G, 631I/D, miR-938G>A, and miR-1302-3C>T, are associated with the risk of idiopathic recurrent pregnancy loss (RPL). Blood samples were collected from 388 patients with idiopathic RPL (at least two consecutive spontaneous abortions) and 227 control participants. We found the miR-604 AG and AG + GG genotypes of miR-604, the miR-938 GA and GA + AA genotypes of miR-938, and the miR-1302-3CT and CT + TT genotypes of miR-1302-3 are less frequent than the wild-type (WT) genotypes, miR-604AA, miR-938GG, and miR-1302-3CC, respectively, in RPL patients. Using allele-combination multifactor dimensionality reduction (MDR) analysis, we found that eight haplotypes conferred by the miR-604/miR-608/miR-631/miR-938/miR-1302-3 allele combination, A-C-I-G-T, A-C-I-A-C, G-C-I-G-C, G-C-I-G-T, G-G-I-G-C, G-G-I-G-T, G-G-I-A-C, G-G-D-G-C, three from the miR-604/miR-631/miR-938/miR-1302-3 allele combination, A-I-G-T, G-I-G-C, G-I-A-T, one from the miR-604/miR-631/miR-1302-3 allele combination, G-I-C, and two from the miR-604/miR-1302-3 allele combination, G-C and G-T, were less frequent in RPL patients, suggesting protective effects (all p < 0.05). We also identified the miR-604A>G and miR-938G>A polymorphisms within the seed sequence of the mature miRNAs and aligned the seed sequences with the 3′UTR of putative target genes, methylenetetrahydrofolate reductase (MTHFR) and gonadotropin-releasing hormone receptor (GnRHR), respectively. We further found that the binding affinities between miR-604/miR-938 and the 3′UTR of their respective target genes (MTHFR, GnRHR) were significantly different for the common (miR-604A, miR-938G) and variant alleles (miR-604G, miR-938A). These results reveal a significant association between the miR-604A>G and miR-938G>A polymorphisms and idiopathic RPL and suggest that miRNAs can affect RPL in Korean women.