RNA-Seq reveals differentially expressed genes affecting polyunsaturated fatty acids percentage in the Huangshan Black chicken population

Genetic architecture and key regulatory genes of fatty acid composition in Gushi chicken breast muscle determined by GWAS and WGCNA

BACKGROUND

Fatty acids composition in poultry muscle is directly related to its tenderness, flavour, and juiciness, whereas its genetic mechanisms have not been elucidated. In this study, the genetic structure and key regulatory genes of the breast muscle fatty acid composition of local Chinese chicken, Gushi-Anka F2 resource population by integrating genome-wide association study (GWAS) and weighted gene co-expression network analysis (WGCNA) strategies. GWAS was performed based on 323,306 single nucleotide polymorphisms (SNPs) obtained by genotyping by sequencing (GBS) method and 721 chickens from the Gushi-Anka F2 resource population with highly variable fatty acid composition traits in the breast muscle. And then, according to the transcriptome data of the candidate genes that were obtained and phenotypic data of fatty acid composition traits in breast muscle of Gushi chickens at 14, 22, and 30 weeks of age, we conducted a WGCNA.

RESULTS

A total of 128 suggestive significantly associated SNPs for 11 fatty acid composition traits were identified and mapped on chromosomes (Chr) 2, 3, 4, 5, 13, 17, 21, and 27. Of these, the two most significant SNPs were Chr13:5,100,140 (P = 4.56423e-10) and Chr13:5,100,173 (P = 4.56423e-10), which explained 5.6% of the phenotypic variation in polyunsaturated fatty acids (PUFA). In addition, six fatty acid composition traits, including C20:1, C22:6, saturated fatty acid (SFA), unsaturated fatty acids (UFA), PUFA, and average chain length (ACL), were located in the same QTL intervals on Chr13. We obtained 505 genes by scanning the linkage disequilibrium (LD) regions of all significant SNPs and performed a WGCNA based on the transcriptome data of the above 505 genes. Combining two strategies, 9 hub genes (ENO1, ADH1, ASAH1, ADH1C, PIK3CD, WISP1, AKT1, PANK3, and C1QTNF2) were finally identified, which could be the potential candidate genes regulating fatty acid composition traits in chicken breast muscle.

CONCLUSION

The results of this study deepen our understanding of the genetic mechanisms underlying the regulation of fatty acid composition traits, which is helpful in the design of breeding strategies for the subsequent improvement of fatty acid composition in poultry muscle.

Integrated Transcriptomics Profiling in Chahua and Digao Chickens' Breast for Assessment Molecular Mechanism of Meat Quality Traits

Meat quality traits are an important economic trait and remain a major argument, from the producer to the consumer. However, there are a few candidate genes and pathways of chicken meat quality traits that were reported for chicken molecular breeding. The purpose of the present study is to identify the candidate genes and pathways associated with meat quality underlying variations in meat quality. Hence, transcriptome profiles of breast tissue in commercial Digao (DG, 5 male) and Chahua (CH, 5 male) native chicken breeds were analyzed at the age of 100 days. The results found 3525 differentially expressed genes (DEGs) in CH compared to DG with adjusted p-values of ≤0.05 and log2FC ≥ 0.1 FDR ≤ 0.05. Functional analysis of GO showed that the DEGs are mainly involved in the two types of processes of meat quality, such as positive regulation of the metabolic process, extracellular structure organization, collagen trimer, cellular amino acid metabolic process, cellular amino acid catabolic process, and heme binding. Functional analysis of KEGG showed that the DEGs are mainly involved in the two types of processes of meat quality, such as oxidative phosphorylation, carbon metabolism, valine, leucine, and isoleucine degradation, and fatty acid degradation. Many of the DEGs are well known to be related to meat quality, such as COL28A1, COL1A2, MB, HBAD, HBA1, ACACA, ACADL, ACSL1, ATP8A1, CAV1, FADS2, FASN, DCN, CHCHD10, AGXT2, ALDH3A2, and MORN4. Therefore, the current study detected multiple pathways and genes that could be involved in the control of the meat quality traits of chickens. These findings should be used as an essential resource to improve the accuracy of selection for meat traits in chickens using marker-assisted selection based on differentially expressed genes.

