Knockout of TRDMT1 methyltransferase affects DNA methylome in glioblastoma cells.

PURPOSE: We have previously shown that TRDMT1 methyltransferase is a regulator of chemotherapy-associated responses in glioblastoma cells. Despite the fact that glioblastoma, a common and malignant brain tumor, is widely characterized in terms of genetic and epigenetic markers, there are no data on TRDMT1-related changes in 5-methylcytosine pools in the genome. In the present study, the effect of TRDMT1 gene knockout (KO) on DNA methylome was analyzed. METHODS: CRISPR-based approach was used to obtain TRDMT1 KO glioblastoma cells. Total 5-methylcytosine levels in DNA, DNMT1 pools and DNMT activity were studied using ELISA. Reduced representation bisulfite sequencing (RRBS) was considered to comprehensively evaluate DNA methylome in glioblastoma cells with TRDMT1 KO. RESULTS: TRDMT1 KO cells were characterized by decreased levels of total 5-methylcytosine in DNA and DNMT1, and DNMT activity. RRBS-based methylome analysis revealed statistically significant differences in methylation-relevant DMS-linked genes in control cells compared to TRDMT1 KO cells. TRDMT1 KO-associated changes in DNA methylome may affect the activity of several processes and pathways such as telomere maintenance, cell cycle and longevity regulating pathway, proteostasis, DNA and RNA biology. CONCLUSIONS: TRDMT1 may be suggested as a novel modulator of gene expression by changes in DNA methylome that may affect cancer cell fates during chemotherapy. We postulate that the levels and mutation status of TRDMT1 should be considered as a prognostic marker and carefully monitored during glioblastoma progression.

Trichostatin A-Mediated Epigenetic Modulation Predominantly Triggers Transcriptomic Alterations in the Ex Vivo Expanded Equine Chondrocytes.

Epigenetic mechanisms of gene regulation are important for the proper differentiation of cells used for therapeutic and regenerative purposes. The primary goal of the present study was to investigate the impacts of 5-aza-2' deoxycytidine (5-AZA-dc)- and/or trichostatin A (TSA)-mediated approaches applied to epigenomically modulate the ex vivo expanded equine chondrocytes maintained in monolayer culture on the status of chondrogenic cytodifferentiation at the transcriptome level. The results of next-generation sequencing of 3' mRNA-seq libraries on stimulated and unstimulated chondrocytes of the third passage showed no significant influence of 5-AZA-dc treatment. Chondrocytes stimulated with TSA or with a combination of 5-AZA-dc+TSA revealed significant expressional decline, mainly for genes encoding histone and DNA methyltransferases, but also for other genes, many of which are enriched in canonical pathways that are important for chondrocyte biology. The TSA- or 5-AZA-dc+TSA-induced upregulation of expanded chondrocytes included genes that are involved in histone hyperacetylation and also genes relevant to rheumatoid arthritis and inflammation. Chondrocyte stimulation experiments including a TSA modifier also led to the unexpected expression incrementation of genes encoding HDAC3, SIRT2, and SIRT5 histone deacetylases and the MBD1 CpG-binding domain protein, pointing to another function of the TSA agent besides its epigenetic-like properties. Based on the transcriptomic data, TSA stimulation seems to be undesirable for chondrogenic differentiation of passaged cartilaginous cells in a monolayer culture. Nonetheless, obtained transcriptomic results of TSA-dependent epigenomic modification of the ex vivo expanded equine chondrocytes provide a new source of data important for the potential application of epigenetically altered cells for transplantation purposes in tissue engineering of the equine skeletal system.

Expression of cytokines in dairy cattle mammary gland parenchyma during chronic staphylococcal infection.

