Analysis of four complete linkage sequence variants within a novel lncRNA located in a growth QTL on chromosome 1 related to growth traits in chickens

Genomic analysis of Nigerian indigenous chickens reveals their genetic diversity and adaptation to heat-stress

Indigenous poultry breeds from Africa can survive in harsh tropical environments (such as long arid seasons, excessive rain and humidity, and extreme heat) and are resilient to disease challenges, but they are not productive compared to their commercial counterparts. Their adaptive characteristics are in response to natural selection or to artificial selection for production traits that have left selection signatures in the genome. Identifying these signatures of positive selection can provide insight into the genetic bases of tropical adaptations observed in indigenous poultry and thereby help to develop robust and high-performing breeds for extreme tropical climates. Here, we present the first large-scale whole-genome sequencing analysis of Nigerian indigenous chickens from different agro-climatic conditions, investigating their genetic diversity and adaptation to tropical hot climates (extreme arid and extreme humid conditions). The study shows a large extant genetic diversity but low level of population differentiation. Using different selection signature analyses, several candidate genes for adaptation were detected, especially in relation to thermotolerance and immune response (e.g., cytochrome P450 2B4-like, TSHR, HSF1, CDC37, SFTPB, HIF3A, SLC44A2, and ILF3 genes). These results have important implications for conserving valuable genetic resources and breeding improvement of chickens for thermotolerance.

Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology

BACKGROUND

The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues.

RESULTS

A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network.

CONCLUSIONS

These validated results show significant improvement over current bovine genome annotations.

Association between the copy number variation of CCSER1 gene and growth traits in Chinese Capra hircus (goat) populations

Recently, Coiled-coil serine-rich protein 1 (CCSER1) gene is reported to be related to economic traits in livestock, and become a hotspot. In our study, we detected CCSER1 gene CNV in 693 goats from six breeds (GZB, GZW, AN, BH, HG, TH) by quantitative real-time PCR (qPCR) and the association analysis between the types of CNV and growth traits. Then, CCSER1 gene expression pattern was discovered in seven tissues from NB goats. Our results showed that the CCSER1 gene copy numbers were distributed differently in the aforementioned six breeds. The type of CCSER1 gene CNV was significantly associated with body weight and heart girth traits in GZW goat, in which individuals with deletion type were dominant in body weight trait (P < 0.05), while the normal type individuals were more advantageous in heart girth trait (P < 0.01); and there was a significant association with heart girth in TH goat (P < 0.05), which normal type was the dominant one. The expression profile revealed that CCSER1 gene has the highest level in the lung, followed by the small intestine and heart. In conclusion, our result is dedicated to an in-depth study of the novel CCSER1 gene CNV site and to provide essential information for Chinese goats molecular selective breeding in the future.

JustRNA: a database of plant long noncoding RNA expression profiles and functional network

Long noncoding RNAs (lncRNAs) are regulatory RNAs involved in numerous biological processes. Many plant lncRNAs have been identified, but their regulatory mechanisms remain largely unknown. A resource that enables the investigation of lncRNA activity under various conditions is required because the co-expression between lncRNAs and protein-coding genes may reveal the effects of lncRNAs. This study developed JustRNA, an expression profiling resource for plant lncRNAs. The platform currently contains 1 088 565 lncRNA annotations for 80 plant species. In addition, it includes 3692 RNA-seq samples derived from 825 conditions in six model plants. Functional network reconstruction provides insight into the regulatory roles of lncRNAs. Genomic association analysis and microRNA target prediction can be employed to depict potential interactions with nearby genes and microRNAs, respectively. Subsequent co-expression analysis can be employed to strengthen confidence in the interactions among genes. Chromatin immunoprecipitation sequencing data of transcription factors and histone modifications were integrated into the JustRNA platform to identify the transcriptional regulation of lncRNAs in several plant species. The JustRNA platform provides researchers with valuable insight into the regulatory mechanisms of plant lncRNAs. JustRNA is a free platform that can be accessed at http://JustRNA.itps.ncku.edu.tw.

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.

Genome-wide identification evolution and expression of vestigial-like gene family in chicken

