Copy Number Variation in SOX6 Contributes to Chicken Muscle Development

SOX6 AU controls myogenesis by cis-modulation of SOX6 in cattle

Long non-coding RNAs (lncRNAs) have been shown to be involved in the regulation of skeletal muscle development through multiple mechanisms. The present study revealed that the lncRNA SOX6 AU (SRY-box transcription factor 6 antisense upstream) is reverse transcribed from upstream of the bovine sex-determining region Y (SRY)-related high-mobility-group box 6 (SOX6) gene. SOX6 AU was significantly differentially expressed in muscle tissue among different developmental stages in Xianan cattle. Subsequently, knockdown and overexpression experiments discovered that SOX6 AU promoted primary skeletal muscle cells proliferation, apoptosis, and differentiation in bovine. The overexpression of SOX6 AU in bovine primary skeletal muscle cells resulted in 483 differentially expressed genes (DEGs), including 224 upregulated DEGs and 259 downregulated DEGs. GO functional annotation analysis showed that muscle development-related biological processes such as muscle structure development and muscle cell proliferation were significantly enriched. KEGG pathway analysis revealed that the PI3K/AKT and MAPK signaling pathways were important pathways for DEG enrichment. Notably, we found that SOX6 AU inhibited the mRNA and protein expression levels of the SOX6 gene. Moreover, knockdown of the SOX6 gene promoted the proliferation and apoptosis of bovine primary skeletal muscle cells. Finally, we showed that SOX6 AU promoted the proliferation and apoptosis of bovine primary skeletal muscle cells by cis-modulation of SOX6 in cattle. This work illustrates our discovery of the molecular mechanisms underlying the regulation of SOX6 AU in the development of beef.

Genome-wide association study of copy number variations with shank traits in a F2 crossbred chicken population

Copy number variations (CNVs) are large-scale changes in the DNA sequence that can affect the genetic structure and phenotype of an organism. The purpose of this study was to investigate the existing CNVs and their associations with the shank diameter (ShD) and shank length (ShL) traits using data from an F2 crossbred chicken population. To carry out the study, 312 chickens were genotyped using the Illumina 60k SNP Beadchip. The shank traits of the birds were measured from day 1 to 12 weeks of age. penncnv and cnvruler tools were used to find copy numbers and regions with copy number changes (CNVR), respectively. The CNVRanger package was used to perform a genome-wide association study between shank traits and CNVs. Gene ontology research in CNVRs was carried out using the david database. In this investigation, 966 CNVs and 606 regions with copy number changes were discovered. The copy number states and variations were randomly distributed along the length of the autosomal chromosomes. Weeks 1-4, 9 and 12 of growth revealed a significant association of copy number variations with shank traits, false discovery rate (FDR-corrected p-value < 0.01), and the majority of CNVs that were statistically significant were found on chromosomes 1-3. These CNV segments are nearby genes such as KCNJ12, FGF6 and MYF5, which are fundamental to growth and development. In addition, gene set analyses revealed terms related to muscle physiology, regulation of cellular processes and potassium channels.

MeRIP sequencing reveals the regulation of N6-methyladenosine in muscle development between hypertrophic and leaner broilers

Meat production performance is the most important economic trait in broilers, and skeletal muscle, as the largest organ in animals, is directly related to meat production during embryonic and postnatal growth and development. N6-Methyladenosine (m6A) is a chemical modification occurs on RNA adenosine that has been reported to participate in a variety of biological processes in all species. However, there are still few reports on the regulatory role of muscle growth and development in poultry after birth. This study aims to reveal the distribution of m6A modification sites in chicken pectoralis major muscle after birth and find out the regulatory relationship between m6A and muscle development. As representatives of leaner (Xinghua chicken [XH]) and hypertrophic (White Recessive Rock chicken [WRR]) broilers, there are significant differences in body weight, muscle fiber diameter, and muscle fiber cross-sectional area between XH and WRR chickens. RNA sequencing detected a total of 397 differentially expressed genes (DEG) in the pectoralis major muscle of XH and WRR chicken, and these DEGs were mainly enriched in catalytic activity and metabolic pathways. MeRIP sequencing results showed that among all 6,476 differentially modified m6A peaks, about 90% peaks (5,823) were differentially down regulated in XH chickens. The joint analysis of the mRNA and MeRIP sequencing data found 145 DEGs with differential m6A peak, ALKBH5 as a m6A demethylase, was also included. The highly expression of ALKBH5 in the muscle tissue of poultry and differential expression between XH and WRR chickens suggest that ALKBH5 may play a crucial role in regulating muscle development. Our results revealed that there were significant differences in growth rate, body weight, muscle fiber diameter, and fiber cross-section area between WRR and XH chicken, as well as significant differences in m6A methylation level and muscle metabolism level.