The study aim was to determine the expression of genes potentially related to chronic mastitis at the mRNA and protein levels, viz. chemokine C-C motif receptor 1 (CCR1), C-C motif chemokine ligand 2 (CCL2), C-C motif chemokine ligand 5 (CXCL5), tumor necrosis factor α (TNFα), interleukin 1ß (IL-1ß), interleukin 6 (IL-6), interleukin 8 (IL-8), interleukin 18 (IL-18), in bovine mammary gland parenchyma. The study examines the differences in expression of selected genes between cows with chronic mastitis caused by coagulase-positive (CoPS) or coagulase-negative staphylococci (CoNS) and those with healthy udders (H). Samples were collected from the udder quarters from 40 Polish Holstein-Friesian cows; 54 of these samples were chosen for analysis based on microbiological analysis of milk taken two days before slaughter. They were categorized into three groups: CoPS (N = 27), CoNS (N = 14) and H (N = 13). The RNA expression was analyzed by RT-qPCR and protein concentration by ELISA. No differences in the mRNA levels of seven genes (TNFα, IL-18, CCR1, IL-1ß, CCL2, IL-8, IL-6) and four proteins (TNFα, IL-18, CCR1, IL-1ß) were identified between the CoPS and H groups. Higher transcript levels of CXCL5 (p ≤ 0.05) gene were noted in CoPS than in H. Compared to H, higher concentrations of IL-8 and CXCL5 (p ≤ 0.05) were observed in CoPS (0.05 < p < 0.1) and CCL2 (0.05 < p < 0.1) in CoNS, while lower levels of Il-6 were found in CoPS. This may suggest that during chronic mastitis the organism stops producing pro-inflammatory cytokines, probably to protect the host tissues against their damage during prolonged infection.

MicroRNA expression profile in bovine mammary gland parenchyma infected by coagulase-positive or coagulase-negative staphylococci.

MicroRNAs (miRNAs) are short, non-coding RNAs, 21-23 nucleotides in length which are known to regulate biological processes that greatly impact immune system activity. The aim of the study was to compare the miRNA expression in non-infected (H) mammary gland parenchyma samples with that of glands infected with coagulase-positive staphylococci (CoPS) or coagulase-negative staphylococci (CoNS) using next-generation sequencing. The miRNA profile of the parenchyma was found to change during mastitis, with its profile depending on the type of pathogen. Comparing the CoPS and H groups, 256 known and 260 potentially new miRNAs were identified, including 32 that were differentially expressed (p ≤ 0.05), of which 27 were upregulated and 5 downregulated. Comparing the CoNS and H groups, 242 known and 171 new unique miRNAs were identified: 10 were upregulated (p ≤ 0.05), and 2 downregulated (p ≤ 0.05). In addition, comparing CoPS with H and CoNS with H, 5 Kyoto Encyclopedia of Genes and Genomes pathways were identified; in both comparisons, differentially-expressed miRNAs were associated with the bacterial invasion of epithelial cells and focal adhesion pathways. Four gene ontology terms were identified in each comparison, with 2 being common to both immune system processes and signal transduction. Our results indicate that miRNAs, especially miR-99 and miR-182, play an essential role in the epigenetic regulation of a range of cellular processes, including immunological systems bacterial growth in dendritic cells and disease pathogenesis (miR-99), DNA repair and tumor progression (miR-182).

Transcripts and protein levels of CSN1S1 and CSN3 genes in dairy cattle mammary gland secretory tissue during chronic staphylococcal infection.

Our objective was to determine the influence of chronic coagulase-positive staphylococci (CoPS) or coagulase-negative staphylococci (CoNS) infection on the mRNA and protein levels of two main milk proteins responsible for cheese curd quantity and quality, alpha-S1-casein (CSN1S1) and kappa-casein (CSN3). Measurements were made in cow mammary parenchyma with a prevalence of secretory tissue (MGST). Samples of MGST were collected from the separate quarters and divided into CoPS, CoNS and bacteria-free (H) groups according to the microbiological status of the quarter milk. No differences in CSN1S1 and CSN3 mRNA level were found between groups, however, CSN1S1 protein level was significantly higher in the H group than the CoNS group, and CSN3 protein level was significantly higher in H than CoPS group. Hence, while the CSN1S1 and CSN3 genes appear to be constitutively expressed at the mRNA level in dairy cow MGST during mastitis, CoNS infection negatively affected CSN1S1 protein level, and CoPS infection negatively affected CSN3 protein level. The lack of change at the mRNA level suggests that staphylococcal infection may affect the post-transcriptional or post-translational modifications.

Corn dried distillers grains with solubles (cDDGS) in the diet of pigs change the expression of adipose genes that are potential therapeutic targets in metabolic and cardiovascular diseases.

