Influence of breed and environment on leukocyte telomere length in cattle

High milk yield is associated with reduced longevity in high-producing dairy cattle breeds. Pre-term culling leads to high replacement heifer demand and economic losses for the dairy industry. Selection for this trait is limited because of low heritability and difficulties in phenotype measurement. Telomeres are elements found at the ends of chromosomes, consisting of repetitive DNA sequences, several thousand base pairs in length, coupled with nucleoprotein complexes. Eventually, in humans and most other animals, telomere length reduces with age. When telomeric DNA is truncated to a critical length, cell ageing, cell cycle arrest, and apoptosis are induced. As a result, telomere length can be considered as a predictor of health risks and an individual's lifespan. The leukocyte telomere length may be used as a proxy phenotype of productive lifespan to improve cattle selection. Our objectives were to assess the effects of breed and breed group (dairy vs. beef) on the leukocyte telomere length and to estimate the effect of cold climate on this trait in Kalmyk cattle populations from the South (Rostov Oblast) and Far North (Republic of Sakha) regions of Russia. The leukocyte telomere lengths were estimated computationally from whole-genome resequencing data. We leveraged data on leukocyte telomere length, sex, and age of 239 animals from 17 cattle breeds. The breed factor had a significant effect on leukocyte telomere length across our sample. There was no difference in leukocyte telomere length between dairy and beef groups. The population factor had a significant effect on leukocyte telomere length in Kalmyk animals. In conclusion, we found that breed, but not breed group (dairy vs. beef), was significantly associated with leukocyte telomere length in cattle. Residence in colder climates was associated with longer leukocyte telomere length in Kalmyk breed cattle.

Association of three single nucleotide polymorphisms in the LPIN1 gene with milk production traits in cows of the Yaroslavl breed

Lipin-1 is a member of the evolutionarily conserved family of proteins and is expressed predominantly in adipose tissue and skeletal muscle. On the one hand, lipin-1 is an enzyme that catalyzes the dephosphorylation of phosphatidic acid to diacylglycerol (DAG) and thus participates in the metabolic pathways of biosynthesis of storage lipids in the cell, membrane phospholipids, and intracellular signaling molecules. On the other hand, lipin-1 is able to be transported from the cytoplasm to the nucleus and is a coactivator of lipid metabolism gene transcription. It was shown, using the analysis of single nucleotide polymorphism (SNP) associations, that the lipin-1 coding gene (LPIN1) is a promising candidate gene for milk production traits in Holstein and Brown Swiss cows. However, it is unclear how much of its effect depends on the breed. The Yaroslavl dairy cattle breed was created in the 18-19 centuries in Russia by breeding northern Great Russian cattle, which were short and poor productive, but well adapted to local climatic conditions and bad food base. It was shown by whole genome genotyping and sequencing that the Yaroslavl breed has unique genetics compared to Russian and other cattle breeds. The aim of the study was to assess the frequency of alleles and genotypes of three SNPs in the LPIN1 gene and to study the association of these SNPs with milk production traits in Yaroslavl cows. Blood samples from 142 cows of the Yaroslavl breed were obtained from two farms in the Yaroslavl region. Genotyping of SNPs was carried out by polymerase chain reaction-restriction fragment length polymorphism method. Associations of SNPs with 305-day milk yield, fat yield, fat percentages, protein yield, and protein percentages were studied from the first to the fourth lactation. Statistical tests were carried out using a mixed linear model, taking into account the relationship between individuals. We identified three SNPs - rs110871255, rs207681322 and rs109039955 with a frequency of a rare allele of 0.042-0.261 in Yaroslavl cows. SNP rs110871255 was associated with fat yield during the third and fourth lactations. SNP rs207681322 was associated with milk yield for the second, third and fourth lactations, as well as protein yield for the third lactation. Thus, we identified significant associations of SNPs rs207681322 and rs110871255 in the LPIN1 gene with a number of milk production traits during several lactations in Yaroslavl cows.

Candidate genes for domestication and resistance to cold climate according to whole genome sequencing data of Russian cattle and sheep breeds