Vestigial-like (Vgll) genes are widespread in vertebrates and play an important role in muscle development. In this study, we used bioinformatics methods to systematically identify the chicken VGLL family in the whole genome and investigated its evolutionary history and gene structure features. Tissue expression spectra combined with real-time PCR data were used to analyze the organizational expression pattern of the genes. Based on the maximum likelihood method, a phylogenetic tree of the VGLL family was constructed, and 94 VGLL genes were identified in 24 breeds, among which four VGLL family genes were identified in the chicken genome. Ten motifs were detected in the VGLL genes, and the analysis of introns combined with gene structure revealed that the family was conserved during evolution. Tissue expression analysis suggested that the expression profiles of the VGLL family genes in 16 tissues differed between LU Shi and AA broilers. In addition, a single gene (VGLL2) showed increased expression in chickens at embryonic days 10-16 and was involved in the growth and development of skeletal muscle in chickens in the embryonic stage. In summary, VGLL genes are involved in chicken muscle growth and development, which provides useful information for subsequent functional studies of VGLL genes.

Two cSNPs sites in the fatty acid-binding protein 4 (FABP4) gene and their association analysis with body measurement data in five Chinese cattle breeds

Single nucleotide polymorphisms (SNPs) include the transition and transversion of a single base. In this study, we found two cSNPs (rs110757796 and rs110652478) in the FABP4 (fatty acid-binding protein 4) gene and tested their population genetic parameters in five Chinese cattle breeds. FABP4 is a key marker molecule for lipid production. It plays a crucial role in the growth and development of animals. Thus, we also analyze the association between the two cSNPs of FABP4 gene and body measurement data of Chinese cattle. Our results were indicated that rs110757796 was significantly associated with the chest width in Chinese cattle groups (p < 0.05). In addition, we tested the spatiotemporal expression profile of the bovine FABP4 gene and effect of genetic variation on its expression. It was found that bovine FABP4 gene has tissue-differential expression. Then, the genetic variation located on the exon has a certain effect on the expression of FABP4 gene in bovine muscle. Overall, our results provide that FABP4 could as a candidate gene to improve the progress of cattle molecular breeding.

RNA sequencing reveals lncRNA-mediated non-mendelian inheritance of feather growth change in chickens

Background

Long non-coding RNAs (lncRNAs) play an essential role in biological processes. However, the expression patterns of lncRNAs that regulate the non-Mendelian inheritance feather phenotypes remain unknown.

Objective

This study aimed to compare the expression profiles of lncRNAs in the follicles of the late-feathering cocks (LC) and late-feathering hens (LH) that followed genetic rules and the early-feathering hen (EH) and early-feathering cock (EC) that did not conform to the genetic laws.

Methods

We performed RNA sequencing and investigated the differentially expressed lncRNAs (DElncRNAs) between the early- and late-feathering chickens, which function by cis-acting or participate in the competing endogenous RNA (ceRNA) network.

Results

A total of 53 upregulated and 43 downregulated lncRNAs were identified in EC vs. LC, and 58 upregulated and 109 downregulated lncRNAs were identified in EH vs. LH. The target mRNAs regulated by lncRNAs in cis were enriched in the pentose phosphate pathway, TGF-β signaling pathway and Jak-STAT signaling pathway in EC vs. LC and were associated with the TGF-β signaling pathway, Wnt signaling pathway, p53 signaling pathway and Jak-STAT signaling pathway in EH vs. LH. In addition, the lncRNA-mediated ceRNA regulatory pathways of hair follicle formation were mainly enriched in the TGF-β signaling pathway, Wnt signaling pathway, melanogenesis, and calcium signaling pathways. The levels of ENSGALG00000047626 were significantly higher in the late-feathering chickens than in the early-feathering chickens, which regulated the expression of SSTR2 by gga-miR-1649-5p.

Conclusion

This study provides a novel molecular mechanism of lncRNA’s response to the feather rate that does not conform to the genetic laws in chickens.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13258-022-01304-2.

The copy number variation of DMBT1 gene effects body traits in two Chinese cattle breeds

