A short insertion mutation disrupts genesis of miR-16 and causes increased body weight in domesticated chicken

Identification of key module and hub genes affecting broiler body weight through weighted gene co-expression network analysis

Body weight (BW) is an important economic trait in chickens. The hypothalamus serves as a central regulator of appetite and energy balance, and extensive research has demonstrated its pivotal role in regulating BW. However, the molecular network of the hypothalamus regulating BW traits in chickens needs to be further illuminated. In the present study, 200 1-day-old male 817 broilers were reared to 50 d of age, and BW were recorded. 20 birds with the lowest BW were classified as the low body weight group (L-BWG), and 20 birds with the highest BW were classified as the high body weight group (H-BWG). 18 hypothalamic tissue samples were collected, including 5 from the L-BWG, 5 from the H-BWG, and 8 from the middle weight range, and were analyzed using RNA-seq and weighted gene co-expression network analysis (WGCNA). Among the 18 RNA-seq samples, 5 samples from the L-BWG and 5 from the H-BWG were selected for differential expression gene analysis. Compared with the L-BWG, 195 and 1,241 genes were upregulated and downregulated in the H-BWG, respectively. The WGCNA analysis classified all co-expressed genes in the hypothalamus of 817 broilers into 20 modules. Among these modules, the pink module was identified as significantly negatively (r = -0.81, P = 4×10-5) associated with BW. Furthermore, several genes, including Wnt family member 6 (WNT6), growth differentiation factor 11 (GDF11), bone morphogenetic protein 4 (BMP4), and erb-b2 receptor tyrosine kinase 4 (ERBB4), involved in "regulation of developmental process" and "response to growth factor," were identified as hub genes that contribute to the regulation of BW. These results provide valuable information for further understanding of the gene expression and regulation affecting BW traits and will contribute to the molecular breeding of chickens in the future.

Host genetics and gut microbiota synergistically regulate feed utilization in egg-type chickens

BACKGROUND

Feed efficiency is a crucial economic trait in poultry industry. Both host genetics and gut microbiota influence feed efficiency. However, the associations between gut microbiota and host genetics, as well as their combined contributions to feed efficiency in laying hens during the late laying period, remain largely unclear.

METHODS

In total, 686 laying hens were used for whole-genome resequencing and liver transcriptome sequencing. 16S rRNA gene sequencing was conducted on gut chyme (duodenum, jejunum, ileum, and cecum) and fecal samples from 705 individuals. Bioinformatic analysis was performed by integrating the genome, transcriptome, and microbiome to screen for key genetic variations, genes, and gut microbiota associated with feed efficiency.

RESULTS

The heritability of feed conversion ratio (FCR) and residual feed intake (RFI) was determined to be 0.28 and 0.48, respectively. The ileal and fecal microbiota accounted for 15% and 10% of the FCR variance, while the jejunal, cecal, and fecal microbiota accounted for 20%, 11%, and 10% of the RFI variance. Through SMR analysis based on summary data from liver eQTL mapping and GWAS, we further identified four protein-coding genes, SUCLA2, TNFSF13B, SERTM1, and MARVELD3, that influence feed efficiency in laying hens. The SUCLA2 and TNFSF13B genes were significantly associated with SNP 1:25664581 and SNP rs312433097, respectively. SERTM1 showed significant associations with rs730958360 and 1:33542680 and is a potential causal gene associated with the abundance of Corynebacteriaceae in feces. MARVELD3 was significantly associated with the 1:135348198 and was significantly correlated with the abundance of Enterococcus in ileum. Specifically, a lower abundance of Enterococcus in ileum and a higher abundance of Corynebacteriaceae in feces were associated with better feed efficiency.

CONCLUSIONS

This study confirms that both host genetics and gut microbiota can drive variations in feed efficiency. A small portion of the gut microbiota often interacts with host genes, collectively enhancing feed efficiency. Therefore, targeting both the gut microbiota and host genetic variation by supporting more efficient taxa and selective breeding could improve feed efficiency in laying hens during the late laying period.

Regulatory Effects of 198-bp Structural Variants in the GSTA2 Promoter Region on Adipogenesis in Chickens

Molecular breeding accelerates animal breeding and improves efficiency by utilizing genetic mutations. Structural variations (SVs), a significant source of genetic mutations, have a greater impact on phenotypic variation than SNPs. Understanding SV functional mechanisms and obtaining precise information are crucial for molecular breeding. In this study, association analysis revealed significant correlations between 198-bp SVs in the GSTA2 promoter region and abdominal fat weight, intramuscular fat content, and subcutaneous fat thickness in chickens. High expression of GSTA2 in adipose tissue was positively correlated with the abdominal fat percentage, and different genotypes of GSTA2 exhibited varied expression patterns in the liver. The 198-bp SVs regulate GSTA2 expression by binding to different transcription factors. Overexpression of GSTA2 promoted preadipocyte proliferation and differentiation, while interference had the opposite effect. Mechanistically, the 198-bp fragment contains binding sites for transcription factors such as C/EBPα that regulate GSTA2 expression and fat synthesis. These SVs are significantly associated with chicken fat traits, positively influencing preadipocyte development by regulating cell proliferation and differentiation. Our work provides compelling evidence for the use of 198-bp SVs in the GSTA2 promoter region as molecular markers for poultry breeding and offers new insights into the pivotal role of the GSTA2 gene in fat generation.

The roles of miRNAs in adult skeletal muscle satellite cells

