Next Article in Journal
MiR-371a-5p Positively Associates with Hepatocellular Carcinoma Malignancy but Sensitizes Cancer Cells to Oxaliplatin by Suppressing BECN1-Dependent Autophagy
Previous Article in Journal
The Performance and Mechanism of Sludge Reduction by the Bioaugmentation Approach
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Life
Volume 12
Issue 10
10.3390/life12101650
life-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two Retrotransposon Elements in Intron of Porcine BMPR1B Is Associated with Phenotypic Variation

by
Chenglin Chi
1,
Jia He
1,
Zhanyu Du
1,
Yao Zheng
1,
Enrico D’Alessandro
2,
Cai Chen
1,
Ali Shoaib Moawad
1,3,
Emmanuel Asare
1,
Chengyi Song
1 and
Xiaoyan Wang
1,*
1
College of Animal Science & Technology, Yangzhou University, Yangzhou 225009, China
2
Department of Veterinary Science, Division of Animal Production, University of Messina, 98168 Messina, Italy
3
Department of Animal Production, Faculty of Agriculture, Kafrelsheikh University, Kafrelsheikh 33516, Egypt
*
Author to whom correspondence should be addressed.
Life 2022, 12(10), 1650; https://doi.org/10.3390/life12101650
Submission received: 31 August 2022 / Revised: 28 September 2022 / Accepted: 16 October 2022 / Published: 20 October 2022
(This article belongs to the Special Issue Regulation Mechanism of Gene Expression Mediated by Transposon)
Browse Figures
Review Reports Versions Notes

Abstract

:
It has been established that through binding to bone morphogenetic proteins (BMPs), bone morphogenetic protein receptor I B (BMPR1B) can mediate transforming growth factor β (TGF-β) signal transduction, and is involved in the regulation of several biological processes, such as bone and muscle formation and homeostasis, as well as folliculogenesis. Also known as FecB, BMPR1B has been reported as the major gene for sheep prolificacy. A number of previous studies have analyzed the relationship between single nucleotide polymorphisms (SNPs) in this gene and its related performance. In recent years, with the illustration of the effect of retrotransposon insertion on the expression of the proximal genes or phenotypic variation, retrotransposon insertion polymorphisms (RIPs) have been used as a novel type of molecular marker in the evaluation of evolution, population structure and breeding of plant and domestic animals. In this study, the RIPs in porcine BMPR1B gene were excavated, and thereafter verified using a comparative genome and polymerase chain reaction (PCR). The potential effects of phenotype, gene expression and functions related to RIPs were also explored. The results showed that 13 distinct RIPs were identified in introns of porcine BMPR1B. Among these, only BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 displayed a close relationship with the growth traits of Large White pigs. Moreover, the total number of BMPR1B-SINE+/+-RIP9 individuals born was found to be significantly higher than that of SINE−/− (p < 0.05). These two RIPs showed an obvious distribution pattern among Chinese indigenous breeds and Western commercial breeds. The expression of BMPR1B in ovaries of adult BMPR1B-SINE+/+-RIP9 Sushan pigs was found to be significantly higher in comparison to those of BMPR1B-SINE−/−-RIP9 (p < 0.05). SINE insertion of BMPR1B-SINE-RIP9 and LINE insertion of BMPR1B-LINE-RIP13 were observed to significantly increase the activity of Octamer binding transcription factor 4 (OCT4) minipromoter in CHO and C2C12 cells (p < 0.01). Therefore, these two RIPs could serve as useful molecular markers for modulating the growth or reproductive traits in assisted selection of pig breeding, while the mechanisms of the insertion function should be studied further.

1. Introduction

Bone morphogenetic protein receptor I B (BMPR1B) is a member of the bone morphogenic protein (BMP) receptor family, which can effectively bind BMPs to mediate transforming growth factor β (TGF-β)-induced signal transduction events [1], and is involved in the regulation of important biological processes, such as bone and muscle formation and homeostasis [2,3], folliculogenesis [4], etc. The protein encoded by BMPR1B gene contains four conserved residues, which can play different roles in defining BMP2, 6, 7 and GDF5 ligand affinity [5]. The variation and expression change of the BMPR1B have been related to the human ovarian insufficiency and defects of folliculogenesis [6,7], Pierre Robin sequence [8], isolated acromesomelic dysplasia [9], neuroblastoma [10], growth plate defects [11] and cancers [2].
For domestic animals, the point mutation of A746G in the coding region of sheep BMPR1B, also known as FecB, was reported to cause a significant increase in the ovulation rate and Booroola phenotype in 2001 [12,13,14]. Since then, a number of studies have been performed to decipher the relationship between this causal mutation and the prolificacy of sheep or goats [15,16,17]. The expression difference among these mutation allele individuals has been identified, in order to further elucidate the mechanism of the A746G mutation on biological processes [18,19,20]. The gene editing technology on this point mutation was used to successfully obtain the FecB-mutated sheep [21]. However, the transgenic pigs generated by introducing 746G in BMPR1B showed a normal litter size performance [22]. In addition, other point mutation and indels in the BMPR1B gene were found to significantly affect the litter size of sheep [23,24]. Overall, this gene has been thought as the major gene involved in regulating sheep prolificacy, and conformed connecting with pig prolificacy [25]. Up to now, there have only been few studies reported on the retrotransposon insertion polymorphisms (RIPs) of BMPR1B in pig breeding.
Retrotransposons, which are the dominant type of transposable elements (TEs), can effectively contribute to one-third to one-half of the mammalian genome [26]. They are an interesting genome component that can mobilize and generate new copies of themselves, as they are reversely transcribed into the genome using the pattern of “copy and paste”. The cis-regulatory sequences contained in TEs could recruit the host-encoded factor, such as transcription factors (TFs), and ensure their amplification within the new loci of host [27]. Therefore, TE insertion can affect the proximal genes by functioning as enhancers, promoters, silencers and boundary elements contributing to the genome architecture [28,29].
According to the presence of long terminal repeats (LTRs) or non-LTR, retrotransposons can be further classified into two types. Endogenous retroviruses (ERVs) are the main part of LTRs. Non-LTR retrotransposons have been divided into long interspersed elements (LINEs) and short-interspersed elements (SINEs). Based on the retrotransposon insertion polymorphism, several types of molecular markers have been developed in plant breeding [30,31]. Recently, RIPs were also identified as a causative mutation of the phenotype, or closely related to the performance of domestic animals [32,33,34,35,36,37,38,39]. Considering the importance of the BMPR1B protein in several biological process, the RIPs in BMPR1B were excavated and verified. The potential relationship of these RIPs with pig growth, as well as the reproductive performance and the mechanism of these insertions, were identified to find the effective and novel molecular markers that could be used in pig breeding.

2. Materials and Methods

2.1. Animals and DNA Isolation

