Innovative construction of the first reliable catalogue of bovine circular RNAs.

The aim of this study was to compare the circular transcriptome of divergent tissues in order to understand: i) the presence of circular RNAs (circRNAs) that are not exonic circRNAs, i.e. originated from backsplicing involving known exons and, ii) the origin of artificial circRNA (artif_circRNA), i.e. circRNA not generated in-vivo. CircRNA identification is mostly an in-silico process, and the analysis of data from the BovReg project (https://www.bovreg.eu/) provided an opportunity to explore new ways to identify reliable circRNAs. By considering 117 tissue samples, we characterized 23,926 exonic circRNAs, 337 circRNAs from 273 introns (191 ciRNAs, 146 intron circles), 108 circRNAs from small non-coding genes and nearly 36.6K circRNAs classified as other_circRNAs. Furthermore, for 63 of those samples we analysed in parallel data from total-RNAseq (ribosomal RNAs depleted prior to library preparation) with paired mRNAseq (library prepared with poly(A)-selected RNAs). The high number of circRNAs detected in mRNAseq, and the significant number of novel circRNAs, mainly other_circRNAs, led us to consider all circRNAs detected in mRNAseq as artificial. This study provided evidence of 189 false entries in the list of exonic circRNAs: 103 artif_circRNAs identified by total RNAseq/mRNAseq comparison using two circRNA tools, 26 probable artif_circRNAs, and 65 identified by deep annotation analysis. Extensive benchmarking was performed (including analyses with CIRI2 and CIRCexplorer-2) and confirmed 94% of the 23,737 reliable exonic circRNAs. Moreover, this study demonstrates the effectiveness of a panel of highly expressed exonic circRNAs (5-8%) in analysing the tissue specificity of the bovine circular transcriptome.

Genomic analyses of withers height and linear conformation traits in German Warmblood horses using imputed sequence-level genotypes.

BACKGROUND: Body conformation, including withers height, is a major selection criterion in horse breeding and is associated with other important traits, such as health and performance. However, little is known about the genomic background of equine conformation. Therefore, the aim of this study was to use imputed sequence-level genotypes from up to 4891 German Warmblood horses to identify genomic regions associated with withers height and linear conformation traits. Furthermore, the traits were genetically characterised and putative causal variants for withers height were detected. RESULTS: A genome-wide association study (GWAS) for withers height confirmed the presence of a previously known quantitative trait locus (QTL) on Equus caballus (ECA) chromosome 3 close to the LCORL/NCAPG locus, which explained 16% of the phenotypic variance for withers height. An additional significant association signal was detected on ECA1. Further investigations of the region on ECA3 identified a few promising candidate causal variants for withers height, including a nonsense mutation in the coding sequence of the LCORL gene. The estimated heritability for withers height was 0.53 and ranged from 0 to 0.34 for the conformation traits. GWAS identified significantly associated variants for more than half of the investigated conformation traits, among which 13 showed a peak on ECA3 in the same region as withers height. Genetic parameter estimation revealed high genetic correlations between these traits and withers height for the QTL on ECA3. CONCLUSIONS: The use of imputed sequence-level genotypes from a large study cohort led to the discovery of novel QTL associated with conformation traits in German Warmblood horses. The results indicate the high relevance of the QTL on ECA3 for various conformation traits, including withers height, and contribute to deciphering causal mutations for body size in horses.

Meta-analysis of six dairy cattle breeds reveals biologically relevant candidate genes for mastitis resistance.

BACKGROUND: Mastitis is a disease that incurs significant costs in the dairy industry. A promising approach to mitigate its negative effects is to genetically improve the resistance of dairy cattle to mastitis. A meta-analysis of genome-wide association studies (GWAS) across multiple breeds for clinical mastitis (CM) and its indicator trait, somatic cell score (SCS), is a powerful method to identify functional genetic variants that impact mastitis resistance. RESULTS: We conducted meta-analyses of eight and fourteen GWAS on CM and SCS, respectively, using 30,689 and 119,438 animals from six dairy cattle breeds. Methods for the meta-analyses were selected to properly account for the multi-breed structure of the GWAS data. Our study revealed 58 lead markers that were associated with mastitis incidence, including 16 loci that did not overlap with previously identified quantitative trait loci (QTL), as curated at the Animal QTLdb. Post-GWAS analysis techniques such as gene-based analysis and genomic feature enrichment analysis enabled prioritization of 31 candidate genes and 14 credible candidate causal variants that affect mastitis. CONCLUSIONS: Our list of candidate genes can help to elucidate the genetic architecture underlying mastitis resistance and provide better tools for the prevention or treatment of mastitis, ultimately contributing to more sustainable animal production.

Cows with diverging haplotypes show differences in differential milk cell count, milk parameters and vaginal temperature after S. aureus challenge but not after E. coli challenge.

BACKGROUND: In dairy cattle, mastitis causes high financial losses and impairs animal well-being. Genetic selection is used to breed cows with reduced mastitis susceptibility. Techniques such as milk cell flow cytometry may improve early mastitis diagnosis. In a highly standardized in vivo infection model, 36 half-sib cows were selected for divergent paternal Bos taurus chromosome 18 haplotypes (Q vs. q) and challenged with Escherichia coli for 24 h or Staphylococcus aureus for 96 h, after which the samples were analyzed at 12 h intervals. Vaginal temperature (VT) was recorded every three minutes. The objective of this study was to compare the differential milk cell count (DMCC), milk parameters (fat %, protein %, lactose %, pH) and VT between favorable (Q) and unfavorable (q) haplotype cows using Bayesian models to evaluate their potential as improved early indicators of differential susceptibility to mastitis. RESULTS: After S. aureus challenge, compared to the Q half-sibship cows, the milk of the q cows exhibited higher PMN levels according to the DMCC (24 h, p < 0.001), a higher SCC (24 h, p < 0.01 and 36 h, p < 0.05), large cells (24 h, p < 0.05) and more dead (36 h, p < 0.001) and live cells (24 h, p < 0.01). The protein % was greater in Q milk than in q milk at 0 h (p = 0.025). In the S. aureus group, Q cows had a greater protein % (60 h, p = 0.048) and fat % (84 h, p = 0.022) than q cows. Initially, the greater VT of S. aureus-challenged q cows (0 and 12-24 h, p < 0.05) reversed to a lower VT in q cows than in Q cows (48-60 h, p < 0.05). Additionally, the following findings emphasized the validity of the model: in the S. aureus group all DMCC subpopulations (24 h-96 h, p < 0.001) and in the E. coli group nearly all DMCC subpopulations (12 h-24 h, p < 0.001) were higher in challenged quarters than in unchallenged quarters. The lactose % was lower in the milk samples of E. coli-challenged quarters than in those of S. aureus-challenged quarters (24 h, p < 0.001). Between 12 and 18 h, the VT was greater in cows challenged with E. coli than in those challenged with S. aureus (3-h interval approach, p < 0.001). CONCLUSION: This in vivo infection model confirmed specific differences between Q and q cows with respect to the DMCC, milk component analysis results and VT results after S. aureus inoculation but not after E. coli challenge. However, compared with conventional milk cell analysis monitoring, e.g., the global SCC, the DMCC analysis did not provide refined phenotyping of the pathogen response.

Sequence-based GWAS meta-analyses for beef production traits.

BACKGROUND: Combining the results of within-population genome-wide association studies (GWAS) based on whole-genome sequences into a single meta-analysis (MA) is an accurate and powerful method for identifying variants associated with complex traits. As part of the H2020 BovReg project, we performed sequence-level MA for beef production traits. Five partners from France, Switzerland, Germany, and Canada contributed summary statistics from sequence-based GWAS conducted with 54,782 animals from 15 purebred or crossbred populations. We combined the summary statistics for four growth, nine morphology, and 15 carcass traits into 16 MA, using both fixed effects and z-score methods. RESULTS: The fixed-effects method was generally more informative to provide indication on potentially causal variants, although we combined substantially different traits in each MA. In comparison with within-population GWAS, this approach highlighted (i) a larger number of quantitative trait loci (QTL), (ii) QTL more frequently located in genomic regions known for their effects on growth and meat/carcass traits, (iii) a smaller number of genomic variants within the QTL, and (iv) candidate variants that were more frequently located in genes. MA pinpointed variants in genes, including MSTN, LCORL, and PLAG1 that have been previously associated with morphology and carcass traits. We also identified dozens of other variants located in genes associated with growth and carcass traits, or with a function that may be related to meat production (e.g., HS6ST1, HERC2, WDR75, COL3A1, SLIT2, MED28, and ANKAR). Some of these variants overlapped with expression or splicing QTL reported in the cattle Genotype-Tissue Expression atlas (CattleGTEx) and could therefore regulate gene expression. CONCLUSIONS: By identifying candidate genes and potential causal variants associated with beef production traits in cattle, MA demonstrates great potential for investigating the biological mechanisms underlying these traits. As a complement to within-population GWAS, this approach can provide deeper insights into the genetic architecture of complex traits in beef cattle.

Metabogenomic analysis to functionally annotate the regulatory role of long non-coding RNAs in the liver of cows with different nutrient partitioning phenotype.

Long non-coding RNAs (lncRNAs) hold gene regulatory potential, but require substantial further functional annotation in livestock. Applying two metabogenomic approaches by combining transcriptomic and metabolomic analyses, we aimed to identify lncRNAs with potential regulatory function for divergent nutrient partitioning of lactating crossbred cows and to establish metabogenomic interaction networks comprising metabolites, genes and lncRNAs. Through correlation analysis of lncRNA expression with transcriptomic and metabolomic data, we unraveled lncRNAs that have a putative regulatory role in energy and lipid metabolism, the urea and tricarboxylic acid cycles, and gluconeogenesis. Especially FGF21, which correlated with a plentitude of differentially expressed genes, differentially abundant metabolites, as well as lncRNAs, suggested itself as a key metabolic regulator. Notably, lncRNAs in close physical proximity to coding-genes as well as lncRNAs with natural antisense transcripts appear to perform a fine-tuning function in gene expression involved in metabolic pathways associated with different nutrient partitioning phenotypes.

A 50-kb deletion disrupting the RSPO2 gene is associated with tetradysmelia in Holstein Friesian cattle.

BACKGROUND: Tetradysmelia is a rare genetic disorder that is characterized by an extremely severe reduction of all limb parts distal of the scapula and pelvic girdle. We studied a Holstein Friesian backcross family with 24 offspring, among which six calves displayed autosomal recessive tetradysmelia. In order to identify the genetic basis of the disorder, we genotyped three affected calves, five dams and nine unaffected siblings using a Bovine Illumina 50 k BeadChip and sequenced the whole genome of the sire. RESULTS: Pathological examination of four tetradysmelia cases revealed a uniform and severe dysmelia of all limbs. Applying a homozygosity mapping approach, we identified a homozygous region of 10.54 Mb on chromosome 14 (Bos taurus BTA14). Only calves that were diagnosed with tetradysmelia shared a distinct homozygous haplotype for this region. We sequenced the whole genome of the cases' sire and searched for heterozygous single nucleotide polymorphisms (SNPs) and small variants on BTA14 that were uniquely present in the sire and absent from 3102 control whole-genome sequences of the 1000 Bull Genomes Project, but none were identified in the 10.54-Mb candidate region on BTA14. Therefore, we subsequently performed a more comprehensive analysis by also considering structural variants and detected a 50-kb deletion in the targeted chromosomal region that was in the heterozygous state in the cases' sire. Using PCR, we confirmed that this detected deletion segregated perfectly within the family with tetradysmelia. The deletion spanned three exons of the bovine R-spondin 2 (RSPO2) gene, which encode three domains of the respective protein. R-spondin 2 is a secreted ligand of leucine-rich repeats containing G protein-coupled receptors that enhance Wnt signalling and is involved in a broad range of developmental processes during embryogenesis. CONCLUSIONS: We identified a 50-kb deletion on BTA14 that disrupts the coding sequence of the RSPO2 gene and is associated with bovine tetradysmelia. To our knowledge, this is the first reported candidate causal mutation for tetradysmelia in a large animal model. Since signalling pathways involved in limb development are conserved across species, the observed inherited defect may serve as a model to further elucidate fundamental pathways of limb development.

Development and evaluation of a milk protein transcript depletion method for differential transcriptome analysis in mammary gland tissue.

BACKGROUND: In the mammary gland transcriptome of lactating dairy cows genes encoding milk proteins are highly abundant, which can impair the detection of lowly expressed transcripts and can bias the outcome in global transcriptome analyses. Therefore, the aim of this study was to develop and evaluate a method to deplete extremely highly expressed transcripts in mRNA from lactating mammary gland tissue. RESULTS: Selective RNA depletion was performed by hybridization of antisense oligonucleotides targeting genes encoding the caseins (CSN1S1, CSN1S2, CSN2 and CSN3) and whey proteins (LALBA and PAEP) within total RNA followed by RNase H-mediated elimination of the respective transcripts. The effect of the RNA depletion procedure was monitored by RNA sequencing analysis comparing depleted and non-depleted RNA samples from Escherichia coli (E. coli) challenged and non-challenged udder tissue of lactating cows in a proof of principle experiment. Using RNase H-mediated RNA depletion, the ratio of highly abundant milk protein gene transcripts was reduced in all depleted samples by an average of more than 50% compared to the non-depleted samples. Furthermore, the sensitivity for discovering transcripts with marginal expression levels and transcripts not yet annotated was improved. Finally, the sensitivity to detect significantly differentially expressed transcripts between non-challenged and challenged udder tissue was increased without leading to an inadvertent bias in the pathogen challenge-associated biological signaling pathway patterns. CONCLUSIONS: The implementation of selective RNase H-mediated RNA depletion of milk protein gene transcripts from the mammary gland transcriptome of lactating cows will be highly beneficial to establish comprehensive transcript catalogues of the tissue that better reflects its transcriptome complexity.

Unmapped reads from cattle RNAseq data: A source for missing and misassembled sequences in the reference assemblies and for detection of pathogens in the host.

Usually, reads from transcriptome sequencing data unmapped to the target species' reference genome are disregarded. A recent RNAseq project on the new fatal disease Bovine Neonatal Pancytopenia had indicated an unexplained immune response signature to a double-stranded RNA virus. To unravel its background, contigs were de novo assembled from unmapped RNAseq reads and aligned against the bovine genome assemblies and multispecies NCBI databases. Lack of genuine virus sequence contigs rejected the hypothesis of a live virus being causal for the unexplained immune response. Alignment data also demonstrated incomplete bovine reference genome assemblies. In addition, we found that several parasite and virus genome reference assemblies in NCBI were contaminated with bovine DNA and confirmed recombination of bovine DNA into BVD virus strains. Exploring unmapped reads can extract useful biological information regarding the presence of microorganisms and can highlight issues with reference genome assemblies of host and pathogen species.

A novel RNAseq-assisted method for MHC class I genotyping in a non-model species applied to a lethal vaccination-induced alloimmune disease.

BACKGROUND: MHC class I genotyping is essential for a wide range of biomedical, immunological and biodiversity applications. Whereas in human a comprehensive MHC class I allele catalogue is available, respective data in non-model species is scarce in spite of decades of research. RESULTS: Taking advantage of the new high-throughput RNA sequencing technology (RNAseq), we developed a novel RNAseq-assisted method (RAMHCIT) for MHC class I typing at nucleotide level. RAMHCIT is performed on white blood cells, which highly express MHC class I molecules enabling reliable discovery of new alleles and discrimination of closely related alleles due to the high coverage of alleles with reads. RAMHCIT is more comprehensive than previous methods, because no targeted PCR pre-amplification of MHC loci is necessary, which avoids preselection of alleles as usually encountered, when amplification with MHC class I primers is performed prior to sequencing. In addition to allele identification, RAMHCIT also enables quantification of MHC class I expression at allele level, which was remarkably consistent across individuals. CONCLUSIONS: Successful application of RAMHCIT is demonstrated on a data set from cattle with different phenotype regarding a lethal, vaccination-induced alloimmune disease (bovine neonatal pancytopenia), for which MHC class I alleles had been postulated as causal agents.

Genomic prediction of unordered categorical traits: an application to subpopulation assignment in German Warmblood horses.

BACKGROUND: Categorical traits without ordinal representation of classes do not qualify for threshold models. Alternatively, the multinomial problem can be assessed by a sequence of independent binary contrasts using schemes such as one-vs-all or one-vs-one. Class probabilities can be arrived at by normalization or pair-wise coupling strategies. We assessed the predictive ability of whole-genome regression models and support vector machines for the classification of horses into four German Warmblood breeds. RESULTS: Prediction accuracies of leave-one-out cross-validation were high and ranged from 0.75 to 0.97 depending on the binary classifier and breeds incorporated in the training. An analysis of the population structure using eigenvectors of the genomic relationship matrix revealed clustering of individuals beyond the given breed labels. Admixture between two breeds became apparent which had substantial impact on the prediction accuracies between those two breeds and also influenced the contrasts between other breeds. CONCLUSIONS: Genomic prediction of unordered categorical traits was successfully applied to subpopulation assignment of German Warmblood horses. The applied methodology is a straightforward extension of existing binary threshold models for genomic prediction.

Detection of genetic variants affecting cattle behaviour and their impact on milk production: a genome-wide association study.

Behaviour traits of cattle have been reported to affect important production traits, such as meat quality and milk performance as well as reproduction and health. Genetic predisposition is, together with environmental stimuli, undoubtedly involved in the development of behaviour phenotypes. Underlying molecular mechanisms affecting behaviour in general and behaviour and productions traits in particular still have to be studied in detail. Therefore, we performed a genome-wide association study in an F2 Charolais × German Holstein cross-breed population to identify genetic variants that affect behaviour-related traits assessed in an open-field and novel-object test and analysed their putative impact on milk performance. Of 37,201 tested single nucleotide polymorphism (SNPs), four showed a genome-wide and 37 a chromosome-wide significant association with behaviour traits assessed in both tests. Nine of the SNPs that were associated with behaviour traits likewise showed a nominal significant association with milk performance traits. On chromosomes 14 and 29, six SNPs were identified to be associated with exploratory behaviour and inactivity during the novel-object test as well as with milk yield traits. Least squares means for behaviour and milk performance traits for these SNPs revealed that genotypes associated with higher inactivity and less exploratory behaviour promote higher milk yields. Whether these results are due to molecular mechanisms simultaneously affecting behaviour and milk performance or due to a behaviour predisposition, which causes indirect effects on milk performance by influencing individual reactivity, needs further investigation.

Establishment of a cloning-free CRISPR/Cas9 protocol to generate large deletions in the bovine MDBK cell line.

The CRISPR/Cas9 technique applied to modify the cattle genome has value in increasing animal health and welfare. Here, we established a simple, fast, and efficient cloning-free CRISPR/Cas9 protocol for large deletions of genomic loci in the frequently used model bovine MDBK cell line. The main advantages of our protocol are as follows: (i) pre-screening of the sgRNA efficiency with a fast and simple cleavage assay, (ii) reliable detection of genomic edits primarily by PCR and confirmed by DNA sequencing, and (iii) single cell sorting with FACS providing specific genetic information from modified cells of interest. Therefore, our method could be successfully applied in different studies, including functional validation of any genetic or regulatory elements.

Cryptosporidium parvum infection alters the intestinal mucosa transcriptome in neonatal calves: implications for immune function.

One of the leading causes of infectious diarrhea in newborn calves is the apicomplexan protozoan Cryptosporidium parvum (C. parvum). However, little is known about its immunopathogenesis. Using next generation sequencing, this study investigated the immune transcriptional response to C. parvum infection in neonatal calves. Neonatal male Holstein-Friesian calves were either orally infected (N = 5) or not (CTRL group, N = 5) with C. parvum oocysts (gp60 subtype IIaA15G2R1) at day 1 of life and slaughtered on day 7 after infection. Total RNA was extracted from the jejunal mucosa for short read. Differentially expressed genes (DEGs) between infected and CTRL groups were assessed using DESeq2 at a false discovery rate < 0.05. Infection did not affect plasma immunohematological parameters, including neutrophil, lymphocyte, monocyte, leucocyte, thrombocyte, and erythrocyte counts as well as hematocrit and hemoglobin concentration on day 7 post infection. The immune-related DEGs were selected according to the UniProt immune system process database and were used for gene ontology (GO) and pathway enrichment analysis using Cytoscape (v3.9.1). Based on GO analysis, DEGs annotated to mucosal immunity, recognizing and presenting antigens, chemotaxis of neutrophils, eosinophils, natural killer cells, B and T cells mediated by signaling pathways including toll like receptors, interleukins, tumor necrosis factor, T cell receptor, and NF-KB were upregulated, while markers of macrophages chemotaxis and cytosolic pattern recognition were downregulated. This study provides a holistic snapshot of immune-related pathways induced by C. parvum in calves, including novel and detailed feedback and feedforward regulatory mechanisms establishing the crosstalk between innate and adaptive immune response in neonate calves, which could be utilized further to develop new therapeutic strategies.

A selective block of nuclear actin export stabilizes the giant nuclei of Xenopus oocytes.

Actin is a major cytoskeletal element and is normally kept cytoplasmic by exportin 6 (Exp6)-driven nuclear export. Here, we show that Exp6 recognizes actin features that are conserved from yeast to human. Surprisingly however, microinjected actin was not exported from Xenopus laevis oocyte nuclei, unless Exp6 was co-injected, indicating that the pathway is inactive in this cell type. Indeed, Exp6 is undetectable in oocytes, but is synthesized from meiotic maturation onwards, which explains how actin export resumes later in embryogenesis. Exp6 thus represents the first example of a strictly developmentally regulated nuclear transport pathway. We asked why Xenopus oocytes lack Exp6 and observed that ectopic application of Exp6 renders the giant oocyte nuclei extremely fragile. This effect correlates with the selective disappearance of a sponge-like intranuclear scaffold of F-actin. These nuclei have a normal G2-phase DNA content in a volume 100,000 times larger than nuclei of somatic cells. Apparently, their mechanical integrity cannot be maintained by chromatin and the associated nuclear matrix, but instead requires an intranuclear actin-scaffold.

Monitoring the immune response to vaccination with an inactivated vaccine associated to bovine neonatal pancytopenia by deep sequencing transcriptome analysis in cattle.

Bovine neonatal pancytopenia (BNP) is a new fatal, alloimmune/alloantibody mediated disease of new-born calves induced by ingestion of colostrum from cows, which had been vaccinated with a specific vaccine against the Bovine Virus Diarrhoea Virus (BVDV). The hypothesis of pathogenic MHC class I molecules in the vaccine had been put up, but no formal proof of specific causal MHC class I alleles has been provided yet. However, the unique features of the vaccine obviously result in extremely high specific antibody titres in the vaccinated animals, but apparently also in further molecules inducing BNP. Thus, a comprehensive picture of the immune response to the vaccine is essential. Applying the novel approach of next generation RNA sequencing (RNAseq), our study provides a new holistic, comprehensive analysis of the blood transcriptome regulation after vaccination with the specific BVDV vaccine. Our RNAseq approach identified a novel cytokine-like gene in the bovine genome that is highly upregulated after vaccination. This gene has never been described before in any other species and might be specific to ruminant immune response. Furthermore, our data revealed a very coordinated immune response to double-stranded (ds) RNA or a dsRNA analogue after vaccination with the inactivated single-stranded (ss) RNA vaccine. This would suggest either a substantial contamination of the vaccine with dsRNA from host cells after virus culture or a dsRNA analogue applied to the vaccine. The first option would highlight the potential risks associated with virus culture on homologous cells during vaccine production; the latter option would emphasise the potential risks associated with immune stimulating adjuvants used in vaccine production.

Improving the annotation of the cattle genome by annotating transcription start sites in a diverse set of tissues and populations using Cap Analysis Gene Expression sequencing.

Understanding the genomic control of tissue-specific gene expression and regulation can help to inform the application of genomic technologies in farm animal breeding programs. The fine mapping of promoters [transcription start sites (TSS)] and enhancers (divergent amplifying segments of the genome local to TSS) in different populations of cattle across a wide diversity of tissues provides information to locate and understand the genomic drivers of breed- and tissue-specific characteristics. To this aim, we used Cap Analysis Gene Expression (CAGE) sequencing, of 24 different tissues from 3 populations of cattle, to define TSS and their coexpressed short-range enhancers (<1 kb) in the ARS-UCD1.2_Btau5.0.1Y reference genome (1000bulls run9) and analyzed tissue and population specificity of expressed promoters. We identified 51,295 TSS and 2,328 TSS-Enhancer regions shared across the 3 populations (dairy, beef-dairy cross, and Canadian Kinsella composite cattle from 2 individuals, 1 of each sex, per population). Cross-species comparative analysis of CAGE data from 7 other species, including sheep, revealed a set of TSS and TSS-Enhancers that were specific to cattle. The CAGE data set will be combined with other transcriptomic information for the same tissues to create a new high-resolution map of transcript diversity across tissues and populations in cattle for the BovReg project. Here we provide the CAGE data set and annotation tracks for TSS and TSS-Enhancers in the cattle genome. This new annotation information will improve our understanding of the drivers of gene expression and regulation in cattle and help to inform the application of genomic technologies in breeding programs.

From the comparative study of a circRNA originating from an mammalian ATXN2L intron to understanding the genesis of intron lariat-derived circRNAs.

Circular intronic RNAs (ciRNAs) are still unexplored regarding mechanisms for their emergence. We considered the ATXN2L intron lariat-derived circular RNA (ciRNA-ATXN2L) as an opportunity to conduct a cross-species examination of ciRNA genesis. To this end, we investigated 207 datasets from 4 tissues and from 13 mammalian species. While in eight species, ciRNA-ATXN2L was never detected, in pigs and rabbits, ciRNA-ATXN2L was expressed in all tissues and sometimes at very high levels. Bovine tissues were an intermediate case and in macaques and cats, only ciRNA-ATXN2L traces were detected. The pattern of ciRNA-ATXN2L restricted to only five species is not related to a particular evolution of intronic sequences. To empower our analysis, we considered 221 additional introns including 80 introns where a lariat-derived ciRNA was previously described. The primary driver of micro-ciRNA genesis (< 155 nt as ciRNA-ATXN2L) appears to be the absence of a canonical "A" (i.e. a "tnA" located in the usual branching region) to build the lariat around this adenosine. The balance between available "non canonical-A" (no ciRNA genesis) and "non-A" (ciRNA genesis) for use as a branch point to build the lariat could modify the expression level of ciRNA-ATXN2L. In addition, the rare localization of the 2'-5' bond in an open RNA secondary structure could also negatively affect the lifetime of ciRNAs (macaque ciRNA-ATXN2L). Our analyses suggest that ciRNA-ATXN2L is likely a functionless splice remnant. This study provides a better understanding of the ciRNAs origin, especially drivers for micro ciRNA genesis.

Replacement of microsatellite markers by imputed medium-density SNP arrays for parentage control in German warmblood horses.

In horses, parentage control is currently performed based on an internationally standardized panel of 17 microsatellite (MS) markers comprising 12 mandatory and five optional markers. Unlike MS, single nucleotide polymorphism (SNP) profiles support a wider portfolio of genomic applications, including parentage control. A transition to SNP-based parentage control is favorable, but requires additional efforts for ensuring generation-overlapping availability of marker genotypes of the same type. To avoid double genotyping of either parents or offspring for changing to SNP technology and enable efficient transition, we tested whether MS genotypes used for parentage control could be reliably imputed from a medium-density SNP panel in German warmblood horses. Imputation accuracy was tested in a tenfold cross-validation with two approaches: within breed (option A) and across breeds (option B). Average imputation accuracies of 97.98% (A) and 96.17% (B) were achieved, respectively. Due to interbreed differences in genotyping rates, five MS markers of low genotyping rate (GTR; < 90%) could be imputed with higher accuracy within breed (98.18%) than across breeds (90.73%). MS markers with high GTR performed homogeneously well in option B (98.44%) and showed slightly lower accuracy in option A (97.90%). Among these markers, AHT5 proved to be problematic for imputation regardless of the approach, revealing accuracies of 86.40% (A) and 88.70% (B). Better results for MS markers with high GTR and savings in computational processing justified the choice of option B for routine implementation. To date, more than 9500 horses have undergone the new parentage control based on imputed MS genotypes.