Satellite cells are bona fide muscle stem cells that are indispensable for successful post-natal muscle growth and regeneration after severe injury. These cells also participate in adult muscle adaptation in several capacities. MicroRNAs (miRNAs) are post-transcriptional regulators of mRNA that are implicated in several aspects of stem cell function. There is evidence to suggest that miRNAs affect satellite cell behavior in vivo during development and myogenic progenitor behavior in vitro, but the role of miRNAs in adult skeletal muscle satellite cells is less studied. In this review, we provide evidence for how miRNAs control satellite cell function with emphasis on satellite cells of adult skeletal muscle in vivo. We first outline how miRNAs are indispensable for satellite cell viability and control the phases of myogenesis. Next, we discuss the interplay between miRNAs and myogenic cell redox status, senescence, and communication to other muscle-resident cells during muscle adaptation. Results from recent satellite cell miRNA profiling studies are also summarized. In vitro experiments in primary myogenic cells and cell lines have been invaluable for exploring the influence of miRNAs, but we identify a need for novel genetic tools to further interrogate how miRNAs control satellite cell behavior in adult skeletal muscle in vivo.

MicroRNA control of the myogenic cell transcriptome and proteome: the role of miR-16

MicroRNAs (miRs) control stem cell biology and fate. Ubiquitously expressed and conserved miR-16 was the first miR implicated in tumorigenesis. miR-16 is low in muscle during developmental hypertrophy and regeneration. It is enriched in proliferating myogenic progenitor cells but is repressed during differentiation. The induction of miR-16 blocks myoblast differentiation and myotube formation, whereas knockdown enhances these processes. Despite a central role for miR-16 in myogenic cell biology, how it mediates its potent effects is incompletely defined. In this investigation, global transcriptomic and proteomic analyses after miR-16 knockdown in proliferating C2C12 myoblasts revealed how miR-16 influences myogenic cell fate. Eighteen hours after miR-16 inhibition, ribosomal protein gene expression levels were higher relative to control myoblasts and p53 pathway-related gene abundance was lower. At the protein level at this same time point, miR-16 knockdown globally upregulated tricarboxylic acid (TCA) cycle proteins while downregulating RNA metabolism-related proteins. miR-16 inhibition induced specific proteins associated with myogenic differentiation such as ACTA2, EEF1A2, and OPA1. We extend prior work in hypertrophic muscle tissue and show that miR-16 is lower in mechanically overloaded muscle in vivo. Our data collectively point to how miR-16 is implicated in aspects of myogenic cell differentiation. A deeper understanding of the role of miR-16 in myogenic cells has consequences for muscle developmental growth, exercise-induced hypertrophy, and regenerative repair after injury, all of which involve myogenic progenitors.

Molecular Characterization, Expression Profile, and A 21-bp Indel within the ASB9 Gene and Its Associations with Chicken Production Traits

A growing number of studies have shown that members of the ankyrin repeat and suppressors of cytokine signaling (SOCS) box-containing protein (ASB) family are extensively involved in biological processes such as cell growth, tissue development, insulin signaling, ubiquitination, protein degradation, and skeletal muscle membrane protein formation, while the specific biological role of ankyrin-repeat and SOCS box protein 9 (ASB9) remains unclear. In this study, a 21 bp indel in the intron of ASB9 was identified for the first time in 2641 individuals from 11 different breeds and an F₂ resource population, and differences were observed among individuals with different genotypes (II, ID, and DD). An association study of a cross-designed F₂ resource population revealed that the 21-bp indel was significantly related to growth and carcass traits. The significantly associated growth traits were body weight (BW) at 4, 6, 8, 10, and 12 weeks of age; sternal length (SL) at 4, 8, and 12 weeks of age; body slope length (BSL) at 4, 8, and 12 weeks of age; shank girth (SG) at 4 and 12 weeks of age; tibia length (TL) at 12 weeks of age; and pelvic width (PW) at 4 weeks of age (p < 0.05). This indel was also significantly correlated with carcass traits including semievisceration weight (SEW), evisceration weight (EW), claw weight (CLW), breast muscle weight (BMW), leg weight (LeW), leg muscle weight (LMW), claw rate (CLR), and shedding weight (ShW) (p < 0.05). In commercial broilers, the II genotype was the dominant genotype and underwent extensive selection. Interestingly, the ASB9 gene was expressed at significantly higher levels in the leg muscles of Arbor Acres broilers than those of Lushi chickens, while the opposite was true for the breast muscles. In summary, the 21-bp indel in the ASB9 gene significantly influenced the expression of the ASB9 gene in muscle tissue and was associated with multiple growth and carcass traits in the F₂ resource population. These findings suggested that the 21-bp indel within the ASB9 gene could be used in marker-assisted selection breeding for traits related to chicken growth.

QPCTL Affects the Daily Weight Gain of the F2 Population and Regulates Myogenic Cell Proliferation and Differentiation in Chickens

Molecular breeding can accelerate the process of animal breeding and improve the breeding efficiency. To date, many Indel molecular markers have been identified in livestock and poultry, but how Indels affect economic traits is not well understood. For molecular breeding, it is crucial to reveal the mechanism of action of Indels and to provide more accurate information. The purpose of this study was to investigate how the 52/224-bp multiallelic Indels of the chicken QPCTL promoter area affect the daily weight gain of chickens and the potential regulatory mechanism of the QPCTL gene. The analysis was conducted by association analysis, qPCR, dual-fluorescence assay and Western blotting. The results showed that Indels in the QPCTL promoter region were significantly associated with the daily weight gain in chickens and that QPCTL expression showed a decreasing trend in embryonic breast muscle tissues. Furthermore, QPCTL expression was significantly higher in breast muscle tissues of the AC genotype than in those of the AB and BB genotypes. Based on the transcriptional activity results, the pGL3-C vector produced more luciferase activity than pGL3-A and pGL3-B. In addition, overexpression of QPCTL promoted chicken primary myoblast (CPM) proliferation and inhibited differentiation. The results of this study suggest that Indels in the promoter region of the QPCTL gene may regulate the proliferation and differentiation of CPMs by affecting the expression of QPCTL, which ultimately affects the growth rate of chickens. These Indels have important value for the molecular breeding of chickens, and QPCTL can be used as a candidate gene to regulate and improve chicken growth and development.

