Genetic architecture modulates diet-induced hepatic mRNA and miRNA expression profiles in Diversity Outbred mice

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Finally, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility     

Obesity and genetics regulate microRNAs in islets, liver, and adipose of diabetic mice

Type 2 diabetes results from severe insulin resistance coupled with a failure of β cells to compensate by secreting sufficient insulin. Multiple genetic loci are involved in the development of diabetes, although the effect of each gene on diabetes susceptibility is thought to be small. MicroRNAs (miRNAs) are noncoding 19–22-nucleotide RNA molecules that potentially regulate the expression of thousands of genes. To understand the relationship between miRNA regulation and obesity-induced diabetes, we quantitatively profiled approximately 220 miRNAs in pancreatic islets, adipose tissue, and liver from diabetes-resistant (B6) and diabetes-susceptible (BTBR) mice. More than half of the miRNAs profiled were expressed in all three tissues, with many miRNAs in each tissue showing significant changes in response to genetic obesity. Furthermore, several miRNAs in each tissue were differentially responsive to obesity in B6 versus BTBR mice, suggesting that they may be involved in the pathogenesis of diabetes. In liver there were approximately 40 miRNAs that were downregulated in response to obesity in B6 but not BTBR mice, indicating that genetic differences between the mouse strains play a critical role in miRNA regulation. In order to elucidate the genetic architecture of hepatic miRNA expression, we measured the expression of miRNAs in genetically obese F2 mice. Approximately 10% of the miRNAs measured showed significant linkage (miR-eQTLs), identifying loci that control miRNA abundance. Understanding the influence that obesity and genetics exert on the regulation of miRNA expression will reveal the role miRNAs play in the context of obesity-induced type 2 diabetes.     

Using hippocampal microRNA expression differences between mouse inbred strains to characterise miRNA function

Micro-RNAs (miRNAs) are short, single-stranded, noncoding RNAs that are involved in the regulation of protein-coding genes at the level of messenger RNA (mRNA). They are involved in the regulation of numerous traits, including developmental timing, apoptosis, immune function, and neuronal development. To better understand how the expression of the miRNAs themselves is regulated, we looked for miRNA expression differences among four mouse inbred strains, A/J, BALB/cJ, C57BL/6J, and DBA/2J, in one tissue, the hippocampus. A total of 166 miRNA RT-PCR assays were used to screen RNA pools for each strain. Twenty miRNA species that were markedly different between strains were further investigated using eight individual samples per strain, and 11 miRNAs showed significant differences across strains (p < 0.05). This is the first observation of miRNA expression differences across inbred mice strains. We conducted an in silico correlation analysis of the expression of these differentially expressed miRNAs with phenotype data and mRNA expression to better characterise the effects of these miRNAs on both phenotype and the regulation of mRNA expression. This approach has allowed us to nominate miRNAs that have potential roles in anxiety, exploration, and learning and memory.     

cis-acting determinants of basal and lipid-regulated apolipoprotein A-IV expression in mice

The levels of apolipoprotein A-IV (apoA-IV) mRNA are regulated by dietary lipid in the liver of both the mouse and rat. Thirteen different inbred mouse strains were fed a high lipid diet, and the effect on apoA-IV liver mRNA levels was examined. It was found that each strain responded in one of two ways. Mice of four strains had higher liver apoA-IV mRNA levels as compared with syngeneic mice fed a normal chow diet. Mice of the other nine strains had decreased liver apoA-IV mRNA levels as compared with syngeneic mice fed a normal chow diet. Using F1 hybrids between mice from BALB/c, C3H, and C57BL/6 and between 129 and C57BL/6, as well as recombinant inbred strains derived from a cross between BALB/c and C57BL/6, we have shown that both the normal level of liver apoA-IV mRNA in the chow-fed mice and the lipid-dependent regulation of apoA-IV mRNA levels are controlled by cis-acting genetic elements. The apoA-IV mRNA levels in mice fed a normal diet varied dramatically among strains, with the largest difference (90-fold) being between the 129/J inbred strain and the C57BL/6J strain. In addition, we have examined the expression of apoA-IV during mouse development. ApoA-IV mRNA is expressed early in mouse liver (16 days postcoitum), whereas others have shown previously that rat liver apoA-IV mRNA is undetectable until 14 days after birth. ApoA-IV mRNA levels in the intestine and apoA-I mRNA levels in the liver and intestine, by contrast, mirror the pattern seen in the rat.     

