Whole-exome sequencing identified a novel mutation of AURKC in a Chinese family with macrozoospermia

Purpose Macrozoospermia is a rare sperm morphologic abnormality associated with male infertility and is characterized by a high percentage of spermatozoa with large irregular heads. The aim of this study was to identify the genetic cause of an infertile male with macrozoospermia from a consanguineous family. Methods Whole-exome sequencing (WES) was performed using peripheral blood genomic DNA from the patient and his parents. Results WES analysis of the patient with macrozoospermia from a consanguineous family allowed the identification of a novel homozygous missense variant in the AURKC gene (c.269G>A). Bioinformatics analysis also suggested this variant a pathogenic mutation. Quantitative real-time PCR analysis showed that the mRNA level of AURKC is significantly decreased in the patient compared with his father. Moreover, no embryos were available for transfer after ICSI. Conclusions These results further support the important role of AURKC in male infertility and guide the practitioner in optimal decision making for patients with macrozoospermia. Electronic supplementary material The online version of this article (10.1007/s10815-018-1374-3) contains supplementary material, which is available to authorized users.


Introduction
Approximately 70 million couples worldwide suffer from infertility [1], and approximately half of these cases are due to male factors [2]. Macrozoospermia is a rare sperm morphologic abnormality associated with male infertility and is characterized by a high percentage of spermatozoa with large irregular heads [3,4]. This syndrome, first reported in 1977, affects < 1% of infertile men [3]. It is considered to be an autosomal recessive type of teratozoospermia that results in male infertility [5]. Currently, the most relevant single-gene defect that has been identified in a patient with macrozoospermia is a mutation in Aurora Kinase C (AURKC) [5][6][7].
Here, we report a novel homozygous missense variant in the AURKC gene in an infertile male from a consanguineous family identified via whole-exome sequencing. In addition, this variant has a high probability of pathogenicity according to in silico analysis. ICSI was also performed and failed to generate any embryos suitable for transfer. These results further support the important role of AURKC in male infertility and guide the practitioner in optimal decision-making for patients with macrozoospermia.

Patient
The proband was a 27-year-old treated for infertility at Anhui Provincial Hospital. The results of the patient's semen tests are described in Table 1. The results revealed that the patient suffered from macrozoospermia (close to 100% large-headed spermatozoa with a sperm count of 1 M/ml) ( Table 1). The parents of the patient were first-degree cousins (Fig. 1). The patient exhibited normal erection and ejaculation, and reported having regular sexual intercourse 2-3 times/week; however, his wife had been unable to get pregnant since they were married in 2015. The patient had no history of unhealthy activity or contact with adverse chemicals. Physical examination revealed normally developed male external genitalia, normal bilateral testicular size, and no abnormality in the bilateral spermatic veins upon palpation. The proband did not have primary microcephaly or respiratory disease. The chromosomal karyotype of the patient was normal (46; XY), and no deletion was found in the Y chromosome. The patient hormone levels were normal.
This project was approved by the ethics committee of Anhui Medical University. All study members gave their written informed consent before sample collection.
Genomic DNA extraction and whole-exome sequencing (WES) The genomic DNA from the patient and his parents was extracted from peripheral blood samples using a QIAamp DNA blood midi kit (Qiagen, Hilden, Germany) according to the manufacturer's protocol.
Genomic DNA from the proband and his parents was subjected to WES. WES was performed by the WuXi NextCODE in Shanghai on a HiSeq2000 sequencing platform (Illumina, San Diego, CA, USA). After the adaptors were removed, the WES raw reads were aligned to the reference genome Hg19 using the Burrows-Wheeler aligner, followed by removal of PCR duplicates. Variants including single-nucleotide polymorphisms and indels were identified using SAMtools and annotated by ANNOVAR software. A candidate gene was considered a variant that fulfilled the following criteria: (i) missense, nonsense, frame shift, and splice site variants, (ii) absent or rare (frequency below 1%) in the two databases (dbSNP, 1000G), and (iii) homozygous variants in the patient and heterozygous variants in his parents.

