Minja Zorc

Catalog of MicroRNA Seed Polymorphisms in Vertebrates

MicroRNAs (miRNAs) are a class of non-coding RNA that plays an important role in posttranscriptional regulation of mRNA. Evidence has shown that miRNA gene variability might interfere with its function resulting in phenotypic variation and disease susceptibility. A major role in miRNA target recognition is ascribed to complementarity with the miRNA seed region that can be affected by polymorphisms. In the present study, we developed an online tool for the detection of miRNA polymorphisms (miRNA SNiPer) in vertebrates (http://www.integratomics-time.com/miRNA-SNiPer) and generated a catalog of miRNA seed region polymorphisms (miR-seed-SNPs) consisting of 149 SNPs in six species. Although a majority of detected polymorphisms were due to point mutations, two consecutive nucleotide substitutions (double nucleotide polymorphisms, DNPs) were also identified in nine miRNAs. We determined that miR-SNPs are frequently located within the quantitative trait loci (QTL), chromosome fragile sites, and cancer susceptibility loci, indicating their potential role in the genetic control of various complex traits. To test this further, we performed an association analysis between the mmu-miR-717 seed SNP rs30372501, which is polymorphic in a large number of standard inbred strains, and all phenotypic traits in these strains deposited in the Mouse Phenome Database. Analysis showed a significant association between the mmu-miR-717 seed SNP and a diverse array of traits including behavior, blood-clinical chemistry, body weight size and growth, and immune system suggesting that seed SNPs can indeed have major pleiotropic effects. The bioinformatics analyses, data and tools developed in the present study can serve researchers as a starting point in testing more targeted hypotheses and designing experiments using optimal species or strains for further mechanistic studies.

Genome-wide and species-wide in silico screening for intragenic MicroRNAs in human, mouse and chicken.

MicroRNAs (miRNAs) are non-coding RNAs (ncRNAs) involved in regulation of gene expression. Intragenic miRNAs, especially those exhibiting a high degree of evolutionary conservation, have been shown to be coordinately regulated and/or expressed with their host genes, either with synergistic or antagonistic correlation patterns. However, the degree of cross-species conservation of miRNA/host gene co-location is not known and co-expression information is incomplete and fragmented among several studies. Using the genomic resources (miRBase and Ensembl) we performed a genome-wide in silico screening (GWISS) for miRNA/host gene pairs in three well-annotated vertebrate species: human, mouse, and chicken. Approximately half of currently annotated miRNA genes resided within host genes: 53.0% (849/1,600) in human, 48.8% (418/855) in mouse, and 42.0% (210/499) in chicken, which we present in a central publicly available Catalog of intragenic miRNAs (http://www.integratomics-time.com/miR-host/catalog). The miRNA genes resided within either protein-coding or ncRNA genes, which include long intergenic ncRNAs (lincRNAs) and small nucleolar RNAs (snoRNAs). Twenty-seven miRNA genes were found to be located within the same host genes in all three species and the data integration from literature and databases showed that most (26/27) have been found to be co-expressed. Particularly interesting are miRNA genes located within genes encoding for miRNA silencing machinery (DGCR8, DICER1, and SND1 in human and Cnot3, Gdcr8, Eif4e, Tnrc6b, and Xpo5 in mouse). We furthermore discuss a potential for phenotype misattribution of miRNA host gene polymorphism or gene modification studies due to possible collateral effects on miRNAs hosted within them. In conclusion, the catalog of intragenic miRNAs and identified 27 miRNA/host gene pairs with cross-species conserved co-location, co-expression, and potential co-regulation, provide excellent candidates for further functional annotation of intragenic miRNAs in health and disease.

Obesity gene atlas in mammals.

Obesity in humans has increased at an alarming rate over the past two decades and has become one of the leading public health problems worldwide. Studies have revealed a large number of genes/markers that are associated with obesity and/or obesity-related phenotypes, indicating an urgent need to develop a central database for helping the community understand the genetic complexity of obesity. In the present study, we collected a total of 1,736 obesity associated loci and created a freely available obesity database, including 1,515 protein-coding genes and 221 microRNAs (miRNAs) collected from four mammalian species: human, cattle, rat, and mouse. These loci were integrated as orthologs on comparative genomic views in human, cattle, and mouse. The database and genomic views are freely available online at: http://www.integratomics-time.com/fat_deposition. Bioinformatics analyses of the collected data revealed some potential novel obesity related molecular markers which represent focal points for testing more targeted hypotheses and designing experiments for further studies. We believe that this centralized database on obesity and adipogenesis will facilitate development of comparative systems biology approaches to address this important health issue in human and their potential applications in animals.

Genome-wide in silico screening for microRNA genetic variability in livestock species

MicroRNAs are a class of non-coding RNAs that post-transcriptionally regulate target gene expression. Previous studies have shown that microRNA gene variability can interfere with its function, resulting in phenotypic variation. Polymorphisms within microRNA genes present a source of novel biomarkers for phenotypic traits in animal breeding. However, little is known about microRNA genetic variability in livestock species, which is also due to incomplete data in genomic resource databases. Therefore, the aim of this study was to perform a genome-wide in silico screening of genomic sources and determine the genetic variability of microRNA genes in livestock species using mirna sniper 3.0 (http://www.integratomics-time.com/miRNA-SNiPer/), a new version of our previously developed tool. By examining Ensembl and miRBase genome builds, it was possible to design a tool-based generated search of 16 genomes including four livestock species: pig, horse, cattle and chicken. The analysis revealed 65 polymorphisms located within mature microRNA regions in these four species, including 28% within the seed region in cattle and chicken. Polymorphic microRNA genes in cattle and chicken were further examined for mapping to quantitative trait loci regions associated with production and health traits. The developed bioinformatics tool enables the analysis of polymorphic microRNA genes and prioritization of potential regulatory polymorphisms and therefore contributes to the development of microRNA-based biomarkers in livestock species. The assembled catalog and the developed tool can serve the animal science community to efficiently select microRNA SNPs for further quantitative and molecular genetic evaluations of their phenotypic effects and causal associations with livestock production traits.

Pivotal role of the muscle-contraction pathway in cryptorchidism and evidence for genomic connections with cardiomyopathy pathways in RASopathies

BackgroundCryptorchidism is the most frequent congenital disorder in male children; however the genetic causes of cryptorchidism remain poorly investigated. Comparative integratomics combined with systems biology approach was employed to elucidate genetic factors and molecular pathways underlying testis descent.MethodsLiterature mining was performed to collect genomic loci associated with cryptorchidism in seven mammalian species. Information regarding the collected candidate genes was stored in MySQL relational database. Genomic view of the loci was presented using Flash GViewer web tool (http://gmod.org/wiki/Flashgviewer/). DAVID Bioinformatics Resources 6.7 was used for pathway enrichment analysis. Cytoscape plug-in PiNGO 1.11 was employed for protein-network-based prediction of novel candidate genes. Relevant protein-protein interactions were confirmed and visualized using the STRING database (version 9.0).ResultsThe developed cryptorchidism gene atlas includes 217 candidate loci (genes, regions involved in chromosomal mutations, and copy number variations) identified at the genomic, transcriptomic, and proteomic level. Human orthologs of the collected candidate loci were presented using a genomic map viewer. The cryptorchidism gene atlas is freely available online: http://www.integratomics-time.com/cryptorchidism/. Pathway analysis suggested the presence of twelve enriched pathways associated with the list of 179 literature-derived candidate genes. Additionally, a list of 43 network-predicted novel candidate genes was significantly associated with four enriched pathways. Joint pathway analysis of the collected and predicted candidate genes revealed the pivotal importance of the muscle-contraction pathway in cryptorchidism and evidence for genomic associations with cardiomyopathy pathways in RASopathies.ConclusionsThe developed gene atlas represents an important resource for the scientific community researching genetics of cryptorchidism. The collected data will further facilitate development of novel genetic markers and could be of interest for functional studies in animals and human. The proposed network-based systems biology approach elucidates molecular mechanisms underlying co-presence of cryptorchidism and cardiomyopathy in RASopathies. Such approach could also aid in molecular explanation of co-presence of diverse and apparently unrelated clinical manifestations in other syndromes.

Genetic Variability of MicroRNA Genes in 15 Animal Species.

MicroRNAs (miRNA) are a class of non-coding RNAs important in posttranscriptional regulation of target genes. Previous studies have proven that genetic variability of miRNA genes (miR-SNP) has an impact on phenotypic variation and disease susceptibility in human, mice and some livestock species. MicroRNA gene polymorphisms could therefore represent biomarkers for phenotypic traits also in other animal species. We upgraded our previously developed tool miRNA SNiPer to the version 4.0 which enables the search of miRNA genetic variability in 15 animal genomes: http://www.integratomics-time.com/miRNA-SNiPer. Genome-wide in silico screening (GWISS) of 15 genomes revealed that based on the current database releases, miRNA genes are most polymorphic in cattle, followed by human, fruitfly, mouse, chicken, pig, horse, and sheep. The difference in the number of miRNA gene polymorphisms between species is most probably not due to a biological reason and lack of genetic variability in some species, but to different stage of sequencing projects and differences in development of genomic resource databases in different species. Genome screening revealed several interesting genomic hotspots. For instance, several multiple nucleotide polymorphisms (MNPs) are present within mature seed region in cattle. Among miR-SNPs 46 are present on commercial whole-genome SNP chips: 16 in cattle, 26 in chicken, two in sheep and two in pig. The update of the miRNA SNiPer tool and the generated catalogs will serve researchers as a starting point in designing projects dealing with the effects of genetic variability of miRNA genes.

MicroRNA Polymorphisms in Cancer: A Literature Analysis

Single nucleotide polymorphisms (SNPs) located in microRNA (miRNA) genes (miR-SNPs) have attracted increasing attention in recent years due to their involvement in the development of various types of cancer. Therefore, a systematic review on this topic was needed. From 55 scientific publications we collected 20 SNPs, which are located within 18 miRNA encoding genes and have been associated with 16 types of cancer. Among 20 miRNA gene polymorphisms 13 are located within the premature miRNA region, five within mature, and two within mature seed miRNA region. We graphically visualized a network of miRNA-cancer associations which revealed miRNA genes and cancer types with the highest number of connections. Our study showed that, despite a large number of variations currently known to be located within miRNA genes in humans, most of them have not yet been tested for association with cancer. MicroRNA SNPs collected in this study represent only 0.43% of known miRNA gene variations (20/4687). Results of the present study will be useful to researchers investigating the clinical use of miRNAs, such as the roles of miRNAs as diagnostic markers and therapeutic targets.

Catalog of genetic variants within mature microRNA seed regions in chicken

MicroRNA (miRNA) is a class of noncoding RNA important in posttranscriptional regulation of target genes. The regulation mechanism requires complementarity between target mRNA and the miRNA region responsible for their recognition and binding, also called the seed region. It has been estimated that each miRNA targets approximately 200 genes and genetic variability of miRNA genes has been associated with phenotypic variation and disease susceptibility in humans, livestock species, and model organisms. Polymorphisms in miRNA genes especially within the seed region could therefore represent biomarkers for phenotypic traits important in livestock animals. Using the updated Version 5.0 of our previously developed bioinformatics tool miRNA SNiPer we assembled polymorphic miRNA genes in chicken. Out of 740 miRNA genes 263 were polymorphic, among them 77 had SNPs located within the mature region, and 29 of them within the miRNA seed region. Because several polymorphisms in databases result from sequencing errors, we performed experimental validation of polymorphisms located within 4 selected miRNA genes in chicken (gga-mir-1614, −1644, −1648, and −1657). We confirmed the presence of nine polymorphisms and identified 3 additional novel polymorphisms within primary miRNA regions in chicken representing 3 layer-type breeds, one layer-type hybrid, and one meat-type intercrossed population. The developed catalog of mir-SNPs in chicken can serve researchers as a starting point for association studies dealing with poultry production traits and designing functional experiments.

Immune related gene expression in worker honey bee (Apis mellifera carnica) pupae exposed to neonicotinoid thiamethoxam and Varroa mites (Varroa destructor)

Varroa destructor is one of the most common parasites of honey bee colonies and is considered as a possible co-factor for honey bee decline. At the same time, the use of pesticides in intensive agriculture is still the most effective method of pest control. There is limited information about the effects of pesticide exposure on parasitized honey bees. Larval ingestion of certain pesticides could have effects on honey bee immune defense mechanisms, development and metabolic pathways. Europe and America face the disturbing phenomenon of the disappearance of honey bee colonies, termed Colony Collapse Disorder (CCD). One reason discussed is the possible suppression of honey bee immune system as a consequence of prolonged exposure to chemicals. In this study, the effects of the neonicotinoid thiamethoxam on honey bee, Apis mellifera carnica, pupae infested with Varroa destructor mites were analyzed at the molecular level. Varroa-infested and non-infested honey bee colonies received protein cakes with or without thiamethoxam. Nurse bees used these cakes as a feed for developing larvae. Samples of white-eyed and brown-eyed pupae were collected. Expression of 17 immune-related genes was analyzed by real-time PCR. Relative gene expression in samples exposed only to Varroa or to thiamethoxam or simultaneously to both Varroa and thiamethoxam was compared. The impact from the consumption of thiamethoxam during the larval stage on honey bee immune related gene expression in Varroa-infested white-eyed pupae was reflected as down-regulation of spaetzle, AMPs abaecin and defensin-1 and up-regulation of lysozyme-2. In brown-eyed pupae up-regulation of PPOact, spaetzle, hopscotch and basket genes was detected. Moreover, we observed a major difference in immune response to Varroa infestation between white-eyed pupae and brown-eyed pupae. The majority of tested immune-related genes were upregulated only in brown-eyed pupae, while in white-eyed pupae they were downregulated.