Adaptative response to changes in pyruvate metabolism on the epigenetic landscapes and transcriptomics of bovine embryos.

The epigenetic reprogramming that occurs during the earliest stages of embryonic development has been described as crucial for the initial events of cell specification and differentiation. Recently, the metabolic status of the embryo has gained attention as one of the main factors coordinating epigenetic events. In this work, we investigate the link between pyruvate metabolism and epigenetic regulation by culturing bovine embryos from day 5 in the presence of dichloroacetate (DCA), a pyruvate analog that increases the pyruvate to acetyl-CoA conversion, and iodoacetate (IA), which inhibits the glyceraldehyde-3-phosphate dehydrogenase (GAPDH), leading to glycolysis inhibition. After 8 h of incubation, both DCA and IA-derived embryos presented higher mitochondrial membrane potential. Nevertheless, in both cases, lower levels of acetyl-CoA, ATP-citrate lyase and mitochondrial membrane potential were found in blastocysts, suggesting an adaptative metabolic response, especially in the DCA group. The metabolic alteration found in blastocysts led to changes in the global pattern of H3K9 and H3K27 acetylation and H3K27 trimethylation. Transcriptome analysis revealed that such alterations resulted in molecular differences mainly associated to metabolic processes, establishment of epigenetic marks, control of gene expression and cell cycle. The latter was further confirmed by the alteration of total cell number and cell differentiation in both groups when compared to the control. These results corroborate previous evidence of the relationship between the energy metabolism and the epigenetic reprogramming in preimplantation bovine embryos, reinforcing that the culture system is decisive for precise epigenetic reprogramming, with consequences for the molecular control and differentiation of cells.

Transcriptome annotation of 17 porcine tissues using nanopore sequencing technology.

The annotation of animal genomes plays an important role in elucidating molecular mechanisms behind the genetic control of economically important traits. Here, we employed long-read sequencing technology, Oxford Nanopore Technology, to annotate the pig transcriptome across 17 tissues from two Yorkshire littermate pigs. More than 9.8 million reads were obtained from a single flow cell, and 69 781 unique transcripts at 50 108 loci were identified. Of these transcripts, 16 255 were found to be novel isoforms, and 22 344 were found at loci that were novel and unannotated in the Ensembl (release 102) and NCBI (release 106) annotations. Novel transcripts were mostly expressed in cerebellum, followed by lung, liver, spleen, and hypothalamus. By comparing the unannotated transcripts to existing databases, there were 21 285 (95.3%) transcripts matched to the NT database (v5) and 13 676 (61.2%) matched to the NR database (v5). Moreover, there were 4324 (19.4%) transcripts matched to the SwissProt database (v5), corresponding to 11 356 proteins. Tissue-specific gene expression analyses showed that 9749 transcripts were highly tissue-specific, and cerebellum contained the most tissue-specific transcripts. As the same samples were used for the annotation of cis-regulatory elements in the pig genome, the transcriptome annotation generated by this study provides an additional and complementary annotation resource for the Functional Annotation of Animal Genomes effort to comprehensively annotate the pig genome.

Prediction of transcript isoforms in 19 chicken tissues by Oxford Nanopore long-read sequencing.

To identify and annotate transcript isoforms in the chicken genome, we generated Nanopore long-read sequencing data from 68 samples that encompassed 19 diverse tissues collected from experimental adult male and female White Leghorn chickens. More than 23.8 million reads with mean read length of 790 bases and average quality of 18.2 were generated. The annotation and subsequent filtering resulted in the identification of 55,382 transcripts at 40,547 loci with mean length of 1,700 bases. We predicted 30,967 coding transcripts at 19,461 loci, and 16,495 lncRNA transcripts at 15,512 loci. Compared to existing reference annotations, we found ∼52% of annotated transcripts could be partially or fully matched while ∼47% were novel. Seventy percent of novel transcripts were potentially transcribed from lncRNA loci. Based on our annotation, we quantified transcript expression across tissues and found two brain tissues (i.e., cerebellum and cortex) expressed the highest number of transcripts and loci. Furthermore, ∼22% of the transcripts displayed tissue specificity with the reproductive tissues (i.e., testis and ovary) exhibiting the most tissue-specific transcripts. Despite our wide sampling, ∼20% of Ensembl reference loci were not detected. This suggests that deeper sequencing and additional samples that include different breeds, cell types, developmental stages, and physiological conditions, are needed to fully annotate the chicken genome. The application of Nanopore sequencing in this study demonstrates the usefulness of long-read data in discovering additional novel loci (e.g., lncRNA loci) and resolving complex transcripts (e.g., the longest transcript for the TTN locus).

