Genetic variants of the hypoxia-inducible factor 3 alpha subunit (Hif3a) gene in the Fat and Lean mouse selection lines.

BACKGROUND: Adipose tissue hypoxia and members of the hypoxia-inducible factor alpha (HIFA) are involved in development of obesity. However, the mechanism and functions of HIF3A, one of three HIFA paralogs, in fat deposition have not been sufficiently studied. METHODS AND RESULTS: In the present study, we investigated whether Hif3a sequence variants are associated with divergent fat deposition in mouse selection lines for fatness and leanness. Sequencing and RFLP were used to analyse sequence variants within Hif3a. To identify candidate regulatory variants, we performed literature screening and used databases and bioinformatics tools like Ensembl, MethPrimer, TargetScanMouse, miRDB, PolyAsite, RISE, LncRRIsearch, RNAfold, PredictProtein, CAIcal, and switches.ELM Resource. There are 90 sequence variants in Hif3a between the two mouse lines. While most Fat line variants locate within intronic regions, Lean line variants are mainly in 3' UTR. We constructed a map of Hif3a potential regulatory regions and identified 39 regulatory variants by integrating data on constrained and regulatory elements, CpGs, and miRNAs and lncRNAs binding sites. Moreover, 3' UTR and two exonic variants may influence mRNA stability, translation rate and protein functionality. We propose as priority candidates for further functional studies a missense (rs37398126) and synonymous (rs37739792) variants, and intronic (rs47471302) variant that overlap conserved element in promoter region and predicted lncRNAs binding site. CONCLUSION: The results indicate a potential involvement of Hif3a in fat deposition. Additionally, approach used in the present study may serve as a general guideline for constructing an integrative gene map for prioritizing candidate gene variants with phenotypic effects.

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.

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.

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.

Establishment and characterization of proliferating primary cultures of equine epidermal keratinocytes.

Skin-derived tissue cultures are a useful model to study molecular mechanisms of skin renewal and pathogenesis of dermal diseases. Horses often suffer from skin diseases, skin trauma and problems with proper wound healing, which could be improved by in vitro grown keratinocyte grafts. Herein we describe establishment and characterization of equine skin-derived primary cell cultures, using enzymatic and explant methods. The established cell lines of primary equine keratinocytes (peK) maintained high proliferative capacity for over five passages and expressed different epithelial/keratinocyte-specific markers. Characterization of the primary culture was performed in parallel with localization studies of the markers in the skin histological sections, using commercially available antibodies. Relative expression of typical differentiation stage-specific markers was determined in the established cell lines, using RT-qPCR. Basal (proliferating) keratinocytes were the predominant cell type in the established cell lines, but low expression of post-mitotic keratinocyte markers was also detected. Differences in marker expression were observed neither between the peK originating from two different animals nor between the peK established with two different methods (enzymatically or by explanting). The described methods in combination with the suggested characterization and differentiation markers are suitable for establishment of proliferating peK and evaluation of their differentiation status.

Correction to: Genome-wide SNP data unveils the globalization of domesticated pigs.

An amendment to this paper has been published and can be accessed via the original article.

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.

Genetic diversity of autochthonous pig breeds analyzed by microsatellite markers and mitochondrial DNA D-loop sequence polymorphism.

The evaluation of the genetic structure of autochthonous pig breeds is very important for conservation of local pig breeds and preservation of diversity. In this study, 18 microsatellite loci were used to detect genetic relationship between autochthonous pig breeds [Black Slavonian (BS), Turopolje pig (TP), and Croatian wild boar] and to determine phylogenetic relationship among Croatian autochthonous pig breeds and certain Asian and European pigs using the mitochondrial DNA (mtDNA) D-loop sequence polymorphism. Relatively high degree of genetic variation was found between the observed populations. The analysis of mtDNA showed that haplotypes of the studied pig populations are different from the other European and Chinese haplotypes. BS pigs showed some similarities with Mangalitsa and Duroc breeds. The genetic distances of TP can be explained by high degree of inbreeding during the past century. Despite the European origin of Croatian pig breeds with some impact of Chinese breeds in the past, the results of present study show that genetic diversity is still pronounced within investigated breeds. Furthermore, the genetic diversity is even more pronounced between Croatian breeds and other European and Chinese pig breeds. Thus, conservation of Croatian pig breeds will contribute to overall genetic diversity preservation of pig breeds.

