A Single-Cell Transcriptome of Bovine Milk Somatic Cells.

The production of milk by dairy cows far exceeds the nutritional needs of the calf and is vital for the economical use of dairy cattle. High milk yield is a unique production trait that can be effectively enhanced through traditional selection methods. The process of lactation in cows serves as an excellent model for studying the biological aspects of lactation with the aim of exploring the mechanistic base of this complex trait at the cellular level. In this study, we analyzed the milk transcriptome at the single-cell level by conducting scRNA-seq analysis on milk samples from two Holstein Friesian cows at mid-lactation (75 and 93 days) using the 10× Chromium platform. Cells were pelleted and fat was removed from milk by centrifugation. The cell suspension from each cow was loaded on separate channels, resulting in the recovery of 9313 and 14,544 cells. Library samples were loaded onto two lanes of the NovaSeq 6000 (Illumina) instrument. After filtering at the cell and gene levels, a total of 7988 and 13,973 cells remained, respectively. We were able to reconstruct different cell types (milk-producing cells, progenitor cells, macrophages, monocytes, dendritic cells, T cells, B cells, mast cells, and neutrophils) in bovine milk. Our findings provide a valuable resource for identifying regulatory elements associated with various functions of the mammary gland such as lactation, tissue renewal, native immunity, protein and fat synthesis, and hormonal response.

Selection Signatures Reveal Candidate Genes for the Cornish Rex Breed-Specific Phenotype.

Many coat color, behavioral and morphological traits are specific and fixed across cat breeds, with several variants influencing these traits being common among different breeds. In the domestic cat, rexoid mutations have been documented in several breeds. In the Cornish Rex, four bp deletion in the LPAR6 gene has been found to cause a frame shift and a premature stop codon. In addition to the rexoid coat, Cornish Rex cats also have a characteristic head, ear shape and body type. Analysis of the selection signatures in the Cornish Rex genome revealed several regions that are under selective pressure. One of these is located in CFA B4, in the region where the ALX1 gene is located. The ALX1 gene in Burmese cats disrupts the cranial morphogenesis and causes brachycephaly in the heterozygous state. In our study, we confirmed the presence of a deletion in LPAR6 in 20 Cornish Rex and in four F1 hybrids between Cornish Rex and domestic cat. However, we did not confirm the presence of the deletion in ALX1 in Cornish Rex cats. Genome-wide selection signature analysis was performed using ROH islands and integrated haplotype score (iHS) statistics based on publicly available SNP array data of 11 Cornish Rex cats. The selection signatures were detected on chromosomes A1, A3, C2, B1, B4 and D1.

Genome-Wide Analysis of Milk Production Traits and Selection Signatures in Serbian Holstein-Friesian Cattle.

To improve the genomic evaluation of milk-related traits in Holstein-Friesian (HF) cattle it is essential to identify the associated candidate genes. Novel SNP-based analyses, such as the genetic mapping of inherited diseases, GWAS, and genomic selection, have led to a new era of research. The aim of this study was to analyze the association of each individual SNP in Serbian HF cattle with milk production traits and inbreeding levels. The SNP 60 K chip Axiom Bovine BovMDv3 was deployed for the genotyping of 334 HF cows. The obtained genomic results, together with the collected phenotypic data, were used for a GWAS. Moreover, the identification of ROH segments was performed and served for inbreeding coefficient evaluation and ROH island detection. Using a GWAS, a polymorphism, rs110619097 (located in the intron of the CTNNA3 gene), was detected to be significantly (p < 0.01) associated with the milk protein concentration in the first lactation (adjusted to 305 days). The average genomic inbreeding value (FROH) was 0.079. ROH islands were discovered in proximity to genes associated with milk production traits and genomic regions under selection pressure for other economically important traits of dairy cattle. The findings of this pilot study provide useful information for a better understanding of the genetic architecture of milk production traits in Serbian HF dairy cows and can be used to improve lactation performances in Serbian HF cattle breeding programs.

Identification of potentially pathogenic variants for autism spectrum disorders using gene-burden analysis.

Gene- burden analyses have lately become a very successful way for the identification of genes carrying risk variants underlying the analysed disease. This approach is also suitable for complex disorders like autism spectrum disorder (ASD). The gene-burden analysis using Testing Rare Variants with Public Data (TRAPD) software was conducted on whole exome sequencing data of Slovenian patients with ASD to determine potentially novel disease risk variants in known ASD-associated genes as well as in others. To choose the right control group for testing, principal component analysis based on the 1000 Genomes and ASD cohort samples was conducted. The subsequent protein structure and ligand binding analysis usingI-TASSER package were performed to detect changes in protein structure and ligand binding to determine a potential pathogenic consequence of observed mutation. The obtained results demonstrate an association of two variants-p.Glu198Lys (PPP2R5D:c.592G>A) and p.Arg253Gln (PPP2R5D:c.758G>A) with the ASD. Substitution p.Glu198Lys (PPP2R5D:c.592G>A) is a variant, previously described as pathogenic in association with ASD combined with intellectual disability, whereas p.Arg253Gln (PPP2R5D:c.758G>A) has not been described as an ASD-associated pathogenic variant yet. The results indicate that the filtering process was suitable and could be used in the future for detection of novel pathogenic variants when analysing groups of ASD patients.

Genetic Diversity of Montenegrin Local Sheep Breeds Based on Microsatellite Markers.

The Montenegrin sheep population mostly consists of local breeds and their crossbreeds that are very valuable from their genome preservation point of view. The aim of this study was the investigation of the genetic diversity of seven Montenegrin sheep breeds (Jezeropivska-JP, Sora-SOR, Zetska zuja-ZZ, Bardoka-BAR, Sjenicka-SJ, Ljaba-Lj, and Piperska zuja-PIP) using 18 microsatellite sets of markers. The genotyping was done for 291 samples from seven populations using the multiplex amplification of sequences with polymerase chain reaction (PCR). The parameters of genetic diversity were estimated using several software tools. In total, 243 alleles were found, with a range of 6 to 25 by locus. The mean observed heterozygosity (Ho), polymorphism information content (PIC), and Fis values (fixation index) per marker were 0.728, 0.781, and -0.007, respectively. The mean number of alleles per breed varied from 4.889 in ZZ to 10.056 in JP. The highest Ho was estimated for JP (0.763) and the lowest for ZZ (0.640). The genetic structure showed close relations between SOR and JP, and both of them with SJ, while ZZ, LJ, and PIP were more distanced. This study provides useful indicators for the development of further in-depth studies and the creation of appropriate conservation programs.

The Consequences of Mitochondrial T10432C Mutation in Cika Cattle: A "Potential" Model for Leber's Hereditary Optic Neuropathy.

While mitogenome mutations leading to pathological manifestations are rare, more than 200 such mutations have been described in humans. In contrast, pathogenic mitogenome mutations are rare in domestic animals and have not been described at all in cattle. In the small local Slovenian cattle breed Cika, we identified (next-generation sequencing) two cows with the T10432C mitogenome mutation in the ND4L gene, which corresponds to the human T10663C mutation known to cause Leber's hereditary optic neuropathy (LHON). Pedigree analysis revealed that the cows in which the mutation was identified belong to two different maternal lineages with 217 individual cows born between 1997 and 2020. The identified mutation and its maternal inheritance were confirmed by Sanger sequencing across multiple generations, whereas no single analysis revealed evidence of heteroplasmy. A closer clinical examination of one cow with the T10432C mutation revealed exophthalmos, whereas histopathological examination revealed retinal ablations, subretinal oedema, and haemorrhage. The results of these analyses confirm the presence of mitochondrial mutation T10432C with homoplasmic maternal inheritance as well as clinical and histopathological signs similar to LHON in humans. Live animals with the mutation could be used as a suitable animal model that can improve our understanding of the pathogenesis of LHON and other mitochondriopathies.

Genetic diversity and population structure of six autochthonous pig breeds from Croatia, Serbia, and Slovenia.

BACKGROUND: The importance of local breeds as genetic reservoirs of valuable genetic variation is well established. Pig breeding in Central and South-Eastern Europe has a long tradition that led to the formation of several local pig breeds. In the present study, genetic diversity parameters were analysed in six autochthonous pig breeds from Slovenia, Croatia and Serbia (Banija spotted, Black Slavonian, Turopolje pig, Swallow-bellied Mangalitsa, Moravka and Krskopolje pig). Animals from each of these breeds were genotyped using microsatellites and single nucleotide polymorphisms (SNPs). The results obtained with these two marker systems and those based on pedigree data were compared. In addition, we estimated inbreeding levels based on the distribution of runs of homozygosity (ROH) and identified genomic regions under selection pressure using ROH islands and the integrated haplotype score (iHS). RESULTS: The lowest heterozygosity values calculated from microsatellite and SNP data were observed in the Turopolje pig. The observed heterozygosity was higher than the expected heterozygosity in the Black Slavonian, Moravka and Turopolje pig. Both types of markers allowed us to distinguish clusters of individuals belonging to each breed. The analysis of admixture between breeds revealed potential gene flow between the Mangalitsa and Moravka, and between the Mangalitsa and Black Slavonian, but no introgression events were detected in the Banija spotted and Turopolje pig. The distribution of ROH across the genome was not uniform. Analysis of the ROH islands identified genomic regions with an extremely high frequency of shared ROH within the Swallow-bellied Mangalitsa, which harboured genes associated with cholesterol biosynthesis, fatty acid metabolism and daily weight gain. The iHS approach to detect signatures of selection revealed candidate regions containing genes with potential roles in reproduction traits and disease resistance. CONCLUSIONS: Based on the estimation of population parameters obtained from three data sets, we showed the existence of relationships among the six pig breeds analysed here. Analysis of the distribution of ROH allowed us to estimate the level of inbreeding and the extent of homozygous regions in these breeds. The iHS analysis revealed genomic regions potentially associated with phenotypic traits and allowed the detection of genomic regions under selection pressure.