Unveiling the Influence of Copy Number Variations on Genetic Diversity and Adaptive Evolution in China's Native Pig Breeds via Whole-Genome Resequencing

Copy number variations (CNVs) critically influence individual genetic diversity and phenotypic traits. In this study, we employed whole-genome resequencing technology to conduct an in-depth analysis of 50 pigs from five local swine populations [Rongchang pig (RC), Wuzhishan pig (WZS), Tibetan pig (T), Yorkshire (YL) and Landrace (LR)], aiming to assess their genetic potential and explore their prospects in the field of animal model applications. We identified a total of 96,466 CNVs, which were subsequently integrated into 7112 non-redundant CNVRs, encompassing 1.3% of the swine genome. Functional enrichment analysis of the genes within these CNVRs revealed significant associations with sensory perception, energy metabolism, and neural-related pathways. Further selective scan analyses of the local pig breeds RC, T, WZS, along with YL and LR, uncovered that for the RC variety, the genes PLA2G10 and ABCA8 were found to be closely related to fat metabolism and cardiovascular health. In the T breed, the genes NCF2 and CSGALNACT1 were associated with immune response and connective tissue characteristics. As for the WZS breed, the genes PLIN4 and CPB2 were primarily linked to fat storage and anti-inflammatory responses. In summary, this research underscores the pivotal role of CNVs in fostering the diversity and adaptive evolution of pig breeds while also offering valuable insights for further exploration of the advantageous genetic traits inherent to China's local pig breeds. This facilitates the creation of experimental animal models tailored to the specific characteristics of these breeds, contributing to the advancement of livestock and biomedical research.

Genome-wide CNV analysis reveals candidate genes associated with milk production traits in water buffalo (Bubalus bubalis)

Buffaloes are vital contributors to the global dairy industry. Understanding the genetic basis of milk production traits in buffalo populations is essential for breeding programs and improving productivity. In this study, we conducted whole-genome resequencing on 387 buffalo genomes from 29 diverse Asian breeds, including 132 river buffaloes, 129 swamp buffaloes, and 126 crossbred buffaloes. We identified 36,548 copy number variant (CNVs) spanning 133.29 Mb of the buffalo genome, resulting in 2,100 copy number variant regions (CNVRs), with 1,993 shared CNVRs being found within the studied buffalo types. Analyzing CNVRs highlighted distinct genetic differentiation between river and swamp buffalo subspecies, verified by evolutionary tree and principal component analyses. Admixture analysis grouped buffaloes into river and swamp categories, with crossbred buffaloes displaying mixed ancestry. To identify candidate genes associated with milk production traits, we employed 3 approaches. First, we used Vst-based population differentiation, revealing 11 genes within CNVRs that exhibited significant divergence between different buffalo breeds, including genes linked to milk production traits. Second, expression quantitative loci (eQTL) analysis revealed differential expression of CNVR-driven genes (DECGs) associated with milk production traits. Notably, known milk production-related genes were among these DECGs, validating their relevance. Last, a genome-wide association study (GWAS) identified 3 CNVRs significantly linked to peak milk yield. Our study provides comprehensive genomic insights into buffalo populations and identifies candidate genes associated with milk production traits. These findings facilitate genetic breeding programs aimed at increasing milk yield and improving quality in this economically important livestock species.

Large-scale genomic and transcriptomic analyses elucidate the genetic basis of high meat yield in chickens

INTRODUCTION

Investigating the genetic markers and genomic signatures related to chicken meat production by combing multi-omics methods could provide new insights into modern chicken breeding technology systems.

OBJECT

Chicken is one of the most efficient and environmentally friendly livestock, especially the fast-growing white-feathered chicken (broiler), which is well known for high meat yield, but the underlying genetic basis is poorly understood.

METHOD

We generated whole-genome resequencing of three purebred broilers (n = 748) and six local breeds/lines (n = 114), and sequencing data of twelve chicken breeds (n = 199) were obtained from the NCBI database. Additionally, transcriptome sequencing of six tissues from two chicken breeds (n = 129) at two developmental stages was performed. A genome-wide association study combined with cis-eQTL mapping and the Mendelian randomization was applied.