It is known that different species of animals, when living in the same environmental conditions, can form similar phenotypes. The study of the convergent evolution of several species under the influence of the same environmental factor makes it possible to identify common mechanisms of genetic adaptation. Local cattle and sheep breeds have been formed over thousands of years under the influence of domestication, as well as selection aimed at adaptation to the local environment and meeting human needs. Previously, we identified a number of candidate genes in genome regions potentially selected during domestication and adaptation to the climatic conditions of Russia, in local breeds of cattle and sheep using whole genome genotyping data. However, these data are of low resolution and do not reveal most nucleotide substitutions. The aim of the work was to create, using the whole genome sequencing data, a list of genes associated with domestication, selection and adaptation in Russian cattle and sheep breeds, as well as to identify candidate genes and metabolic pathways for selection for cold adaptation. We used our original data on the search for signatures of selection in the genomes of Russian cattle (Yakut, Kholmogory, Buryat, Wagyu) and sheep (Baikal, Tuva) breeds. We used the HapFLK, DCMS, FST and PBS methods to identify DNA regions with signatures of selection. The number of candidate genes in potentially selective regions was 946 in cattle and 151 in sheep. We showed that the studied Russian cattle and sheep breeds have at least 10 genes in common, apparently involved in the processes of adaptation/selection, including adaptation to a cold climate, including the ASTN2, PM20D1, TMEM176A, and GLIS1 genes. Based on the intersection with the list of selected genes in at least two Arctic/Antarctic mammal species, 20 and 8 genes, have been identified in cattle and sheep, respectively, that are potentially involved in cold adaptation. Among them, the most promising for further research are the ASPH, NCKAP5L, SERPINF1, and SND1 genes. Gene ontology analysis indicated the existence of possible common biochemical pathways for adaptation to cold in domestic and wild mammals associated with cytoskeleton disassembly and apoptosis.

Compilation and functional classification of telomere length-associated genes in humans and other animal species

Telomeres are the terminal regions of chromosomes that ensure their stability while cell division. Telomere shortening initiates cellular senescence, which can lead to degeneration and atrophy of tissues, so the process is associated with a reduction in life expectancy and predisposition to a number of diseases. An accelerated rate of telomere attrition can serve as a predictor of life expectancy and health status of an individual. Telomere length is a complex phenotypic trait that is determined by many factors, including the genetic ones. Numerous studies (including genome-wide association studies, GWAS) indicate the polygenic nature of telomere length control. The objective of the present study was to characterize the genetic basis of the telomere length regulation using the GWAS data obtained during the studies of various human and other animal populations. To do so, a compilation of the genes associated with telomere length in GWAS experiments was collected, which included information on 270 human genes, as well as 23, 22, and 9 genes identified in the cattle, sparrow, and nematode, respectively. Among them were two orthologous genes encoding a shelterin protein (POT1 in humans and pot-2 in C. elegans). Functional analysis has shown that telomere length can be influenced by genetic variants in the genes encoding: (1) structural components of telomerase; (2) the protein components of telomeric regions (shelterin and CST complexes); (3) the proteins involved in telomerase biogenesis and regulating its activity; (4) the proteins that regulate the functional activity of the shelterin components; (5) the proteins involved in telomere replication and/or capping; (6) the proteins involved in the alternative telomere lengthening; (7) the proteins that respond to DNA damage and are responsible for DNA repair; (8) RNA-exosome components. The human genes identified by several research groups in populations of different ethnic origins are the genes encoding telomerase components such as TERC and TERT as well as STN1 encoding the CST complex component. Apparently, the polymorphic loci affecting the functions of these genes may be the most reliable susceptibility markers for telomere-related diseases. The systematized data about the genes and their functions can serve as a basis for the development of prognostic criteria for telomere length-associated diseases in humans. Information about the genes and processes that control telomere length can be used for marker-assisted and genomic selection in the farm animals, aimed at increasing the duration of their productive lifetime.

Comparative analysis of allele frequencies for DNA polymorphisms associated with disease and economically important traits in the genomes of Russian and foreign cattle breeds