Genome-wide SNP data unveils the globalization of domesticated pigs.

BACKGROUND: Pigs were domesticated independently in Eastern and Western Eurasia early during the agricultural revolution, and have since been transported and traded across the globe. Here, we present a worldwide survey on 60K genome-wide single nucleotide polymorphism (SNP) data for 2093 pigs, including 1839 domestic pigs representing 122 local and commercial breeds, 215 wild boars, and 39 out-group suids, from Asia, Europe, America, Oceania and Africa. The aim of this study was to infer global patterns in pig domestication and diversity related to demography, migration, and selection. RESULTS: A deep phylogeographic division reflects the dichotomy between early domestication centers. In the core Eastern and Western domestication regions, Chinese pigs show differentiation between breeds due to geographic isolation, whereas this is less pronounced in European pigs. The inferred European origin of pigs in the Americas, Africa, and Australia reflects European expansion during the sixteenth to nineteenth centuries. Human-mediated introgression, which is due, in particular, to importing Chinese pigs into the UK during the eighteenth and nineteenth centuries, played an important role in the formation of modern pig breeds. Inbreeding levels vary markedly between populations, from almost no runs of homozygosity (ROH) in a number of Asian wild boar populations, to up to 20% of the genome covered by ROH in a number of Southern European breeds. Commercial populations show moderate ROH statistics. For domesticated pigs and wild boars in Asia and Europe, we identified highly differentiated loci that include candidate genes related to muscle and body development, central nervous system, reproduction, and energy balance, which are putatively under artificial selection. CONCLUSIONS: Key events related to domestication, dispersal, and mixing of pigs from different regions are reflected in the 60K SNP data, including the globalization that has recently become full circle since Chinese pig breeders in the past decades started selecting Western breeds to improve local Chinese pigs. Furthermore, signatures of ongoing and past selection, acting at different times and on different genetic backgrounds, enhance our insight in the mechanism of domestication and selection. The global diversity statistics presented here highlight concerns for maintaining agrodiversity, but also provide a necessary framework for directing genetic conservation.

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.

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.

Isolation, characterization and differentiation of cells expressing pluripotent/multipotent markers from adult human ovaries.

Pluripotent stem cells are still generally accepted not to exist in adult human ovaries, although increasing studies confirm the presence of pluripotent/multipotent stem cells in adult mammalian ovaries, including those of humans. The aim of this study is to isolate, characterize and differentiate in vitro stem cells that originate from the adult human ovarian cortex and that express markers of pluripotency/multipotency. After enzymatic degradation of small ovarian cortex biopsies retrieved from 18 women, ovarian cell cultures were successfully established from 17 and the formation of cell colonies was observed. The presence of cells/colonies expressing some markers of pluripotency (alkaline phosphatase, surface antigen SSEA-4, OCT4, SOX-2, NANOG, LIN28, STELLA), germinal lineage (DDX4/VASA) and multipotency (M-CAM/CD146, Thy-1/CD90, STRO-1) was confirmed by various methods. Stem cells from the cultures, including small round SSEA-4-positive cells with diameters of up to 4 µm, showed a relatively high degree of plasticity. We were able to differentiate them in vitro into various types of somatic cells of all three germ layers. However, these cells did not form teratoma when injected into immunodeficient mice. Our results thus show that ovarian tissue is a potential source of stem cells with a pluripotent/multipotent character for safe application in regenerative medicine.

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.

Identification of goat mammary stem/progenitor cells.