Sanger sequencing validation
The mutation in AURKC in the proband and his parents was validated using Sanger sequencing. We amplified the PCR products of exon 3 of the AURKC gene using specific primes (the forward primer is 5′-AACCAGGATTCGAGTGTCTG-3′, and the reverse primer is 5′-CAATCTCCAGGTAG ACGATGGAG-3′). Then, the PCR products were sequenced on an ABI 3730XL automated sequencer (Applied Biosystems, Forster City, CA, USA).

RNA extraction and Q-PCR
RNA extraction was carried out on the whole blood using TRI REAGENT®BD (Molecular Research Center) using the manufacturer's protocol. Reverse transcription was performed in patient and his parents with 5 μl of extracted RNA. Hybridization of the oligo-dT was performed by incubating for 30 min at 42°C and quenching on ice with the following mix: 5 μl of RNA, 10 μl of 2Xsupermix (10 mM, Pharmacia), 1 μl of gDNA (0.5 mM, Roche diagnostics), and 4 μl of H 2 O. Then 2 μl of the obtained cDNA mix was used for the quantitative PCR (Q-PCR) using a StepOne-PlusTM Real-Time PCR System (Life Technologies) with Power SYBR Green PCR Master Mix (Life Technologies) according to the manufacturer's protocol. Quantification of the fold change in gene expression was determined by the relative quantification method (2 −ΔΔCT ) using the gapdh gene as a reference. Data are shown as the average fold increase standard error of the mean. Primers are described in supplementary file 1.

ICSI, embryo, and evaluation of embryo quality
The patient and his wife underwent intracytoplasmic sperm injection (ICSI) at Anhui Provincial Hospital. Briefly, embryo culture was performed using Vitrolife G-SERIES™ culture media (Vitrolife, Goteborg, Sweden) according to the manufacturer's instructions. The patient's wife had undergone one ICSI cycle. Detailed results are described in the BResults^section.

Sequence alignment of AURKC protein
Sequence alignment of the AURKC protein in different species was conducted using ClustalX2.1. The number of each species was as follows: Homo sapiens (NP_001015878.1), Mus

WES analysis of a patient with macrozoospermia
To identify the genetic cause of the macrozoospermia, we performed WES using peripheral blood genomic DNA from the patient and his parents to identify the putative pathogenic mutation. Considering the family history of consanguinity, we focused on homozygous mutations. After of the exclusions of frequent variants and application of technical and biological filters (see BMaterials and methods^section and Fig. 2A), a limited list of homozygous mutations was established (Supplementary Table 1). To determine whether any of these 36 genes may be related to macrozoospermia, we first examined the tissue expression patterns for each gene. Among these homozygous variant genes, only two genes exhibited testisenriched expression (Fig. S1). One of the identified genes was AURKC, and a new mutation (c.269G>A, GenBank accession number, NM_001015878) in exon 3 was identified in the patient by WES. Accordingly, the change in amino acids was determined to be Arg90Gln. Another identified gene was gametogenetin (GGN), and a mutation (c.148T>C, GenBank accession number, NM_152657.3) in exon 3 was identified in the patient by WES. Accordingly, the change in amino acids was determined to be Trp50Arg.

Validation by Sanger sequencing
We next performed Sanger sequencing to verify the homozygous variant in the AURKC and GGN genes in both the patient and his parents. Exon 3 of the AURKC gene was amplified from genomic DNA of the patient and his parents using PCR. The homozygous missense variant in AURKC was validated in the patient, and his parents had heterozygous alleles (Fig. 2B).

Detrimental effects of the identified variant in AURKC
Unfortunately, we could not get access to testicular biopsy samples from the affected patient and his parents; therefore, we set up Q-PCR to study the effect of the variation. Q-PCR analysis demonstrated that mRNA levels of AURKC were down-regulated significantly in the patient (Fig. 2C).

