• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    BMPR-IB gene disruption causes severe limb deformities in pigs

    2022-06-07 10:50:08QiangYangChuanMinQiaoWeiWeiLiuHaoYunJiangQiQiJingYaYaLiaoJunRenYuYunXing
    Zoological Research 2022年3期

    Qiang Yang, Chuan-Min Qiao, Wei-Wei Liu, Hao-Yun Jiang, Qi-Qi Jing, Ya-Ya Liao, Jun Ren, Yu-Yun Xing,*

    1 State Key Laboratory of Pig Genetic Improvement and Production Technology, Jiangxi Agricultural University, Nanchang, Jiangxi 330045, China

    ABSTRACT In an attempt to generate g.A746G substitution in the BMPR-IB gene, we unexpectedly obtained BMPR-IB homozygous knockout piglets (BMPR-IB-/-) and heterogeneous knockout piglets with one copy of the A746G mutation (BMPR-IB-/746G) via CRISPR/Cas9 editing.Polymerase chain reaction (PCR) and sequencing revealed complex genomic rearrangements in the target region.All BMPR-IB-disrupted piglets showed an inability to stand and walk normally.Both BMPR-IB-/- and BMPR-IB-/746G piglets exhibited severe skeletal dysplasia characterized by distorted and truncated forearms(ulna, radius) and disordered carpal, metacarpal, and phalangeal bones in the forelimbs.The piglets displayed more severe deformities in the hindlimbs by visual inspection, including fibular hemimelia,enlarged tarsal bone, and disordered toe joint bones.Limb deformities were more profound in BMPR-IB-/-piglets than in the BMPR-IB-/746G piglets.Proteomic analysis identified 139 differentially expressed proteins (DEPs) in the hindlimb fibula of BMPR-IB-/746G piglets compared to the wild-type (WT)controls.Most DEPs are involved in skeletal or embryonic development and/or the TGF-β pathway and tumor progression.Gene Ontology (GO) and protein domain enrichment analysis suggested alterations in these processes.Of the top 50 DEPs, a large proportion, e.g., C1QA, MYO1H, SRSF1,P3H1, GJA1, TCOF1, RBM10, SPP2, MMP13, and PHAX, were significantly associated with skeletal development.Our study provides novel findings on the role of BMPR-IB in mammalian limb development.

    Keywords: BMPR-IB; A746G; Pigs; Limb deformities

    INTRODUCTION

    Bone morphogenetic proteins (BMPs), members of the transforming growth factor β (TGF-β) family, play important roles in the formation of bone and cartilage and the development of other organs, such as muscle, kidney, and blood vessels (Katagiri & Watabe, 2016).BMPs transduce their signals through type I and type II serine-threonine kinase receptors (BMPRI and BMPRII) (Miyazono et al., 2010).Perturbations of BMP signaling via BMPRI have been linked to multiple diseases in bone, cartilage, and muscles (Lin et al.,2016).As a type of BMPRI, activin-like kinase 6 (ALK6, also calledBMPR-IB) is a critical regulator of chondrogenesis and osteogenesis (Lin et al., 2016).Patients with missense (I200K and R486W) or deletion (del359-366) mutations inBMPR-IBsuffer severe limb deformations, including short stature, fibula aplasia, severe brachydactyly, and ulnar deviation of thehands, which are mainly caused by chondrodysplasia during skeletal development (Demirhan et al., 2005; Lehmann et al.,2003).In addition, null mutations ofBMPR-IBin mice result in limb abnormalities, demonstrating thatBMPR-IBis required for chondrocyte proliferation, differentiation, and maturation(Baur et al., 2000; Yi et al., 2000).In vitrostudies have shown that continuous expression of activeBMPR-IBinduces mineralized bone matrix formation, while inhibition of endogenousBMPR-IBblocksBMP2-induced osteoblast differentiation and mineralized bone matrix formation,suggesting thatBMPR-IBis required for osteoblast differentiation and bone formation (Chen et al., 1998).Mice expressing truncated dominant-negativeBMPR-IBin target osteoblasts exhibit impaired postnatal bone formation (Zhao et al., 2002).Thus, osteoblasticBMPR-IBappears to play a necessary role during postnatal bone modeling and remodeling.

    TheBMPR-IBgene also affects prolificacy in sheep and plays a vital role in the control of follicular growth and development (Davis et al., 2006; Reader et al., 2012).The A746G mutation inBMPR-IBis reported to be highly associated with increased ovulation rates and litter size in sheep (Mulsant et al., 2001; Wilson et al., 2001).As theBMPR-IBA746G mutation was not detected in any pig breeds in our collection, we introduced the 746 GG mutation into the porcine genome via traditional transgenic technology in a previous study (Zhao et al., 2016), aiming to improve reproductive performance.Though the resulting transgenic boar exhibited stronger spermatogenic ability (Zhao et al.,2016), an important functional gene,NAGLU, was disrupted due to random insertion (Yang et al., 2018).

    To avoid random integration of exogenous vectors, we used the CRISPR/Cas9 genome-editing platform (Cong et al., 2013)to edit AA746 into 746GG of theBMPR-IBgene in porcine fetal fibroblasts (PFFs) using a polymerase chain reaction(PCR)-amplified donor, with cloned pigs then generated by somatic cell nuclear transfer (SCNT).Unexpectedly, all nine cloned piglets from four deliveries showed severe limb deformities and were unable to walk and stand normally.Based on these ubiquitous limb deformities, we proposed that accidental gene disruption may have occurred.In this study,we first genotyped the target region across the donor sequence and found that theBMPR-IBgene was disrupted in these piglets.We then dissected the limb phenotypes and investigated molecular regulation in these pigs.

    MATERIALS AND METHODS

    Animals and ethics statement

    The PFFs were isolated from one-month-old large white pig fetuses using 200 U/mL collagenase type IV (Sigma-Aldrich,USA).Surrogate sows were housed individually according to standard procedures.All experiments involving animals were conducted in accordance with the Guidelines for the Care and Use of Laboratory Animals formulated by the Ministry of Agriculture and Rural Affairs of the People’s Republic of China.The study was approved by the Ethics Committee of Jiangxi Agricultural University.All animal operations were performed under anesthesia to minimize suffering.Sixteen cloned piglets (nine individuals with deformed limbs and seven normal individuals) aged 3-5 days were used in this experiment.

    Construction of CRISPR/Cas9 vector and double-stranded DNA donor template

    The pSpCas9(BB)-2A-GFP (PX458, #48138) plasmid used in this study was purchased from Addgene (http://www.addgene.org/CRISPR/, USA).Single-guide RNA (sgRNA) targeting exon 8 of porcineBMPR-IBwas designed using the Benchling online tool (https://www.benchling.com/).Here, 5'-TCATTGCTGCAGACATCAAA-3' was selected as the sgRNA,and the corresponding protospacer adjacent motif (PAM)sequence was GGG (Figure 1A).Paired synthesized oligonucleotides containing theBMPR-IB-sgRNA sequence(Supplementary Table S1) were ordered from Sangon Biotech Co., Ltd.(China); the DNA pair were annealed and ligated into theBbsI-digested PX458 vector to generate the recombinant plasmid.The recombinant plasmid DNA was transformed into Trans5ɑ competent cells (TransGen, China), and then extracted using an EndoFree?Plasmid Maxi Kit (Qiagen,Germany).

    Figure 1 Generation of BMPR-IB-modified PFFs and piglets

    Linear double-stranded donor DNA was obtained by bridgeoverlap-extension PCR.First, we used primer pair two(Supplementary Table S1) to amplify the region harboring the sgRNA sequence targetingBMPR-IBin large white pigs by routine PCR, and DNA of the individual with the GGC sequence at the PAM site was selected as a template for the next amplification.Two PCR assays using primer pairs three and four (Supplementary Table S1) were then performed in a thermocycler under the following conditions: 94 °C for 2 min;35 cycles at 98 °C for 10 s, 68 °C for 3 min; and 72 °C for 10 min.We mixed the two PCR products at a ratio of 1:1, and then conducted bridging PCR using primer pair five(Supplementary Table S1).The bridging PCR products were subjected to a final round of amplification using a touch-down PCR protocol (primer pair six for detection, Supplementary Table S1) under the following conditions: 94 °C for 5 min; 26 cycles at 94 °C for 30 s, 68 °C (-0.5 °C/cycle) for 45 s, 72 °C for 2 min; 14 cycles at 94 °C for 30 s, 55 °C for 45 s, 72 °C for 2 min; and 72 °C for 10 min.Finally, the amplified PCR product (donor DNA containing A746G mutation) was gel purified using a QIAquick PCR Purification Kit (Qiagen,Germany).

    PFF transfection and selection

    To obtainBMPR-IB746G mutation cell clones, 25 μg of Cas9-sgRNA plasmid, 25 μg of purified donor DNA, and 12.5 μg of PX459 v2 plasmid were co-transfected into 3×106PFF cells using the BTX ECM 2001 (USA) electroporation system (200 V, 1 ms, 3 pulses, 1 repeat).The electroporated PFFs were transferred into a 10 cm Petri Dish with growth medium containing Dulbecco’s Modified Eagle Medium (DMEM), 15%fetal bovine serum (FBS), 100 IU/mL penicillin and 100 μg/mL streptomycin, and SCR7 (Xcessbio, USA) to a final concentration of 1 μmol/L.After 30 h, PFFs were selected using 3.5 μg/mL puromycin (Sigma, Japan) for two days.The cells were then plated into 40 Petri Dishes (10 cm) at various cell densities for an additional 6-8 days of culture at 37 °C.The single-cell colonies in the Petri Dishes were collected and seeded in 24-well plates.After reaching 90% confluency, the cells in each colony were passaged in 6-well plates and subcultured at 37 °C for 48 h.About 20% of each colony was digested (56 °C, 60 min; 95 °C, 10 min) in 10 μL of lysis buffer(0.5% NP40 and 2 μg/μL of Proteinase K) to extract DNA, and the remaining cells were stored in liquid nitrogen for SCNT.Primer pair seven (Supplementary Table S1) was first used to identify the A746G locus in cell clones.The sgRNA sequence region was amplified using primer pair two (Supplementary Table S1) for A746G positive clones.The PCR products were sequenced on a 3130XL Genetic Analyzer (Applied Biosystem, USA).Colonies carrying 746GG were selected as donor cells for SCNT, and cells carrying AA746 were used as controls.

    SCNT and embryo transplantation

    Five cell colonies (Figure 1E) were used as nuclear donors to produce cloned pigs via SCNT, as described previously (Gong et al., 2004).Briefly, cumulus-oocyte complexes werecollected and matured at 38.5 °C for 20 h in maturation medium comprised of M199 (Gibco, USA) supplemented with 10% FBS, 0.01 U/mL basal follicle stimulating hormone(bFSH), 0.01 U/mL basal luteinizing hormone (bLH), 1 μg/mL estradiol, and 1% (v:v) penicillin/streptomycin.The first polar body was then aspirated from the mature oocytes using a glass pipette, and, finally, the donor cells were fused with enucleated oocytes using BTX ECM 2001 (USA) electrofusion equipment.The reconstructed embryos were cultured in embryo-development medium at 38.5 °C for 24 h and then transferred surgically into the oviducts of estrus-synchronized pig surrogates (each recipient receiving approximate 350 embryos).Pregnancy status was determined by ultrasonography 30 days after embryo transfer.All cloned piglets were delivered by natural birth.

    Genotyping cloned piglets

    Genomic DNA was extracted from the ear tissues of cloned piglets using a cell/tissue genomic DNA extraction kit(Generay Biotech, Shanghai, China).Primer pair eight(Supplementary Table S1) was used for long-range PCR amplification across the donor sequence inBMPR-IB.The PCR protocol included 25 μL of 2×Gflex PCR Buffer (TaKaRa,Japan), 1 μL of each forward and reverse primer (10 μmol/L),100 ng of genomic DNA, and ddH2O to a final volume of 50 μL.Cycling parameters were: 94 °C for 1 min; 30 cycles at 98 °C for 10 s, 60 °C for 15 s, 68 °C for 5.5 min; and 68 °C for 10 min.After sequencing the PCR products, the mutant piglets were all found to beBMPR-IBdisrupted, eitherBMPR-IB-/-orBMPR-IB-/746G, as described below.

    To detect whether the CRISPR plasmid vectors were integrated into the cloned piglet genomes, DNA of the cloned piglets was used as a template for PCR amplification with primer pairs nine and ten, respectively (Supplementary Table S1).The touchdown PCR conditions were 94 °C for 5 min; 26 cycles at 94 °C for 30 s, 68 °C (-0.5 °C/cycle) for 30 s, 72 °C for 45 s; 14 cycles at 94 °C for 30 s, 55 °C for 30 s, 72 °C for 45 s; and 72 °C for 10 min.To detect whether off-target mutations existed in theBMPR-IB-disrupted piglets, the top 15 predicted off-target sites (OTS) were selected using online software (https://www.benchling.com/).Amplicons were subjected to Sanger sequencing.The primers used are listed in Supplementary Table S2.

    Phenotype analysis

    OneBMPR-IB-/-, twoBMPR-IB-/746G, and two wild-type (WT)cloned piglets were used for skeletal phenotype analyses.We used the Digital Diagnost system (Philips, Netherlands) to take X-ray pictures of the whole-body skeletons of cloned piglets.The images were taken at 50 KV with 3 mA exposure.We further analyzed the anatomical structure of limbs in theBMPR-IB-disrupted and WT cloned piglets.Briefly, we used a scalpel, scissors, and tweezers to peel and remove the skin,muscle, and related adhesion tissues of the limbs.The limb skeletons of theBMPR-IB-/-(n=1),BMPR-IB-/746G(n=2), and WT (n=2) cloned piglets were then collected and stored at-80 °C.Both forelimbs and hindlimbs of each piglet were analyzed using the Micro-CT system (Nemo NMC-100,Pingseng Healthcare Inc., China) at a resolution of 50 μm,voltage of 90 kV, and current of 60 μA.Quantification of trabecular bone was assessed by measuring a 2 mm section starting 2 mm distal to the growth plate of the forelimb radius and hindlimb tibia.For cortical bone, 3 mm sections above the midshaft of the forelimb radius and hindlimb tibia were analyzed.Based on threshold segmentation and threedimensional (3D) measurements, quantitative analyses of trabecular and cortical bone (BV/TV, ratio of segmented bone volume to total volume; BS/TV, ratio of segmented bone surface to total volume; Tb.Sp, mean distance between trabeculae, assessed using direct 3D methods; Tb.N, average number of trabeculae per unit length, a key parameter for trabecular bone architecture; Tb.BMD, bone density measure,reflecting strength of bones as represented by calcium content; Ct.BMD, cortical bone mineral density; Ct.Th,average cortical thickness; Ct.ar/Tt.ar, cortical area fraction) in the radius and tibia were performed using Avatar v1.6.5(Pingseng Healthcare Inc., China).

    Quantitative reverse transcription PCR (qRT-PCR) and western blot analyses

    Total RNA was extracted from the liver, kidney, testicle,forelimb and hindlimb cartilage, forelimb ulna, and hindlimb fibula ofBMPR-IB-disrupted (BMPR-IB-/-,n=1;BMPR-IB-/746G,n=3) and WT (n=3) piglets using TRIzol Reagent (Invitrogen,USA).The mRNA was reverse transcribed into cDNA using a PrimeScriptTMRT Reagent Kit with gDNA Eraser (Takara,Japan).qRT-PCR was performed using the ABI 7 900 HT Fast Real-Time PCR System (Applied Biosystems, USA), and the thermal parameters were 50 °C for 2 min; 95 °C for 10 min; 40 cycles at 95 °C for 15 s, 60 °C for 50 s; 95 °C for 15 s; 60 °C for 15 s; and 95 °C for 15 s.The 2-??Ctformula was used to determine relative gene expression, which was normalized to the level ofGAPDHmRNA.Each reaction was performed in technical quadruplicate.Primer pairs eleven and twelve for qRT-PCR are listed in Supplementary Table S1.

    Total protein was extracted from the liver and kidney ofBMPR-IB-/-piglets (n=1) and forelimb ulna and hindlimb fibula ofBMPR-IB-/746G(n=3) and WT piglets (n=3) using a protein extraction kit (Applygen, China).Protein concentration was determined using the bicinchoninic acid (BCA) assay, and 30 μg of protein from each sample was resolved by 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDSPAGE) (Angle Gene, China) and transferred onto polyvinylidene fluoride (PVDF) membranes (Millipore, USA).After blocking with Quick Block? Blocking Buffer (Beyotime,China) for 1 h at room temperature, the membranes were washed in Tris Buffered Saline with Tween 20 (TBST) and incubated overnight at 4 °C with primary antibodies, including:rabbit polyclonal anti-BMPR-IB antibody (1:500; Sino Biological, China), mouse monoclonal antibodies for C1QA,SRSF1, P3H1, GJA1, TCOF1, RBM10, MMP13, and PHAX(1:100; Santa Cruz Biotechnology, USA), and β-actin (1:1 000;Abcam, UK).The membranes were washed and incubated with horseradish peroxidase (HRP)-conjugated secondary anti-rabbit or anti-mouse antibodies (1:2 000; Abcam, UK) for 1 h at room temperature.β-actin was used as the loading control.The bands were visualized with an ECL chemiluminescent kit (Beyotime, China) and scanned usingthe ChemiDocTMMP imaging system (Bio-Rad, USA).Band intensities were quantified using ImageJ (Rawak Software,Germany).The relative amounts of BMPR-IB, C1QA, SRSF1,P3H1, GJA1, TCOF1, RBM10, MMP13, and PHAX were calculated after normalization to β-actin.

    4D label-free quantitative proteomic analysis

    Proteomic analysis was conducted by PTM Bio Co., Ltd.(China).Hindlimb fibula samples from theBMPR-IB-/746G(n=3)and WT piglets (n=3) were ground into powder in liquid nitrogen and thoroughly homogenized in 8 mol/L urea containing 1% protease inhibitor cocktail.After insoluble debris was removed by centrifugation at 12 000gfor 10 min at 4 °C, the supernatant was collected and protein concentration was determined using a BCA protein assay kit (Beyotime,China).For digestion, the protein solution was reduced with 5 mmol/L dithiothreitol for 30 min at 56 °C and alkylated with 11 mmol/L iodoacetamide for 15 min at room temperature in the dark.The protein samples were then diluted with 100 mmol/L triethylammonium bicarbonate (TEAB) to a urea concentration of 2 mol/L.Finally, trypsin was added at a 1:50 trypsin-to-protein mass ratio for the first digestion overnight and a 1:100 trypsin-to-protein mass ratio for the second digestion of 4 h.

    The tryptic peptides were dissolved in 0.1% formic acid and 2% acetonitrile (solvent A) and loaded onto a home-made reversed-phase analytical column (15 cm long, 75 μm i.d.).The gradient was increased from 6% to 23% solvent B (0.1%formic acid in 90% acetonitrile) over 38 min, 23% to 35% in 14 min, and to 80% in 4 min, then held at 80% for the last 4 min,all at a constant flow rate of 550 nL/min on an EASY-nLC 1 000 ultra-performance liquid chromatography (UPLC) system(Thermo Fisher Scientific, USA).The peptides were subjected to nanospray ionization (NSI) source followed by tandem mass spectrometry (MS/MS) in a Q ExactiveTMPlus (Thermo Fisher Scientific, USA) coupled to the UPLC system.The electrospray voltage applied was 2.3 kV.The m/z scan range was 400 to 1 200 for full scan, and intact peptides were detected using Orbitrap at a resolution of 60 000.Peptides were then selected for MS/MS using normalized collision energy (NCE) setting as 28 and the fragments were detected using Orbitrap at a resolution of 17 500.Standard datadependent acquisition (DDA) procedures were implemented to detect and quantify peptides, including one MS acquisition at a mass/charge ratio (m/z) of 400-1 500, with the top 20 intense precursor ions subjected to MS/MS scans with 15.0 s dynamic exclusion.Automatic gain control (AGC) was set at 5×104and fixed first mass was set to 100 m/z.

    The resulting MS/MS data were processed using the MaxQuant search engine (v.1.5.2.8).Tandem mass spectra were searched against the Sus_scrofa_9823 database concatenated with the reverse decoy database.Trypsin/P was specified as the cleavage enzyme allowing up to four missing cleavages.Mass tolerance for the precursor ions was set to 20×10-6in the first search and 5×10-6in the main search, and mass tolerance for the fragment ions was set to 0.02 Da.Carbamidomethyl on Cys was specified as fixed modification and acetylation modification and oxidation on Met were specified as variable modifications.The false discovery rate(FDR) was adjusted to <1% and minimum score for modified peptides was set to >40.

    Functional enrichment analysis

    Proteins were classified by Gene Ontology (GO) annotation into three categories: biological process, cellular component,and molecular function.For each category, two-tailed Fisher’s exact test was employed to test differentially expressed protein (DEP) enrichment against all identified proteins.GO terms with a correctedP<0.05 were considered significant.For each DEP, the InterPro (http://www.ebi.ac.uk/interpro/)database was searched and two-tailed Fisher’s exact test was employed to test its (domain) enrichment against all identified proteins.Protein domains with a correctedP<0.05 were considered significant.

    Target proteomic analysis

    Parallel reaction monitoring (PRM) analysis was performed to verify and quantify the DEPs.One or two unique peptides of each target protein were selected from proteomic measurements based on thresholds of detection frequencies>50%, missed cleavage=0, andP<0.05.Finally, a spectral library of 10 selected proteins represented by 18 peptides was created.Proteins extracted from the hindlimb fibula ofBMPRIB-/746G(n=3) and WT piglets (n=3) were alkylated and digested, as described in “4D label-free quantitative proteomic analysis”.Protein (1 μg) was injected into the liquid chromatography-tandem mass spectrometry (LC-MS/MS)system for the PRM assay (2.1 kV electrospray voltage and 35 000 resolution (AGC target 3×106, maximum injection time 50 ms)).The isolation window for MS/MS was set at 2.0 m/z.Peptide settings were: enzyme, trypsin [KR/P]; peptide length,8-25; variable modification, carbamidomethyl on Cys and oxidation on Met; and max variable modifications, 3.Transition settings: precursor charges, 2 and 3; ion charges, 1 and 2; ion types, b, y, and p; product ions, ion 3 to last ion; and ion match tolerance, 0.02 Da.Skyline software (v3.6) was used for relative quantification for PRM study (Henderson et al.,2018).

    Statistical analysis

    All data are presented as mean±standard deviation (SD).Statistical analysis was performed using an independent twotailed Student’st-test.A null probability ofP<0.05 was considered statistically significant.SPSS v17.0 was used for all analyses.

    RESULTS

    CRISPR/Cas9-mediated BMPR-IB A746G mutation in PFFs

    To convert A746 of theBMPR-IBgene to 746G, we designed one sgRNA and 1 977 bp linear double-stranded donor(template) DNA for homology-directed repair (HDR)-mediated genome editing by CRISPR/Cas9 (Figure 1A).The linear double-stranded donor DNA, obtained from bridging PCR(Supplementary Figure S1A-C), encompassed the target point mutation (A746G) and sgRNA sequence, and was confirmed by Sanger sequencing (Supplementary Figure S1D).The donor DNA carried an altered PAM (GGG to GGC, Figure 1A;Supplementary Figure S1D).The Cas9/sgRNA expressionvector and template DNA were both co-transfected into an early passage of PFFs.A total of 113 single-cell clones were obtained after 48 h of puromycin selection and subsequent subculture.These clones were genotyped via short-range PCR and Sanger sequencing.Of these 113 clones, 39 contained the A746G mutation (16 homozygotes and 23 heterozygotes), 30 contained homozygous or heterozygous deletions, and 44 were WT (Supplementary Table S3).

    Generation of BMPR-IB-edited pigs

    We selected twoBMPR-IB/746GG and threeBMPR-IB/AA746 cell colonies identified by short-range PCRs for SCNT,obtaining approximately 3 500 embryos.All embryos were surgically transferred to 10 surrogate sows, six of which carried to term and gave birth to 16 cloned piglets (Figure 1E).The genotypes of the cloned piglets were determined by Sanger sequencing of long-range PCR products across the 1 977 bp donor sequence.As shown in Figure 1B-F, seven cloned piglets carried the expected 746G mutation at the target locus, but all had a 1 365 bp heterozygous deletion of theBMPR-IBgene (referred to as “BMPR-IB-/746G”).Two piglets had a 2 bp deletion on one allele and a 2 431 bp deletion on the other allele of theBMPR-IBgene (referred to as “BMPR-IB-/-”).Piglets carrying either of the two genotypes were considered asBMPR-IB-disrupted pigs.

    To detect whether the CRISPR plasmids were integrated into the host genome, we genotyped the Cas9 and sgRNA domains, and no integrations were found in any of the cloned piglets (Supplementary Figure S2).As off-target effects are a major concern when using the CRISPR/Cas9 system, we performed PCR amplification and sequencing for the top 15 potential OTS.No off-target effects were found in theBMPRIBmutant piglets (Supplementary Table S4).

    Phenotype characterization of BMPR-IB-disrupted pigs

    TheBMPR-IB-disrupted piglets were unable to stand normally,with their splayed limbs only able to move in a “paddling”motion (Video 1: https://figshare.com/s/966c57aee3cf5147 604c).Close examination revealed different malformed shapes of the forelimb and hindlimb.As shown in Figure 2A,the forelimbs ofBMPR-IB-disrupted piglets (BMPR-IB-/746GorBMPR-IB-/-) were severely distorted towards the abdomen,and the wrist joints could not be flexed or extended normally.In the hindlimbs ofBMPR-IB-disrupted piglets, the knee joints were moderately rigid and the ankles and toes were severely bent backwards.

    Radiographic and anatomical examinations showed marked skeletal hypoplasia of the limbs inBMPR-IB-disrupted piglets.In the forelimbs of theBMPR-IB-/746GandBMPR-IB-/-piglets,the ulna and radius were severely deformed and shortened.The carpal, metacarpal, and phalangeal bones ofBMPR-IB-disrupted piglets were disorderedly arranged and partially absent.Moreover, bone in the wrist joint was hardened,rendering the joint unable to flex (Figure 2B, C;Supplementary Figure S3).In the hindlimbs ofBMPR-IB-/746Gpiglets, the fibula was largely absent, tarsal bones were enlarged and partially absent, and phalangeal bones showed disordered arrangement.In addition, the ankle joint was swollen, with a large joint gap and collapsed third proximal phalanx (os compedale) (Figure 2B, D; Supplementary Figure S3).Of note, theBMPR-IB-/-piglets had nearly the same hindlimb deformity phenotypes as theBMPR-IB-/746Gpiglets,except for complete loss of the fibula inBMPR-IB-/-piglets(Figure 2B, D).

    Figure 2 Phenotypic characterization of BMPR-IB-disrupted piglets

    High-resolution micro-CT scans generated clearer visualization of limb bone morphology and skeletal arrangement, as well as quantitative bone morphometry(Figures 3, 4).Quantitative analysis of the forelimb radius(Figure 3C) revealed a 42% and 15% reduction in BV/TV in theBMPR-IB-/-andBMPR-IB-/746Gpiglets, respectively,compared with the WT piglets.Furthermore, trabecular BS/TV was 28% lower in theBMPR-IB-/-group than in the WT group,with no significant difference observed in theBMPR-IB-/746Ggroup.In addition, compared with the WT piglets, Tb.Sp was 24% and 61% larger in theBMPR-IB-/746GandBMPR-IB-/-groups, respectively.BMPR-IBdisruption did not appear to have a negative effect on Tb.BMD or Ct.BMD.In addition, no significant effects were observed in Ct.Th or Ct.ar/Tt.ar due to the large variations in WT piglets.Regarding the hindlimb tibia(Figure 4C), theBMPR-IB-/746GandBMPR-IB-/-groups had distinctly lower BV/TV (-26% and -35%, respectively)compared to the WT group.Trabecular BS/TV was 11%higher in theBMPR-IB-/-group compared to the WT, while no significant difference was found in theBMPR-IB-/746Ggroup.Furthermore, Tb.BMD was lower in both groups (-8% and-11%, respectively), but there was no significant difference in Ct.BMD.There was a 35% reduction in Ct.Th in theBMPR-IB-/-group, while the decrease in Ct.Th was not significant in theBMPR-IB-/746Ggroup (Figure 4C).In general, theBMPR-IB-/-group showed poorer morphometric deterioration than theBMPR-IB-/746Ggroup.

    Analysis of BMPR-IB expression in cloned piglets

    The qRT-PCR results showed thatBMPR-IBgene expression was lost in all examined tissues of theBMPR-IB-/-piglets.In theBMPR-IB-/746Gpiglets,BMPR-IBexpression was normal in the forelimb cartilage and ulna but was significantly elevated(P<0.05 orP<0.01) in the hindlimb cartilage and fibula, as well as the liver, kidney, and testes (Figure 5A).

    Western blot analysis confirmed the loss ofBMPR-IBin the liver and kidney ofBMPR-IB-/-piglets (Figure 5B, C).The western blots showed no significant difference in BMPR-IB protein expression in the forelimb ulna between theBMPR-IB-/746Gand WT piglets.However, BMPR-IB protein expression in the hindlimb fibula in theBMPR-IB-/746Gpiglets was significantly higher (P<0.05) than that in the WT piglets(Figure 5B, C).

    Figure 3 Micro-CT analyses of forelimb bones of piglets

    Figure 4 Micro-CT analyses of hindlimb bones of piglets

    Identification of DEPs

    Hindlimb fibula samples fromBMPR-IB-/746G(n=3) and WT piglets (n=3) were analyzed using the LC-MS/MS platform.A total of 6 412 proteins were identified and 5 957 proteins were quantified.Based on cut-off criteria (P<0.05, number of peptides>1, fold-change>1.5), a total of 139 proteins with one or more unique peptides were significantly differentially expressed in the hindlimb fibula between the two groups,including 51 down-regulated and 88 up-regulated in theBMPR-IB-/746Gpiglets (Figure 6A; Supplementary Table S5).

    GO and protein domain enrichment analysis

    To obtain a better mechanistic understanding of the protein networks that may be related to limb deformity formation inBMPR-IBdisrupted piglets, we used GO analysis to categorize DEPs.The DEPs in the hindlimb fibula were primarily involved in bone development regulation, such as regulation of bone remodeling, positive regulation of bone resorption, positive regulation of bone remodeling, regulation of bone mineralization, osteoblast differentiation, and regulation of biomineral tissue development (Figure 6B).

    Protein domain enrichment analysis was performed to identify functional domains of the DEPs.As shown in Figure 6C, secreted phosphoprotein 24 (Spp-24), high mobility group (HMG) box, chromo shadow domain, transferrin, and serpin (serine protease inhibitor) were significantly enriched in the hindlimb fibula.

    Analysis of target proteins using PRM

    To verify the changes in differential expression identified by 4D label-free quantitative proteomic analysis, PRM was employed using the same samples.Ten up-regulated proteins with potential involvement in bone development regulation were selected and a panel of seven proteins with one or two unique peptides were successfully detected and quantified.All seven proteins were significantly up-regulated in the hindlimb fibula ofBMPR-IB-/746Gpiglets (P<0.05,BMPR-IB-/746G/WT ratio>1) (Table 1).Thus, the PRM results were highly consistent with those from the 4D label-free quantitative phase(Table 1).

    Validation of expression levels of selected DEPs

    To validate the 4D label-free quantitative proteomic results,the expression levels of several DEPs in the hindlimb fibula were detected via western blotting.As shown in Figure 7, the protein expression levels of C1QA, SRSF1, TCOF1, RBM10,and MMP13 in the hindlimb fibula were significantly higher in theBMPR-IB-/746Gpiglets (P<0.01 orP<0.05) than in the WTpiglets, and GJA1 expression was higher in theBMPR-IB-/746Gpiglets than in the WT piglets, but with a weakly significant effect (P=0.096).PHAX protein expression in the hindlimb fibula was significantly lower inBMPR-IB-/746Gpiglets (P<0.01)than in WT piglets (Figure 7).Thus, the expression levels of the selected DEPs detected by western blotting were highly consistent with the 4D label-free quantitative proteomic data(Figure 7; Supplementary Table S5).

    Table 1 PRM-verified DEPs in hindlimb fibula of BMPR-IB-/746G piglets compared with WT individuals

    Figure 5 qRT-PCR and western blot analyses of BMPR-IB expression levels

    DISCUSSION

    In this study, we originally planned to produce cloned pigs with aBMPR-IB746GGmutation.We isolated cell clones carrying this mutation via CRISPR/Cas9 techniques using bridge PCR-amplified products as the HDR template.However, our primary design had two significant defects.First, we attempted to avoid re-cutting the repair template by changing GGG(PAM) to GGC (Figure 1A).However, we overlooked the fact that an extra NGG PAM was introduced by the A746G mutation (Figure 1A), enabling Cas9 to re-cleave the repair template.Second, we implemented short-range PCR(Supplementary Figure S1) rather than long-range PCR across the whole donor region to identify isolated cell clones.This led to unknown genome rearrangements in theBMPR-IBgene in the selected colonies (Figure 1F).These two factors resulted in the accidental production of piglets carrying a disruptedBMPR-IBgene, which, to the best of our knowledge,is the first report on the generation ofBMPR-IBgene disruption in pigs.Recent evidence indicates that complex rearrangements are frequently observed in CRISPR/Cas9 editing (Alanis-Lobato et al., 2021; Canaj et al., 2019;Skryabin et al., 2020).In this study, theBMPR-IB-/746Gpiglets carried a discontinuous 1 365 bp deletion with short fragments of intron 8 and exon 9 retained in the gene (Figure 1F).RegardingBMPR-IB-/-piglets, a 2 431 bp segment was deleted in one allele and a 2 bp segment was deleted in the other (Figure 1F).Both the 1 365 bp and 2 431 bp deletions spanned the donor region.Our study indicated that multiple PCR procedures may be a better option for detecting repair donor and flanking regions.This is because complex rearrangements can occur during genome editing, and a particular pair of primers may amplify only one allele of the target gene in cases where the other allele contains a deletion.

    In our study,BMPR-IB-/746Gpiglets harboring compoundmutations (g.746G and 1 365 bp del, Figure 1E) in theBMPRIBgene displayed phenotypes typically observed in limbs with skeletal dysplasia.The g.746A>G mutation in theBMPR-IBgene has been previously identified in domestic sheep and is suggested to increase ovulation rates and litter size (Chu et al., 2007; Mahdavi et al., 2014; Reader et al., 2012; Roy et al.,2011; Zhou et al., 2018).In the current study,BMPR-IBmRNA and protein expression levels in five of the seven examined tissues were significantly higher inBMPR-IB-/746Gpiglets than in WT piglets (P<0.01 orP<0.05; Figure 5), suggesting that the expression levels of single intactBMPR-IBallele (carrying the746Gvariant without deletion) inBMPR-IB-/746Gindividuals exceeded the sum of the transcript levels of the two alleles in the WT individuals.We speculate two possible reasons for this phenomenon.First, the A746G mutation may be located on acis-acting element that regulatesBMPR-IBexpression, and it may alter the function of this element, thereby affecting gene expression.Second, genetic compensation may be induced due to deleterious mutations (El-Brolosy & Stainier, 2017).Indeed, bioinformatic analysis in this study indicated that the A746G mutation in the highly conserved intracellular kinase signaling domain of the BMP-IB receptor is likely to damage its structure and function, and intriguingly, cause distinct alterations in secondary structures in pigs and sheep(Supplementary Figure S4).As the function of the mutant BMPR-IB protein was greatly reduced but not necessarilycompletely lost, theBMPR-IB-/746Gindividuals may exhibited a genetic compensation response, in which the intact mutant allele and its translated protein were overexpressed to compensate for its functional deficits.This may also explain whyBMPR-IB-/746Gpiglets are anatomically and pathologically similar toBMPR-IB-/-individuals but show less severe deformities in skeletal development (Figures 2D, 3, 4).

    Figure 6 DEPs in hindlimb fibula between BMPR-IB-/746G and WT groups

    Figure 7 Western blot analyses of DEP expression levels in hindlimb fibula between BMPR-IB-/746G and WT groups

    Disruption ofBMPR-IBin mice and humans often results in severe limb abnormalities.For example, in BMPR-IB-/-mice,the proximal and middle phalanges are reduced and fused,the radius, ulna, and tibia lengths are normal, but the metacarpals/metatarsals are shorter, and several carpal/tarsal bones are affected (Yi et al., 2000).In the current study,however, mutant piglets not only exhibited malformation in the carpal, metacarpal, and phalangeal bones, but the ulna and radius were also severely distorted and shortened (Figures 2B, C, 3A).In humans, both I200K and R486W missense mutations are reported to cause brachydactyly type A2(Lehmann et al., 2003).Furthermore, loss-of-function mutation(del359-366) can cause acromesomelic chondrodysplasia and genital anomalies, with hypoplasia of the femoral neck and head (Demirhan et al., 2005).The severe skeletal defects in ourBMPR-IBmutant (BMPR-IB-/746GorBMPR-IB-/-) piglets,such as loss of fibula, abnormally developed tarsal bones,disorderedly arranged phalangeal bones, and severely distorted and shortened ulna and radius, are analogous to many clinical phenotypes.We also found that the forelimb wrist joint inBMPR-IBmutant piglets was hardened and could not be flexed, the hindlimb ankle joint inBMPR-IB-/746Gpiglets was malformed, and the third proximal phalanx (os compedale) was collapsed (Figure 2C, D).However, no defects were found in other parts of the skeleton.

    As the hindlimb fibula was completely missing inBMPR-IB-/-piglets (Figures 2B, D, 4A), we obtained fibula samples fromBMPR-IB-/746Gand WT piglets for proteomic analysis.In total,139 DEPs (88 up-regulated, 51 down-regulated) were identified (Figure 6A; Supplementary Table S5).PRM analysis was performed to verify the quantitative data of seven randomly selected proteins (Table 1).Most DEPs were significantly involved in one or multiple functions of skeletal development, embryonic development, TGF-β pathway, and immune system, especially tumor progression (Supplementary Table S5).These results are consistent with the known functions of theBMPR-IBgene (Dituri et al., 2019; Rahman et al., 2015).A large proportion of genes in the top 50 DEPs were significantly associated with skeletal or embryonic development (Supplementary Table S5), includingC1QA,MYO1H,SRSF1,P3H1,GJA1,TCOF1,RBM10,SPP2,MMP13, andPHAX.C1QAis a subunit chain of C1q, which is implicated in osteoclast development from monocytes (Teo et al., 2012).Polymorphisms inMYO1Hare associated with mandibular retrognathism (Arun et al., 2016) and sagittal and vertical craniofacial skeletal patterns (Cunha et al., 2019).SRSF1is a prototype member of the serine/arginine (SR) rich family of splicing proteins and may modulate pattern formation(including cartilage formation) by inhibiting transcription of tissue-specific genes during embryogenesis (Lee et al., 2016).P3H1, also known asLEPRE1, forms a complex with cyclophilin B in the endoplasmic reticulum (Vranka et al.,2010).P3H1-null mice display abnormalities in fibrillarcollagen-rich tissues, including bones (Vranka et al., 2010),and mutations inP3H1can cause non-lethal (Takagi et al.,2012) and lethal recessive (Cabral et al., 2012) osteogenesis imperfections.Mutations inGJA1can lead to human oculodentodigital dysplasia, including face, eye, tooth, and limb deformities (Paznekas et al., 2009; Sargiannidou et al.,2021).TCOF1encodes the nucleolar phosphoprotein treacle,and its mutation is responsible for Treacher Collins syndrome(Dai et al., 2016; Grzanka & Piekie?ko-Witkowska, 2021).Mutations inRBM10can lead to splicing changes that affect mouse palate development (Rodor et al., 2017).SPP2is a highly phosphorylated and glycosylated sialoprotein and a prominent component of mineralized extracellular matrices in bone.Full-length SPP2 (24 kD) inhibits BMP-induced bone formation, while its degradation product (18.5 kD) (designated spp18.5) appears to be pro-osteogenic (Brochmann et al.,2009; Sintuu et al., 2008).We speculate that the biological function of SPP2 may be finely regulated by proteolysis, and this process may be altered inBMPR-IB-disrupted pigs.MMP13is considered to play a crucial role in bone formation and remodeling and is expressed in both terminal hypertrophic chondrocytes in the growth plate and in osteoblasts (Li et al.,2017; Stickens et al., 2004).MMP13is up-regulated at the onset of osteoarthritis (Li et al., 2017), andMMP13-deficient mice show abnormal development of the skeletal growth plate(Stickens et al., 2004).A research suggests thatPHAXis associated with Pierre Robin syndrome, which is characterized by congenital micrognathia, glossoptosis and airway obstruction (Ansari et al., 2014).Here, GO analysis showed noticeable enrichment in biological pathways related to bone and embryonic development, including regulation of bone remodeling, positive regulation of bone resorption,positive regulation of bone remodeling, regulation of bone mineralization, osteoblast differentiation, and regulation of biomineral tissue development.

    In summary, theBMPR-IB-disrupted piglets obtained via CRISPR-Cas9 and SCNT exhibited walking difficulties and severe developmental deformities in the forelimbs and hindlimbs.TheBMPR-IB-/746GandBMPR-IB-/-piglets showed similar phenotypes.We identified 139 DEPs in the hindlimb fibula ofBMPR-IB-/746Gpiglets with limb deformities, most of which are involved in skeletal and (or) embryonic development.Our study provides novel insights intoBMPR-IBdeficiency in a large mammal.

    SUPPLEMENTARY DATA

    Supplementary data to this article can be found online.

    COMPETING INTERESTS

    The authors declare that they have no competing interests.

    AUTHORS’ CONTRIBUTIONS

    Conceptualization, Y.Y.X.and Q.Y.; formal analysis and experiment, Q.Y.and C.M.Q.; investigation, W.W.L.and H.Y.J.; methodology, Q.Q.J.and Y.Y.L.; resources, Y.Y.X.and J.R.; visualization, Q.Y.and C.M.Q; writing-original draft, Q.Y.;writing-review & editing, Y.Y.X.All authors read and approved the final version of the manuscript.

    ACKNOWLEDGEMENTS

    We would like to thank Prof.Lu-Sheng Huang at Jiangxi Agricultural University for kind suggestions during the experimental design process.

    女人爽到高潮嗷嗷叫在线视频| 午夜福利欧美成人| 亚洲欧美精品综合一区二区三区| 中文字幕av电影在线播放| 亚洲成人手机| 在线十欧美十亚洲十日本专区| 在线观看免费视频日本深夜| 少妇精品久久久久久久| 日本精品一区二区三区蜜桃| 视频区欧美日本亚洲| 香蕉丝袜av| 精品福利永久在线观看| 国产精品一区二区在线观看99| 国产亚洲欧美在线一区二区| 精品国产乱码久久久久久男人| 69av精品久久久久久 | 午夜福利视频精品| 精品久久久精品久久久| 又黄又粗又硬又大视频| 精品久久久精品久久久| 成人黄色视频免费在线看| 男女床上黄色一级片免费看| 色在线成人网| 人人妻,人人澡人人爽秒播| 国产日韩欧美在线精品| 欧美精品高潮呻吟av久久| 久久亚洲精品不卡| 五月开心婷婷网| 免费看十八禁软件| 成人国产av品久久久| 少妇猛男粗大的猛烈进出视频| 老司机亚洲免费影院| 侵犯人妻中文字幕一二三四区| 亚洲成人免费电影在线观看| 妹子高潮喷水视频| 亚洲三区欧美一区| 在线观看免费视频日本深夜| 777久久人妻少妇嫩草av网站| 一区二区三区激情视频| 久久久精品区二区三区| 国产一区二区三区视频了| 老熟妇仑乱视频hdxx| 大型黄色视频在线免费观看| 精品一区二区三卡| 免费在线观看完整版高清| 黑丝袜美女国产一区| 中文欧美无线码| 99精品久久久久人妻精品| 狂野欧美激情性xxxx| 大陆偷拍与自拍| 中文字幕人妻丝袜一区二区| 国产av一区二区精品久久| 国产精品香港三级国产av潘金莲| 老熟妇仑乱视频hdxx| 美女午夜性视频免费| 久久久精品免费免费高清| 久久精品熟女亚洲av麻豆精品| 国产免费现黄频在线看| 亚洲精华国产精华精| 午夜精品久久久久久毛片777| 免费在线观看黄色视频的| 精品午夜福利视频在线观看一区 | 亚洲成人免费电影在线观看| 亚洲专区字幕在线| 成人av一区二区三区在线看| 精品亚洲乱码少妇综合久久| 久久天堂一区二区三区四区| 日本黄色日本黄色录像| 亚洲精品中文字幕一二三四区 | 日本一区二区免费在线视频| 欧美老熟妇乱子伦牲交| 亚洲精品粉嫩美女一区| 亚洲全国av大片| 精品福利观看| 丝瓜视频免费看黄片| 国产精品久久久人人做人人爽| 大陆偷拍与自拍| 日韩有码中文字幕| 久久人人爽av亚洲精品天堂| 黑人巨大精品欧美一区二区mp4| 亚洲av成人不卡在线观看播放网| 国产不卡一卡二| 国产欧美亚洲国产| 国产一区二区三区在线臀色熟女 | 成人18禁高潮啪啪吃奶动态图| 国产av国产精品国产| 亚洲精品成人av观看孕妇| 在线观看人妻少妇| 精品人妻在线不人妻| 国产精品一区二区精品视频观看| 午夜日韩欧美国产| 欧美成人午夜精品| 午夜福利欧美成人| 他把我摸到了高潮在线观看 | 午夜久久久在线观看| 久久久久久久久免费视频了| netflix在线观看网站| 一级a爱视频在线免费观看| 欧美+亚洲+日韩+国产| 男人舔女人的私密视频| 国产免费av片在线观看野外av| 少妇裸体淫交视频免费看高清 | 免费观看av网站的网址| 国产成人一区二区三区免费视频网站| 在线播放国产精品三级| 高清欧美精品videossex| 精品久久久精品久久久| 在线观看免费视频日本深夜| 搡老岳熟女国产| videosex国产| 欧美亚洲 丝袜 人妻 在线| 夜夜夜夜夜久久久久| 国产日韩一区二区三区精品不卡| 一本久久精品| 一边摸一边抽搐一进一出视频| av线在线观看网站| a级片在线免费高清观看视频| 亚洲黑人精品在线| 在线观看免费午夜福利视频| 一级毛片电影观看| 黑人操中国人逼视频| 悠悠久久av| 99国产精品免费福利视频| 亚洲第一青青草原| 亚洲国产精品一区二区三区在线| 一区二区av电影网| 最近最新中文字幕大全电影3 | 香蕉国产在线看| 成人18禁高潮啪啪吃奶动态图| 五月天丁香电影| 免费黄频网站在线观看国产| 黄频高清免费视频| 久9热在线精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品亚洲精品国产色婷小说| 午夜福利视频在线观看免费| av视频免费观看在线观看| 亚洲一码二码三码区别大吗| 在线av久久热| 搡老熟女国产l中国老女人| 国产精品久久久久久精品古装| 日本五十路高清| 国产一卡二卡三卡精品| 纵有疾风起免费观看全集完整版| 亚洲精品国产一区二区精华液| 啦啦啦视频在线资源免费观看| 精品国产亚洲在线| avwww免费| 国产成人精品在线电影| 69精品国产乱码久久久| 91老司机精品| 下体分泌物呈黄色| 中文字幕制服av| 91成年电影在线观看| 久久免费观看电影| 成人三级做爰电影| 精品高清国产在线一区| 超碰97精品在线观看| 美女高潮喷水抽搐中文字幕| 日本av免费视频播放| 黄色丝袜av网址大全| 国产精品99久久99久久久不卡| 丁香六月天网| 99在线人妻在线中文字幕 | 最黄视频免费看| 精品熟女少妇八av免费久了| 免费一级毛片在线播放高清视频 | 国产在线视频一区二区| 精品熟女少妇八av免费久了| 国精品久久久久久国模美| 日韩免费av在线播放| 亚洲成人国产一区在线观看| 亚洲精品中文字幕一二三四区 | 成人18禁高潮啪啪吃奶动态图| 美女主播在线视频| 国产熟女午夜一区二区三区| 久久中文字幕人妻熟女| tocl精华| 热re99久久国产66热| 侵犯人妻中文字幕一二三四区| 国产成人系列免费观看| 1024香蕉在线观看| 黄频高清免费视频| 色婷婷av一区二区三区视频| 人人妻人人添人人爽欧美一区卜| 天堂俺去俺来也www色官网| 午夜免费成人在线视频| 国产精品香港三级国产av潘金莲| 国产亚洲精品久久久久5区| 三级毛片av免费| 最新美女视频免费是黄的| 国产精品98久久久久久宅男小说| tube8黄色片| 美女高潮喷水抽搐中文字幕| 一区二区三区精品91| 久久午夜综合久久蜜桃| 一本久久精品| 王馨瑶露胸无遮挡在线观看| 91精品国产国语对白视频| av线在线观看网站| 一二三四社区在线视频社区8| 国产一区二区三区在线臀色熟女 | 亚洲中文日韩欧美视频| 999久久久精品免费观看国产| 日本欧美视频一区| 啦啦啦 在线观看视频| 黄片大片在线免费观看| 国产成人精品久久二区二区免费| 日韩精品免费视频一区二区三区| 啪啪无遮挡十八禁网站| 国产片内射在线| 热re99久久精品国产66热6| 午夜福利免费观看在线| 99热国产这里只有精品6| 亚洲国产中文字幕在线视频| 一区二区三区乱码不卡18| 无遮挡黄片免费观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲伊人久久精品综合| 一区二区三区国产精品乱码| 国产日韩欧美视频二区| 亚洲,欧美精品.| 久久久精品免费免费高清| 国产成人精品无人区| av又黄又爽大尺度在线免费看| 99九九在线精品视频| 三上悠亚av全集在线观看| 91麻豆精品激情在线观看国产 | 国产视频一区二区在线看| 精品一区二区三卡| 亚洲七黄色美女视频| 国产精品偷伦视频观看了| 中文字幕最新亚洲高清| 老司机在亚洲福利影院| 自线自在国产av| av有码第一页| 免费人妻精品一区二区三区视频| 久久久国产一区二区| 女性被躁到高潮视频| 午夜福利乱码中文字幕| a在线观看视频网站| 高清在线国产一区| 下体分泌物呈黄色| 久久毛片免费看一区二区三区| 日韩人妻精品一区2区三区| 欧美日韩福利视频一区二区| 免费在线观看日本一区| 久久久久久人人人人人| 人人妻人人爽人人添夜夜欢视频| 欧美日韩成人在线一区二区| 伦理电影免费视频| 日韩一区二区三区影片| a级片在线免费高清观看视频| 最近最新免费中文字幕在线| 大片电影免费在线观看免费| 老司机亚洲免费影院| 精品人妻1区二区| 最近最新中文字幕大全免费视频| 成人三级做爰电影| 国产精品98久久久久久宅男小说| 午夜福利欧美成人| 日韩人妻精品一区2区三区| av福利片在线| 亚洲色图综合在线观看| 免费在线观看影片大全网站| 在线天堂中文资源库| 两性夫妻黄色片| 夫妻午夜视频| 在线观看免费视频日本深夜| 午夜福利免费观看在线| 中文字幕色久视频| 国产精品久久久久成人av| 国产真人三级小视频在线观看| 免费看十八禁软件| 淫妇啪啪啪对白视频| 三上悠亚av全集在线观看| 欧美激情高清一区二区三区| 法律面前人人平等表现在哪些方面| 另类精品久久| 国产男女超爽视频在线观看| 首页视频小说图片口味搜索| 9热在线视频观看99| 国产欧美日韩精品亚洲av| 90打野战视频偷拍视频| 在线观看免费视频网站a站| 国产成人影院久久av| 午夜久久久在线观看| 在线观看www视频免费| 91av网站免费观看| 可以免费在线观看a视频的电影网站| 欧美精品亚洲一区二区| 日韩欧美三级三区| 亚洲伊人久久精品综合| 美女国产高潮福利片在线看| av视频免费观看在线观看| 国产欧美日韩一区二区精品| 亚洲av第一区精品v没综合| 美女主播在线视频| 午夜日韩欧美国产| 丁香六月天网| 这个男人来自地球电影免费观看| 日日爽夜夜爽网站| 精品视频人人做人人爽| 国产在线一区二区三区精| 黄片小视频在线播放| 91国产中文字幕| 久久青草综合色| 丁香六月天网| 母亲3免费完整高清在线观看| 高清av免费在线| 国产亚洲一区二区精品| 国产在线视频一区二区| 日本wwww免费看| 老司机深夜福利视频在线观看| 亚洲第一av免费看| 久久人妻福利社区极品人妻图片| 黑人欧美特级aaaaaa片| 久久久久久久大尺度免费视频| 色94色欧美一区二区| 成人国语在线视频| 亚洲中文字幕日韩| 在线永久观看黄色视频| 欧美日韩成人在线一区二区| 真人做人爱边吃奶动态| 国产男女内射视频| 在线观看舔阴道视频| 久久精品国产a三级三级三级| 亚洲av日韩在线播放| 欧美久久黑人一区二区| 亚洲国产成人一精品久久久| 成在线人永久免费视频| 精品欧美一区二区三区在线| 成人国产一区最新在线观看| 99精国产麻豆久久婷婷| 黄色丝袜av网址大全| 少妇粗大呻吟视频| 国产精品久久久久久人妻精品电影 | 国产成人av教育| 午夜91福利影院| 国产又爽黄色视频| 欧美成狂野欧美在线观看| 真人做人爱边吃奶动态| 国产成人系列免费观看| 最新的欧美精品一区二区| 色播在线永久视频| 国产xxxxx性猛交| 91国产中文字幕| 无限看片的www在线观看| 他把我摸到了高潮在线观看 | 亚洲国产成人一精品久久久| 久久人妻熟女aⅴ| 成人18禁高潮啪啪吃奶动态图| 成人国产一区最新在线观看| 人成视频在线观看免费观看| 青青草视频在线视频观看| 欧美日韩亚洲高清精品| 精品国产乱码久久久久久男人| 欧美日韩亚洲国产一区二区在线观看 | 纵有疾风起免费观看全集完整版| 制服人妻中文乱码| 国产在线一区二区三区精| 在线观看免费日韩欧美大片| 欧美日韩av久久| 天天添夜夜摸| 夫妻午夜视频| 亚洲少妇的诱惑av| 男人操女人黄网站| 国产精品98久久久久久宅男小说| 视频区欧美日本亚洲| 两个人看的免费小视频| 日韩视频在线欧美| 一区在线观看完整版| 国产精品熟女久久久久浪| avwww免费| 18禁黄网站禁片午夜丰满| 777久久人妻少妇嫩草av网站| 免费高清在线观看日韩| 天堂8中文在线网| 国产真人三级小视频在线观看| 无人区码免费观看不卡 | 中亚洲国语对白在线视频| 亚洲国产欧美网| 精品少妇黑人巨大在线播放| 国产高清videossex| 国产在线观看jvid| 国产麻豆69| 在线观看免费视频日本深夜| 久久精品国产亚洲av高清一级| 免费一级毛片在线播放高清视频 | 在线观看免费午夜福利视频| 侵犯人妻中文字幕一二三四区| 亚洲欧美一区二区三区黑人| 天天躁日日躁夜夜躁夜夜| 少妇被粗大的猛进出69影院| 久久久精品国产亚洲av高清涩受| 欧美精品亚洲一区二区| 老司机影院毛片| 啦啦啦在线免费观看视频4| 美女高潮喷水抽搐中文字幕| 少妇的丰满在线观看| 好男人电影高清在线观看| 黄片大片在线免费观看| 免费在线观看视频国产中文字幕亚洲| 啦啦啦 在线观看视频| 国产成人欧美在线观看 | 亚洲中文字幕日韩| 91av网站免费观看| av国产精品久久久久影院| 精品一区二区三卡| 99香蕉大伊视频| 下体分泌物呈黄色| 午夜福利免费观看在线| 制服人妻中文乱码| 国产又爽黄色视频| 日本撒尿小便嘘嘘汇集6| 999精品在线视频| 亚洲黑人精品在线| 老司机亚洲免费影院| 人妻一区二区av| 国产激情久久老熟女| 韩国精品一区二区三区| 国产精品 国内视频| 又黄又粗又硬又大视频| 深夜精品福利| 老司机福利观看| 日韩欧美国产一区二区入口| a级毛片在线看网站| 一本综合久久免费| 日本精品一区二区三区蜜桃| 成人三级做爰电影| 老熟女久久久| 桃红色精品国产亚洲av| 欧美性长视频在线观看| 国产一区二区三区综合在线观看| 国产成人啪精品午夜网站| 精品亚洲成国产av| 欧美日韩一级在线毛片| 国产精品免费视频内射| 欧美日韩亚洲综合一区二区三区_| 高清视频免费观看一区二区| 国产一区二区三区在线臀色熟女 | 国产91精品成人一区二区三区 | 精品一区二区三卡| av福利片在线| 一级黄色大片毛片| 亚洲精品中文字幕一二三四区 | 色94色欧美一区二区| 精品第一国产精品| 日韩精品免费视频一区二区三区| 国产成人一区二区三区免费视频网站| 午夜两性在线视频| 51午夜福利影视在线观看| 99久久99久久久精品蜜桃| 国产欧美亚洲国产| 亚洲情色 制服丝袜| 亚洲第一av免费看| 国产精品秋霞免费鲁丝片| 国产一区二区三区视频了| 这个男人来自地球电影免费观看| 激情在线观看视频在线高清 | 在线看a的网站| 在线 av 中文字幕| 国产三级黄色录像| 午夜福利视频精品| 纵有疾风起免费观看全集完整版| 久久久久久免费高清国产稀缺| 波多野结衣一区麻豆| 国产男靠女视频免费网站| 国产精品久久电影中文字幕 | 香蕉丝袜av| 国产精品一区二区精品视频观看| 麻豆国产av国片精品| www.自偷自拍.com| 欧美另类亚洲清纯唯美| 人人妻人人爽人人添夜夜欢视频| 久久精品国产综合久久久| cao死你这个sao货| 黄色 视频免费看| 亚洲成a人片在线一区二区| 国产高清激情床上av| 人人妻人人添人人爽欧美一区卜| 欧美日韩亚洲综合一区二区三区_| 色老头精品视频在线观看| 久久性视频一级片| videos熟女内射| 男女高潮啪啪啪动态图| av有码第一页| 9热在线视频观看99| 成人黄色视频免费在线看| 精品亚洲成国产av| 亚洲国产中文字幕在线视频| 黄色视频不卡| 欧美乱码精品一区二区三区| 久久精品91无色码中文字幕| 亚洲精品中文字幕在线视频| 成人国语在线视频| 免费在线观看完整版高清| 亚洲午夜理论影院| 国产免费福利视频在线观看| 欧美日韩黄片免| 性少妇av在线| 热99久久久久精品小说推荐| 国产人伦9x9x在线观看| 热99久久久久精品小说推荐| 欧美大码av| 国产又色又爽无遮挡免费看| 女人久久www免费人成看片| 电影成人av| 亚洲精品在线美女| 视频区图区小说| 在线观看人妻少妇| 大码成人一级视频| 久久久欧美国产精品| 黄片播放在线免费| 亚洲精品粉嫩美女一区| 久久 成人 亚洲| 国产精品av久久久久免费| 色老头精品视频在线观看| 视频在线观看一区二区三区| 建设人人有责人人尽责人人享有的| 两性夫妻黄色片| 激情视频va一区二区三区| 桃红色精品国产亚洲av| 亚洲色图av天堂| av网站免费在线观看视频| 亚洲一区中文字幕在线| 亚洲av欧美aⅴ国产| 精品午夜福利视频在线观看一区 | 视频区图区小说| 侵犯人妻中文字幕一二三四区| 欧美激情高清一区二区三区| 水蜜桃什么品种好| 啪啪无遮挡十八禁网站| 9热在线视频观看99| 99riav亚洲国产免费| 激情视频va一区二区三区| 在线观看免费午夜福利视频| 好男人电影高清在线观看| 欧美成狂野欧美在线观看| 亚洲第一欧美日韩一区二区三区 | 国产精品成人在线| 另类精品久久| 久9热在线精品视频| 国产日韩欧美亚洲二区| 99热国产这里只有精品6| 日韩制服丝袜自拍偷拍| 国产极品粉嫩免费观看在线| 好男人电影高清在线观看| 中文字幕人妻熟女乱码| 亚洲精品在线观看二区| 亚洲国产av影院在线观看| 大型黄色视频在线免费观看| 无人区码免费观看不卡 | 如日韩欧美国产精品一区二区三区| www日本在线高清视频| 亚洲黑人精品在线| 亚洲三区欧美一区| av福利片在线| 十八禁高潮呻吟视频| 国产区一区二久久| 天天影视国产精品| 丰满人妻熟妇乱又伦精品不卡| 亚洲av美国av| 日本av手机在线免费观看| 成人免费观看视频高清| 女性被躁到高潮视频| 久久九九热精品免费| 亚洲精品中文字幕在线视频| bbb黄色大片| tube8黄色片| 午夜久久久在线观看| 日本精品一区二区三区蜜桃| 日日夜夜操网爽| 啦啦啦 在线观看视频| 国产免费视频播放在线视频| 国产精品一区二区在线观看99| 亚洲av成人一区二区三| 国产av精品麻豆| 欧美性长视频在线观看| 1024视频免费在线观看| 久久久久久免费高清国产稀缺| bbb黄色大片| 丁香欧美五月| 久久人人97超碰香蕉20202| 国产成人免费观看mmmm| 老汉色av国产亚洲站长工具| 九色亚洲精品在线播放| 亚洲免费av在线视频| 久久久久精品人妻al黑| 国产精品1区2区在线观看. | 免费久久久久久久精品成人欧美视频| 91成年电影在线观看| 国产亚洲一区二区精品| 777米奇影视久久| 91成年电影在线观看| 久久久精品区二区三区| 欧美另类亚洲清纯唯美| 天天添夜夜摸| 丁香六月欧美| 亚洲av日韩在线播放| 国产精品久久久久久精品古装| 午夜久久久在线观看| 国产精品av久久久久免费| 久久婷婷成人综合色麻豆| 免费av中文字幕在线| 少妇裸体淫交视频免费看高清 | 男男h啪啪无遮挡| 亚洲五月婷婷丁香| 国产精品香港三级国产av潘金莲| netflix在线观看网站| 后天国语完整版免费观看| 午夜福利一区二区在线看| 国产精品久久久久成人av| 人人妻人人澡人人看| 国产伦理片在线播放av一区|