Biallelic variants in MAD2L1BP (p31comet) cause female infertility characterized by oocyte maturation arrest

  1. Lingli Huang
  2. Wenqing Li
  3. Xingxing Dai
  4. Shuai Zhao
  5. Bo Xu
  6. Fengsong Wang
  7. Ren-Tao Jin
  8. Lihua Luo
  9. Limin Wu
  10. Xue Jiang
  11. Yu Cheng
  12. Jiaqi Zou
  13. Caoling Xu
  14. Xianhong Tong
  15. Heng-Yu Fan  Is a corresponding author
  16. Han Zhao  Is a corresponding author
  17. Jianqiang Bao  Is a corresponding author
  1. Reproductive and Genetic Hospital, The First Affiliated Hospital of USTC, Division of Life Sciences and Medicine, University of Science and Technology of China, China
  2. Center for Reproductive Medicine, Department of Obstetrics and Gynecology, Anhui Provincial Hospital Affiliated to Anhui Medical University, China
  3. Hefei National Laboratory for Physical Sciences at Microscale, Biomedical Sciences and Health Laboratory of Anhui Province, University of Science and Technology of China (USTC), China
  4. Life Sciences Institute, Zhejiang University, China
  5. International Institutes of Medicine, the Fourth Affiliated Hospital of Zhejiang University School of Medicine, China
  6. Hospital for Reproductive Medicine, State Key Laboratory of Reproductive Medicine and Offspring Health, Key Laboratory of Reproductive Endocrinology of Ministry of Education, Shandong Key Laboratory of Reproductive Medicine, Shandong Provincial Clinical Research Center for Reproductive Health, Shandong University, China
  7. School of Life Science, Anhui Medical University, China
4 figures, 4 tables and 2 additional files

Figures

Figure 1 with 2 supplements
Identification of pathogenic variants in MAD2L1BP from three unrelated families.

(A) Pedigrees of three female patients with oocyte maturation arrest at the MI stage from three unrelated families. The first case was from a consanguineous family, while the other two were not. The patients carried homozygous c.853C>T [p.R285*] and c.541C>T [p.R181*] variants from Family 1 and Family 3, respectively, while the patient from Family 2 had a compound heterozygous mutation c.518delT and c.21–94G>A as indicated. The chromatograms below show the Sanger sequencing results of the PCR-amplified fragments containing the respective variants in each family. (B) Phylogenetic conservation of the identified amino acid mutations in MAD2L1BP. The positions of all mutations are indicated in the schematic genomic structure of MAD2L1BP. The MAD2-interacting site and the seatbelt self-interaction region for MAD2L1BP are labeled in the secondary structure of MAD2L1BP protein at the bottom. All identified mutations residing in Exon 4 of MAD2L1BP generated the premature STOP codon, likely resulting in truncated proteins. (C) The representative morphology of normal and affected individual oocytes. The normal MI and MII oocytes are shown in the top panel. The red arrow points to the first polar body (PB1). Oocytes retrieved from Family 1 (II-1) and Family 2 (II-1) were arrested at the MI stage. The asterisks indicate the oocytes with a magnified view in the right panel. MI: metaphase I; MII: metaphase II. Scale bar, 80 μm.

Figure 1—figure supplement 1
Multi-sequence alignment of amino acids for MAD2L1BP (p31comet) (A) and its interacting partner MAD2 (B) orthologs from eutherian species with Jalview.

Pathogenic mutations in affected individuals reported in this study are marked in red box.

Figure 1—figure supplement 2
Semen analysis of the individual II-2 in Family 1 with a homozygous mutation (c.853C>T [p.R285*]) in MAD2L1BP.

(A–C) Photomicrographs of the ejaculated semen smear stained by papanicolaou. Few normally shaped spermatozoa can be observed in the ejaculated semen from the individual II-2 in Family 1. (D) The control sperm morphology by papanicolaou staining from a fertile man. Scale bars, 10 μm. See Supplementary file 1A for semen parameters of the individual II-2 in Family 1.

Figure 2 with 2 supplements
Adverse impacts of MAD2L1BP variants on mRNA expression and protein function.

