RNA Sequencing of Sperm from Healthy Cattle and Horses Reveals the Presence of a Large Bacterial Population.

RNA molecules within ejaculated sperm can be characterized through whole-transcriptome sequencing, enabling the identification of pivotal transcripts that may influence reproductive success. However, the profiling of sperm transcriptomes through next-generation sequencing has several limitations impairing the identification of functional transcripts. In this study, we explored the nature of the RNA sequences present in the sperm transcriptome of two livestock species, cattle and horses, using RNA sequencing (RNA-seq) technology. Through processing of transcriptomic data derived from bovine and equine sperm cell preparations, low mapping rates to the reference genomes were observed, mainly attributed to the presence of ribosomal RNA and bacteria in sperm samples, which led to a reduced sequencing depth of RNAs of interest. To explore the presence of bacteria, we aligned the unmapped reads to a complete database of bacterial genomes and identified bacteria-associated transcripts which were characterized. This analysis examines the limitations associated with sperm transcriptome profiling by reporting the nature of the RNA sequences among which bacterial RNA was found. These findings can aid researchers in understanding spermatozoal RNA-seq data and pave the way for the identification of molecular markers of sperm performance.

Transcriptome profiling revealed that key rumen epithelium functions change in relation to short-chain fatty acids and rumen epithelium-attached microbiota during the weaning transition.

This study aims to characterize the functional changes of the rumen epithelium associated with ruminal short-chain fatty acid (SCFA) concentration and epithelium-attached microbes during the weaning transition in dairy calves. Ruminal SCFA concentrations were determined, and transcriptome and microbiota profiling in biopsied rumen papillae were obtained from Holstein calves before and after weaning using RNA- and amplicon sequencing. Metabolic pathway analysis showed that pathways related to SCFA metabolism and cell apoptosis were up- and down-regulated postweaning, respectively. Functional analysis showed that genes related to SCFA absorption, metabolism, and protective roles against oxidative stress were positively correlated with ruminal SCFA concentrations. The relative abundance of epithelium-attached Rikenellaceae RC9 gut group and Campylobacter was positively correlated with genes involved in SCFA absorption and metabolism, suggesting that these microbes can cooperatively affect host functions. Future research should examine the contribution of attenuated apoptosis on rumen epithelial functional shifts during the weaning transition.

Unravelling transmission ratio distortion across the bovine genome: identification of candidate regions for reproduction defects.

BACKGROUND: Biological mechanisms affecting gametogenesis, embryo development and postnatal viability have the potential to alter Mendelian inheritance expectations resulting in observable transmission ratio distortion (TRD). Although the discovery of TRD cases have been around for a long time, the current widespread and growing use of DNA technologies in the livestock industry provides a valuable resource of large genomic data with parent-offspring genotyped trios, enabling the implementation of TRD approach. In this research, the objective is to investigate TRD using SNP-by-SNP and sliding windows approaches on 441,802 genotyped Holstein cattle and 132,991 (or 47,910 phased) autosomal SNPs. RESULTS: The TRD was characterized using allelic and genotypic parameterizations. Across the whole genome a total of 604 chromosomal regions showed strong significant TRD. Most (85%) of the regions presented an allelic TRD pattern with an under-representation (reduced viability) of carrier (heterozygous) offspring or with the complete or quasi-complete absence (lethality) for homozygous individuals. On the other hand, the remaining regions with genotypic TRD patterns exhibited the classical recessive inheritance or either an excess or deficiency of heterozygote offspring. Among them, the number of most relevant novel regions with strong allelic and recessive TRD patterns were 10 and 5, respectively. In addition, functional analyses revealed candidate genes regulating key biological processes associated with embryonic development and survival, DNA repair and meiotic processes, among others, providing additional biological evidence of TRD findings. CONCLUSIONS: Our results revealed the importance of implementing different TRD parameterizations to capture all types of distortions and to determine the corresponding inheritance pattern. Novel candidate genomic regions containing lethal alleles and genes with functional and biological consequences on fertility and pre- and post-natal viability were also identified, providing opportunities for improving breeding success in cattle.

Correlation scan: identifying genomic regions that affect genetic correlations applied to fertility traits.