Copy number variations (CNVs) belong to mutations in the genome level with loci in the region of genic or intergenic. It is through different effects (such as position effect and dose effect) that influence complex traits and diseases. Deleted in Malignant Brain Tumors 1 (DMBT1) gene is a member of the scavenger receptor cysteine-rich super family. In cattle, this gene has been associated with the susceptibility to bovine tuberculosis. In this study, a new CNV was found in DMBT1 gene of Chinese cattle breeds and tested in two different Chinese cattle breeds (Jiaxian red and Pinan) for frequency distribution analysis. Besides, the body size data such as body length, body height, chest girth, chest width, rump length, and rump girth for Jiaxian (JX) and Pinan (PN) cattle were collected and associated with the newly identified CNV. The CNV was significantly associated with the body length and chest girth of JX cattle, and the rump length of PN cattle (P < 0.05). Furthermore, the expression profile of the DMBT1 gene was tested in calves' tissues and the myoblasts differentiation. It was found that the DMBT1 gene expression was high in tuberculosis susceptible tissues (liver and lungs) at the calf stage and high in myoblast early differentiation. These tests were done using the qPCR method. As the result, the CNV of DMBT1 gene could be used as a candidate marker for bovine growth and health in marker-assisted selection (MAS) breeding.

Copy number variation of the CCDC39 gene is associated with growth traits in Chinese cattle

Background

Copy number variation (CNV) has become an essential part of genetic structural variation. Coiled‐coil domain containing 39 (CCDC39) is a gene that related to the growth and development of organs and tissues. It is identified that it has a CNV region by animal genome resequencing.

Objective

In this study, we detected the phenotypic traits and different distributions of CCDC39 gene copy numbers in five Chinese cattle breeds (Qinchuan (QC) cattle, Yunling (YL) cattle, Xianan (XN) cattle, Pinan (PN) cattle and Jiaxian (JX) cattle).

Methods

Five hundred and six cattle were randomly selected for CNV distribution detection. Blood samples were taken and genomic DNA was extracted. Different tissues were obtained from adult (n = 3) XN cattle, including heart, liver, kidney, skeletal muscle and lung. The genome qPCR experiment was performed with SYBR Green in triplicate. CDNA qPCR was used to detect the expression level of CCDC39 in different tissues and varieties. Using SPSS v20.0 software, the relationship between CCDC39 CNV and the growth traits of PN, XN, QC, NY and YL cattle breeds was analyzed by one‐way analysis of variance (ANOVA).

Results

The results showed that the expression of CCDC39 in lung was higher than that in other tissues. The expression in liver and kidney was similar, but the expression in heart and muscle was less. It can be seen that the duplication type of QC cattle CCDC39 CNV is higher than the deletion or normal in the height at hip cross. The normal type of PN cattle in body length and hip width was better than duplication and deletion (p < 0.05). In XN cattle, the deletion type of CNV had superior growth characteristics in heart girth and cannon bone circumference compared with the duplication type and the normal type (p < 0.05).

Conclusion

The study revealed a significant association between CNV of CCDC39 gene and growth traits in different Chinese cattle breeds.

CNV analysis of VAMP7 gene reveals variation associated with growth traits in Chinese cattle

Copy number variant (CNV), a common genetic polymorphism, is closely related to the phenotypic variation traits of organisms. Vesicle-associated membrane protein 7 gene (VAMP7) codes a protein, which is a member of the SNARE proteins family and plays an important role in the process of intracellular vesicle transport. In this study, a total of four cattle breeds (Yunling cattle, Xianan cattle, Pinan cattle, Jiaxian red cattle) were used to investigate the copy numbers, and we found an association relationship between CNV of VAMP7 gene and growth traits of cattle by SPSS 20.0 software. The results showed that the CNV type of VAMP7 gene in four cattle breeds had the same distribution, Duplication type occupies a dominant position among the four varieties. In Yunling cattle, the Duplication type of VAMP7 is significantly related to the height at the hip cross (p < 0.05), Individuals with Duplication type commonly have less performance on growth and development, which indicates that the Duplication type of the VAMP7 gene may have a negative effect on cattle growth. Individuals with the other two CNV types may become the breeding direction of the VAMP7 gene. This study provided a new perspective and basic material for the molecular genetics of the CNV of the VAMP7 gene, and also promoted the breeding progress of Chinese local cattle.

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.

Transcriptome analysis reveals the potential roles of long non-coding RNAs in feed efficiency of chicken