Investigation of clinical characteristics and genome associations in the 'UK Lipoedema' cohort

Lipoedema is a chronic adipose tissue disorder mainly affecting women, causing excess subcutaneous fat deposition on the lower limbs with pain and tenderness. There is often a family history of lipoedema, suggesting a genetic origin, but the contribution of genetics is currently unclear. A tightly phenotyped cohort of 200 lipoedema patients was recruited from two UK specialist clinics. Objective clinical characteristics and measures of quality of life data were obtained. In an attempt to understand the genetic architecture of the disease better, genome-wide single nucleotide polymorphism (SNP) genotype data were obtained, and a genome wide association study (GWAS) was performed on 130 of the recruits. The analysis revealed genetic loci suggestively associated with the lipoedema phenotype, with further support provided by an independent cohort taken from the 100,000 Genomes Project. The top SNP rs1409440 (ORmeta ≈ 2.01, Pmeta ≈ 4 x 10-6) is located upstream of LHFPL6, which is thought to be involved with lipoma formation. Exactly how this relates to lipoedema is not yet understood. This first GWAS of a UK lipoedema cohort has identified genetic regions of suggestive association with the disease. Further replication of these findings in different populations is warranted.

Transcriptome-metabolome analysis reveals how sires affect meat quality in hybrid sheep populations

Crossbreeding improves and enhances meat quality and is widely used in sheep production; however, the molecular mechanisms underlying the meat quality of various crossbred sheep remain unknown. In this study, male Southdown, Suffolk and Hu sheep were crossbred with female Hu sheep, and the transcriptomes and metabolomes of the longissimus dorsi muscle of the F1 generation were sequenced to explore how different sire breeds affect meat quality. The results showed that 631 differentially expressed genes and 119 significantly altered metabolites contributed to muscle development characteristics and meat quality-related diversity (P < 0.05). These genes and metabolites were significantly enriched in lipid metabolism pathways, including arachidonic acid metabolism and PPAR signaling. Several candidate genes were associated with muscle growth, such as MYLK3, MYL10, FIGN, MYH8, MYOM3, LMCD1, and FLRT1. Among these, MYH8 and MYL10 participated in regulating muscle growth and development and were correlated with meat quality-related fatty acid levels (|r| > 0.5 and p < 0.05). We selected mRNA from four of these genes to verify the accuracy of the sequencing data via qRT-PCR. Our findings provide further insight into the key genes and metabolites involved in muscle growth and meat quality in hybrid sheep populations.

Whole transcriptome RNA Sequencing Reveals the Global Molecular Responses and circRNA/lncRNA-miRNA-mRNA ceRNA Regulatory Network in Chicken Fat Deposition

Fat deposition is a vital factor affecting the economics of poultry production. Numerous studies on fat deposition have been done. However, the molecular regulatory mechanism is still unclear. In the present study, the whole-transcriptome RNA sequencing in abdominal fat, back skin, and liver both high- and low-abdominal fat groups was used to uncover the competitive endogenous RNA (ceRNA) regulation network related to chicken fat deposition. The results showed that differentially expressed (DE) genes in abdominal fat, back skin, liver were 1207(784 mRNAs, 330 lncRNAs, 41 circRNAs, 52 miRNAs), 860 (607 mRNAs, 166 lncRNAs, 26 circRNAs, 61 miRNAs), and 923 (501 mRNAs, 262 lncRNAs, 15 circRNAs, 145 miRNAs), respectively. The ceRNA regulatory network analysis indicated that the fatty acid metabolic process, monocarboxylic acid metabolic process, carboxylic acid metabolic process, glycerolipid metabolism, fatty acid metabolism, and peroxisome proliferator-activated receptor (PPAR) signaling pathway took part in chicken fat deposition. Meanwhile, we scan the important genes, FADS2, HSD17B12, ELOVL5, AKR1E2, DGKQ, GPAM, PLIN2, which were regulated by gga-miR-460b-5p, gga-miR-199-5p, gga-miR-7470-3p, gga-miR-6595-5p, gga-miR-101-2-5p. While these miRNAs were competitive combined by lncRNAs including MSTRG.18043, MSTRG.7738, MSTRG.21310, MSTRG.19577, and circRNAs including novel_circ_PTPN2, novel_circ_CTNNA1, novel_circ_PTPRD. This finding provides new insights into the regulatory mechanism of mRNA, miRNA, lncRNA, and circRNA in chicken fat deposition.