The genetic makeup of a breed including its genetic differences from other breeds determines its appearance and characteristics, including economically important traits and resistance to pathologies. To date, many loci controlling significant phenotypes have been identified, which is successfully used in the world practice of marker-assisted selection to improve breed properties. The aim of this study was a comparative analysis of frequencies for known causative nucleotide substitutions, insertions and deletions associated with disease and economically important traits in Russian and foreign cattle breeds. As a result, we identified frequencies of these DNA polymorphisms in the populations of Russian cattle breeds, compared them with those of foreign populations of the same breed, as well as other foreign breeds. Our results indicate similarities in frequencies for most of such alleles within breeds (populations of Russian and foreign breeding), as well as the relationship between the causative allele prevalence and the presence of phenotypic traits under the effect. We also found an excess of some undesirable alleles in the Russian cattle populations, which should be paid attention to when designing breeding programs. We found that the alleles increasing fertility in the Hereford breed have a higher frequency in the Russian Hereford population compared to the foreign counterpart. Interestingly, unlike for the European breeds, for Asian Turano-Mongolian Wagyu and Yakut cattle, there was a less clear link between phenotypic traits and frequencies of known causative alleles. Our work points to specific genetic variants that could be used to improve and/or maintain the performance of certain cattle breeds bred in the Russian Federation.

[Signatures of selection and candidate genes for adaptation to extreme environmental factors in the genomes of Turano-Mongolian cattle breeds]

Changes in the environment force populations of organisms to adapt to new conditions, either through phenotypic plasticity or through genetic or epigenetic changes. Signatures of selection, such as specific changes in the frequency of alleles and haplotypes, as well as the reduction or increase in genetic diversity, help to identify changes in the cattle genome in response to natural and artificial selection, as well as loci and genetic variants directly affecting adaptive and economically important traits. Advances in genetics and biotechnology enable a rapid transfer of unique genetic variants that have originated in local cattle breeds in the process of adaptation to local environments into the genomes of cosmopolitan high-performance breeds, in order to preserve their outstanding performance in new environments. It is also possible to use genomic selection approach to increase the frequency of already present adaptive alleles in cosmopolitan breeds. The review examines recent work on the origin and evolution of Turano-Mongolian cattle breeds, adaptation of Turano-Mongolian cattle to extreme environments, and summarizes available information on potential candidate genes for climate adaptation of Turano-Mongolian breeds, including cold resistance genes, immune response genes, and high-altitude adaptation genes. The authors conclude that the current literature data do not provide preference to one of the two possible scenarios of Turano-Mongolian breed origins: as a result of the domestication of a wild aurochs at East Asia or as a result of the migration of taurine proto-population from the Middle East. Turano-Mongolian breeds show a high degree of adaptation to extreme climatic conditions (cold, heat, lack of oxygen in the highlands) and parasites (mosquitoes, ticks, bacterial and viral infections). As a result of high-density genotyping and sequencing of genomes and transcriptomes, prospective candidate genes and genetic variants involved in adaptation to environmental factors have recently been identified.

Copy number variants in genomes of local sheep breeds from Russia

Copy number variants (CNVs) are genomic structural variations that contribute to many adaptive and economically important traits in livestock. In this study, we detected CNVs in 354 animals from 16 Russian indigenous sheep breeds and analysed their possible functional roles. Our analysis of the entire sample set resulted in 4527 CNVs forming 1450 CNV regions (CNVRs). When constructing CNVRs for individual breeds, a total of 2715 regions ranging from 88 in Groznensk to 337 in Osetin breeds were identified. To make interbreed CNVR frequency comparison possible, we also identified core CNVRs using CNVs with overlapping chromosomal locations found in different breeds. This resulted in 137 interbreed CNVRs with frequency >15% in at least one breed. Functional enrichment analysis of genes affected by CNVRs in individual breeds revealed 12 breeds with significant enrichments in olfactory perception, PRAME family proteins, and immune response. Function of genes affected by interbreed and breed-specific CNVRs revealed candidates related to domestication, adaptation to high altitudes and cold climates, reproduction, parasite resistance, milk and meat qualities, wool traits, fat storage, and fat metabolism. Our work is the first attempt to uncover and characterise the CNV makeup of Russian indigenous sheep breeds. Further experimental and functional validation of CNVRs would help in developing new and improving existing sheep breeds.