RESULT

We identified > 17 million high-quality SNPs, of which 21.74% were newly identified, based on 21 chicken breeds/lines. A total of 163 protein-coding genes underwent positive selection in purebred broilers, and 83 genes were differentially expressed between purebred broilers and local chickens. Notably, muscle development was proven to be the major difference between purebred broilers and local chickens, or ancestors, based on genomic and transcriptomic evidence from multiple tissues and stages. The MYH1 gene family showed the top selection signatures and muscle-specific expression in purebred broilers. Furthermore, we found that the causal gene SOX6 influenced breast muscle yield and also related to myopathy occurrences. A refined haplotype was provided, which had a significant effect on SOX6 expression and phenotypic changes.

CONCLUSION

Our study provides a comprehensive atlas comprising the typical genomic variants and transcriptional characteristics for muscle development and suggests a new regulatory target (SOX6-MYH1s axis) for breast muscle yield and myopathy, which could aid in the development of genome-scale selective breeding aimed at high meat yield in broiler chickens.

Sox6 Differentially Regulates Inherited Myogenic Abilities and Muscle Fiber Types of Satellite Cells Derived from Fast- and Slow-Type Muscles

Adult skeletal muscle is primarily divided into fast and slow-type muscles, which have distinct capacities for regeneration, metabolism and contractibility. Satellite cells plays an important role in adult skeletal muscle. However, the underlying mechanisms of satellite cell myogenesis are poorly understood. We previously found that Sox6 was highly expressed in adult fast-type muscle. Therefore, we aimed to validate the satellite cell myogenesis from different muscle fiber types and investigate the regulation of Sox6 on satellite cell myogenesis. First, we isolated satellite cells from fast- and slow-type muscles individually. We found that satellite cells derived from different muscle fiber types generated myotubes similar to their origin types. Further, we observed that cells derived from fast muscles had a higher efficiency to proliferate but lower potential to self-renew compared to the cells derived from slow muscles. Then we demonstrated that Sox6 facilitated the development of satellite cells-derived myotubes toward their inherent muscle fiber types. We revealed that higher expression of Nfix during the differentiation of fast-type muscle-derived myogenic cells inhibited the transcription of slow-type isoforms (MyH7B, Tnnc1) by binding to Sox6. On the other hand, Sox6 activated Mef2C to promote the slow fiber formation in slow-type muscle-derived myogenic cells with Nfix low expression, showing a different effect of Sox6 on the regulation of satellite cell development. Our findings demonstrated that satellite cells, the myogenic progenitor cells, tend to develop towards the fiber type similar to where they originated. The expression of Sox6 and Nfix partially explain the developmental differences of myogenic cells derived from fast- and slow-type muscles.

Transcriptome sequencing analysis of the role of miR-499-5p and SOX6 in chicken skeletal myofiber specification

MicroRNAs (miRNAs) might play critical roles in skeletal myofiber specification. In a previous study, we found that chicken miR-499-5p is specifically expressed in slow-twitch muscle and that its potential target gene is SOX6. In this study, we performed RNA sequencing to investigate the effects of SOX6 and miR-499-5p on the modulation and regulation of chicken muscle fiber type and its regulatory mechanism. The expression levels of miR-499-5p and SOX6 demonstrated opposing trends in different skeletal muscles and were associated with muscle fiber type composition. Differential expression analysis revealed that miR-499-5p overexpression led to significant changes in the expression of 297 genes in chicken primary myoblasts (CPMs). Myofiber type-related genes, including MYH7B and CSRP3, showed expression patterns similar to those in slow-twitch muscle. According to functional enrichment analysis, differentially expressed genes were mostly associated with muscle development and muscle fiber-related processes. SOX6 was identified as the target gene of miR-499-5p in CPM using target gene mining and luciferase reporter assays. SOX6 knockdown resulted in upregulation of the slow myosin genes and downregulation of fast myosin genes. Furthermore, protein-protein interaction network analysis revealed that MYH7B and RUNX2 may be the direct targets of SOX6. These results indicated that chicken miR-499-5p may promote slow-twitch muscle fiber formation by repressing SOX6 expression. Our study provides a dataset that can be used as a reference for animal meat quality and human muscle disease studies.

Genome-wide detection of copy number variation in American mink using whole-genome sequencing