Embryonic transcriptome unravels mechanisms and pathways underlying embryonic development with respect to muscle growth, egg production, and plumage formation in native and broiler chickens

Background: Muscle development, egg production, and plumage colors are different between native and broiler chickens. The study was designed to investigate why improved Aseel (PD4) is colorful, stronger, and grew slowly compared with the control broiler (CB).

Methods: A microarray was conducted using the 7th-day embryo (7EB) and 18th-day thigh muscle (18TM) of improved Aseel and broiler, respectively. Also, we have selected 24 Gallus gallus candidate reference genes from NCBI, and total RNA was isolated from the broiler, improved Aseel embryo tissues, and their expression profiles were studied by real-time quantitative PCR (qPCR). Furthermore, microarray data were validated with qPCR using improved Aseel and broiler embryo tissues.

Results: In the differential transcripts screening, all the transcripts obtained by microarray of slow and fast growth groups were screened by fold change ≥ 1 and false discovery rate (FDR) ≤ 0.05. In total, 8,069 transcripts were differentially expressed between the 7EB and 18TM of PD4 compared to the CB. A further analysis showed that a high number of transcripts are differentially regulated in the 7EB of PD4 (6,896) and fewer transcripts are differentially regulated (1,173) in the 18TM of PD4 compared to the CB. On the 7th^- and 18th-day PD4 embryos, 3,890, 3,006, 745, and 428 transcripts were up- and downregulated, respectively. The commonly up- and downregulated transcripts are 91 and 44 between the 7th- and 18th-day of embryos. In addition, the best housekeeping gene was identified. Furthermore, we validated the differentially expressed genes (DEGs) related to muscle growth, myostatin signaling and development, and fatty acid metabolism genes in PD4 and CB embryo tissues by qPCR, and the results correlated with microarray expression data.

Conclusion: Our study identified DEGs that regulate the myostatin signaling and differentiation pathway; glycolysis and gluconeogenesis; fatty acid metabolism; Jak-STAT, mTOR, and TGF-β signaling pathways; tryptophan metabolism; and PI3K-Akt signaling pathways in PD4. The results revealed that the gene expression architecture is present in the improved Aseel exhibiting embryo growth that will help improve muscle development, differentiation, egg production, protein synthesis, and plumage formation in PD4 native chickens. Our findings may be used as a model for improving the growth in Aseel as well as optimizing the growth in the broiler.

miR-30a-3p can inhibit the proliferation and promote the differentiation of chicken primary myoblasts

1. Chicken muscle is an important factor in meat quality and its development is controlled by a complex regulatory network.2. The following study examined the expression of miR-30a-3p in Gushi chicken breast muscle tissue and found that it was differentially expressed at different embryonic stages, reaching a peak in the 14-day-old embryo (E14).3. The effect of miR-30a-3p on chicken primary myoblasts (CPMs) was explored. Results from both cell counting kit-8 (CCK-8) and 5-ethynyl-2'-deoxyuridine (EdU) showed that this can inhibit the proliferation of myoblasts, and through cell cycle experiments, the inhibition of myoblast proliferation was found, which may be due to G0/G1 arrest in the cell cycle.4. The effect of miR-30a-3p on the differentiation of myoblasts was studied. The results showed that miR-30a-3p can promote the expression of MYOD, myogenin (MYOG), and myosin heavy chain (MYHC) genes to promote the differentiation of myoblasts. Through MYHC protein immunofluorescence experiments, it was found that miR-30a-3p can effectively increase the area of myotubes.5. Finally, mRNA transcriptome data was analysed, which showed that miR-30a-3p has 51 potential target genes. Among them, forkhead box O3 (FOXO3), ankyrin repeat domain 1 (ANKRD1), and insulin-induced 1 (INSIG1) genes were differentially expressed at different developmental stages and were enriched in Gene Ontology (GO) terms, such as cell differentiation and cellular developmental process. The data showed that tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein gamma (YWHAG), BUB1 mitotic checkpoint serine/threonine kinase (BUB1), and growth arrest and DNA damage-inducible 45 (GADD45) genes were enriched in the cell cycle pathway.6. It can be speculated that miR-30a-3p plays roles through these genes in myoblast development. This research provides information for further improving knowledge of the chicken muscle development regulation network.

Genome-Wide Association Studies Provide Insight Into the Genetic Determination for Hyperpigmentation of the Visceral Peritoneum in Broilers