Large-Scale Multiplexing Permits Full-Length Transcriptome Annotation of 32 Bovine Tissues From a Single Nanopore Flow Cell.

A comprehensive annotation of transcript isoforms in domesticated species is lacking. Especially considering that transcriptome complexity and splicing patterns are not well-conserved between species, this presents a substantial obstacle to genomic selection programs that seek to improve production, disease resistance, and reproduction. Recent advances in long-read sequencing technology have made it possible to directly extrapolate the structure of full-length transcripts without the need for transcript reconstruction. In this study, we demonstrate the power of long-read sequencing for transcriptome annotation by coupling Oxford Nanopore Technology (ONT) with large-scale multiplexing of 93 samples, comprising 32 tissues collected from adult male and female Hereford cattle. More than 30 million uniquely mapping full-length reads were obtained from a single ONT flow cell, and used to identify and characterize the expression dynamics of 99,044 transcript isoforms at 31,824 loci. Of these predicted transcripts, 21% exactly matched a reference transcript, and 61% were novel isoforms of reference genes, substantially increasing the ratio of transcript variants per gene, and suggesting that the complexity of the bovine transcriptome is comparable to that in humans. Over 7,000 transcript isoforms were extremely tissue-specific, and 61% of these were attributed to testis, which exhibited the most complex transcriptome of all interrogated tissues. Despite profiling over 30 tissues, transcription was only detected at about 60% of reference loci. Consequently, additional studies will be necessary to continue characterizing the bovine transcriptome in additional cell types, developmental stages, and physiological conditions. However, by here demonstrating the power of ONT sequencing coupled with large-scale multiplexing, the task of exhaustively annotating the bovine transcriptome - or any mammalian transcriptome - appears significantly more feasible.

Sperm morphometry is affected by increased inbreeding in the Retinta cattle breed: A molecular approach.

The effect of inbreeding depression on sperm motility is well documented, but its influence on sperm morphometry has been scarcely examined to date. Here, we combined the use of computer-assisted sperm morphometry analysis (CASMA) with a SNP-based genomic approach to determine and characterize the effect of inbreeding on the sperm shape of a highly inbred cattle population. We determined seven morphometric parameters on frozen-thawed sperm samples of 57 Retinta bulls: length (L, µm), width (W, µm), area (A, µm2 ), perimeter (P, µm), ellipticity (ELI; L/W), elongation (L-W)/(L + W) and perimeter-to-area shape factor (p2a; P2 /4 × π × A). The comparison of highly inbred (HI) and lowly inbreed (LI) individuals based on runs of homozygosity (ROH) inbreeding values (F ROH ) showed no differences between groups. An additional two-step unsupervised sperm subpopulation analysis based on morphometric parameters showed significant differences in the abundance of different sperm subpopulations between groups (p < 0.05). This analysis revealed that HI bulls harbored a higher percentage of narrow-head sperm as opposed to the higher percentage of large- and round-headed sperm detected in LI. A further genomic characterization revealed 23 regions differentially affected by inbreeding in both groups, detecting six genes (SPAG6, ARMC3, PARK7, VAMP3, DYNLRB2, and PHF7) previously related to different spermatogenesis-associated processes.

Transcription initiation mapping in 31 bovine tissues reveals complex promoter activity, pervasive transcription, and tissue-specific promoter usage.

Characterizing transcription start sites is essential for understanding the regulatory mechanisms that control gene expression. Recently, a new bovine genome assembly (ARS-UCD1.2) with high continuity, accuracy, and completeness was released; however, the functional annotation of the bovine genome lacks precise transcription start sites and contains a low number of transcripts in comparison to human and mouse. By using the RAMPAGE approach, this study identified transcription start sites at high resolution in a large collection of bovine tissues. We found several known and novel transcription start sites attributed to promoters of protein-coding and lncRNA genes that were validated through experimental and in silico evidence. With these findings, the annotation of transcription start sites in cattle reached a level comparable to the mouse and human genome annotations. In addition, we identified and characterized transcription start sites for antisense transcripts derived from bidirectional promoters, potential lncRNAs, mRNAs, and pre-miRNAs. We also analyzed the quantitative aspects of RAMPAGE to produce a promoter activity atlas, reaching highly reproducible results comparable to traditional RNA-seq. Coexpression networks revealed considerable use of tissue-specific promoters, especially between brain and testicle, which expressed several genes in common from alternate loci. Furthermore, regions surrounding coexpressed modules were enriched in binding factor motifs representative of each tissue. The comprehensive annotation of promoters in such a large collection of tissues will substantially contribute to our understanding of gene expression in cattle and other mammalian species, shortening the gap between genotypes and phenotypes.