Resequencing and signatures of selection scan in two Siberian native sheep breeds point to candidate genetic variants for adaptation and economically important traits

Russian sheep breeds represent an important economic asset by providing meat and wool, whilst being adapted to extreme climates. By resequencing two Russian breeds from Siberia: Tuva (n = 20) and Baikal (n = 20); and comparing them with a European (UK) sheep outgroup (n = 14), 41 million variants were called, and signatures of selection were identified. High-frequency missense mutations on top of selection peaks were found in genes related to immunity (LOC101109746) in the Baikal breed and wool traits (IDUA), cell differentiation (GLIS1) and fat deposition (AADACL3) in the Tuva breed. In addition, genes found under selection owing to haplotype frequency changes were related to wool traits (DSC2), parasite resistance (CLCA1), insulin receptor pathway (SOCS6) and DNA repair (DDB2) in the Baikal breed, and vision (GPR179) in the Tuva breed. Our results present candidate genes and SNPs for future selection programmes, which are necessary to maintain and increase socioeconomic gain from Siberian breeds.

Whole genome studies of origin, selection and adaptation of the Russian cattle breeds

Our review presents several recent studies on the genetic history and signatures of selection in genomes of the native Russian cattle breeds. Most of these works are not easily accessible for the Russian-speaking audience. We describe the origins of appearance of the Russian cattle breeds from the genetics perspective. We point to the links between most of the Russian breeds with the taurine breeds of the European origin and for some Russian breeds with the breeds of the Asian origin. We describe major phylogenetic clusters of the Russian breeds and point to those that still maintain their unique genetics, meaning that their preservation is a priority. In addition, we review the results of the search for signatures of selection in genomes of the Russian cattle breeds. Some unique signatures of selection present in the genomes of so-called “turano-mongolian” cattle (i. e. the Yakut cattle) are described which allowed the Yakut cattle to adapt to harsh environments found above the Polar Circle. Signatures of selection which could help other cattle breeds of the Russian origin to adapt to various climatic condition of the Russian Federation are reviewed. The Russian cattle genomes also contain known signatures of selection related to cattle domestication about 8–10 thousand years ago. The most profound ones include genes related to changes of the coat colour. This phenotype in many cases could be related to the distinction of the first domesticated populations and lead to the formation of so-called land races (primitive breeds). Whole-genome association studies of Russian cattle breeds pointed to a novel gene which could be related to the “white-faced” phenotype and to a gene which is related to body temperature support under the acute cold stress. The data presented in our review could be used for identification of genetic markers to focus on in future efforts on designing new highly productive cattle breeds adapted to climates of the Russian Federation and other countries with similar climates.

Genome-wide association study for body weight in cattle populations from Siberia

Body weight is a complex trait in cattle associated with commonly used commercial breeding measurements related to growth. Although many quantitative trait loci (QTL) for body weight have been identified in cattle so far, searching for genetic determinants in different breeds or environments is promising. Therefore, we carried out a genome-wide association study (GWAS) in two cattle populations from the Russian Federation (Siberian region) using the GGP HD150K array containing 139 376 single nucleotide polymorphism (SNP) markers. Association tests for 107 550 SNPs left after filtering revealed five statistically significant SNPs on BTA5, considering a false discovery rate of less than 0.05. The chromosomal region containing these five SNPs contains the CCND2 gene, which was previously associated with average daily weight gain and body mass index in US beef cattle populations and in humans respectively. Our study is the first GWAS for body weight in beef cattle populations from the Russian Federation. The results provided here suggest that, despite the existence of breed- and species-specific QTL, the genetic architecture of body weight could be evolutionarily conserved in mammals.

Genes related to the white face colour pattern in eight Russian cattle breeds