Seasonal and sex-dependent gene expression in emu (Dromaius novaehollandiae) fat tissues

Emu (Dromaius novaehollandiae) farming has been gaining wide interest for fat production. Oil rendered from this large flightless bird’s fat is valued for its anti-inflammatory and antioxidant properties for uses in therapeutics and cosmetics. We analyzed the seasonal and sex-dependent differentially expressed (DE) genes involved in fat metabolism in emus. Samples were taken from back and abdominal fat tissues of a single set of four male and four female emus in April, June, and November for RNA-sequencing. We found 100 DE genes (47 seasonally in males; 34 seasonally in females; 19 between sexes). Seasonally DE genes with significant difference between the sexes in gene ontology terms suggested integrin beta chain-2 (ITGB2) influences fat changes, in concordance with earlier studies. Six seasonally DE genes functioned in more than two enriched pathways (two female: angiopoietin-like 4 (ANGPTL4) and lipoprotein lipase (LPL); four male: lumican (LUM), osteoglycin (OGN), aldolase B (ALDOB), and solute carrier family 37 member 2 (SLC37A2)). Two sexually DE genes, follicle stimulating hormone receptor (FSHR) and perilipin 2 (PLIN2), had functional investigations supporting their influence on fat gain and loss. The results suggested these nine genes influence fat metabolism and deposition in emus.

Seasonal and sexual variation in mRNA expression of selected adipokine genes affecting fat deposition and metabolism of the emu (Dromaius novaehollandiae)

Emus are farmed for fat production. Oil rendered from their back and abdominal fat pads has good anti-oxidant and anti-inflammatory properties and has ingredients that promote cell growth. Our objective is to examine the mRNA expression of 7 emu adipokine genes (eFABP4, eSCD1, eAdipoQ, eAdipoR1, eAdipoR2, eLEP and eLepR) to identify gene markers that may help improve emu fat production. Back and abdominal fat tissues from 11 adult emus were biopsied at four time points (April, June, August and November). Total RNA was isolated and cDNA was synthesized. Gene specific primers were designed for partial cloning fragments to amplify the open reading frame of the 7 genes. eLEP was not expressed in emu fat tissue. Nucleotides and amino acids sequences of the 6 expressed gene were compared with homologs from other species and phylogenetic relationships established. Seasonal mRNA expression of each gene was assessed by quantitative RT-PCR and differential expression analysed by the 2^−ΔΔC_T method. The 6 expressed genes showed seasonal variation in expression and showed association of expression level with back fat adiposity. More whole-genome scanning studies are needed to develop novel molecular markers that can be applied to improve fat production in emus.

Molecular Regulation of Lipogenesis, Adipogenesis and Fat Deposition in Chicken

In the poultry industry, excessive fat deposition is considered an undesirable factor, affecting feed efficiency, meat production cost, meat quality, and consumer’s health. Efforts to reduce fat deposition in economically important animals, such as chicken, can be made through different strategies; including genetic selection, feeding strategies, housing, and environmental strategies, as well as hormone supplementation. Recent investigations at the molecular level have revealed the significant role of the transcriptional and post-transcriptional regulatory networks and their interaction on modulating fat metabolism in chickens. At the transcriptional level, different transcription factors are known to regulate the expression of lipogenic and adipogenic genes through various signaling pathways, affecting chicken fat metabolism. Alternatively, at the post-transcriptional level, the regulatory mechanism of microRNAs (miRNAs) on lipid metabolism and deposition has added a promising dimension to understand the structural and functional regulatory mechanism of lipid metabolism in chicken. Therefore, this review focuses on the progress made in unraveling the molecular function of genes, transcription factors, and more notably significant miRNAs responsible for regulating adipogenesis, lipogenesis, and fat deposition in chicken. Moreover, a better understanding of the molecular regulation of lipid metabolism will give researchers novel insights to use functional molecular markers, such as miRNAs, for selection against excessive fat deposition to improve chicken production efficiency and meat quality.

