Associations between variants of bone morphogenetic protein 7 gene and growth traits in chickens

Artificial selection footprints in indigenous and commercial chicken genomes

BACKGROUND

Although many studies have been done to reveal artificial selection signatures in commercial and indigenous chickens, a limited number of genes have been linked to specific traits. To identify more trait-related artificial selection signatures and genes, we re-sequenced a total of 85 individuals of five indigenous chicken breeds with distinct traits from Yunnan Province, China.

RESULTS

We found 30 million non-redundant single nucleotide variants and small indels (< 50 bp) in the indigenous chickens, of which 10 million were not seen in 60 broilers, 56 layers and 35 red jungle fowls (RJFs) that we compared with. The variants in each breed are enriched in non-coding regions, while those in coding regions are largely tolerant, suggesting that most variants might affect cis-regulatory sequences. Based on 27 million bi-allelic single nucleotide polymorphisms identified in the chickens, we found numerous selective sweeps and affected genes in each indigenous chicken breed and substantially larger numbers of selective sweeps and affected genes in the broilers and layers than previously reported using a rigorous statistical model. Consistent with the locations of the variants, the vast majority (~ 98.3%) of the identified selective sweeps overlap known quantitative trait loci (QTLs). Meanwhile, 74.2% known QTLs overlap our identified selective sweeps. We confirmed most of previously identified trait-related genes and identified many novel ones, some of which might be related to body size and high egg production traits. Using RT-qPCR, we validated differential expression of eight genes (GHR, GHRHR, IGF2BP1, OVALX, ELF2, MGARP, NOCT, SLC25A15) that might be related to body size and high egg production traits in relevant tissues of relevant breeds.

CONCLUSION

We identify 30 million single nucleotide variants and small indels in the five indigenous chicken breeds, 10 million of which are novel. We predict substantially more selective sweeps and affected genes than previously reported in both indigenous and commercial breeds. These variants and affected genes are good candidates for further experimental investigations of genotype-phenotype relationships and practical applications in chicken breeding programs.

The 12-bp indel in the SMAD family member 2 gene is associated with goat growth traits

SMAD family member 2 (SMAD2) is a member of the TGFβ signaling pathway and functions as an essential regulator in the processes of development, cell proliferation, and bone formation. A previous observation reported that a 12-bp deletion of this gene affected the litter size in goats. However, according to our knowledge, no study has reported an association between this polymorphism and goat body measurement traits. The purpose of this study was to investigate the association of the insertion/deletion (indel) within the SMAD2 gene with the growth traits of goats. The indel polymorphism was found to be significantly associated with chest width and bust (p < 0.05), while cannon circumference was significantly the strongest compared to other traits (p < 0.01) and individuals with the DD genotypes were more dominant genotypes than other genotypes. In summary, we found evidence that the 12-bp indel within the SMAD2 gene could improve goat body measurement traits, paving the way for marker-assisted selection in the field of goat genetics and breeding.

Molecular Characterization and Function of Bone Morphogenetic Protein 7 (BMP7) in the Pacific Abalone, Haliotis discus hannai

Bone morphogenetic proteins (BMPs) play important roles in a lot of biological processes, such as bone development, cell proliferation, cell differentiation, growth, etc. However, the functions of abalone BMP genes are still unknown. This study aimed to better understand the characterization and biological function of BMP7 of Haliotis discus hannai (hdh-BMP7) via cloning and sequencing analysis. The coding sequence (CDS) length of hdh-BMP7 is 1251 bp, which encodes 416 amino acids including a signal peptide (1-28 aa), a transforming growth factor-β (TGF-β) propeptide (38-272 aa), and a mature TGF-β peptide (314-416 aa). The analysis of expression showed that hdh-BMP7 mRNA was widely expressed in all the examined tissues of H. discus hannai. Four SNPs were related to growth traits. The results of RNA interference (RNAi) showed that the mRNA expression levels of hdh-BMPR I, hdh-BMPR II, hdh-smad1, and hdh-MHC declined after hdh-BMP7 was silenced. After RNAi experiment for 30 days, the shell length, shell width, and total weight were found to be reduced in H. discus hannai (p < 0.05). The results of real-time quantitative reverse transcription PCR revealed that the hdh-BMP7 mRNA was lower in abalone of the S-DD-group than in the L-DD-group. Based on these data, we hypothesized that BMP7 gene has a positive role in the growth of H. discus hannai.