BACKGROUND

Copy number variations (CNVs) represent a major source of genetic diversity and contribute to the phenotypic variation of economically important traits in livestock species. In this study, we report the first genome-wide CNV analysis of American mink using whole-genome sequence data from 100 individuals. The analyses were performed by three complementary software programs including CNVpytor, DELLY and Manta.

RESULTS

A total of 164,733 CNVs (144,517 deletions and 20,216 duplications) were identified representing 5378 CNV regions (CNVR) after merging overlapping CNVs, covering 47.3 Mb (1.9%) of the mink autosomal genome. Gene Ontology and KEGG pathway enrichment analyses of 1391 genes that overlapped CNVR revealed potential role of CNVs in a wide range of biological, molecular and cellular functions, e.g., pathways related to growth (regulation of actin cytoskeleton, and cAMP signaling pathways), behavior (axon guidance, circadian entrainment, and glutamatergic synapse), lipid metabolism (phospholipid binding, sphingolipid metabolism and regulation of lipolysis in adipocytes), and immune response (Wnt signaling, Fc receptor signaling, and GTPase regulator activity pathways). Furthermore, several CNVR-harbored genes associated with fur characteristics and development (MYO5A, RAB27B, FGF12, SLC7A11, EXOC2), and immune system processes (SWAP70, FYN, ORAI1, TRPM2, and FOXO3).

CONCLUSIONS

This study presents the first genome-wide CNV map of American mink. We identified 5378 CNVR in the mink genome and investigated genes that overlapped with CNVR. The results suggest potential links with mink behaviour as well as their possible impact on fur quality and immune response. Overall, the results provide new resources for mink genome analysis, serving as a guideline for future investigations in which genomic structural variations are present.

Assessment of linkage disequilibrium patterns between structural variants and single nucleotide polymorphisms in three commercial chicken populations

BACKGROUND

Structural variants (SV) are causative for some prominent phenotypic traits of livestock as different comb types in chickens or color patterns in pigs. Their effects on production traits are also increasingly studied. Nevertheless, accurately calling SV remains challenging. It is therefore of interest, whether close-by single nucleotide polymorphisms (SNPs) are in strong linkage disequilibrium (LD) with SVs and can serve as markers. Literature comes to different conclusions on whether SVs are in LD to SNPs on the same level as SNPs to other SNPs. The present study aimed to generate a precise SV callset from whole-genome short-read sequencing (WGS) data for three commercial chicken populations and to evaluate LD patterns between the called SVs and surrounding SNPs. It is thereby the first study that assessed LD between SVs and SNPs in chickens.

RESULTS

The final callset consisted of 12,294,329 bivariate SNPs, 4,301 deletions (DEL), 224 duplications (DUP), 218 inversions (INV) and 117 translocation breakpoints (BND). While average LD between DELs and SNPs was at the same level as between SNPs and SNPs, LD between other SVs and SNPs was strongly reduced (DUP: 40%, INV: 27%, BND: 19% of between-SNP LD). A main factor for the reduced LD was the presence of local minor allele frequency differences, which accounted for 50% of the difference between SNP - SNP and DUP - SNP LD. This was potentially accompanied by lower genotyping accuracies for DUP, INV and BND compared with SNPs and DELs. An evaluation of the presence of tag SNPs (SNP in highest LD to the variant of interest) further revealed DELs to be slightly less tagged by WGS SNPs than WGS SNPs by other SNPs. This difference, however, was no longer present when reducing the pool of potential tag SNPs to SNPs located on four different chicken genotyping arrays.

CONCLUSIONS

The results implied that genomic variance due to DELs in the chicken populations studied can be captured by different SNP marker sets as good as variance from WGS SNPs, whereas separate SV calling might be advisable for DUP, INV, and BND effects.

Identification and Tissue-Specific Characterization of Novel SHOX-Regulated Genes in Zebrafish Highlights SOX Family Members Among Other Genes