Although the genetic correlations between complex traits have been estimated for more than a century, only recently we have started to map and understand the precise localization of the genomic region(s) that underpin these correlations. Reproductive traits are often genetically correlated. Yet, we don't fully understand the complexities, synergism, or trade-offs between male and female fertility. In this study, we used reproductive traits in two cattle populations (Brahman; BB, Tropical Composite; TC) to develop a novel framework termed correlation scan (CS). This framework was used to identify local regions associated with the genetic correlations between male and female fertility traits. Animals were genotyped with bovine high-density single nucleotide polymorphisms (SNPs) chip assay. The data used consisted of ~1000 individual records measured through frequent ovarian scanning for age at first corpus luteum (AGECL) and a laboratory assay for serum levels of insulin growth hormone (IGF1 measured in bulls, IGF1b, or cows, IGF1c). The methodology developed herein used correlations of 500-SNP effects in a 100-SNPs sliding window in each chromosome to identify local genomic regions that either drive or antagonize the genetic correlations between traits. We used Fisher's Z-statistics through a permutation method to confirm which regions of the genome harboured significant correlations. About 30% of the total genomic regions were identified as driving and antagonizing genetic correlations between male and female fertility traits in the two populations. These regions confirmed the polygenic nature of the traits being studied and pointed to genes of interest. For BB, the most important chromosome in terms of local regions is often located on bovine chromosome (BTA) 14. However, the important regions are spread across few different BTA's in TC. Quantitative trait loci (QTLs) and functional enrichment analysis revealed many significant windows co-localized with known QTLs related to milk production and fertility traits, especially puberty. In general, the enriched reproductive QTLs driving the genetic correlations between male and female fertility are the same for both cattle populations, while the antagonizing regions were population specific. Moreover, most of the antagonizing regions were mapped to chromosome X. These results suggest regions of chromosome X for further investigation into the trade-offs between male and female fertility. We compared the CS with two other recently proposed methods that map local genomic correlations. Some genomic regions were significant across methods. Yet, many significant regions identified with the CS were overlooked by other methods.

Unrevealing functional candidate genes for bovine fertility through RNA sequencing meta-analysis and regulatory elements networks of co-expressed genes and lncRNAs.

The high genetic heterogeneity and environmental effects of subfertility in livestock species make the elucidation of the genetic mechanisms associated with reproductive efficiency a difficult task. Network and co-expression network meta-analyses were applied alongside genetic variant calling and long non-coding RNA (lncRNA) characterization to identify functionally relevant target genes and regulatory subnetworks associated with fertility in dairy cattle. In total, 505 lncRNAs (441 previously annotated in the bovine reference genome ARS-UCD 1.2 and 64 novel lncRNAs) were identified. Seven differentially expressed genes between high-fertile (HF) and sub-fertile (SF) Holstein cows were identified in the network meta-analysis (CA5A, ENSBTAG00000051149, ENSBTAG00000003272, DEFB7, DIO2, TRPV3, and COL4A4). Additionally, seven functional candidate differentially co-expressed (DcoExp) modules with a differential regulatory pattern (|z-score|>2) were identified between HF and SF cows. The functional candidate genes and DcoExp modules identified were associated with fertility relevant processes such as the regulation of embryonic implantation and proliferation, interaction and molecule transfer between the fetus and the cow, and the immune system. These results help to better understand the genetic mechanisms associated with reproductive efficiency in dairy cattle through the identification of potential biomarkers and genetic variants associated with differentially expressed regulatory gene and lncRNAs regulatory element networks.

Small genetic variation affecting mRNA isoforms associated with marbling and meat color in beef cattle.

The aim of this study was to identify mRNA isoforms and small genetic variants that may be affecting marbling and beef color in Nellore cattle. Longissimus thoracis muscle samples from 20 bulls with different phenotypes (out of 80 bulls set) for marbling (moderate (n = 10) and low (n = 10) groups) and beef color (desirable (n = 10) and undesirable (n = 9) group) traits were used to perform transcriptomic analysis using RNA sequencing. Fourteen and 15 mRNA isoforms were detected as differentially expressed (DE) (P-value ≤ 0.001) between divergent groups for marbling and meat color traits, respectively. Some of those DE mRNA isoforms have shown sites of splicing modified by small structural variants as single nucleotide variant (SNV), insertion, and/or deletion. Enrichment analysis identified metabolic pathways, such as O2/CO2 exchange in erythrocytes, tyrosine biosynthesis, and phenylalanine degradation. The results obtained suggest potential key regulatory genes associated with these economically important traits for the beef industry and for the consumer.