Hyperpigmentation of the visceral peritoneum (HVP) has been becoming one of the most challenging problems in yellow-feathered chicken production, which seriously affected chicken carcass quality traits. Detecting which genes dominantly impact pigmentation in the peritoneum tissues is of great benefit to the genetic improvement of HVP. To investigate the genetic mechanism of HVP in yellow-feathered broilers, genome-wide association studies (GWASs) were conducted in the F₂ generation of a cross broiler population with 395 birds. A total of 115,706 single-nucleotide polymorphisms (SNPs) of 122,415 were retained to identify quantitative trait loci (QTL) associated to HVP in chicken. The GWAS results based on the logistic mixed model (LMM) revealed that a narrow genomic location on chromosomes 1 (49.2–51.3 Mb) was significantly associated (p ≤ 4.32 × 10⁻⁷) with HVP, which contained 23 SNP makers related to 14 functional genes (MFNG, POLDIP3, POLR2F, PICK1, PDXP, SGSM3, RANGAP1, MYH9, RPL3, GALP3, LGALS1, MICALL1, ATF4, and CYP2D6). Four highly associated (p < 10⁻⁵) haplotype blocks of 0.80 kb (two SNPs), 0.06 kb (two SNPs), 0.95 kb (two SNPs), and 0.03 kb (two SNPs) were identified with two, two, four, and four distinct haplotypes, respectively. As a melanoma-associated gene, CYP2D6 were also possibly involved in the development of HVP occurring in chicken with two significant variations (rs314284996 and rs317955795) in the promoter regions. Further tests revealed that the expression of CYP2D6 was obviously higher in the visceral peritoneum tissue of chicken with HVP than that in the normal group (p < 0.05). Our results provide a novel clue to understand the genetic mechanism of HVP generation in chicken, and the mapped QTL or candidate genes might serve for genomic selection to improve carcass quality in the yellow-feathered chicken industry.

Genetic effect of an InDel in the promoter region of the NUDT15 and its effect on myoblast proliferation in chickens

Background

Molecular breeding accelerates the speed of animal breeding. Screening molecular markers that can affect economic traits through genome-wide association studies (GWAS) can provide a theoretical basis for molecular breeding. At present, a large number of molecular markers have been screened in poultry research, but few reports on how molecular markers affect economic traits exist. It is particularly important to reveal the action mechanisms of molecular markers, which can provide more accurate information for molecular breeding.

Results

The aim of this study was to investigate the relationships between two indels (NUDT15-indel-2777 and NUDT15-indel-1673) in the promoter region of NUDT15 and growth and carcass traits in chickens and to explore the regulatory mechanism of NUDT15. Significant differences were found in genotype and allele frequencies among commercial broilers, commercial laying hens and dual-purpose chickens. The results of association analyses showed that these two indel loci could significantly affect growth traits, such as body weight, and carcass traits. Tissue expression profiling at E12 showed that the expression of NUDT15 was significantly higher in skeletal muscle, and time-expression profiling of leg muscle showed that the expression of NUDT15 in myoblasts was significantly higher in the E10 and E12 proliferation stages than in other stages. Promoter activity analysis showed that pro-1673-I and pro-1673-D significantly inhibited promoter activity, and the promoter activity of pro-1673-D was significantly lower than that of pro-1673-I. In addition, when NUDT15 was overexpressed or underwent interference in chicken primary myoblasts (CPMs), NUDT15 could inhibit the proliferation of CPMs.

Conclusion

The results suggest that the studied indels in the promoter region of NUDT15 may regulate the proliferation of CPMs by affecting NUDT15 expression, ultimately affecting the growth and carcass traits of chickens. These indel polymorphisms may be used together as molecular markers for improving economic traits in chickens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08362-6.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

A 104-bp Structural Variation of the ADPRHL1 Gene Is Associated With Growth Traits in Chickens

Analyzing marker-assisted breeding is an important method utilized in modern molecular breeding. Recent studies have determined that a large number of molecular markers appear to explain the impact of "lost heritability" on human height. Therefore, it is necessary to locate molecular marker sites in poultry and investigate the possible molecular mechanisms governing their effects. In this study, we found a 104-bp insertion/deletion polymorphism in the 5'UTR of the ADPRHL1 gene through resequencing. In cross-designed F₂ resource groups, the indel was significantly associated with weight at 0, 2, 4, 6, and 10 weeks and a number of other traits [carcass weight (CW), semi-evisceration weight (SEW), evisceration weight (EW), claw weight (CLW), wings weight (DWW), gizzard weight (GW), pancreas weight (PW), chest muscle weight (CMW), leg weight (LW), leg muscle weight (LMW), shedding Weight (SW), liver rate (LR), and leg muscle rate (LMR)] (P < 0.05). In brief, the insertion-insertion (II) genotype was significantly associated with the greatest growth traits and meat quality traits, whereas the values associated with the insertion-deletion (ID) genotype were the lowest in the F₂ reciprocal cross chickens. The mutation sites were genotyped in 4,526 individuals from 12 different chicken breeds and cross-designed F₂ resource groups. The II genotype is the most important genotype in commercial broilers, and the I allele frequency observed in these breeds is relatively high. Deletion mutations tend to be fixed in commercial broilers. However, there is still considerable great potential for breeding in dual-purpose chickens and commercial laying hens. A luciferase reporter assay showed that the II genotype of the ADPRHL1 gene possessed 2.49-fold higher promoter activity than the DD genotype (P < 0.05). We hypothesized that this indel might affect the transcriptional activity of ADPRHL1, thereby affecting the growth traits of chickens. These findings may help to elucidate the function of the ADPRHL1 gene and facilitate enhanced reproduction in the chicken industry.

Centennial Review: Metabolic microRNA - shifting gears in the regulation of metabolic pathways in poultry