Gametes from stem cells: Status and applications in animal reproduction.

In vitro gamete differentiation could revolutionize animal production by decreasing generation intervals, increasing the number of gametes per animal and facilitating the dissemination of elite genetics. In addition, it could help to develop new strategies for the conservation of endangered species. The recent in vitro reconstitution of germ cell development in mice has inspired researchers to invest their best efforts into reproducing this achievement in livestock species. With this goal in mind, multiple differentiation approaches and cell sources have been evaluated. The degree of success in these evaluations varies according to the species and the stage of development studied, but, in general, partially positive results have been obtained. Evidence suggests that although functional gametes with true reproductive potential are still to be obtained, it is a matter of time before this goal is achieved.

In vitro breeding: application of embryonic stem cells to animal production†.

Embryonic stem cells (ESCs) are derived from the inner cell mass of preimplantation blastocysts. For decades, attempts to efficiently derive ESCs in animal livestock species have been unsuccessful, but this goal has recently been achieved in cattle. Together with the recent reconstitution of the germ cell differentiation processes from ESCs in mice, these achievements open new avenues for the development of promising technologies oriented toward improving health, animal production, and the environment. In this article, we present a strategy that will notably accelerate genetic improvement in livestock populations by reducing the generational interval, namely in vitro breeding (IVB). IVB combines genomic selection, a widely used strategy for genetically improving livestock, with ESC derivation and in vitro differentiation of germ cells from pluripotent stem cells. We also review the most recent findings in the fields on which IVB is based. Evidence suggests this strategy will be soon within reach.

Runs of homozygosity in a selected cattle population with extremely inbred bulls: Descriptive and functional analyses revealed highly variable patterns.

The analysis of runs of homozygosity (ROH), using high throughput genomic data, has become a valuable and frequently used methodology to characterize the genomic and inbreeding variation of livestock and wildlife animal populations. However, this methodology has been scarcely used in highly inbred domestic animals. Here, we analyzed and characterized the occurrence of ROH fragments in highly inbred (HI; average pedigree-based inbreeding coefficient FPED = 0.164; 0.103 to 0.306) and outbred Retinta bulls (LI; average FPED = 0.008; 0 to 0.025). We studied the length of the fragments, their abundance, and genome distribution using high-density microarray data. The number of ROH was significantly higher in the HI group, especially for long fragments (>8Mb). In the LI group, the number of ROH continuously decreased with fragment length. Genome-wide distribution of ROH was highly variable between samples. Some chromosomes presented a larger number of fragments (BTA1, BTA19, BTA29), others had longer fragments (BTA4, BTA12, BTA17), while other ones showed an increased ROH accumulation over specific loci (BTA2, BTA7, BTA23, BTA29). Similar differences were observed in the analysis of 12 individuals produced by a similar inbred event (FPED3 = 0.125). The correlation between the fraction of the genome covered by ROH (FROH) and FPED was high (0.79), suggesting that ROH-based estimations are indicative of inbreeding levels. On the other hand, the correlation between FPED and the microsatellite-based inbreeding coefficient (FMIC) was only moderate (r = 0.44), suggesting that STR-based inbreeding estimations should be avoided. Similarly, we found a very low correlation (r = -0.0132) between recombination rate and ROH abundance across the genome. Finally, we performed functional annotation analyses of genome regions with significantly enriched ROH abundance. Results revealed gene clusters related to pregnancy-associated proteins and immune reaction. The same analysis performed for regions enriched with recently formed ROH (> 8 Mb) showed gene clusters related to flagellum assembly. In both cases, the processes were related to male and female reproductive functions, which may partially explain the reduced fertility associated with inbred populations.

Genetic Variation in FABP4 and Evaluation of Its Effects on Beef Cattle Fat Content.