Functional Differentiation of BMP7 Genes in Zebrafish: bmp7a for Dorsal-Ventral Pattern and bmp7b for Melanin Synthesis and Eye Development

Bone morphogenetic protein 7 (BMP7) belongs to the transforming growth factor β (TGF-β) family, which not only induces cartilage and bone formation, but also regulates eye development and melanoma tumorigenesis in mammals. In teleosts, BMP7 differentiates into two subtypes, bmp7a and bmp7b, which have clearly differentiated structures. To fully understand the functional differentiation of bmp7a and bmp7b in fish species, we successfully constructed bmp7a and bmp7b gene deletion mutants in zebrafish using CRISPR/Cas9-mediated gene editing technology. Our results showed that bmp7a mutation caused abnormal development of the embryo’s dorsal-ventral pattern that led to death; bmp7b mutation induced growth inhibition and increased melanin production in the skin and eye of mutants. Histological analysis revealed that melanin in the retina of the eyes in bmp7b mutants increased, and behavioral observation showed that the vision and sensitivity to food of the mutants were reduced. Transcriptome analysis of the skin and eye tissues showed that the expression changes of wnt7ba and gna14 in bmp7b mutants might promote the increase of melanin. Additionally, the eye transcriptome analysis indicated that changes in the structure of the eyes in bmp7b mutants led to defects in phototransduction, and seven DEGs (rgs9a, rgs9b, rcvrn2, guca1d, grk1b, opn1mw4, and gc2) were identified as key candidate genes that affected the photonic response of the eyes. The study revealed the functional differentiation of bmp7a and bmp7b in teleosts and the first report about the inhibitory effect of bmp7b on melanogenesis may provide useful information for the future research on human melanoma-related diseases.

Associations of CALM1 and DRD1 polymorphisms, and their expression levels, with Taihang chicken egg-production traits

Egg production by hens is an important reproductive performance index in the poultry industry. To investigate the effects of the CALM1 and DRD1 genes on egg production in chicken, their mRNA expression and single nucleotide polymorphisms (SNP) levels were investigated, and bioinformatics and egg-production association analyses were performed. Three SNPs (g.44069941G > A and g.44069889A > G in CALM1 and g.10742639C > T in DRD1) were detected in the exons and introns of CALM1 and DRD1 in 400 Taihang chickens. Among them, g.44069941G > A was significantly associated with Taihang chicken egg production on the 500th day (p < 0.05), whereas g.10742639C > T was significantly associated with the 300th day (p < 0.05). The expression levels of CALM1 and DRD1 in ovarian tissues of a high-yielding Taihang group were greater than in a low-yielding group (p < 0.05). The bioinformatics analysis revealed that the mutations influenced the mRNA secondary structures of CALM1 and DRD1. This study provides new insights into the potential effects of CALM1 and DRD1 polymorphisms on chicken egg production. The two SNPs g.44069941G > A and g.10742639C > T are potential molecular markers for improving the reproductive traits of Taihang chicken.

A Potential Role of Bone Morphogenetic Protein 7 in Shell Formation and Growth in the Razor Clam Sinonovacula constricta

Bone morphogenetic proteins (BMPs) not only play essential roles in bone development but also are involved in embryonic growth, organogenesis cell proliferation and differentiation. However, the previous studies on the functions of shellfish BMPs genes are still very limited. To better understand its molecular structure and biological function, BMP7 of the razor clam Sinonovacula constricta (Sc-BMP7) was cloned and characterized in this study. The full length of Sc-BMP7 is 2252 bp, including an open reading frame (ORF) of 1257 bp encoding 418 amino acids. The protein sequence included a signal peptide (1–32 aa), a prodomain (38–270 aa) and a TGF-β domain (317–418 aa). The quantitative expression of eleven adult tissues showed that Sc-BMP7 was significantly higher expressed in the gill, foot, and mantle (P < 0.05), but lower in hemocytes and hepatopancreas. In the early development stages, low expression was detected in the stages of unfertilized mature eggs, fertilized eggs, 4-cell embryos, blastula, gastrulae, whereas it increased after the stage of trochophore and demonstrated the highest expression in umbo larvae (P < 0.01). In shell repair experiment, Sc-BMP7 showed increasing expression level after 12 h. The higher expression of Sc-BMP7 was detected while Ca²⁺ concentration was reduced in seawater. After inhibiting Sc-BMP7 expression using RNA interference (RNAi) technology, expression of Sc-BMP7 mRNA and protein were significantly down-regulated (P < 0.05) in the central zone of mantle (nacre formation related tissue) and the pallial zone of mantle (prismatic layer formation related tissue). Association analysis identified two shared SNPs in exon of Sc-BMP7 gene from 246 individuals of two groups. These results indicated that BMP7 might be involved in shell formation and growth. These results would contribute to clarify the role of Sc-BMP7 in the regulation of growth and shell formation, and provide growth-related markers for molecular marker assisted breeding of this species.