BACKGROUND: Corn dried distillers grains with solubles (cDDGS) are a byproduct of biofuel and alcohol production. cDDGS have been used in pig feed for many years, because they are readily available and rich in protein, fiber, unsaturated fatty acids and phytosterols. However, feed mixtures too high in cDDGS result in the worsening of backfat quality. We performed RNA-sequencing analysis of backfat from crossbred pigs fed different diets. The diets were isoenergetic but contained different amounts of cDDGS and various sources of fats. The animals were divided into four dietary groups during the two months of experimentation: group I (control (-cDDGS+rapeseed oil)), group II (+cDDGS+rapeseed oil), group III (+cDDGS+beef tallow), and group IV (+cDDGS+coconut oil). The aim of the present experiment was to evaluate changes in the backfat transcriptome of pigs fed isoenergetic diets that differed in cDDGS presence. RESULTS: Via DESeq2 software, we identified 93 differentially expressed genes (DEGs) between groups I and II, 13 between groups I and III, and 125 between groups I and IV. DEGs identified between group I (-cDDGS+rapeseed oil) and group II (+cDDGS+rapeseed oil) were highly overrepresented in several KEGG pathways: metabolic pathways (FDR < 1.21e-06), oxidative phosphorylation (FDR < 0.00189), fatty acid biosynthesis (FDR < 0.00577), Huntington's disease (FDR < 0.00577), fatty acid metabolism (FDR < 0.0112), Parkinson's disease (FDR < 0.0151), non-alcoholic fatty liver disease (NAFLD) (FDR < 0.016), Alzheimer's disease (FDR < 0.0211) and complement and coagulation cascades (FDR < 0.02). CONCLUSIONS: We observed that the addition of cDDGS positively affects the expression of several genes that have been recently proposed as potential targets for the treatment of obesity, diabetes, cardiovascular disease, and Alzheimer's disease (e.g., FASN, AACS, ALAS1, HMGCS1, and VSIG4). Thus, our results support the idea of including cDDGS into the diets of companion animals and humans and encourage research into the bioactive ingredients of cDDGS.

Mutations in Leptin (LEP) Gene Are Associated with Carcass and Meat Quality Traits in Crossbreed Rabbits.

Leptin is a hormone synthesized and secreted primarily in adipose cells that help to regulate energy balance. This study examined the associations of single nucleotide polymorphisms in the rabbit leptin gene with growth traits, slaughter traits and physicochemical parameters of New Zealand White (NZW) and Belgian Giant Grey (BGG) crossbreed rabbits. In total, 320 crossbreed animals were genotyped for polymorphisms within exon 2-g.16081633T>C, intron 1_2-g.16081420C>T, and within UTR-g.16079636C>G for association analysis. Identified polymorphisms within rabbits leptin gene showed significant differences for dissectible fat percentage in carcass and dissectible fat weight in intermediate part (g.16081633T>C). Moreover, meat traits like protein content (g.16081633T>C; g.16079636C>G), intramuscular fat content (g.16081633T>C; g.16079636C>G, g.16081420C>T), dry matter (g.16081420C>T), ash (g.16081420C>T), water (g.16081420C>T), and cohesiveness (g.16081420C>T, g.16079636C>G) were affected by polymorphisms in leptin gene. We conclude that polymorphism in the rabbit leptin gene influences important carcass and meat traits of NZW × BGG crossbreeds. Therefore, polymorphisms identified in this study may be used in selection as a meat trait markers.

The use of SNP markers for cattle breed identification.

A potential application of single nucleotide polymorphisms (SNPs) in animal husbandry and production is identification of the animal breed. In this study, using chosen marker selection methods and genotypic data obtained with the use of Illumina Bovine SNP50 BeadChip for individuals belonging to ten cattle breeds, the reduced panels containing the most informative SNP markers were developed. The suitability of selected SNP panels for the effective and reliable assignment of the studied individuals to the breed of origin was checked by three allocation algorithms implemented in GeneClass 2. The studied breeds set included both Polish-native breeds under the genetic resources conservation programs and highly productive breeds with a global range. For all of the tested marker selection methods ("delta" and two FST-based variants), two separate methodological approaches of marker assortment were used and three marker panels were created with 96, 192, and 288 SNPs respectively, to determine the minimum number of markers required for effective differentiation of the studied breeds. Moreover, the usefulness of the most effective panels of markers to assess the population structure and genetic diversity of the analyzed breeds was examined. The conducted analyses showed the possibility of using SNP subsets from medium-density genotypic microarrays to distinguish breeds of cattle kept in Poland and to analyze their genetic structure.

