Identification of quantitative trait transcripts for growth traits in the large scales of liver and muscle samples

RNA Sequencing of the Pituitary Gland and Association Analyses Reveal PRKG2 as a Candidate Gene for Growth and Carcass Traits in Chinese Ningdu Yellow Chickens

Growth and carcass traits are of great economic importance to the chicken industry. The candidate genes and mutations associated with growth and carcass traits can be utilized to improve chicken growth. Therefore, the identification of these genes and mutations is greatly importance. In this study, a total of 17 traits related to growth and carcass were measured in 399 Chinese Ningdu yellow chickens. RNA sequencing (RNA-seq) was performed to detect candidate genes using 12 pituitary gland samples (six per group), which exhibited extreme growth and carcass phenotypes: either a high live weight and carcass weight (H group) or a low live weight and carcass weight (L group). A differential expression analysis, utilizing RNA-seq, between the H and L groups identified 428 differentially expressed genes (DEGs), including 110 up-regulated genes and 318 down-regulated genes. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of the identified genes showed a significant enrichment of 158 GO terms and two KEGG pathways, including response to stimulus and neuroactive ligand-receptor interaction, respectively. Furthermore, RNA-seq data, qRT–PCR, and quantitative trait transcript (QTT) analysis results suggest that the PRKG2 gene is an important candidate gene for growth and carcass traits of Chinese Ningdu yellow chickens. More specifically, association analyses of a single nucleotide polymorphism (SNP) in PRKG2 and growth and carcass traits showed that the SNP rs16400745 was significantly associated with 12 growth and carcass traits (P < 0.05), such as carcass weight (P = 9.68E-06), eviscerated weight (P = 3.04E-05), and semi-eviscerated weight (P = 2.14E-04). Collectively, these results provide novel insights into the genetic basis of growth in Chinese Ningdu yellow chickens and the SNP rs16400745 reported here could be incorporated into the selection programs involving this breed.

Unravelling the effect of environment on the genome of Sarda breed ewes using Runs of Homozygosity

Natural adaptation and artificial selection have shaped the genome of modern livestock breeds. Among SNP-based metrics that are used to detect signatures of selection at genome-wide level, runs of homozygosity (ROH) are getting increasing popularity. In this paper, ROH distribution and features of a sample of 823 Sarda breed ewes farmed at different levels of altitude are analysed to investigate the effect of the environment on the patterns of homozygosity. A total of 46,829 (33,087 unique) ROH were detected. OAR2 exhibited the largest average number of ROH per animal. The most frequent ROH (OAR27, 38.9-44.2 Mb) was shared by 327. ROH length was statistically affected (p < 0.001) by both the altitude and temperature of the place where the flock was located. The highest probability of a SNP falling in a ROH was observed for hill ewes, whereas the smallest one for mountain. A total of 457 SNP exceeded the 99th percentile of the ROH count per SNP distribution and were considered significant. These markers mapped in eight chromosomes and they clustered into 17 ROH islands, where 80 candidate genes were mapped. Results of this study highlighted differences in the ROH distribution and features among sheep farmed in flocks located at different levels of altitude, confirming the role of environmental adaptability in shaping the genome of this breed.

The Contribution of Copy Number Variants and Single Nucleotide Polymorphisms to the Additive Genetic Variance of Carcass Traits in Cattle

The relative contributions of both copy number variants (CNVs) and single nucleotide polymorphisms (SNPs) to the additive genetic variance of carcass traits in cattle is not well understood. A detailed understanding of the relative importance of CNVs in cattle may have implications for study design of both genomic predictions and genome-wide association studies. The first objective of the present study was to quantify the relative contributions of CNV data and SNP genotype data to the additive genetic variance of carcass weight, fat, and conformation for 945 Charolais, 923 Holstein-Friesian, and 974 Limousin sires. The second objective was to jointly consider SNP and CNV data in a least absolute selection and shrinkage operator (LASSO) regression model to identify genomic regions associated with carcass weight, fat, and conformation within each of the three breeds separately. A genomic relationship matrix (GRM) based on just CNV data did not capture any variance in the three carcass traits when jointly evaluated with a SNP-derived GRM. In the LASSO regression analysis, a total of 987 SNPs and 18 CNVs were associated with at least one of the three carcass traits in at least one of the three breeds. The quantitative trait loci (QTLs) corresponding to the associated SNPs and CNVs overlapped with several candidate genes including previously reported candidate genes such as MSTN and RSAD2, and several potential novel candidate genes such as ACTN2 and THOC1. The results of the LASSO regression analysis demonstrated that CNVs can be used to detect associations with carcass traits which were not detected using the set of SNPs available in the present study. Therefore, the CNVs and SNPs available in the present study were not redundant forms of genomic data.

An Integrative Analysis of Transcriptome and GWAS Data to Identify Potential Candidate Genes Influencing Meat Quality Traits in Pigs