SHOX deficiency causes a spectrum of clinical phenotypes related to skeletal dysplasia and short stature, including Léri-Weill dyschondrosteosis, Langer mesomelic dysplasia, Turner syndrome, and idiopathic short stature. SHOX controls chondrocyte proliferation and differentiation, bone maturation, and cellular growth arrest and apoptosis via transcriptional regulation of its direct target genes NPPB, FGFR3, and CTGF. However, our understanding of SHOX-related pathways is still incomplete. To elucidate the underlying molecular mechanisms and to better understand the broad phenotypic spectrum of SHOX deficiency, we aimed to identify novel SHOX targets. We analyzed differentially expressed genes in SHOX-overexpressing human fibroblasts (NHDF), and confirmed the known SHOX target genes NPPB and FGFR among the most strongly regulated genes, together with 143 novel candidates. Altogether, 23 genes were selected for further validation, first by whole-body characterization in developing shox-deficient zebrafish embryos, followed by tissue-specific expression analysis in three shox-expressing zebrafish tissues: head (including brain, pharyngeal arches, eye, and olfactory epithelium), heart, and pectoral fins. Most genes were physiologically relevant in the pectoral fins, while only few genes were also significantly regulated in head and heart tissue. Interestingly, multiple sox family members (sox5, sox6, sox8, and sox18) were significantly dysregulated in shox-deficient pectoral fins together with other genes (nppa, nppc, cdkn1a, cdkn1ca, cyp26b1, and cy26c1), highlighting an important role for these genes in shox-related growth disorders. Network-based analysis integrating data from the Ingenuity pathways revealed that most of these genes act in a common network. Our results provide novel insights into the genetic pathways and molecular events leading to the clinical manifestation of SHOX deficiency.

Genome-wide detection of CNVs and their association with performance traits in broilers

Background

Copy number variations (CNVs) are a major type of structural genomic variants that underlie genetic architecture and phenotypic variation of complex traits, not only in humans, but also in livestock animals. We identified CNVs along the chicken genome and analyzed their association with performance traits. Genome-wide CNVs were inferred from Affymetrix® high density SNP-chip data for a broiler population. CNVs were concatenated into segments and association analyses were performed with linear mixed models considering a genomic relationship matrix, for birth weight, body weight at 21, 35, 41 and 42 days, feed intake from 35 to 41 days, feed conversion ratio from 35 to 41 days and, body weight gain from 35 to 41 days of age.

Results

We identified 23,214 autosomal CNVs, merged into 5042 distinct CNV regions (CNVRs), covering 12.84% of the chicken autosomal genome. One significant CNV segment was associated with BWG on GGA3 (q-value = 0.00443); one significant CNV segment was associated with BW35 (q-value = 0.00571), BW41 (q-value = 0.00180) and BW42 (q-value = 0.00130) on GGA3, and one significant CNV segment was associated with BW on GGA5 (q-value = 0.00432). All significant CNV segments were verified by qPCR, and a validation rate of 92.59% was observed. These CNV segments are located nearby genes, such as KCNJ11, MyoD1 and SOX6, known to underlie growth and development. Moreover, gene-set analyses revealed terms linked with muscle physiology, cellular processes regulation and potassium channels.

Conclusions

Overall, this CNV-based GWAS study unravels potential candidate genes that may regulate performance traits in chickens. Our findings provide a foundation for future functional studies on the role of specific genes in regulating performance in chickens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07676-1.

Genome-wide comparative analyses reveal selection signatures underlying adaptation and production in Tibetan and Poll Dorset sheep

The identification of genome-wide selection signatures can provide insights on the mechanisms of natural and/or artificial selection and uncover genes related to biological functions and/or phenotypes. Tibetan sheep are an important livestock in Tibet, providing meat and wool for Tibetans who are renown for breeding livestock that adapt well to high altitudes. Using whole-genome sequences with an effective sequencing depth of 5×, we investigated the genomic diversity and structure and, identified selection signatures of White Tibetan, Oula and Poll Dorset sheep. We obtained 30,163,679 Single Nucleotide Polymorphisms (SNPs) and 5,388,372 indels benchmarked against the ovine Oar_v4.0 genome assembly. Next, using F_ST, ZHp and XP-EHH approaches, we identified selection signatures spanning a set of candidate genes, including HIF1A, CAPN3, PRKAA1, RXFP2, TRHR and HOXA10 that are associated with pathways and GO categories putatively related to hypoxia responses, meat traits and disease resistance. Candidate genes and GO terms associated with coat color were also identified. Finally, quantification of blood physiological parameters, revealed higher levels of mean corpuscular hemoglobin measurement and mean corpuscular hemoglobin concentration in Tibetan sheep compared with Poll Dorset, suggesting a greater oxygen-carrying capacity in the Tibetan sheep and thus better adaptation to high-altitude hypoxia. In conclusion, this study provides a greater understanding of genome diversity and variations associated with adaptive and production traits in sheep.