Analyses of lambing dates in sheep breeds using von Mises distribution.

Regular changes in the environment and biological responses generate seasonal patterns in the reproduction in small ruminants. Breeding seasonality is a significant constraint impacting efficiency of lamb production. However, seasonality-related traits present a special peculiarity from a statistical point of view being circular data (day of year running 1:365). Recently, circular mixed models have been developed on the basis of the von Mises distribution and were applied to analyse lambing day distribution recorded from five major Canadian sheep breeds (Rideau Arcott, Romanov, Dorset, Suffolk and Polypay). In a simulation study, the linear model was not able to capture the variance components simulated under the circular paradigm; however, the von Mises model evidenced its ability to infer the variance components of simulated circular records. Using real data of sheep, mostly negligible variances were observed for additive genetic effect when using a linear model on circular data values. In contrast, when using the von Mises model, genetic variances were different across breeds, and it raises the possibility to delay the peak of reproduction and to change the seasonality of the ewes. However, a large variance was captured by flock-year effects emphasizing the strong influence of management in lambing seasons for Canadian sheep populations. Finally, the results suggest the potential of using the von Mises model to analyse circular data, and further research is needed for better understand the complexity of this trait and the von Mises models.

Genome-wide association study between copy number variants and hoof health traits in Holstein dairy cattle.

Genome-wide association studies based on SNP have been completed for multiple traits in dairy cattle; however, copy number variants (CNV) could add genomic information that has yet to be harnessed. The objectives of this study were to identify CNV in genotyped Holstein animals and assess their association with hoof health traits using deregressed estimated breeding values as pseudophenotypes. A total of 23,256 CNV comprising 1,645 genomic regions were identified in 5,845 animals. Fourteen genomic regions harboring structural variations, including 9 deletions and 5 duplications, were associated with at least 1 of the studied hoof health traits. This group of traits included digital dermatitis, interdigital dermatitis, heel horn erosion, sole ulcer, white line lesion, sole hemorrhage, and interdigital hyperplasia; no regions were associated with toe ulcer. Twenty candidate genes overlapped with the regions associated with these traits including SCART1, NRXN2, KIF26A, GPHN, and OR7A17. In this study, an effect on infectious hoof lesions could be attributed to the PRAME (Preferentially Expressed Antigen in Melanoma) gene. Almost all genes detected in association with noninfectious hoof lesions could be linked to known metabolic disorders. The knowledge obtained considering information of associated CNV to the traits of interest in this study could improve the accuracy of estimated breeding values. This may further increase the genetic gain for these traits in the Canadian Holstein population, thus reducing the involuntary animal losses due to lameness.

The value of incorporating carcass trait phenotypes in terminal sire selection indexes to improve carcass weight and quality of heavy lambs.

Genetic selection for carcass traits is paramount to maximize the profitability and long-term sustainability of any meat-producing livestock species. The main objectives of this research were to evaluate the efficiency of indicator traits for the genetic improvement of lamb carcass traits and to determine the value of including carcass traits into terminal sire selection indexes for the Canadian sheep industry. The carcass traits included hot carcass weight (HCW), fat depth at the GR site (FATGR) and average carcass conformation score (AVGCONF), and were measured on heavy lambs (slaughter age less than 365 days and HCW greater than 16.3 kg) in commercial abattoirs. Growth traits were found to be moderately efficient indicator traits for the genetic improvement of HCW but selection on ultrasound traits was necessary to substantially improve the carcass quality traits (FATGR and AVGCONF). Economic selection indexes were designed by adding various combinations of carcass traits into the Canadian Sheep Genetic Evaluation System terminal indexes. Records measured on individuals and progeny were assumed to be the sources of information for live animal and carcass traits, respectively. The changes in index accuracy, efficiency and expected correlated response were used to assess the value of their inclusion. HCW was found to have a large economic value, and its inclusion into terminal selection indexes was expected to substantially increase their accuracy (0.08-0.12 points) and efficiency (20%-30%). However, further including FATGR (measured 110 mm from the carcass midline over the 12th rib) and AVGCONF had little impact on the accuracy (≤0.03) and efficiency (1%-7%) of the proposed indexes. Thus, the inclusion of carcass traits into the existing terminal selection indexes could be beneficial for the genetic improvement of HCW, but further research is needed to determine optimal methods of increasing carcass fatness and muscularity.

