Copy number variations (CNVs) and karyotyping analysis in males with azoospermia and oligospermia

Background Considering the essential roles that genetic factors play in azoospermia and oligospermia, this study aims to identify abnormal chromosomes using karyotyping and CNVs and elucidate the associated genes in patients. Methods A total of 1157 azoospermia and oligospermia patients were recruited, of whom, 769 and 674 underwent next-generation sequencing (NGS) to identify CNVs and routine G-band karyotyping, respectively. Results First, 286 patients were co-analyzed using CNV sequencing (CNV-seq) and karyotyping. Of the 725 and 432 patients with azoospermia and oligospermia, 33.8% and 48.9% had abnormal karyotypes and CNVs, respectively. In particular, 47,XXY accounted for 44.18% and 26.33% of abnormal karyotypes and CNVs, respectively, representing the most frequent genetic aberration in azoospermia and oligospermia patients. Nevertheless, big Y and small Y accounted for 7.46% and 16.67% of abnormal karyotypes, respectively. We also identified high-frequency CNVs-loci, such as Xp22.31 and 2p24.3, in azoospermia and oligospermia patients. Conclusion Sex chromosome and autosomal CNV loci, such as Xp22.31 and 2p24.3, as well as the associated genes, such as VCX and NACAP9, could be candidate spermatogenesis genes. The high-frequency abnormal karyotypes, CNV loci, and hot genes represent new targets for future research. Supplementary Information The online version contains supplementary material available at 10.1186/s12920-023-01652-2.

Approximately 10% of male infertility patients are categorized as non-obstructive azoospermia (NOA) due to failure in spermatogenesis.Microdeletions and chromosomal abnormalities of the Y chromosome are the key genetic factors in NOA [9], with incidences of 9.7% and 13% in males with azoospermia, respectively [10].Chromosomal abnormalities that cause male infertility clinically manifest as azoospermia, oligospermia, teratozoospermia, and globozoospermia.That is, numerical and structural changes in sex chromosomes can cause impaired spermatogenesis, whereas autosomal structural abnormalities, including balanced reciprocal translocations on chromosomes 1, 6, 12, and 22, can also lead to azoospermia or severe oligospermia [11].
Small deletions of the long arm of the Y chromosome (Yq) are not visible through conventional banding karyotype.Nevertheless, deletions in the azoospermia factor (AZF) region of the Yq arm among men manifest as azoospermia and oligospermia.Four spermatogenesis loci have been characterized in the Yq11 head, middle, and terminus, namely, AZFa, AZFb, AZFc, and AZFd [2,12].However, more than 4000 genes contribute to human spermatogenesis [13].
The microdeletion and microduplication of chromosomal fragments can be detected by analyzing the CNVs [14].CNV-seq-based on NGS technology for low-depth whole genome sequencing, which detects chromosomal aneuploidy, chromosomal CNVs over 0.1 Mb, polyploidy, and mosaicism as low as 5% [15,16].Moreover, CNVseq has a low associated cost, short run period, high resolution, high accuracy, and requires a small amount of input DNA.Meanwhile, karyotyping, with a maximum resolution of 3 Mb, is the gold standard for identifying chromosomal abnormalities to detect whole and partial chromosomal aneuploidies, unbalanced translocations, polyploidies, and mosaicism [16,17].
This study analyzed the NGS-based CNV-seq and karyotyping in males with azoospermia and oligospermia to identify novel microdeletion and microduplication targets, as well as genes and abnormal chromosomal karyotypes associated with azoospermia and oligospermia.

Inclusion criteria
Healthy couples that had not become pregnant after one year of sexual intercourse without contraception, and who showed abnormal semen counts upon semen examination were included in this study.

Exclusion criteria
Adult males with a history of genital tract obstruction or dysfunction (varicocele and obstructive azoospermia) or congenital defects in the structure of the urogenital system were excluded.Adult male patients who had received radiotherapy for reproductive system carcinoma or chemotherapy for malignancies were also excluded.

Semen analysis
All sperm samples were tested using standard procedures in the laboratory department of Tongji Hospital in Wuhan, China.Semen samples were collected through masturbation after 3-5 days of abstinence, and semen characteristics were evaluated within 1 h after ejaculation.Semen analysis was performed in accordance with the World Health Organization (WHO) standard protocol for diagnosis of azoospermia (absence of sperm in semen) and oligospermia (sperm counts < 15 × 10 6 ml − 1 ).