Over 20 yr ago, a small noncoding class of RNA termed microRNA (miRNA) that was able to recognize sequences in mRNAs and inhibit their translation was discovered in Caenorhabditis elegans. In the intervening years, miRNA have been discovered in most eukaryotes and are now known to regulate the majority of protein-coding genes. It has been discovered that disruption of miRNA function often leads to the development of pathological conditions. One physiological system under extensive miRNA-mediated regulation is metabolism. Metabolism is one of the most dynamic of biological networks within multiple organs, including the liver, muscle, and adipose tissue, working in concert to respond to ever-changing nutritional cues and energy demands. Therefore, it is not surprising that miRNA regulate virtually all aspects of eukaryotic metabolism and have been linked to metabolic disorders, such as obesity, fatty liver diseases, and diabetes, just to name a few. Chickens, and birds in general, face their own unique metabolic challenges, particularly after hatching, when their metabolism must completely transform from using lipid-rich yolk to carbohydrate-rich feed as fuel in a very short period of time. Furthermore, commercial poultry breeds have undergone extensive selection over the last century for more desirable production traits, which has resulted in numerous metabolic consequences. Here, we review the current knowledge of miRNA-mediated regulation of metabolic development and function in chickens.

Genetic regulation and heritability of miRNA and mRNA expression link to phosphorus utilization and gut microbiome

Improved utilization of phytates and mineral phosphorus (P) in monogastric animals contributes significantly to preserving the finite resource of mineral P and mitigating environmental pollution. In order to identify pathways and to prioritize candidate genes related to P utilization (PU), the genomic heritability of 77 and 80 trait-dependent expressed miRNAs and mRNAs in 482 Japanese quail were estimated and eQTL (expression quantitative trait loci) were detected. In total, 104 miR-eQTL (microRNA expression quantitative traits loci) were associated with SNP markers (false discovery rate less than 10%) including 41 eQTL of eight miRNAs. Similarly, 944 mRNA-eQTL were identified at the 5% False discovery rate threshold, with 573 being cis-eQTL of 36 mRNAs. High heritabilities of miRNA and mRNA expression coincide with highly significant eQTL. Integration of phenotypic data with transcriptome and microbiome data of the same animals revealed genetic regulated mRNA and miRNA transcripts (SMAD3, CAV1, ENNPP6, ATP2B4, miR-148a-3p, miR-146b-5p, miR-16-5p, miR-194, miR-215-5p, miR-199-3p, miR-1388a-3p) and microbes (Candidatus Arthromitus, Enterococcus) that are associated with PU. The results reveal novel insights into the role of mRNAs and miRNAs in host gut tissue functions, which are involved in PU and other related traits, in terms of the genetic regulation and inheritance of their expression and in association with microbiota components.

An in silico approach to identify and prioritize miRNAs target sites polymorphisms in colorectal cancer and obesity

Colorectal cancer (CRC) and obesity are linked clinical entities with a series of complex processes being engaged in their development. MicroRNAs (miRNAs) participate in these processes through regulating CRC and obesity‐related genes. This study aimed to develop an in silico approach to systematically identify and prioritize miRNAs target sites polymorphisms in obesity and CRC. Data from genome‐wide association studies (GWASs) were used to retrieve CRC and obesity‐associated variants. The polymorphisms that were resided in experimentally verified or computationally predicted miRNA target sites were retrieved and prioritized using a range of bioinformatics analyses. We found 6284 CRC and 38931 obesity unique variants. For CRC 33 haplotypes variants in 134 interactions were in miRNA targetome, while for obesity we found more than 935 unique interactions. Functionally prioritized SNPs revealed that, SNPs in 153 obesity and 50 CRC unique interactions were have disruptive effects on miRNA:mRNA integration by changing on target RNA secondary structure. Structural accessibility of target sites were decreased in 418 and 103 unique interactions and increased in 516 and 79 interactions, for obesity and CRC, respectively. The miRNA:mRNA hybrid stability was increased in 127 and 17 unique interactions and decreased in 33 and 24 interactions for the effect of obesity and CRC SNPs, respectively. In this study, seven SNPs with 15 interactions and three SNPs with four interactions were prioritized for obesity and CRC, respectively. These SNPs could be used for future studies for finding potential biomarkers for diagnoses, prognosis, or treatment of CRC and obesity.

Genome-wide association study reveals the genetic determinism of growth traits in a Gushi-Anka F2 chicken population

Chicken growth traits are economically important, but the relevant genetic mechanisms have not yet been elucidated. Herein, we performed a genome-wide association study to identify the variants associated with growth traits. In total, 860 chickens from a Gushi-Anka F₂ resource population were phenotyped for 68 growth and carcass traits, and 768 samples were genotyped based on the genotyping-by-sequencing (GBS) method. Finally, 734 chickens and 321,314 SNPs remained after quality control and removal of the sex chromosomes, and these data were used to carry out a GWAS analysis. A total of 470 significant single-nucleotide polymorphisms (SNPs) for 43 of the 68 traits were detected and mapped on chromosomes (Chr) 1-6, -9, -10, -16, -18, -23, and -27. Of these, the significant SNPs in Chr1, -4, and -27 were found to be associated with more than 10 traits. Multiple traits shared significant SNPs, indicating that the same mutation in the region might have a large effect on multiple growth or carcass traits. Haplotype analysis revealed that SNPs within the candidate region of Chr1 presented a mosaic pattern. The significant SNPs and pathway enrichment analysis revealed that the MLNR, MED4, CAB39L, LDB2, and IGF2BP1 genes could be putative candidate genes for growth and carcass traits. The findings of this study improve our understanding of the genetic mechanisms regulating chicken growth and carcass traits and provide a theoretical basis for chicken breeding programs.

Multiple ancestral haplotypes harboring regulatory mutations cumulatively contribute to a QTL affecting chicken growth traits