Ear or blood samples were collected from the twelve pig breeds or population including Duroc, Landrace, Large White, Sujiang, Diannan small ear, Erhualian, Wuzhishan, Bama, Tibetan, Meishan, Fengjing and wild boars. Genomic DNA was extracted using a standard phenol-chloroform extraction procedure and stored at −20 °C. DNA pool from 12 breeds or population were prepared for RIPs verification, with 6 individuals for each breed. A certain number of breeds, including Duroc, Landrace, Large White, Erhualian, Wuzhishan, Sushan and Sujiang were prepared for analyzing the population diversity. The correlation analysis was performed between the different RIP genotypes and the performance records of thickness of backfat, age at 30 kg bodyweight, age at 100 kg bodyweight, first litter size, number born alive, litter weight of Large White pigs. Both the origin of breeds and number have been shown in Table S1. Among them, Large White, Duroc and Landrace are three lean type commercial breeds. Fengjing, Meishan, Erhualian and Wuzhishan are four fat type indigenous breeds. Sujiang is a cross breed which originated from Duroc, Jiangquhai and Fengjing. Sushan is also a cross breed developed from Meishan, Erhualian, Landrace and Large White in Jiangsu Province [40].

2.2. RIPs Excavation

The sequence of BMPR1B (Gene ID: 396691) with the flanking region (5 kb 5′ upstream and 3 kb 3′ downstream) were downloaded from NCBI (https://ncbi.nlm.nih.gov/, accessed on 17 May 2021). As the bait, BMPR1B sequences from the fourteen assembled non-reference genomes (Bama, Pietrain, Jinhua, cross-breed of Yorkshire/Landrace/Duroc, Landrace, Hampshire, Wuzhishan, Berkshire, Large White, Tibetan, Meishan, Bamei, Gottingen and Rongchang) were downloaded from the NCBI database (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 19 May 2021). All the sequences were aligned to excavate structural variations using the ClustalX (Version 2.0). In addition, large structural variations (more than 50 bp and presenting more than one population) were retained for the retrotransposon annotation using RepeatMasker (Versions 4.0.7). The structural variations with over 60% similarity to the sequences in RepeatMasker were considered as the predicted retrotransposons and used for further identification by PCR.
The conservatism of BMPR1B genes from the different species, including pig (ENSSSCG00000029621), Chacoan Peccary (ENSCWAG00000015161), human (ENSG00000138696), cattle (ENSBTAG00000002081) and horse (ENSECAG00000012140), were analyzed using mVISTA (http://genome.lbl.gov/vista/mvista/submit.shtml, accessed on 21 May 2021).

2.3. RIPs Identification by PCR

The primers spanning the predicted retrotransposons were designed using Oligo7. The primer sequences and annealing temperature are shown in Table S2. All the primers in this study were synthesized by Beijing Liuhe Huada Gene Technology Co. Ltd. (Beijing, China). PCR was carried out using ABI PRISM 3100 Genetic Analyzer (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA, USA), and the reaction system was prepared as described previously by using the 12 DNA pools as template [38].

2.4. Population Diversity Based on RIPs

BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 were used to illustrate the genetic diversity of the represented breeds, including the lean-type commercial breeds (Large white, Landrace and Duroc), fat-type indigenous breeds (Fengjing, Meishan, Erhualian and Wuzhishan) and Jiangsu cross breeds (Sujiang and Sushan). Genotype frequency, allele frequency and Hardy-Weinberg test of RIPs were analyzed by POPGENE 1.32 [41]. Polymorphic information content (PIC) was calculated by the following formula:
P I C = 1 − ∑ i = 1 m P i 2 − ∑ i = 1 m − 1 ∑ j = i + 1 m 2 P i 2 P j 2

2.5. qPCR

The primers (Table S2) for the pig BMPR1B (mRNA XM_021100423.1) and GAPDH (mRNA XM_003126531.5) were designed based on the sequences available in GenBank Using Oligo7 (accessed on 3 June 2021).
Four Sushan individuals for each genotype of BMPR1B-SINE-RIP9 and for LINE+/+ and LINE+/− individuals of BMPR1B-LINE-RIP13 were slaughtered to obtain the ovary, backfat and longissimus dorsi. Thereafter, the tissue RNA was extracted using TRIzol (Tiangen, Beijing, China), and the quantity and quality of total RNA were monitored using a spectrophotometer (NanoPhotometer N60 Touch, Implen Gmbh, Munich, Germany). The cDNA was prepared using Takara Kit (Takara, Dalian, China). A real-time quantitative polymerase chain reaction (qPCR) was performed in a 25μL reaction mixture, containing 12.5 μL TB Green® premix Ex TaqTM II (2×, Biomedical, Dalian, China), 1 μL each forward and reverse primer (10 ng/µL), 8.5 µL ddH2O, and 2 µL cDNA (100 ng/µL). The reaction conditions used were as follows: pre-denaturation at 95 °C for 30 s; denaturation at 95 °C for 3 s; and extension at 60 °C for 30 s (40 cycles). The relative expression levels of BMPR1B mRNA were analyzed using the 2−ΔΔCt method, and the GAPDH gene was used as the internal reference.

2.6. Dual-Luciferase Reporter Assay

The insertion fragments with the double restriction enzyme cutting site of BMPR1B-SINE-RIP9 (273 bp) and BMPR1B-LINE-RIP13 (804 bp) were cloned from the Sushan genomic DNA (the primers have been listed in Table S2), and verified by sequencing. The fragments were then inserted into the pGL3-basic vectors (Promega, Madison, WI, USA). Furthermore, Oct4 mini promoters were cloned from the pTol2-Oct4-mCherry vector [42] and inserted into pGL3-basic vectors with or without the insertion fragment to evaluate the potential enhancer or depressor activity. The correct recombinant plasmids were designated as pGL3-SINE-Oct4-basic and pGL3-LINE-Oct4-basic. Using a Lipofectamine 3000 reagent (Invitrogen, Carlsbad, CA, USA), 2 × 104 CHO and C2C12 cells were plated and transfected with the plasmids. The cells were collected after 48 h to evaluate the luciferase activity using the dual luciferase reporter system (Promega, Madison, WI, USA) according to the manufacturer’s instructions. The luciferase activity was calculated by the following formula:
Luciferase   activity = Firefly   luciferase   activity Renilla   luciferase   activity

2.7. Statistical Analysis

All the data has been shown as mean ± standard error, and p < 0.05 was considered a significant difference. The results were processed by SPSS17.0 (SPSS, Chicago, IL, USA) using one-way analysis of variance with Duncan’s test. The association analysis was also calculated by using one-way ANOVA in the SPSS software. The general linear model was used to determine the possible association between RIPs and growth and reproductive traits. Equation 1 for growth trait and 2 for reproductive traits used were as follows: (1) yijk= u + gi + aj + sk + eijk, (2) yijk = u + gi + sk + eik, where yijk was the phenotypic observation, u was the mean, gi was the fixed effect of genotype, aj was the fixed effect of gender, sk was the fixed effect of born weight, and eijk or eik was the residual effect.

2.8. Animal Welfare

All the treatments and protocols involving animals in this study were strictly performed in accordance with the guidelines of the Animal Experiment Ethics Committee of Yangzhou University (approval number: YZUDWSY2018-12).

3. Results

3.1. Thirteen RIPs Were Verified by PCR in the Pig BMPR1B Gene

By RepeatMasker annotation, 38 SINE, 17 LINE and 6 LTR insertion polymorphisms were predicted in the BMPR1B gene among the reference genome, and another 15 non-reference genomes. The PCR was employed to identify these predicted RIPs by using the specific primers (Table S1), and 13 RIPs were verified in DNA pool (Figure 1A). They were all distributed in the introns of BMRP1B gene (Figure 1B). Most of the insertion loci were relatively not conserved (Figure 1C). The details, including the insertion sequence size, direction and subfamily of retrotransposon, have been listed in Table 1. These RIPs were named as BMPR1B-SINE-RIP1, BMPR1B-SINE-RIP2, BMPR1B-SINE-RIP3, BMPR1B-SINE-RIP4, BMPR1B-SINE-RIP5, BMPR1B-SINE-RIP6, BMPR1B-SINE-RIP7, BMPR1B-ERV-RIP8, BMPR1B-SINE-RIP9, BMPR1B-SINE-RIP10, BMPR1B-LINE-RIP11, BMPR1B-ERV-RIP12 and BMPR1B-LINE-RIP13, which corresponded to labels 1–13 in Figure 1B.

3.2. Association between RIPs and Phenotype of Large White Pigs

The association of 13 RIPs with the phenotype of Large White pigs was analyzed, and BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 were found to be closely associated with the growth and reproductive traits of Large White pigs. As shown in Table 2, the age at 100 kg body weight of SINE−/− individuals for BMPR1B-SINE-RIP9 was observed to be significantly lower than those of SINE+/− individuals (p < 0.01). The age at 30 kg body weight and 100 kg body weight of LINE–/– individuals for BMPR1B-LINE-RIP13 was noted to be significantly lower than those of LINE+/+ individuals (p < 0.05). The first litter size of SINE+/+ individuals for BMPR1B-SINE-RIP9 of Large White pigs was significantly higher as compared to those of SINE−/− individuals (p < 0.05) (Table 3).

3.3. BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 Distribution in the Different Pig Breeds

The potential distribution and genetic parameters of BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 in lean-type commercial breeds (Large White, Landrace and Duroc), fat-type indigenous breeds (Fengjing, Erhualian, Meishan, Jinhua and Wuzhishan) and cross breeds in Jiangsu province (Sujiang and Sushan) are shown in Table 4. The SINE+ frequency of BMPR1B-SINE-RIP9 was dominant in the lean-type breeds (Large White and Landrace) and cross breeds (Sushan and Sujiang) while the SINE− was dominant in fat-type breeds (Wuzhishan). The LINE+ frequency of BMPR1B-LINE-RIP13 was obviously found to be higher in lean-type breeds (Duroc, Large White and Landrace) in comparison to cross breeds (Sushan and Sujiang) and fat-type breeds (Jinhua and Erhualian). Except for Large White pigs, all RIPs in the different breeds were in Hardy-Weinberg equilibrium. BMPR1B-SINE-RIP9 in fat-type indigenous breeds (Wuzhishan) and Sujiang were low polymorphic, while lean-type breeds (Duroc, Large White and Landrace) and Sushan were medium polymorphic. BMPR1B-LINE-RIP13 in the cross pigs (Sushan and Sujiang) and fat-type breeds (Jinhua and Erhualian) were low polymorphic, but in lean-type breeds (Duroc, Large White and Landrace) they were medium polymorphic.

3.4. The Expression Pattern of Pig BMPR1B in Sushan Pigs for Two RIPs

The expression in ovary, longissimus dorsi and backfat of the different genotype individuals of six-month old Sushan pigs for BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 was shown in Figure 2. For BMPR1B-SINE-RIP9, there was significant difference in the expression of BMPR1B among SINE+/+ and SINE−/− individuals in ovaries (p < 0.05). Two genotypes, LINE+/− and LINE−/−, were found in Sushan pigs for BMPR1B-LINE-RIP13. There was no significant difference noted between these two genotype individuals for the expression of BMPR1B in the three tissues of Sushan pigs (p > 0.05).

3.5. SINE and LINE Insertion Might Act as an Enhancer in CHO and C2C12 Cells

The insertion fragments, which were 273 bp and 692 bp for BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13, respectively, were amplified from Large White genomic DNA and then inserted into the pGL3-Oct4-basic vector. The recombinant plasmids were named as BMPR1B-SINE9-Oct4-Luc+ and BMPR1B-LINE13-Oct4-Luc+ (Figure 3C). These two plasmids were transfected into CHO and C2C12 cells, and the pGL3-Oct4-basic plasmid was used as a negative control (Figure 3D). The luciferase activity of BMPR1B-SINE9-Oct4-Luc+ and BMPR1B-LINE13-Oct4-Luc+ was found to be significantly higher than that of pGL3-Oct4-basic negative control (p < 0.01) in CHO and C2C12, Therefore, the 273 bp fragment and 692 bp fragment might act as an enhancer to affect the Oct4 promoter activity in CHO and C2C12 cells.

4. Discussion

BMPR1B is a type 1B receptor and a pivotal member of the BMP receptor family. The protein encoded by the BMPR1B gene mainly participates in the regulation of the BMP signaling pathway. BMPs are secreted ligands in the TGFβ superfamily, and can perform various important roles, including neurogenesis, organogenesis, pregnancy, cancer, etc [43]. All porcine BMP genes were screened to find out useful RIPs, but nothing was found within the gene and upstream or downstream flanking in the previous reports. The similar results were also observed in GH, Leptin gene, which is necessary for maintaining the different physiological processes including growth, development and reproduction. The studies have shown that 85–90% of mouse and human protein coding genes contain transposon elements sequences in their introns [44]. Thus, it is possible that essential genes might have developed measures to control the retrotransposon mobility in order to ensure proper physiological function [45]. However, their receptors could present some RIPs, which lead to the sequence diversification of organisms. All identified RIPs were found in the intron of the porcine BMPR1B gene. In our previous study, only 0.63% of retrotransposons were inserted into the exon of protein coding genes, while 35.1% retrotransposons were inserted into the intron of protein coding genes [46]. The considerable transposon elements’ under-representation in the coding regions could be attributed to a strong negative selection to avoid repercussions in the protein structure [47].
A number of studies have reported that BMPR1B is the major gene for determining the litter size of sheep. As a new-type molecular marker, 13 RIPs were used to analyze the correlation with the growth and reproductive traits of Large White pigs in this study. It was found that only BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 were closely related to the growth traits of Large White pigs, while BMPR1B-SINE-RIP9 was closely linked to the first litter size of Large White pigs. Using high throughput sequencing technology, it has been reported that the selective sweep that included BMPR1B (chr8:134.05–134.15 Mb and chr8: 123.3–147.4 Mb), which overlapped with the QTL to BMPR1B, significantly affected the number of piglets born [25,48]. BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 were within or very close to the selective sweep for litter size on the Chr8. Thereafter, the distribution of these two RIPs were evaluated in the representative pigs of three different type pigs, including fat-type, lean-type and cross breeds. Chinese indigenous pigs are famous for their high prolificacy, good meat quality, slow growth rate and low lean percentage in comparison with commercial pigs originating from Western nations [49,50,51,52,53]. There were obvious different distribution patterns noticed among these three types of pigs. Some nucleotide variants in the intragenic and intergenic regions indicated population characteristics among the different originated pig breeds [54,55,56,57]. BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 showed the obvious variable distribution pattern among the different type of breeds, which may be attributed to the difference of artificial selection concerned with these pigs. Therefore, we focused our attention on these two RIPs.
BMPR1B defects in mice have been shown to cause infertility [58], and high expression of BMPR1B in ovaries has been related to highly prolific sheep and goats [20,59]. In this study, the SINE-insertion individuals exhibited higher expression of BMPR1B in ovary and greater first litter size in comparison to individuals without SINE insertion. The function of BMPR1B was diverse, and the protein contributed to the growth and development of the mammal. The two RIPs were also found to regulate the growth rate of Large White pigs. However, the expression difference was not found in back fat and muscle among these genotypes. It is possible that the expression of BMPR1B in the tissues of brain was relatively high compared with other tissues, including the ovary [20,23], which made it play a major role through modulating the hypothalamic–pituitary axis to affect the growth and development.
To explore the potential mechanism of retrotransposon insertion affecting the expression of BMPR1B, the dual-luciferase reporter assay was employed to evaluate the potential activities of insertion sequences. Retrotransposon can disperse vast amounts of promoters, enhancers, and even repressive elements [60,61]. In this study, the SINE and LINE insertion significantly increased the dual-luciferase activity in CHO and C2C12 cells. Since SINE insert is present upstream of CDS domain, it could act as an enhancer to increase the expression of BMPR1B in ovaries and affect the reproductive traits of Large White pigs. However, the mechanisms through which BMPR1B-LINE-RIP13 can affect growth traits is still unclear, but the function of LINE insertion should be further explored. Overall, combined analysis of the relationship between these two RIPs and related traits of Large White pigs could lead to useful molecular markers for assisted selection in pig breeding.

5. Conclusions

Thirteen distinct RIPs were identified in introns of porcine BMPR1B by comparative genome and PCR. Only BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 displayed a close relationship with the growth rate or first litter size of Large White pigs. These two RIPs exhibited an obvious distribution pattern that the SINE+ frequency of BMPR1B-SINE-RIP9 and the LINE+ frequency of BMPR1B-LINE-RIP13 was dominant in the commercial breeds in comparison to Chinese indigenous breeds, which may be attributed to the difference of artificial selection concerned with these pigs. The expression of BMPR1B in ovaries was altered significantly by the SINE insertion of BMPR1B-SINE-RIP9 in adult Sushan pigs (p < 0.05). SINE insertion of BMPR1B-SINE-RIP9 and LINE insertion of BMPR1B-LINE-RIP13 significantly increased the activity of OCT4 minipromoter in both CHO and C2C12 cells (p < 0.01). Therefore, these two RIPs could serve as useful molecular markers for predicting the growth or reproductive traits in the assisted selection of pig breeding.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life12101650/s1, Table S1. Number and origin of pig breeds for BMPR1B RIPs detection. Table S2. The primers for PCR, vectors construction and q-PCR.

Author Contributions

Designing research, C.C. (Chenglin Chi), J.H., Z.D., Y.Z., E.D., C.C. (Cai Chen), A.S.M., E.A., C.S. and X.W.; performing, analyzing, and interpreting all experiments, C.C. (Chenglin Chi), X.W., C.C. (Cai Chen) and C.S.; writing the paper, C.C. (Chenglin Chi) and X.W.; and approving the manuscript, E.D., A.S.M., C.S. and X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant numbers 31872977 to Chengyi Song and 32002146 to Cai Chen); the Independent Innovation Fund Project of Agricultural Science and Technology in Jiangsu Province (CX(19)2016) to Chengyi Song; the Revitalization of Seed Industry (JBGS) in Jiangsu Province (JBGS(2021)028) to Chengyi Song; the Priority Academic Program Development of Jiangsu Higher Education Institutions; and the High-End Talent Support Program of Yangzhou University to Chengyi Song.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Animal Experiment Ethics Committee of Yangzhou University (no. NSFC2020-dkxy-02, 27 March 2020).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yamaji, N.; Celeste, A.J.; Thies, R.S.; Song, J.J.; Bernier, S.M.; Goltzman, D.; Lyons, K.M.; Nove, J.; Rosen, V.; Wozney, J.M. A mammalian serine/threonine kinase receptor specifically binds BMP-2 and BMP-4. Biochem. Biophys. Res. Commun. 1994, 205, 1944–1951. [Google Scholar] [CrossRef] [PubMed]
  2. Yoon, B.S.; Ovchinnikov, D.A.; Yoshii, I.; Mishina, Y.; Behringer, R.R.; Lyons, K.M. Bmpr1a and Bmpr1b have overlapping functions and are essential for chondrogenesis in vivo. Proc. Natl. Acad. Sci. USA 2005, 102, 5062–5067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Sartori, R.; Sandri, M. BMPs and the muscle–bone connection. Bone 2015, 80, 37–42. [Google Scholar] [CrossRef] [PubMed]
  4. Souza, C.J.; Campbell, B.K.; McNeilly, A.S.; Baird, D.T. Bone morphogenetic proteins and folliculogenesis: Lessons from the Booroola mutation. Reprod. Suppl. 2003, 61, 361–370. [Google Scholar] [CrossRef] [PubMed]
  5. Evangelista, W.; Yeh, L.C.; Gmyrek, A.; Lee, J.C.; Lee, J.C. Long-Range Communication Network in the Type 1B Bone Morphogenetic Protein Receptor. Biochemistry 2015, 54, 7079–7088. [Google Scholar] [CrossRef]
  6. Renault, L.; Patiño, L.C.; Magnin, F.; Delemer, B.; Young, J.; Laissue, P.; Binart, N.; Beau, I. BMPR1A and BMPR1B missense mutations cause primary ovarian insufficiency. J. Clin. Endocrinol. Metab. 2020, 105, e1449–e1457. [Google Scholar] [CrossRef]
  7. França, M.M.; Mendonca, B.B. Genetics of ovarian insufficiency and defects of folliculogenesis. Best Pract. Res. Clin. Endocrinol. Metab. 2021, 36, 101594. [Google Scholar] [CrossRef]
  8. Yang, Y.; Yuan, J.; Yao, X.; Zhang, R.; Yang, H.; Zhao, R.; Guo, J.; Jin, K.; Mei, H.; Luo, Y.; et al. BMPR1B mutation causes Pierre Robin sequence. Oncotarget 2017, 8, 25864–25871. [Google Scholar] [CrossRef] [Green Version]
  9. Khan, S.; Basit, S.; Khan, M.A.; Muhammad, N.; Ahmad, W. Genetics of human isolated acromesomelic dysplasia. Eur. J. Med. Genet. 2016, 59, 198–203. [Google Scholar] [CrossRef]
  10. Alshangiti, A.M.; Wyatt, S.L.; McCarthy, E.; Collins, L.M.; Hegarty, S.V.; Sullivan, A.M.; O′Keeffe, G.W. Association of distinct type 1 bone morphogenetic protein receptors with different molecular pathways and survival outcomes in neuroblastoma. Neuronal Signal 2020, 4, NS20200006. [Google Scholar] [CrossRef]
  11. Long, F.; Shi, H.; Li, P.; Guo, S.; Ma, Y.; Wei, S.; Li, Y.; Gao, F.; Gao, S.; Wang, M.; et al. A SMOC2 variant inhibits BMP signaling by competitively binding to BMPR1B and causes growth plate defects. Bone 2021, 142, 115686. [Google Scholar] [CrossRef] [PubMed]
  12. Mulsant, P.; Lecerf, F.; Fabre, S.; Schibler, L.; Monget, P.; Lanneluc, I.; Pisselet, C.; Riquet, J.; Monniaux, D.; Callebaut, I.; et al. Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Merino ewes. Proc. Natl. Acad. Sci. USA 2001, 98, 5104–5109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Souza, C.J.H.; MacDougall, C.; Campbell, B.K.; McNeilly, A.S.; Baird, D.T. The Booroola (FecB) phenotype is associated with a mutation in the bone morphogenetic receptor type 1 B (BMPR1B) gene. J. Endocrinol. 2001, 169, R1–R6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wilson, T.; Wu, X.Y.; Juengel, J.L.; Ross, I.K.; Lumsden, J.M.; Lord, E.A.; Dodds, K.G.; Walling, G.A.; McEwan, J.C.; O′Connell, A.R.; et al. Highly prolific Booroola sheep have a mutation in the intracellular kinase domain of bone morphogenetic protein IB receptor (ALK-6) that is expressed in both oocytes and granulosa cells. Biol. Reprod. 2001, 64, 1225–1235. [Google Scholar] [CrossRef] [PubMed]
  15. Fabre, S.; Pierre, A.; Mulsant, P.; Bodin, L.; Di Pasquale, E.; Persani, L.; Monget, P.; Monniaux, D. Regulation of ovulation rate in mammals: Contribution of sheep genetic models. Reprod. Biol. Endocrinol. 2006, 4, 20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Hua, G.H.; Yang, L.G. A review of research progress of FecB gene in Chinese breeds of sheep. Anim. Reprod. Sci. 2009, 116, 1–9. [Google Scholar] [CrossRef]
  17. Li, D.; Zhang, L.; Wang, Y.; Chen, X.; Li, F.; Yang, L.; Cui, J.; Li, R.; Cao, B.; An, X.; et al. FecB mutation and litter size are associated with a 90-base pair deletion in BMPR1B in East Friesian and Hu crossbred sheep. Anim. Biotechnol. 2022, 5, 1–10. [Google Scholar] [CrossRef]
  18. Tang, J.; Hu, W.; Di, R.; Liu, Q.; Wang, X.; Zhang, X.; Zhang, J.; Chu, M. Expression analysis of the prolific candidate genes, BMPR1B, BMP15, and GDF9 in Small Tail Han ewes with three fecundity (FecB gene) genotypes. Animals 2018, 8, 166. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, S.; Guo, X.; He, X.; Di, R.; Zhang, X.; Zhang, J.; Wang, X.; Chu, M. Insight into Pituitary lncRNA and mRNA at Two Estrous Stages in Small Tail Han Sheep with Different FecB Genotypes. Front. Endocrinol. 2021, 12, 789564. [Google Scholar] [CrossRef]
  20. Wen, Y.-L.; Guo, X.-F.; Ma, L.; Zhang, X.-S.; Zhang, J.-L.; Zhao, S.-G.; Chu, M.-X. The expression and mutation of BMPR1B and its association with litter size in small-tail Han sheep (Ovis aries). Arch. Anim. Breed. 2021, 64, 211–221. [Google Scholar] [CrossRef]
  21. Zhou, S.; Ding, Y.; Liu, J.; Liu, Y.; Zhao, X.; Li, G.; Zhang, C.; Li, C.; Wang, Y.; Kalds, P.; et al. Highly efficient generation of sheep with a defined FecB(B) mutation via adenine base editing. Genet. Sel. Evol. 2020, 52, 35. [Google Scholar] [CrossRef] [PubMed]
  22. Zhao, X.; Yang, Q.; Zhao, K.; Jiang, C.; Ren, D.; Xu, P.; He, X.; Liao, R.; Jiang, K.; Ma, J.; et al. Production of Transgenic Pigs with an Introduced Missense Mutation of the Bone Morphogenetic Protein Receptor Type IB Gene Related to Prolificacy. Asian-Australas. J. Anim. Sci. 2016, 29, 925–937. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Jia, J.; Chen, Q.; Gui, L.; Jin, J.; Li, Y.; Ru, Q.; Hou, S. Association of polymorphisms in bone morphogenetic protein receptor-1B gene exon-9 with litter size in Dorset, Mongolian, and Small Tail Han ewes. Asian-Australas. J. Anim. Sci. 2019, 32, 949–955. [Google Scholar] [CrossRef] [PubMed]
  24. Li, H.; Xu, H.; Akhatayeva, Z.; Liu, H.; Lin, C.; Han, X.; Lu, X.; Lan, X.; Zhang, Q.; Pan, C. Novel indel variations of the sheep FecB gene and their effects on litter size. Gene 2021, 767, 145176. [Google Scholar] [CrossRef] [PubMed]
  25. Li, W.-T.; Zhang, M.-M.; Li, Q.-G.; Tang, H.; Zhang, L.-F.; Wang, K.-J.; Zhu, M.-Z.; Lu, Y.-F.; Bao, H.-G.; Zhang, Y.-M. Whole-genome resequencing reveals candidate mutations for pig prolificacy. Proc. R. Soc. B Biol. Sci. 2017, 284, 20172437. [Google Scholar] [CrossRef] [Green Version]
  26. Platt, R.N.; Vandewege, M.W.; Ray, D.A. Mammalian transposable elements and their impacts on genome evolution. Chromosome Res. 2018, 26, 25–43. [Google Scholar] [CrossRef] [Green Version]
  27. Hermant, C.; Torres-Padilla, M.E. TFs for TEs: The transcription factor repertoire of mammalian transposable elements. Genes Dev. 2021, 35, 22–39. [Google Scholar] [CrossRef]
  28. Sundaram, V.; Wysocka, J. Transposable elements as a potent source of diverse cis-regulatory sequences in mammalian genomes. Philos. Trans. R. Soc. B 2020, 375, 20190347. [Google Scholar] [CrossRef] [Green Version]
  29. Li, D.; He, M.; Tang, Q.; Tian, S.; Zhang, J.; Li, Y.; Wang, D.; Jin, L.; Ning, C.; Zhu, W. Comparative 3D genome architecture in vertebrates. BMC Biol. 2022, 20, 99. [Google Scholar] [CrossRef]
  30. Kalendar, R.; Schulman, A.H. Transposon-based tagging: IRAP, REMAP, and iPBS. In Molecular Plant Taxonomy; Springer: Berlin/Heidelberg, Germany, 2014; Volume 1115, pp. 233–255. [Google Scholar]
  31. Amiteye, S. Basic concepts and methodologies of DNA marker systems in plant molecular breeding. Heliyon 2021, 7, e08093. [Google Scholar] [CrossRef]
  32. Wang, Z.; Qu, L.; Yao, J.; Yang, X.; Li, G.; Zhang, Y.; Li, J.; Wang, X.; Bai, J.; Xu, G.; et al. An EAV-HP insertion in 5′ Flanking region of SLCO1B3 causes blue eggshell in the chicken. PLoS Genet. 2013, 9, e1003183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Ballif, B.; Ramirez, C.; Carl, C.; Sundin, K.; Krug, M.; Zahand, A.; Shaffer, L.; Flores-Smith, H. The PMEL Gene and Merle in the Domestic Dog: A Continuum of Insertion Lengths Leads to a Spectrum of Coat Color Variations in Australian Shepherds and Related Breeds. Cytogenet. Genome Res. 2018, 156, 22–34. [Google Scholar] [CrossRef] [PubMed]
  34. Li, J.; Davis, B.W.; Jern, P.; Dorshorst, B.J.; Siegel, P.B.; Andersson, L. Characterization of the endogenous retrovirus insertion in CYP19A1 associated with henny feathering in chicken. Mob. DNA 2019, 10, 38. [Google Scholar] [CrossRef] [PubMed]
  35. Chen, J.; Dalirsefat, S.B.; Han, D.; Dong, X.; Hua, G.; Zheng, X.; Xia, T.; Shao, T.; Deng, X.; Wu, C. An EAV-HP insertion in the 5 flanking region of SLCO1B3 is associated with its tissue-expression profile in blue-eggshell Yimeng chickens (Gallus gallus)-ScienceDirect. Poult. Sci. 2020, 99, 6371–6377. [Google Scholar] [CrossRef] [PubMed]
  36. Chen, C.; Zheng, Y.; Wang, M.; Murani, E.; D’Alessandro, E.; Moawad, A.S.; Wang, X.; Wimmers, K.; Song, C. SINE Insertion in the Intron of Pig GHR May Decrease Its Expression by Acting as a Repressor. Animals 2021, 11, 1871. [Google Scholar] [CrossRef]
  37. Liang, D.; Zhao, P.; Si, J.; Fang, L.; Pairo-Castineira, E.; Hu, X.; Xu, Q.; Hou, Y.; Gong, Y.; Liang, Z. Genomic analysis revealed a convergent evolution of LINE-1 in coat color: A case study in water buffaloes (Bubalus bubalis). Mol. Biol. Evol. 2021, 38, 1122–1136. [Google Scholar] [CrossRef]
  38. Wang, X.; Chen, Z.; Murani, E.; D′Alessandro, E.; An, Y.; Chen, C.; Li, K.; Galeano, G.; Wimmers, K.; Song, C. A 192 bp ERV fragment insertion in the first intron of porcine TLR6 may act as an enhancer associated with the increased expressions of TLR6 and TLR1. Mob. DNA 2021, 12, 20. [Google Scholar] [CrossRef]
  39. Du, Z.; D′Alessandro, E.; Zheng, Y.; Wang, M.; Chen, C.; Wang, X.; Song, C. Retrotransposon Insertion Polymorphisms (RIPs) in Pig Coat Color Candidate Genes. Animals 2022, 12, 969. [Google Scholar] [CrossRef]
  40. Wang, X.; D′Alessandro, E.; Chi, C.; Moawad, A.S.; Zong, W.; Chen, C.; Song, C. Genetic Evaluation and Population Structure of Jiangsu Native Pigs in China Revealed by SINE Insertion Polymorphisms. Animals 2022, 12, 1345. [Google Scholar] [CrossRef]
  41. Yeh, F.; Yang, R.; Boyle, T. POPGENE Version 1.32: Microsoft Windows–Based Freeware for Population Genetic Analysis, Quick User Guide; Center for International Forestry Research, University of Alberta: Edmonton, AB, Canada, 1999; pp. 1–29. [Google Scholar]
  42. Chan, S.; Shen, D.; Sang, Y.; Wang, S.; Wang, Y.; Chen, C.; Gao, B.; Song, C. Development of enhancer-trapping and-detection vectors mediated by the Tol2 transposon in zebrafish. PeerJ 2019, 7, e6862. [Google Scholar] [CrossRef]
  43. Wozney, J.M.; Rosen, V.; Celeste, A.J.; Mitsock, L.M.; Whitters, M.J.; Kriz, R.W.; Hewick, R.M.; Wang, E.A. Novel regulators of bone formation: Molecular clones and activities. Science 1988, 242, 1528–1534. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, Y.; Romanish, M.T.; Mager, D.L. Distributions of transposable elements reveal hazardous zones in mammalian introns. PLoS Comput. Biol. 2011, 7, e1002046. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Malicki, M.; Iliopoulou, M.; Hammann, C. Retrotransposon Domestication and Control in Dictyostelium discoideum. Front. Microbiol. 2017, 8, 1869. [Google Scholar] [CrossRef]
  46. Chen, C.; Wang, W.; Wang, X.; Shen, D.; Wang, S.; Wang, Y.; Gao, B.; Wimmers, K.; Mao, J.K.; Song, C. Retrotransposons evolution and impact on lncRNA and protein coding genes in pigs. Mob. DNA 2019, 10, 19. [Google Scholar] [CrossRef] [PubMed]
  47. Philippsen, G.S.; DeMarco, R. Impact of transposable elements in the architecture of genes of the human parasite Schistosoma mansoni. Mol. Biochem. Parasitol. 2019, 228, 27–31. [Google Scholar] [CrossRef]
  48. Hernandez, S.; Finlayson, H.; Ashworth, C.; Haley, C.; Archibald, A. A genome-wide linkage analysis for reproductive traits in F 2 L arge W hite× M eishan cross gilts. Anim. Genet. 2014, 45, 191–197. [Google Scholar] [CrossRef] [Green Version]
  49. Legault, C. Selection of breeds, strains and individual pigs for prolificacy. J. Reprod. Fertil. Suppl. 1985, 33, 151–166. [Google Scholar] [CrossRef]
  50. Liu, X.; Xiong, X.; Yang, J.; Zhou, L.; Yang, B.; Ai, H.; Ma, H.; Xie, X.; Huang, Y.; Fang, S. Genome-wide association analyses for meat quality traits in Chinese Erhualian pigs and a Western Duroc×(Landrace× Yorkshire) commercial population. Genet. Sel. Evol. 2015, 47, 44. [Google Scholar] [CrossRef] [Green Version]
  51. Liu, Y.; Yang, X.; Jing, X.; He, X.; Wang, L.; Liu, Y.; Liu, D. Transcriptomics analysis on excellent meat quality traits of skeletal muscles of the Chinese indigenous min pig compared with the large white breed. Int. J. Mol. Sci. 2017, 19, 21. [Google Scholar] [CrossRef] [Green Version]
  52. Ma, X.; Li, P.; Zhu, M.; He, L.; Sui, S.; Gao, S.; Su, G.; Ding, N.; Huang, Y.; Lu, Z. Genome-wide association analysis reveals genomic regions on Chromosome 13 affecting litter size and candidate genes for uterine horn length in Erhualian pigs. Animal 2018, 12, 2453–2461. [Google Scholar] [CrossRef]
  53. Chen, Q.; Zhang, W.; Cai, J.; Ni, Y.; Xiao, L.; Zhang, J. Transcriptome analysis in comparing carcass and meat quality traits of Jiaxing Black Pig and Duroc× Duroc× Berkshire× Jiaxing Black Pig crosses. Gene 2022, 808, 145978. [Google Scholar] [CrossRef] [PubMed]
  54. Zhang, Y.; Du, H.; Chen, J.; Yang, G.; Zhang, X. Porcine growth differentiation factor 9 gene polymorphisms and their associations with litter size. J. Genet. Genom. 2008, 35, 163–169. [Google Scholar] [CrossRef]
  55. Tao, H.; Mei, S.; Sun, X.; Peng, X.; Zhang, X.; Ma, C.; Wang, L.; Hua, L.; Li, F. Associations of TCF12, CTNNAL1 and WNT10B gene polymorphisms with litter size in pigs. Anim. Reprod. Sci. 2013, 140, 189–194. [Google Scholar] [CrossRef] [PubMed]
  56. Charoensook, R.; Gatphayak, K.; Brenig, B.; Knorr, C. Genetic diversity analysis of Thai indigenous pig population using microsatellite markers. Asian-Australas. J. Anim. Sci. 2019, 32, 1491. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Lee, S.H.; Seo, D.W.; Cho, E.S.; Choi, B.H.; Kim, Y.M.; Hong, J.K.; Han, H.D.; Jung, Y.B.; Kim, D.J.; Choi, T.J. Genetic diversity and ancestral study for Korean Native Pigs using 60K SNP chip. Animals 2020, 10, 760. [Google Scholar] [CrossRef] [PubMed]
  58. Mishina, Y.; Suzuki, A.; Ueno, N.; Behringer, R.R. Bmpr encodes a type I bone morphogenetic protein receptor that is essential for gastrulation during mouse embryogenesis. Genes Dev. 1995, 9, 3027–3037. [Google Scholar] [CrossRef] [Green Version]
  59. Pan, Z.Y.; Di, R.; Tang, Q.Q.; Jin, H.H.; Chu, M.X.; Huang, D.W.; He, J.N.; Liu, Q.Y.; Hu, W.P.; Wang, X.Y.; et al. Tissue-specific mRNA expression profiles of GDF9, BMP15, and BMPR1B genes in prolific and non-prolific goat breeds. Czech J. Anim. Sci. 2015, 60, 452–458. [Google Scholar] [CrossRef] [Green Version]
  60. Chuong, E.B.; Elde, N.C.; Feschotte, C. Regulatory activities of transposable elements: From conflicts to benefits. Nat. Rev. Genet. 2017, 18, 71–86. [Google Scholar] [CrossRef] [Green Version]
  61. Bourque, G.; Burns, K.H.; Gehring, M.; Gorbunova, V.; Seluanov, A.; Hammell, M.; Imbeault, M.; Izsvak, Z.; Levin, H.L.; Macfarlan, T.S.; et al. Ten things you should know about transposable elements. Genome Biol. 2018, 19, 199. [Google Scholar] [CrossRef]
Life 12 01650 g001 550
Figure 1. RIPs identification and distribution in BMPR1B gene. (A). RIPs identification by PCR and electrophoresis in DNA pool of twelve population. M: DL2000 marker; 1: Duroc, 2: Landrace, 3: Large White, 4: Bama, 5: Diannan small-ear, 6: Fengjing, 7: Sujiang, 8: Wuzhishan, 9: Tibetan, 10: Meishan, 11: Erhualian, 12: Wild boars; (B). RIPs Insertion loci diagram in BMPR1B; and (C). Sequence alignment for different species in BMPR1B.
Figure 1. RIPs identification and distribution in BMPR1B gene. (A). RIPs identification by PCR and electrophoresis in DNA pool of twelve population. M: DL2000 marker; 1: Duroc, 2: Landrace, 3: Large White, 4: Bama, 5: Diannan small-ear, 6: Fengjing, 7: Sujiang, 8: Wuzhishan, 9: Tibetan, 10: Meishan, 11: Erhualian, 12: Wild boars; (B). RIPs Insertion loci diagram in BMPR1B; and (C). Sequence alignment for different species in BMPR1B.
Life 12 01650 g001
Life 12 01650 g002 550
Figure 2. The expression pattern of pig BMPR1B in different genotype Sushan pigs of two RIPs. (A). Genotypes of BMPR1B-SINE-RIP9 in the Sushan pigs, M: DL2000 marker; (B). Genotypes of BMPR1B-LINE-RIP13 in the Sushan pigs, M: DL2000 marker; (C). Influence of SINE insertion on expression of BMPR1B gene; and (D). Effect of LINE insertion on expression of BMPR1B gene. The values are shown as mean ± SD, and * showed p < 0.05. “+/+” means Homozygous insert. “+/−” means heterozygous insert. “−/−” means no insert. The values are shown as mean ± SD, and * showed p < 0.05.
Figure 2. The expression pattern of pig BMPR1B in different genotype Sushan pigs of two RIPs. (A). Genotypes of BMPR1B-SINE-RIP9 in the Sushan pigs, M: DL2000 marker; (B). Genotypes of BMPR1B-LINE-RIP13 in the Sushan pigs, M: DL2000 marker; (C). Influence of SINE insertion on expression of BMPR1B gene; and (D). Effect of LINE insertion on expression of BMPR1B gene. The values are shown as mean ± SD, and * showed p < 0.05. “+/+” means Homozygous insert. “+/−” means heterozygous insert. “−/−” means no insert. The values are shown as mean ± SD, and * showed p < 0.05.
Life 12 01650 g002
Life 12 01650 g003 550
Figure 3. SINE and LINE insertion may act as a potential enhancer in CHO and C2C12 cells. (A). Result of target sequence and transposon reference sequence alignment; (B). The electrophoretogram of double-enzyme digestion for the recombinant vector, M: DL5000 marker; (C). A schematic diagram of the recombinant vector structure; and (D). A schematic diagram of the results of enhancer activity verification. The values shown are mean ± SD, and ** showed p < 0.01.
Figure 3. SINE and LINE insertion may act as a potential enhancer in CHO and C2C12 cells. (A). Result of target sequence and transposon reference sequence alignment; (B). The electrophoretogram of double-enzyme digestion for the recombinant vector, M: DL5000 marker; (C). A schematic diagram of the recombinant vector structure; and (D). A schematic diagram of the results of enhancer activity verification. The values shown are mean ± SD, and ** showed p < 0.01.
Life 12 01650 g003
Table
Table 1. The details of retrotransposon insertion in the BMPR1B gene.
Table 1. The details of retrotransposon insertion in the BMPR1B gene.
LociInsertion in Exon/IntronLength (bp)DirectionTypeChromosome (Sus scrofa11.1)
BMPR1B-SINE-RIP1Intron1286−SINEA1125,030,204–125,030,205
BMPR1B-SINE-RIP2Intron1307+SINEA1125,020,478–125,020,478
BMPR1B-SINE-RIP3Intron1259−SINEA3124,975,543–124,975,544
BMPR1B-SINE-RIP4Intron1283+SINEA1124,934,102–124,934,103
BMPR1B-SINE-RIP5Intron2304+SINEA1124,833,678–124,833,679
BMPR1B-SINE-RIP6Intron2291+SINEA1124,831,282–124,831,283
BMPR1B-SINE-RIP7Intron2285−SINEA1124,829,563–124,829,564
BMPR1B-ERV-RIP8Intron2289+ERVIII124,804,706–124,804,995
BMPR1B-SINE-RIP9Intron2273−SINEA1124,800,210–124,800,483
BMPR1B-SINE-RIP10Intron2315+SINEA1124,788,936–124,788,937
BMPR1B-LINE-RIP11Intron5134−L1D11124,607,615–124,607,749
BMPR1B-ERV-RIP12Intron5565−ERVIII124,606,427–124,606,428
BMPR1B-LINE-RIP13Intron5692+L1D21124,595,110–124,595,802
Note: “+” indicates forward insertion, “−” indicates reverse insertion.
Table
Table 2. Correlation between BMPR1B-SINE-RIP9, BMPR1B-LINE-RIP13 and growth traits of Large White pigs.
Table 2. Correlation between BMPR1B-SINE-RIP9, BMPR1B-LINE-RIP13 and growth traits of Large White pigs.
RIPsGenotypeNumberThickness of Back Fat/mmAge at 30 kg Body Weight/dAge at 100 kg Body Weight/d
BMPR1B-SINE-RIP9SINE+\+16810.51 ± 0.2074.99 ± 0.71162.14 ± 0.55 a
SINE+\−28610.87 ± 0.1774.91 ± 0.52163.71 ± 0.52 A
SINE−\−7511.25 ± 0.3474.89 ± 1.27160.80 ± 0.90 b
BMPR1B-LINE-RIP13LINE+/+25910.65 ± 0.1774.92 ± 0.52 a163.25 ± 0.56 a
LINE+/−11510.67 ± 0.2874.36 ± 0.94 a161.97 ± 0.76 ab
LINE−/−2111.62 ± 0.6270.07 ± 2.58 b159.26 ± 1.48 b
Note: +/+: homozygote with retrotransposon insertion; +/−: heterozygote with retrotransposon insertion; −/−: homozygote without retrotransposon insertion. The same superscript letters in the same column (a and ab or b and ab) mean that the difference between groups was not significant. Different superscript lowercase letters (a and b) indicated a significant difference between groups (p < 0.05). Different superscript capital letters (A and b) indicated an extremely significant difference between groups (p < 0.01).
Table
Table 3. Correlation between BMPR1B-SINE-RIP9 and reproductive performance of Large White pigs.
Table 3. Correlation between BMPR1B-SINE-RIP9 and reproductive performance of Large White pigs.
GenotypeNumberFirst Litter SizeNumber Born AliveBorn Weight/kg
SINE+/+9910.47 ± 0.28 a9.86 ± 0.2613.42 ± 0.36
SINE+/−11810.03 ± 0.23 ab9.49 ± 0.2412.73 ± 0.33
SINE−/−269.19 ± 0.57 b8.88 ± 0.6311.85 ± 0.90
Note: +/+: homozygote with retrotransposon insertion; +/−: heterozygote with retrotransposon insertion; −/−: homozygote without retrotransposon insertion. The same superscript letters in the same column (a and ab or b and ab) mean that the difference between groups was not significant. Different superscript lowercase letters (a and b) indicated a significant difference between groups (p < 0.05).
Table
Table 4. Distribution of RIPs in different pig breeds.
Table 4. Distribution of RIPs in different pig breeds.
Polymorphic SiteBreedsNumberGenotype FrequencyAllele FrequencyHardy-Weinberg EquilibriumPolymorphic Information Content
+/++/−−/−+−
BMPR1B-SINE-RIP9Large White52931.7654.0614.1858.7941.210.010.367
Landrace2450.0037.5012.5068.7531.250.530.337
Sushan2433.3341.6725.0054.1745.830.430.373
Sujiang2180.9519.050.0090.489.520.630.157
Wuzhishan244.1720.8375.0014.5885.420.420.218
BMPR1B-LINE-RIP13Large White39566.5729.115.3280.1319.870.090.268
Duroc2462.5029.178.3377.0822.920.390.291
Landrace2420.8337.5041.6739.5860.420.290.364
Sushan32015.6384.387.8192.190.630.134
Sujiang214.7619.0576.1914.2985.710.310.215
Jinhua24016.6783.338.3391.670.660.141
Erhualian2404.1795.832.0897.920.920.040
Note: +/+: homozygote with retrotransposon insertion; +/−: heterozygote with retrotransposon insertion; −/−: homozygote without retrotransposon insertion; +: allele with retrotransposon insertion; and −: allele without retrotransposon insertion.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chi, C.; He, J.; Du, Z.; Zheng, Y.; D’Alessandro, E.; Chen, C.; Moawad, A.S.; Asare, E.; Song, C.; Wang, X. Two Retrotransposon Elements in Intron of Porcine BMPR1B Is Associated with Phenotypic Variation. Life 2022, 12, 1650. https://doi.org/10.3390/life12101650

AMA Style

Chi C, He J, Du Z, Zheng Y, D’Alessandro E, Chen C, Moawad AS, Asare E, Song C, Wang X. Two Retrotransposon Elements in Intron of Porcine BMPR1B Is Associated with Phenotypic Variation. Life. 2022; 12(10):1650. https://doi.org/10.3390/life12101650

Chicago/Turabian Style

Chi, Chenglin, Jia He, Zhanyu Du, Yao Zheng, Enrico D’Alessandro, Cai Chen, Ali Shoaib Moawad, Emmanuel Asare, Chengyi Song, and Xiaoyan Wang. 2022. "Two Retrotransposon Elements in Intron of Porcine BMPR1B Is Associated with Phenotypic Variation" Life 12, no. 10: 1650. https://doi.org/10.3390/life12101650

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop