Less than 2% of mammalian genomes code for proteins, but ‘the majority of its bases can be found in primary transcripts’ – a phenomenon termed the pervasive transcription, which was first reported in 2007. Even though most of the transcripts do not code for proteins, they play a variety of biological functions, with regulation of gene expression appearing as the most common one. Those transcripts are divided into two groups based on their length: small non-coding RNAs, which are maximally 200 bp long, and long non-coding RNAs (lncRNAs), which are longer than 200 nucleotides. The advances in next-generation sequencing methods provided a new possibility of investigating the full set of RNA molecules in the cell. In this review, we summarized the current state of knowledge on lncRNAs in three major livestock species – Sus scrofa, Bos taurus and Gallus gallus, based on the literature and the content of biological databases. In the NONCODE database, the largest number of identified lncRNA transcripts is available for pigs, but cattle have the largest number of lncRNA genes. Poultry is represented by less than a half of records. Genomic annotation of lncRNAs showed that the majority of them are assigned to introns (pig, poultry) or intergenic (cattle). The comparison with well-annotated human and mouse genomes indicates that such annotation is a result of lack of proper lncRNA annotation data. Since lncRNAs play an important role in genomic studies, their characterization in farm animals’ genomes is critical in bridging the gap between genotype and phenotype.

Mastitis is an inflammatory disease of the mammary gland, which has a significant economic impact and is an animal welfare concern. This work examined the association between single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) with the incidence of clinical mastitis (CM). Using information from 16 half-sib pairs of Holstein-Friesian cows (32 animals in total) we searched for genomic regions that differed between a healthy (no incidence of CM) and a mastitis-prone (multiple incidences of CM) half-sib. Three cows with average sequence depth of coverage below 10 were excluded, which left 13 half-sib pairs available for comparisons. In total, 191 CNV regions were identified, which were deleted in a mastitis-prone cow, but present in its healthy half-sib and overlapped in at least nine half-sib pairs. These regions overlapped with exons of 46 genes, among which APP (BTA1), FOXL2 (BTA1), SSFA2 (BTA2), OTUD3 (BTA2), ADORA2A (BTA17), TXNRD2 (BTA17) and NDUFS6 (BTA20) have been reported to influence CM. Moreover, two duplicated CNV regions present in nine healthy individuals and absent in their mastitis-affected half-sibs overlapped with exons of a cholinergic receptor nicotinic α 10 subunit on BTA15 and a novel gene (ENSBTAG00000008519) on BTA27. One CNV region deleted in nine mastitis-affected sibs overlapped with two neighbouring long non-coding RNA sequences located on BTA12. Single nucleotide polymorphisms with differential genotypes between a healthy and a mastitis-affected sib included 17 polymorphisms with alternate alleles in eight affected and healthy half-sib families. Three of these SNPs were located introns of genes: MET (BTA04), RNF122 (BTA27) and WRN (BTA27). In summary, structural polymorphisms in form of CNVs, putatively play a role in susceptibility to CM. Specifically, sequence deletions have a greater effect on reducing resistance against mastitis, than sequence duplications have on increasing resistance against the disease.

The number of studies of Copy Number Variation in cattle has increased in recent years. This has been prompted by the increased availability of data on polymorphisms and their relationship with phenotypes. In addition, livestock species are good models for some human phenotypes. In the present study, we described the landscape of CNV driven genetic variation in a large population of 146 individuals representing 13 cattle breeds, using whole genome DNA sequence. A highly significant variation among all individuals and within each breed was observed in the number of duplications (P < 10−15) and in the number of deletions (P < 10−15). We also observed significant differences between breeds for duplication (P = 0.01932) and deletion (P = 0.01006) counts. The same variation CNV length - inter-individual and inter-breed differences were significant for duplications (P < 10−15) and deletions (P < 10−15). Moreover, breed-specific variants were identified, with the largest proportion of breed-specific duplications (9.57%) found for Fleckvieh and breed-specific deletions found for Brown Swiss (5.00%). Such breed-specific CNVs were predominantly located in intragenic regions, however in Simmental, one deletion present in five individuals was found in the coding sequence of a novel gene ENSBTAG00000000688 on chromosome 18. In Brown Swiss, Norwegian Red and Simmental breed-specific deletions were located within KIT and MC1R genes, which are responsible for a coat colour. The functional annotation of coding regions underlying the breed-specific CNVs showed that in Norwegian Red, Guernsey, and Simmental significantly under- and overrepresented GO terms were related to chemical stimulus involved in sensory perception of smell and the KEGG pathways for olfactory transduction. In addition, specifically for the Norwegian Red breed, the dopaminergic synapse KEGG pathway was significantly enriched within deleted parts of the genome. The CNV landscape in Bos taurus genome revealed by this study was highly complex, with inter-breed differences, but also a significant variation within breeds. The former, may explain some of the phenotypic differences among analysed breeds, and the latter contributes to within-breed variation available for selection.