Feed efficiency is an important economic trait and reduces the production costs per unit of animal product. Up to now, few studies have conducted transcriptome profiling of liver tissue in feed efficiency-divergent chickens (Ross vs native breeds). Also, molecular mechanisms contributing to differences in feed efficiency are not fully understood, especially in terms of long non-coding RNAs (lncRNAs). Hence, transcriptome profiles of liver tissue in commercial and native chicken breeds were analyzed. RNA-Seq data along with bioinformatics approaches were applied and a series of lncRNAs and target genes were identified. Furthermore, protein-protein interaction network construction, co-expression analysis, co-localization analysis of QTLs and functional enrichment analysis were used to functionally annotate the identified lncRNAs. In total, 2,290 lncRNAs were found (including 1,110 annotated, 593 known and 587 novel), of which 53 (including 39 known and 14 novel), were identified as differentially expressed genes between two breeds. The expression profile of lncRNAs was validated by RT-qPCR. The identified novel lncRNAs showed a number of characteristics similar to those of known lncRNAs. Target prediction analysis showed that these lncRNAs have the potential to act in cis or trans mode. Functional enrichment analysis of the predicted target genes revealed that they might affect the differences in feed efficiency of chicken by modulating genes associated with lipid metabolism, carbohydrate metabolism, growth, energy homeostasis and glucose metabolism. Some gene members of significant modules in the constructed co-expression networks were reported as important genes related to feed efficiency. Co-localization analysis of QTLs related to feed efficiency and the identified lncRNAs suggested several candidates to be involved in residual feed intake. The findings of this study provided valuable resources to further clarify the genetic basis of regulation of feed efficiency in chicken from the perspective of lncRNAs.

Transcriptome sequencing analysis for the identification of stable lncRNAs associated with bovine Staphylococcus aureus mastitis

Background

Staphylococcus aureus (S. aureus) mastitis is one of the most difficult diseases to treat in lactating dairy cows worldwide. S. aureus with different lineages leads to different host immune responses. Long non-coding RNAs (lncRNAs) are reported to be widely involved in the progress of inflammation. However, no research has identified stable lncRNAs among different S. aureus strain infections. In addition, folic acid (FA) can effectively reduce inflammation, and whether the inflammatory response caused by S. aureus can be reduced by FA remains to be explored.

Methods

lncRNA transcripts were identified from Holstein mammary gland tissues infected with different concentrations of S. aureus (in vivo) and mammary alveolar cells (Mac-T cells, in vitro) challenged with different S. aureus strains. Differentially expressed (DE) lncRNAs were evaluated, and stable DE lncRNAs were identified in vivo and in vitro. On the basis of the gene sequence conservation and function conservation across species, key lncRNAs with the function of potentially immune regulation were retained for further analysis. The function of FA on inflammation induced by S. aureus challenge was also investigated. Then, the association analysis between these keys lncRNA transcripts and hematological parameters (HPs) was carried out. Lastly, the knockdown and overexpression of the important lncRNA were performed to validate the gene function on the regulation of cell immune response.

Results

Linear regression analysis showed a significant correlation between the expression levels of lncRNA shared by mammary tissue and Mac-T cells (P < 0.001, R² = 0.3517). lncRNAs PRANCR and TNK2–AS1 could be regarded as stable markers associated with bovine S. aureus mastitis. Several HPs could be influenced by SNPs around lncRNAs PRANCR and TNK2–AS1. The results of gene function validation showed PRANCR regulates the mRNA expression of SELPLG and ITGB2 within the S. aureus infection pathway and the Mac-T cells apoptosis. In addition, FA regulated the expression change of DE lncRNA involved in toxin metabolism and inflammation to fight against S. aureus infection.

Conclusions

The remarkable association between SNPs around these two lncRNAs and partial HP indicates the potentially important role of PRANCR and TNK2–AS1 in immune regulation. Stable DE lncRNAs PRANCR and TNK2–AS1 can be regarded as potential targets for the prevention of bovine S. aureus mastitis. FA supplementation can reduce the negative effect of S. aureus challenge by regulating the expression of lncRNAs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00639-2.

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.

Systematic screening of long intergenic noncoding RNAs expressed during chicken embryogenesis