Use of gene expression profile to identify potentially relevant transcripts to myofibrillar fragmentation index trait.

The myofibrillar fragmentation index (MFI) is an indicative trait for meat tenderness. Longissimus thoracis muscle samples from the 20 most extreme bulls (out of 80 bulls set) for MFI (high (n = 10) and low (n = 10) groups) trait were used to perform transcriptomic analysis, using RNA Sequencing (RNA-Seq). An average of 24.616 genes was expressed in the Nellore muscle transcriptome analysis. A total of 96 genes were differentially expressed (p value ≤ 0.001) between the two groups of divergent bulls for MFI. The HEBP2 and BDH1 genes were overexpressed in animals with high MFI. The MYBPH and MYL6, myosin encoders, were identified. The differentially expressed genes were related to increase mitochondria efficiency, especially in cells under oxidative stress conditions, and these also were related to zinc and calcium binding, membrane transport, and muscle constituent proteins, such as actin and myosin. Most of those genes were involved in metabolic pathways of oxidation-reduction, transport of lactate in the plasma membrane, and muscle contraction. This is the first study applying MFI phenotypes in transcriptomic studies to identify and understand differentially expressed genes for beef tenderness. These results suggest that differences detected in gene expression between high and low MFI animals are related to reactive mechanisms and structural components of oxidative fibers under the condition of cellular stress. Some genes may be selected as positional candidate genes to beef tenderness, MYL6, MYBPH, TRIM63, TRIM55, TRIOBP, and CHRNG genes. The use of MFI phenotypes could enhance results of meat tenderness studies.

Gene expression profiles of bovine genital ridges during sex determination and early differentiation of the gonads†.

Most current knowledge of sex determination in mammals has emerged from mouse and human studies. To investigate the molecular regulation of the sex determination process in cattle, we used an RNA sequencing strategy to analyze the transcriptome landscape of male and female bovine fetal gonads collected in vivo at key developmental stages: before, during, and after SRY gene activation on fetal days D35 (bipotential gonad formation), D39 (peak SRY expression), and D43 (early gonad differentiation). Differentially expressed genes (DEGs) were identified in male vs. female germinal ridges and among group genes showing similar expression profiles during the three periods. There were 143, 96, and 658 DEG between males and female fetuses at D35, D39, and D43, respectively. On D35, genes upregulated in females were enriched in translation, nuclear export, RNA localization, and mRNA splicing events, whereas those upregulated in males were enriched in cell proliferation regulation and male sex determination terms. In time-course experiments, 767 DEGs in males and 545 DEGs in females were identified between D35 vs. D39, and 3157 DEGs in males and 2008 in females were identified between D39 vs. D43. Results highlight unique aspects of sex determination in cattle, such as the expression of several Y chromosome genes (absent in mice and humans) before SRY expression and an abrupt increase in the nuclear expression of SOX10 (instead of SOX9 expression in the Sertoli cell cytoplasm as observed in mice) during male determination and early differentiation.

Intermittent leucine pulses during continuous feeding alters novel components involved in skeletal muscle growth of neonatal pigs.

When neonatal pigs continuously fed formula are supplemented with leucine pulses, muscle protein synthesis and body weight gain are enhanced. To identify the responsible mechanisms, we combined plasma metabolomic analysis with transcriptome expression of the transcriptome and protein catabolic pathways in skeletal muscle. Piglets (n = 23, 7-day-old) were fed continuously a milk replacement formula via orogastric tube for 21 days with an additional parenteral infusion (800 µmol kg-1 h-1) of either leucine (LEU) or alanine (CON) for 1 h every 4 h. Plasma metabolites were measured by liquid chromatography-mass spectrometry. Gene and protein expression analyses of longissimus dorsi muscle were performed by RNA-seq and Western blot, respectively. Compared with CON, LEU pigs had increased plasma levels of leucine-derived metabolites, including 4-methyl-2-oxopentanoate, beta-hydroxyisovalerate, ß-hydroxyisovalerylcarnitine, and 3-methylglutaconate (P ≤ 0.05). Leucine pulses downregulated transcripts enriched in the Kyoto Encyclopedia of Genes and Genomes terms "spliceosome," "GAP junction," "endocytosis," "ECM-receptor interaction," and "DNA replication". Significant correlations were identified between metabolites derived from leucine catabolism and muscle genes involved in protein degradation, transcription and translation, and muscle maintenance and development (P ≤ 0.05). Further, leucine pulses decreased protein expression of autophagic markers and serine/threonine kinase 4, involved in muscle atrophy (P ≤ 0.01). In conclusion, results from our studies support the notion that leucine pulses during continuous enteral feeding enhance muscle mass gain in neonatal pigs by increasing protein synthetic activity and downregulating protein catabolic pathways through concerted responses in the transcriptome and metabolome.