FABP4 is a protein primarily expressed in adipocytes and macrophages that plays a key role in fatty acid trafficking and lipid hydrolysis. FABP4 gene polymorphisms have been associated with meat quality traits in cattle, mostly in Asian breeds under feedlot conditions. The objectives of this work were to characterize FABP4 genetic variation in several worldwide cattle breeds and evaluate possible genotype effects on fat content in a pasture-fed crossbred (Angus-Hereford-Limousin) population. We re-sequenced 43 unrelated animals from nine cattle breeds (Angus, Brahman, Creole, Hereford, Holstein, Limousin, Nelore, Shorthorn, and Wagyu) and obtained 22 single nucleotide polymorphisms (SNPs) over 3,164 bp, including four novel polymorphisms. Haplotypes and linkage disequilibrium analyses showed a high variability. Five SNPs were selected to perform validation and association studies in our crossbred population. Four SNPs showed well-balanced allele frequencies (minor frequency > 0.159), and three showed no significant deviations from Hardy-Weinberg proportions. SNPs showed significant effects on backfat thickness and fatty acid composition (P < 0.05). The protein structure of one of the missense SNPs was analyzed to elucidate its possible effect on fat content in our studied population. Our results revealed a possible blockage of the fatty acid binding site by the missense mutation.

Genetic characterisation of PPARG, CEBPA and RXRA, and their influence on meat quality traits in cattle.

BACKGROUND: Peroxisome proliferator-activated receptor gamma (PPARG), CCAAT/enhancer binding protein alpha (CEBPA) and retinoid X receptor alpha (RXRA) are nuclear transcription factors that play important roles in regulation of adipogenesis and fat deposition. The objectives of this study were to characterise the variability of these three candidate genes in a mixed sample panel composed of several cattle breeds with different meat quality, validate single nucleotide polymorphisms (SNPs) in a local crossbred population (Angus - Hereford - Limousin) and evaluate their effects on meat quality traits (backfat thickness, intramuscular fat content and fatty acid composition), supporting the association tests with bioinformatic predictive studies. RESULTS: Globally, nine SNPs were detected in the PPARG and CEBPA genes within our mixed panel, including a novel SNP in the latter. Three of these nine, along with seven other SNPs selected from the Single Nucleotide Polymorphism database (SNPdb), including SNPs in the RXRA gene, were validated in the crossbred population (N = 260). After validation, five of these SNPs were evaluated for genotype effects on fatty acid content and composition. Significant effects were observed on backfat thickness and different fatty acid contents (P < 0.05). Some of these SNPs caused slight differences in mRNA structure stability and/or putative binding sites for proteins. CONCLUSIONS: PPARG and CEBPA showed low to moderate variability in our sample panel. Variations in these genes, along with RXRA, may explain part of the genetic variation in fat content and composition. Our results may contribute to knowledge about genetic variation in meat quality traits in cattle and should be evaluated in larger independent populations.

Effectiveness of single-nucleotide polymorphisms to investigate cattle rustling.

Short tandem repeats (STR)s have been the eligible markers for forensic animal genetics, despite single-nucleotide polymorphisms (SNP)s became acceptable. The technology, the type, and amount of markers could limit the investigation in degraded forensic samples. The performance of a 32-SNP panel genotyped through OpenArrays(TM) (real-time PCR based) was evaluated to resolve cattle-specific forensic cases. DNA from different biological sources was used, including samples from an alleged instance of cattle rustling. SNPs and STRs performance and repeatability were compared. SNP call rate was variable among sample type (average = 80.18%), while forensic samples showed the lowest value (70.94%). The repeatability obtained (98.7%) supports the used technology. SNPs had better call rates than STRs in 12 of 20 casework samples, while forensic index values were similar for both panels. In conclusion, the 32-SNPs used are as informative as the standard bovine STR battery and hence are suitable to resolve cattle rustling investigations.

Characterization of the bovine gene LIPE and possible influence on fatty acid composition of meat.