Combining Random Forests and a Signal Detection Method Leads to the Robust Detection of Genotype-Phenotype Associations

Genome wide association studies (GWAS) are a well established methodology to identify genomic variants and genes that are responsible for traits of interest in all branches of the life sciences. Despite the long time this methodology has had to mature the reliable detection of genotype-phenotype associations is still a challenge for many quantitative traits mainly because of the large number of genomic loci with weak individual effects on the trait under investigation. Thus, it can be hypothesized that many genomic variants that have a small, however real, effect remain unnoticed in many GWAS approaches. Here, we propose a two-step procedure to address this problem. In a first step, cubic splines are fitted to the test statistic values and genomic regions with spline-peaks that are higher than expected by chance are considered as quantitative trait loci (QTL). Then the SNPs in these QTLs are prioritized with respect to the strength of their association with the phenotype using a Random Forests approach. As a case study, we apply our procedure to real data sets and find trustworthy numbers of, partially novel, genomic variants and genes involved in various egg quality traits.

POMC gene expression, polymorphism, and the association with reproduction traits in chickens

Reproduction trait is one of the most important economic traits in poultry industry. This study was aimed to investigate the mRNA expression levels, single nucleotide polymorphisms (SNP) of POMC gene, and the association with reproduction traits in chickens. Five SNP (g.958 G > A, g.1374 G > C, g.1393 G > A, g.1817 C > T, and g.1918G > A) were detected in introns of POMC gene in 317 Zhenning yellow chickens. Association analysis revealed that g.958 G > A and g.1817 C > T showed significantly associations with fertilization rate, hatching rate of hatching eggs, and hatching rate of fertilized eggs in chickens. Simultaneously, g.1374 G > C and g.1918G > A were both associated with egg weight at 300 D of age (P < 0.05). The SNP of g.958 G > A, g.1393 G > A, and g.1817 C > T were all associated with E2 hormone levels (P < 0.05). The result of mRNA expression levels in different tissues showed that POMC mRNA expression level in the pituitary was higher than those in the other tissues and varied in different genotypes. In conclusion, the results in this study provided new evidences that polymorphisms of the POMC gene have potential effects on reproduction traits in chickens. The 5 SNP detected in this study could be potential markers for improving reproduction traits in chickens.

Genome-Wide Association Study of Growth and Body-Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea) Using ddRAD Sequencing

Large yellow croaker (Larimichthys crocea) is an economically important marine fish species of China. Due to overfishing and marine pollution, the wild stocks of this croaker have collapsed in the past decades. Meanwhile, the cultured croaker is facing the difficulties of reduced genetic diversity and low growth rate. To explore the molecular markers related to the growth traits of croaker and providing the related SNPs for the marker-assisted selection, we used double-digest restriction-site associated DNA (ddRAD) sequencing to dissect the genetic bases of growth traits in a cultured population and identify the SNPs that associated with important growth traits by GWAS. A total of 220 individuals were genotyped by ddRAD sequencing. After quality control, 27,227 SNPs were identified in 220 samples and used for GWAS analysis. We identified 13 genome-wide significant associated SNPs of growth traits on 8 chromosomes, and the beta P of these SNPs ranged from 0.01 to 0.86. Through the definition of candidate regions and gene annotation, candidate genes related to growth were identified, including important regulators such as fgf18, fgf1, nr3c1, cyp8b1, fabp2, cyp2r1, ppara, and ccm2l. We also identified SNPs and candidate genes that significantly associated with body shape, including bmp7, col1a1, col11a2, and col18a1, which are also economically important traits for large yellow croaker aquaculture. The results provided insights into the genetic basis of growth and body shape in large yellow croaker population and would provide reliable genetic markers for molecular marker-assisted selection in the future. Meanwhile, the result established a basis for our subsequent fine mapping and related gene study.