The number of cases of pancreatic cancers in 2019 in Poland was 3852 (approx. 2% of all cancers). The course of the disease is very fast, and the average survival time from the diagnosis is 6 months. Only <2% of patients live for 5 years from the diagnosis, 8% live for 2 years, and almost half live for only about 3 months. A family predisposition to pancreatic cancer occurs in about 10% of cases. Several oncogenes in which somatic changes lead to the development of tumours, including genes BRCA1/2 and PALB2, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1, are involved in pancreatic cancer. Between 4% and 10% of individuals with pancreatic cancer will have a mutation in one of these genes. Six percent of patients with pancreatic cancer have NTRK pathogenic fusion. The pathogenesis of pancreatic cancer can in many cases be characterised by homologous recombination deficiency (HRD)—cell inability to effectively repair DNA. It is estimated that from 24% to as many as 44% of pancreatic cancers show HRD. The most common cause of HRD are inactivating mutations in the genes regulating this DNA repair system, mainly BRCA1 and BRCA2, but also PALB2, RAD51C and several dozen others.

Thirty-two whole genome DNA sequences of cows were analyzed to evaluate inter-individual variability in the distribution and length of copy number variations (CNV) and to functionally annotate CNV breakpoints. The total number of deletions per individual varied between 9,731 and 15,051, whereas the number of duplications was between 1,694 and 5,187. Most of the deletions (81%) and duplications (86%) were unique to a single cow. No relation between the pattern of variant sharing and a family relationship or disease status was found. The animal-averaged length of deletions was from 5,234 to 9,145 bp and the average length of duplications was between 7,254 and 8,843 bp. Highly significant inter-individual variation in length and number of CNV was detected for both deletions and duplications. The majority of deletion and duplication breakpoints were located in intergenic regions and introns, whereas fewer were identified in noncoding transcripts and splice regions. Only 1.35 and 0.79% of the deletion and duplication breakpoints were observed within coding regions. A gene with the highest number of deletion breakpoints codes for protein kinase cGMP-dependent type I, whereas the T-cell receptor α constant gene had the most duplication breakpoints. The functional annotation of genes with the largest incidence of deletion/duplication breakpoints identified 87/112 Kyoto Encyclopedia of Genes and Genomes pathways, but none of the pathways were significantly enriched or depleted with breakpoints. The analysis of Gene Ontology (GO) terms revealed that a cluster with the highest enrichment score among genes with many deletion breakpoints was represented by GO terms related to ion transport, whereas the GO term cluster mostly enriched among the genes with many duplication breakpoints was related to binding of macromolecules. Furthermore, when considering the number of deletion breakpoints per gene functional category, no significant differences were observed between the “housekeeping” and “strong selection” categories, but genes representing the “low selection pressure” group showed a significantly higher number of breakpoints.

In Bos taurus the universality of the reference genome is biased towards genetic variation represented by only two related individuals representing the same Hereford breed. Therefore, results of genetic analyses based on this reference may not be reliable. The 1000 Bull Genomes resource allows for identification of breed-specific polymorphisms and for the construction of breed-specific reference genomes. Whole-genome sequences or 936 bulls allowed us to construct seven breed specific reference genomes of Bos taurus for Angus, Brown Swiss, Fleckvieh, Hereford, Jersey, Limousin and Simmental. In order to identify breed-specific variants all detected SNPs were filtered within-breed to satisfy criteria of the number of missing genotypes not higher than 7% and the alternative allele frequency equal to unity. The highest number of breed-specific SNPs was identified for Jersey (130,070) and the lowest—for the Simmental breed (197). Such breed-specific polymorphisms were annotated to coding regions overlapping with 78 genes in Angus, 140 in Brown Swiss, 132 in Fleckvieh, 100 in Hereford, 643 in Jersey, 10 in Limousin and no genes in Simmental. For most of the breeds, the majority of breed-specific variants from coding regions was synonymous. However, most of Fleckvieh-specific and Hereford-specific polymorphisms were missense mutations. Since the identified variants are characteristic for the analysed breeds, they form the basis of phenotypic differences observed between them, which result from different breeding programmes. Breed-specific reference genomes can enhance the accuracy of SNP driven inferences such as Genome-wide Association Studies or SNP genotype imputation.