Insights into the regulation of human CNV-miRNAs from the view of their target genes
© Wu et al.; licensee BioMed Central Ltd. 2012
Received: 22 April 2012
Accepted: 7 December 2012
Published: 18 December 2012
microRNAs (miRNAs) represent a class of small (typically 22 nucleotides in length) non-coding RNAs that can degrade their target mRNAs or block their translation. Recent research showed that copy number alterations of miRNAs and their target genes are highly prevalent in cancers; however, the evolutionary and biological functions of naturally existing copy number variable miRNAs (CNV-miRNAs) among individuals have not been studied extensively throughout the genome.
In this study, we comprehensively analyzed the properties of genes regulated by CNV-miRNAs, and found that CNV-miRNAs tend to target a higher average number of genes and prefer to synergistically regulate the same genes; further, the targets of CNV-miRNAs tend to have higher variability of expression within and between populations. Finally, we found the targets of CNV-miRNAs are more likely to be differentially expressed among tissues and developmental stages, and participate in a wide range of cellular responses.
Our analyses of CNV-miRNAs provide new insights into the impact of copy number variations on miRNA-mediated post-transcriptional networks. The deeper interpretation of patterns of gene expression variation and the functional characterization of CNV-miRNAs will help to broaden the current understanding of the molecular basis of human phenotypic diversity.
KeywordsCopy number variation miRNA Expression variation HapMap ethnic population
miRNAs are a class of small non-coding RNAs, which act through binding in a sequence-specific manner to the 3′UTR of target genes. Each miRNA can potentially regulate many transcripts and at least one-third of human genes are estimated to be miRNA targets. miRNAs participate in posttranscriptional gene regulation by repressing the expression of their target genes through inhibition of translation or cleavage of mRNAs[2–6]. miRNAs also contribute to genetic buffering of the gene expression variation, and play an important role in maintaining the identity of mature tissues through a feed-forward loop regulatory architecture[7, 8], such as the relationship between miR-9a and E(spl) in Drosophila[9, 10] and the regulation of E2F1 by miR-17 in human.
A primary goal in medical and evolutionary genomics is to understand the genetic mechanisms of natural variation in gene expression[12–16]. The structure of the human genome is highly variable and the copy number variations (CNVs) refer to alterations of genomic segments of more than 1,000 nucleotides that are present at significant frequencies within a population[17–19]. Many studies showed that CNVs can expand dosage variation of the associated genes, leading to the under-representation of dosage-sensitive protein-coding functional units such as transcription factors and members of protein complexes[20, 21]. CNVs can be discovered by cytogenetic techniques, such as fluorescent in situ hybridization, comparative genomic hybridization, array comparative genomic hybridization, and next-generation sequencing[22–24]. In humans, more than 30,000 genomic regions with segmental duplications have been recognized by systematic comparative genomic hybridizations on the DNA of healthy human subjects; however, the CNVs of other animals were far less studied (see http://projects.tcag.ca/variation). For example, only about 2,000 CNVs have been identified in Pan troglodytes and about 4,000 CNVs in inbred Mus musculus[26, 27].
Recent studies revealed a high frequency in copy number abnormality of miRNA processing genes, such as Dicer1 and Argonature2, in breast and ovarian cancers[28, 29]. Although copy number alterations of miRNAs and their regulatory genes were frequently investigated in oncogenesis[28–30], the evolutionary and functional impact of CNV-miRNAs on the human genome has not been studied extensively. Based on the human genomic structure variations, Marcinkowska et al. recently detected about 30% miRNAs located in the human CNV-regions, indicating that non-coding RNAs also have potential functional variants.
In this study, we comprehensively analyzed the properties of genes regulated by CNV-miRNAs and explored the potential involvement of CNV-miRNAs in the expression variability of their targets within and between populations. Our analysis revealed significant functional differences between the targets of CNV-miRNAs and the targets of non-CNV-miRNAs. The involvement of CNV-miRNAs in a wide range of cellular responses provided us with valuable information of the impact of CNVs on the post-transcriptional network.
Characterization of the regulation of CNV-miRNAs from the view of their target genes
We first compiled the genes regulated by CNV-miRNAs using the targets from TargetScan5.1, which predicts miRNA targets based on sequence complementarities, sequence context information and binding energy. Because of its high confidence, TargetScan5.1 has been widely used in a variety of “omics” studies (see Methods)[32–34]. From among the miRNA-Target associations that were obtained, the representative miRNA for each family with the lowest total context score was presented, but all other miRNAs from the same family were considered to target the same gene at the same target sites. To study the non-redundant miRNA binding sites directly, we replaced the miRNAs by their miRNA-family ID. Finally, 63,428 regulatory relationships were constructed comprising 541 miRNA-families and 9,174 targets (see Additional file1).
According to the study by Marcinkowska et al., a total of 209 miRNAs were found to locate in the human CNV-regions. These miRNAs belong to 172 families (see Additional file2); the remaining 369 miRNA-families had no members in the CNV-regions. In the following analysis, these two types were referred to as CNV-miRNAs and non-CNV-miRNAs, respectively.
We investigated target genes of the non-CNV-miRNAs and CNV-miRNAs and classified them into three groups (see Additional file3). The first group contains a total of 1,134 target genes that are regulated exclusively by CNV-miRNAs, 823 of the genes are regulated by one CNV-miRNA, 211 by two CNV-miRNAs, 67 by three CNV-miRNAs, 22 by four CNV-miRNAs, and 11 by ≥ 5 CNV-miRNAs. The second group contains a total of 5,710 target genes that are regulated by non-CNV-miRNAs and at least one CNV-miRNA. The third group consists of 2,330 target genes that are regulated exclusively by non-CNV-miRNAs, 1,408 of the genes are regulated by one non-CNV-miRNA, 515 by two non-CNV-miRNAs, 207 by three non-CNV-miRNAs, 95 by four non-CNV-miRNAs and 105 by ≥ 5 non-CNV-miRNAs.
Simulation analysis to explore the target-recognition preference of CNV-miRNAs and non-CNV-miRNAs
The numberof regulatory miRNAs
Mean of 1,000 simulations
Std. dev of 1,000 simulations
Genes regulated exclusively by CNV miRNAs
Genes regulated exclusively by non-CNV miRNAs
Target genes of CNV-miRNAs tend to be differentially expressed among individuals within a population
Intuitively, CNVs of miRNA genes can dramatically change their dosage, and this would then affect the expression levels of the target genes in the corresponding individuals[5, 15]. Recently, a series of genome-wide gene expression profiles have been measured in four HapMap ethnic populations, CEU (U.S. residents with Northern and Western European ancestry), YRI (Yoruba people of Ibadan, Nigeria), CHB (Chinese Han in Beijing) and JPT (Japanese from Tokyo). We calculated the coefficient of variation (CV) for each protein-coding gene across individuals in the four populations to quantify the within-population expression variability of each of the genes (see Methods). Briefly, the CV is the ratio of the standard deviation of gene’s expression to its mean intensity, which is considered to be an unbiased and comprehensive metric to measure the regulation diversity at the expression level among individuals (see Additional file4).
Target genes of CNV-miRNAs are more likely to be differentially expressed between populations
A good study has demonstrated that the within-population expression variability of genes can influence the propensity of their differential expression levels between populations. Here, some CNV-miRNAs may live in different populations; thus, the genes targeted by these CNV-miRNAs are likely to be differentially expressed among individuals within a population and also between different populations.
Using the method described above, we identified genes that were differentially expressed in at least one of the four ethnic populations (see Additional file6). As shown in Figure4B, a similar number of genes were differentially expressed among six population pairs selected from the four ethnic populations. We then investigated whether genes targeted by CNV-miRNAs were over-represented in these differentially expressed genes. As shown in Figure4C, the proportion of differentially expressed genes was 15.7% for targets regulated exclusively by non-CNV-miRNAs, 17.4% for targets regulated by both CNV-miRNAs and non-CNV-miRNAs (p=0.060, Chi-square, two-tail test), the proportion increased further to 21.7% for targets regulated exclusively by CNV-miRNAs (p=0.001, Chi-square, two-tail test).
Target genes of CNV-miRNAs tend to be differentially expressed across tissues and developmental stages
Functional differences between target genes regulated exclusively by CNV-miRNAs and target genes regulated exclusively by non-CNV-miRNAs
It is interesting to know whether or not the orthologs of human CNV-miRNAs were also located in CNV-regions of other animals. We compiled the available CNVs of Pan troglodytes and Mus musculus[26, 27], and then intersected the location of their miRNAs with the coordinates of the CNVs. The results showed that only 21 and eight miRNA-families have members located in CNV-regions in Pan troglodytes and inbred Mus musculus, respectively (see Additional file8). Hence, the human genome contained the highest proportion of CNV-miRNAs, making it the best model to detect the mechanisms and function of CNV-miRNAs.
Animal genomes have the characteristics of dynamics and plasticity, giving them the ability to adapt to changing environmental conditions. Mobile and evolving elements such as telomeres, transposons, and copy number variants have been studied in investigations into the potential effect of environment on genomes. For example, Haasl and Payseur designed a mathematical model to study microsatellite variations, such as the expected distribution of repeat sizes, and the expected squared difference in repeat size among samples; their simulations revealed that microsatellites, especially triplet repeats, provided adaptation facilitators for beneficial evolution of genomes. miRNAs are relatively newly discovered genomic elements, but their post-transcriptional regulation is present early on in metazoan evolution. The number of miRNAs in a genome correlates with the morphological complexity of the animal, indicating that they play roles in evolutionary changes of body structure. It is now widely accepted that an increase in the complexity of gene regulatory mechanisms, at both the genomic and transcriptomic level, drives the appearance of more complex organisms. Two distinct mechanisms of increasing complexity of gene expression, namely, the co-evolution between CNVs and miRNAs, have been recently recognized and studied. Marcinkowska et al. compared the fractions of miRNA loci and the fraction of genome covered by CNVs, and reported that the CNV purification effect was insignificant. Felekkis et al. demonstrated that the number of distinct miRNA types and the average number of miRNA binding sites in genes in CNV regions were significantly higher than genes in non-CNV regions. In this study, we proposed the miRNA-target recognition may play important roles in escape from purification of the CNV-miRNAs that target the same genes. Further analysis revealed that “targeting by CNV-miRNAs” seems to be favored and that the target genes participate in a wide-range of cellular responses to environmental factors. For target genes regulated by one miRNA, CNV-miRNAs tend to target a higher average number of genes than non-CNV-miRNAs. From an evolutionary viewpoint, if the CNV-miRNAs were deleterious and only remained in the genome because they were difficult to remove, then we might expect them to have a tendency to target, on average, a lesser number of genes than non-CNV-miRNAs; furthermore, if the CNV-miRNAs were neutral and their retention attributed to random genetic drift, the CNV-miRNAs and non-CNV-miRNAs should target a similar average number of genes. Therefore, some CNV-miRNAs seems to be beneficial to the organism and “targeting by CNV-miRNAs” may provide positive selective pressure to their target genes.
Our analyses revealed pervasive impacts of CNV on the miRNA-mediated post-transcription regulatory network. Previous studies demonstrated that miRNAs preferentially regulated the hubs of protein interaction and metabolic networks. We here propose that the CNV of miRNAs may fluctuate the dosage balance of signal transduction pathways, metabolic flux or protein complexes[53, 54], leading eventually to individuals of the same population or different populations having different susceptibility to diseases. Although it is difficult to identify these CNV-miRNAs without a comprehensive investigation of health risks among human populations, recent experimental studies have discovered CNV-causing dysregulation of miRNAs that confirmed their roles in disease occurrence. In one study, next-generation sequencing technology was used to explore CNV as a potential mechanism of miRNA mis-expression, the affected miRNA loci were consistently found to be either lost or gained, and their candidate mRNA targets were coordinately dysregulated; the authors demonstrated the structure variation of the miRNA loci clearly characterized the pre-invasive stage of breast cancer. In another study, genetic networks were inferred from miRNA expression in normal and cancer tissues, and cancer networks built from disjointed sub-networks were found to accompany miRNA copy number alterations, such as the amplification of the hsa-miR-17/92 family, the deletion of the hsa-miR-143/145 cluster, and the physical alteration of the hsa-miR-204/30 at the DNA copy number level. The results of these studies clearly demonstrate the feasibility of using the dysregulation of CNV-miRNAs as biological markers for disease screening; indicating that CNV-miRNAs and their targets should be given more attention in studies of human health.
To the best of our knowledge, this is the first genome-wide integrative analysis among human CNVs, miRNAs, their targets and expression variations. Our results will pave the way for future studies for the functional characterization of CNV-miRNAs. This study reveals more clear roles of CNV-miRNAs and is valuable for studying the impact of CNVs on human health.
Compilation of human miRNA target genes
The miRNAs and their predicted targets were taken from TargetScan (http://www.targetscan.org version 5.1)[32, 33]. Targets with a total context score of −0.3 or lower were ignored, where the score quantitatively measure the overall target efficacy. A total of 9,174 targets with at least one conserved 7-mer or 8-mer were selected as reliable miRNA targets (see Additional file1).
Analysis of human gene expression data
The microarray-based gene expression profiles were derived from lymphoblastic cell lines of 270 HapMap individuals (http://www.sanger.ac.uk/humgen/genevar, GSE6536), including 90 samples of YRI (Yoruba people of Ibadan, Nigeria), 90 samples of CEU (U.S. residents with northern and western European ancestry), 45 samples of CHB (Chinese Han in Beijing) and 45 samples of JPT (Japanese from Tokyo)[60, 61]. The annotation table was retrieved from http://www.ncbi.nlm.nih.gov/projects/geo/query/acc.cgi?acc=GPL2507. The RefSeq identifiers were transformed to Ensembl Gene ID through BioMart. Finally, the expression profiles of 16,686 human genes (including 8,636 miRNA targets) across four HapMap populations were complied.
The following formulas were adopted to calculate the coefficient of variation (CV) of gene i in each ethnic population.
The mean intensity M i calculated by
The standard deviation σ i calculated by,
The coefficient of variation CV i calculated by
Where j=1 to n, n represents the number of samples in a population, S ij represents the expression signal of gene i in sample j. Greater CV implies higher expression variability of a gene across individuals within the corresponding population (see Additional file4).
Calculation of MAFs of SNPs in UTRs of human genes
Minor allele frequency (MAF) refers to the frequency at which the less common allele occurs in a given population. SNPs with a minor allele frequency of 5% or greater were targeted by the HapMap project and have been widely employed in Genome Wide Association Studies for complex traits (GWAS)[62, 63].
Where N aa represents the count of individuals who are homozygous for allele1, N Aa represents the count of individuals who are heterozygous, N aa represents the count of individuals who are homozygous for allele2.
Compilation of DERs of human genes
The differential expression ratios (DER) of human genes were obtained from the study by Chen et al. (FitSNPs, http://fitsnps.stanford.edu/download.php). Briefly, the authors downloaded 476 human GEO datasets from the NCBI Gene Expression Omnibus and categorized each GEO dataset into 24 types of comparisons, such as disease state, cell type, metabolism and so on. A total of 4,877 subset-versus-subset comparisons were performed to identify differentially expressed genes with a cutoff of q value ≤ 0.05 by SAM package. For each human gene, the count of GEO datasets in which it was differentially expressed was divided by the count of its measured GEO.
The gene symbols and EntrezGene IDs were transformed to their Ensembl gene IDs using the BioMart program.The Ensembl genes with available DERs were then intersected with the genes that were used for TargetScan5.1 prediction. Finally, the DER values of 9,784 genes that are not regulated by miRNAs and 8,979 target genes of miRNAs were obtained.
Functional analysis of human genes based on gene ontology
The Gene Ontology (GO) has developed three structured controlled vocabularies to describe gene products in terms of their associated biological processes, cellular components and molecular functions. The human gene association file was downloaded from http://www.geneontology.org/gene-associations/. For each GO term, the proportion of annotated genes was compared between the genes regulated exclusively by CNV-miRNAs and the genes regulated exclusively by non-CNV-miRNAs. The p-value was estimated by Fisher’s exact two-tailed test, and a cutoff of p ≤ 0.05 was used to identify the over-represented or under-represented GO terms among the genes that are regulated exclusively by CNV-miRNAs.
The project was started and completed in Dalian Institute of chemical Physics. Computations were performed on a Linux cluster with 50 nodes (Intel 5130, 2.0 GHz CPU, 4G memory, Laboratory of Molecular Modeling and Design, Dalian Institute of Chemical Physics, Chinese Academy of Sciences). Perl (http://perl.org) and R (http://www.r-project.org/) scripts were used for analysis, and can be obtained on request.
Copy number variation
miRNA that is located in copy number variation regions
miRNA that is not located in copy number variation regions
U.S. residents with northern and western European ancestry
Yoruba people of Ibadan, Nigeria
Chinese Han in Beijing
Japanese from Tokyo
The coefficient of variation ratio
Minor allele frequency
This work was supported by funding from “Hundred Talents Program” of Chinese Academy of Sciences and State key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences.
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.View ArticlePubMedGoogle Scholar
- He L, Hannon GJ: MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004, 5: 522-531. 10.1038/nrg1379.View ArticlePubMedGoogle Scholar
- Rosero S, Bravo-Egana V, Jiang Z, Khuri S, Tsinoremas N, Klein D, Sabates E, Correa-Medina M, Ricordi C, Domínguez-Bendala J, Diez J, Pastori RL: MicroRNA signature of the human developing pancreas. BMC Genomics. 2010, 11: 509-10.1186/1471-2164-11-509.PubMed CentralView ArticlePubMedGoogle Scholar
- Ding XC, Grosshans H: Repression of C. elegans microRNA targets at the initiation level of translation requires GW182 proteins. EMBO J. 2009, 28: 213-222. 10.1038/emboj.2008.275.PubMed CentralView ArticlePubMedGoogle Scholar
- Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM: Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005, 433: 769-773. 10.1038/nature03315.View ArticlePubMedGoogle Scholar
- Vivek J, Mark L, David DF M, Yang YH: Identification of microRNA-mRNA modules using microarray data. BMC Genomics. 2011, 12: 138-10.1186/1471-2164-12-138.View ArticleGoogle Scholar
- Yu Z, Jian Z, Shen SH, Purisima E, Wang E: Global analysis of microRNA target gene expression reveals that miRNA targets are lower expressed in mature mouse and drosophila tissues than in the embryos. Nucleic Acids Res. 2007, 35: 152-164.PubMed CentralView ArticlePubMedGoogle Scholar
- Hornstein E, Shomron N: Canalization of development by microRNAs. Nat Genet. 2006, 38: S20-S24. 10.1038/ng1803.View ArticlePubMedGoogle Scholar
- Li Y, Wang F, Lee JA, Gao FB: MicroRNA-9a ensures the precise specification of sensory organ precursors in Drosophila. Genes Dev. 2006, 20: 2793-2805. 10.1101/gad.1466306.PubMed CentralView ArticlePubMedGoogle Scholar
- Cohen SM, Brennecke J, Stark A: Denoising feedback loops by thresholding – a new role for microRNAs. Genes Dev. 2006, 20: 2769-2772. 10.1101/gad.1484606.View ArticlePubMedGoogle Scholar
- O’Donnell KA, Wentzel EA, Zeller KI, Dang CV, Mendell JT: c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 2005, 435: 839-843. 10.1038/nature03677.View ArticlePubMedGoogle Scholar
- Morley M, Molony CM, Weber TM, Devlin JL, Ewens KG, Spielman RS, Cheung VG: Genetic analysis of genome-wide variation in human gene expression. Nature. 2004, 430: 743-747. 10.1038/nature02797.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheung VG, Spielman RS, Ewens KG, Weber TM, Morley M, Burdick JT: Mapping determinants of human gene expression by regional and genome-wide association. Nature. 2005, 437: 1365-1369. 10.1038/nature04244.PubMed CentralView ArticlePubMedGoogle Scholar
- GuhaThakurta D, Xie T, Anand M, Edwards SW, Li G, Wang SS, Schadt EE: Cis-regulatory variations: a study of SNPs around genes showing cis-linkage in segregating mouse populations. BMC Genomics. 2006, 7: 235-10.1186/1471-2164-7-235.PubMed CentralView ArticlePubMedGoogle Scholar
- Henrichsen CN, Chaignat E, Reymond A: Copy number variants, diseases and gene expression. Hum Mol Genet. 2009, 18 (R1): R1-R8. 10.1093/hmg/ddp011.View ArticlePubMedGoogle Scholar
- Pickrell JK, Marioni JC, Pai AA, Degner JF, Engelhardt BE, Nkadori E, Veyrieras JB, Stephens M, Gilad Y, Pritchard JK: Understanding mechanisms underlying human gene expression variation with RNA sequencing. Nature. 2010, 464: 768-772. 10.1038/nature08872.PubMed CentralView ArticlePubMedGoogle Scholar
- Wong KK, deLeeuw RJ, Dosanjh NS, Kimm LR, Cheng Z, Horsman DE, MacAulay C, Ng RT, Brown CJ, Eichler EE, Lam WL: A comprehensive analysis of common copy-number variations in the human genome. Am J Hum Genet. 2007, 80: 91-104. 10.1086/510560.PubMed CentralView ArticlePubMedGoogle Scholar
- Bonaglia MC, Giorda R, Beri S, De Agostini C, Novara F, Fichera M, Grillo L, Galesi O, Vetro A, Ciccone R, Bonati MT, Giglio S, Guerrini R, Osimani S, Marelli S, Zucca C, Grasso R, Borgatti R, Mani E, Motta C, Molteni M, Romano C, Greco D, Reitano S, Baroncini A, Lapi E, Cecconi A, Arrigo G, Patricelli MG, Pantaleoni C, D’Arrigo S, Riva D, Sciacca F, Dalla Bernardina B, Zoccante L, Darra F, Termine C, Maserati E, Bigoni S, Priolo E, Bottani A, Gimelli S, Bena F, Brusco A, di Gregorio E, Bagnasco I, Giussani U, Nitsch L, Politi P, Martinez-Frias ML, Martínez-Fernández ML, Martínez Guardia N, Bremer A, Anderlid BM, Zuffardi O: Molecular mechanisms generating and stabilizing terminal 22q13 deletions in 44 subjects with Phelan/McDermid Syndrome. PLoS Genet. 2011, 7: e1002173-10.1371/journal.pgen.1002173.PubMed CentralView ArticlePubMedGoogle Scholar
- Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P, Fitzgerald T, Hu M, Ihm CH, Kristiansson K, Macarthur DG, Macdonald JR, Onyiah I, Pang AW, Robson S, Stirrups K, Valsesia A, Walter K, Wei J, Tyler-Smith C, Carter NP, Lee C, Scherer SW, Hurles ME, Wellcome Trust Case Control Consortium: Origins and functional impact of copy number variation in the human genome. Nature. 2010, 464: 704-712. 10.1038/nature08516.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang RT, Sangtae A, Park CC, Khan AH, Kenneth L, Smith DJ: Effects of genome-wide copy number variation on expression in mammalian cells. BMC Genomics. 2011, 12: 562-10.1186/1471-2164-12-562.PubMed CentralView ArticlePubMedGoogle Scholar
- Woodwark C, Bateman A: The characterization of three types of genes that overlie copy number variable regions. PLoS One. 2011, 6 (5): e14814-10.1371/journal.pone.0014814.PubMed CentralView ArticlePubMedGoogle Scholar
- Korbel JO, Urban AE, Affourtit JP, Godwin B, Grubert F, Simons JF, Kim PM, Palejev D, Carriero NJ, Du L, Taillon BE, Chen Z, Tanzer A, Saunders AC, Chi J, Yang F, Carter NP, Hurles ME, Weissman SM, Harkins TT, Gerstein MB, Egholm M, Snyder M: Paired-end mapping reveals extensive structural variation in the human genome. Science. 2007, 318: 420-426. 10.1126/science.1149504.PubMed CentralView ArticlePubMedGoogle Scholar
- Sudmant PH, Kitzman JO, Antonacci F, Alkan C, Malig M, Tsalenko A, Sampas N, Bruhn L, Shendure J, Eichler EE, 1000 Genomes Project: Diversity of human copy number variation and multicopy genes. Science. 2010, 330: 641-646. 10.1126/science.1197005.PubMed CentralView ArticlePubMedGoogle Scholar
- Mills RE, Walter K, Stewart C, Handsaker RE, Chen K, Alkan C, Abyzov A, Yoon SC, Ye K, Cheetham RK, Chinwalla A, Conrad DF, Fu Y, Grubert F, Hajirasouliha I, Hormozdiari F, Iakoucheva LM, Iqbal Z, Kang S, Kidd JM, Konkel MK, Korn J, Khurana E, Kural D, Lam HY, Leng J, Li R, Li Y, Lin CY, Luo R, Mu XJ, Nemesh J, Peckham HE, Rausch T, Scally A, Shi X, Stromberg MP, Stütz AM, Urban AE, Walker JA, Wu J, Zhang Y, Zhang ZD, Batzer MA, Ding L, Marth GT, McVean G, Sebat J, Snyder M, Wang J, Ye K, Eichler EE, Gerstein MB, Hurles ME, Lee C, McCarroll SA, Korbel JO, 1000 Genomes Project: Mapping copy number variation by population-scale genome sequencing. Nature. 2011, 470: 59-65. 10.1038/nature09708.PubMed CentralView ArticlePubMedGoogle Scholar
- Perry GH, Yang F, Marques-Bonet T, Murphy C, Fitzgerald T, Lee AS, Hyland C, Stone AC, Hurles ME, Tyler-Smith C, Eichler EE, Carter NP, Lee C, Redon R: Copy number variation and evolution in humans and chimpanzees. Genome Res. 2008, 18: 1698-1710. 10.1101/gr.082016.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Cutler G, Marshall LA, Chin N, Baribault H, Kassner PD: Significant gene content variation characterizes the genomes of inbred mouse strains. Genome Res. 2007, 17: 1743-1754. 10.1101/gr.6754607.PubMed CentralView ArticlePubMedGoogle Scholar
- Agam A, Yalcin B, Bhomra A, Cubin M, Webber C, Holmes C, Flint J, Mott R: Elusive copy number variation in the mouse genome. PLoS One. 2010, 5 (9): e12839-10.1371/journal.pone.0012839.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang L, Huang J, Yang N, Greshock J, Megraw MS, Giannakakis A, Liang S, Naylor TL, Barchetti A, Ward MR, Yao G, Medina A, O’brien-Jenkins A, Katsaros D, Hatzigeorgiou A, Gimotty PA, Weber BL, Coukos G: microRNAs exhibit high frequency genomic alterations in human cancer. Proc Natl Acad Sci USA. 2006, 103: 9136-9141. 10.1073/pnas.0508889103.PubMed CentralView ArticlePubMedGoogle Scholar
- Lionetti M, Agnelli L, Mosca L, Fabris S, Andronache A, Todoerti K, Ronchetti D, Deliliers GL, Neri A: Integrative high-resolution microarray analysis of human myeloma cell lines reveals deregulated miRNA expression associated with allelic imbalances and gene expression profiles. Genes Chromosomes Cancer. 2009, 48: 521-531. 10.1002/gcc.20660.View ArticlePubMedGoogle Scholar
- Maire G, Martin JW, Yoshimoto M, Chilton-MacNeill S, Zielenska M, Squire JA: Analysis of miRNA-gene expression-genomic profiles reveals complex mechanisms of microRNA deregulation in osteosarcoma. Cancer Genet. 2011, 204: 138-146. 10.1016/j.cancergen.2010.12.012.View ArticlePubMedGoogle Scholar
- Marcinkowska M, Szymanski M, Krzyzosiak WJ, Kozlowski P: Copy number variation of microRNA genes in the human genome. BMC Genomics. 2011, 12: 183-10.1186/1471-2164-12-183.PubMed CentralView ArticlePubMedGoogle Scholar
- Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120: 15-20. 10.1016/j.cell.2004.12.035.View ArticlePubMedGoogle Scholar
- Chen K, Rajewsky N: Natural selection on human microRNA binding sites inferred from SNP data. Nat Genet. 2006, 38: 1452-1456. 10.1038/ng1910.View ArticlePubMedGoogle 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: 91-105. 10.1016/j.molcel.2007.06.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Fay JC, Wyckoff GJ, Wu CI: Positive and negative selection on the human genome. Genetics. 2001, 158: 1227-1234.PubMed CentralPubMedGoogle Scholar
- Nielsen R, Hellmann I, Hubisz M, Bustamante C, Clark AG: Recent and ongoing selection in the human genome. Nat Rev Genet. 2007, 8: 857-868.PubMed CentralView ArticlePubMedGoogle Scholar
- Felekkis K, Voskarides K, Dweep H, Sticht C, Gretz N, Deltas C: Increased number of microRNA target sites in genes encoded in CNV regions, Evidence for an evolutionary genomic interaction. Mol Biol Evol. 2011, 28: 2421-2424. 10.1093/molbev/msr078.View ArticlePubMedGoogle Scholar
- Kaern M, Elston TC, Blake WJ, Collins JJ: Stochasticity in gene expression: from theories to phenotypes. Nat Rev Gene. 2005, 6: 451-464. 10.1038/nrg1615.View ArticleGoogle Scholar
- Hartl D: A Primer of Population Genetics. 2000, Sunderland, MA, USA: Sinauer Associates, Inc., 3Google Scholar
- Smedley D, Haider S, Ballester B, Holland R, London D, Thorisson G, Kasprzyk A: BioMart-biological queries made easy. BMC Genomics. 2009, 10: 22-10.1186/1471-2164-10-22.PubMed CentralView ArticlePubMedGoogle Scholar
- The International HapMap Consortium: Integrating common and rare genetic variation in diverse human populations. Nature. 2010, 467: 52-58. 10.1038/nature09298.View ArticleGoogle Scholar
- Li J, Liu Y, Kim T, Min R, Zhang Z: Gene expression variability within and between human populations and implications toward disease susceptibility. PLoS Comput Biol. 2010, 6 (8): e1000910-10.1371/journal.pcbi.1000910.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen R, Morgan AA, Dudley J, Deshpande T, Li L, Kodama K, Chiang AP, Butte AJ: FitSNPs: highly differentially expressed genes are more likely to have variants associated with disease. Genome Biol. 2008, 9: R170-10.1186/gb-2008-9-12-r170.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen R, Li L, Butte AJ: AILUN: reannotating gene expression data automatically. Nat Methods. 2007, 4: 879-10.1038/nmeth1107-879.PubMed CentralView ArticlePubMedGoogle Scholar
- Day-Richter J, Harris MA, Haendel M, Lewis S, Gene Ontology OBO-Edit Working Group: OBO-Edit–an ontology editor for biologists. Bioinformatics. 2007, 23: 2198-2200. 10.1093/bioinformatics/btm112.View ArticlePubMedGoogle Scholar
- Haasl RJ, Payseur BA: The number of alleles at a microsatellite defines the allele frequency spectrum and facilitates fast accurate estimation of theta. Mol Biol Evol. 2010, 12: 2702-2715.View ArticleGoogle Scholar
- Sempere LF, Cole CN, McPeek MA, Peterson KJ: The phylogenetic distribution of metazoan microRNAs: insights into evolutionary complexity and constraint. J Exp Zool B Mol Dev Evol. 2006, 306: 575-588.View ArticlePubMedGoogle Scholar
- Heimberg AM, Sempere LF, Moy VN, Donoghue PC, Peterson KJ: MicroRNAs and the advent of vertebrate morphological complexity. Proc Natl Acad Sci USA. 2008, 105 (8): 2946-2950. 10.1073/pnas.0712259105.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu CI, Shen Y, Tang T: Evolution under canalization and the dual roles of microRNAs–A hypothesis. Genome Res. 2009, 19 (5): 734-743. 10.1101/gr.084640.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou J, Lemos B, Dopman EB, Hartl DL: Copy-number variation: the balance between gene dosage and expression in drosophila melanogaster. Genome Biol Evol. 2011, 3: 1014-1024. 10.1093/gbe/evr023.PubMed CentralView ArticlePubMedGoogle Scholar
- Liang H, Li WH: MicroRNA regulation of human protein–protein interaction network. RNA. 2007, 13 (9): 1402-1408. 10.1261/rna.634607.PubMed CentralView ArticlePubMedGoogle Scholar
- Tibiche C, Wang E: MicroRNA regulatory patterns on the human metabolic network. The Open Systems Biology Journal. 2008, 1: 1-8.View ArticleGoogle Scholar
- Veitia RA: Gene dosage balance in cellular pathways: implications for dominance and gene duplicability. Genetics. 2004, 168: 569-574. 10.1534/genetics.104.029785.PubMed CentralView ArticlePubMedGoogle Scholar
- Veitia RA, Bottani S, Birchler JA: Cellular reactions to gene dosage imbalance: genomic, transcriptomic and proteomic effects. Trends Genet. 2008, 24: 390-397. 10.1016/j.tig.2008.05.005.View ArticlePubMedGoogle Scholar
- Knight JC: Human Genetic Diversity: Functional Consequences for Health and Disease. 2009, Oxford, UK: Oxford University Press, 1View ArticleGoogle Scholar
- Bethany Noelle Hannafon: An integrated analysis of the coordinated dysregulation of microRNAs and their targets in pre-invasive breast cancer. PhD thesis. 2010, Boston UniversityGoogle Scholar
- Volinia S, Galasso M, Costinean S, Tagliavini L, Gamberoni G, Drusco A, Marchesini J, Mascellani N, Sana ME, Abu Jarour R, Desponts C, Teitell M, Baffa R, Aqeilan R, Iorio MV, Taccioli C, Garzon R, Di Leva G, Fabbri M, Catozzi M, Previati M, Ambs S, Palumbo T, Garofalo M, Veronese A, Bottoni A, Gasparini P, Harris CC, Visone R, Pekarsky Y, de la Chapelle A, Bloomston M, Dillhoff M, Rassenti LZ, Kipps TJ, Huebner K, Pichiorri F, Lenze D, Cairo S, Buendia MA, Pineau P, Dejean A, Zanesi N, Rossi S, Calin GA, Liu CG, Palatini J, Negrini M, Vecchione A, Rosenberg A, Croce CM: Reprogramming of miRNA networks in cancer and leukemia. Genome Res. 2010, 20 (5): 589-599. 10.1101/gr.098046.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP: The impact of microRNAs on protein output. Nature. 2008, 455: 64-71. 10.1038/nature07242.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu X, Song Y: Preferential regulation of miRNA targets by environmental chemicals in the human genome. BMC Genomics. 2011, 12: 244-10.1186/1471-2164-12-244.PubMed CentralView ArticlePubMedGoogle Scholar
- Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, Redon R, Bird CP, de Grassi A, Lee C, Tyler-Smith C, Carter N, Scherer SW, Tavaré S, Deloukas P, Hurles ME, Dermitzakis ET: Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007, 315: 848-853. 10.1126/science.1136678.PubMed CentralView ArticlePubMedGoogle Scholar
- Stranger BE, Nica AC, Forrest MS, Dimas A, Bird CP, Beazley C, Ingle CE, Dunning M, Flicek P, Koller D, Montgomery S, Tavaré S, Deloukas P, Dermitzakis ET: Population genomics of human gene expression. Nat Genet. 2007, 39: 1217-1224. 10.1038/ng2142.PubMed CentralView ArticlePubMedGoogle Scholar
- Serre D, Gurd S, Ge B, Sladek R, Sinnett D, Harmsen E, Bibikova M, Chudin E, Barker DL, Dickinson T, Fan JB, Hudson TJ: Differential allelic expression in the human genome: a robust approach to identify genetic and epigenetic cis-acting mechanisms regulating gene expression. PLoS Genet. 2008, 4 (2): e1000006-10.1371/journal.pgen.1000006.PubMed CentralView ArticlePubMedGoogle Scholar
- Spencer CC, Su Z, Donnelly P, Marchini J: Designing genome-wide association studies: sample size, power, imputation, and the choice of genotyping chip. PLoS Gene. 2009, 5 (5): e1000477-10.1371/journal.pgen.1000477.View ArticleGoogle Scholar