Goat mammary gland epithelial cells have been used to establish primary and permanent cell lines, but to date, no data have been available regarding mammary stem cells (MaSCs) in this species. The detection and characterization of goat MaSCs is an important task for a better understanding of the cyclic character of mammary gland development, which will also offer the potential for manipulation of lactation yield and persistency. The objective of the present study was to demonstrate that a subpopulation of goat MaSCs resides in the goat mammary gland. Mammary tissue from lactating Saanen goats (Capra hircus) was dissociated and processed to a single-cell suspension. Using an in vitro colony-forming assay, we demonstrated that distinct colony types, which expressed specific lineage markers, arose from unipotent progenitors. Using two different growth media, we showed that the frequencies of caprine clonogenic progenitors differed according to growth conditions. Goat epithelial cells were transplanted under the kidney capsule of nonobese diabetic/severe combined immunodeficient (NOD/SCID) mice, where they formed organized, bilayered structures. Our results indicate the presence of goat MaSCs in the caprine mammary gland. To our knowledge, these data represent the first description of the tissue hierarchy of the goat mammary gland and demonstrate the regenerative potential of adult goat MaSCs.

A simplified PCR-based method for detection of gray coat color allele in horse.

Coat color of gray horses is associated with a 4.6-kb duplication, which can be determined using PCR amplification of about 5-kb DNA fragment. In practice, this means that amplification might fail frequently. Therefore, a novel genetic screening method based on amplification of the 246 bp DNA fragment has been developed.

Impaired Neurodevelopmental Genes in Slovenian Autistic Children Elucidate the Comorbidity of Autism With Other Developmental Disorders.

Autism spectrum disorders (ASD) represent a phenotypically heterogeneous group of patients that strongly intertwine with other neurodevelopmental disorders (NDDs), with genetics playing a significant role in their etiology. Whole exome sequencing (WES) has become predominant in molecular diagnostics for ASD by considerably increasing the diagnostic yield. However, the proportion of undiagnosed patients still remains high due to complex clinical presentation, reduced penetrance, and lack of segregation analysis or clinical information. Thus, reverse phenotyping, where we first identified a possible genetic cause and then determine its clinical relevance, has been shown to be a more efficient approach. WES was performed on 147 Slovenian pediatric patients with suspected ASD. Data analysis was focused on identifying ultrarare or "single event" variants in ASD-associated genes and further expanded to NDD-associated genes. Protein function and gene prioritization were performed on detected clinically relevant variants to determine their role in ASD etiology and phenotype. Reverse phenotyping revealed a pathogenic or likely pathogenic variant in ASD-associated genes in 20.4% of patients, with subsequent segregation analysis indicating that 14 were de novo variants and 1 was presumed compound heterozygous. The diagnostic yield was further increased by 2.7% by the analysis of ultrarare or "single event" variants in all NDD-associated genes. Protein function analysis established that genes in which variants of unknown significance (VUS) were detected were predominantly the cause of intellectual disability (ID), and in most cases, features of ASD as well. Using such an approach, variants in rarely described ASD-associated genes, such as SIN3B, NR4A2, and GRIA1, were detected. By expanding the analysis to include functionally similar NDD genes, variants in KCNK9, GNE, and other genes were identified. These would probably have been missed by classic genotype-phenotype analysis. Our study thus demonstrates that in patients with ASD, analysis of ultrarare or "single event" variants obtained using WES with the inclusion of functionally similar genes and reverse phenotyping obtained a higher diagnostic yield despite limited clinical data. The present study also demonstrates that most of the causative genes in our cohort were involved in the syndromic form of ASD and confirms their comorbidity with other developmental disorders.

Large-scale mitogenome sequencing reveals consecutive expansions of domestic taurine cattle and supports sporadic aurochs introgression.