(A–B) Quantitative PCR (qPCR) assays comparing the mRNA expression levels of MAD2L1BP in peripheral blood samples from Family 1 and Family 2 as indicated. The experiments were performed in technical triplicates. (C) qPCR assay showing the mRNA expression levels of MAD2 in individual blood samples from Family 1. The experiments were performed in technical triplicates. (D) The ectopic protein expression levels of Myc-MAD2, Flag-MAD2L1BP WT(wild-type), and Flag-MAD2L1BP MUT (p.R285*) by immunoblotting following transfection of individual plasmids into 293T cells in culture. Full-length and truncated MAD2L1BP protein can be seen in the panel. (E) Co-immunoprecipitation (Co-IP) assays demonstrating that the truncated MAD2L1BP MUT (p.R285*) lost interaction with MAD2, as compared with MAD2L1BP WT, when co-transfected into 293T cells in vitro. (F) Ribbon and topological diagrams showing the structure of MAD2L1BP WT (PDB ID: 2QYF) and the predicted MAD2L1BP MUT (p.R285*). Notably, MAD2L1BP MUT (p.R285*) lacks a C-terminal seatbelt configuration (highlighted in red) resulting in structural instability.

Figure 2—source data 1

Blots for detecting Myc-MAD2, Flag-MAD2L1BP WT, Flag-MAD2L1BP MUT protein expression.

https://cdn.elifesciences.org/articles/85649/elife-85649-fig2-data1-v2.zip
Figure 2—source data 2

Co-immunoprecipitation (Co-IP) blots for detecting interaction between MAD2 and MAD2L1BP WT or MAD2L1BP MUT.

https://cdn.elifesciences.org/articles/85649/elife-85649-fig2-data2-v2.zip
Figure 2—source data 3

The original file of blots for Figure 2D.

https://cdn.elifesciences.org/articles/85649/elife-85649-fig2-data3-v2.zip
Figure 2—source data 4

The original file of blots for Figure 2E.

https://cdn.elifesciences.org/articles/85649/elife-85649-fig2-data4-v2.zip
Figure 2—figure supplement 1
Phylogenetic tree of MAD2L1BP and mRNA expression pattern of MAD2L1BP and MAD2 in human follicles and mouse oocytes.

(A) Evolution of the MAD2L1BP gene orthologs in metazoan species. Sequence conservation analysis was performed with Jalview software, and processed on Evolview website. (B–C) mRNA expression levels of MAD2L1BP (B) and MAD2L1(MAD2) (C) in developing human follicles and mouse oocytes. RNA-seq data was downloaded from GSE71434 and GSE107746. Reads were mapped to the mouse genome (mm10) or human genome (hg38) using STAR. Relative expression quantification of genes and transcripts was performed using RSEM. PN, postnatal; GV, germinal vesicle; MII, metaphase II; rep, replicate. O, oocyte; G, granulosa cell.

Figure 2—figure supplement 2
Minigene assay showing the aberrant alternative splicing of MAD2L1BP exons associated with the variant c.21-94G>A from the Family 2.

(A) Schematic diagram of minigene splicing assay. The G→A substitution in intron 1 produced a novel acceptor site for Exon 2, leading to the generation of a premature termination codon (PTC). HEK293T cells were transfected with minigene constructs (WT or c.21-94G>A), RT-PCR was performed with the primers as indicated. (B) Different splicing patterns in WT and c.21-94G>A group were shown by the agarose gel. (C) Representative Sanger sequencing chromatogram validating the excessively splicing of the exon 2 in MAD2L1BP mRNA with c.21-94G>A and thus caused deletion of large fragment of MAD2L1BP protein.

Figure 3 with 1 supplement
cRNA microinjection of Mad2l1bp accelerated or rescued the meiotic division in mouse germinal vesicle (GV) oocytes.

(A) A schematic diagram depicting the time points of distinctive events through meiotic progression in mouse oocytes. cRNAs were microinjected into GV oocytes after 44 hr of PMSG priming. (B) Representative morphology of oocytes after microinjection with cRNAs encoding mock GFP, mouse Mad2l1bp WT, and Mad2l1bp Mut (equivalent to human p.R285*). Non-injected group was treated similarly to the other three groups without cRNA microinjection. Images were taken 14 hr following release from milrinone inhibitor. Scale bar, 100 μm. (C and D) Kinetic recordings showing the percentages of oocytes with GV breakdown (GVBD) (C), and the oocyte polar body extrusion (PBE) rate (D) through the time-lapse imaging experiment. The total numbers of oocytes microinjected were labeled as indicated. The experiments were performed in biological triplicates. Data are presented as mean ± SEM. (E) Representative morphology of oocytes after cRNA co-microinjection of Mad2 supplemented with WT or Mad2l1bp Mut (equivalent to human p.R285*) into the mouse GV oocytes. Images were taken 14 hr after release from milrinone. Scale bar, 100 μm. (F) Comparison of the percentages of oocytes with GVBD and PBE among the four microinjection groups in (E). The total numbers of oocytes microinjected were labeled as indicated. The experiments were performed in biological triplicates. The statistics were performed with Student’s t-test. Data are presented as mean ± SEM. ‘***‘ indicates p<0.001.