Long noncoding RNAs (lncRNAs) have emerged as important regulators of many biological processes, including embryogenesis and development. To provide a systematic analysis of lncRNAs expressed during chicken embryogenesis, we used Iso-Seq and RNA-Seq to identify potential lncRNAs at embryonic stages from d 1 to d 8 of incubation: sequential stages covering gastrulation, somitogenesis, and organogenesis. The data characterized an expanded landscape of lncRNAs, yielding 45,410 distinct lncRNAs (31,282 genes). Amongst these, a set of 13,141 filtered intergenic lncRNAs (lincRNAs) transcribed from 9803 lincRNA gene loci, of which, 66.5% were novel, were further analyzed. These lincRNAs were found to share many characteristics with mammalian lincRNAs, including relatively short lengths, fewer exons, lower expression levels, and stage-specific expression patterns. Functional studies motivated by "guilt-by-association" associated individual lincRNAs with specific GO functions, providing an important resource for future studies of lincRNA function. Most importantly, a weighted gene co-expression network analysis suggested that genes of the brown module were specifically associated with the day 2 stage. LincRNAs within this module were co-expressed with proteins involved in hematopoiesis and lipid metabolism. This study presents the systematic identification of lincRNAs in developing chicken embryos and will serve as a powerful resource for the study of lincRNA functions.

Genome-Wide Association Study and Transcriptome Differential Expression Analysis of the Feather Rate in Shouguang Chickens

The feather rate phenotype in chicks, including early-feathering and late-feathering phenotypes, are widely used as a sexing system in the poultry industry. The objective of this study was to obtain candidate genes associated with the feather rate in Shouguang chickens. In the present study, we collected 56 blood samples and 12 hair follicle samples of flight feathers from female Shouguang chickens. Then we identified the chromosome region associated with the feather rate by genome-wide association analysis (GWAS). We also performed RNA sequencing and analyzed differentially expressed genes between the early-feathering and late-feathering phenotypes using HISAT2, StringTie, and DESeq2. We identified a genomic region of 10.0–13.0 Mb of chromosome Z, which is statistically associated with the feather rate of Shouguang chickens at one-day old. After RNA sequencing analysis, 342 differentially expressed known genes between the early-feathering (EF) and late-feathering (LF) phenotypes were screened out, which were involved in epithelial cell differentiation, intermediate filament organization, protein serine kinase activity, peptidyl-serine phosphorylation, retinoic acid binding, and so on. The sperm flagellar 2 gene (SPEF2) and prolactin receptor (PRLR) gene were the only two overlapping genes between the results of GWAS and differential expression analysis, which implies that SPEF2 and PRLR are possible candidate genes for the formation of the chicken feathering phenotype in the present study. Our findings help to elucidate the molecular mechanism of the feather rate in chicks.

LncRNAs and their regulatory networks in breast muscle tissue of Chinese Gushi chickens during late postnatal development

BACKGROUND

Chicken skeletal muscle is an important economic product. The late stages of chicken development constitute the main period that affects meat production. LncRNAs play important roles in controlling the epigenetic process of growth and development. However, studies on the role of lncRNAs in the late stages of chicken breast muscle development are still lacking. In this study, to investigate the expression characteristics of lncRNAs during chicken muscle development, 12 cDNA libraries were constructed from Gushi chicken breast muscle samples from 6-, 14-, 22-, and 30-week-old chickens.

RESULTS

A total of 1252 new lncRNAs and 1376 annotated lncRNAs were identified. Furthermore, 53, 61, 50, 153, 117, and 78 DE-lncRNAs were found in the W14 vs. W6, W22 vs. W14, W22 vs. W6, W30 vs. W6, W30 vs. W14, and W30 vs. W22 comparison groups, respectively. After GO enrichment analysis of the DE-lncRNAs, several muscle development-related GO terms were found in the W22 vs. W14 comparison group. Moreover, it was found that the MAPK signaling pathway was one of the most frequently enriched pathways in the different comparison groups. In addition, 12 common target DE-miRNAs of DE-lncRNAs were found in different comparison groups, some of which were muscle-specific miRNAs, such as gga-miR-206, gga-miR-1a-3p, and miR-133a-3p. Interestingly, the precursors of four newly identified miRNAs were found to be homologous to lncRNAs. Additionally, we found some ceRNA networks associated with muscle development-related GO terms. For example, the ceRNA networks contained the DYNLL2 gene with 12 lncRNAs that targeted 2 miRNAs. We also constructed PPI networks, such as IGF-I-EGF and FZD6-WNT11.

CONCLUSIONS

This study revealed, for the first time, the dynamic changes in lncRNA expression in Gushi chicken breast muscle at different periods and revealed that the MAPK signaling pathway plays a vital role in muscle development. Furthermore, MEF2C and its target lncRNA may be involved in muscle regulation through the MAPK signaling pathway. This research provided valuable resources for elucidating posttranscriptional regulatory mechanisms to promote the development of chicken breast muscles after hatching.