The evaluation of the usefulness of pedigree verification-dedicated SNPs for breed assignment in three polish cattle populations.

The breed assignment in cattle is one of the issues of molecular genetics which needs further testing and development. Although several statistical approaches have been developed to enable such application, the obtained results strongly depend on specific populations differentiation and power of markers discrimination or their informativeness. Currently, all breeding animals are being tested for parentage with the use of panel of 12 microsatellite markers, which in near future probably will be replaced by about 100 single nucleotide polymorphisms (SNPs). Despite the fact that SNPs are mainly bi-allelic, the multilocus genotypes can reach the level of polymorphism of a panel of microsatellite markers. In this study we attempted to determine the breed of origin of 741 cattle by using 120 SNPs dedicated for parentage testing and included in the BovineSNP50 BeadChip genotyping assay (Illumina). The applied Bayesian and frequency-based methods allowed such differentiation, however, the reliability of the results was not completely satisfying, suggesting that the studied markers are not the best tool for breed assignment.

The evaluation of bovine SNP50 beadchip assay performance in Polish red cattle breed.

According to increasing interest in the use of high density SNP (single nucleotide polymorphism) genotyping assays for genome-wide genetic studies in farm animals, there is a need to assess the usefulness of currently available genomic tools for application in different breeds. The performance of the assays may differ between the breeds because of discrepancies in allele frequencies of the polymorphisms or differences in linkage disequilibrium patterns. In this study we attempted to test the performance of the Bovine SNP50 v2 genotyping assay (Illumina) for population genetic and other applications in the Polish Red cattle breed. We found that 37,977 of the 53,438 autosomal markers included in the assay give high quality genomic information and can be used for different applications concerning this breed. The remaining markers were denoted as "of limited use" or redundant because of their weak performance, deviation from Hardy-Weinberg equilibrium, low minor allele frequency, high linkage disequilibrium with other neighboring markers or location on sex chromosomes.

Methylome and transcriptome data integration reveals aberrantly regulated genes in equine sarcoids.

DNA methylation is a key mechanism in transcription regulation, and aberrant methylation is a common and important mechanism in tumor initiation, maintenance, and progression. To find genes that are aberrantly regulated by altered methylation in horse sarcoids, we used reduced representation bisulfite sequencing (RRBS) accompanied by RNA sequencing (RNA-Seq) for methylome (whole genome DNA methylation sequencing) and transcriptome profiling, respectively. We found that the DNA methylation level was generally lower in lesion samples than in controls. In the analyzed samples, a total of 14,692 differentially methylated sites (DMSs) in the context of CpG (where cytosine and guanine are separated by a phosphate), and 11,712 differentially expressed genes (DEGs) were identified. The integration of the methylome and transcriptome data suggests that aberrant DNA methylation may be involved in the deregulation of expression of the 493 genes in equine sarcoid. Furthermore, enrichment analysis of the genes demonstrated the activation of multiple molecular pathways related to extracellular matrix (ECM), oxidative phosphorylation (OXPHOS), immune response, and disease processes that can be related to tumor progression. The results provide further insight into the epigenetic alterations in equine sarcoids and provide a valuable resource for follow-up studies to identify biomarkers for predicting susceptibility to this common condition in horses.

Methylation Status of Gene Bodies of Selected microRNA Genes Associated with Neoplastic Transformation in Equine Sarcoids.