In depth studies of quantitative trait loci (QTL) can provide insights to the genetic architectures of complex traits. A major effect QTL at the distal end of chicken chromosome 1 has been associated with growth traits in multiple populations. This locus was fine-mapped in a fifteen-generation chicken advanced intercross population including 1119 birds and explored in further detail using 222 sequenced genomes from 10 high/low body weight chicken stocks. We detected this QTL that, in total, contributed 14.4% of the genetic variance for growth. Further, nine mosaic precise intervals (Kb level) which contain ancestral regulatory variants were fine-mapped and we chose one of them to demonstrate the key regulatory role in the duodenum. This is the first study to break down the detail genetic architectures for the well-known QTL in chicken and provides a good example of the fine-mapping of various of quantitative traits in any species.

A 31-bp indel in the 5' UTR region of GNB1L is significantly associated with chicken body weight and carcass traits

BACKGROUND

G-protein subunit beta 1 like (GNB1L) encodes a G-protein beta-subunit-like polypeptide. Chicken GNB1L is upregulated in the breast muscle of high feed efficiency chickens, and its expression is 1.52-fold that in low feed efficiency chickens. However, no report has described the effects of GNB1L indels on the chicken carcass and growth traits.

RESULTS

This study identified a 31-bp indel in the 5' untranslated region (UTR) of GNB1L and elucidated the effect of this gene mutation on the carcass and growth traits in chickens. The 31-bp indel showed a highly significant association with the body weight at 8 different stages and was significantly correlated with daily gains at 0 to 4 weeks and 4 to 8 weeks. Similarly, the mutation was significantly associated with small intestine length, breast width, breast depth and breast muscle weight. Moreover, DD and ID were superior genotypes for chicken growth and carcass traits.

CONCLUSIONS

These results show that the 31-bp indel of GNB1L significantly affects chicken body weight and carcass traits and can serve as a candidate molecular marker for chicken genetics and breeding programs.

Endocrine and genetic factors affecting egg laying performance in chickens: a review

1. Egg-laying performance reflects the overall reproductive performance of breeding hens. The genetic traits for egg-laying performance have low or medium heritability, and, depending on the period involved, usually ranges from 0.16 to 0.64. Egg-laying in chickens is regulated by a combination of environmental, endocrine and genetic factors. 2. The main endocrine factors that regulate egg-laying are gonadotropin-releasing hormone (GnRH), prolactin (PRL), follicle-stimulating hormone (FSH) and luteinising hormone (LH). 3. In the last three decades, many studies have explored this aspect at a molecular genetic level. Recent studies identified 31 reproductive hormone-based candidate genes that were significantly associated with egg-laying performance. With the development of genome-sequencing technology, 64 new candidate genes and 108 single nucleotide polymorphisms (SNPs) related to egg-laying performance have been found using genome-wide association studies (GWAS), providing novel insights into the molecular genetic mechanisms governing egg production. At the same time, microRNAs that regulate genes responsible for egg-laying in chickens were reviewed. 4. Research on endocrinological and genetic factors affecting egg-laying performance will greatly improve the reproductive performance of chickens and promote the protection, development, and utilisation of poultry. This review summarises studies on the endocrine and genetic factors of egg-laying performance in chickens from 1972 to 2019.

A 51 bp indel polymorphism within the PTH1R gene is significantly associated with chicken growth and carcass traits

Parathyroid hormone (PTH) is a crucial regulator of calcium homeostasis and bone remodeling, and the parathyroid hormone 1 receptor (PTH1R) belongs to a class II G-protein-coupled receptor. PTH activates PTH1R, which mediates catabolic and anabolic processes in the skeleton. However, the functional mechanism of PTH1R has not been thoroughly elucidated in organisms. This study identified a 51 bp indel mutation in the first intron of the PTH1R gene and elucidated the effect of this gene mutation on the growth and carcass traits in chickens. The results indicated that the 51 bp indel was significantly associated with subcutaneous fat thickness, abdominal fat weight, body weight and daily gain over 4-8 weeks. Furthermore, we found that PTH1R gene expression was highest in the kidney and liver tissues, and it showed a trend of decreasing in leg and breast muscle tissues at different embryonic stages. In addition, we examined the expression of the three genotypes of the PTH1R gene in the liver, breast muscle and abdominal fat and found that the II genotype was significantly higher than the DD and ID genotypes. In summary, these findings suggest that the PTH1R gene can serve as a potential molecular marker for chicken breeding.

Antagonistic effect of selenium on lead-induced neutrophil apoptosis in chickens via miR-16-5p targeting of PiK3R1 and IGF1R

Environmental contamination by heavy metals, such as lead (Pb), can lead to severe immune dysfunction. MicroRNAs (miRNAs) are involved in regulating immunity. Whether Pb can regulate neutrophil apoptosis through miRNA, and whether selenium (Se) can antagonize this response are still unknown. We treated neutrophils with 12.5 μM (CH₃OO)₂Pb and 1 μM Na₂SeO₃ for 3 h, after which apoptosis was evaluated using acrideine orange/ethidium bromide (AO/EB) dual fluorescent staining and flow cytometry. The results showed that neutrophil apoptosis was significantly increased following Pb exposure, and that this response was prevented upon Se addition. Pb up-regulates miR-16-5p and leads to the subsequent down-regulation of the target genes phosphoinositide-3-kinase regulatory subunit 1 (PiK3R1), insulin-like growth factor 1 receptor (IGF1R), and phosphatidylinositol 3 kinase (Pi3K)-protein kinase B (AKT), followed by activation of the tumor protein P53 (P53)-B-cell lymphoma-2 (Bcl-2)/Bcl-2-Associated X protein (Bax)-cytochrome c (Cytc)-Caspase 9 (mitochondrial apoptotic pathway) and the tumor necrosis factor receptor superfamily member 6 (Fas)-Fas-associated death domain protein (Fadd)-Caspase 8 (death receptor pathway). Pb also triggered oxidative stress and indirectly activated the mitochondrial apoptotic pathway. We conclude that miR-16-5p plays a key role in the apoptosis of neutrophils exposed to Pb by down-regulating the expression of PiK3R1 and IGFR1, thereby activating the mitochondrial apoptotic pathway and death receptor pathway. Se can prevent Pb-induced apoptosis.

