Introduction: Individual risk assessment of assisted reproductive
technologies is essential for personalized treatment strategies. Genetic
and genomic indicators of the response to stress by cells could provide
individual prognostic indicators for in vitro fertilization (IVF)
success. Such indicators include the copy number of ribosomal genes (rDNA), which
modulates the level of protein synthesis, and the abundance of mitochondrial DNA
(mtDNA), which provides the cell with energy, while the content of telomere
repeats (TRs) indicate the biological age. Materials and Methods: The
contents of the three repeats in DNA isolated from blood leukocytes of 40 women
before and after ovarian stimulation were assayed prior to IVF. Then, we divided
the women into a successful IVF group, IVF+ (N = 17, 7 cases of twins), and a
group of failed cases, IVF– (N = 23). The control group included 17 non-pregnant women with natural childbirth in the past. The
nonradioactive quantitative hybridization (NQH) method was applied to assay the
genome repeat contents. Results: The number of rDNA copies in the IVF+
group was significantly higher than in the IVF– group (p
The growing number of infertile couples is a global trend that currently affects approximately 10–15% of couples worldwide. Accordingly, the number of women who resort to in vitro fertilization (IVF) is increasing. Infertility is a very complex multifactorial condition, meaning different couples have different chances of success. The individual outcome depends on the health, infertility causes, age, emotional state, lifestyle, eggs and sperm quality, and many other factors relating to the couple. The IVF procedure often fails because of unaccountable reasons. Even in cases where there is a successful conception, an increased risk of pregnancy termination (stopping the development of the embryo, spontaneous abortion, miscarriage) remains. Therefore, uncovering novel factors for infertility and IVF failure is necessary to determine the feasibility of using the patient’s own reproductive cells and making an adequate choice of precision-assisted reproductive technology.
IVF is an effective and physiologically stressful procedure . Therefore, indicators of the effectiveness of the cellular response to stress can be considered potential predictors of the success of the IVF procedure. Such indicators include the cellular protein synthesis level, biological age, and ability to sustain a stress response with an adequate ATP (adenosine triphosphate) supply.
Protein synthesis is limited by, among other factors, the number of available ribosomes. Consequently, the number of ribosomes depends on the ribosomal gene copy number in the genome [2, 3]. Previously, we conducted a pilot study to show that the number of ribosomal gene (rDNA) copies in the blood leukocytes of the female allowed for the prediction of IVF success . The number of rDNA copies in the genomes of women with several unsuccessful IVF attempts was significantly lower than in the genomes of women who succeeded in getting pregnant. The ribosomal repeat harbors the genes for 18S, 5.8S, and 28S ribosomal RNA (rRNA) and belongs to the class of moderate tandem repeats located in the p-regions of the human acrocentric chromosomes. There is still no clear answer in the literature to the variability in rDNA copy number under stress. There is evidence for both stable rDNA content in the genome and variation of this parameter under stress [5, 6, 7]. The variability of rDNA content in the genome of women during IVF has, to our knowledge, never before been studied.
Telomere repeats (TRs), which also belong to the class of moderate tandem repeats, are situated at the ends of the chromosomes. A decrease in the number of TR copies in the genome also termed telomere length shortening, is associated with aging and oxidative stress [8, 9, 10, 11]. TR is considered one of the cellular aging markers. The IVF success is deemed as depending on age-specific changes in women and the level of oxidative stress. However, published data on the link between the TR content and IVF success are contradictory [12, 13, 14, 15].
Mitochondrial repeats (mtDNA) are a component of mitochondria. The mitochondria provide the cell with energy, which is stored in the form of ATP molecules. The ATP synthesis level depends on the number of mitochondria in the cell, i.e., on the mtDNA copy number. Under stress, the mtDNA abundance changes [16, 17], meaning a lack of energy may be a cause of IVF failure. Some published reports disprove a positive association between the IVF result and the mtDNA copy number in the female genome, whereas some studies have corroborated the relationship [18, 19, 20, 21, 22, 23].
The aim of this study was to investigate the possibility that a successful pregnancy could be predicted in IVF patients by analyzing the contents of three genome repeats in DNA derived from blood leukocytes—ribosomal, telomere, and mitochondrial—as well as monitoring the changes in their contents before and after ovarian stimulation (preparation for IVF).
A total of 40 women (aged 28 to 40) were observed by embryologists for the IVF (In Vitro Fertilization) procedure. Among them, the IVF procedures ended with failure in 23 patients (IVF– group), while in 17 patients, the procedures ended in pregnancy (IVF+ group). Exclusion criteria from the study were (1) aged over 40 years; (2) diseases or pathological conditions accompanied by impaired fertility and/or reproductive function; (3) severe chronic diseases. Blood sampling from women who underwent IVF was conducted prior to beginning hormone therapy (sample 1) and 2–4 weeks later, immediately before the IVF procedure (sample 2).
The control group consisted of 17 non-pregnant women (aged 28 to 40), who were prophylactically observed in the same clinic and had no previous problems with conception or pregnancy. Blood sampling was performed on women in the control group only once in the first 10 days of their menstrual cycle.
Venous blood was collected into test tubes containing heparin. Erythrocyte lysis was performed, and leukocyte sediment was isolated. A solution containing 0.04 M EDTA, 2% sodium lauryl sarcosylate, and 150 µg/mL of RNase A (Sigma, St. Louis, Missouri 68178, USA) was added to the precipitate for 45 minutes at 37 °C. Then, the cell lysate was treated with proteinase K (200 µg/mL, Promega, Madison, Wisconsin 53711, USA) for 24 hours at 37 °C. Protein extraction was carried out using equal volumes of a mixture of phenol/chloroform/isoamyl alcohol (25:24:1), phenol, and a mixture of chloroform/isoamyl alcohol (24:1). Phenol was stabilized with 8-hydroxyquinoline. DNA was precipitated by adding 1/10 volume of 3 M sodium acetate (pH 5.2) and 2.5 volumes of ice-cold ethyl alcohol. DNA was collected by centrifugation for 15 minutes at 4 °C and washed with 70% ethanol before being dried and dissolved in water.
The repeat contents in the genomic DNA were assessed using the nonradioactive
quantitative hybridization (NQH) technique, which has previously been described
in detail [5, 7, 24]. Membranes (Optitran BA-S85, GE Healthcare) were soaked in a
20-fold SSC buffer pH 7.0. The DNA samples with a concentration of
~75 ng/µL were denatured in 0.1 M NaOH solution for
10 min. Then, the solution was neutralized by adding an equal volume of 20-fold
SSC buffer at pH 4.0. The denatured DNA samples were applied to activated
membranes in the form of dots (1.5 µL) in 4 repetitions. The
membranes were dried at 80 °C for 1.5 hours. Prehybridization was
carried out in a buffer (5-fold SSC buffer, 5 mg/mL BSA, 5 mg/mL ficoll, 5 mg/mL
polyvinylpyrrolidone, 0.5% SDS, 50% formamide) at 45 °C for 20
min. Hybridization was carried out under the same conditions in the buffer with a
biotin-labeled DNA probe. Nonspecific sorption was blocked using a buffer (0.05 M
Tris–HCl, pH 7.5, 0.5 M NaCl, 0.1% milk powder, 0.1% gelatin) for 10 minutes
with constant stirring. Then, streptavidin conjugate with alkaline phosphatase
(Merck) was added to the buffer “AP 7.5 oligo” (0.05 M Tris–HCl, 0.5 M NaCl,
0.05 M MgCl
DNA probes were labeled with biotin using nick translation. The pBR322–rDNA plasmid was used for rDNA analysis. The pBR322 vector contains an insert of the rDNA EcoR1 fragment (5836 b, between –515 and 5321 nucleotides from the transcription start site, “GenBank” No. U13369). Ribosomal DNA content is represented as the number of rDNA copies per diploid genome. Biotinylated mtDNA was used for mtDNA quantification. MtDNA content is represented as a number of mtDNA copies per diploid genome. TR quantification was performed using biotin-labeled oligonucleotide biotin-(TTAGGG) 7 (synthesized by Syntol, RF). TR content is presented in units of mass: pg TR/µg DNA.
Statistical data processing was performed using Microsoft Office Excel,
StatGraph, and StatPlus2007 Professional software. Nonparametric Mann–Whitney
(U-test) and Kolmogorov–Smirnov tests were applied. The differences were
considered statistically significant at p
The contents of the three genome repeats (ribosomal, mitochondrial, and telomeric) were determined from the blood leukocytes in each group. Then, the IVF procedure was performed. Next, the women were divided into two groups, according to the outcome: IVF+ (pregnancy occurred, N = 17) and IVF– (negative result, N = 23). Of the 17 women who underwent a successful IVF, 7 gave birth to twins. Fig. 1 shows the abundance of the three repeats in the genomes from leukocytes in the three groups: IVF+, IVF–, and C. Groups were compared by nonparametric statistics (U-test). Table 1 represents the descriptive statistics for the parameters studied in each group.
Comparisons of the groups by the contents of the three repeats: rDNA (A), TR (B), and mtDNA (C) in the genomic DNA; Mann–Whitney test.
|rDNA, copy number||IVF||40||388||113||215–673||374||0.29|
|mtDNA, copy number||IVF||40||345||236||40–917||253||0.69|
|TR, pg/µg, DNA||IVF||40||259||110||68–462||252||0.43|
The IVF+ and IVF– groups differed significantly (p
The IVF+ group differed significantly from the other three groups (p
In the groups of women who experienced problems with conception and had
resorted to IVF, we found significantly greater telomere repeat variability than
in the control C group. The maximum mean values of the TR content were observed
in the IVF+ group (306
Fig. 2A–C shows the pairwise dependencies linking the analyzed parameters. The table in Fig. 2D shows the results of the correlation analysis (only significant indicators are shown). For the IVF+ group, we found a positive correlation between the number of mtDNA copies and the TR content. Indeed, the highest was the TR content, whereby the longer the TRs were, the more mitochondria the cell contained. A similar correlation was observed across the entire sample.
An interesting trend was revealed (Fig. 2A–C), whereby for the 10 women in the IVF+ group with the highest mtDNA copy counts (more than 500) and elevated TR contents (more than 240 pg/µg of DNA), 7 gave birth to twins. Additionally, the rDNA CN ranged from 438 to 576 in these women.
Correlations between the three repeat copy numbers. (A–C) Pairwise correlations between the three repeat contents. The dotted lines surround an area of points corresponding to the IVF– group. (D) Three DNA repeat correlation analyses.
The IVF procedure was performed 2–4 weeks after the first blood sampling. Preparation for the IVF procedure included hormonal ovarian stimulation. Fig. 3A shows data reflecting the changes in the three genome repeat contents during ovarian stimulation. For the IVF+ and IVF– groups, significant changes were observed in the TR and mtDNA contents. The content of these repeats in the DNA of some participants increased, while in others, it did not change or even decrease. For the IVF+ group, a negative correlation was observed between the initial repeat contents in the DNA (sample 1) and the relative content change after preparation for IVF (sample 2). During the observation period, the content of the repeats in the leukocytes increased in the women who had initially contained low repeat copy numbers. In the sample leukocytes with higher TRs and mtDNA contents (1), a decrease was observed in the repeat copy numbers in the sample (2) (Fig. 3B,C).
Changes in the three genome repeat copy numbers. (A) Changes in the three genome repeat contents in women undergoing IVF during a 2–4-week preparation period. Blood samples (1) were taken during the first inspection. Blood samples (2) were also taken shortly before the IVF procedure. (B) The dependence of mtDNA content changes on the mtDNA content in the sample (1). (C) The dependence of TR changes on the TR content in the sample (1). The data analyses on correlation are presented in the boxes.
The ribosomal repeats in samples (1) and (2) varied in number, yet were within the experimental error, i.e., remained stable, suggesting it is a constant trait over a lifetime.
In the present study, we analyzed, for the first time, the copy numbers of the three repeats in the female genome in blood leukocytes. All three repeats perform essential cellular functions, such as modulating the level of protein synthesis (rDNA), controlling the biological age (TR), and supplying energy to the cell (mtDNA). A distinctive feature of the female genomes in women who had not been able to conceive naturally and had resorted to IVF was a higher variability in the contents of the three repeats compared to the control group (Table 1).
The maximum differences (p
The IVF+ and IVF– groups also differed in mtDNA CN (copy number), although the
differences were less pronounced (p
Previously, a number of authors have explored the possibility of using telomere
repeat numbers to predict the IVF outcome since failed IVF attempts were found to
be associated with a lower TR number in the genome [11, 12, 14, 15]. However, these
results contradict the data presented in the previous study . In our study,
we found no significant differences according to TR content in the IVF+ and IVF–
An important requirement for an index to reflect the chances of a successful IVF treatment is its stability over long periods of a lifespan. The TR and mtDNA contents cannot be deemed as stable genetic traits. Cells are heterogeneous in the abundance of these repeats over the same cell culture. The amount of mtDNA can increase or decrease under conditions of stress of varying intensities [16, 17]. The TR amount in blood cells is also influenced by psycho-emotional (affective) stress . Oxidative stress reduces the TR content in some cells in the cell pool, which manifests as a TR count decrease in the DNA isolated from the cells. However, the exaggerated clearance of cells with short telomeres via apoptosis following the induction of the stress response impacts the results of an increase in TR abundance.
We found high variability in the TR and mtDNA contents in the DNA isolated from the leukocytes of the same women before and after (2–4 weeks) preparation for the IVF procedure (Fig. 3). Preparation for IVF had a modulating effect on the contents of these repeats in DNA and was most pronounced in the IVF+ group. Moreover, changes in the TR and mtDNA repeat copy numbers can be attributed to both hormonal ovarian stimulation and the affective stress experienced by the women during this difficult period.
Unlike TR and mtDNA, the genomic ribosomal repeat content was stable and did not depend on the hormonal preparation for IVF (Fig. 3). Previously, we obtained the same result (variations in the TR content, albeit at a stable number of rDNA copies) in blood leukocyte DNA isolated from medical students, when studying the response to the affective stress caused by passing their exams . The study of replicative aging of cultured skin fibroblasts also showed that the rDNA content stability was in contrast to a decrease in the TR count [7, 31]. We found few cases of exception in aging cultures when analyzing late passages since these cells had lost hypermethylated rDNA copies, which were functionally inactive, and located outside the nucleolar structures, meaning they were a problem for replication. This process resulted in a reduction in the rDNA CN. However, the number of such cells in the total cell pool did not exceed 5–10% .
Our study has three limitations. Firstly, the sample size was restricted. Secondly, only Caucasian patients were involved. Finally, we explored the three repeats in the genomes of lymphocytes, without studying other tissues and organs, including reproductive. Thus, the question remains open as to what extent our findings can be extrapolated to additional tissues and organs and to different ethnicities. Our previous decades-long studies of phenotypic manifestations relating to the rDNA copy numbers suggest that the ribosomal genes, as universal “housekeeping” genes, manifest similarly in every cell and tissue (except for malignant), and their dosage effects do not depend on the racial origin. Moreover, we have already verified this assumption in recent studies on the different effects of rDNA copy numbers. For instance, we found elevated rDNA copy numbers in the lymphocytes of Caucasian patients with schizophrenia, a multifactorial disease, by still as yet obscure causes . Recently, these findings were corroborated by Sen Li and co-authors in East Asian (Japanese) patients, from a direct assay of their brain cells .
This pilot study has shown that the very low number of ribosomal repeat copies in the DNA isolated from blood leukocytes can be considered a stable and effective predictor of the failure of the IVF procedure. The combination of high mitochondrial DNA and telomere repeat content along with a moderately elevated ribosomal repeat number can be a potential predictor of multiple pregnancies. The major limitation of our pilot study was the use of a small sample size. Therefore, to validate these conclusions, further studies on larger samples and racially different populations (e.g., Asians) are warranted.
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
ESE, NNV and VLI designed the study; MPK, PAK and MAK selected and described patients and conducted blood sampling, ESE and EVK worked with cell cultures and performed DNA isolation; TAS determined the ribosomal repeat count and performed fluorescence microscopy; ONA determined the mtDNA and telomere repeat abundance; LNP and NNV wrote the text; LNP translated the text to English language; EVK and VLI statistically processed the data obtained; LNP, PEU, SIK and SVK made substantial contributions to data analysis and interpretation; SIK and SVK supervised the research. All authors contributed to editorial changes in the manuscript. All authors read and approved the final manuscript. All authors have participated sufficiently in the work and agreed to be accountable for all aspects of the work.
The study was carried out in accordance with the latest version of the Declaration of Helsinki and was approved by the Independent Interdisciplinary Ethics Committee on Ethical Review for Clinical Studies (Protocol No. 12 dated October 27, 2021). All participants provided written informed consent to participate in the study after the procedures had been completely explained.
FGFF-2022-0007/State assignment of the Ministry of Science and Higher Education (changes in TR and mtDNA repeat content) supported this research.
The authors declare no conflict of interest.
Publisher’s Note: IMR Press stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.