Dominique Rocha

The myostatin protein is a potent negative regulator of skeletal muscle growth encoded by the MSTN gene. MSTN loss‐of‐function variants lead to a particular cattle phenotype characterized by an increase in skeletal muscle mass, known as “double muscling” or “double muscled”. However, most of the MSTN causal variants that have been linked to this phenotype lack experimental validation. This is the case, for example, for the five missense MSTN variants reported to be causal according to the Online Mendelian Inheritance in Animals. RNA splicing plays a major role in regulating gene expression; therefore, exploring the effects of variants on RNA splicing may provide relevant information on their functional impact. Here, we have set up a full‐length gene assay (FLGA) to functionally assess MSTN splicing variants, and we have used it to test the five missense variants plus a well‐described deep intronic splicing variant as a positive control. We also evaluated the performances of SpliceAI and Pangolin, two deep learning‐based splice predictors, to identify potential splicing effects of these six variants. Our FLGA system performed well and showed that none of the missense variants has an effect on splicing, unlike the positive control. For each variant, splicing program predictions were perfectly concordant with the effect observed in the FLGA. We have produced a relevant and powerful assay to analyze MSTN splicing variants in cattle. SpliceAI and Pangolin may be efficiently used to screen large datasets of MSTN variants and sort the best candidates prior to experimental validation using an FLGA.

Marbling is one of the most important beef quality traits. An association between a non‐synonymous variant located in EGFLAM (EGF‐like fibronectin and laminin G) and marbling in Hanwoo cattle has recently been published. We therefore investigated the association between this SNP (rs109436056 SNP) and marbling in Limousin cattle. A total of 355 animals were phenotyped for marbling and genotyped for this SNP. Significant association (p < 0.05) was observed, in which genotype CC exhibited higher marbling. This SNP could be used for the genetic improvement of marbling in Limousin cattle.

Standardised analysis pipelines are an important part of FAIR bioinformatics research. Over the last decade, there has been a notable shift from point-and-click pipeline solutions such as Galaxy towards command-line solutions such as Nextflow and Snakemake. We report on recent developments in the nf-core and Nextflow frameworks that have led to widespread adoption across many scientific communities. We describe how adopting nf-core standards enables faster development, improved interoperability, and collaboration with the >8,000 members of the nf-core community. The recent development of Nextflow Domain-Specific Language 2 (DSL2) allows pipeline components to be shared and combined across projects. The nf-core community has harnessed this with a library of modules and subworkflows that can be integrated into any Nextflow pipeline, enabling research communities to progressively transition to nf-core best practices. We present a case study of nf-core adoption by six European research consortia, grouped under the EuroFAANG umbrella and dedicated to farmed animal genomics. We believe that the process outlined in this report can inspire many large consortia to seek harmonisation of their data analysis procedures.

Mitochondrial DNA sequences are frequently transferred into the nuclear genome, generating nuclear mitochondrial DNA sequences (NUMTs). Here, we analysed, for the first time, NUMTs in the domestic yak genome. We obtained 499 alignment matches covering 340.2 kbp of the yak nuclear genome. After a merging step, we identified 167 NUMT regions with a total length of ~ 503 kbp, representing 0.02% of the nuclear genome. We discovered copies of all mitochondrial regions and found that most NUMT regions are intergenic or intronic and mostly untranscribed. 98 different NUMT regions from domestic yak showed high homology with cow and/or wild yak genomes, suggesting selection or hybridization between domestic/wild yak and cow. To rule out the possibility that the identified NUMTs could be artifacts of the domestic yak genome assembly, we validated experimentally five NUMT regions by PCR amplification. As NUMT regions show high similarity to the mitochondrial genome can potentially pose a risk to domestic yak DNA mitochondrial studies, special care is therefore needed to select primers for PCR amplification of mitochondrial DNA sequences.

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.

Chicken is a valuable model for understanding fundamental biology and vertebrate evolution, as well as a major source of nutrient-dense and lean-protein-enriched food globally. Although it is the first non-mammalian amniote genome to be sequenced, the chicken genome still lacks a systematic characterization of functional variation. Here, through integrating bulk RNA-Seq of 7,015 samples, single-cell RNA-Seq of 127,598 cells and 2,869 whole-genome sequences, we presented an atlas of regulatory variants across 28 chicken tissues and 43 cell types, including millions of regulatory effects on primary expression (protein-coding genes, lncRNA and exon) and post-transcriptional modifications (alternative splicing and 3’ untranslated region alternative polyadenylation). We highlighted distinct molecular mechanisms underlying these regulatory variants and their context dependence and utility in interpreting genome-wide associations of 108 chicken complex traits. Finally, we illustrated lineage-specific features of gene regulation between chickens and mammals, and demonstrated how this resource can assist with cross-species gene mapping of complex traits.

Chicken is a valuable model for understanding fundamental biology, vertebrate evolution and diseases, as well as a major source of nutrient-dense and lean-protein-enriched food globally. Although it is the first non-mammalian amniote genome to be sequenced, the chicken genome still lacks a systematic characterization of functional impacts of genetic variants. Here, through integrating 7,015 RNA-Seq and 2,869 whole-genome sequence data, the Chicken Genotype- Tissue Expression (ChickenGTEx) project presents the pilot reference of regulatory variants in 28 chicken tissue transcriptomes, including millions of regulatory effects on primary expression (including protein-coding genes, lncRNA and exon) and post-transcriptional modifications (alternative splicing and 3’ untranslated region alternative polyadenylation). We explored the tissue-sharing and context-specificity of these regulatory variants, their underlying molecular mechanisms of action, and their utility in interpreting adaptation and genome-wide associations of 108 chicken complex traits. Finally, we illustrated shared and lineage-specific features of gene regulation between chickens and mammals, and demonstrated how the ChickenGTEx resource can further assist with translating genetic findings across species. One-Sentence Summary: The ChickenGTEx provides a multi-tissue reference of regulatory variants for chicken genetics and genomics, functional genomics, precision breeding, veterinary medicine, vertebrate evolution and even human biomedicine.

DGAT1 is playing a major role in fat metabolism and triacylglyceride synthesis. Only two DGAT1 loss-of-function variants altering milk production traits in cattle have been reported to date, namely p.M435L and p.K232A. The p.M435L variant is a rare alteration and has been associated with skipping of exon 16 which results in a non-functional truncated protein, and the p.K232A-containing haplotype has been associated with modifications of the splicing rate of several DGAT1 introns. In particular, the direct causality of the p.K232A variant in decreasing the splicing rate of the intron 7 junction was validated using a minigene assay in MAC-T cells. As both these DGAT1 variants were shown to be spliceogenic, we developed a full-length gene assay (FLGA) to re-analyse p.M435L and p.K232A variants in HEK293T and MAC-T cells. Qualitative RT-PCR analysis of cells transfected with the full-length DGAT1 expression construct carrying the p.M435L variant highlighted complete skipping of exon 16. The same analysis performed using the construct carrying the p.K232A variant showed moderate differences compared to the wild-type construct, suggesting a possible effect of this variant on the splicing of intron 7. Finally, quantitative RT-PCR analyses of cells transfected with the p.K232A-carrying construct did not show any significant modification on the splicing rate of introns 1, 2 and 7. In conclusion, the DGAT1 FLGA confirmed the p.M435L impact previously observed in vivo, but invalidated the hypothesis whereby the p.K232A variant strongly decreased the splicing rate of intron 7.