Molecular characterization and an 80-bp indel polymorphism within the prolactin receptor (PRLR) gene and its associations with chicken growth and carcass traits

The prolactin receptor (PRLR), a type I cytokine receptor, must bind prolactin (PRL) to act on target cells to mediate various physiological functions, including reproduction and lactation. This study identified an 80-bp insertion/deletion (indel) polymorphism in the 3'-untranslated region (3'-UTR) of the chicken PRLR gene in 3736 individuals from 15 breeds and analyzed its associations with growth and carcass traits in an F₂ resource population. The results of the association analysis indicated that the 80-bp indel polymorphism was significantly (P < 0.05) or very significantly (P < 0.01) associated with multiple growth and carcass traits, such as body weight, leg weight, and shank length. In addition, we found that during the breeding process of commercial laying hens and commercial broilers, the 80-bp indel locus was artificially selected for the II genotype. Together, our findings reveal that this 80-bp indel polymorphism has potential as a new molecular marker for marker-assisted selection of chicken growth and carcass traits.

Evaluation of a Chicken 600K SNP genotyping array in non-model species of grouse

The use of single nucleotide polymorphism (SNP) arrays to generate large SNP datasets for comparison purposes have recently become an attractive alternative to other genotyping methods. Although most SNP arrays were originally developed for domestic organisms, they can be effectively applied to wild relatives to obtain large panels of SNPs. In this study, we tested the cross-species application of the Affymetrix 600K Chicken SNP array in five species of North American prairie grouse (Centrocercus and Tympanuchus genera). Two individuals were genotyped per species for a total of ten samples. A high proportion (91%) of the total 580 961 SNPs were genotyped in at least one individual (73–76% SNPs genotyped per species). Principal component analysis with autosomal SNPs separated the two genera, but failed to clearly distinguish species within genera. Gene ontology analysis identified a set of genes related to morphogenesis and development (including genes involved in feather development), which may be primarily responsible for large phenotypic differences between Centrocercus and Tympanuchus grouse. Our study provided evidence for successful cross-species application of the chicken SNP array in grouse which diverged ca. 37 mya from the chicken lineage. As far as we are aware, this is the first reported application of a SNP array in non-passerine birds, and it demonstrates the feasibility of using commercial SNP arrays in research on non-model bird species.

Bidirectional Selection for Body Weight on Standing Genetic Variation in a Chicken Model

Experimental populations of model organisms provide valuable opportunities to unravel the genomic impact of selection in a controlled system. The Virginia body weight chicken lines represent a unique resource to investigate signatures of selection in a system where long-term, single-trait, bidirectional selection has been carried out for more than 60 generations. At 55 generations of divergent selection, earlier analyses of pooled genome resequencing data from these lines revealed that 14.2% of the genome showed extreme differentiation between the selected lines, contained within 395 genomic regions. Here, we report more detailed analyses of these data exploring the regions displaying within- and between-line genomic signatures of the bidirectional selection applied in these lines. Despite the strict selection regime for opposite extremes in body weight, this did not result in opposite genomic signatures between the lines. The lines often displayed a duality of the sweep signatures, where an extended region of homozygosity in one line, in contrast to mosaic pattern of heterozygosity in the other line. These haplotype mosaics consisted of short, distinct haploblocks of variable between-line divergence, likely the results of a complex demographic history involving bottlenecks, introgressions and moderate inbreeding. We demonstrate this using the example of complex haplotype mosaicism in the growth1 QTL. These mosaics represent the standing genetic variation available at the onset of selection in the founder population. Selection on standing genetic variation can thus result in different signatures depending on the intensity and direction of selection.

Upregulated gga-miR-16-5p Inhibits the Proliferation Cycle and Promotes the Apoptosis of MG-Infected DF-1 Cells by Repressing PIK3R1-Mediated the PI3K/Akt/NF-κB Pathway to Exert Anti-Inflammatory Effect

Mycoplasma gallisepticum (MG) mainly infects chickens to initiate chronic respiratory disease (CRD). microRNAs (miRNAs) play vital roles according to previously reported studies. Our previous study showed that gga-miR-16-5p, in MG-infected lungs of chicken embryo, was upregulated by Illumina sequencing. The study aimed to reveal what role gga-miR-16-5p plays in CRD progression. gga-miR-16-5p was upregulated in MG-infected fibroblast cells (DF-1). Phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1) was demonstrated as the target gene of gga-miR-16-5p. Furthermore, PIK3R1 expression was lower in MG-infected groups than it in noninfected controls measured by qPCR. Additionally, overexpressed gga-miR-16-5p could downregulate PIK3R1 and phosphorylated serine/threonine kinase (p-Akt) to express protein, whereas there is an opposite effect on inhibition. Overexpressed gga-miR-16-5p resulted in decreased activity of tumor necrosis factor alpha (TNF-α) and the nuclear factor-kappaB (NF-κB) by qPCR. Furthermore, overexpressed gga-miR-16-5p restricted cell multiplication, cycle progression, and increased apoptosis of MG-infected DF-1 cells, whereas inhibited gga-miR-16-5p led to the opposite effect. Collectively, upregulated gga-miR-16-5p could decrease multiplication, cycle progression, and increase apoptosis of MG-infected DF-1 cells, at least partly through directly targeting PIK3R1 and inhibiting PI3K/Akt/NF-κB pathway to exert an anti-inflammatory effect. Our results will provide more experimental evidence to bring pathogenesis of MG infection to light.