Horses are of great importance in recreation, livestock production, as working animals in poorly developed countries, and for equine-assisted therapy. Equine sarcoids belong to the most commonly diagnosed tumors in this species. They may cause discomfort, pain, and can lead to the permanent impairment of motor function. The molecular bases of their formation are still under investigation. Our previous studies revealed altered microRNA (miRNA) expression and DNA methylation levels in sarcoid tumors. Abnormal patterns of methylation may be responsible for changes in gene expression levels, including microRNAs. Recently, the DNA methylation of gene bodies has also been shown to have an impact on gene expression. Thus, the aim of the study was to investigate the methylation pattern of gene bodies of chosen miRNAs identified in sarcoid tissue (miR-101, miR-10b, miR-200a, and miR-338-3p), which have also been established to play roles in neoplastic transformation. To this end, we applied qRT-PCR, Bisulfite Sequencing PCR (BSP), and Mquant methods. As a result, we identified the statistically significant downregulation of pri-mir-101-1, pri-mir-10b, and pri-mir-200a in the sarcoid samples in comparison to the control. The DNA methylation analysis revealed their hypermethylation. This suggests that DNA methylation may be one mechanism responsible for the downregulation of theses miRNAs. However, the identified differences in the methylation levels are not very high, which implies that other mechanisms may also underlie the downregulation of the expression of these miRNAs in equine sarcoids. For the first time, the results obtained shed light on microRNA expression regulation by gene body methylation in equine sarcoids and provide bases for further deeper studies on other mechanisms influencing the miRNA repertoire.

RNA 5-methylcytosine status is associated with DNMT2/TRDMT1 nuclear localization in osteosarcoma cell lines.

Osteosarcoma (OS) is a pediatric malignant bone tumor with unsatisfying improvements in survival rates due to limited understanding of OS biology and potentially druggable targets. The present study aims to better characterize osteosarcoma U-2 OS, SaOS-2, and MG-63 cell lines that are commonly used as in vitro models of OS. We focused on evaluating the differences in cell death pathways, redox equilibrium, the activity of proliferation-related signaling pathways, DNA damage response, telomere maintenance, DNMT2/TRDMT1-based responses and RNA 5-methylcytosine status. SaOS-2 cells were characterized by higher levels of superoxide and nitric oxide that promoted AKT and ERK1/2 activation thus modulating cell death pathways. OS cell lines also differed in the levels and localization of DNA repair regulator DNMT2/TRDMT1. SaOS-2 cells possessed the lowest levels of total, cytoplasmic and nuclear DNMT2/TRDMT1, whereas in MG-63 cells, the highest levels of nuclear DNMT2/TRDMT1 were associated with the most pronounced status of RNA 5-methylcytosine. In silico analysis revealed potential phosphorylation sites at DNMT2/TRDMT1 that may be related to the regulation of DNMT2/TRDMT1 localization. We postulate that redox homeostasis, proliferation-related pathways and DNMT2/TRDMT1-based effects can be modulated as a part of anti-osteosarcoma strategy reflecting diverse phenotypic features of OS cells.

Another lesson from unmapped reads: in-depth analysis of RNA-Seq reads from various horse tissues.

In recent years, a vast amount of sequencing data has been generated and large improvements have been made to reference genome sequences. Despite these advances, significant portions of reads still do not map to reference genomes and these reads have been considered as junk or artificial sequences. Recent studies have shown that these reads can be useful, e.g., for refining reference genomes or detecting contaminating microorganisms present in the analyzed biological samples. A special case of this is RNA sequencing (RNA-Seq) reads that come from tissue transcriptomes. Unmapped reads from RNA-Seq have received much less attention than those from whole-genome sequencing. In particular, in the horse, an analysis of unmapped RNA reads has not been performed yet. Thus, in this study, we analyzed the unmapped reads originating from the RNA-Seq performed through the Functional Annotation of Animal Genomes (FAANG) project in the horse, using eight different tissues from two mares. We demonstrated that unmapped reads from RNA-Seq could be easily assembled into transcripts relating to many important genes present in the sequences of other mammals. Large portions of these transcripts did not have coding potential and, thus, can be considered as non-coding RNA. Moreover, reads that were not mapped to the reference genome but aligned to the entries in NCBI database of horse proteins were enriched for biological processes that largely correspond to the functions of organ from which RNA was isolated and thus are presumably true transcripts of genes associated with cell metabolism in those tissues. In addition, a portion of reads aligned to the common pathogenic or neutral microbiota, of which the most common was Brucella spp. These data suggest that unmapped reads can be an important target for in-depth analysis that may substantially enrich results of initial RNA-Seq experiments for various tissues and organs.