Complete mitochondrial genome and phylogenetic analysis of Huangshan Black chicken (Gallus gallus)

In this study, we sequenced the complete mitochondrial genome (mitogenome) of the Huangshan Black chicken (HBC). Results showed that the complete HBC mitogenome was 16,785 bp in size, comprising 22 transfer RNA genes, 2 ribosomal RNA genes, 13 protein-coding genes, and 1 non-coding control region (D-loop). The overall nucleotide composition was 32.5% for C, 30.3% for A, 23.7% for T, and 13.5% for G. Phylogenetic analysis showed that the HBC mitogenome was clustered with Xianju chicken, which belonged to the haplogroup D2. Our results therefore demonstrate that the origin of HBC corresponds to haplogroup D2 distribution and might have at least one maternal lineage originated from Southeast Asia.

Prospecting genomic regions associated with milk production traits in Egyptian buffalo

The objectives of the current study were to detect putative genomic loci and to identify candidate genes associated with milk production traits in Egyptian buffalo. A total number of 161 479 daily milk yield (DMY) records and 60 318 monthly measures for fat and protein percentages (FP and PP, respectively), along with fat and protein yields (FY and PY, respectively) from 1670 animals were used. Genotyping was performed using Axiom® Buffalo Genotyping 90 K array. Genome-wide association study (GWAS) for each trait was performed using PLINK. After Bonferroni correction, 47 SNPs were associated with one or more milk production traits. These SNPs were distributed over 36 quantitative trait loci (QTL) and located on 20 buffalo chromosomes (BBU). For the 47 SNPs, one was overlapped for three traits (DMY, FY, and PY), six were associated with two traits (one for PP and PY and five for FY and PY) while the rest were associated with only one trait. Out of 36 identified QTL, eleven were overlapped with previously reported loci in buffalo and/or cattle populations. Some of these SNPs are placed within or close to potential candidate genes, for example: TPD52, ZBTB10, RALYL and SNX16 on BBU15, ADGRD1 on BBU17, ESRRG on BBU5 and GRIP1 on BBU4. This is the first reported study between genome-wide markers and milk components in Egyptian buffalo. Our findings provide useful information to explore the genetic mechanisms and relevant genes contributing to the variation in milk production traits. Further confirmation studies with larger population size are necessary to validate the findings and detect the causal genetic variants.

Isolation and Identification of Bovine Preadipocytes and Screening of MicroRNAs Associated with Adipogenesis

Simple Summary

Promoting fat deposition in beef cattle has been a focus of modern animal breeding research. However, previous researchers have not examined the mechanism of adipogenesis in much detail. MicroRNAs (miRNAs) are small noncoding RNAs that play a pivotal role in adipogenesis. In this study, to explore the molecular regulatory mechanism of adipocyte differentiation and formation, bovine preadipocytes were isolated and induced into adipocytes, and then the expression patterns of miRNAs between preadipocytes and adipocytes were detected through RNA sequencing. Deep sequence analysis has identified 78, 71, and 48 novel miRNAs and 497, 491, and 524 known miRNAs in the preadipocytes, and 44, 54, and 47 novel miRNAs and 519, 522, and 504 known miRNAs in the adipocytes. Among the annotated miRNAs, 131 bovine miRNAs were upregulated in adipocytes, and 119 bovine miRNAs were downregulated in adipocytes, such as bta-miR-3604, bta-miR-23b-3p, bta-miR-26a, and bta-miR-129-3p. Bovine target gene prediction results of these miRNAs show that numerous genes are associated with lipid metabolism. These results can provide both technical support and a research basis for promoting bovine adipocyte fat deposition.

Abstract