Post-genotyping optimization of dataset formation could affect genetic diversity parameters: an example of analyses with alpine goat breeds.

BACKGROUND: Local breeds retained unique genetic variability important for adaptive potential especially in light of challenges related to climate change. Our first objective was to perform, for the first time, a genome-wide diversity characterization using Illumina GoatSNP50 BeadChip of autochthonous Dreznica goat breed from Slovenia, and five and one local breeds from neighboring Austria and Italy, respectively. For optimal conservation and breeding programs of endangered local breeds, it is important to detect past admixture events and strive for preservation of purebred representatives of each breed with low or without admixture. In the second objective, we hence investigated the effect of inclusion or exclusion of outliers from datasets on genetic diversity and population structure parameters. RESULTS: Distinct genetic origin of the Dreznica goat was demonstrated as having closest nodes to Austrian and Italian breeds. A phylogenetic study of these breeds with other goat breeds having SNP data available in the DRYAD repository positioned them in the alpine, European and global context. Swiss breeds clustered with cosmopolitan alpine breeds and were closer to French and Spanish breeds. On the other hand, the Dreznica goat, Austrian and Italian breeds were closer to Turkish breeds. Datasets where outliers were excluded affected estimates of genetic diversity parameters within the breed and increased the pairwise genetic distances between most of the breeds. Alpine breeds, including Dreznica, Austrian and Italian goats analyzed here, still exhibit relatively high levels of genetic variability, homogeneous genetic structure and strong geographical partitioning. CONCLUSIONS: Genetic diversity analyses revealed that the Slovenian Dreznica goat has a distinct genetic identity and is closely related to the neighboring Austrian and Italian alpine breeds. These results expand our knowledge on phylogeny of goat breeds from easternmost part of the European Alps. The here employed outlier test and datasets optimization approaches provided an objective and statistically powerful tool for removal of admixed outliers. Importance of this test in selecting the representatives of each breed is warranted to obtain more objective diversity parameters and phylogenetic analysis. Such parameters are often the basis of breeding and management programs and are therefore important for preserving genetic variability and uniqueness of local rare breeds.

Genomic Characterization of the Istrian Shorthaired Hound.

Istrian shorthaired hound is an old indigenous Croatian dog breed with historical traces of its origin, which date back to the 14th century. Due to its intelligence and great hunting abilities, it is considered an excellent hunting dog. Despite its ancient origin, there is no data on genetic diversity, population structure, and degree of inbreeding that could be used for advanced management and conservation of this breed. Our study aimed to provide a high-resolution population structure of the Istrian shorthaired hound using a 220K HD SNP array, to compare the obtained data with the genealogical records and to place the breed in a broader context of world dog populations. Relatively high population size and low inbreeding coefficient estimated from genealogical data indicate a preserved genetic diversity in this breed. The principle component analysis, the NeighborNet network, and TreeMix were used to determine the genetic relationship between the Istrian shorthaired hound and other breeds. The Istrian shorthaired hound was found to be genetically related to Italian hunting dogs sharing the same branch with the Segugio Italiano a Pelo Raso and Segugio Italiano a Pelo Forte. The ADMIXTURE analysis indicated that the Istrian shorthaired hound could be involved in the development of some other hunting dog breeds. The estimated effective population size (Ne) based on SNP data was similar to Ne calculated from genealogical data indicating the absence of bottlenecks and well-balanced use of breeding animals. The low genomic inbreeding coefficient, together with the higher number of short runs of homozygosity, observed in the Istrian shorthaired hound, confirms the ancient origin of the breed based on historical documents. The analysis of selective sweeps identified genomic regions with the strongest selection signals in the vicinity of the genes associated with cognitive performance and behavior. Genome analysis proved to be a useful tool for estimating population parameters and can be implemented in the conservation plan for this indigenous breed.