Polymorphisms of AMY1A gene and their association with growth, carcass traits and feed intake efficiency in chickens

Investigations on the association between chicken traits and genetic variations can provide basic information to improve production performance in chickens. In our previous work, we genotyped 450 male chickens with a 600 K SNP array [1] and found that several SNPs in the genomic regions of the amylase alpha 1A (AMY1A) gene were significantly associated with feed intake efficiency and carcass traits. Given the lower accuracy of the SNP array, we performed direct sequencing with male and female chickens to further test chicken AMY1A polymorphisms and investigate their association with 17 traits in chickens. The results showed that 7 SNPs in the 5' flanking region, exon, intron and 3' UTR (3' untranslated region) of AMY1A, were significantly associated with daily gain (DG), average daily feed intake (ADFI), leg muscle weight (LMW) and abdominal fat (AF) (p < 0.05). Additionally, the haplotypes based on three SNPs, rs15910189, rs314354067 and rs316026696, showed significant associations with DG (p < 0.01), ADFI and AF (p < 0.05). To better understand the transcriptional regulation of AMY1A, we cloned its 5' flanking region and found that the SNPs rs316436216 and rs314213090 which might change the transcriptional regulator binding sites, were in the suppressor and enhancer regions, respectively. In addition, luciferase assays revealed that the SNP rs314613110 in the 3' UTR influenced the binding of the miRNA gga-miR-1764-3p. To validate whether there is any copy number variation in AMY1A in our population, we performed a genome-wide assessment of CNVs through whole-genome resequencing data. However, no CNV was found in AMY1A in our population, which is different from the increased copy number of AMY1A found in humans who consume a high-starch diet. Therefore, the present study provides substantial evidence for the association of AMY1A polymorphisms with growth traits and feed intake efficiency, which might contribute to chicken breeding programs.

Emerging roles of Sox6 in the renal and cardiovascular system

The function of Sex-determining Region Y (SRY)-related high-mobility-group box (Sox) family of transcription factors in cell fate decisions during embryonic development are well-established. Accumulating evidence indicates that the Sox family of transcription factors are fundamental in adult tissue homeostasis, regeneration, and physiology. The SoxD subfamily of genes are expressed in various cell types of different organs during embryogenesis and adulthood and have been involved in cell-fate determination, cellular proliferation and survival, differentiation, and terminal maturation in a number of cell lineages. The dysregulation in the function of SoxD proteins (i.e. Sox5, Sox6, Sox13, and Sox23) have been implicated in different disease conditions such as chondrodysplasia, cancer, diabetes, hypertension, autoimmune diseases, osteoarthritis among others. In this minireview, we present recent developments related to the transcription factor Sox6, which is involved in a number of diseases such as diabetic nephropathy, adipogenesis, cardiomyopathy, inflammatory bowel disease, and cancer. Sox6 has been implicated in the regulation of renin expression and JG cell recruitment in mice during sodium depletion and dehydration. We provide a current perspective of Sox6 research developments in last five years, and the implications of Sox6 functions in cardiovascular physiology and disease conditions.

Genomic variations and signatures of selection in Wuhua yellow chicken

Wuhua yellow chicken (WHYC) is an important traditional yellow-feathered chicken from China, which is characterized by its white tail feathers, white flight feathers, and strong disease resistance. However, the genomic basis of these unique traits associated with WHYC is poorly understood. In this study, whole-genome resequencing was performed with an average coverage of 20.77-fold to investigate heritable variation and identify selection signals in WHYC. Reads were mapped onto the chicken reference genome (Galgal5) with a coverage of 85.95%. After quality control, 11,953,471 single nucleotide polymorphisms and 1,069,574 insertion/deletions were obtained. In addition, 41,408 structural variants and 33,278 copy number variants were found. Comparative genomic analysis of WHYC and other yellow-feathered chicken breeds showed that selected regions were enriched in genes involved in transport and catabolism, immune system, infectious diseases, signal transduction, and signaling molecules and interactions. Several genes associated with disease resistance were also identified, including IFNA, IFNB, CD86, IL18, IL11RA, VEGFC, and ATG10. Furthermore, our results suggest that PMEL and TYRP1 may contribute to the white feather coloring in WHYC. These findings can improve our understanding of the genetic characteristics of WHYC and may contribute to future breed improvement.

Genome-wide genetic structure and selection signatures for color in 10 traditional Chinese yellow-feathered chicken breeds

BACKGROUND

Yellow-feathered chickens (YFCs) have a long history in China. They are well-known for the nutritional and commercial importance attributable to their yellow color phenotype. Currently, there is a huge paucity in knowledge of the genetic determinants responsible for phenotypic and biochemical properties of these iconic chickens. This study aimed to uncover the genetic structure and the molecular underpinnings of the YFCs trademark coloration.

RESULTS

The whole-genomes of 100 YFCs from 10 major traditional breeds and 10 Huaibei partridge chickens from China were re-sequenced. Comparative population genomics based on autosomal single nucleotide polymorphisms (SNPs) revealed three geographically based clusters among the YFCs. Compared to other Chinese indigenous chicken genomes incorporated from previous studies, a closer genetic proximity within YFC breeds than between YFC breeds and other chicken populations is evident. Through genome-wide scans for selective sweeps, we identified RALY heterogeneous nuclear ribonucleoprotein (RALY), leucine rich repeat containing G protein-coupled receptor 4 (LGR4), solute carrier family 23 member 2 (SLC23A2), and solute carrier family 2 member 14 (SLC2A14), besides the classical beta-carotene dioxygenase 2 (BCDO2), as major candidates pigment determining genes in the YFCs.

CONCLUSION

We provide the first comprehensive genomic data of the YFCs. Our analyses show phylogeographical patterns among the YFCs and potential candidate genes giving rise to the yellow color trait of the YFCs. This study lays the foundation for further research on the genome-phenotype cross-talks that define important poultry traits and for formulating genetic breeding and conservation strategies for the YFCs.

Associations of ORMDL1 gene copy number variations with growth traits in four Chinese sheep breeds

Copy number variations (CNVs) are gains and losses of genomic sequence of more than 50 bp between two individuals of a species. Also, CNV is considered to be one of the main elements affecting the phenotypic diversity and evolutionary adaptation of animals. ORMDL sphingolipid biosynthesis regulator 1 (ORMDL1) is a protein-coding gene associated with diseases and development. In our study, the polymorphism of ORMDL1 gene copy numbers in four Chinese sheep breeds (abbreviated CK, HU, STH, and LTH) was detected. In addition, we analyzed the transcriptional expression level of ORMDL1 gene in different tissues of sheep and examined the association of ORMDL1 CNV with growth traits. The statistical analysis revealed that ORMDL1 CNV was remarkably correlated with body height, heart girth, and circumference of cannon bone in HU sheep (P<0.05), and there are significant effects on body weight, body height, body length, chest depth, and height of hip cross in STH sheep (P<0.05). In conclusion, our results provide a basis for the relationship between CNV of ORMDL1 gene and sheep growth traits, suggesting that ORMDL1 CNV may be considered a promising marker for the molecular breeding of Chinese sheep.

Identification of Copy Number Variation in Domestic Chicken Using Whole-Genome Sequencing Reveals Evidence of Selection in the Genome

Simple Summary

Chickens have been bred for meat and egg production as a source of animal protein. With the increase of productivity as the main purpose of domestication, factors such as metabolism and immunity were boosted, which are detectable signs of selection on the genome. This study focused on copy number variation (CNV) to find evidence of domestication on the genome. CNV was detected from whole-genome sequencing of 65 chickens including Red Jungle Fowl, broilers, and layers. After that, CNV region, the overlapping region of CNV between individuals, was made to identify which genomic regions showed copy number differentiation. The 663 domesticated-specific CNV regions were associated with various functions such as metabolism and organ development. Also, by performing population differentiation analyses such as clustering analysis and ANOVA test, we found that there are a lot of genomic regions with different copy number patterns between broilers and layers. This result indicates that different genetic variations can be found, depending on the purpose of artificial selection and provides considerations for future animal breeding.

Abstract

Copy number variation (CNV) has great significance both functionally and evolutionally. Various CNV studies are in progress to find the cause of human disease and to understand the population structure of livestock. Recent advances in next-generation sequencing (NGS) technology have made CNV detection more reliable and accurate at whole-genome level. However, there is a lack of CNV studies on chickens using NGS. Therefore, we obtained whole-genome sequencing data of 65 chickens including Red Jungle Fowl, Cornish (broiler), Rhode Island Red (hybrid), and White Leghorn (layer) from the public databases for CNV region (CNVR) detection. Using CNVnator, a read-depth based software, a total of 663 domesticated-specific CNVRs were identified across autosomes. Gene ontology analysis of genes annotated in CNVRs showed that mainly enriched terms involved in organ development, metabolism, and immune regulation. Population analysis revealed that CN and RIR are closer to each other than WL, and many genes (LOC772271, OR52R1, RD3, ADH6, TLR2B, PRSS2, TPK1, POPDC3, etc.) with different copy numbers between breeds found. In conclusion, this study has helped to understand the genetic characteristics of domestic chickens at CNV level, which may provide useful information for the development of breeding systems in chickens.

Genome-wide detection of copy number variations in polled yak using the Illumina BovineHD BeadChip

Background

Copy number variations (CNVs), which are genetic variations present throughout mammalian genomes, are a vital source of phenotypic variation that can lead to the development of unique traits. In this study we used the Illunima BovineHD BeadChip to conduct genome-wide detection of CNVs in 215 polled yaks.

Results

A total of 1066 CNV regions (CNVRs) were detected with a total length of 181.6 Mb, comprising ~ 7.2% of the bovine autosomal genome. The size of these CNVRs ranged from 5.53 kb to 1148.45 kb, with an average size of 170.31 kb. Eight out of nine randomly chosen CNVRs were successfully validated by qPCR. A functional enrichment analysis of the CNVR-associated genes indicated their relationship to a number of molecular adaptations that enable yaks to thrive at high altitudes. One third of the detected CNVRs were mapped to QTLs associated with six classes of economically important traits, indicating that these CNVRs may play an important role in variations of these traits.

Conclusions

Our genome-wide yak CNV map may thus provide valuable insights into both the molecular mechanisms of high altitude adaptation and the potential genomic basis of economically important traits in yak.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5759-1) contains supplementary material, which is available to authorized users.