Identification of Ovine Serum miRNAs Following Bacterial Lipopolysaccharide Challenge.

Host-pathogen interactions are complex and influenced by host genetic and epigenetic modifications. Recently, the significance of microRNAs (miRNAs) in pathogenic infection and the regulation of immune response has been highlighted. However, information on miRNAs' role in the course of inflammation is still very limited in small ruminants. The present study was intended to identify changes in the expression of circulatory miRNAs post-lipopolysaccharide (LPS)-challenge. In this study, young ewes (n = 18) were challenged with Escherichia coli LPS (400 ng/kg i.v.) and blood samples were collected for serum miRNA isolation at two-time points; prior to challenge (T0), and 4 h (T4) post-challenge, reflecting the peak cortisol response. A total of 91 miRNAs were profiled, including 84 miRNAs on a commercial ovine miRNA-PCR array, and seven individual miRNAs. Forty five miRNAs were differentially expressed (DE) with 35 being up-regulated (Fold regulation, FR > 2) and 10 being down-regulated (FR < 1, p < 0.05) at T4. Among the up-regulated miRNAs, 14 were significantly (p < 0.05) induced, including oar-miRs: 369-3p, 495-3p, 376a-3p, 543-3p, 668-3p, 329a-3p, 655-3p, 411a-5p, and 154a-3p, which were located on ovine chromosome 18 forming four miRNA clusters within 10 kb. The elevated miRNAs belonged to different functional classes, playing roles in activating the hypothalamic-pituitary-adrenal axis; increasing cell survival and differentiation; and inducing inflammatory responses and targeted PI3K-Akt and MAPK signaling and chemokine signaling pathways. In summary, these results reveal the dynamic nature of ovine serum miRNAs during LPS-induced stress and highlight the potential role of identified miRNA-clusters on chromosome 18 to understand the regulation of the acute-phase response. Some of these identified circulating miRNAs may also serve as stress biomarkers for livestock in the future.

About the existence of common determinants of gene expression in the porcine liver and skeletal muscle.

BACKGROUND: The comparison of expression QTL (eQTL) maps obtained in different tissues is an essential step to understand how gene expression is genetically regulated in a context-dependent manner. In the current work, we have compared the transcriptomic and eQTL profiles of two porcine tissues (skeletal muscle and liver) which typically show highly divergent expression profiles, in 103 Duroc pigs genotyped with the Porcine SNP60 BeadChip (Illumina) and with available microarray-based measurements of hepatic and muscle mRNA levels. Since structural variation could have effects on gene expression, we have also investigated the co-localization of cis-eQTLs with copy number variant regions (CNVR) segregating in this Duroc population. RESULTS: The analysis of differential expresssion revealed the existence of 1204 and 1490 probes that were overexpressed and underexpressed in the gluteus medius muscle when compared to liver, respectively (|fold-change| > 1.5, q-value < 0.05). By performing genome scans in 103 Duroc pigs with available expression and genotypic data, we identified 76 and 28 genome-wide significant cis-eQTLs regulating gene expression in the gluteus medius muscle and liver, respectively. Twelve of these cis-eQTLs were shared by both tissues (i.e. 42.8% of the cis-eQTLs identified in the liver were replicated in the gluteus medius muscle). These results are consistent with previous studies performed in humans, where 50% of eQTLs were shared across tissues. Moreover, we have identified 41 CNVRs in a set of 350 pigs from the same Duroc population, which had been genotyped with the Porcine SNP60 BeadChip by using the PennCNV and GADA softwares, but only a small proportion of these CNVRs co-localized with the cis-eQTL signals. CONCLUSION: Despite the fact that there are considerable differences in the gene expression patterns of the porcine liver and skeletal muscle, we have identified a substantial proportion of common cis-eQTLs regulating gene expression in both tissues. Several of these cis-eQTLs influence the mRNA levels of genes with important roles in meat (CTSF) and carcass quality (TAPT1), lipid metabolism (TMEM97) and obesity (MARC2), thus evidencing the practical importance of dissecting the genetic mechanisms involved in their expression.

Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles.

