1. Subjects
The patient, a 27-year-old woman, had been married for over a year and presented with infertility despite not using contraception. She had a history of Polycystic Ovary Syndrome (PCOS) and a menstrual cycle ranging from 23 to 37 days, with a 3-day menstrual period, light flow, and occasional dysmenorrhea. The patient had undergone ovulation stimulation treatment at the Reproductive Medicine Center of Guizhou Provincial People's Hospital, resulting in dominant follicle maturation and ovulation. However, she remained infertile after regular intercourse. She had not undergone in vitro fertilization (IVF) or intracytoplasmic sperm injection (ICSI) and denied consanguinity in her marriage.
The patient's family history revealed four women with a history of infertility: II3, III5, III8, and IV1. In contrast, III3 and III7 had children (see Fig. 1a for details).
In this study, peripheral blood samples were collected from the proband (IV1), her father (III2), mother (III1), and aunt for genetic testing. The remaining family members declined to participate. The study was approved by the Ethics Committee of Guizhou Provincial People's Hospital, and all participants provided written informed consent to participate.
2. Whole Exome Sequencing (WES) and Bioinformatics Analysis
Whole exome sequencing (P039-Exome) was performed on the proband with a target mean depth of 104.99 and target region coverage of 98.79. The sequenced reads were aligned to the UCSC hg19 human reference genome using the BWA algorithm (https://www.sentieon.com/) within Sentieon software (https://www.sentieon.com/). Single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) were identified using Sentieon’s variant calling tool. The detected variants were annotated using ANNOVAR software (http://annovar.openbioinformatics.org/en/latest/), which incorporates data from HGMD® and ClinVar (https://www.ncbi.nlm.nih.gov/clinvar). The pathogenicity of the identified variants was predicted using bioinformatics tools, including SIFT (http://sift.jcvi.org/), PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/), MutationTaster (https://www.mutationtaster.org/), and GERP++ (http://mendel.stanford.edu/SidowLab/downloads/gerp/index.html).
3. Sanger Sequencing
To validate the TUBB8 gene mutation identified by WES, Sanger sequencing was performed on the proband, her parents, and her aunt to confirm the presence of potential pathogenic variations. The forward and reverse primers used for sequencing were 5'-CCGTGGCATCCTGATATTGC-3' and 5'-TGACTGTGGCTGAGCTTACC-3', respectively. PCR amplification products were sequenced using a 3130XL sequencer via Sanger sequencing.
4. Conservation Analysis and Protein Structure Prediction
To assess the conservation of the mutant amino acid sequences, we utilized the MEGA11 software tool (https://www.megasoftware.net/dload_win_beta). The three-dimensional structures of both wild-type (WT) and mutant (MUT) proteins were predicted using Phyre2 software (http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgiid=index). Additionally, ChimeraX software (https://www.cgl.ucsf.edu/chimerax/) was used to analyze the structural and physicochemical differences between the wild-type and mutant proteins at the mutation sites.
5. Immunofluorescence
Plasmids TUBB8-WT and TUBB8-A1039G-MUT were constructed using a 3xFlag-pcDNA3.1 vector, expressing the FLAG tag fusion protein. These plasmids were transiently transfected into HeLa cells using Lip3000 (HB-LF3, Hanhen Bio, China). After a 48-hour incubation post-transfection, cells were fixed with 4% formaldehyde at 4°C for 30 minutes, then washed with PBST. Permeabilization was performed using 0.2% Triton X-100 for 20 minutes at room temperature, followed by blocking with 5% BSA for 1–2 hours. Immunostaining was conducted using Mouse anti-FLAG tag primary antibody (1:200 dilution, Abclonal, AE005), 594-conjugated Goat anti-Mouse IgG (H + L) (1:200 dilution, AS054, ABclonal, China), or FITC-conjugated Goat anti-Mouse IgG (H + L) (1:300 dilution, AS001, ABclonal, China). Nuclei were stained with DAPI (10 µg/ml, C0065, Solarbio, China) for 10 minutes. The cellular microtubule structure and protein localization were observed using a confocal microscope (LSM980, ZEISS, Germany), and fluorescence intensity was analyzed using ImageJ software.
6. RT-PCR
Plasmids TUBB8-WT and TUBB8-A1039G-MUT were constructed using the pEGFP-N1 empty plasmid and verified by Western blotting. HEK293T cells or HeLa cells were seeded in six-well plates at a density of 4 × 10^5 cells per well, 24 hours before transfection. The plasmids were transfected into the cells using Lip3000 (HB-LF3, Hanhen Bio, China), and RNA was extracted 36 hours post-transfection. Reverse transcription to cDNA was performed using the RT Master Mix for qPCR II (gDNA digester plus) (HY-K0511A-1, MCE, China), followed by real-time quantitative PCR using the TB Green® Premix Ex Taq™ (Tli RNaseH Plus) (RR420A, TAKARA, Japan) on a real-time PCR instrument (CFX OPUS 96, BIO-RAD, USA). The forward and reverse primers for TUBB8 were: TUBB8-F: 5'-TTCAGGCCAGACAACTTCATCTT-3' and TUBB8-R: 5'-CTCCTTTCTGACAACGTCCATCA-3'. The internal control primers were: ACTB-F: 5'-CCACGAAACTACCTTCAACTCCATC-3' and ACTB-R: 5'-AGTGATCTCCTTCTGCATCCTGTC-3'. The Q-PCR reaction mixture (25 µL) consisted of cDNA (2 µL), forward primer (10 µM, 1 µL), reverse primer (10 µM, 1 µL), 2× TB Ex Taq premix (12.5 µL), and sterile water (8.5 µL). The relative gene expression levels were calculated using the 2^−ΔΔCt method (ΔΔCt = ΔCt(treatment group) - ΔCt(control group), ΔCt = Ct(target gene) - Ct(reference gene)).
7. Western Blot
Plasmids TUBB8-WT and TUBB8-A1039G-MUT were constructed using the pEGFP-N1 empty vector and verified by Western blotting. HEK293T or HeLa cells were seeded in six-well plates at a density of 4 × 10^5 cells per well, 24 hours before transfection. The plasmids were transfected into the cells using Lip3000 (HB-LF3, Hanhen Bio, China). Protein extraction was performed 48 hours post-transfection, and Western blotting was conducted. The primary antibodies used were Mouse anti-GFP tag (AE012, Abclonal, China) and Mouse anti-GAPDH (AC033, Abclonal, China), both diluted to 1:4000 and 1:40,000, respectively. The secondary antibody, Goat anti-Mouse IgG HRP (AS003, Abclonal, China), was diluted to 1:4000. Protein bands were visualized, and signals were analyzed using ImageJ software.
8. Statistical Analysis
Data were analyzed using GraphPad Prism software. Values are presented as the mean ± standard error of the mean (SEM), and statistical significance was determined using the t-test. A P-value of less than 0.05 was considered statistically significant.