Differential microRNA response to a high-cholesterol, high-fat diet in livers of low and high LDL-C baboons
© Karere et al.; 2012
Received: 21 October 2011
Accepted: 30 June 2012
Published: 18 July 2012
Dysregulation of microRNA (miRNA) expression has been implicated in molecular genetic events leading to the progression and development of atherosclerosis. We hypothesized that miRNA expression profiles differ between baboons with low and high serum low-density lipoprotein cholesterol (LDL-C) concentrations in response to diet, and that a subset of these miRNAs regulate genes relevant to dyslipidemia and risk of atherosclerosis.
Using Next Generation Illumina sequencing methods, we sequenced hepatic small RNA libraries from baboons differing in their LDL-C response to a high-cholesterol, high-fat (HCHF) challenge diet (low LDL-C, n = 3; high LDL-C, n = 3), resulting in 517 baboon miRNAs: 490 were identical to human miRNAs and 27 were novel. We compared miRNA expression profiles from liver biopsies collected before and after the challenge diet and observed that HCHF diet elicited expression of more miRNAs compared to baseline (chow) diet for both low and high LDL-C baboons. Eighteen miRNAs exhibited differential expression in response to HCHF diet in high LDL-C baboons compared to 10 miRNAs in low LDL-C baboons. We used TargetScan/Base tools to predict putative miRNA targets; miRNAs expressed in high LDL-C baboons had significantly more gene targets than miRNAs expressed in low LDL-C responders. Further, we identified miRNA isomers and other non-coding RNAs that were differentially expressed in response to the challenge diet in both high LDL-C and low LDL-C baboons.
We sequenced and annotated baboon liver miRNAs from low LDL-C and high LDL-C responders using high coverage Next Gen sequencing methods, determined expression changes in response to a HCHF diet challenge, and predicted target genes regulated by the differentially expressed miRNAs. The identified miRNAs will enrich the database for non-coding small RNAs including the extent of variation in these sequences. Further, we identified other small non-coding RNAs differentially expressed in response to diet. Our discovery of differentially expressed baboon miRNAs in response to a HCHF diet challenge that differ by LDL-C phenotype is a fundamental step in understating the role of non-coding RNAs in dyslipidemia.
Cardiovascular disease (CVD), a leading cause of mortality in the United States, is commonly due to the development of atherosclerosis. Atherosclerosis is a chronic inflammatory disease affecting medium and large arteries, resulting from a complex interaction between genotype and environmental factors including diet, which influences dyslipidemia. Clinical and epidemiological studies indicate that atherosclerosis is positively correlated with serum LDL-C and inversely correlated with high-density lipoprotein cholesterol (HDL-C)[2–5], demonstrating that increased levels of serum LDL-C is a major risk factor for atherosclerosis. However, the exact mechanism by which LDL-C promotes atherosclerosis remains an active aspect for research. Early studies by Brown and Goldstein demonstrated that oxidized LDL-C (ox-LDL-C) is a prerequisite for initiation of atherosclerotic lesions, leading to the development of intermediate and complex plaques[7, 8]. Subsequent studies have confirmed that ox-LDL-C is an essential atherosclerotic risk factor, triggering arterial cellular modification[9–11]. Lowering LDL-C reduces the risk for atherosclerosis in humans, and modification of diet and lifestyle, and administration of drugs is a major strategy to control CVD.
Apart from the physiological characteristics, understanding the genetic mechanisms underlying dyslipidemia and subsequent development of atherosclerosis is fundamental to therapeutic interventions. Many genetic factors influencing risk factors for atherosclerosis have been identified, including genes underlying the vulnerability of arterial intimal layer to ox-LDL-C[10, 14]. The expression of genes important for development of atherosclerosis is modulated by a variety of elements, including microRNAs (miRNAs).
miRNAs are small endogenous non-coding RNAs that are evolutionarily conserved across many species and the biogenesis of miRNAs has been described[15, 16]. miRNAs down-regulate gene expression by translational repression, degradation and deadenylation[17, 18]. Occasionally miRNAs can also up-regulate gene expression[19, 20]. These small non-coding RNAs have been implicated in a range of biological processes, including cell growth, differentiation, apoptosis, and cholesterol metabolism[15, 21]. These processes are important for the initiation and progression of atherosclerosis. Consequently, dysregulation of miRNA function may lead to atherosclerosis. For example, miR-365 is up-regulated by ox-LDL-C, leading to reduced expression of BCL2, an antiapoptotic gene, in arterial endothelial cells. miR-155 and miR-146a are associated with accumulation of ox-LDL-C in monocytes[23, 24]. It has been shown that miR-92a controls angiogenesis of atherosclerotic plaques in mice, leading to destabilization and rupture. miRNAs, let-7 and miR-17-92 cluster, are known to target thrombospondin-1, an inhibitor of angiogenesis. Furthermore, miR-210, -15b, -16, and -20b are implicated in down-regulation of vascular endothelial growth factor[27–29], a key gene in cell migration in the atherosclerotic lesion. Studies also indicate that the up-regulation of miR-335 and −122 is associated with lipid metabolism in obese mice[30, 31] and chimpanzee. This suggests that expressed miRNAs play a fundamental role in dyslipidemia, a major risk factor for atherosclerosis. Importantly miRNAs have been shown to exhibit temporal and spatial expression patterns, emphasizing the need to define cell–type–specific miRNA expression profiles, which require expression-sensitive tools to capture low-abundance small RNAs.
Although the baboon is a well-characterized model for human biomedical studies, including CVD, very little is known about baboon (Papio hamadryas) miRNAs. Discovering liver miRNAs, identifying diet responsive miRNAs that differ between low and high LDL-C responders, and understanding how miRNAs regulate gene expression in the baboon will provide insights on the role of miRNAs in dyslipidemia. In this study, we have used Next-Generation sequencing methods to perform deep sequencing of small RNA libraries derived from liver biopsies collected before and after a HCHF challenge diet from six half-sib baboons (low LDL-C, n = 3; high LDL-C, n = 3), differing in their LDL-C response to dietary fat and cholesterol. The sequence reads were aligned to human genome; notably the current draft of baboon genome is an unannotated draft assembly which does not provide chromosomal coordinates or transcript identification as with the human genome sequence. However, based on physical maps, the baboon and human genomes exhibit considerable sequence conservation and synteny. We identified differentially expressed miRNAs in response to the challenge diet in the two LDL-C phenotypes and used in-silico methods to identify miRNA targets.
Annotation of small RNAs and identification of novel miRNA genes
Baboon liver small RNA libraries: number of clusters, sequence tags and their expression levels
No. of Clusters
Mapped unique tags
% tags mapped to human miRNAs
% miRNA expression
miRNAs expressed in baboon livers
Identical baboon-human miRNAs
Novel baboon miRNAs
Common to low and high LDL-C
Frequency of miRNA expression for each LDL-C phenotype and each diet (n = 3 per group)
Identical baboon-human miRNAs (mean 397.50 ± 19.82)
Novel baboon miRNAs (mean 4.08 ± 1.98)
miRNA expression categorized by diet for each LDL-C phenotype
Identical baboon-human miRNAs
Differential expression profiling of miRNAs responsive to dietary fat
There was a statistical difference (Fisher’s exact test; p < 0.05) between the number of miRNAs differentially expressed in high LDL-C baboons in response to HCHF diet (n = 18) compared to 10 miRNAs for low LDL-C responders. Two miRNA families are shared between the two baboon phenotypes. In this study miR-29 family members were significantly down-regulated in response to HCHF diet in both low and high LDL-C baboons. Interestingly, the expression of some miRNA gene families showed polycistronic regulation, such that a subset of members of one family was either down-regulated in one phenotype or up-regulated in another. For example, a subset of miRNA let-7 family members was down-regulated in low LDL-C baboons but up-regulated in high LDL-C baboons (Additional file1). In the Additional file3: Table S1 lists the miRNAs that were differentially expressed, p-value, fold change and the complete mature miRNA sequences.
Moreover, some miRNAs were strictly diet-specific; these miRNAs, although not differentially expressed, exhibited peculiar characteristics in which they were expressed either in response to chow diet but not HCHF diet, or vice versa (Additional file3: Table S2). In addition some miRNAs were polymorphic between low and high LDL-C baboons. For example, miR-181c and 302a were expressed in low LDL-C baboon livers while 181a and 320b were detected in high LDL-C baboon livers.
We observed that a substantial portion of sequence tags perfectly mapped to other non-coding RNAs including snoRNA, siRNA, snRNAs, and tRNAs in Rfam database. These RNAs accounted for 20-33% of the total expressed tags in the two phenotypic groups (Figures 2b and 2C). Moreover, 31% (1,064/3,456) of other non-coding RNAs in high LDL-C baboons were differentially expressed between chow and HCHF diets (False Discovery Rate (FDR) ≤ 0.05) (Additional file4). These small RNAs were not differentially expressed in low LDL-C baboon livers in response to the HCHF diet.
Targets of differentially expressed miRNAs in baboon livers
Number of predicted and validated gene targets for the 28 differentially expressed miRNA families
miRNAs families differentially expressed in response to HCHF diet, p ≤ 0.05
Human gene targets
Rhesus gene targets
Overlap between rhesus and human genes
Experimentally validated targets
Novel baboon miRNAs categorized by diet and LDL-C. The numbers of predicted target genes, and the 2-7mer seed region for each novel miRNA, are shown
No. of targets
No. of targets
Total unique targets
Total unique targets
No. of targets
No. of targets
Total unique targets
Total unique targets
No. of targets
No. of targets
Total unique targets
Total unique targets
Baboon and dyslipidemia
High LDL-C is the major risk factor for atherosclerosis, the leading cause of CVD. Underlying the pathogenesis of atherosclerosis are complex networks of genes regulated by a plethora of elements including miRNAs. In this study we investigated differential expression of hepatic miRNAs in baboons differing in serum LDL-C concentrations in response to HCHF diet, and predicted miRNA targets. Liver is a major organ involved in lipid metabolism and regulation of serum cholesterol including LDL-C, and hepatic inflammation is associated with early atherosclerosis. To understand the liver microRNAome relevant to LDL-C phenotype and dietary fat and cholesterol, we performed high coverage Next-Generation sequencing, annotated the baboon liver small RNA transcriptome and quantified expression of small RNAs. The baboon is a well-characterized model for dyslipidemia and atherosclerosis[35, 36], and is very similar to humans both physiologically and genetically. In addition to genetic background, we are able to control diet and environment in a nonhuman primate, which is not feasible in humans. For this study, animals exhibiting discordant LDL-C serum measures (low and high LDL-C responders) were challenged for 7 weeks with a HCHF diet. Liver samples were collected from six half-sibs discordant for LDL-C before and after the challenge diet. By analyzing half-sib baboons discordant for the quantitative trait (< 2 S.D. difference), we minimized genetic variation due to genetic background and increased the likelihood that we would discover miRNAs that are diet responsive and that differ by LDL-C phenotype.
Identification of baboon liver miRNAs
The distribution of non-coding RNA sequences identified in this study is in concordance with previous studies identifying miRNAs using deep sequencing methods. The length distribution of sequence reads indicates that 22-nt fragments were the most abundant (Figure2), corresponding to the average size of mature miRNAs. We identified various categories of small non-coding RNAs from baboon liver samples and observed that miRNAs were the most abundant (62%), (Figure2a), in concordance with previous studies[38, 39].
Various methods were applied in the identification of baboon miRNAs. The first study to identify baboon miRNAs employed a combined approach involving genomic computation and microarray expression. In our present study, we identified 490 miRNAs identical to human miRNAs. This number far exceeds the 227 liver miRNAs identified by genomic computational prediction. Moreover, the total number of expressed baboon liver miRNAs identified by sequencing increased by 21% that detected by microarray. In addition we identified 27 novel baboon miRNAs with corresponding complimentary (‘star’) sequences (Table 2). The presence of star miRNAs is compelling evidence for the DICER-like processing from a miRNA hairpin and stringency measures applied in the analysis. The identified baboon miRNAs will contribute to the species-specific miRNA database. These results demonstrate the power of parallel and high-throughput Next-Generation sequencing in elucidating differential expression profiles of miRNAs in baboon livers.
Next Generation (NG) deep sequencing has become a method of choice for discovering miRNA sequences and expression. NG technology allows detection of low copy number miRNAs in the transcriptome and discovery of novel small RNAs, with excellent reproducibility. Although QRT-PCR is the method of choice for validating array-based RNA expressions, it is a relative measure of expression normalized against expression of a ‘house-keeping’ gene. In contrast, expression profiling by deep sequencing is normalized against the whole transcriptome providing more robust data in a high coverage sequencing experiment than what can be attained by QRT-PCR. For example, in this study the average read counts for baboon liver small RNAs was approximately 800,000. Moreover, NG sequencing enables detection of variation in mature miRNAs and RNA editing. Thus, deep sequencing remains an extremely sensitive approach for identifying miRNA expression levels and for detecting novel and variant miRNA sequences.
Expression of small non-coding RNAs
Expression profiles of miRNAs associated with metabolic disorder risk factors
Low LDL-C (Chow vs HCHF)
High LDL-C (Chow vs HCHF)
Biological processes associated with disease risk factors
Demethylation of MMP-2/MMP-9 genes involved in cell migration in aorta
Regulator of vascular smooth muscle cell function
Angiogenesis in hypertrophy
Angiogenesis in hypertrophy
Proliferation of vascular smooth muscle cells
Cell proliferation in carotid artery
Repolarization and diastolic activity
LDL-C oxidization accumulation in monocytes
Oxidized LDL-C accumulation in monocytes
Regulation of pro-inflammation in aorta
Although miR-33a/b were not differentially expressed in response to the HCHF challenge diet both miRNAs were detected in low and high LDL-C livers, however, at differing amplitudes. miR-33a was more abundant in low LDL-C while miR-33b was more abundant in high LDL-C livers. These findings are consistent with previous studies that indicate miR-33a is highly expressed in low sterol condition whereas miR-33b is more abundant in cholesterol-filled cells. miR-33 family is transcribed from the introns of sterol-regulatory element-binding factor (SREBF) isoforms and is linked to cholesterol homeostasis.
In this study, we identified differentially expressed miRNAs that are predominantly engulfed by HDL particles. miR-877 and 191 were differentially expressed in low and high LDL-C livers, respectively. The results are consistent with findings by Vickers and colleagues, indicating that miR-877 is abundantly expressed in normal individuals while miR-191 is up-regulated in familial hypercholesterolemia (FH) human subjects. Mutation in LDL receptor (LDLR) leads to FH and severe development of atherosclerotic lesions in homozygous individuals. In addition we found that miR-19 was differentially expressed in high LDL-C baboons and not in low LDL-C responders. miR-19 is abundant in LDLR mutant mice. Together, these findings indicate that these miRNAs are potentially involved in regulation of dyslipidemia and HDL particles may be involved in their transport to recipient cells as an anti-atherogenic mechanism.
During the biogenesis of miRNAs, a duplex RNA is generated. Only a single strand (mature miRNA) is loaded into RNA-induced silencing complex (RISC) while a near complementary strand, miRNA star (miRNA*), is degraded[54, 55]. Recent studies based on deep sequencing have uncovered the presence of miRNA* in different systems but in low abundance. This indicates the power of deep sequencing in detecting low copy number miRNAs previously thought to be degraded. In our study we confirmed presence of miRNAs* (Additional file1). In addition we discovered that some miRNAs* exhibit higher expression than the mature miRNAs, are inversely regulated between low and high LDL-C responders, or are detected in some samples where mature miRNAs are absent. For example, miR-222* was up-regulated in response to HCHF diet in low LDL-C baboon livers but down-regulated in high LDL-C responders while mature sequences were down-regulated in both LDL-C phenotypes. These observations indicate that processing of miRNAs is regulated and that regulation of processing can affect gene expressions.
Although miRNA families have a consensus 5’ end ‘seed’ region, little is reported on upstream polymorphisms, which may affect the binding efficiency of the miRNA on mRNA UTRs, in relation to the phenotype of interest. In this study we identified some miRNAs that have nucleotide polymorphisms between low and high LDL-C responders including miR-181 and 302. miR-181 variants have identical ‘seed’ region; however, variant 181a has a G allele insertion flanking the 8-mer seed region in comparison to 181c. miR-181 has a consensus binding site in humans, chimpanzee and rhesus macaque but not mice and rats. In comparison, 302a and c have G/A substitutions at the 3’ end, contributing to differing total context scores and binding efficiency. Mice, rats, rhesus macaque, chimpanzee and humans have a consensus binding sequence for miR-302. In contrast rabbits, which are often used to model dyslipidemia and atherosclerosis do not have sites for miR-181 or 302.
Furthermore, differentially expressed miRNAs in baboon livers (Figure3) have been implicated to have multiple roles in different types of diseases. For example miR-16 is reportedly over-expressed in lung cancer, down-regulated in prostate cancer and up-regulated in rheumatoid arthritis. Down-regulation of miR-34b is associated with tumorogenesis suppression and Parkinson’s disease and its suppression via methylation confers anti-apoptotic effects on colorectal cancer. miR-34b targets include TENC1 (tensin like C1 domain containing phosphatase), which has been implicated in apoptotic events in vitro. These observations suggest that miRNAs have different roles in distinct tissues and that the effects on gene expression should be interpreted in the context of cell state.
Our results indicate that a significant proportion of expressed small non-coding RNAs include other miRNA-sized small non-coding RNAs such as snoRNAs, siRNAs, snRNAs and tRNAs. Non-coding RNAs (n = 3,455) were expressed in high LDL-C responders while 4,645 were identified in low LDL-C livers. We observed that 31% (1,064/3,455) of other non-coding RNAs in high LDL-C baboons were differentially expressed (FDR ≤ 0.05) in response to HCHF diet. Non-coding RNAs were not differentially expressed in low LDL-C baboons. A total of 2,277 RNAs were expressed in both baboon phenotypes including 99% of those differentially expressed in high LDL-C livers. The functional role of these small RNAs in dyslipidemia is not clear. However, prior reports indicate that deep sequencing generates a considerable number of sequence tags perfectly mapping to other known small non-coding RNAs. At least one study identified snoRNAs with miRNA-like function and DICER processing signatures. It is plausible these RNAs may have functional roles in dyslipidemia.
miRNAs are often expressed in gene clusters, and their transcription is thought to be closely linked[16, 66]. In this study we observed that miRNA clusters and cluster members were coordinately regulated. For example, the majority of let-7 family members are down-regulated in low LDL-C responders but up-regulated in high LDL-C responders. However, not all members of the family exhibited a coordinated pattern, suggesting an independent regulation event. A recent study confirmed that some miRNAs show differential expression in distinct cells due to post-transcriptional regulation associated with variation in the levels of DICER, an enzyme involved in biosynthesis of mature miRNA. It is plausible that the expressions of some let-7 cluster members are dependent on cell context or independent of cluster transcription unit. Recent studies have uncovered the association of clustered miRNAs and complex diseases including B-cell leukemia[67, 68] and testicular carcinoma[69–71]. These miRNA clusters tend to occur in fragile regions of genomes and may be species-specific[16, 72, 73].
miRNA gene targets and atherosclerosis
Many of the predicted targets of expressed miRNA family members, with identical seed regions, overlapped considerably. Thus, it is logical to focus on identical seed regions or the miRNA family names when considering targets of miRNAs; miRNAs with same seed region are classified as one family. In this study, we sought to identify genes targeted by expressed miRNA families in low and high LDL-C baboons using in silico methods.
For differentially expressed miRNAs responsive to the HCHF diet, we observed a substantial (> 90%) overlap between targets identified in human and rhesus macaque genomes (Table 5). Because the genome assembly of baboon is in draft form, we used the rhesus genome to identify predicted gene targets. Baboon is closely related to rhesus macaque with 98% sequence identity. This conservation of a large percentage of miRNA target sites between human and rhesus and the high sequence identity between rhesus and baboon suggests that our findings in baboons may translate to humans given that our data constitute sequences that are identical between human and baboon genomes. Nonetheless, a recent study has shown that evolutionary changes in miRNA binding sites, may inhibit or enhance the effect of miRNAs on gene expression regulation. The study by Richardson et al. revealed that the PLIN4 3’ UTR has undergone nucleotide substitution allowing the binding of miR-522 for human individuals with A but not G allele. The miR-522 binding site with A allele was not found in Neanderthal and non-human primates suggesting a recent evolutionally change. Until an annotated baboon genome sequence is available we will not be able to assess conservation of miRNA targets sites between baboon and human and baboon and rhesus.
In this study, we observed many-genes-to-many-miRNAs regulation modules, suggesting that many genes are coordinately and cooperatively regulated by miRNAs. Genes targeted by multiple miRNAs are tightly regulated and may show graded expression in response to expression of different miRNAs. In addition, genes tightly regulated by redundant miRNAs may be an indication of their central role in the cell. We observed that, on average, about 2,000 genes are targeted by a single miRNA family, which presents daunting challenges for the development of a pharmaceutical agent with specific target. Moreover miRNA target prediction revealed that some genes are cooperatively targeted by miRNAs differing in their response to HCHF diet; i.e., one miRNA was up-regulated while the other was down-regulated in response to the HCHF diet, but both miRNAs have the same target. For example, LDLRAP1 (low density lipoprotein receptor adapter protein 1, which interacts with LDL-C receptor) is targeted by miR-27a and -26a that are up- and down-regulated, respectively, in response to the HCHF diet. The molecular interactions of these miRNAs with LDLRAP1 are not known, but it is plausible this phenomenon may be an inherent mechanism to balance gene expression levels. miR-27a and 26a are differentially regulated and implicated in adipogenesis[76, 77].
We estimate that differentially expressed miRNAs in our study targeted 9% of the total genes implicated in atherosclerosis[78, 79]. The following miRNA targets are included among those potentially associated with atherosclerosis: 1) LDLR, a cell surface protein involved in receptor-mediated endocytosis of LDL-C. Mutations in this gene cause the autosomal dominant disorder familial hypercholesterolemia and increased carotid artery intima-media thickness. 2) VLDLR (very low density lipoprotein receptor) plays important roles in VLDL triglyceride metabolism. 3) LIPA (lipase A precursor) catalyzes the hydrolysis of cholesterol esters and triglycerides in lysosomes, and mutations in the gene is associated with impaired endothelial function in CVD. Other targets of interests include PTEN (phosphatase and tensin homolog), which functions as a tumor suppressor by negatively regulating intracellular levels of phosphatidylinositol-3,4,5-trisphosphate, substrates for AKT (tyrosine kinase B protein) in the intracellular AKT/PKB signaling pathway, and TENC1 (tensin-like C1 domain containing phosphatase), which is thought to act in a similar manner as PTEN in depriving substrates for AKT. Further, we observed that ACVR1B (activin A receptor, type IB isoform c precursor) is considerably targeted by both novel and known miRNAs. ACVR1B is well-conserved in primates including human, rhesus, and orangutan, but relatively less-conserved in chimpanzee. ACVR1B is a receptor for activin growth factor, which belongs to TGF-beta (transforming growth factor-beta) superfamily. TGF-beta is involved and implicated in inflammatory processes, such as atherosclerosis in arteries. Altogether this study identifies potential atherogenic candidate genes, targeted by differentially expressed baboon miRNAs. Identifying the roles of miRNAs and their target genes and signaling pathways in atherosclerosis will be critical for future research. miRNAs represent new regulatory factors for cardiovascular biology and are an emerging novel class of therapeutic targets for CVD.
Although identification of miRNA gene targets is a key step in elucidating miRNA functional mechanisms, comprehensive identification of targets still presents a daunting challenge. There are many bioinformatics tools available to predict miRNA targets; however, consistency among different platforms is limited. In this study, we employed TargetScan, which is the most widely used tool in target prediction. Because of the limitation of the high-throughput in-silico prediction, experimental validations of the miRNA-mRNA interaction are necessary. The experiments include reporter gene assay and loss- or gain-of-function of miRNA to measure target repression. However, these validation experiments are cumbersome, and, as a result, there are only a few validated miRNA targets. For example, our results generated using IPA indicate only 1% of the total miRNA targets have been validated. Recently, a high-throughput and sensitive technique for identification and quantification of miRNA targets has been reported. The technique involves in-situ co-immunoprecipitation of mRNA associated with RISC, loaded with a specific miRNA, and subsequent deep sequencing. The technique, RISC-Seq, is a modification of a previous approach that utilized microarray profiling and has been further optimized. There is no doubt the technique provides important insight into determining the full spectrum of miRNA-mRNA interactions. However, the method is not yet fully optimized, and the modifications suggested by Matkovich et al. present challenges due to over-expression of Argonaute2 protein in RISC complex, which may induce out-of-context cell state. Second, the technique does not address all outcomes of functions of miRNAs; a substantial number of miRNAs act by silencing the translation machinery and epigenetic modifications. Thus, until the RISC-Seq approach is fully optimized, in-silico prediction and basic experimental validations will remain useful approaches in deciphering the miRNA-mRNA interactions.
We have sequenced baboon liver small RNA samples from half-sibs discordant for serum LDL-C and were fed chow and HCHF diets. We identified 517 baboon miRNAs including 490 that are identical to human and 27 novel miRNAs. The list of miRNAs reported in this study will contribute to the public database for small RNAs. Further, we identified 30 differentially expressed miRNAs in response to a HCHF challenge diet in two distinct baboon phenotypes previously characterized for extreme responses to dietary fat. Six of these differentially expressed miRNAs have been previously identified as regulators of genes that play roles in CVD. We also identified miRNA variants and other non-coding RNAs differentially expressed in response to the HCHF diet. Analysis of target genes of the differentially expressed miRNAs showed that 9% of the currently annotated CVD-related genes are predicted to be regulated by these miRNAs. Future studies that include gene sequencing in these samples to identify all expressed genes and gene variants in combination with gene/miRNA interactions in these LDL-C responders will take us closer to understanding the genetic and epigenetic networks that underlie LDL-C response to dietary fat and cholesterol, thus providing potential central genetic targets for miRNA-based therapeutic agents to regulate LDL-C.
Baboons (Papio hamadryas) were maintained by the Southwest National Primate Research Center (SNPRC) at the Texas Biomedical Research Institute (Texas Biomed), which is accredited by the Association for Assessment and Accreditation of Laboratory Animal Care International. Experimental protocols were approved by the Institutional Animal Care and Use Committee and were conducted under the supervision of Southwest National Primate Research Center staff veterinarians.
Selection of baboon sib-pairs discordant for LDL-C
Based on phenotypic and genotypic analysis of 951 pedigreed baboons, we identified three half-sib-pairs with contrasting phenotypes for LDL-C as described. The sib-pairs differed by at least two standard deviations for LDL-C serum concentrations. In addition, members of each selected sibling pair were discordant for at least 1 marker loci within the support interval for an LDL-C QTL on baboon chromosome 4.
High-cholesterol, high-fat (HCHF) dietary challenge
The mean (mg/dL) and standard deviation of the LDL-C levels in low and high LDL-C baboons on chow and HCHF diets
High LDL-C Responders
Low LDL-C Responders
Liver biopsy collection
Baboons were sedated with ketamine (10 mg/kg), given atropine (0.025 mg/kg) and intubated. Anesthesia was induced and maintained with isoflurane (1-2%). Blood pressure was measured by automated arm cuff (Collin) and oxygen saturation, heart rate, and respiration were monitored by pulse oximetry. Biopsies were collected from the left lobe of the liver using Temno Evolution Needle, gauge 14 (CareFusion, San Diego, CA). During post biopsy recovery analgesia was provided in the form of Stadol, 0.15 mg/kg, bid, for 3 days and ampicillin, 25 mg/day for 10 days.
Small RNA library preparation
miRNA-enriched fractions were isolated from 10 mg of baboon liver (n = 12) using the mirVanaTM miRNA Isolation Kit (Ambion) according to the manufacturer’s protocol and stored at −80°C until further use. Small RNA libraries were prepared from the enriched fractions with Illumina Small RNA Prep Kit v1.5 following the manufacturer’s instructions. Briefly, the Illumina 3’ adapter (5’-ATCTCGTATGCCGTCTTCTGCTTGT) was ligated to 2 μg of enriched small RNAs using T4 RNA truncated ligase (New England BioLabs) and 15% PEG at 22°C for 2 hrs. Then, the Illumina 5’ adapter (5’-GTTCAGAGTTCTACAGTCCGACGATC) was ligated onto the 3’ adapter-linked RNA using T4 RNA ligase (New England BioLabs) for 1 hr at 20°C. The 5’ and 3’ adapter-linked RNA was converted to a single-stranded cDNA using SuperScript II Reverse Transcriptase (Invitrogen) and Illumina’s RT-Primer (5’-CAAGCAGAAGACGGCATACGA) following the manufacturer’s protocol. The cDNA was PCR-amplified with Hotstart Phusion DNA Polymerase using the following profile: denaturation at 98°C for 30s, 15 cycles of 98°C for 10s, 60°C for 30s, 72°C for 15 s, and a final extension at 72°C for 10 min.
The PCR product was purified on a 6% TBE urea polyacrylamide gel (Invitrogen). The gel was stained with ethidium bromide and visualized under UV light. The gel band corresponding to 93–100 bp was excised, passed through a bottom-perforated 0.5 ml tube and eluted in 1X Elution Buffer following the manufacturer’s protocol. Notably, while the average length of a miRNA is 22 nucleotides (nts), the excised band corresponds to the total length of the adapter-ligated miRNAs. The resulting gel slurry was spun in a Spin-X filter and RNA precipitated with 2ul of glycogen, 10ul of 3 M-ammonium acetate and 325 μl of 100% cold ethanol. The mixture was centrifuged, and the resulting pellet washed with 70% ethanol. The dried RNA pellet was suspended in 10 μl of Resuspension Buffer (Illumina).
Next generation sequencing of small RNAs
Small RNA sequencing was performed using Illumina’s Genome Analyzer (GAIIx). Diluted cDNA (10 nmol/L) was amplified, sequenced, and single-read clusters generated using cBot Cluster Generation Kit and TruSeq SBS Kit v5 (36 cycles). For the control sample, we used PhiX Control Kit v3 (Illumina).
Analysis of sequence reads
Reads filtering and adapter trimming
miRNA sequence reads were analyzed using mirTools, an interactive web-based tool for analyzing deep sequencing reads. Briefly, the small RNA sequence reads were filtered to exclude low-quality reads and 3’/5’ adapter- and polyA-sequences. The reads were trimmed off the 3’ adapter sequences using a custom perl script. The resulting clean full-length reads were formatted into a non-redundant FASTA file. Each sequence tag in the file contains the number of reads represented by the tag.
Small RNA annotation
The resulting sequence tags (18–30 nts) were mapped onto the human genome sequence (hg18) using the SOAP2 program embedded in mirTools. Subsequently, the reads were aligned against various databases including miRBase, Rfam, RepeatMasker, and the reference genome transcriptome database. Consequently the sequence tags were classified as a known miRNA, a degraded fragment from non-coding RNA, or RNA derived from repeat or coding sequences. Sequences that did not conform to any of the categories but mapped to a unique location on the human genome were termed “unclassified.”
Discovery of novel miRNAs
The prediction of novel miRNAs was performed by miRDeep software and associated RNAfold software, incorporated in mirTools using default settings. Unclassified tags were used to predict candidate novel miRNAs. Using RNAfold software, a sequence fragment containing 100 nts flanking each side of the locus aligning to an unclassified tag was excised from the genome sequence. RNAfold software folds and performs computational assessment of the fragment to ensure that there is a minimum of 14 base pairings between the mature and star sequences of the hairpin. Because false positive hairpins can be formed from different types of small RNAs in the genome, miRDeep software performs further stringent analysis to test whether reads aligning to a potential miRNA precursor (pre-miRNA) are consistent with the Dicer processing signature to simulate miRNA biogenesis in the cell; presence of 2-nt 3’ overhangs and positional alignment of the reads on mature, star, or loop of the precursor sequence. Extracted sequences that do not have a Dicer signature and fold into a hairpin were discarded. An unclassified sequence tag is considered novel miRNA if it aligns to a candidate pre-miRNA. Sequences with dominantly abundant reads were denoted as mature miRNA and its complementary sequence as miRNA*.
Detection of differential expression of small RNAs
miRNAs and other non-coding RNAs differentially expressed between chow and HCHF samples were determined by assessing the read counts as described. A sequence tag consists of one or more reads, which reflects the number of times the tag was sequenced, hence the expression level. We normalized an expression level of a sequence tag, within and between samples, by dividing the absolute read count by the total read count of the sample and multiplied by 1,000,000. The normalized count corresponds to the relative expression level of each tag. Normalized expression data generated from animals on chow and challenge diets were used to determine differential expression levels in each phenotype. A two-tailed student t-test was performed to test for significant expression levels for miRNAs expressed on both diets and Welch’s t-test, assuming unequal variance, was performed for miRNAs with expression detected in only one of the two groups. A small non-coding RNA was considered differentially expressed between samples when the p-value was ≤ 0.05, and multiple testing adjustments applied, where appropriate, using False Discovery Rate.
We identified potential targets of diet responsive differentially expressed miRNAs in low LDL-C and high LDL-C responders using TargetScan/Base Release 5.1 (http://www.targetscan.org/) and Ingenuity Pathway Analysis (http://www.ingenuity.com) tools. We identified novel miRNA targets using TargetScan, employing the 7-mer seed-region of each miRNA sequence to customize prediction of targets.
High-density lipoprotein cholesterol
Low-density lipoprotein cholesterol
B-cell lymphoma protein 2 beta isoform
Short oligonucleotide alignment program
Small interfering RNAs
- sno RNAs:
Small nucleolar RNAs
Small nuclear RNAs
RNA-induced silencing complex
- TENC1 :
Tensin like C1 domain containing phosphatase
- VLDLR :
Very low density lipoprotein receptor
- LDLR :
Low density lipoprotein receptor
- LIPA :
Lipase A precursor
- PTEN :
Phosphatase and tensin homolog
- AKT :
Tyrosine kinase B protein
- TENC1 :
Tensin-like C1 domain containing phosphatase
- ACVR1B :
Activin A receptor, type IB isoform c precursor
- TGF-beta :
Transforming growth factor-beta
False Discovery Rate.
Acknowledgements and funding
We thank Dr. Kimberly D. Spradling for assistance in data analysis and manuscript proofreading, Ms. Jane VandeBerg for assistance in measurements of serum LDL-C concentrations in baboons, and Mr. Roy Garcia for sequencing the small RNA libraries. This work was supported by National Institutes of Health grants P01 HL028972-27, P01 HL028972-27 S1 and P51 RR013986. This investigation was conducted in part in facilities constructed with support from Research Facilities Improvement Program Grant Number C06 RR013556 and C06 RR015456 from the National Center for Research Resources, National Institutes of Health. Genesio Karere is a recipient of a research grant in 2009 from the Southwest Foundation Forum Research Support.
- Ross R: The pathogenesis of atherosclerosis: a perspective for the 1990s. Nature. 1993, 362 (6423): 801-809.PubMedGoogle Scholar
- McGill HC, McMahan CA, Kruski AW, Mott GE: Relationship of lipoprotein cholesterol concentrations to experimental atherosclerosis in baboons. Arteriosclerosis. 1981, 1 (1): 3-12.PubMedGoogle Scholar
- Amarenco P, Goldstein LB, Szarek M, Sillesen H, Rudolph AE, Callahan A, Hennerici M, Simunovic L, Zivin JA, Welch KM, SPARCL Investigators: Effects of intense low-density lipoprotein cholesterol reduction in patients with stroke or transient ischemic attack: the Stroke Prevention by Aggressive Reduction in Cholesterol Levels (SPARCL) trial. Stroke. 2007, 38 (12): 3198-3204.PubMedGoogle Scholar
- Amarenco P, Labreuche J, Touboul PJ: High-density lipoprotein-cholesterol and risk of stroke and carotid atherosclerosis: a systematic review. Atherosclerosis. 2008, 196 (2): 489-496.PubMedGoogle Scholar
- Assmann G, Schulte H: Relation of high-density lipoprotein cholesterol and triglycerides to incidence of atherosclerotic coronary artery disease (the PROCAM experience). Prospective Cardiovascular Munster study. Am J Cardiol. 1992, 70 (7): 733-737.PubMedGoogle Scholar
- Brown MS, Goldstein JL: Lipoprotein metabolism in the macrophage: implications for cholesterol deposition in atherosclerosis. Annu Rev Biochem. 1983, 52: 223-261.PubMedGoogle Scholar
- McGill HC, McMahan CA, Herderick EE, Tracy RE, Malcom GT, Zieske AW, Strong JP: Effects of coronary heart disease risk factors on atherosclerosis of selected regions of the aorta and right coronary artery. PDAY Research Group. Pathobiological Determinants of Atherosclerosis in Youth. Arterioscler Thromb Vasc Biol. 2000, 20 (3): 836-845.PubMedGoogle Scholar
- Witztum JL, Steinberg D: Role of oxidized low density lipoprotein in atherogenesis. J Clin Invest. 1991, 88 (6): 1785-1792.PubMed CentralPubMedGoogle Scholar
- Jialal I, Devaraj S: The role of oxidized low density lipoprotein in atherogenesis. J Nutr. 1996, 126 (4 Suppl): 1053S-1057S.PubMedGoogle Scholar
- Daugherty A, Dunn JL, Rateri DL, Heinecke JW: Myeloperoxidase, a catalyst for lipoprotein oxidation, is expressed in human atherosclerotic lesions. J Clin Invest. 1994, 94 (1): 437-444.PubMed CentralPubMedGoogle Scholar
- Schiopu A, Frendeus B, Jansson B, Soderberg I, Ljungcrantz I, Araya Z, Shah PK, Carlsson R, Nilsson J, Fredrikson GN: Recombinant antibodies to an oxidized low-density lipoprotein epitope induce rapid regression of atherosclerosis in apobec-1(−/−)/low-density lipoprotein receptor(−/−) mice. J Am Coll Cardiol. 2007, 50 (24): 2313-2318.PubMedGoogle Scholar
- Kushwaha RS, McGill HC: Diet, plasma lipoproteins and experimental atherosclerosis in baboons (Papio sp.). Hum Reprod Update. 1998, 4 (4): 420-429.PubMedGoogle Scholar
- Monraats PS, Fang Y, Pons D, Pires NM, Pols HA, Zwinderman AH, de Maat MP, Doevendans PA, DeWinter RJ, Tio RA, Waltenberger J, Frants RR, Quax PH, van der Laarse A, van der Wall EE, Uitterlinden AG, Jukema JW: Vitamin D receptor: a new risk marker for clinical restenosis after percutaneous coronary intervention. Expert Opin Ther Targets. 2010, 14 (3): 243-251.PubMedGoogle Scholar
- Zeibig S, Li Z, Wagner S, Holthoff HP, Ungerer M, Bultmann A, Uhland K, Vogelmann J, Simmet T, Gawaz M, Munch G: Effect of the oxLDL binding protein Fc-CD68 on plaque extension and vulnerability in atherosclerosis. Circ Res. 2011, 108 (6): 695-703.PubMedGoogle Scholar
- Ambros V: The functions of animal microRNAs. Nature. 2004, 431 (7006): 350-355.PubMedGoogle Scholar
- Karere GM, Glenn JP, VandeBerg JL, Cox LA: Identification of baboon microRNAs expressed in liver and lymphocytes. J Biomed Sci. 2010, 17: 54-PubMed CentralPubMedGoogle Scholar
- Galasso M, Elena Sana M, Volinia S: Non-coding RNAs: a key to future personalized molecular therapy?. Genome Med. 2010, 2 (2): 12-PubMed CentralPubMedGoogle Scholar
- Wu L, Fan J, Belasco JG: MicroRNAs direct rapid deadenylation of mRNA. Proc Natl Acad Sci U S A. 2006, 103 (11): 4034-4039.PubMed CentralPubMedGoogle Scholar
- Orom UA, Nielsen FC, Lund AH: MicroRNA-10a binds the 5'UTR of ribosomal protein mRNAs and enhances their translation. Mol Cell. 2008, 30 (4): 460-471.PubMedGoogle Scholar
- Ritchie W, Rajasekhar M, Flamant S, Rasko JE: Conserved expression patterns predict microRNA targets. PLoS Comput Biol. 2009, 5 (9): e1000513-PubMed CentralPubMedGoogle Scholar
- Guo S, Lu J, Schlanger R, Zhang H, Wang JY, Fox MC, Purton LE, Fleming HH, Cobb B, Merkenschlager M, Golub TR, Scadden DT: MicroRNA miR-125a controls hematopoietic stem cell number. Proc Natl Acad Sci U S A. 2010, 107 (32): 14229-14234.PubMed CentralPubMedGoogle Scholar
- Qin B, Xiao B, Liang D, Xia J, Li Y, Yang H: MicroRNAs expression in ox-LDL treated HUVECs: MiR-365 modulates apoptosis and Bcl-2 expression. Biochem Biophys Res Commun. 2011, 410 (1): 127-133.PubMedGoogle Scholar
- Yang K, He YS, Wang XQ, Lu L, Chen QJ, Liu J, Sun Z, Shen WF: MiR-146a inhibits oxidized low-density lipoprotein-induced lipid accumulation and inflammatory response via targeting toll-like receptor 4. FEBS Lett. 2011, 585 (6): 854-860.PubMedGoogle Scholar
- Huang RS, Hu GQ, Lin B, Lin ZY, Sun CC: MicroRNA-155 silencing enhances inflammatory response and lipid uptake in oxidized low-density lipoprotein-stimulated human THP-1 macrophages. J Investig Med. 2010, 58 (8): 961-967.PubMedGoogle Scholar
- Bonauer A, Carmona G, Iwasaki M, Mione M, Koyanagi M, Fischer A, Burchfield J, Fox H, Doebele C, Ohtani K, Chavakis E, Potente M, Tjwa M, Urbich C, Zeiher AM, Dimmeler S: MicroRNA-92a controls angiogenesis and functional recovery of ischemic tissues in mice. Science. 2009, 324 (5935): 1710-1713.PubMedGoogle Scholar
- Kuehbacher A, Urbich C, Dimmeler S: Targeting microRNA expression to regulate angiogenesis. Trends Pharmacol Sci. 2008, 29 (1): 12-15.PubMedGoogle Scholar
- Cascio S, D'Andrea A, Ferla R, Surmacz E, Gulotta E, Amodeo V, Bazan V, Gebbia N, Russo A: miR-20b modulates VEGF expression by targeting HIF-1 alpha and STAT3 in MCF-7 breast cancer cells. J Cell Physiol. 2010, 224 (1): 242-249.PubMedGoogle Scholar
- Fasanaro P, D'Alessandra Y, Di Stefano V, Melchionna R, Romani S, Pompilio G, Capogrossi MC, Martelli F: MicroRNA-210 modulates endothelial cell response to hypoxia and inhibits the receptor tyrosine kinase ligand Ephrin-A3. J Biol Chem. 2008, 283 (23): 15878-15883.PubMed CentralPubMedGoogle Scholar
- Liu XD, Wu X, Yin YL, Liu YQ, Geng MM, Yang HS, Blachier F, Wu GY: Effects of dietary L: -arginine or N-carbamylglutamate supplementation during late gestation of sows on the miR-15b/16, miR-221/222. 2011, Amino Acids: VEGFA and eNOS expression in umbilical veinGoogle Scholar
- Esau C, Davis S, Murray SF, Yu XX, Pandey SK, Pear M, Watts L, Booten SL, Graham M, McKay R, Subramaniam A, Propp S, Lollo BA, Freier S, Bennett CF, Bhanot S, Monia BP: miR-122 regulation of lipid metabolism revealed by in vivo antisense targeting. Cell Metab. 2006, 3 (2): 87-98.PubMedGoogle Scholar
- Nakanishi N, Nakagawa Y, Tokushige N, Aoki N, Matsuzaka T, Ishii K, Yahagi N, Kobayashi K, Yatoh S, Takahashi A, Suzuki H, Urayama O, Yamada N, Shimano H: The up-regulation of microRNA-335 is associated with lipid metabolism in liver and white adipose tissue of genetically obese mice. Biochem Biophys Res Commun. 2009, 385 (4): 492-496.PubMedGoogle Scholar
- Lanford RE, Hildebrandt-Eriksen ES, Petri A, Persson R, Lindow M, Munk ME, Kauppinen S, Orum H: Therapeutic silencing of microRNA-122 in primates with chronic hepatitis C virus infection. Science. 2010, 327 (5962): 198-201.PubMed CentralPubMedGoogle Scholar
- Cox LA, Mahaney MC, Vandeberg JL, Rogers J: A second-generation genetic linkage map of the baboon (Papio hamadryas) genome. Genomics. 2006, 88 (3): 274-281.PubMedGoogle Scholar
- Kleemann R, Verschuren L, van Erk MJ, Nikolsky Y, Cnubben NH, Verheij ER, Smilde AK, Hendriks HF, Zadelaar S, Smith GJ, Kaznacheev V, Nikolskaya T, Melnikov A, Hurt-Camejo E, van der Greef J, van Ommen B, Kooistra T: Atherosclerosis and liver inflammation induced by increased dietary cholesterol intake: a combined transcriptomics and metabolomics analysis. Genome Biol. 2007, 8 (9): R200-PubMed CentralPubMedGoogle Scholar
- McGill HC, McMahan CA, Kruski AW, Kelley JL, Mott GE: Responses of serum lipoproteins to dietary cholesterol and type of fat in the baboon. Arteriosclerosis. 1981, 1 (5): 337-344.PubMedGoogle Scholar
- McGill HC, McMahan CA, Zieske AW, Tracy RE, Malcom GT, Herderick EE, Strong JP: Association of Coronary Heart Disease Risk Factors with microscopic qualities of coronary atherosclerosis in youth. Circulation. 2000, 102 (4): 374-379.PubMedGoogle Scholar
- Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75 (5): 843-854.PubMedGoogle Scholar
- Vaz C, Ahmad HM, Sharma P, Gupta R, Kumar L, Kulshreshtha R, Bhattacharya A: Analysis of microRNA transcriptome by deep sequencing of small RNA libraries of peripheral blood. BMC Genomics. 2010, 11: 288-PubMed CentralPubMedGoogle Scholar
- Morin RD, O'Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu AL, Zhao Y, McDonald H, Zeng T, Hirst M, Eaves CJ, Marra MA: Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res. 2008, 18 (4): 610-621.PubMed CentralPubMedGoogle Scholar
- t Hoen PA, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, de Menezes RX, Boer JM, van Ommen GJ, den Dunnen JT: Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36 (21): e141-Google Scholar
- Felli N, Fontana L, Pelosi E, Botta R, Bonci D, Facchiano F, Liuzzi F, Lulli V, Morsilli O, Santoro S, Valtieri M, Calin GA, Liu CG, Sorrentino A, Croce CM, Peschle C: MicroRNAs 221 and 222 inhibit normal erythropoiesis and erythroleukemic cell growth via kit receptor down-modulation. Proc Natl Acad Sci U S A. 2005, 102 (50): 18081-18086.PubMed CentralPubMedGoogle Scholar
- Maxwell KN, Soccio RE, Duncan EM, Sehayek E, Breslow JL: Novel putative SREBP and LXR target genes identified by microarray analysis in liver of cholesterol-fed mice. J Lipid Res. 2003, 44 (11): 2109-2119.PubMedGoogle Scholar
- Cheng C, Tempel D, Den Dekker WK, Haasdijk R, Chrifi I, Bos FL, Wagtmans K, van de Kamp EH, Blonden L, Biessen EA, Moll F, Pasterkamp G, Serruys PW, Schulte-Merker S, Duckers HJ: Ets2 determines the inflammatory state of endothelial cells in advanced atherosclerotic lesions. Circ Res. 2011, 109 (4): 382-395.PubMedGoogle Scholar
- Leeper NJ, Raiesdana A, Kojima Y, Chun HJ, Azuma J, Maegdefessel L, Kundu RK, Quertermous T, Tsao PS, Spin JM: MicroRNA-26a is a novel regulator of vascular smooth muscle cell function. J Cell Physiol. 2011, 226 (4): 1035-1043.PubMed CentralPubMedGoogle Scholar
- Quintavalle M, Elia L, Condorelli G, Courtneidge SA: MicroRNA control of podosome formation in vascular smooth muscle cells in vivo and in vitro. J Cell Biol. 2010, 189 (1): 13-22.PubMed CentralPubMedGoogle Scholar
- Ji R, Cheng Y, Yue J, Yang J, Liu X, Chen H, Dean DB, Zhang C: MicroRNA expression signature and antisense-mediated depletion reveal an essential role of MicroRNA in vascular neointimal lesion formation. Circ Res. 2007, 100 (11): 1579-1588.PubMedGoogle Scholar
- Yang L, Duan R, Chen D, Wang J, Chen D, Jin P: Fragile X mental retardation protein modulates the fate of germline stem cells in Drosophila. Hum Mol Genet. 2007, 16 (15): 1814-1820.PubMedGoogle Scholar
- Li T, Cao H, Zhuang J, Wan J, Guan M, Yu B, Li X, Zhang W: Identification of miR-130a, miR-27b and miR-210 as serum biomarkers for atherosclerosis obliterans. Clin Chim Acta. 2011, 412 (1–2): 66-70.PubMedGoogle Scholar
- Fang Y, Shi C, Manduchi E, Civelek M, Davies PF: MicroRNA-10a regulation of proinflammatory phenotype in athero-susceptible endothelium in vivo and in vitro. Proc Natl Acad Sci U S A. 2010, 107 (30): 13450-13455.PubMed CentralPubMedGoogle Scholar
- Rayner KJ, Esau CC, Hussain FN, McDaniel AL, Marshall SM, van Gils JM, Ray TD, Sheedy FJ, Goedeke L, Liu X, Khatsenko OG, Kaimal V, Lees CJ, Fernandez-Hernando C, Fisher EA, Temel RE, Moore KJ: Inhibition of miR-33a/b in non-human primates raises plasma HDL and lowers VLDL triglycerides. Nature. 2011, 478 (7369): 404-407.PubMed CentralPubMedGoogle Scholar
- Fernandez-Hernando C, Moore KJ: MicroRNA modulation of cholesterol homeostasis. Arterioscler Thromb Vasc Biol. 2011, 31 (11): 2378-2382.PubMed CentralPubMedGoogle Scholar
- Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT: MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. 2011, 13 (4): 423-433.PubMed CentralPubMedGoogle Scholar
- Rader DJ, Cohen J, Hobbs HH: Monogenic hypercholesterolemia: new insights in pathogenesis and treatment. J Clin Invest. 2003, 111 (12): 1795-1803.PubMed CentralPubMedGoogle Scholar
- Friedlander MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, Rajewsky N: Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotechnol. 2008, 26 (4): 407-415.PubMedGoogle Scholar
- Filipowicz W, Bhattacharyya SN, Sonenberg N: Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?. Nat Rev Genet. 2008, 9 (2): 102-114.PubMedGoogle Scholar
- Zhang C, Zheng H, Yu Q, Yang P, Li Y, Cheng F, Fan J, Liu E: A practical method for quantifying atherosclerotic lesions in rabbits. J Comp Pathol. 2010, 142 (2–3): 122-128.PubMedGoogle Scholar
- Navarro A, Diaz T, Gallardo E, Vinolas N, Marrades RM, Gel B, Campayo M, Quera A, Bandres E, Garcia-Foncillas J, Ramirez J, Monzo M: Prognostic implications of miR-16 expression levels in resected non-small-cell lung cancer. J Surg Oncol. 2011, 103 (5): 411-415.PubMedGoogle Scholar
- Porkka KP, Ogg EL, Saramaki OR, Vessella RL, Pukkila H, Lahdesmaki H, van Weerden WM, Wolf M, Kallioniemi OP, Jenster G, Visakorpi T: The miR-15a-miR-16-1 locus is homozygously deleted in a subset of prostate cancers. Genes Chromosomes Cancer. 2011, 50 (7): 499-509.PubMedGoogle Scholar
- Feng ZT, Li J, Ren J, Lv Z: Expression of miR-146a and miR-16 in peripheral blood mononuclear cells of patients with rheumatoid arthritis and their correlation to the disease activity. Nan Fang Yi Ke Da Xue Xue Bao. 2011, 31 (2): 320-323.PubMedGoogle Scholar
- Kalimutho M, Di Cecilia S, Del Vecchio Blanco G, Roviello F, Sileri P, Cretella M, Formosa A, Corso G, Marrelli D, Pallone F, Federici G, Bernardini S: Epigenetically silenced miR-34b/c as a novel faecal-based screening marker for colorectal cancer. Br J Cancer. 2011, 104 (11): 1770-1778.PubMed CentralPubMedGoogle Scholar
- Minones-Moyano E, Porta S, Escaramis G, Rabionet R, Iraola S, Kagerbauer B, Espinosa-Parrilla Y, Ferrer I, Estivill X, Marti E: MicroRNA profiling of Parkinson's disease brains identifies early downregulation of miR-34b/c which modulate mitochondrial function. Hum Mol Genet. 2011, 20 (15): 3067-3078.PubMedGoogle Scholar
- Toyota M, Suzuki H, Sasaki Y, Maruyama R, Imai K, Shinomura Y, Tokino T: Epigenetic silencing of microRNA-34b/c and B-cell translocation gene 4 is associated with CpG island methylation in colorectal cancer. Cancer Res. 2008, 68 (11): 4123-4132.PubMedGoogle Scholar
- Hafizi S, Ibraimi F, Dahlback B: C1-TEN is a negative regulator of the Akt/PKB signal transduction pathway and inhibits cell survival, proliferation, and migration. FASEB J. 2005, 19 (8): 971-973.PubMedGoogle Scholar
- Jung CH, Hansen MA, Makunin IV, Korbie DJ, Mattick JS: Identification of novel non-coding RNAs using profiles of short sequence reads from next generation sequencing data. BMC Genomics. 2010, 11: 77-PubMed CentralPubMedGoogle Scholar
- Ender C, Krek A, Friedlander MR, Beitzinger M, Weinmann L, Chen W, Pfeffer S, Rajewsky N, Meister G: A human snoRNA with microRNA-like functions. Mol Cell. 2008, 32 (4): 519-528.PubMedGoogle Scholar
- Sempere LF, Freemantle S, Pitha-Rowe I, Moss E, Dmitrovsky E, Ambros V: Expression profiling of mammalian microRNAs uncovers a subset of brain-expressed microRNAs with possible roles in murine and human neuronal differentiation. Genome Biol. 2004, 5 (3): R13-PubMed CentralPubMedGoogle Scholar
- Agueli C, Cammarata G, Salemi D, Dagnino L, Nicoletti R, La Rosa M, Messana F, Marfia A, Bica MG, Coniglio ML, Pagano M, Fabbiano F, Santoro A: 14q32/miRNA clusters loss of heterozygosity in acute lymphoblastic leukemia is associated with up-regulation of BCL11a. Am J Hematol. 2010, 85 (8): 575-578.PubMedGoogle Scholar
- Patz M, Pallasch CP, Wendtner CM: Critical role of microRNAs in chronic lymphocytic leukemia: overexpression of the oncogene PLAG1 by deregulated miRNAs. Leuk Lymphoma. 2010, 51 (8): 1379-1381.PubMedGoogle Scholar
- Voorhoeve PM, le Sage C, Schrier M, Gillis AJ, Stoop H, Nagel R, Liu YP, van Duijse J, Drost J, Griekspoor A, Zlotorynski E, Yabuta N, De Vita G, Nojima H, Looijenga LH, Agami R: A genetic screen implicates miRNA-372 and miRNA-373 as oncogenes in testicular germ cell tumors. Cell. 2006, 124 (6): 1169-1181.PubMedGoogle Scholar
- Murray MJ, Saini HK, van Dongen S, Palmer RD, Muralidhar B, Pett MR, Piipari M, Thornton CM, Nicholson JC, Enright AJ, Coleman N: The two most common histological subtypes of malignant germ cell tumour are distinguished by global microRNA profiles, associated with differential transcription factor expression. Mol Cancer. 2010, 9: 290-PubMed CentralPubMedGoogle Scholar
- Palmer RD, Murray MJ, Saini HK, van Dongen S, Abreu-Goodger C, Muralidhar B, Pett MR, Thornton CM, Nicholson JC, Enright AJ, Coleman N: Children's Cancer and Leukaemia Group: Malignant germ cell tumors display common microRNA profiles resulting in global changes in expression of messenger RNA targets. Cancer Res. 2010, 70 (7): 2911-2923.PubMed CentralPubMedGoogle Scholar
- Altuvia Y, Landgraf P, Lithwick G, Elefant N, Pfeffer S, Aravin A, Brownstein MJ, Tuschl T, Margalit H: Clustering and conservation patterns of human microRNAs. Nucleic Acids Res. 2005, 33 (8): 2697-2706.PubMed CentralPubMedGoogle Scholar
- Bentwich I: Prediction and validation of microRNAs and their targets. FEBS Lett. 2005, 579 (26): 5904-5910.PubMedGoogle Scholar
- Freemerman AJ, Wright RM, Flickinger CJ, Herr JC: Cloning and sequencing of baboon and cynomolgus monkey intra-acrosomal protein SP-10: homology with human SP-10 and a mouse sperm antigen (MSA-63). Mol Reprod Dev. 1993, 34 (2): 140-148.PubMedGoogle Scholar
- Richardson K, Louie-Gao Q, Arnett DK, Parnell LD, Lai CQ, Davalos A, Fox CS, Demissie S, Cupples LA, Fernandez-Hernando C, Ordovas JM: The PLIN4 variant rs8887 modulates obesity related phenotypes in humans through creation of a novel miR-522 seed site. PLoS One. 2011, 6 (4): e17944-PubMed CentralPubMedGoogle Scholar
- Kajimoto K, Naraba H, Iwai N: MicroRNA and 3 T3-L1 pre-adipocyte differentiation. RNA. 2006, 12 (9): 1626-1632.PubMed CentralPubMedGoogle Scholar
- Lin Q, Gao Z, Alarcon RM, Ye J, Yun Z: A role of miR-27 in the regulation of adipogenesis. FEBS J. 2009, 276 (8): 2348-2358.PubMedGoogle Scholar
- Rainwater DL, Cox LA, Rogers J, VandeBerg JL, Mahaney MC: Localization of multiple pleiotropic genes for lipoprotein metabolism in baboons. J Lipid Res. 2009, 50 (7): 1420-1428.PubMed CentralPubMedGoogle Scholar
- Skogsberg J, Lundstrom J, Kovacs A, Nilsson R, Noori P, Maleki S, Kohler M, Hamsten A, Tegner J, Bjorkegren J: Transcriptional profiling uncovers a network of cholesterol-responsive atherosclerosis target genes. PLoS Genet. 2008, 4 (3): e1000036-PubMed CentralPubMedGoogle Scholar
- Yamakawa K, Yanagi H, Saku K, Sasaki J, Okafuji T, Shimakura Y, Kawai K, Tsuchiya S, Takada K, Naito S: Family studies of the LDL receptor gene of relatively severe hereditary hypercholesterolemia associated with Achilles tendon xanthomas. Hum Genet. 1991, 86 (5): 445-449.PubMedGoogle Scholar
- Pauciullo P, Giannino A, De Michele M, Gentile M, Liguori R, Argiriou A, Carlotto A, Faccenda F, Mancini M, Bond MG, De Simone V, Rubba P: Increased carotid artery intima-media thickness is associated with a novel mutation of low-density lipoprotein receptor independently of major cardiovascular risk factors. Metabolism. 2003, 52 (11): 1433-1438.PubMedGoogle Scholar
- Yuan G, Liu Y, Sun T, Xu Y, Zhang J, Yang Y, Zhang M, Cianflone K, Wang DW: The therapeutic role of very low-density lipoprotein receptor gene in hyperlipidemia in type 2 diabetic rats. Hum Gene Ther. 2011, 22 (3): 302-312.PubMedGoogle Scholar
- Wild PS, Zeller T, Schillert A, Szymczak S, Sinning CR, Deiseroth A, Schnabel RB, Lubos E, Keller T, Eleftheriadis MS, Bickel C, Rupprecht HJ, Wilde S, Rossmann H, Diemert P, Cupples LA, Perret C, Erdmann J, Stark K, Kleber ME, Epstein SE, Voight BF, Kuulaasma K, Li M, Schafer AS, Klopp N, Braund PS, Sager HB, Demissie S, Proust C, Konig IR, Wichmann HE, Reinhard W, Hoffmann MM, Virtamo J, Burnett MS, Siscovick D, Wiklund PG, Qu L, El Mokthari NE, Thompson JR, Peters A, Smith AV, Yon E, Baumert J, Hengstenberg C, Marz W, Amouyel P, Devaney J, Schwartz SM, Saarela O, Mehta NN, Rubin D, Silander K, Hall AS, Ferrieres J, Harris TB, Melander O, Kee F, Hakonarson H, Schrezenmeir J, Gudnason V, Elosua R, Arveiler D, Evans A, Rader DJ, Illig T, Schreiber S, Bis JC, Altshuler D, Kavousi M, Witteman JC, Uitterlinden AG, Hofman A, Folsom AR, Barbalic M, Boerwinkle E, Kathiresan S, Reilly MP, O'Donnell CJ, Samani NJ, Schunkert H, Cambien F, Lackner KJ, Tiret L, Salomaa V, Munzel T, Ziegler A, Blankenberg S: A Genome-wide Association Study Identifies LIPA as a Susceptibility Gene for Coronary Artery Disease. Circ Cardiovasc Genet. 2011, 4 (4): 403-U203.PubMed CentralPubMedGoogle Scholar
- Hafizi S, Gustafsson A, Oslakovic C, Idevall-Hagren O, Tengholm A, Sperandio O, Villoutreix BO, Dahlback B: Tensin2 reduces intracellular phosphatidylinositol 3,4,5-trisphosphate levels at the plasma membrane. Biochem Biophys Res Commun. 2010, 399 (3): 396-401.PubMedGoogle Scholar
- Toma I, McCaffrey TA: Transforming growth factor-beta and atherosclerosis: interwoven atherogenic and atheroprotective aspects. Cell Tissue Res. 2011, 1007/s00441-011-1189-3. PMID: 21626289Google Scholar
- Ding J, Wang DZ: "RISCing" the heart: In vivo identification of cardiac microRNA targets by RISCome. Circ Res. 2011, 108 (1): 3-5.PubMedGoogle Scholar
- Xiao F, Zuo Z, Cai G, Kang S, Gao X, Li T: miRecords: an integrated resource for microRNA-target interactions. Nucleic Acids Res. 2009, 37 (Database issue): D105-10.PubMed CentralPubMedGoogle Scholar
- Hanina SA, Mifsud W, Down TA, Hayashi K, O'Carroll D, Lao K, Miska EA, Surani MA: Genome-wide identification of targets and function of individual MicroRNAs in mouse embryonic stem cells. PLoS Genet. 2010, 6 (10): e1001163-PubMed CentralPubMedGoogle Scholar
- Karginov FV, Conaco C, Xuan Z, Schmidt BH, Parker JS, Mandel G, Hannon GJ: A biochemical approach to identifying microRNA targets. Proc Natl Acad Sci U S A. 2007, 104 (49): 19291-19296.PubMed CentralPubMedGoogle Scholar
- Matkovich SJ, Van Booven DJ, Eschenbacher WH, Dorn GW: RISC RNA sequencing for context-specific identification of in vivo microRNA targets. Circ Res. 2011, 108 (1): 18-26.PubMed CentralPubMedGoogle Scholar
- Cox LA, Birnbaum S, Mahaney MC, Rainwater DL, Williams JT, VandeBerg JL: Identification of promoter variants in baboon endothelial lipase that regulate high-density lipoprotein cholesterol levels. Circulation. 2007, 116 (10): 1185-1195.PubMedGoogle Scholar
- Kammerer CM, Rainwater DL, Cox LA, Schneider JL, Mahaney MC, Rogers J, VandeBerg JL: Locus controlling LDL cholesterol response to dietary cholesterol is on baboon homologue of human chromosome 6. Arterioscler Thromb Vasc Biol. 2002, 22 (10): 1720-1725.PubMedGoogle Scholar
- Rainwater DL, Kammerer CM, Hixson JE, Carey KD, Rice KS, Dyke B, VandeBerg JF, Slifer SH, Atwood LD, McGill HC, Vandeberg JL: Two major loci control variation in beta-lipoprotein cholesterol and response to dietary fat and cholesterol in baboons. Arterioscler Thromb Vasc Biol. 1998, 18 (7): 1061-1068.PubMedGoogle Scholar
- Singh AT, Rainwater DL, Kammerer CM, Sharp RM, Poushesh M, Shelledy WR, VandeBerg JL: Dietary and genetic effects on LDL size measures in baboons. Arterioscler Thromb Vasc Biol. 1996, 16 (12): 1448-1453.PubMedGoogle Scholar
- Zhu E, Zhao F, Xu G, Hou H, Zhou L, Li X, Sun Z, Wu J: mirTools: microRNA profiling and discovery based on high-throughput sequencing. Nucleic Acids Res. 2010, 38 (Web Server issue): W392-7.PubMed CentralPubMedGoogle Scholar
- Li R, Li Y, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment program. Bioinformatics. 2008, 24 (5): 713-714.PubMedGoogle Scholar
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36 (Database issue): D154-8.PubMed CentralPubMedGoogle Scholar
- Griffiths-Jones S: Annotating non-coding RNAs with Rfam. Curr Protoc Bioinformatics. 2005, 0.1002/0471250953.bi1205s9. PMID: 18428745Google Scholar
- Tarailo-Graovac M, Chen N: Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009, 10.1002/0471250953.bi0410s25. PMID: 19274634Google Scholar
- Hofacker IL, Stadler PF: Memory efficient folding algorithms for circular RNA secondary structures. Bioinformatics. 2006, 22 (10): 1172-1176.PubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100 (16): 9440-9445.PubMed CentralPubMedGoogle Scholar
- Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007, 27 (1): 91-105.PubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.