The contribution of domestic cattle in human societies is enormous, making cattle, along with other essential benefits, the economically most important domestic animal in the world today. To expand existing knowledge on cattle domestication and mitogenome diversity, we performed a comprehensive complete mitogenome analysis of the species (802 sequences, 114 breeds). A large sample was collected in South-east Europe, an important agricultural gateway to Europe during Neolithization and a region rich in cattle biodiversity. We found 1725 polymorphic sites (810 singletons, 853 parsimony-informative sites and 57 indels), 701 unique haplotypes, a haplotype diversity of 0.9995 and a nucleotide diversity of 0.0015. In addition to the dominant T3 and several rare haplogroups (Q, T5, T4, T2 and T1), we have identified maternal line in Austrian Murbodner cattle that possess surviving aurochs' mitochondria haplotype P1 that diverged prior to the Neolithization process. This is convincing evidence for rare female-mediated adaptive introgression of wild aurochs into domesticated cattle in Europe. We revalidated the existing haplogroup classification and provided Bayesian phylogenetic inference with a more precise estimated divergence time than previously available. Occasionally, classification based on partial mitogenomes was not reliable; for example, some individuals with haplogroups P and T5 were not recognized based on D-loop information. Bayesian skyline plot estimates (median) show that the earliest population growth began before domestication in cattle with haplogroup T2, followed by Q (~10.0-9.5 kyBP), whereas cattle with T3 (~7.5 kyBP) and T1 (~3.0-2.5 kyBP) expanded later. Overall, our results support the existence of interactions between aurochs and cattle during domestication and dispersal of cattle in the past, contribute to the conservation of maternal cattle diversity and enable functional analyses of the surviving aurochs P1 mitogenome.

How Can We Advance Integrative Biology Research in Animal Science in 21st Century? Experience at University of Ljubljana from 2002 to 2022.

In this perspective analysis, we strive to answer the following question: how can we advance integrative biology research in the 21st century with lessons from animal science? At the University of Ljubljana, Biotechnical Faculty, Department of Animal Science, we share here our three lessons learned in the two decades from 2002 to 2022 that we believe could inform integrative biology, systems science, and animal science scholarship in other countries and geographies. Cultivating multiomics knowledge through a conceptual lens of integrative biology is crucial for life sciences research that can stand the test of diverse biological, clinical, and ecological contexts. Moreover, in an era of the current COVID-19 pandemic, animal nutrition and animal science, and the study of their interactions with human health (and vice versa) through integrative biology approaches hold enormous prospects and significance for systems medicine and ecosystem health.

Mycoplasma gallisepticum and Mycoplasma synoviae express a cysteine protease CysP, which can cleave chicken IgG into Fab and Fc.

Major poultry pathogens M. gallisepticum and M. synoviae share a gene encoding a putative cysteine protease CysP similar to papain cysteine protease (C1A subfamily). Comparison of the cysP gene sequences of 18 M. synoviae and 10 M. gallisepticum strains sequenced in this study showed polymorphisms, including deletions. Seven M. synoviae strains, including the type strain WVU 1853, had a 39 bp deletion in the 3' end of the cysP gene. In the same cysP region, all M. gallisepticum strains showed a deletion of 66 bp. Immunoblot analysis with specific antibodies demonstrated that M. synoviae strains expressed CysP, which was approximately 65 kDa. Both M. synoviae and M. gallisepticum were able to digest chicken IgG (cIgG). Incubation of cIgG (∼170 kDa) with M. synoviae or M. gallisepticum cells (∼15 h at 37 °C) resulted in a papain-like cleavage pattern of cIgG and fragments corresponding to the antigen-binding fragment of IgG (Fab, ∼45 kDa) and the crystallizable region fragment (Fc) of the IgG heavy chain (dimer of ∼60 kDa). Iodoacetamide (50 mM) prevented cleavage of cIgG by both Mycoplasma species. Following site-directed mutagenesis (eight TGA codons were changed to TGG) the cysP gene of M. synoviae ULB 925 was expressed as a His-tagged protein in a cell-free system. Purified recombinant CysP (rCysP; ∼67 kDa, pI∼8) cleaved cIgG into Fab and Fc fragments. This indicates that CysP is responsible for the cIgG cleavage caused by M. synoviae and, probably, by M. gallisepticum. This is the first evidence to our knowledge that mycoplasmas have enzymes that can cleave the host IgG and indicates a novel strategy used by M. gallisepticum and M. synoviae for prolonged survival despite the antibody response of their host.