Figure 3—source data 1

cRNA microinjection of full-length or truncated Mad2l1bp uncovered their discordant roles in driving the extrusion of polar body 1 (PB1) in mouse oocytes.

https://cdn.elifesciences.org/articles/85649/elife-85649-fig3-data1-v2.zip
Figure 3—source data 2

Co-injection of Mad2 cRNAs with Mad2l1bp cRNAs (WT vs p.R285* mutant) indicated the Mad2l1bp variant lost its function in vivo as compared with its WT counterpart.

https://cdn.elifesciences.org/articles/85649/elife-85649-fig3-data2-v2.zip
Figure 3—figure supplement 1
Mad2l1bp variant lost its function in silencing the spindle assembly checkpoint (SAC) in oocytes.

(A) Immunofluorescence co-staining of BUB3 (mitotic checkpoint complex [MCC]) and CREST (kinetochore) using chromosomes spreads from oocytes co-injected with Mad2l1bp WT and oocytes co-injected with Mad2l1bp Mut, at 16 hr after release from milrinone. Schematic diagrams of bivalents at metaphase I (MI) or univalent chromosomes at MII are shown on the right side to aid interpretation. Scale bar, 5 μm. (B) Calculation of the percentages of BUB3 (+) bivalents, BUB3 (-) bivalents, BUB3 (-) univalent, and BUB3 (+) univalent for oocytes at 16 hr after release from milrinone. The numbers of oocytes in each group were indicated. (+), positive; (-), negative. (C) Comparison of the average BUB3 signal intensities among oocytes co-injected with Mad2l1bp WT and Mad2l1bp Mut.

Figure 4 with 1 supplement
Meiotic rescue of the metaphase I (MI)-arrested frozen human oocytes by MAD2L1BP cRNA microinjection.

(A) The resumption of polar body extrusion through MAD2L1BP cRNA microinjection into the frozen oocytes of the patient in Family 1 (II-1). The extrusion of the first polar body indicates the completion of meiosis I of the oocytes. In the control group, the retrieved MI oocytes remained arrested at MI after sham microinjection (upper panel). Six oocytes from the patient in Family 1 (II-1) were microinjected with full-length MAD2L1BP cRNAs, of which four successfully finished extrusion of polar body 1 (PB1) (lower panel). The red arrows point to PB1. Scale bar, 80 μm. (B) A summarized working model depicting how MAD2L1BP(p31comet) deficiency causes human oocyte meiotic arrest. The mitotic checkpoint complex (MCC) is well known to comprise four core members: BUBR1, BUB3, CDC20, and MAD2. MAD2L1BP(p31comet) is known as the core adaptor that bridges the interaction between MAD2 and TRIP13 (Yang et al., 2007; Alfieri et al., 2018). The MAD2L1BP (p31comet)-mediated interplay between MAD2 and TRIP13, known as the MAD2•MAD2L1BP(p31comet)•TRIP13 axis, disassembles the MCC signaling, which consequently drives the meiotic progression during oocyte meiosis. Top panel: The spindle assembly checkpoint (SAC) is switched ‘ON’ until all the chromosome kinetochores are correctly connected to the spindles. Bottom panel: In the presence of WT MAD2L1BP, SAC is turned ‘OFF’ through the recruitment of MAD2-MAD2L1BP-TRIP13 signaling to silence MCC, leading to the meiotic cell-cycle progression. By contrast, in the presence of mutated MAD2L1BP, MCC is unable to be silenced, resulting in persistent SAC activation and oocyte meiotic arrest.

Figure 4—figure supplement 1
Modified Smart-seq2 analysis revealed that MAD2L1BP variant impaired human oocyte meiotic progression.

(A) Bar graph showing the average numbers of transcripts detected in two MI_Control and two MI_MAD2L1BPR285* oocytes (TPM >1). Data are presented as mean ± SEM, n=2. (B) Volcano plot of Smart-seq2 data obtained from human oocytes in MI_Control and MI_MAD2L1BPR285* groups. ‘Up’ and ‘Down’ represent the up-regulated and down-regulated transcripts in MI_MAD2L1BPR285* oocytes. Cutoff: fold change (FC) ≥2, adjusted p≤0.05. (C–D) Gene Ontology (GO) enrichment analysis of down-regulated (C) and up-regulated (D) genes in arrested MAD2L1BPR285* oocytes compared with healthy oocytes. (E) The expression levels of mitotic checkpoint complex (MCC) components and TRIP13 in both groups. Data are presented as mean ± SEM, n=2. Statistical analyses were performed using R software. Detailed p-value as indicated. (F) The expression levels of selected gene transcripts relevant to cell division. Data are presented as mean ± SEM, n=2. Statistical analyses were performed using R software. Detailed p-value as indicated. (G) The expression levels of selected genes related to mitochondria. Data are presented as mean ± SEM, n=2. Statistical analyses were performed using R software. Detailed p-value as indicated.