Genetic Variability and Population Structure of Polish Konik Horse Maternal Lines Based on Microsatellite Markers.

The aim of the conservation programme is to maintain the population size of endangered livestock breeds of less economic importance at a level that ensures the survival of the breed, the preservation of genetic diversity, and the preservation of as many pedigree lines as possible. The Polish Konik, a native Polish primitive-type horse breed and is one of the breeds included in such a programme in Poland. Presently, there are only 16 (of the 35 maternal lines known in 1962), some of which are endangered. We examined the genetic variability and structure of the Polish Konik maternal lines (176 individuals) on the basis of the pedigree data and 17 microsatellite markers (STRs) from parentage testing. The overall mean number of alleles was 7.647 (±0.411), the effective number of alleles was 3.935 (±0.271), the mean number of alleles for which the frequency was equal to or lower than 5% was 4.471 (±0.286), and the mean information index was 1.506 (±0.087). The structure of the population and admixture patterns were calculated with the Structure and Structure Harvester software. The structural analysis indicated three likely genetic clusters; as the most optimal K value was estimated as 3, with ∆K of 15.4188. The F-statistics results indicated a low level of inbreeding (average inbreeding coefficient FIT was 0.0188, coefficient of differentiation FST was 0.0304, and mean inbreeding index value FIS was -0.0119). Variability monitoring should be carried out in order to avoid inbreeding depression, while breeding strategies should be designed to prevent the decrease of genetic variability in the Polish horse breed and to sustain the active female lines.

The Identification of a Novel Fucosidosis-Associated FUCA1 Mutation: A Case of a 5-Year-Old Polish Girl with Two Additional Rare Chromosomal Aberrations and Affected DNA Methylation Patterns.

Fucosidosis is a rare neurodegenerative autosomal recessive disorder, which manifests as progressive neurological and psychomotor deterioration, growth retardation, skin and skeletal abnormalities, intellectual disability and coarsening of facial features. It is caused by biallelic mutations in FUCA1 encoding the α-L-fucosidase enzyme, which in turn is responsible for degradation of fucose-containing glycoproteins and glycolipids. FUCA1 mutations lead to severe reduction or even loss of α-L-fucosidase enzyme activity. This results in incomplete breakdown of fucose-containing compounds leading to their deposition in different tissues and, consequently, disease progression. To date, 36 pathogenic variants in FUCA1 associated with fucosidosis have been documented. Among these are three splice site variants. Here, we report a novel fucosidosis-related 9-base-pair deletion (NG_013346.1:g.10233_10241delACAGGTAAG) affecting the exon 3/intron 3 junction within a FUCA1 sequence. This novel pathogenic variant was identified in a five-year-old Polish girl with a well-defined pattern of fucosidosis symptoms. Since it is postulated that other genetic, nongenetic or environmental factors can also contribute to fucosidosis pathogenesis, we performed further analysis and found two rare de novo chromosomal aberrations in the girl's genome involving a 15q11.1-11.2 microdeletion and an Xq22.2 gain. These abnormalities were associated with genome-wide changes in DNA methylation status in the epigenome of blood cells.

FISH mapping of six genes responsible for development of the nervous and skeletal systems on donkey (Equus asinus) chromosomes.

The results obtained in the present study made it possible to place selected markers responsible for development of the nervous and skeletal systems on the physical map of the donkey genome. Fluorescence in situ hybridization (FISH) was used to localize genes such as GDF5 (15q13), FRZB (4q23.1), TWIST (1q31), PAX6 (20q25), SALL1 (24q15) and SHH (1q35) on donkey chromosomes. The identification of their localization confirmed previously proposed homologies using ZOO-FISH technique, except for FRZB and SALL1 genes. This suggests that they were affected by rearrangements that changed their localization compared to horse, and in the case of the SALL1 gene also compared to human.

PRINS detection of 18S rDNA in pig, red fox and Chinese raccoon dog, and centromere DNA in horse.