Exposure of honey bee larvae to thiamethoxam and its interaction with Nosema ceranae infection in adult honey bees.

During their lifetime honey bees (Apis mellifera) rarely experience optimal conditions. Sometimes, a simultaneous action of multiple stressors, natural and chemical, results in even greater effect than of any stressor alone. Therefore, integrative investigations of different factors affecting honey bees have to be carried out. In this study, adult honey bees exposed to thiamethoxam in larval and/or adult stage and infected with Nosema ceranae were examined. Newly emerged bees from colonies, non-treated or treated with thiamethoxam, were organized in six groups and kept in cages. Thiamethoxam treated bees were further exposed to either thiamethoxam or Nosema (groups TT and TN), or simultaneously to both (group TTN). Newly emerged bees from non-treated colonies were exposed to Nosema (group CN). From both, treated and non-treated colonies two groups were organized and further fed only with sugar solution (groups C and TC). Here, we present the expression profile of 19 genes in adult worker honey bees comprising those involved in immune, detoxification, development and apoptosis response. Results showed that gene expression patterns changed with time and depended on the treatment. In group TC at the time of emergence the majority of tested genes were downregulated, among which nine were significantly altered. The same gene pattern was observed on day six, where the only significantly upregulated gene was defensin-1. On day nine most of analyzed genes in all experimental groups showed upregulation compared to control group, where upregulation of antimicrobial peptide genes abaecin, defensin-1 and defensin-2 was significant in groups TT and TTN. On day 15 we observed a similar pattern of expression in groups TC and TT exposed to thiamethoxam only, where most of the detoxification genes were downregulated. Additionally RNA loads of Nosema and honey bee viruses were recorded. We detected a synergistic interaction of thiamethoxam and Nosema, reflected in lowest honey bee survival.

Population structure and genetic history of Tibetan Terriers.

BACKGROUND: Tibetan Terrier is a popular medium-sized companion dog breed. According to the history of the breed, the western population of Tibetan Terriers includes two lineages, Lamleh and Luneville. These two lineages derive from a small number of founder animals from the native Tibetan Terrier population, which were brought to Europe in the 1920s. For almost a century, the western population of Tibetan Terriers and the native population in Tibet were reproductively isolated. In this study, we analysed the structure of the western population of Tibetan Terriers, the original native population from Tibet and of different crosses between these two populations. We also examined the genetic relationships of Tibetan Terriers with other dog breeds, especially terriers and some Asian breeds, and the within-breed structure of both Tibetan Terrier populations. RESULTS: Our analyses were based on high-density single nucleotide polymorphism (SNP) array (Illumina HD Canine 170 K) and microsatellite (18 loci) genotypes of 64 Tibetan Terriers belonging to different populations and lineages. For the comparative analysis, we used 348 publicly available SNP array genotypes of dogs from other breeds. We found that the western population of Tibetan Terriers and the native Tibetan Terriers clustered together with other Asian dog breeds, whereas all other terrier breeds were grouped into a separate group. We were also able to differentiate the western Tibetan Terrier lineages (Lamleh and Luneville) from the native Tibetan Terrier population. CONCLUSIONS: Our results reveal the relationships between the western and native populations of Tibetan Terriers and support the hypothesis that Tibetan Terrier belongs to the group of ancient dog breeds of Asian origin, which are close to the ancestors of the modern dog that were involved in the early domestication process. Thus, we were able to reject the initial hypothesis that Tibetan Terriers belong to the group of terrier breeds. The existence of this native population of Tibetan Terriers at its original location represents an exceptional and valuable genetic resource.

TLR2 polymorphism (rs650082970) is associated with somatic cell count in goat milk.

Pathogens invading the mammary gland are recognized through a range of pattern recognition receptors (PRRs), residing on the plasma membrane of mammary epithelial cells. Toll-like receptor 2 (TLR2) signalling is responsible for recognition of Gram-positive bacteria, which are the most common mastitis-causing pathogens in goats. Somatic cell counts (SCC) in milk are routinely determined in goat dairy flocks and serve as an indicator of milk quality, which is highly correlated to intramammary infections. Recently, a single nucleotide polymorphism of the TLR2 was suggested to be associated with SCC in goat milk. To further test the suggested association, we genotyped 61 Slovenian Alpine goats included in the dataset. The effect of the genotype was analysed using the general linear model (GLM) procedure of SAS/STAT software. We found the TLR2 genotypes significantly (p = 0.0007) associated with milk SCC. Animals with the A/G genotype had significantly (p ≤ 0.05) lower SCC value in milk compared to the G/G genotype. Our data suggest that the A allele is the minor one and is associated with lower milk SCC. In the current study, we provide a validated PCR-RFLP based genotyping assay for the TLR2 SNP (rs650082970) and confirm its association with milk SCC. Further studies to confirm the association on a larger number of animals of different breeds and to explain functional consequences of the polymorphism in relation to SCC are encouraged.