Tables

Table 1
Clinical characteristics of affected individuals and their retrieved oocytes.
Age (years)Duration of infertility (years)IVF/ICSI cyclesProtocolTotal number of oocytes retrievedGV oocytesMI oocytesPB1 oocytesFertilized oocytesCleaved embryos
Family 1 (II-1)
2221(IVF)Modified ultra-long11011000
2(IVF)Mild stimulation615000
3(IVF)PPOS13013000
Family 2 (II-1)
3161(IVF)PPOS909000
Family 3 (II-1)
3491(IVF)Natural101000
2(IVF)Short404000
  1. GV, geminal vesicle; MI, metaphase I; PB1, first polar body; PPOS, progestin-primed ovarian stimulation.

Table 2
Summary of the variants identified after filtering by whole-exome sequencing (WES) in patients from three families (F1: II-1, F2: II-1, and F3: II-1).
F1: II-1F2: II-1F3: II-1
Total variants134,159392,39266,914
After excluding variants reported in dbSNP, 1000 genomes, EVS, ExAC, and gnomAD (MAF <0.01)361914,02322,832
Exonic nonsynonymous or splicing variants, or coding indels2735719660
Homozygous or compound heterozygous (excluding X and Y chromosomes)121411
Homozygous653
In homozygous region >2 Mb1
Table 2—source data 1

Homozygous and compound heterozygous variants identified by whole-exome sequencing (WES) that survived filtering in patient F1: II-1, F2: II-1, and F3: II-1.