Study population
We recruited 1157 male infertility patients with azoospermia or oligospermia from the obstetrics and andrology department of Tongji Hospital, Wuhan, China, between February 2017 and July 2019.All CNV-seq and karyotype analyses of infertile male samples were performed at the perinatal laboratory of Tongji Hospital (Wuhan, China).

Chromosome karyotype analysis of peripheral blood
Approximately 2 ml of peripheral blood was collected in heparin vacutainers for chromosomal karyotype analysis.For each sample, 0.5 ml of peripheral blood was cultured in 5 ml complete medium (RPMI 1640 media + 15% fetal calf serum) for 72 h.The basic cytogenetic method included chromosome harvesting, slide preparation, staining, and banding; subsequently, chromosome number and banding patterns were evaluated.Chromosome preparation involved arresting the cells, harvesting and preparing a single cell suspension, hypotonic treatment, and fixation.The culture process was performed according to standard procedures [18].G-banded metaphase chromosomes were obtained by treating the cells with trypsin and staining with Giemsa stain.In all cases, at least 30 cells in metaphase were analyzed.In mosaicismsuspected cases, 50 cells were counted.The karyotypes were interpreted according to the International System for Human Cytogenetic Nomenclature (ISHCN) recommendations.An abnormal clone was defined as exhibiting structural aberration, or chromosomal gain in at least two metaphases, and chromosome deletion observed in at least three metaphases.

Sample preparation and copy number variation sequencing (CNV-seq)
For CNV-seq analysis, approximately 1-2 ml of peripheral blood was collected in EDTA vacutainers and sent to the perinatal laboratory of Tongji Hospital, Wuhan, China.DNA was extracted from the peripheral blood using a DNeasy Blood and Tissue Kit (Qiagen, Germany).The DNA was quantified using a NanoDrop spectrophotometer (Thermo Fisher Scientific, United States) (concentration > 50 ng/µl; OD260/280 > 1.8; OD260/230 > 1.5).DNA samples were sent to Beijing Berry Genomics Co., Ltd., for CNV sequencing using the NextSeq CN500 platform (Illumina Inc.).The annotation and interpretation were carried out based on the guidelines of the American College of Medical Genetics and Genomics.The results were determined by referring to the hg19 version of the human genome and the latest data available from the Database of Genomic Variants (DGV), DECIPHER, Online Mendelian Inheritance in Man (OMIM), University of California Santa Cruz (UCSC), and PubMed.