miR-429-3p/LPIN1 Axis Promotes Chicken Abdominal Fat Deposition via PPARγ Pathway

To explore the regulatory mechanism of abdominal fat deposition in broilers, 100-day-old Sanhuang chickens (n = 12) were divided into high-fat and low-fat groups, according to the abdominal fat ratio size. Total RNA isolated from the 12 abdominal fat tissues was used for miRNA and mRNA sequencing. Results of miRNA and mRNA sequencing revealed that miR-429-3p was highly expressed in high-fat chicken whereas LPIN1 expression was downregulated. Further, we determined that miR-429-3p promoted preadipocyte proliferation and differentiation, whereas LPIN1 exerted an opposite effect. Notably, we found that the miR-429-3p/LPIN1 axis facilitated PPARγ pathway activation, which is closely associated with the progression of adipogenesis. In conclusion, our results provide evidence that a novel miR-429-3p/LPIN1 axis is involved in the regulation of adipogenesis, which may have a guiding role in the improvement of breeding for abdominal fat traits in broiler chickens.

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.

circFLT1 and lncCCPG1 Sponges miR-93 to Regulate the Proliferation and Differentiation of Adipocytes by Promoting lncSLC30A9 Expression

Although many circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs) have been discovered in adipocytes, their precise functions and molecular mechanisms remain poorly understood. Based on existing circRNA and lncRNA sequencing data of bovine adipocytes, we screened for the differential expression of circFLT1 and lncCCPG1 in preadipocytes and adipocytes and further analyzed their function and regulation during adipogenesis. The overexpression of circFLT1 and lncCCPG1 together facilitated adipocyte differentiation and suppressed proliferation. Computationally, the RNA hybrid showed that circFLT1 and lncCCPG1 had multiple potential binding sites with miR-93. Additionally, luciferase reporting experiments verified that circFLT1 and lncCCPG1 may interact with miR-93. We also demonstrated that overexpressed miR-93 effectively suppresses the expression of lncSLC30A9. Signaling pathway enrichment analysis, luciferase activity assay, and expression analysis revealed that lncSLC30A9 inhibits proliferation by inhibiting the expression of AKT protein and promotes differentiation by recruiting the FOS protein to the promoter of peroxisome proliferator-activated receptor gamma (PPARG). In sum, our results elucidate the regulatory mechanisms of circFLT1 and lncCCPG1 as miR-93 sponges in bovine adipocytes.

Molecular characterization and a duplicated 31-bp indel within the LDB2 gene and its associations with production performance in chickens

Many studies have shown that the LDB2 gene plays a regulatory role in retinal development and the cell cycle, but its biological role remains unclear. In this study, a 31-bp indel in the LDB2 gene was found for the first time on the basis of 2797 individuals from 10 different breeds, which led to different genotypes among individuals (II, ID and DD). Among these genotypes, DD was the most dominant. Association analysis of an F₂ resource population crossed with the Gushi (GS) chicken and Anka chicken showed that the DD genotype conferred a significantly greater semi-evisceration weight (SEW, 1108.665 g ± 6.263), evisceration weight (EW, 927.455 g ± 5.424), carcass weight (CW, 1197.306 g ± 6.443), breast muscle weight (BMW, 71.05 g ± 0.574), and leg muscle weight (LMW, 100.303 g ± 0.677) than the ID genotype (SEW, 1059.079 g ± 16.86; EW, 879.459 g ± 14.446; CW, 1141.821 g ± 17.176; BMW, 67.164 g ± 1.523; and LMW, 96.163 g ± 1.823). In addition, LDB2 gene expression in different breeds was significantly higher in the breast muscles and leg muscles than in other tissues. The expression level in the breast muscle differed significantly among stages of GS chicken development, with the highest expression observed at 6 weeks. The expression levels in the pectoral muscles differed significantly among Ross 308 genotypes. In summary, we studied the relationships between a 31-bp indel in the LDB2 gene and economic traits in chickens. The indel was significantly correlated with multiple growth and carcass traits in the F₂ resource population and affected the expression of the LDB2 gene in muscle tissue. In short, our study revealed that the LDB2 gene 31-bp indel can be used as a potential genetic marker for molecular breeding.