Understanding the genetic factors behind meat quality traits is of great significance to animal breeding and production. We previously conducted a genome-wide association study (GWAS) for meat quality traits in a White Duroc × Erhualian F2 pig population using Illumina porcine 60K SNP data. Here, we further investigate the functional candidate genes and their network modules associated with meat quality traits by integrating transcriptomics and GWAS information. Quantitative trait transcript (QTT) analysis, gene expression QTL (eQTL) mapping, and weighted gene co-expression network analysis (WGCNA) were performed using the digital gene expression (DGE) data from 493 F2 pig's muscle and liver samples. Among the quantified 20,108 liver and 23,728 muscle transcripts, 535 liver and 1,014 muscle QTTs corresponding to 416 and 721 genes, respectively, were found to be significantly (p < 5 × 10^-4) correlated with 22 meat quality traits measured on longissiums dorsi muscle (LM) or semimembranosus muscle (SM). Transcripts associated with muscle glycolytic potential (GP) and pH values were enriched for genes involved in metabolic process. There were 42 QTTs (for 32 genes) shared by liver and muscle tissues, of which 10 QTTs represent GP- and/or pH-related genes, such as JUNB, ATF3, and PPP1R3B. Furthermore, a genome-wide eQTL mapping revealed a total of 3,054 eQTLs for all annotated transcripts in muscle (p < 2.08 × 10^-5), including 1,283 cis-eQTLs and 1771 trans-eQTLs. In addition, WGCNA identified five modules relevant to glycogen metabolism pathway and highlighted the connections between variations in meat quality traits and genes involved in energy process. Integrative analysis of GWAS loci, eQTL, and QTT demonstrated GALNT15/GALNTL2 and HTATIP2 as strong candidate genes for drip loss and pH drop from postmortem 45 min to 24 h, respectively. Our findings provide valuable insights into the genetic basis of meat quality traits and greatly expand the number of candidate genes that may be valuable for future functional analysis and genetic improvement of meat quality.

Whole-exome sequencing identifies novel mutations in ABC transporter genes associated with intrahepatic cholestasis of pregnancy disease: a case-control study

Background

Intrahepatic cholestasis of pregnancy (ICP) can cause premature delivery and stillbirth. Previous studies have reported that mutations in ABC transporter genes strongly influence the transport of bile salts. However, to date, their effects are still largely elusive.

Methods

A whole-exome sequencing (WES) approach was used to detect novel variants. Rare novel exonic variants (minor allele frequencies: MAF < 1%) were analyzed. Three web-available tools, namely, SIFT, Mutation Taster and FATHMM, were used to predict protein damage. Protein structure modeling and comparisons between reference and modified protein structures were performed by SWISS-MODEL and Chimera 1.14rc, respectively.

Results

We detected a total of 2953 mutations in 44 ABC family transporter genes. When the MAF of loci was controlled in all databases at less than 0.01, 320 mutations were reserved for further analysis. Among these mutations, 42 were novel. We classified these loci into four groups (the damaging, probably damaging, possibly damaging, and neutral groups) according to the prediction results, of which 7 novel possible pathogenic mutations were identified that were located in known functional genes, including ABCB4 (Trp708Ter, Gly527Glu and Lys386Glu), ABCB11 (Gln1194Ter, Gln605Pro and Leu589Met) and ABCC2 (Ser1342Tyr), in the damaging group. New mutations in the first two genes were reported in our recent article. In addition, compared to the wild-type protein structure, the ABCC2 Ser1342Tyr-modified protein structure showed a slight change in the chemical bond lengths of ATP ligand-binding amino acid side chains. In placental tissue, the expression level of the ABCC2 gene in patients with ICP was significantly higher (P < 0.05) than that in healthy pregnant women. In particular, the patients with two mutations in ABC family genes had higher average values of total bile acids (TBA), aspartate transaminase (AST), direct bilirubin (DBIL), total cholesterol (CHOL), triglycerides (TG) and high-density lipoprotein (HDL) than the patients who had one mutation, no mutation in ABC genes and local controls.

Conclusions

Our present study provide new insight into the genetic architecture of ICP and will benefit the final identification of the underlying mutations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12884-021-03595-x.

Signatures of positive selection underlying beef production traits in Korean cattle breeds

The difference in the breeding programs and population history may have diversely shaped the genomes of Korean native cattle breeds. In the absence of phenotypic data, comparisons of breeds that have been subjected to different selective pressures can aid to identify genomic regions and genes controlling qualitative and complex traits. In this study to decipher genetic variation and identify evidence of divergent selection, 3 Korean cattle breeds were genotyped using the recently developed high-density GeneSeek Genomic Profiler F250 (GGP-F250) array. The three Korean cattle breeds clustered according to their coat color phenotypes and breeding programs. The Heugu breed reliably showed smaller effective population size at all generations considered. Across the autosomal chromosomes, 113 and 83 annotated genes were identified from Hanwoo-Chikso and Hanwoo-Heugu comparisons, respectively of which 16 genes were shared between the two pairwise comparisons. The most important signals of selection were detected on bovine chromosomes 14 (24.39–25.13 Mb) and 18 (13.34–15.07 Mb), containing genes related to body size, and coat color (XKR4, LYN, PLAG1, SDR16C5, TMEM68, CDH15, MC1R, and GALNS). Some of the candidate genes are also associated with meat quality traits (ACSF3, EIF2B1, BANP, APCDD1, and GALM) and harbor quantitative trait locus (QTL) for beef production traits. Further functional analysis revealed that the candidate genes (DBI, ACSF3, HINT2, GBA2, AGPAT5, SCAP, ELP6, APOB, and RBL1) were involved in gene ontology (GO) terms relevant to meat quality including fatty acid oxidation, biosynthesis, and lipid storage. Candidate genes previously known to affect beef production and quality traits could be used in the beef cattle selection strategies.