Genes involved in CNVs
Microdeletion and microduplication of chromosomal fragments detected through CNV-seq may include changes in genetic function and affect gene expression in testicular or other tissues, leading to spermatogenic failure.Hence, the PubMed public data platform of the National Library of Medicine and National Center for Biotechnology Information (NCBI) through Genome Data Viewer was used to search genes of the CNV loci (https://www.ncbi.nlm.nih.gov/genome/guide/human/).

Analysis process
This study recruited 1157 males with infertility, including 725 with azoospermia and 432 with oligospermia.The average ages of patients with azoospermia and oligospermia were 32 ± 0.4 vs. 31 ± 0.6 years (p = 0.57), respectively.The simplified analysis process is shown in Fig. 1.Samples for 769 and 674 patients were evaluated via CNV-seq or chromosome karyotyping, respectively, while those for 286 patients were assessed by both methods.Moreover, 26.92% (77/286) of patients had abnormal CNVs and normal chromosome karyotyping, whereas CNVs were not detected in 7.34% (21/286) of patients that had abnormal chromosome karyotyping (Supplemental Table 1).
Abnormal CNVs were detected in 376 patients, of which 52 carried duplication and deletion CNVs, while 324 patients had either duplication or deletion CNVs.No statistical differences were observed between azoospermia and oligospermia patients in the total deletion and duplication CNVs or in the deletion and duplication CNVs of different chromosomes (p > 0.05; Supplementary Table 2).

Analysis of abnormal chromosomal karyotyping in males with azoospermia and oligospermia
Among the males with infertility, 446/674 cases had normal chromosomal karyotyping, while 228 had abnormal karyotyping.Of these 228 cases, the most frequent abnormal karyotyping was 47,XXY, 103/228 (45.18%) in 100 azoospermia cases and 3 oligospermia cases.Small Y (38/228, 20 azoospermia cases, 18 oligospermia cases), big Y (17/228, 5 azoospermia cases, 12 oligospermia cases), abnormal chromosome Y (25/228), and abnormal chromosome X (2/228) were also of interest.The distribution of abnormal chromosomal karyotyping in azoospermia and oligospermia patients is shown in Supplemental Table 3.The 25 cases with karyotyping related to abnormal chromosome Y are shown in Table 1.

High-frequency duplication CNVs
Abnormal CNVs were detected in 376 of the 769 patients.Among these 376, 341 had duplication CNVs, while the others had deletion CNVs (131 patients).Regarding the * represent chromosome X accompanied by others chromosomal microdeletion and microduplication, was showed in supplemental Table 5.All genes of CNVs covered was showed in supplemental material of excel 1 (sheet 1) deletion and microdeletion CNVs and involved genes of chromosome Y are shown in Table 4.All genes covered by CNV loci are shown in the Supplemental material of Excel 3.

Discussion
The WHO estimates that male factors have a crucial role in half of all fertility cases [19].Male infertility is a complex multifactorial pathological condition that presents as azoospermia and oligospermia.The causes of spermatogenic defects are associated with the environment, lifestyle, nutrition, and genetic factors [20].In fact, in azoospermia, approximately 25% of cases can be associated with genetic factors; however, the genetic factors of male infertility are highly complex [21].
Somatic cells, germ cells, and endocrine signals are essential for the development and maturation of sperm cells.Both sex chromosomal and autosomal genes are involved in the regulation of spermatogenesis.Cytogenetic karyotyping and NGS-based CNV-seq technologies for chromosomal analysis have been widely used in clinical practice.Chromosomal karyotype analysis is the first choice for genetic testing of infertility in males, especially males with azoospermia and oligospermia.The reason is that the karyotype analysis is simple and convenient, with high accuracy and low cost.In this study, Klinefelter's syndrome (KS, 47,XXY) was identified as a common chromosomal abnormality in males with infertility (103/674, 15.28%).However, a 45,X/46,XY karyotype was also observed in two males with azoospermia.Chromosomal polymorphisms are primarily considered * represent chromosome 2 accompanied by others chromosomal microdeletion and microduplication, was showed in supplemental Table 6.All genes of CNVs covered was showed in supplemental material of excel 2 (sheet 1) to have no obvious clinical phenotypes, however, other studies have reported that they are related to infertility and recurrent spontaneous abortion (RSA) [6].Big Y and small Y are considered Y chromosomal polymorphisms associated with RSA.Hence, we speculate that big Y and small Y could be related to azoospermia and oligospermia.Other chromosome Y abnormalities, such as 47,XYY and 46,X,+mar, and structural aberrations of chromosome Y (inversions, translocations, deletions, and additions), as well as Yqh + and Yqh-were also observed in males with azoospermia and oligospermia.We found a rare case of a male with azoospermia who had a 46,XX karyotype.This testicular disorder of sex development is a rare genetic syndrome, causing male infertility due to the absence of the AZF region in chromosome Yq [22].Sex chromosomes and autosomal genes participate in spermatogenesis regulation [20].Meanwhile, 47,XXY (KS) is the main genetic factor of male infertility, accounting for 10-15% of oligospermia and NOA cases [23].In our study, we detected 99 KS patients through CNV-seq, 77 of which were 47,XXY; 72/77 cases showed azoospermia, and 5 exhibited oligospermia.Furthermore, 22 cases (22/99) were 47,XXY accompanied by other chromosomal microdeletions and microduplications, except for one case of 47,XXY with dup(18) (p11.32) which manifested as oligospermia, while the other 21 cases manifested as azoospermia (Supplemental Table 4).Chromosome X is enriched in genes associated with spermatogenesis.Infertile males reportedly have Xp22.31,Xp22.1, Xp11.23, and Xq28 CNV loci [24].In this study, 24 patients with chromosome X microduplication CNVs were detected, of which 7 were accompanied by other chromosomal microdeletions and microduplications (Table 2).Moreover, 5 of 6 patients with dup(X)(p22.31)exhibited azoospermia.As such, we speculate that the Xp22.31 loci might be associated with male fertility.The regions involved VCX, VCX3A, GPC3, PUDP, SOX3, and HSFX2, among others.X-linked/Y-linked (VCX/Y) genes encode a group of proteins expressed exclusively in male germ cells among normal tissue.The family members include VCX, VCX2, VCX3A, VCX3B, VCY, and VCY1B.All copies of VCX and VCY are transcribed exclusively in the testis, likely within germ cells [25].Moreover, the VCX/Y protein family might participate in ribosome assembly regulation during spermatogenesis [26].VCX may act by mediating the mitochondria-dependent apoptosis pathway or p53-Bax pathway to prompt cell apoptosis and inhibit cell growth, delaying cell progression in G1 to S transition, and leading to spermatogenesis impairment and failure [27].Some genes have significance in spermatogenesis, while the functions of others are unknown.For example, SOX3 is important for regulating embryonic development and shows homology to the SRY gene.There is enrichment in male-specific genes on the X chromosome [20].Therefore, the microdeletion and microduplication of CNVs on the X chromosome are responsible for azoospermia and oligospermia.The severity of the clinical phenotype is related to the length of deletion and duplication bases, the location of the region, and whether other chromosomal microdeletions and microduplications are present.In these patients, we detected that the most frequent autosomal CNVs was microduplication on chromosome 2 (30 cases).Moreover, 7/30 cases had dup(2)(p24.3)(chromosome 2: 13,520,000-14,340,000), involving NACAP9, LOC105373438, and LINC00276, among others (Table 3).We considered that the CNV loci of 2p24.The microdeletions in the AZF region on chromosome Y are well-characterized causes of male infertility involving azoospermia and oligospermia [9,13].The AZF region contains many critical genes associated with spermatogenesis, germ cell development, and spermatocyte maturation, such as USP9Y, KDM5D, and DAZ [29].Microdeletion loci on chromosome Y, such as Yp11.2,Yq11.21q11.221,and Yq11.223, may be associated with spermatogenesis.Hence, the associated genes, including NLGN4Y, TTTY5, and SRY, require further functional analyses to clarify their relationship with spermatogenesis.With the rapid development of sequencing technology, genomics has provided new tools to better understand the genetics of male infertility.In particular, deletion and microdeletion of the AZF region of chromosome Y have been found to involve genes related to spermatogenesis.Moreover, microdeletion and microduplication of chromosome X, autosomal genes, involved genes, and abnormal chromosomal karyotypes are associated with spermatogenesis and warrant further research.
The limitations of our study include its small sample size and single research center design.Moreover, we did not perform CFTR analysis nor aCGH and Taqman assays to verify detected CNV loci.Hence, since CNVseq cannot detect single base mutations and certain microdeletions or microduplications, exome sequencing is needed in future studies.

Conclusion
Various genetic factors contribute to spermatogenic dysfunction.Deletion of the AZF region of the Yq arm primarily manifests as azoospermia and oligospermia.Meanwhile, small Y and big Y may be related to azoospermia and oligospermia.CNV-seq detected that the AZF region of chromosome Y deletion primarily manifested as azoospermia and oligospermia.Various microdeletion and microduplication CNV loci of chromosome Y and the covered genes require further investigation to elucidate their roles in spermatogenesis.Nevertheless, our CNV-seq analysis detected important CNV loci, such as Xp22.31 and 2p24.3, and the associated genes (VCX3A and NACAP9), which might serve as candidate spermatogenesis CNV loci and hot genes.

Data Availability
The data that support the findings of this study are not openly available due to reasons of sensitivity and are available from the corresponding author

Fig. 1
Fig. 1 Simplified flow diagram of analysis process in this study

Table 1
Distribution of abnormal chromosome Y in males with azoospermia and oligospermia Fig. 2 Frequency of CNVs on different chromosomes in males with azoospermia and oligospermia.Green, red, and blue represent total CNVs, microduplication CNVs, and microdeletion CNVs, respective

Table 3 )
. Other chromosomal microdeletion and microduplications accompanied by chromosome 2 aberrations are shown in Supplemental Table6.All genes covered by CNV loci are presented in the Supplemental material of Excel 2.

Table 2
Details regarding 24 males with azoospermia and oligospermia with chromosome X microduplication CNVs and the involved genes

Table 3
Details regarding 30 males with azoospermia and oligospermia with chromosome 2 microduplication CNVs and the involved genes

Table 4
Details of 35 males with azoospermia and oligospermia with chromosome Y microdeletion CNVs and the involved genes Inheritance in Man This study was funded by National Natural Science Foundation of China and Beijing Kanghua Western and Chinese Medicine Development Foundation (Grant No. 81902299 & 2201102997).