The elucidation of the mechanisms of preadipocyte differentiation and fat accumulation in adipocytes is a major work in beef cattle breeding. As important post-transcriptional regulators, microRNAs (miRNAs) take part in cell proliferation, differentiation, apoptosis, and fat metabolism through binding seed sites of targeting mRNAs. The aim of this study was to isolate and identify bovine preadipocytes and screen miRNAs associated with adipogenesis. Bovine preadipocytes were isolated from subcutaneous fatty tissue and induced to differentiate into adipocytes. Verification of preadipocytes and adipocytes was performed by qRT-PCR (real-time quantitative reverse transcription PCR), Oil Red O staining, and immunofluorescence staining. Total RNA was extracted for small RNA sequencing. The sequencing data showed that 131 miRNAs were highly expressed in adipocytes, and 119 miRNAs were highly expressed in preadipocytes. Stem–loop qPCR (stem–loop quantitative real-time PCR) results showed that the expression patterns of 11 miRNAs were consistent with the sequencing results (miR-149-5p, miR-24-3p, miR-199a-5p, miR-33a, etc.). According to KEGG pathway and Gene Ontology (GO) analyses, multiple predicted target genes were associated with lipid metabolism. In summary, this study provides a protocol of isolating bovine preadipocytes and screening various differently expressed miRNAs during preadipocyte differentiation.

Comparison of carcass and meat quality traits between lean and fat Pekin ducks

OBJECTIVE

According to market demand, meat duck breeding mainly includes 2 breeding directions: lean Pekin duck (LPD) and fat Pekin duck (FPD). The aim of the present study was to compare carcass and meat quality traits between 2 strains, and to provide basic data for guidelines of processing and meat quality improvement.

METHODS

A total of 62 female Pekin ducks (32 LPDs and 30 FPDs) were slaughtered at the age of 42 days. The live body weight and carcass traits were measured and calculated. Physical properties of breast muscle were determined by texture analyzer and muscle fibers were measured by paraffin sections. The content of inosine monophosphate (IMP), intramuscular fat (IMF) and fatty acids composition were measured by high-performance liquid chromatography, Soxhlet extraction method and automated gas chromatography respectively.

RESULTS

The results showed that the bodyweight of LPDs was higher than that of FPDs. FPDs were significantly higher than LPDs in subcutaneous fat thickness, subcutaneous fat weight, subcutaneous fat percentage, abdominal fat percentage and abdominal fat shear force (p<0.01). LPDs were significantly higher than FPDs in breast muscle thickness, breast muscle weight, breast muscle rate and breast muscle shear force (p<0.01). The muscle fiber average area and fiber diameter of LPDs were significantly higher than those of FPDs (p<0.01). The muscle fiber density of LPDs was significantly lower than that of FPDs (p<0.01). The IMF of LPDs in the breast muscle was significantly higher than that in the FPDs (p<0.01). There was no significant difference between the 2 strains in IMP content (p>0.05). The polyunsaturated fatty acid content of LPDs was significantly higher than that of FPDs (p<0.01), and FPDs had higher saturated fatty acid and monounsaturated fatty acid levels (p<0.05).

CONCLUSION

Long-term breeding work resulted in vast differences between the two strains Pekin ducks. This study provides a reference for differences between LPD and FPD that manifest as a result of long-term selection.

Discovery of genomic variations by whole-genome resequencing of the North American Araucana chicken

Gallus gallus (chicken) is phenotypically diverse, with over 60 recognized breeds, among the myriad species within the Aves lineage. Domestic chickens have been under artificial selection by humans for thousands of years for agricultural purposes. The North American Araucana (NAA) breed arose as a cross between the Chilean “Collonocas” that laid blue eggs and was rumpless and the “Quetros” that had unusual tufts but with tail. NAAs were introduced from South America in the 1940s and have been kept as show birds by enthusiasts since then due to several distinctive traits: laying eggs with blue eggshells, characteristic ear-tufts, a pea comb, and rumplessness. The population has maintained variants for clean-faced and tufted, as well as tailed and rumplessness traits making it advantageous for genetic studies. Genome resequencing of six NAA chickens with a mixture of these traits was done to 71-fold coverage using Illumina HiSeq 2000 paired-end reads. Trimmed and concordant reads were mapped to the Gallus_gallus-5.0 reference genome (galGal5), generated from a female Red Junglefowl (UCD001). To identify candidate genes that are associated with traits of the NAA, their genome was compared with the Korean Araucana, Korean Domestic and White Leghorn breeds. Genomic regions with significantly reduced levels of heterogeneity were detected on five different chromosomes in NAA. The sequence data generated confirm the identity of variants responsible for the blue eggshells, pea comb, and rumplessness traits of NAA and propose one for ear-tufts.