Generation of Retinaldehyde for Retinoic Acid Biosynthesis

The concentration of all-trans-retinoic acid, the bioactive derivative of vitamin A, is critically important for the optimal performance of numerous physiological processes. Either too little or too much of retinoic acid in developing or adult tissues is equally harmful. All-trans-retinoic acid is produced by the irreversible oxidation of all-trans-retinaldehyde. Thus, the concentration of retinaldehyde as the immediate precursor of retinoic acid has to be tightly controlled. However, the enzymes that produce all-trans-retinaldehyde for retinoic acid biosynthesis and the mechanisms responsible for the control of retinaldehyde levels have not yet been fully defined. The goal of this review is to summarize the current state of knowledge regarding the identities of physiologically relevant retinol dehydrogenases, their enzymatic properties, and tissue distribution, and to discuss potential mechanisms for the regulation of the flux from retinol to retinaldehyde.

Mice lacking the epidermal retinol dehydrogenases SDR16C5 and SDR16C6 display accelerated hair growth and enlarged meibomian glands

Retinol dehydrogenases catalyze the rate-limiting step in the biosynthesis of retinoic acid, a bioactive lipid molecule that regulates the expression of hundreds of genes by binding to nuclear transcription factors, the retinoic acid receptors. Several enzymes exhibit retinol dehydrogenase activities in vitro; however, their physiological relevance for retinoic acid biosynthesis in vivo remains unclear. Here, we present evidence that two murine epidermal retinol dehydrogenases, short-chain dehydrogenase/reductase family 16C member 5 (SDR16C5) and SDR16C6, contribute to retinoic acid biosynthesis in living cells and are also essential for the oxidation of retinol to retinaldehyde in vivo Mice with targeted knockout of the more catalytically active SDR16C6 enzyme have no obvious phenotype, possibly due to functional redundancy, because Sdr16c5 and Sdr16c6 exhibit an overlapping expression pattern during later developmental stages and in adulthood. Mice that lack both enzymes are viable and fertile but display accelerated hair growth after shaving and also enlarged meibomian glands, consistent with a nearly 80% reduction in the retinol dehydrogenase activities of skin membrane fractions from the Sdr16c5/Sdr16c6 double-knockout mice. The up-regulation of hair-follicle stem cell genes is consistent with reduced retinoic acid signaling in the skin of the double-knockout mice. These results indicate that the retinol dehydrogenase activities of murine SDR16C5 and SDR16C6 enzymes are not critical for survival but are responsible for most of the retinol dehydrogenase activity in skin, essential for the regulation of the hair-follicle cycle, and required for the maintenance of both sebaceous and meibomian glands.

Characterization of human short chain dehydrogenase/reductase SDR16C family members related to retinol dehydrogenase 10

All-trans-retinoic acid (RA) is a bioactive derivative of vitamin A that serves as an activating ligand for nuclear transcription factors, retinoic acid receptors. RA biosynthesis is initiated by the enzymes that oxidize retinol to retinaldehyde. It is well established that retinol dehydrogenase 10 (RDH10, SDR16C4), which belongs to the 16C family of the short chain dehydrogenase/reductase (SDR) superfamily of proteins, is the major enzyme responsible for the oxidation of retinol to retinaldehyde for RA biosynthesis during embryogenesis. However, several lines of evidence point towards the existence of additional retinol dehydrogenases that contribute to RA biosynthesis in vivo. In close proximity to RDH10 gene on human chromosome 8 are located two genes that are phylogenetically related to RDH10. The predicted protein products of these genes, retinol dehydrogenase epidermal 2 (RDHE2, SDR16C5) and retinol dehydrogenase epidermal 2-similar (RDHE2S, SDR16C6), share 59% and 56% sequence similarity with RDH10, respectively. Previously, we showed that the single ortholog of the human RDHE2 and RDHE2S in frogs, Xenopus laevis rdhe2, oxidizes retinol to retinaldehyde and is essential for frog embryonic development. In this study, we explored the potential of each of the two human proteins to contribute to RA biosynthesis. The results of this study demonstrate that human RDHE2 exhibits a relatively low but reproducible activity when expressed in either HepG2 or HEK293 cells. Expression of the native RDHE2 is downregulated in the presence of elevated levels of RA. On the other hand, the protein encoded by the human RDHE2S gene is unstable when expressed in HEK293 cells. RDHE2S protein produced in Sf9 cells is stable but has no detectable catalytic activity towards retinol. We conclude that the human RDHE2S does not contribute to RA biosynthesis, whereas the low-activity RA-sensitive human RDHE2 may have a role in adjusting the cellular levels of RA in accord with specific physiological conditions.