In silico analysis of the mutation
Bioinformatics analysis of the c.269G>A mutation in the AURKC gene by three online pathogenicity prediction tools (Polyphen-2, SIFT, Mutation taster) suggested that this mutation is most likely a disease-causing mutation (Supplementary Table 2). This variant is a novel mutation that was absent in   Fig. 2 Whole-exome sequencing analysis of the consanguineous family. A Filter strategies used in this study. B Sanger sequencing validation of the affected patient and his parents. The arrow indicates the mutation site. C Q-PCR analysis of mRNA expressions of AURKC in the patient and his father. Messenger RNA expression determined by Q-PCR is calculated as a ratio relative to gapdh and expressed relative to his father. D Sequence alignment of AURKC protein in different species. The arrow refers to the mutation site the ExAC, 1000G, and gnomAD databases. The mutation site at R90 was highly conserved from human to zebrafish, indicating an important role of this site for the function of AURKC protein (Fig. 2D).

ICSI using the patient's sperm
After centrifugation on density gradients and careful examination, a few Bnormal-looking^spermatozoa that fit into an ICSI micropipette were selected. One ICSI cycle was attempted for the patient at Anhui Provincial Hospital. For the ICSI cycles, we collected 10 eggs and 7 MII oocytes; none of the embryos developed to the blastocyst stage.
To date, numerous reports have provided evidence suggesting that failures in chromosome segregation and/or cytokinesis during the first, second, or both meiotic divisions are the major cause of macrozoospermia [15][16][17][18][19]. Aurora kinases are highly evolutionarily conserved kinases that are required for accurate chromosome segregation and cytokinesis during meiosis [20]. AURKC has 7 exons and encodes a 309 amino acid protein in humans. AURKC is a component of the Aurora kinase family and is predominantly expressed in the testis [21,22]. In this study, we found a new missense mutation in exon 3 of the AURKC gene that led to an amino acid change (p.R90Q). The sequence alignment of the AURKC protein showed that this mutation site is conserved among different species (Fig. 2C). We utilized three online pathogenicity predict tools (Polyphen-2, SIFT, Mutation taster) to predict the harmfulness of this variant ( Table 2). The results suggested that this variant has a high probability of pathogenicity.
Previous studies have compared the values of a routine spermogram and spermocytogram between patients with and without an AURKC mutation [8]. The results showed that the proportions of large-headed spermatozoa generally reached far higher values, and the presence of 1% normal spermatozoa is the most discriminant parameter in patients with AURKC mutations [8]. Consistent with previous studies, the patient with missense variant in the AURKC gene (c.269G>A) in this study showed 96% large-headed spermatozoa. However, the proportions of normal spermatozoa reached 4%. One explanation could be that in addition to macrozoospermia, this patient studied here also has a low sperm number (1.71 × 10 6 ). Variable scoring of the normal spermatozoa features between operators and laboratories could also influence the results as some spermatozoa with a slight morphology defect in head or flagella could be considered normal [8].
Numerous publications have described a failure of pregnancy in patients with macrozoospermia [5,7,12,[16][17][18][23][24][25]. In addition, further studies have recommended that systematic genetic screening of the AURKC gene be performed when a patient has more than 30% enlarged head spermatozoa. If a mutation in the AURKC gene is found, ICSI should not be attempted [7,26]. In this study, the patient's wife failed to become pregnant after one ICSI cycle, again demonstrating the correlation between AURKC gene mutation and ICSI outcome. Consequently, ICSI in such a patient is not recommended.
In conclusion, we identified a novel missense mutation (c.269 G>A; p.R90Q) in the AURKC gene. To date, this is the sixth reported variant in the AURKC gene associated with macrozoospermia. This study expands the spectrum of AURKC mutations and helps guide the practitioners to make optimal decisions for patients with macrozoospermia.