MicroRNA Buffering and Altered Variance of Gene Expression in Response to Salmonella Infection

One potential role of miRNAs is to buffer variation in gene expression, although conflicting results have been reported. To investigate the buffering role of miRNAs in response to Salmonella infection in pigs, we sequenced miRNA and mRNA in whole blood from 15 pig samples before and after Salmonella challenge. By analyzing inter-individual variation in gene expression patterns, we found that for moderately and lowly expressed genes, putative miRNA targets showed significantly lower expression variance compared with non-miRNA-targets. Expression variance between highly expressed miRNA targets and non-miRNA-targets was not significantly different. Further, miRNA targets demonstrated significantly reduced variance after challenge whereas non-miRNA-targets did not. RNA binding proteins (RBPs) are significantly enriched among the miRNA targets with dramatically reduced variance of expression after Salmonella challenge. Moreover, we found evidence that targets of young (less-conserved) miRNAs showed lower expression variance compared with targets of old (evolutionarily conserved) miRNAs. These findings point to the importance of a buffering effect of miRNAs for relatively lowly expressed genes, and suggest that the reduced expression variation of RBPs may play an important role in response to Salmonella infection.     

Genetic diversity within honeybee colonies prevents severe infections and promotes colony growth

Multiple mating by social insect queens increases the genetic diversity among colony members, thereby reducing intracolony relatedness and lowering the potential inclusive fitness gains of altruistic workers. Increased genetic diversity may be adaptive, however, by reducing the prevalence of disease within a nest. Honeybees, whose queens have the highest levels of multiple mating among social insects, were investigated to determine whether genetic variation helps to prevent chronic infections. I instrumentally inseminated honeybee queens with semen that was either genetically similar (from one male) or genetically diverse (from multiple males), and then inoculated their colonies with spores of Ascosphaera apis, a fungal pathogen that kills developing brood. I show that genetically diverse colonies had a lower variance in disease prevalence than genetically similar colonies, which suggests that genetic diversity may benefit colonies by preventing severe infections.     

Gene expression profiles in genetically different mice infected with Toxoplasma gondii: ALDH1A2, BEX2, EGR2, CCL3 and PLAU

Toxoplasma gondii can modulate host cell gene expression; however, determining gene expression levels in intermediate hosts after T. gondii infection is not known much. We selected 5 genes (ALDH1A2, BEX2, CCL3, EGR2 and PLAU) and compared the mRNA expression levels in the spleen, liver, lung and small intestine of genetically different mice infected with T. gondii. ALDH1A2 mRNA expressions of both mouse strains were markedly increased at day 1-4 postinfection (PI) and then decreased, and its expressions in the spleen and lung were significantly higher in C57BL/6 mice than those of BALB/c mice. BEX2 and CCR3 mRNA expressions of both mouse strains were significantly increased from day 7 PI and peaked at day 15-30 PI (P<0.05), especially high in the spleen liver or small intestine of C57BL/6 mice. EGR2 and PLAU mRNA expressions of both mouse strains were significantly increased after infection, especially high in the spleen and liver. However, their expression patterns were varied depending on the tissue and mouse strain. Taken together, T. gondii-susceptible C57BL/6 mice expressed higher levels of these 5 genes than did T. gondii-resistant BALB/c mice, particularly in the spleen and liver. And ALDH1A2 and PLAU expressions were increased acutely, whereas BEX2, CCL3 and EGR2 expressions were increased lately. Thus, these demonstrate that host genetic factors exert a strong impact on the expression of these 5 genes and their expression patterns were varied depending on the gene or tissue.     

Regulation by let-7 and lin-4 miRNAs results in target mRNA degradation

MicroRNAs (miRNAs) are approximately 22 nucleotide RNAs that negatively regulate the expression of protein-coding genes. In a present model of miRNA function in animals, miRNAs that form imperfect duplexes with their targets inhibit protein expression without affecting mRNA levels. Here, we report that in C. elegans, regulation by the let-7 miRNA results in degradation of its lin-41 target mRNA, despite the fact that its 3'UTR regulatory sequences can only partially base-pair with the miRNA. Furthermore, lin-14 and lin-28 are targets of the lin-4 miRNA, and we show that the mRNA levels for these protein-coding genes significantly decrease in response to lin-4 expression. This study reveals that mRNAs containing partial miRNA complementary sites can be targeted for degradation in vivo, raising the possibility that regulation at the level of mRNA stability may be more common than previously appreciated for the miRNA pathway.     