Sox Genes Show Spatiotemporal Expression during Murine Tongue and Eyelid Development

The tongue is a critical organ, involved in functions such as speaking, swallowing, mastication, and degustation. Although Sox genes are known to play critical roles in many biological processes, including organogenesis, the expression of the Sox family members during tongue development remains unclear. We therefore performed a comparative in situ hybridization analysis of 17 Sox genes (Sox1–14, 17, 18, and 21) during murine tongue development. Sox2, 4, 6, 8, 9, 10, 11, 12, and 21 were found to be expressed in the tongue epithelium, whereas Sox2, 4–6, 8–11, 13, and 21 showed expression in the mesenchyme of the developing tongue. Expression of Sox1, 4, 6, 8–12, and 21 were observed in the developing tongue muscle. Sox5 and 13 showed expression only at E12, while Sox1 expression was observed only on E18. Sox6, 8, 9, and 12 showed expression at several stages. Although the expression of Sox2, 4, 10, 11, and 21 was detected during all the four stages of tongue development, their expression patterns differed among the stages. We thus identified a dynamic spatiotemporal expression pattern of the Sox genes during murine tongue development. To understand whether Sox genes are involved in the development of other craniofacial organs through similar roles to those in tongue development, we also examined the expression of Sox genes in eyelid primordia, which also contain epithelium, mesenchyme, and muscle. However, expression patterns and timing of Sox genes differed between tongue and eyelid development. Sox genes are thus related to organogenesis through different functions in each craniofacial organ.

Copy number variation in livestock: A mini review

Copy number variation (CNV) is a phenomenon in which sections of the genome, ranging from one kilo base pair (Kb) to several million base pairs (Mb), are repeated and the number of repeats vary between the individuals in a population. It is an important source of genetic variation in an individual which is now being utilized rather than single nucleotide polymorphisms (SNPs), as it covers the more genomic region. CNVs alter the gene expression and change the phenotype of an individual due to deletion and duplication of genes in the copy number variation regions (CNVRs). Earlier, researchers extensively utilized SNPs as the main source of genetic variation. But now, the focus is on identification of CNVs associated with complex traits. With the recent advances and reduction in the cost of sequencing, arrays are developed for genotyping which cover the maximum number of SNPs at a time that can be used for detection of CNVRs and underlying quantitative trait loci (QTL) for the complex traits to accelerate genetic improvement. CNV studies are also being carried out to understand the evolutionary mechanism in the domestication of livestock and their adaptation to the different environmental conditions. The main aim of the study is to review the available data on CNV and its role in genetic variation among the livestock.