Response of adult honey bees treated in larval stage with prochloraz to infection with Nosema ceranae.

Among numerous factors that contribute to honey bee colony losses and problems in beekeeping, pesticides and Nosema ceranae have been often reported. In contrast to insecticides, whose effects on bees have been widely studied, fungicides did not attract considerable attention. Prochloraz, an imidazole fungicide widely used in agriculture, was detected in honey and pollen stored inside hives and has been already proven to alter immune gene expression of honey bees at different developmental stages. The aim of this study was to simulate the realistic conditions of migratory beekeeping, where colonies, both uninfected and infected with N. ceranae, are frequently transported to the vicinity of crop fields treated with prochloraz. We investigated the combined effect of prochloraz and N. ceranae on honey bees that faced fungicide during the larval stage through food consumption and microsporidium infection afterwards. The most pronounced changes in gene expression were observed in newly emerged Nosema-free bees originating from colonies previously contaminated with prochloraz. As exclusively upregulation was registered, prochloraz alone most likely acts as a challenge that induces activation of immune pathways in newly emerged bees. The combination of both stressors (prochloraz and Nosema infection) exerted the greatest effect on six-day-old honey bees. Among ten genes with significantly altered expression, half were upregulated and half downregulated. N. ceranae as a sole stressor had the weakest effects on immune gene expression modulation with only three genes significantly dysregulated. In conclusion, food contaminated with prochloraz consumed in larval stage could present a threat to the development of immunity and detoxification mechanisms in honey bees.

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.

In silico screening of the chicken genome for overlaps between genomic regions: microRNA genes, coding and non-coding transcriptional units, QTL, and genetic variations.

MicroRNAs (miRNAs) are a class of non-coding RNAs involved in posttranscriptional regulation of target genes. Regulation requires complementarity between target mRNA and the mature miRNA seed region, responsible for their recognition and binding. It has been estimated that each miRNA targets approximately 200 genes, and genetic variability of miRNA genes has been reported to affect phenotypic variability and disease susceptibility in humans, livestock species, and model organisms. Polymorphisms in miRNA genes could therefore represent biomarkers for phenotypic traits in livestock animals. In our previous study, we collected polymorphisms within miRNA genes in chicken. In the present study, we identified miRNA-related genomic overlaps to prioritize genomic regions of interest for further functional studies and biomarker discovery. Overlapping genomic regions in chicken were analyzed using the following bioinformatics tools and databases: miRNA SNiPer, Ensembl, miRBase, NCBI Blast, and QTLdb. Out of 740 known pre-miRNA genes, 263 (35.5 %) contain polymorphisms; among them, 35 contain more than three polymorphisms The most polymorphic miRNA genes in chicken are gga-miR-6662, containing 23 single nucleotide polymorphisms (SNPs) within the pre-miRNA region, including five consecutive SNPs, and gga-miR-6688, containing ten polymorphisms including three consecutive polymorphisms. Several miRNA-related genomic hotspots have been revealed in chicken genome; polymorphic miRNA genes are located within protein-coding and/or non-coding transcription units and quantitative trait loci (QTL) associated with production traits. The present study includes the first description of an exonic miRNA in a chicken genome, an overlap between the miRNA gene and the exon of the protein-coding gene (gga-miR-6578/HADHB), and the first report of a missense polymorphism located within a mature miRNA seed region. Identified miRNA-related genomic hotspots in chicken can serve researchers as a starting point for further functional studies and association studies with poultry production and health traits and the basis for systematic screening of exonic miRNAs and missense/miRNA seed polymorphisms in other genomes.

Development of Integrative Map of MicroRNA Gene Regulatory Elements.

Like many other omics types, miRNomics is the topic of increased research interest. Fast accumulation of diverse biological information generates fragmented and sometimes contradictory results. There is an urgent need for development of protocols for systematic assembling, organization and integration of accumulated genomics data. In the present study we developed a protocol for connecting microRNA (miRNA) interactions into an integrative map of miRNA regulatory elements including: upstream regulators, overlapping genomic elements and downstream targets; miRNA-target interactions (MTIs). Data integration (integromics) in miRNomics enables revealing miRNA regulatory networks, obtaining new knowledge, setting more targeted hypotheses, prioritization of candidate regions, and facilitates biomarker development.

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.

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.

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.