miRNA异构体：miRNA序列多样性

微RNA（microRNA,miRNA)是一类抑制基因表达的调控分子,在多种生物学进程中扮演重要角色。近来,基于新一代测序技术获得的小RNA测序数据发现,对于某一miRNA来说,它并不是单一的序列,而是由一系列长度/序列及表达不同的异构微RNA (isomicroRNA, isomiRNA/isomiR)所组成。这些isomiR表达多样且序列多样,甚至引入多样的5′端及种子区域。特定miRNA位点在疾病组织中可具有异常的表达模式,现已证实,部分isomiR具有重要的生物学功能。所关注的经典miRNA序列,其实仅是多重isomiR序列中的1条特殊序列,仅从miRNA角度的研究已不足以揭示小RNA的奥秘,全面的研究应同时在miRNA和isomiR中开展,由此可进一步拓宽miRNA的研究思路。本文主要从isomiR的生物学特点、表达及功能等方面进行介绍。     

Genetic architecture of highly complex chemical resistance traits across four yeast strains

Many questions about the genetic basis of complex traits remain unanswered. This is in part due to the low statistical power of traditional genetic mapping studies. We used a statistically powerful approach, extreme QTL mapping (X-QTL), to identify the genetic basis of resistance to 13 chemicals in all 6 pairwise crosses of four ecologically and genetically diverse yeast strains, and we detected a total of more than 800 loci. We found that the number of loci detected in each experiment was primarily a function of the trait (explaining 46% of the variance) rather than the cross (11%), suggesting that the level of genetic complexity is a consistent property of a trait across different genetic backgrounds. Further, we observed that most loci had trait-specific effects, although a small number of loci with effects in many conditions were identified. We used the patterns of resistance and susceptibility alleles in the four parent strains to make inferences about the allele frequency spectrum of functional variants. We also observed evidence of more complex allelic series at a number of loci, as well as strain-specific signatures of selection. These results improve our understanding of complex traits in yeast and have implications for study design in other organisms.     

Strain-based genetic differences regulate the efficiency of systemic gene delivery as well as expression.

We have characterized the impact of strain-based genetic differences on the efficiency of the intravenous cationic liposome-DNA complex (CLDC)-based gene transfer and expression in mice. We also investigated what steps in the gene delivery and expression pathway appeared responsible for these strain-related differences and whether such differences could be compensated for either by agents that alter host pathways important in CLDC-mediated gene transfer and expression, or by changes in CLDC formulation. We found that different mouse strains can exhibit different expression levels and/or differences in the amount of plasmid DNA delivered to the organs where the DNA is expressed. Furthermore, drug pretreatment or reformulation of the CLDC could improve DNA delivery and/or gene expression in a strain-specific fashion. We conclude that genetic factors critically modify both the tissue deposition and the expression of genetic materials delivered by CLDC. Because manipulation of either the host or the CLDC could at least partially compensate for these strain-related differences, such strategies may be required to effectively use non-viral gene transfer approaches in genetically diverse populations.     

MicroRNA expression in abdominal and gluteal adipose tissue is associated with mRNA expression levels and partly genetically driven

To understand how miRNAs contribute to the molecular phenotype of adipose tissues and related traits, we performed global miRNA expression profiling in subcutaneous abdominal and gluteal adipose tissue of 70 human subjects and characterised which miRNAs were differentially expressed between these tissues. We found that 12% of the miRNAs were significantly differentially expressed between abdominal and gluteal adipose tissue (FDR adjusted p<0.05) in the primary study, of which 59 replicated in a follow-up study of 40 additional subjects. Further, 14 miRNAs were found to be associated with metabolic syndrome case-control status in abdominal tissue and three of these replicated (primary study: FDR adjusted p<0.05, replication: p<0.05 and directionally consistent effect). Genome-wide genotyping was performed in the 70 subjects to enable miRNA expression quantitative trait loci (eQTL) analysis. Candidate miRNA eQTLs were followed-up in the additional 40 subjects and six significant, independent cis-located miRNA eQTLs (primary study: p<0.001; replication: p<0.05 and directionally consistent effect) were identified. Finally, global mRNA expression profiling was performed in both tissues to enable association analysis between miRNA and target mRNA expression levels. We find 22% miRNAs in abdominal and 9% miRNAs in gluteal adipose tissue with expression levels significantly associated with the expression of corresponding target mRNAs (FDR adjusted p<0.05). Taken together, our results indicate a clear difference in the miRNA molecular phenotypic profile of abdominal and gluteal adipose tissue, that the expressions of some miRNAs are influenced by cis-located genetic variants and that miRNAs are associated with expression levels of their predicted mRNA targets.     