The fluorescence in situ hybridization (FISH) technique is widely used in animal cytogenetics. Contrary to FISH procedure, primed in situ DNA synthesis (PRINS) does not require the DNA probe preparation (design, synthesis, gel purification of PCR products and labeling). The PRINS method with primers used as 'DNA probes' is both PCR-sensitive and allows for chromosomal localization of DNA sequences. Here, we show the application of PRINS reaction with one unlabeled oligonucleotide pair to identify 18S rDNA loci in three different animal species: domestic pig (Sus scrofa), red fox (Vulpes vulpes) and Chinese raccoon dog (Nyctereutes procyonoides procyonoides). We present the data of indirect labeling with the digoxigenin-PRINS using two different pairs of primers complementary to centromeric region of horse (Equus caballus) chromosomes. Our new PRINS application may be considered as a useful tool for chromosome investigation in the field of domestic and wild animal genetics and evolution.

Comparative Analysis of CpG Sites and Islands Distributed in Mitochondrial DNA of Model Organisms.

: The information about mtDNA methylation is still limited, thus epigenetic modification remains unclear. The lack of comprehensive information on the comparative epigenomics of mtDNA prompts comprehensive investigations of the epigenomic modification of mtDNA in different species. This is the first study in which the theoretical CpG localization in the mtDNA reference sequences from various species (12) was compared. The aim of the study was to determine the localization of CpG sites and islands in mtDNA of model organisms and to compare their distribution. The results are suitable for further investigations of mtDNA methylation. The analysis involved both strands of mtDNA sequences of animal model organisms representing different taxonomic groups of invertebrates and vertebrates. For each sequence, such parameters as the number, length, and localization of CpG islands were determined with the use of EMBOSS (European Molecular Biology Open Software Suite) software. The number of CpG sites for each sequence was indicated using the newcpgseek algorithm. The results showed that methylation of mtDNA in the analysed species involved mitochondrial gene expression. Our analyses showed that the CpG sites were commonly present in genomic regions including the D-loop, CYTB, ND6, ND5, ND4, ND3, ND2, ND1, COX3, COX2, COX1, ATP6, 16s rRNA, and 12s rRNA. The CpG distribution in animals from different species was diversified. Generally, the number of observed CpG sites of the mitochondrial genome was higher in the vertebrates than in the invertebrates. However, there was no relationship between the frequency of the CpG sites in the mitochondrial genome and the complexity of the analysed organisms. Interestingly, the distribution of the CpG sites for tRNA coding genes was usually cumulated in a larger CpG region in vertebrates. This paper may be a starting point for further research, since the collected information indicates possible methylation regions localized in mtDNA among different species including invertebrates and vertebrates.

Genetic Diversity and Population Structure of Polish Konik Horse Based on Individuals from All the Male Founder Lines and Microsatellite Markers.

The Polish Konik horse is a primitive native breed included in the genetic resource conservation program in Poland. After World War II, intensive breeding work began, aimed at rebuilding this breed. Now, the whole Polish Konik population is represented by six male founder lines (Wicek, Myszak, Glejt I, Goraj, Chochlik and Liliput). Individuals representing all six paternal lineages were selected based on their breeding documentation. We performed a fragment analysis with 17 microsatellite markers (STRs) recommended by the International Society for Animal Genetics (ISAG). The genetic diversity and structure within the paternal lineages and the whole of the studied group were investigated. The average allelic richness was 6.497 for the whole studied group. The fixation index (FST; measure of population differentiation) was low (about 3%), the mean inbreeding coefficient (FIT) was low and close to 0, and the mean inbreeding index value (FIS) was negative. The mean expected heterozygosity was established at 0.7046 and was lower than the observed heterozygosity. The power of discrimination and power of exclusion were 99.9999%. The cumulative parentage exclusion probability equaled 99.9269% when one parental genotype was known and 99.9996% with both parents' genotypic information was available. About 3% of the genetic variation was caused by differences in the breed origin and about 97% was attributed to differences among individuals. Our analysis revealed that there has been no inbreeding in the Polish Konik breed for the studied population. The genetic diversity was high, and its parameters were similar to those calculated for native breeds from other countries reported in the literature. However, due to the small number of breed founders and paternal lineages with unknown representation, the population's genetic diversity and structure should be monitored regularly.