One of the major effects of domestication is change of animal coat colour to up to complete white colour of the whole body. It is possible that white colour of livestock animals had aesthetic significance for humans as well. The first step towards detection of genes and mutations controlling white colouring in animals is the genome-wide association studies. These studies, however, have not been done for the cattle breeds native to the Russian Federation. The aim of this study was therefore to identify genomic intervals and candidate genes that could be responsible for white face colouring in eight Russian cattle breeds. The data on genome-wide genotyping of 131,709 high-quality single nucleotide polymorphisms (SNPs) on 148 animas have been used in the program ­EMMAX. Association analysis has been performed using two related phenotypes: a) the white face with the rest of the body of any colour and b) white face with the rest of the body of different (non-white) colour. In the first case, the only statistically significant marker found was the SNP BovineHD0500019319 located on cattle chromosome (BTA) 5. The same SNP was the most significant within the cluster of three SNPs on BTA5: 68,803,879–69,365,854 associated also with the second phenotype. Five genes were found within this interval in the cattle genome, out of which the most likely functional candidate was SLC41A2, with the SNP BovineHD0500019319 found within its intronic sequence. SLC41A2 encodes a magnesium transporter protein. However, the function of this gene is not well established. Other members of this gene family are the key genes controlling differences in human skin and animal coat colour. Additional significant association signals with the second phenotype have been detected in BTA 1–4, 6–15, 18, 19, 24, 27, and 29. Overall, 37 genomic intervals have been detected associated with white face colouring in eight Russian native cattle breeds.

Status of the cattle genome map

During the last 30 years, the cattle genome map has been expanded from 4 genes linked on chromosome X to over 22,000 genes identified in the cattle genome sequence assembly. This progress has been achieved due to numerous projects on linkage and physical mapping of the cattle genome driven by its agricultural and scientific significance. Indeed, the high-resolution mapping and functional analysis of the genome led to the discovery of major quantitative trait loci (QTL) regions and several quantitative trait nucleotides (QTNs), as well as some disease genes in the cow population. In addition, a comparison of the cattle genome to the genomes of other mammals has revealed its unique features gained during the speciation and adaptation. With the development of non-expensive sequencing techniques, the analysis of the cattle genome will shift towards the identification of differences between breeds or individuals within breeds that account for the unique features of each breed. This approach holds promise for the development of effective tools for the marker assistant selection and disease diagnostics in cattle.

[Role of chromosomal rearrangement and conservative chromosome sites in amniote evolution]

Specific features of evolution of chromosomes in mammals and other amniotes are reviewed. Comparative analysis of chromosome architecture revealed nonrandom distribution of chromosome rearrangement sites in genome, probable role of chromosome rearrangement in adaptation, and evolution-mediated selection of conservative chromosome sites. The chromosome sites stable during evolution are saturated with genes contributing to early organism development. Ruptures of these sites are incompatible with organism survival. Further analysis of chromosome evolution requires more information about completely sequenced genomes.

Discovery, validation and characterization of 1039 cattle single nucleotide polymorphisms

We identified approximately 13 000 putative single nucleotide polymorphisms (SNPs) by comparison of repeat-masked BAC-end sequences from the cattle RPCI-42 BAC library with whole-genome shotgun contigs of cattle genome assembly Btau 1.0. Genotyping of a subset of these SNPs was performed on a panel containing 186 DNA samples from 18 cattle breeds including 43 trios. Of 1039 SNPs confirmed as polymorphic in the panel, 998 had minor allele frequency > or =0.25 among unrelated individuals of at least one breed. When Btau 4.0 became available, 974 of these validated SNPs were assigned in silico to known cattle chromosomes, while 41 SNPs were mapped to unassigned sequence scaffolds, yielding one SNP every approximately 3 Mbp on average. Twenty-four SNPs identified in Btau 1.0 were not mapped to Btau 4.0. Of the 1015 SNPs mapped to Btau 4.0, 959 SNPs had nucleotide bases identical in Btau 4.0 and Btau 1.0 contigs, whereas 56 bases were changed, resulting in the loss of the in silico SNP in Btau 4.0. Because these 1039 SNPs were all directly confirmed by genotyping on the multi-breed panel, it is likely that the original polymorphisms were correctly identified. The 1039 validated SNPs identified in this study represent a new and useful resource for genome-wide association studies and applications in animal breeding.

[Fate of parental mitochondria in embryonic stem hybrid cells]