Characterization of miRNomes in Acute and Chronic Myeloid Leukemia Cell Lines

Myeloid leukemias are highly diverse diseases and have been shown to be associated with microRNA（miRNA） expression aberrations. The present study involved an in-depth miRNome analysis of two human acute myeloid leukemia（AML） cell lines, HL-60 and THP-1, and one human chronic myeloid leukemia（CML） cell line, K562, via massively parallel signature sequencing. mRNA expression profiles of these cell lines that were established previously in our lab facilitated an integrative analysis of miRNA and mRNA expression patterns. miRNA expression profiling followed by differential expression analysis and target prediction suggested numerous miRNA signatures in AML and CML cell lines. Some miRNAs may act as either tumor suppressors or oncomiRs in AML and CML by targeting key genes in AML and CML pathways. Expression patterns of cell type-specific miRNAs could partially reflect the characteristics of K562, HL-60 and THP-1 cell lines, such as actin filament-based processes, responsiveness to stimulus and phagocytic activity. miRNAs may also regulate myeloid differentiation, since they usually suppress differentiation regulators. Our study provides a resource to further investigate the employment of miRNAs in human leukemia subtyping, leukemogenesis and myeloid development. In addition, the distinctive miRNA signatures may be potential candidates for the clinical diagnosis, prognosis and treatment of myeloid leukemias.     

miRNA，siRNA，saRNA与基因表达调控

近年来的研究发现,生物体内存在着大量的非编码RNA（non．codingRNAs,ncRNA）,它们在染色质修饰、基因转录、RNA剪接和mRNA翻译等多种水平上参与了基因表达的调控。ncRNA中的小分子RNA如miRNA能够识别特定的目标mRNA,通过与mRNAs3’非翻译区结合,影响mRNA转录及蛋白质翻译;siRNA是RNA干扰的引发物,能够导致与dsRNA同源的mRNA降解,进而抑制相应基因表达;saRNA是目前最新发现的一种靶向目的基因启动子区的在转录水平激活目的基因表达的dsRNA。miRNA、siRNA和saRNA在生成机制、作用途径等方面关系密切,既区别又相互联系,小分子RNA的研究将是今后分子生物学的研究热点之一。     

Natural quantitative genetic variance in plant growth differs in response to ecologically relevant temperature heterogeneity

Adaptation to large‐scale spatial heterogeneity in the environment accounts for a major proportion of genetic diversity within species. Theory predicts the erosion of adaptive genetic variation on a within‐population level, but considerable genetic diversity is often found locally. Genetic diversity could be expected to be maintained within populations in temporally or spatially variable conditions if genotypic rank orders vary across contrasting microenvironmental settings. Taking advantage of fine‐resolution environmental data, we tested the hypothesis that temperature heterogeneity among years could be one factor maintaining quantitative genetic diversity within a natural and genetically diverse plant population. We sampled maternal families of Boechera stricta, an Arabidopsis thaliana relative, at one location in the central Rocky Mountains and grew them in three treatments that, based on records from an adjacent weather station, simulated hourly temperature changes at the native site during three summers with differing mean temperatures. Treatment had a significant effect on all traits, with 2–3‐fold increase in above‐ and belowground biomass and the highest allocation to roots observed in the treatment simulating the warmest summer on record at the site. Treatment affected bivariate associations between traits, with the weakest correlation between above‐ and belowground biomass in the warmest treatment. The magnitude of quantitative genetic variation for all traits differed across treatments: Genetic variance of biomass was 0 in the warmest treatment, while highly significant diversity was found in average conditions, resulting in broad‐sense heritability of 0.31. Significant genotype × environment interactions across all treatments were found only in root‐to‐shoot ratio. Therefore, temperature variation among summers appears unlikely to account for the observed levels of local genetic variation in size in this perennial species, but may influence family rank order in growth allocation. Our results indicate that natural environmental fluctuations can have a large impact on the magnitude of within‐population quantitative genetic variance.