BACKGROUND: Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles. RESULTS: A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110-114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq. CONCLUSIONS: We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.

Neuropeptidome of the Hypothalamus and Pituitary Gland of Indicine × Taurine Heifers: Evidence of Differential Neuropeptide Processing in the Pituitary Gland before and after Puberty.

Puberty in cattle is regulated by an endocrine axis, which includes a complex milieu of neuropeptides in the hypothalamus and pituitary gland. The neuropeptidome of hypothalamic-pituitary gland tissue of pre- (PRE) and postpubertal (POST) Bos indicus-influenced heifers was characterized, followed by quantitative analysis of 51 fertility-related neuropeptides in these tissues. Comparison of peptide abundances with gene expression levels allowed assessment of post-transcriptional peptide processing. On the basis of classical cleavage, 124 mature neuropeptides from 35 precursor proteins were detected in hypothalamus and pituitary gland tissues of three PRE and three POST Brangus heifers. An additional 19 peptides (cerebellins, PEN peptides) previously reported as neuropeptides that did not follow classical cleavage were also identified. In the pre-pubertal hypothalamus, a greater diversity of neuropeptides (25.8%) was identified relative to post-pubertal heifers, while in the pituitary gland, 38.6% more neuropeptides were detected in the post-pubertal heifers. Neuro-tissues of PRE and POST heifers revealed abundance differences ( p < 0.05) in peptides from protein precursors involved in packaging and processing (e.g., the granin family and ProSAAS) or neuron stimulation (PENK, CART, POMC, cerebellins). On their own, the transcriptome data of the precursor genes could not predict the neuropeptide profile in the exact same tissues in several cases. This provides further evidence of the importance of differential processing of the neuropeptide precursors in the pituitary before and after puberty.

Differential expression of mRNA isoforms in the skeletal muscle of pigs with distinct growth and fatness profiles.

BACKGROUND: The identification of genes differentially expressed in the skeletal muscle of pigs displaying distinct growth and fatness profiles might contribute to identify the genetic factors that influence the phenotypic variation of such traits. So far, the majority of porcine transcriptomic studies have investigated differences in gene expression at a global scale rather than at the mRNA isoform level. In the current work, we have investigated the differential expression of mRNA isoforms in the gluteus medius (GM) muscle of 52 Duroc HIGH (increased backfat thickness, intramuscular fat and saturated and monounsaturated fatty acids contents) and LOW pigs (opposite phenotype, with an increased polyunsaturated fatty acids content). RESULTS: Our analysis revealed that 10.9% of genes expressed in the GM muscle generate alternative mRNA isoforms, with an average of 2.9 transcripts per gene. By using two different pipelines, one based on the CLC Genomics Workbench and another one on the STAR, RSEM and DESeq2 softwares, we have identified 10 mRNA isoforms that both pipelines categorize as differentially expressed in HIGH vs LOW pigs (P-value < 0.01 and ±0.6 log2fold-change). Only five mRNA isoforms, produced by the ITGA5, SEMA4D, LITAF, TIMP1 and ANXA2 genes, remain significant after correction for multiple testing (q-value < 0.05 and ±0.6 log2fold-change), being upregulated in HIGH pigs. CONCLUSIONS: The increased levels of specific ITGA5, LITAF, TIMP1 and ANXA2 mRNA isoforms in HIGH pigs is consistent with reports indicating that the overexpression of these four genes is associated with obesity and metabolic disorders in humans. A broader knowledge about the functional attributes of these mRNA variants would be fundamental to elucidate the consequences of transcript diversity on the determinism of porcine phenotypes of economic interest.

Genetic diversity and signatures of selection in various goat breeds revealed by genome-wide SNP markers.