When hybrid cells are created, not only the nuclear genomes of parental cells unite but their cytoplasm as well. Mitochondrial DNA (mtDNA) is a convenient marker of cytoplasm allowing one to gain insight into the organization of hybrid cell cytoplasm. We analyzed the parental mtDNAs in hybrid cells resulting from fusion of Mus musculus embryonic stem (ES) cells with splenocytes and fetal fibroblasts of DD/c mice or with splenocytes of M. caroli. Identification of the parental mtDNAs in hybrid cells was based on polymorphism among the parental mtDNAs for certain restrictases. We found that intra- and inter-specific ES cell-splenocyte hybrid cells lost entirely or partially mtDNA derived from the somatic partner, whereas ES cell-fibroblast hybrids retained mtDNAs from both parents in similar ratios with a slight bias. The lost of the “somatic” mitochondria by ES-splenocyte hybrids implies non-random segregation of the parental mitochondria as supported by a computer simulation of genetic drift. In contrast, ES cell-fibroblast hybrids show bilateral random segregation of the parental mitochondria judging from analysis of mtDNA in single cells. Preferential segregation of “somatic” mitochondria does not depend on the differences in sequences of the parental mtDNAs but depends on replicative state of the parental cells.

Comparative mapping of cattle chromosome 19: cytogenetic localization of 19 BAC clones

Here we present the results of fluorescent in situ hybridization (FISH) mapping of a set of cattle BAC clones preselected for assignment on cattle chromosome 19 (BTA19). The BAC clones were anchored to human chromosome 17 (HSA17) sequences by BLASTn similarity search of cattle BAC-ends against the human genome sequence (NCBI build 33). Five blocks of homologous synteny were defined in the comparative map of BTA19 and HSA17 built with FISH data and the human genome coordinates. The positions for four evolutionary breakpoints in the bovine and human chromosomes were identified. Comparison of the FISH comparative map with previously published comparative RH, physical, and cytogenetic maps of BTA19 did not reveal major conflicts and allowed for the extension of the boundaries of homology between BTA19 and HSA17. Comparative analysis of HSA17, BTA19, and mouse chromosome 11 (MMU11) demonstrates that most likely mice retain the ancestral organization of the synteny group, and both cattle and human chromosomes underwent several major internal rearrangements after the divergence of Primates, Rodentia, and Cetartiodactyla.

[The structure of a conserved region of porcine genome, represented in human genome by chromosome 17]

Radiation mapping of nine genes (H3F3B, HLR1, MYL4, STAT5B, THRA1, TOP2A, MCP1, NF1, and MPO) to porcine chromosome 12 was carried out. Also, subchromosomal location of the NF1 gene along with the two loci containing the DNA sequences homologous to the DNA of the two human BAC clones was determined. The NF1 position was ascertained via microdissection of chromosome 12 with subsequent PCR amplification of the gene fragment with specific primers. BAC clones were mapped using FISH. Comparative analysis of the gene order in porcine chromosome 12 and in the homologous human chromosome 17 was performed. It was demonstrated that the gene orders in these chromosomes differed relative to the position of the MPO gene.

[Use of PCR markers for mapping swine chromosome 12]

Using PCR analysis of pig-mink and pig-Chinese hamster hybrid cell lines and heterologous and homologous primers of various types, chromosomal and subchromosomal mapping of genes TOP2A, THRA, BRCA1, GAS, HLR1, MYL4, LIS1, MCP1, ENO3, CRYB1, P4HB, STAT5B, and H3F3B to pig chromosome 12 was carried out. The efficiency of using different types of heterologous primers for pig chromosome mapping was compared.

[Chromosome localization and analysis of synteny analysis of some genes in swine, cattle, and sheep (Artiodactyla)]

Using the hybrid cell lines pig-American mink, cow-American mink, and sheep-American mink, the localization of some genes included in a large conservative block localized on human chromosome (chr) 17 was performed by means of electrophoresis of proteins and Southern blot hybridization. Genes NF1, RARA, PRKCA, and ERBB2 were assigned to chr 12 in swine; TK1 and UMPH2, to chr 19 in cattle; and TK1, UMPH2, and PEPA, to chr 11 in sheep. The conserved synteny of these genes in three representatives of the order Artiodactyla was shown.