https://cdn.elifesciences.org/articles/85649/elife-85649-table2-data1-v2.zip
Table 3
MAD2L1BP pathogenic variants observed in the three families.
FamiliesGenomic positioncDNA changeProtein changeMutation typeSIFT*PPH2*Mutation taster*NNsplice*ASSP*gnomAD allelegnomAD homozygotes
1Chr6: 43608202c.853C>Tp.R285*NonsenseNANADNANA3/2485980/248598
2Chr6: 43607867c.518delTp.F173Sfs4*Frameshift deletionNANADNANANot foundNot found
Chr6: 43600837c.21-94G>ASplicingNANANA0.13>0.07.1>2.2111/312502/31250
3Chr6: 43607890c.541C>Tp.R181*NonsenseNANADNANA2/2508500/250850
  1. Abbreviations are as follows: D, deleterious; NA, not available.

  2. *

    Mutation assessment by SIFT, Polyphen-2 (PPH2), Mutation Taster, NNsplice, and ASSP.

  3. Frequency of corresponding mutations in all population of gnomAD.

  4. The variant c.21-94G>A is predicted to introduce an alternative splice acceptor by NNsplice (the increase of the acceptor site score from 0 up to 0.13) and by ASSP (the increase of the acceptor site score from <2.2 up to 7.1).

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Escherichia coli)DH5α competent E. coliThermo FisherCat # 18265017
Genetic reagent (Mus musculus)C57BL/6JJackson LaboratoriesRRID: IMSR_JAX:000664
Cell line (Homo sapiens)HEK293T cellsATCCCat # CRL-3216;
RRID: CVCL_0063
Medium: DMEM+10% FBS
Recombinant DNA reagentpcDNA3.1-human MAD2L1BP isoform1This paperpcDNA3.1-hMAD2L1BP iso1Synthesized from Genewiz. Available from Jianqiang Bao’s lab
Recombinant DNA reagentpcDNA3.1-human MAD2L1BP isoform2This paperpcDNA3.1-hMAD2L1BP iso2Synthesized from Genewiz. Available from Jianqiang Bao’s lab
Recombinant DNA reagentM46-mouse Mad2l1bp full-length-FlagThis paperM46-mMad2l1bp-FlagAvailable from Jianqiang Bao’s lab
Recombinant DNA reagentp3xFLAG-CMV-24-mouse Mad2l1bp truncated (1-254aa, equivalent to human p. R285*)This paperp3xFLAG-mMad2l1bp MutAvailable from Jianqiang Bao’s lab
Recombinant DNA reagentp3xFLAG-CMV-24-mouse Mad2l1bp full-length(1-276aa)This paperp3xFLAG-mMad2l1bp WTAvailable from Jianqiang Bao’s lab
Recombinant DNA reagentpcDNA3.1-mouse Mad2 full-length-MycThis paperpcDNA3.1-mMad2-MycAvailable from Jianqiang Bao’s lab
Recombinant DNA reagentpcDNA-human MAD2L1BP-WT minigeneThis paperpcDNA-hMAD2L1BP-WT minigeneAvailable from Jianqiang Bao’s lab
Recombinant DNA reagentpcDNA-human MAD2L1BP-c.21-94G>A minigeneThis paperpcDNA-hMAD2L1BP-c.21–94G>A minigeneAvailable from Jianqiang Bao’s lab
AntibodyAnti-MAD2L1BP (Mouse monoclonal)Santa CruzCat # sc-134381WB (1:200)
AntibodyAnti-Flag
(Mouse monoclonal)
ProteintechCat # 66008-3-IgWB (1:2000)
AntibodyAnti-Myc
(Mouse monoclonal)
ProteintechCat # 60003-2-IgWB (1:5000)
AntibodyAnti-Beta Actin
(Mouse monoclonal)
ProteintechCat # 66009-1-IgWB (1:10,000)
AntibodyHRP-conjugated goat anti-mouse IgG heavy chain
(Goat polyclonal)
ABclonalCat # AS064WB (1:10,000)
AntibodyHRP-conjugated goat anti-mouse IgG light chain
(Goat polyclonal)
ABclonalCat # AS062WB (1:10,000)
AntibodyAnti-BUB3
(Rabbit polyclonal)
Santa CruzCat # sc-28258IF (1:300)
AntibodyAnti-CREST
(Human, unknown clonality)
Fitzgerald Industries InternationalCat # 90C-CS1058IF (1:100)
AntibodyGoat Anti-Human lgG H&L (FITC)
(Goat polyclonal)
ZENBIOCat # 550022IF (1:500)
AntibodyCoraLite 594 Anti-Rabbit
(Goat polyclonal)
ProteintechCat # SA00013-4IF (1:500)
Software, algorithmR softwareR Project for Statistical ComputingRRID:SCR_001905
http://www.r-project.org/
Software, algorithmDAVIDDAVIDRRID:SCR_001881
https://david.ncifcrf.gov/
Software, algorithmFiji/ImageJFijiRRID:SCR_002285
http://fiji.sc
Software, algorithmJalviewJalviewRRID:SCR_006459
https://www.jalview.org/
Software, algorithmPrismGraphPad PrismRRID:SCR_002798
http://www.graphpad.com/
Software, algorithmPyMOLPyMOLRRID:SCR_000305
http://www.pymol.org/
OtherGenBankNIHRRID:SCR_002760
https://www.ncbi.nlm.nih.gov/genbank/
Database
OthergnomADGenome Aggregation DatabaseRRID:SCR_014964
http://gnomad.broadinstitute.org/
Database
OtherOMIMOMIMRRID:SCR_006437
http://omim.org
Database
OtherASSPASSPhttp://wangcomputing.com/assp/Web resource
OtherEvolViewEvolViewhttp://evolgenius.infoWeb resource
OtherMutation
Taster
MutationTasterRRID:SCR_010777
http://www.mutationtaster.org/
Web resource
OtherNNSpliceNNSplicehttp://www.fruitfly.org/seq_tools/splice.htmlWeb resource
OtherPolyPhen-2PolyPhen: Polymorphism PhenotypingRRID:SCR_013189
http://genetics.bwh.harvard.edu/pph2/
Web resource
OtherSIFTSIFTRRID:SCR_012813
https://sift.bii.a-star.edu.sg/
Web resource

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Lingli Huang
  2. Wenqing Li
  3. Xingxing Dai
  4. Shuai Zhao
  5. Bo Xu
  6. Fengsong Wang
  7. Ren-Tao Jin
  8. Lihua Luo
  9. Limin Wu
  10. Xue Jiang
  11. Yu Cheng
  12. Jiaqi Zou
  13. Caoling Xu
  14. Xianhong Tong
  15. Heng-Yu Fan
  16. Han Zhao
  17. Jianqiang Bao
(2023)
Biallelic variants in MAD2L1BP (p31comet) cause female infertility characterized by oocyte maturation arrest
eLife 12:e85649.
https://doi.org/10.7554/eLife.85649