Chicken ZNF764L gene: mRNA expression profile, alternative splicing analysis and association analysis between first exon indel mutation and economic traits

Zinc finger proteins are a class of transcription factors with finger-like domains and have diverse uses in biological processes, including development, differentiation, and metabolism. In this study, we identified the absence of the 24 bp sequence in the third exon of the zinc finger protein 764-like (ZNF764L) gene that lead to the production of two new transcripts, ZNF764L-SV1 and ZNF764L-SV2, and the sum of the expression levels of the two transcripts is approximately equal the total RNA expression level. Temporal and spatial expression showed that ZNF764L had higher expression during the embryonic stage. Moreover, the research study revealed a 22-bp indel mutation in the first exon region of ZNF764L gene. Statistically significant results (P < 0.05) were encountered for this indel for chicken growth and carcass traits, which include birth weight, chest breadth and body slanting length at 4 weeks of age and subcutaneous fat weight and others. Genetic parameter analysis showed that D is the predominant allele in the commercial chicken population. Gene expression for each genotype showed that birds carrying the II allele had a higher expression level than the other genotypes. These findings enrich the understanding of ZNF764L gene function and enhance reproduction in the chicken industry.

An Evolutionary Genomic Perspective on the Breeding of Dwarf Chickens

The evolutionary history for dwarfism in chickens remains an enigma. Herein, we explore the evolution of the Serama, the smallest breed of chicken. Leveraging comparative population genomics, analyses identify several genes that are potentially associated with the growth and development of bones and muscles. These genes, and in particular both POU1F1 and IGF1, are under strong positive selection. Three allopatric dwarf bantams (Serama, Yuanbao, and Daweishan) with different breeding-histories, form distinct clusters and exhibit unique population structures. Parallel genetic mechanisms underlay their variation in body size. These findings provide insights into the multiple and complex pathways, depending on genomic variation, that chicken can take in response to aviculture selection for dwarfism.

Polygenic Basis and Variable Genetic Architectures Contribute to the Complex Nature of Body Weight -A Genome-Wide Study in Four Chinese Indigenous Chicken Breeds

Body weight (BW) is one of the most important economic traits for animal production and breeding, and it has been studied extensively for its phenotype–genotype associations. While mapping studies have mostly aimed at finding as many loci as possible that contributed to the variation in BW, the role of other factors in its genetic architecture, including their frequencies in the population and their interactions, have been largely overlooked. To comprehensively characterized the genetic architecture of BW, we performed a genome-wide association study (GWAS) both at the single-marker and haplotype level on birds from four indigenous Chinese chicken breeds (Chahua, Silkie, Langshan, and Beard), rather than studying crosses between two founder lines. Additionally, samples from two more breeds (Red Junglefowl and Recessive White) were included to better reflect variable genetic characteristics across populations. Six loci were mapped in this study, revealing the polygenic basis underlying BW. Moreover, by further examining the frequencies of the significantly associated haplotypes in each subpopulation and their effect sizes, most of the loci were found to affect BW in the Beard chicken breed alone. Two loci in GGA9 and GGA27, however, had a common effect on BW across subpopulations, showing that different underlying genetic mechanisms contribute to the phenotypic variability. These findings, particularly the variable genetic architectures found in different loci, improve our understanding of the overall genetic contributions to the large variability in BW among Chinese indigenous chicken breeds. These findings thus will have important implications for future chicken breeding.

miR-16 controls myoblast proliferation and apoptosis through directly suppressing Bcl2 and FOXO1 activities

Myogenesis mainly involves several steps including myoblast proliferation, differentiation, apoptosis and fusion. Except for muscle specific regulators, few miRNAs were proved to coordinate this complex process. Here, we reported that miR-16 inhibited myoblast proliferation and promoted myoblast apoptosis by directly targeting Bcl2 and FOXO1. The expression level of miR-16 was significantly decreased in the hypertrophic pectoral muscle compared to the normal pectoral muscle in chicken. In vitro, elevating miR-16 significantly inhibited myoblast proliferation and promoted myoblast apoptosis, resulting in about 11.2% cells arrested in G1 phase and 12.3% apoptotic cells in the early stage. Bioinformatic and biochemical analyses revealed Bcl2 and FOXO1 as direct targets of miR-16. Consist to the effect of miR-16 on myogenesis, specific inhibition of Bcl2 or FOXO1 significantly suppressed myoblast proliferation and induced myoblast apoptosis, indicating that both Bcl2 and FOXO1 contributed to miR-16 regulatory function in myogenesis. Interestingly, FOXO1, as the core target, mediated multiple growth-related pathways induced by miR-16 such as PI3K-AKT-MAPK and PI3K-AKT-mTOR. Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) revealed that 234 annotated genes bound by FOXO1 in the early-differentiated myoblasts, which were significantly enriched in myogenic proliferation, death and hypotrophy. Altogether, we proposed that miR-16 acted as a coordinated mediator to suppress myogenesis in avian through the control of myoblast proliferation and apoptosis. These findings have provided a novel mechanism whereby miR-16 represses Bcl2 and FOXO1 expression to maintain myoblast growth and skeletal muscle mass.