LIPE is an intracellular neutral lipase, which is capable of hydrolyzing a variety of esters and plays a key role in the mobilization of fatty acids from diacylglycerols. The objectives of this study were to characterize the genetic polymorphism of bovine LIPE gene and to evaluate the possible association between three SNPs in the coding regions of this gene with the fatty acid composition of meat in a cattle population. Forty-three unrelated animals from different cattle breeds were re-sequenced and 21 SNPs were detected over approximately 2600 bp, five of these SNPs were novel. Three SNPs were selected, on the basis of evolutionary conservation, to perform validation and association studies in a crossbred cattle population. Our results may suggest a possible association of SNP1 with contents of oleic acid and total monounsaturated fatty acids (p < 0.01), and SNP2 and SNP3 with Heneicosylic acid content (p < 0.01), may be helpful to improve the quality of meat and improve health.

Assessing the association of single nucleotide polymorphisms in thyroglobulin gene with age of puberty in bulls.

Puberty is a stage of sexual development determined by the interaction of many loci and environmental factors. Identification of genes contributing to genetic variation in this character can assist with selection for early pubertal bulls, improving genetic progress in livestock breeding. Thyroid hormones play an important role in sexual development and spermatogenic function. The objective of this study was to evaluate the association between single nucleotide polymorphisms (SNPs) located in thyroglobulin(TG) gene with age of puberty in Angus bulls. Four SNPs were genotyped in 273 animals using SEQUENOM technology and the association between markers and puberty age was analyzed. Results showed a significant association (P < 0.05) between these markers and puberty age estimated at a sperm concentration of 50 million and a progressive motility of 10%. This is the first report of an association of TG polymorphisms with age of puberty in bulls, and results suggest the importance of thyroidal regulation in bovine sexual development and arrival to puberty.

Comparison of the effectiveness of microsatellites and SNP panels for genetic identification, traceability and assessment of parentage in an inbred Angus herd.

During the last decade, microsatellites (short tandem repeats or STRs) have been successfully used for animal genetic identification, traceability and paternity, although in recent year single nucleotide polymorphisms (SNPs) have been increasingly used for this purpose. An efficient SNP identification system requires a marker set with enough power to identify individuals and their parents. Genetic diagnostics generally include the analysis of related animals. In this work, the degree of information provided by SNPs for a consanguineous herd of cattle was compared with that provided by STRs. Thirty-six closely related Angus cattle were genotyped for 18 STRs and 116 SNPs. Cumulative SNPs exclusion power values (Q) for paternity and sample matching probability (MP) yielded values greater than 0.9998 and 4.32E(-42), respectively. Generally 2-3 SNPs per STR were needed to obtain an equivalent Q value. The MP showed that 24 SNPs were equivalent to the ISAG (International Society for Animal Genetics) minimal recommended set of 12 STRs (MP ∼ 10(-11)). These results provide valuable genetic data that support the consensus SNP panel for bovine genetic identification developed by the Parentage Recording Working Group of ICAR (International Committee for Animal Recording).

Association between GNRHR, LHR and IGF1 polymorphisms and timing of puberty in male Angus cattle.

BACKGROUND: In bovines, there are significant differences within and among beef breeds in the time when bulls reach puberty. Although the timing of puberty is likely to be a multigenic trait, previous studies indicate that there may also be single genes that exert major effects on the timing of puberty within the general population. Despite its economic importance, there are not many SNPs or genetic markers associated with the age of puberty in male cattle. In the present work, we selected three candidate genes, GNRHR, LHR and IGF1, and associated their polymorphisms with the age of puberty in Angus male cattle. RESULTS: After weaning, 276 Angus males were measured every month for weight (W), scrotal circumference (SC), sperm concentration (C) and percentage of motility (M). A total of 4 SNPs, two within GNRHR, one in LHR and one in IGF1 were genotyped using the pyrosequencing technique. IGF1-SnaBI SNP was significant associated (P < 0.01) with age at SC 28 cm, but it were not associated with age at M 10% and C 50 million. Genotype CC exhibited an average age at SC 28 cm of 7 and 11 days higher than CT (p = 0.037) and TT (p = 0.012), respectively. This SNP explained 1.5% of the genetic variance of age of puberty at SC28. LHR-I499L, GNRHR-SNP5 and GNRHR-SNP6 were not associated with any of the measurements. However, GNRHR haplotypes showed a suggestive association with age at SC 28 cm. CONCLUSIONS: The findings presented here could support the hypothesis that IGF1 is a regulator of the arrival to puberty in male calves and is involved in the events that precede and initiate puberty in bull calves. Given that most studies in cattle, as well as in other mammals, were done in female, the present results are the first evidence of markers associated with age at puberty in male cattle.