BACKGROUND: The detection of signatures of selection has the potential to elucidate the identities of genes and mutations associated with phenotypic traits important for livestock species. It is also very relevant to investigate the levels of genetic diversity of a population, as genetic diversity represents the raw material essential for breeding and has practical implications for implementation of genomic selection. A total of 1151 animals from nine goat populations selected for different breeding goals and genotyped with the Illumina Goat 50K single nucleotide polymorphisms (SNP) Beadchip were included in this investigation. RESULTS: The proportion of polymorphic SNPs ranged from 0.902 (Nubian) to 0.995 (Rangeland). The overall mean HO and HE was 0.374 ± 0.021 and 0.369 ± 0.023, respectively. The average pairwise genetic distance (D) ranged from 0.263 (Toggenburg) to 0.323 (Rangeland). The overall average for the inbreeding measures FEH, FVR, FLEUT, FROH and FPED was 0.129, -0.012, -0.010, 0.038 and 0.030, respectively. Several regions located on 19 chromosomes were potentially under selection in at least one of the goat breeds. The genomic population tree constructed using all SNPs differentiated breeds based on selection purpose, while genomic population tree built using only SNPs in the most significant region showed a great differentiation between LaMancha and the other breeds. We hypothesized that this region is related to ear morphogenesis. Furthermore, we identified genes potentially related to reproduction traits, adult body mass, efficiency of food conversion, abdominal fat deposition, conformation traits, liver fat metabolism, milk fatty acids, somatic cells score, milk protein, thermo-tolerance and ear morphogenesis. CONCLUSIONS: In general, moderate to high levels of genetic variability were observed for all the breeds and a characterization of runs of homozygosity gave insights into the breeds' development history. The information reported here will be useful for the implementation of genomic selection and other genomic studies in goats. We also identified various genome regions under positive selection using smoothed FST and hapFLK statistics and suggested genes, which are potentially under selection. These results can now provide a foundation to formulate biological hypotheses related to selection processes in goats.

Parallel responses of human epidermal keratinocytes to inorganic SbIII and AsIII.

SbIII and AsIII are known to exhibit similar chemical properties, but the degree of similarity in their effects on biological systems merits further exploration. Present work compares the responses of human epidermal keratinocytes, a known target cell type for arsenite-induced carcinogenicity, to these metalloids after treatment for a week at environmentally relevant concentrations. Previous work with these cells has shown that arsenite and antimonite have parallel effects in suppressing differentiation, altering levels of several critical enzymes and maintaining colony forming ability. More globally, protein profiling now reveals parallels in SbIII and AsIII effects. The more sensitive technique of transcriptional profiling also shows considerable parallels. Thus, gene expression changes were almost entirely in the same directions for the two treatments, although the degree of change was sometimes significantly different. Inspection of the changes revealed that RYR1 and LRIG1 were among the genes strongly suppressed, consistent with reduced calcium-dependent differentiation and maintenance of EGF-dependent proliferative potential. Moreover, levels of miRNAs in the cells were altered in parallel, with nearly 90% of the 198 most highly expressed ones being suppressed. Among these was miR-203, which is known to decrease proliferative potential. Finally, both SbIII and AsIII were seen to attenuate bone morphogenetic protein 6 induction of dual specificity phosphatases 2 and 14, consistent with maintaining epidermal growth factor receptor signaling. These findings raise the question whether SbIII, like AsIII, could act as a human skin carcinogen.

Single nucleotide polymorphism (SNP) variation of wolves (Canis lupus) in Southeast Alaska and comparison with wolves, dogs, and coyotes in North America.

There is considerable interest in the genetics of wolves (Canis lupus) because of their close relationship to domestic dogs (C. familiaris) and the need for informed conservation and management. This includes wolf populations in Southeast Alaska for which we determined genotypes of 305 wolves at 173662 single nucleotide polymorphism (SNP) loci. After removal of invariant and linked SNP, 123801 SNP were used to quantify genetic differentiation of wolves in Southeast Alaska and wolves, coyotes (C. latrans), and dogs from other areas in North America. There is differentiation of SNP allele frequencies between the species (wolves, coyotes, and dogs), although differentiation is relatively low between some wolf and coyote populations. There are varying levels of differentiation among populations of wolves, including low differentiation of wolves in interior Alaska, British Columbia, and the northern US Rocky Mountains. There is considerable differentiation of SNP allele frequencies of wolves in Southeast Alaska from wolves in other areas. However, wolves in Southeast Alaska are not a genetically homogeneous group and there are comparable levels of genetic differentiation among areas within Southeast Alaska and between Southeast Alaska and other geographic areas. SNP variation and other genetic data are discussed regarding taxonomy and management.