Ribosomal, Satellite III (1q12) and Telomere Repeat Copy Number Variations in Cystic Fibrosis Patients Used as a Model of Permanent Stress and Survivability

E.S. Ershova, E.I. Kondratyeva, L.N. Porokhovnik, A.Y. Voronkova, Yu. L. Melyanovskaya, S.A. Krasovsky, R.V. Veiko, E.K. Zhekaite, M.A. Starinova, N.N. Veiko, S.V. Kostyuk
doi: https://doi.org/10.1101/2024.04.21.24306126

ABSTRACT

Introduction Ribosomal (rDNA), satellite III (f-SatIII) and telomere (TR) tandem repeats perform a variety of functions in the human cell. Copy number variations (CNVs) of these repeats contributes to the global chromatin architecture and genome expression changes in response to stress or pathology. An elevated rDNA abundance and lowered contents of f-SatIII and TR were found in blood leukocytes of patients with schizophrenia. The main question of the study was: are the CNVs of the three repeat types in blood leukocytes a universal phenomenon linked to a patho4logy associated with chronic oxidative stress? Cystic fibrosis (CF), a monogenic disease was chosen as an object of the study.

Materials and Methods We determined the rDNA, f-SatIII and TR content in the blood leukocyte genomes of 545 subjects aged 0.2 to 50 years. The subjects were divided into three groups: Control (HC-group, N = 267), CF group (N= 186) and CF(d) group (severe patients, who died after some time upon sampling, aged 17 to 40 years, N = 92). For each patient, the type of the mutation in CFTR gene had been determined earlier. Non-radioactive quantitative hybridization technique was applied to quantify the repeats.

Results the rDNA abundance was elevated in the DNA of CF and CF(d) groups (565±105 copies per genome, N=278) compared to HC group (445±112 copies, N=267). A patient’s age 3 to 16 years was associated with “severe” mutations in CFTR (98% of cases), increased f-SatIII repeat counts and decreased telomere repeat (TR) contents in genome DNA compared to the age-matched controls. The genomes of deceased patients from CF(d) group also harbored increased numbers of f-SatIII and decreased numbers of TR. Patients above 16 years with a milder course of the disease and relatively low content of “severe” CFTR mutations contained less f-SatIII and more TR in their genomic DNA. A parameter rDNA·(f-SatIII/TR) showed a maximum difference between patients with relatively mild (age 17 to 40) and severe (age 17 to 40) forms of the pathology according to ROC analysis data (AUC = 0.86).

Conclusion Cystic fibrosis was associated with an increase in rDNA abundance and altered f-SatIII and TR contents in the DNA of cases compared to the controls. The severe course of the disease was characterized with high f-SatIII contents and shortened telomeres. Whereas mild CF cases were associated with low contents of f-SatIII and normal or slightly reduced telomere length. The index rDNA· (f-SatIII/TR) might be a predictor of the patient’s life expectancy.

1. INTRODUCTION

The human genome contains a significant number of repetitive sequences. Some of the repeats are clustered into tandem arrays. The tandem repeats play a principle regulation and structuration role in the functioning of the nuclear chromatin. Copy number variants (CNVs) of the tandem repeats have been implicated in many diseases. The tandem repeats demonstrate quantitative polymorphism, diversity in their nucleotide structure, localization on the chromosomes, their functions and stability in the processes of aging and stress response [1–11].

The ribosomal genes (rDNA) are a special type of tandem repeats. Human genome contains 100 to 900 rDNA copies per diploid cell [12–15]. A 44838 nucleotide pair long repeated unit contains a transcribed region, which codes for 47S rRNA, and a non-transcribed intergenic spacer, that includes non-extended repeats – microsatellites and transposons [15, 16]. Processing of 47S rRNA provides RNA synthesis (18S, 5.8S, 28S rRNA). The ribosomal tandem repeats are localized in p-arms of the five acrocentric chromosomes [17]. During interphase, the ribosomal repeats form a special structure termed nucleolus (Fig.1).

Figure 1.
Figure 1. The layout of human genome’s tandem repeats studied. The localization of the repeats in the interphase nucleus is indicated on the base of previous studies [18, 19]. The DNA probes used for hybridization are shown.

The functions of the nucleoli are not limited only to the production of ribosomal subunits [20–22]. The nucleolus is a complex hub, where the processes of ribosome biosynthesis, the cell cycle progress and the response to various types of stresses and viral infections are coordinated. Recent studies have shown that epigenetic status of the ribosomal genes and integrity of the nucleolar structure can modulate the cell homeostasis [23–28]. The discovery of structural and functional links between the nucleolus and the rest of the cell genome has shown that the nucleolus plays a key role in the organization of the nucleus architecture [29]. Copy number variations of the ribosomal repeats change the cell response to DNA damage [30].

A large amount of rDNA in the nucleolus stabilizes the nuclear heterochromatin by blocking the heterochromatin to euchromatin transition [31, 32]. The heterochromatin includes satellite repeats, or satellites. The satellites are tandemly repeated sequences organized into extended arrays and located in the centromeric and pericentromeric regions of the chromosomes. The number of the satellite DNA arrays can vary several times in individual chromatin regions [33]. The satellites play a structural and regulatory role in organizing the chromatin in the nucleus. The critical role of the pericentromeric repeats in maintaining chromosomal integrity during aging has been established [9, 34].

Previously, we studied a fragment of the pericentromeric satellites III (f-SatIII) using the DNA probe PUC1.77 [35] specifically hybridizing to the 1q12 pericentromeric heterochromatin of the first chromosome (Fig.1). A distinctive feature of the f-SatIII is the instability of its content in the DNA of the cell pool under conditions of oxidative stress and aging. Replicative senescence is accompanied by the accumulation of f-SatIII in a fraction of the cells [36]. The cells with high counts of the repeat are unable to divide and do not respond to stress by developing an adaptive response. These cells die when the stress levels rise [18, 37–40]. Transcription of satellite is one of the eventual mechanisms leading to the amplification of this repeat in some cells of a senescent population [41, 42].

The amplification of f-SatIII in senescent cells negatively correlates with a content of the tandem repeat (TTAGGG)n, which forms the telomeres [18, 43]. The telomere regions protect the chromosome ends against agglutination or base pair loss and maintain the stability of the cell cycles [44]. Unlike euchromatin sites of the genomic DNA, repair of single breaks is poorly developed in the telomeres’ DNA [45]. The two major signs of human physiologic aging, namely, loss of genome integrity and telomere shortening, result in disruption of the core architecture and replicative senescence [46]. The telomere shortening is associated with decondensation and transcription of heterochromatin sites [47–52].

We found early an interesting fact: DNA derived from schizophrenia cases contained low copy numbers of f-SatIII repeats together with an increased copy numbers of rDNA (rDNA CN) compared to the controls [43]. A decrease in the content of telomere repeats in patients with schizophrenia was also reported [53, 54]. Schizophrenia is a multifactorial disease that is accompanied by chronic oxidative stress [55]. We believe that oxidative stress is the cause of the observed CNVs of the tandem repeats in schizophrenia patients compared to healthy controls [43].

In order to understand how universal the profile of CNVs of the three repeats is under conditions of chronic oxidative stress induced by the pathology, we analyzed DNA samples isolated from blood leukocytes of patients with cystic fibrosis. Cystic fibrosis (CF) is an autosomal recessive disorder due to mutations in transmembrane conductance regulator (CFTR) gene leading to the abnormality of chloride channels in mucus and sweat producing cells. CFTR dysfunction causes the dysregulation of several signaling pathways to generate an innate oxidative state. An essential deficit of antioxidant molecules and an increase in oxidative stress have been shown in the patients with CF. The sustained imbalance between oxidant and antioxidant species induces chronic inflammation. The role of oxidative stress in the progression of lung injury in the CF patients has been described [56, 57].

Thus, a common feature of the two fundamentally different diseases (schizophrenia and cystic fibrosis) is systemic chronic oxidative stress, which could cause CNVs of the three tandem repeats in the genomes of blood leukocytes derived from patients compared to healthy controls. In order to test this assumption, we quantified the three repeats in DNA samples isolated from blood leukocytes of 267 healthy controls and 278 subjects with cystic fibrosis.

2. MATERIALS AND METHODS

2.1. Subjects

Totally, 278 cystic fibrosis cases aged 0.21 to 49 years were studied. Male to female ratio of the total sample was 1:0.96. The cases were divided into two groups. CF(d) group included 92 patients aged 16.8 up to 38.4 years, with a mean of 24.5±5.0 years, who died within several months after blood samples had been taken. CF group included 186 survived cases aged 0.21 up to 49.6 years, with a mean of 17.4±13.5 years.

CF was diagnosed according to the criteria of the CF clinical care recommendation and the CF national consensus guidelines (2019) [58]. For the assessment of patient’s status and description of the clinical picture, the Russian Federal CF Patient Register data were used. The Register data were collected from anamnesis and outpatient records of patients from Russian CF centers. The Register format corresponded to the European Register of CF patients [59]. The study design and informed consent form were approved by Ethics Committee of Federal State Budgetary Scientific Institution Research Centre for Medical Genetics named after academician N.P.Bochkov (RCMG). Each patient was under dispensary observation at the Research and Clinical Cystic Fibrosis Department of RCMG and was treated at the Cystic Fibrosis Departments of the Moscow Region Research Clinical Institute of Childhood and/or Pulmonology Scientific Research Institute under FMBA (Federal Medical and Biological Agency of the Russian Federation).

In the CF group patients, the following statistically significant correlations with the patient’s age were registered: (1) increasing age of diagnosis; (2) increasing frequency of chronic and intermittent infections caused by Pseudomonas aeruginosa, chronic infections by Burkholderia cepacia, Staphylococcus aureus, Achromobacter spp.; (3) decreasing respiratory function indices such as forced expiratory volume in 1 second (FEV1) and forced vital capacity (FVC). Complications such as allergic bronchopulmonary aspergillosis, diabetes mellitus, pneumothorax, hemoptysis, upper respiratory tract polyposis and osteoporosis prevailed in the groups of adult patients, especially among the deceased, compared to groups of children.

The microbiological status of the bronchopulmonary system was assessed using an algorithm for microbiological diagnosis of chronic lung infection in patients with CF, including application of bacteriological, biochemical, and other phenotypic tests (detection of hemolysis, test for biofilm formation capability) and molecular biology methods [58]. The biomaterials used in the medical inspection of the lower respiratory tract in patients with CF were sputum when coughing, a smear from the throat after coughing, laryngeal or nasopharyngeal aspirate, sputum induced by a hypertensive solution, bronchoalveolar lavage, brush biopsy material during bronchoscopy.

Testing for genetic variants of CFTR gene was conducted according to the standard algorithm [58, 60]. Most ‘mild’ [60] genotypes (34%) were found in the groups of adult patients over 17 compared to the groups of children (2.1%). In CF(d) group, the fraction of patients with ‘mild’ genotype (32%) was significantly (p <0.001) less than in the age-matched CF group (59%).

As healthy controls (HC group), DNA samples were obtained from the collection of the Laboratory of Molecular Biology of RCMG named after N.Bochkov. The dataset included 267 unrelated individuals (51% male) inhabiting Moscow city and the Moscow region. The HC group included healthy persons aged 5 to 50 with a mean age of 22.3 ± 11.5 years, who carried no mutation in CFTR gene and no other diagnosed mutation linked to genetic pathology.

The contents of f-SatIII and TR in human DNA are known to depend on the subject’s age and stress status [37, 40, 43]. In view of this, all patients were divided into five groups:

  1. CF(0-2): N=37, the age was from 0.21 to 2 years, the ‘severe’ genotype CFTR was detected in 97% of patients from this group;

  2. CF (3-16): N=54, the age from 3 to 16 years, the ‘severe’ genotype was detected in 98% of the patients;

  3. CF(d): N=92, the age from 17 to 40 years, the ‘severe’ genotype detected in 68% of the patients. All the cases included in the group died within several months after blood sampling, thus indicating a very severe course of the disease;

  4. CF(17-40): N=83, the age matched the age of patients from CF(d) group, the ‘severe’ genotype was detected in 58%;

  5. CF (41-50): N=12, age from 41 to 50, the ‘severe’ genotype detected in 40%.

The control cohort was divided into three groups, that matched the patient groups by age: HC (3-16), HC (17-40) and HC (41-50).

2.2. DNA Samples

5 mL of blood was collected from the peripheral vein using a syringe flushed with heparin (0.1 mL/5 mL blood) under strict aseptic conditions. Leukocytes were isolated from the blood after erythrocyte lysis. The DNA isolation method was described in detail previously (Chestkov et al. 2018). Briefly, 5 mL of the solution (2% sodium lauryl sarcosylate, 0.04 M EDTA), and 150 μg/mL RNAse A (Sigma, United States) were added to 0.5 mL of the leukocyte suspension (37°C, 45 min). The lysate samples were treated with proteinase K (200 μg /mL, Promega, United States) for 24 h (37°C). The samples were extracted with 5 mL of freshly distilled phenol (stabilized with 8-hydroxyquinoline), 5 mL of phenol/chloroform/isoamyl alcohol (25:24:1), and 5 mL chloroform/isoamyl alcohol (24:1). DNA was precipitated (0.3M sodium acetate, pH 5.2, and 2.5 volume of ethanol, −20°C) and collected by centrifugation (10,000 g, 15 min, 4°C), washed with 70% ethanol, then dried, and dissolved in water. The DNA concentration and purity were determined spectrophotometrically. The final DNA quantification was performed using PicoGreen dsDNA quantification reagent (Molecular Probes, Carlsbad, CA, United States). The DNA concentration in the samples was calculated according to a DNA standard curve. EnSpire equipment (Finland) with excitation and emission wavelengths of 488 and 528 nm was used.

2.3. Non-radioactive Quantitative Hybridization (NQH)

The method was described in detail in a previous publication [61]. Briefly, 2µL (20 ng/µL) of the denatured DNA samples (4-6 dots/sample) were applied to a filter (Optitran BA-S85, GE healthcare). 2µL of six standard samples of the genomic DNA (20 ng/µL) with a known content of the rDNA, f-SatIII and telomere repeats (TR) were applied to the same filter, in order to plot a calibration curve for the dependence of the signal intensity on the contents of the repeats in a particular sample. The filter was heated (80°C in vacuum, 1.5 h). After hybridization had been completed, the membrane filter was treated with a conjugate of streptavidin-alkaline phosphatase (Sigma) and substrates for alkaline phosphatase (NBT + BCIP). The filter was washed, dried, and scanned. In order to measure the hybridization signal, special software was used (“Imager 6,” RCMG, Moscow). The software determined the dot location and calculated the integral dot intensity. Signals from several dots corresponding to the same sample were averaged. The mean and standard error were calculated. The rDNA, f-SatIII and TR contents in the sample were calculated using a calibration curve equation. Relative standard error for NQH was merely 5±2%. Total relative standard error, that includes standard error for DNA extraction, concentration determination and NQH method, was 11±8%.

2.4. The DNA-Probes

Human rDNA pBR322-rDNA probe contains rDNA sequences (5836 bp) cloned into EcoRI site of pBR322 vector. The cloned rDNA fragment covers the positions from −515 to 5321 of the human rDNA (GenBank accession no. U13369), Fig. 1.

F-SatIII probe was a 1.77-kb cloned EcoRI fragment of human satellite DNA [35]. Dr. H. Cook (MRC, Edinburgh, UK) kindly supplied the human chromosome lql2-specific repetitive satellite DNA probe pUC1.77.

The DNA-probes pBR322-rDNA and pUC1.77 were biotinylated with biotin-11-dUTP using a nick translation kit (Biotin NT Labeling Kit, Jena Bioscience GmbH).

For the detection of the human telomere repeat, the probe biotin-(TTAGGG)7 was used. Syntol (Moscow, Russia) performed the synthesis and biotin labeling of the oligo-probe.

2.5. Statistical Analysis

From 4 to 6 parallel dots for one DNA sample were applied in one experiment. All reported results were reproduced at least two times as independent biological replicates. The descriptive statistics for rDNA, f-SatIII and TR are presented in Table 1. The significance of the observed differences was analyzed using non-parametric Mann-Whitney U tests and Kolmogorov– Smirnov statistics. Data were analyzed with StatPlus2007 Professional software (http://www. analystsoft.com/). All p-values were two-sided and considered statistically significant at p < 0.001. Descriptive statistics are also provided by Box Plots, which show the median, lower and upper quartiles, the minimum and maximum values, and outliers. MedCalc (https://www.medcalc.org/manual/roc-curves.php) created a complete sensitivity/specificity report using Receiver Operating Characteristic (ROC) curve analysis. Each point on the ROC curve represents a sensitivity/specificity pair corresponding to a particular decision threshold. The area under the ROC curve (AUC) reflects the parameter difference between two groups (disease/health or health/disease).

View this table:
Table 1. Descriptive statistics for the rDNA, f-SatIII and TR content in the DNA of cystic fibrosis patients (CF) and healthy controls (HC).

3. RESULTS

3.1. Ribosomal repeat CNVs

The cystic fibrosis affected cases differed from the controls by higher values of rDNA CN (Fig. 2Aa1 and Table 1). The patient sample contained no genome with rDNA CN < 350, while 18% of the controls counted less than 350 rDNA copies. The ROC-analysis data (Fig. 2Aa2) corroborated the fact of a significantly higher rDNA content in the genomes of patients with CF. The greatest differences from the controls (AUC=0.87) were observed for the group of deceased patients CF(d).

Figure 2.
Figure 2. Ribosomal DNA (rDNA) abundance in DNA of cystic fibrosis cases and healthy controls.

A. (a1) The rDNA copy numbers in DNA samples from HC, CF and CF(d) groups. A comparison between the groups is reported (U-test). (a2) ROC curves for the groups. The area under the ROC curve (AUC) is a measure of how well a parameter can distinguish between two groups.

B. The dependence of rDNA CN on subject’s age in HC, CF and CF(d) groups. The vertical dotted lines indicate dividing into four age groups.

C. The rDNA CN analysis in different age groups. Box Plot diagrams and group comparisons (U-test) are shown.

The rDNA CN did not depend on the age in the three groups (p>0.1, Fig.2B and Table 2). In each age group, rDNA CN in patient’s DNA were significantly higher than in the controls (Fig. 2C). The maximum coefficients of variation (C.V. = 0.26 – 0.22) and maximum values of rDNA CN (more than 900 copies) were found in groups of affected children and adolescents. In CF (17-40) and CF(d) groups, reduced intervals and coefficients of variation of rDNA CN were found (C.V.= 0.13 – 0.16). No observation of an rDNA content less than 450 and more than 800 copies occurred in patients over 30 years of age.

View this table:
Table 2. Spearman rank correlations for experimentally determined and calculated parameters

3.2. CNVs of f-SatIII repeats

Fig.3Aa1 and Table 1 represent data on f-SatIII repeat abundance in DNA isolated from blood leukocytes of all subjects studied. The f-SatIII contents in DNA samples that belong to HC-group are significantly higher than those in DNA of patients from CF-group. In CF(d) group, the repeat copy numbers in DNA are more than in HC-group (p < 0.005) and significantly exceed the values in CF group. The differences are substantiated by ROC analysis data (Fig.3Aa2).

Fig.3B shows the dependence of f-SatIII content in DNA of the study participants on their age. The f-SatIII contents of HC-group positively correlated with the subject’s age (Table 2, Rs = +0.49, p<10-7, N = 267). The leukocyte DNA of healthy children and adolescents under 16 contained low amounts of f-SatIII (14 ± 2 pg/ng DNA, N=93) compared to subjects at the age of 17 to 50 years (21 ± 6 pg/ng DNA, N=174). The contents of f-SatIII in total DNA samples of the whole CF cohort correlated with the patient’s age to a lesser extent (Table 2, Rs = +0.24, p<10-3, N = 186). In the subgroup of affected children 3 to 16 years of age, 23% of DNA samples carried more repeats, than any sample of the age-matched group of healthy children (Fig.3B).

Figure 3.
Figure 3. Satellite III (f-SatIII) abundance in DNA of cystic fibrosis cases and healthy controls.

A. (a1) The f-SatIII copy numbers in DNA samples from HC, CF and CF(d) groups. A comparison between the groups is reported (U-test). (a2) ROC curves for the groups.

B. The dependence of f-SatIII repeat copy numbers on subject’s age in HC, CF and CF(d) groups. The vertical dotted lines indicate dividing into four age groups.

C. The f-SatIII CN analysis in different age groups. Box Plot diagrams and group comparisons (U-test) are shown.

In figure 3C, one can see the data of f-SatIII quantification in DNA of the three groups (HC, CF and CF(d)) within the different age intervals. The maximum differences in the repeat content between patients with CF and healthy controls were found in age intervals over 16 years. In CF(17-40) and CF(41-50) groups of these age intervals, the coefficients of variation of the parameter were significantly reduced compared to the subgroups CF(d) and HC.

3.3. Telomere repeat CNVs

Fig.4Aa1 and Table 1 represent TR content data for HC, CF and CF(d) groups. The CF cohort differed from the CF(d) subgroup by a higher TR content. The TR content in DNA of cases was considerably lower than in the controls. The differences are corroborated by ROC-analysis results (Fig.4Aa2).

Fig.4B and Fig.4C show a dependence of TR content in DNA on the subject’s age. In HC group, a negative correlation of TR content in DNA with age was observed (Table 2, Rs = −0.32, p<0.0001, N = 267). The TR content in DNA of cases depended non-linearly on the patient’s age and could be described by an equation of the second degree (y = 0.18x2 – 6.8x + 346, R² = 0.2). The contents of the telomere repeats were maximum in the infant CF(0-2) group and considerably decreased in CF(3-16) group. An increase in TR content was observed in patients over 16 years of age. In patients who survived to an age of 40+, the TR content did not differ from the repeat content in the controls. The TR content in DNA samples of CF(d) group was lower than in DNA of age-matched CF cases (Fig.4C).

3.4. Comparison of CNVs of the three repeat types

Fig. 5 (A-C) represents the results of pairwise comparison of the contents of the studied repeats in cases and healthy controls. Table 2 shows the correlation dependencies. Fig.5D displays a three-dimensional plot showing the contents of the three repeats in DNA samples. In the control group, a weak positive correlation was found between the contents of rDNA and f-SatIII (Table 2, Rs =+0.19, p=0.002, N = 267). In the group of died patients with severe disease course CF(d), a positive correlation was found between the contents of rDNA and TR (Table 2, Rs =+0.43, p<0.0001, N = 92).

Figure 4.
Figure 4. Telomere repeat (TR) abundance in DNA of cystic fibrosis cases and healthy controls.

A. (a1) The TR copy numbers in DNA samples from HC, CF and CF(d) groups. A comparison between the groups is reported (U-test). (a2) ROC curves for the groups.

B. The dependence of TR copy numbers on subject’s age in HC, CF and CF(d) group. The vertical dotted lines indicate dividing into four age groups.

C. The TR CN analysis in different age groups. Box Plot diagrams and group comparisons (U-test) are shown.

Figure 5.
Figure 5. A. A plot of f-SatIII content as function of the content of rDNA. B. A plot of TR content as function of the content of rDNA. C. A plot of f-SatIII content as function of the content of TR. D. A plot of f-SatIII content as function of the content of rDNA (X axis) and TR (Y axis).

70% of the DNA samples from CF(d) group appeared within the following ranges of the three parameters: rDNA CN 500 to 650 copies, f-SatIII 6 to 35 pg/ng DNA, and TR 200 to 400 pg/μg DNA (Fig.5A-D). These ranges comprised merely 35% of DNA samples of the group CF and only 16% of samples of the group HC.

64% of DNA samples of HC group were grouped in the following area: 200 to 500 rDNA CN, 6 to 35 pg/ng DNA f-SatIII, and TR more than 300 pg/μg of DNA. This area comprised 19% of DNA samples of CF group and mere 3% of samples belonging to the most affected CF(d) group.

Thus, DNA samples of affected cases and healthy controls significantly differed by the combinations of the contents of the three repeats. The closest grouping of DNA samples in terms of the contents of the three repeats was characteristic of the CF(d) group with a severe course of the disease and a lethal outcome (Fig.5D).

3.5. Analysis of differences in the values of parameters between groups CF and CF(d)

The search for a prognostic marker associated to the disease severity was conducted by comparing parameters (indices) calculated for the patients from CF(d) group aged 17 to 40 years (N=92) and for the age-matched patients from CF(17-40) group (N=83). We analyzed the ratios f-SatIII/rDNA, TR/rDNA and f-SatIII /TR. Indicator rDNA*(f-SatIII/TR) simultaneously takes into account three parameters. The statistics showing the differences in experimental and calculated values of the parameters are listed in Table 3. The greatest differences between the analyzed subgroups CF and CF(d) were observed for the parameters f-SatIII /TR (AUC=0.82) and rDNA*f-SatIII/TR (AUC=0.86, cut-off – 32, specificity – 0.87 and sensitivity – 0.74).

View this table:
Table 3. Comparison of experimentally determined and calculated parameters for CF and CF(d) groups of patients of the same age (17 – 40 y.o.).]

Fig. 6 shows the values of parameters f-SatIII/TR and rDNA*f-SatIII/TR for the three subject groups. The maximum values were found in CF(d). In subgroup CF(3-16), the values were elevated as compared to subgroup CF(0.2-2). The same increase was also found in subgroup HC(17-40) compared to subgroup HC(3-16). In subgroup CF(17-40), index f-SatIII/TR and rDNA*f-SatIII/TR were considerably reduced as compared to subgroup CF (3-16). The analysis of distributions of this parameter corroborated significant differences between groups CF (3-16), CF(d) and the other groups (Fig.6C, data is provided for the index rDNA*f-SatIII/TR).

Figure 6.
Figure 6. The analysis of calculated parameters f-SatIII/TR (A) and rDNA*f-SatIII/TR (B).

(a1, b1) – Relationship between the corresponding parameter and age.

(a2, b2) – The analysis of parameter values in different age groups. Box Plot diagrams and group comparisons (U-test) are shown.

C. The analysis of rDNA*f-SatIII/TR index distributions in different age groups of CF cases using Kolmogorov-Smirnov test.

3.6. Changes in the copy numbers of the three repeats in patients DNA with aging

From the blood leukocytes of 22 patients aged 19 to 40 years, DNA isolation was performed twice with an interval of 5-8 years. DNA was analyzed simultaneously in two samples of blood leukocytes for each patient. Figure 7Aa1-a3 shows the copy numbers of the repeats in two DNA samples of the same patient. The parameter rDNA CN changed in neither patient with the course of time. The parameter f-SatIII increased in 10 patients and decreased in three cases. The parameter TR decreased in 8 cases, while increased in one case only.

Figure 7.
Figure 7. A. Changes of rDNA (a1), f-SatIII (a2), TR (a3) copy numbers and index rDNA*f-SatIII/TR in DNA samples isolated from blood of the same subjects at various time intervals (ΔT, indicated in the caption under X axis). (*)The difference is significant (p< 0.01). B. Distribution of the ratio of the parameter values in the second sample and the first sample obtained from the same subject. The maximum measurement error is indicated (gray box).

Thus, the copy numbers of telomeric and satellite repeats in the DNA of CF patients varied with aging, whereas the ribosomal repeat content remained stable. TR copy numbers tended to decrease, while f-SatIII content tended to increase in DNA with the course of time. Figure 7Aa4 also shows changes in the parameter rDNA*f-SatIII/TR for 22 patients. The dotted line shows the cut-off (32 units) determined via ROC analysis (Fig.6D). The index rDNA*f-SatIII/TR increased in 9 patients over time. In two patients, it reached values of 67 (#20) and 93 (#14), which significantly exceeded the cut-off of 32 units and mignt indicate aggravation in the condition of the patients.

4. DISCUSSION

We quantified three types of tandem repeats of the human genome (Fig.1A) in the same DNA samples isolated from the peripheral blood leukocytes of 278 patients with cystic fibrosis and 267 healthy controls. The studied repeats have different biological functions and characteristics indicating their abundance and localization in the nucleus, their CNVs under stress and in the course of human life or during replicative aging of a cell culture.

Our data have shown that rDNA CN is a stable genetic trait that is formed during the formation of a zygote and remains practically unchanged over the lifespan (Fig.7Aa1). The genomes of patients with cystic fibrosis harbored more rDNA copies compared to the controls (Fig.2). Probably, higher rDNA CN values are essential for survival in the presence of the genetic pathology in the genome of the fetus. The prenatal development with the hereditary pathology requires for an elevated protein synthesis (ribosome biogenesis) that only an elevated rDNA CN can provide [62]. In addition, a large enough number of rDNA copies in the nucleolus stabilizes the heterochromatin of the entire nucleus and reduces the risk of rearrangements under stress induced by a pathology or an external impact [31, 32]. The fact that rDNA copy numbers have been significantly less in DNA samples from soon-to-die patients compared to the survived cases is direct evidence of the selection for this trait resulting in better survivability of patients whose genomes initially carried more copies of ribosomal repeats. The variation in the number of ribosomal repeats over the population occurs during zygote formation and is determined by the combinatorics of the 10 pairs of acrocentric chromosomes in meiosis.

Unlike rDNA, the copy numbers of f-SatIII and TR repeats in human cells changes depending on aging and the endogenous and exogenous stress. Fig. 8 represents a scheme illustrating the changes in the contents of the repeats in primary and cultured human cells, which is based on our previous experimental findings [18, 37, 38, 42]. The cell pool is heterogeneous by the abundance of f-SatIII repeats. The fraction of the repeats in bulk DNA isolated from cells depends on the percentage of cells with different f-SatIII content. Most cells of a young body or a cell culture at early passages (type a, Fig.8) contain low f-SatIII amounts (Fig.3) and long telomeres (Fig.4). The index f-SatIII/ TR for the isolated DNA is minimum (Fig.6). The type a cells prevail in leukocyte pools of patients that belong to groups CF(0-2) and HC(3-16).

Figure 8.
Figure 8. A diagram illustrating changes in the contents of the three repeats in DNA isolated from a cell pool depending on age and/or stress. A, b, c and d are cell types that occur in the cell pool in various proportions depending on age and/or the level of oxidative stress. The diagram is based on the data obtained by the authors earlier.

During normal aging or under stress conditions, a fraction of cells with an elevated f-SatIII content increases in the cell pool [42]. The process that underpins the enlargement of f-SatIII abundance in normal (non-cancer) cells is still not described. For cancer cells, a copy gain scheme has been offered, that includes decondensation of heterochromatin, satellite repeat transcription, reverse transcription of the satellite RNA, and incorporation of DNA/RNA hybrid fragments into the nuclear DNA within pericentromeric loci. The authors suggest that cancer-associated derepression of specific repetitive (satellite) sequences can promote their RNA-driven genomic expansion [41]. The heterochromatin decondensation and transcription activation can promote the telomere shortening. During aging or stress, the cell pool enriches with type b and c cells (Fig.8).

An increase in f-SatIII content and a decrease in TR are observed in DNA isolated from these cells, while rDNA copy number remains stable. The type c cells do not respond to proliferative stimuli, and they are not capable of developing an adaptive response to a stress due to the chromatin rigidity [18]. These cells are eliminated from the pool when stress levels increase. If the fraction of these cells becomes too large and/or the process of their elimination is blocked, the cell population loses its ability to function normally.

In oldest cells or at a high chronic stress level, we observed a decrease in both f-SatIII and TR contents resulting in enrichment of the pool with type d cells (low contents of both telomere and satellite repeats in isolated DNA). The decrease in f-SatIII abundance in these cells might be caused by a block of the satellite transcription process. We showed earlier that the level of satellite III transcription depended non-linearly on the dose of the stress factor. Low doses activated the process, whereas high doses blocked it [18].

The cohort of CF patients was heterogeneous in respect of patient age and disease severity. The highest percentage of severe cases should be expected in CF(0-2) and CF(3-16) groups. The groups comprise 97-98% of “severe” CFTR mutation cases. The presence of “severe” mutations is known to be associated with an unfavorable prognosis and a shorter life expectancy for patients with cystic fibrosis [63, 64]. Older CF groups comprise fewer “severe” CFTR mutation cases. A part of “severe” CF(0-16) patients does not seem to survive and live to adulthood. The copy numbers of f-SatIII and TR depend on two factors: age and stress. The dysregulation of chloride channels results in an elevated level of the oxidative stress [55, 56]. Under the chronic pathology-driven stress the cells age faster. As a result, we observed an increase in f-SatIII/TR index by a factor of 1.5 to 3 in 40% of DNA samples from CF(3-16) group. In controls, this index did not increase before age 17 (Fig.6). Earlier, we observed the same changes in f-SatIII/ TR index during replicative senescence of cultured human skin fibroblasts [42]. Thus, accumulation of type c cells occurred in the population of blood leukocytes in 40% of patients over the age of two (Fig.8). The type c cells are characterized by a high f-SatIII abundance and short telomeres. Telomere shortening has been reported in blood cells of severe CF patients [65].

The comparison of the three age-matched groups [HC, CF(17-40) and CF(d)] has shown that a feature of survived patients older than 16 years with a relatively favorable course of the disease is a very low content of f-SatIII repeats in DNA compared to the age-matched controls (Fig.3) and deceased patients. We found earlier a similar effect of decreased copy numbers of f-SatIII repeats in DNA of some patients with schizophrenia. Like cystic fibrosis, schizophrenia is associated with chronic oxidative stress [54]. The small fraction of cells with an enormously elevated f-SatIII copy numbers in the population of blood leukocytes of the patients under chronic oxidative stress can be explained by several factors.

Firstly, oxidative stress stimulates the process of elimination of defective cells with an enormously high content of f-SatIII, which are not capable of developing an adaptive response to the stress [18]. Secondly, in the case of schizophrenia, we have shown that some compounds used in therapy are capable of suppressing the transcription processes that potentially result in an increase in the satellite content through the reverse transcription-driven copy gain mechanism described for cancer cells and hypothesized for normal cells. In addition, the compounds had stimulated the process of autophagy in cultured cells, which led to the elimination of cells with a high content of the satellite [37, 38]. Certain drugs used to treat complications of cystic fibrosis could also modulate the processes leading to similar changes in the satellite content of the cells of leukocyte pool. A larger fraction of the cells with elevated f-SatIII abundance in the group of deceased patients CF(d) suggests that too many of these cells are formed due to more active repeat transcription and/or disrupted processes of elimination of these cells from the bloodstream in these patients.

When comparing CF(17-40) and CF(d) groups, we found that an index to detect the differences between the groups in the best way is a parameter equal to the product of the f-SatIII/TR index and the content of rDNA (rDNA indicator). Large rDNA CNs are associated with a more severe genetic pathology as they can compensate for stress during prenatal and early postnatal life, while low rDNA CNs are eliminated. However, a severe genetic pathology results in a severe disease course in children. The index f-SatIII/TR reflects a degree of dysfunction in the cell pool due to the severe disease course induced by a severe mutation or other causes.

ROC-analysis determined the optimum cut-off for the index rDNA*f-SatIII/TR as 32 units (Fig.6). Within CF(d) group, the index exceeded this cut-off theshold in 74% of DNA samples. In CF(3-16) subgroup, a fraction of such DNA samples was 51%, while in CF (0-2) and CF(17-50) subgroups, such samples formed merely 9%. The indices rDNA*f-SatIII/TR and f-SatIII/TR might be used as a prognostic marker of the expected lifespan of a cystic fibrosis patient.

The processes that underpin the CNVs of f-SatIII repeats are of great importance for the essential body functions during aging and/or in pathology and require further research. The c type cells (Fig.8) are a ballast for the cell population because they do not respond to various stimuli, including those inducing proliferation or a normal adaptive response. The inability of the cells to respond to stimuli requiring alteration of the genome expression profile is associated with inability to change the structure of chromatin in response to the stimuli [18]. Perhaps, a therapy aimed at reducing the fraction of f-SatIII in the genomes of the cases would help improve the condition of the patients with the severe course of the disease.

5. CONCLUSION

Cystic fibrosis (CF) was associated with higher rDNA abundance and altered f-SatIII and TR contents in the DNA of cases compared to the controls. The severe course of the disease was characterized by high f-SatIII contents and shortened telomeres. The mild course of CF was associated with low contents of f-SatIII and normal or slightly reduced telomere length. The index rDNA*(f-SatIII/TR) might be a predictor of the patient’s life expectancy.

Data Availability

All data produced in the present study are available upon reasonable request to the authors

6. REFERENCES

  1. 1.↵
    Warburton PE, Hasson D, Guillem F, Lescale C, Jin X, Abrusan G. Analysis of the largest tandemly repeated DNA families in the human genome. BMC Genomics. 2008 Nov 7;9:533. doi: 10.1186/1471-2164-9-533.
    OpenUrlCrossRefPubMed
  2. 2.
    Chaisson MJP, Sulovari A, Valdmanis PN, Miller DE, Eichler EE. Advances in the discovery and analyses of human tandem repeats. Emerg Top Life Sci. 2023 Dec 14;7(3):361–381. doi: 10.1042/ETLS20230074.
    OpenUrlCrossRef
  3. 3.
    Miga KH. Centromeric Satellite DNAs: Hidden Sequence Variation in the Human Population. Genes (Basel). 2019 May 8;10(5):352. doi: 10.3390/genes10050352.
    OpenUrlCrossRef
  4. 4.
    Podgornaya OI, Ostromyshenskii DI, Enukashvily NI. Who Needs This Junk, or Genomic Dark Matter. Biochemistry (Mosc). 2018 Apr;83(4):450–466. doi: 10.1134/S0006297918040156.
    OpenUrlCrossRefPubMed
  5. 5.
    Thakur J, Packiaraj J, Henikoff S. Sequence, Chromatin and Evolution of Satellite DNA. Int J Mol Sci. 2021 Apr 21;22(9):4309. doi: 10.3390/ijms22094309.
    OpenUrlCrossRefPubMed
  6. 6.
    Ninomiya K, Yamazaki T, Hirose T. Satellite RNAs: emerging players in subnuclear architecture and gene regulation. EMBO J. 2023 Sep 18;42(18):e114331. doi: 10.15252/embj.2023114331.
    OpenUrlCrossRef
  7. 7.
    Copley KE, Shorter J. Repetitive elements in aging and neurodegeneration. Trends Genet. 2023 May;39(5):381–400. doi: 10.1016/j.tig.2023.02.008.
    OpenUrlCrossRef
  8. 8.
    Brahmachary M, Guilmatre A, Quilez J, Hasson D, Borel C, Warburton P, Sharp AJ. Digital genotyping of macrosatellites and multicopy genes reveals novel biological functions associated with copy number variation of large tandem repeats. PLoS Genet. 2014 Jun 19;10(6):e1004418. doi: 10.1371/journal.pgen.1004418.
    OpenUrlCrossRefPubMed
  9. 9.↵
    Dumbovic G, Forcales SV, Perucho M. Emerging roles of macrosatellite repeats in genome organization and disease development. Epigenetics. 2017 Jul 3;12(7):515–526. doi: 10.1080/15592294.2017.1318235.
    OpenUrlCrossRefPubMed
  10. 10.
    Lower SS, McGurk MP, Clark AG, Barbash DA. Satellite DNA evolution: old ideas, new approaches. Curr Opin Genet Dev. 2018 Apr;49:70–78. doi: 10.1016/j.gde.2018.03.003.
    OpenUrlCrossRefPubMed
  11. 11.↵
    Ugarković Đ, Sermek A, Ljubić S, Feliciello I. Satellite DNAs in Health and Disease. Genes (Basel). 2022 Jun 26;13(7):1154. doi: 10.3390/genes13071154.
    OpenUrlCrossRef
  12. 12.↵
    Parks MM, Kurylo CM, Dass RA, Bojmar L, Lyden D, Vincent CT, Blanchard SC. Variant ribosomal RNA alleles are conserved and exhibit tissue-specific expression. Sci Adv. 2018 Feb 28;4(2):eaao0665. doi: 10.1126/sciadv.aao0665.
    OpenUrlFREE Full Text
  13. 13.
    Agrawal S, Ganley AR. Complete Sequence Construction of the Highly Repetitive Ribosomal RNA Gene Repeats in Eukaryotes Using Whole Genome Sequence Data. Methods Mol Biol. 2016;1455:161–81. doi: 10.1007/978-1-4939-3792-9_13.
    OpenUrlCrossRefPubMed
  14. 14.
    Malinovskaya EM, Ershova ES, Golimbet VE, Porokhovnik LN, Lyapunova NA, Kutsev SI, Veiko NN, Kostyuk SV. Copy Number of Human Ribosomal Genes With Aging: Unchanged Mean, but Narrowed Range and Decreased Variance in Elderly Group. Front Genet. 2018 Aug 7;9:306. doi: 10.3389/fgene.2018.00306.
    OpenUrlCrossRef
  15. 15.↵
    Hori Y, Shimamoto A, Kobayashi T. The human ribosomal DNA array is composed of highly homogenized tandem clusters. Genome Res. 2021 Nov;31(11):1971–1982. doi: 10.1101/gr.275838.121.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    Gonzalez IL, Sylvester JE. Human rDNA: evolutionary patterns within the genes and tandem arrays derived from multiple chromosomes. Genomics. 2001 May 1;73(3):255–63. doi: 10.1006/geno.2001.6540.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    Stults DM, Killen MW, Pierce HH, Pierce AJ. Genomic architecture and inheritance of human ribosomal RNA gene clusters. Genome Res. 2008 Jan;18(1):13–8. doi: 10.1101/gr.6858507.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    Konkova MS, Ershova ES, Savinova EA, Malinovskaya EM, Shmarina GV, Martynov AV, Veiko RV, Zakharova NV, Umriukhin P, Kostyuk GP, Izhevskaya VL, Kutsev SI, Veiko NN, Kostyuk SV. 1Q12 Loci Movement in the Interphase Nucleus Under the Action of ROS Is an Important Component of the Mechanism That Determines Copy Number Variation of Satellite III (1q12) in Health and Schizophrenia. Front Cell Dev Biol. 2020 Jun 5;8:386. doi: 10.3389/fcell.2020.00386.
    OpenUrlCrossRef
  19. 19.↵
    Lai TP, Wright WE, Shay JW. Comparison of telomere length measurement methods. Philos Trans R Soc Lond B Biol Sci. 2018 Mar 5;373(1741):20160451. doi: 10.1098/rstb.2016.0451.
    OpenUrlCrossRefPubMed
  20. 20.↵
    Boisvert FM, van Koningsbruggen S, Navascués J, Lamond AI. The multifunctional nucleolus. Nat Rev Mol Cell Biol. 2007 Jul;8(7):574–85. doi: 10.1038/nrm2184.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.
    Pederson T. The plurifunctional nucleolus. Nucleic Acids Res. 1998 Sep 1;26(17):3871-6. doi: 10.1093/nar/26.17.3871.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    Pederson T, Tsai RY. In search of nonribosomal nucleolar protein function and regulation. J Cell Biol. 2009 Mar 23;184(6):771–6. doi: 10.1083/jcb.200812014.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    Ni C, Buszczak M. The homeostatic regulation of ribosome biogenesis. Semin Cell Dev Biol. 2023 Feb 28;136:13–26. doi: 10.1016/j.semcdb.2022.03.043.
    OpenUrlCrossRef
  24. 24.
    Rubbi CP, Milner J. Disruption of the nucleolus mediates stabilization of p53 in response to DNA damage and other stresses. EMBO J. 2003 Nov 17;22(22):6068–77. doi: 10.1093/emboj/cdg579.
    OpenUrlAbstract/FREE Full Text
  25. 25.
    Bursac S, Brdovcak MC, Donati G, Volarevic S. Activation of the tumor suppressor p53 upon impairment of ribosome biogenesis. Biochim Biophys Acta. 2014 Jun;1842(6):817–30. doi: 10.1016/j.bbadis.2013.08.014.
    OpenUrlCrossRefPubMed
  26. 26.
    Boulon S, Westman BJ, Hutten S, Boisvert FM, Lamond AI. The nucleolus under stress. Mol Cell. 2010 Oct 22;40(2):216–27. doi: 10.1016/j.molcel.2010.09.024.
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.
    Deisenroth C, Zhang Y. Ribosome biogenesis surveillance: probing the ribosomal protein-Mdm2-p53 pathway. Oncogene. 2010 Jul 29;29(30):4253–60. doi: 10.1038/onc.2010.189.
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.↵
    Olson MO. Sensing cellular stress: another new function for the nucleolus? Sci STKE. 2004 Mar 9;2004(224):pe10. doi: 10.1126/stke.2242004pe10.
    OpenUrlAbstract/FREE Full Text
  29. 29.↵
    Chubb JR, Boyle S, Perry P, Bickmore WA. Chromatin motion is constrained by association with nuclear compartments in human cells. Curr Biol. 2002 Mar 19;12(6):439–45. doi: 10.1016/s0960-9822(02)00695-4.
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    Kobayashi T. Regulation of ribosomal RNA gene copy number and its role in modulating genome integrity and evolutionary adaptability in yeast. Cell Mol Life Sci. 2011 Apr;68(8):1395–403. doi: 10.1007/s00018-010-0613-2.
    OpenUrlCrossRefPubMed
  31. 31.↵
    Paredes S, Maggert KA. Ribosomal DNA contributes to global chromatin regulation. Proc Natl Acad Sci U S A. 2009 Oct 20;106(42):17829–34. doi: 10.1073/pnas.0906811106.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    Paredes S, Branco AT, Hartl DL, Maggert KA, Lemos B. Ribosomal DNA deletions modulate genome-wide gene expression: “rDNA-sensitive” genes and natural variation. PLoS Genet. 2011 Apr;7(4):e1001376. doi: 10.1371/journal.pgen.1001376.
    OpenUrlCrossRefPubMed
  33. 33.↵
    Piotrowski A, Bruder CE, Andersson R, Diaz de Ståhl T, Menzel U, Sandgren J, Poplawski A, von Tell D, Crasto C, Bogdan A, Bartoszewski R, Bebok Z, Krzyzanowski M, Jankowski Z, Partridge EC, Komorowski J, Dumanski JP. Somatic mosaicism for copy number variation in differentiated human tissues. Hum Mutat. 2008 Sep;29(9):1118–24. doi: 10.1002/humu.20815.
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.↵
    Bierhoff H, Postepska-Igielska A, Grummt I. Noisy silence: non-coding RNA and heterochromatin formation at repetitive elements. Epigenetics. 2014 Jan;9(1):53–61. doi: 10.4161/epi.26485.
    OpenUrlCrossRefPubMed
  35. 35.↵
    Cooke HJ, Hindley J. Cloning of human satellite III DNA: different components are on different chromosomes. Nucleic Acids Res. 1979 Jul 25;6(10):3177–97. doi: 10.1093/nar/6.10.3177.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    Ershova ES, Malinovskaya EM, Konkova MS, Veiko RV, Umriukhin PE, Martynov AV, Kutsev SI, Veiko NN, Kostyuk SV. Copy Number Variation of Human Satellite III (1q12) With Aging. Front Genet. 2019 Aug 7;10:704. doi: 10.3389/fgene.2019.00704.
    OpenUrlCrossRef
  37. 37.↵
    Porokhovnik LN, Veiko NN, Ershova ES, Kostyuk SV. The Role of Human Satellite III (1q12) Copy Number Variation in the Adaptive Response during Aging, Stress, and Pathology: A Pendulum Model. Genes (Basel). 2021 Sep 28;12(10):1524. doi: 10.3390/genes12101524.
    OpenUrlCrossRef
  38. 38.↵
    Ershova ES, Savinova EA, Kameneva LV, Porokhovnik LN, Veiko RV, Salimova TA, Izhevskaya VL, Kutsev SI, Veiko NN, Kostyuk SV. Satellite III (1q12) Copy Number Variation in Cultured Human Skin Fibroblasts from Schizophrenic Patients and Healthy Controls. Front. Biosci. (Landmark Ed) 2023, 28(8), 191. doi:10.31083/j.fbl2808191.
    OpenUrlCrossRef
  39. 39.
    Ershova ES, Savinova EA, Kameneva LV, Porokhovnik LN, Veiko RV, Salimova TA, Izhevskaya VL, Kutsev SI, Veiko NN, Kostyuk SV. Antipsychotics Affect Satellite III (1q12) Copy Number Variations in the Cultured Human Skin Fibroblasts. Int J Mol Sci. 2023 Jul 10;24(14):11283. doi: 10.3390/ijms241411283.
    OpenUrlCrossRef
  40. 40.↵
    Umriukhin PE, Ershova ES, Filev AD, Agafonova ON, Martynov AV, Zakharova NV, Veiko RV, Porokhovnik LN, Kostyuk GP, Kutsev SI, Veiko NN, Kostyuk SV. The Psychoemotional Stress-Induced Changes in the Abundance of SatIII (1q12) and Telomere Repeats, but Not Ribosomal DNA, in Human Leukocytes. Genes (Basel). 2022 Feb 14;13(2):343. doi: 10.3390/genes13020343.
    OpenUrlCrossRef
  41. 41.↵
    Enukashvily NI, Donev R, Waisertreiger IS, Podgornaya OI. Human chromosome 1 satellite 3 DNA is decondensed, demethylated and transcribed in senescent cells and in A431 epithelial carcinoma cells. Cytogenet Genome Res. 2007;118(1):42–54. doi: 10.1159/000106440.
    OpenUrlCrossRefPubMed
  42. 42.↵
    Bersani F, Lee E, Kharchenko PV, Xu AW, Liu M, Xega K, MacKenzie OC, Brannigan BW, Wittner BS, Jung H, Ramaswamy S, Park PJ, Maheswaran S, Ting DT, Haber DA. Pericentromeric satellite repeat expansions through RNA-derived DNA intermediates in cancer. Proc Natl Acad Sci U S A. 2015 Dec 8;112(49):15148–53. doi: 10.1073/pnas.1518008112.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    Ershova ES, Malinovskaya EM, Golimbet VE, Lezheiko TV, Zakharova NV, Shmarina GV, Veiko RV, Umriukhin PE, Kostyuk GP, Kutsev SI, Izhevskaya VL, Veiko NN, Kostyuk SV. Copy number variations of satellite III (1q12) and ribosomal repeats in health and schizophrenia. Schizophr Res. 2020 Sep;223:199–212. doi: 10.1016/j.schres.2020.07.022.
    OpenUrlCrossRef
  44. 44.↵
    Harley CB. Telomere loss: mitotic clock or genetic time bomb? Mutat Res. 1991 Mar-Nov;256(2-6):271–82. doi: 10.1016/0921-8734(91)90018-7.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    Houben JM, Moonen HJ, van Schooten FJ, Hageman GJ. Telomere length assessment: biomarker of chronic oxidative stress? Free Radic Biol Med. 2008 Feb 1;44(3):235–46. doi: 10.1016/j.freeradbiomed.2007.10.001.
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    Sengupta D, Sengupta K. Lamin A and telomere maintenance in aging: Two to Tango. Mutat Res. 2022 Jul-Dec;825:111788. doi: 10.1016/j.mrfmmm.2022.111788.
    OpenUrlCrossRef
  47. 47.↵
    Toussaint O, Medrano EE, von Zglinicki T. Cellular and molecular mechanisms of stress-induced premature senescence (SIPS) of human diploid fibroblasts and melanocytes. Exp Gerontol. 2000 Oct;35(8):927–45. doi: 10.1016/s0531-5565(00)00180-7.
    OpenUrlCrossRefPubMedWeb of Science
  48. 48.
    Rizvi S, Raza ST, Mahdi F. Telomere length variations in aging and age-related diseases. Curr Aging Sci. 2014;7(3):161–7. doi: 10.2174/1874609808666150122153151.
    OpenUrlCrossRefPubMed
  49. 49.
    D’Mello MJ, Ross SA, Briel M, Anand SS, Gerstein H, Paré G. Association between shortened leukocyte telomere length and cardiometabolic outcomes: systematic review and meta-analysis. Circ Cardiovasc Genet. 2015 Feb;8(1):82–90. doi: 10.1161/CIRCGENETICS.113.000485.
    OpenUrlAbstract/FREE Full Text
  50. 50.
    Wong JY, De Vivo I, Lin X, Fang SC, Christiani DC. The relationship between inflammatory biomarkers and telomere length in an occupational prospective cohort study. PLoS One. 2014 Jan 27;9(1):e87348. doi: 10.1371/journal.pone.0087348. PMID: 24475279; PMCID: PMC3903646.
    OpenUrlCrossRefPubMed
  51. 51.
    Shen J, Terry MB, Liao Y, Gurvich I, Wang Q, Senie RT, Santella RM. Genetic variation in telomere maintenance genes, telomere length and breast cancer risk. PLoS One. 2012;7(9):e44308. doi: 10.1371/journal.pone.0044308.
    OpenUrlCrossRefPubMed
  52. 52.↵
    Barbé-Tuana FM, Parisi MM, Panizzutti BS, Fries GR, Grun LK, Guma FT, Kapczinski F, Berk M, Gama CS, Rosa AR. Shortened telomere length in bipolar disorder: a comparison of the early and late stages of disease. Braz J Psychiatry. 2016 Oct-Dec;38(4):281–286. doi: 10.1590/1516-4446-2016-1910.
    OpenUrlCrossRef
  53. 53.↵
    Solana C, Pereira D, Tarazona R. Early Senescence and Leukocyte Telomere Shortening in schizophrenia: A Role for Cytomegalovirus Infection? Brain Sci. 2018 Oct 18;8(10):188. doi: 10.3390/brainsci8100188.
    OpenUrlCrossRef
  54. 54.↵
    Nguyen TT, Eyler LT, Jeste DV. Systemic Biomarkers of Accelerated Aging in Schizophrenia: A Critical Review and Future Directions. Schizophr Bull. 2018 Feb 15;44(2):398–408. doi: 10.1093/schbul/sbx069.
    OpenUrlCrossRef
  55. 55.↵
    Ermakov EA, Dmitrieva EM, Parshukova DA, Kazantseva DV, Vasilieva AR, Smirnova LP. Oxidative Stress-Related Mechanisms in Schizophrenia Pathogenesis and New Treatment Perspectives. Oxid Med Cell Longev. 2021 Jan 23;2021:8881770. doi: 10.1155/2021/8881770.
    OpenUrlCrossRef
  56. 56.↵
    Causer AJ, Shute JK, Cummings MH, Shepherd AI, Gruet M, Costello JT, Bailey S, Lindley M, Pearson C, Connett G, Allenby MI, Carroll MP, Daniels T, Saynor ZL. Circulating biomarkers of antioxidant status and oxidative stress in people with cystic fibrosis: A systematic review and meta-analysis. Redox Biol. 2020 May;32:101436. doi: 10.1016/j.redox.2020.101436.
    OpenUrlCrossRef
  57. 57.↵
    Olveira C, Padilla A, Dorado A, Contreras V, Garcia-Fuentes E, Rubio-Martin E, Porras N, Doña E, Carmona A, Olveira G. Inflammation and Oxidation Biomarkers in Patients with Cystic Fibrosis: The Influence of Azithromycin. Eurasian J Med. 2017 Jun;49(2):118–123. doi: 10.5152/eurasianjmed.2017.17010.
    OpenUrlCrossRef
  58. 58.↵
    Kondratyeva E.I., Kashirskaya N.Yu., Kapranov N.I., eds. National consensus “Cystic fibrosis: definition, diagnostic criteria, therapy.” 2nd ed. M.: Borges; 2019. Available at: https://mukoviscidoz.org/doc/konsensus/2019/konsensus-2019-bez-rentgenogramm.pdf?ysclid=lc4lwl36qt378677355.
  59. 59.↵
    European Cystic Fibrosis Society. ECFS Patient Registry. [Update: June 20, 2022]. Available at: https://www.ecfs.eu/ecfspr
  60. 60.↵
    Cystic fibrosis genotype-phenotype consortium. Correlation between genotype and phenotype in patients with cystic fibrosis. Engl J Med 1993;329:1308.
    OpenUrl
  61. 61.↵
    Chestkov IV, Jestkova EM, Ershova ES, Golimbet VE, Lezheiko TV, Kolesina NY, Porokhovnik LN, Lyapunova NA, Izhevskaya VL, Kutsev SI, Veiko NN, Kostyuk SV. Abundance of ribosomal RNA gene copies in the genomes of schizophrenia patients. Schizophr Res. 2018 Jul;197:305–314. doi: 10.1016/j.schres.2018.01.001.
    OpenUrlCrossRefPubMed
  62. 62.↵
    Veiko NN, Ershova ES, Veiko RV, Umriukhin PE, Kurmyshev MV, Kostyuk GP, Kutsev SI, Kostyuk SV. Mild cognitive impairment is associated with low copy number of ribosomal genes in the genomes of elderly people. Front Genet. 2022 Sep 19;13:967448. doi: 10.3389/fgene.2022.967448.
    OpenUrlCrossRef
  63. 63.↵
    Mickle JE, Cutting GR. Genotype-phenotype relationships in cystic fibrosis. Med Clin North Am. 2000 May;84(3):597–607. doi: 10.1016/s0025-7125(05)70243-1.
    OpenUrlCrossRefPubMedWeb of Science
  64. 64.↵
    Bareil C, Bergougnoux A. CFTR gene variants, epidemiology and molecular pathology. Arch Pediatr. 2020 Feb;27 Suppl 1:eS8-eS12. doi: 10.1016/S0929-693X(20)30044-0.
    OpenUrlCrossRef
  65. 65.↵
    Glapa-Nowak A, Mutt SJ, Lisowska A, Sapiejka E, Goździk-Spychalska J, Wieczorek-Filipiak M, Drzymała-Czyż S, Nowak JK, Thalmann O, Herzig KH, Walkowiak J. Leukocyte Telomere Length Is Not Reduced in Children and Adults with Cystic Fibrosis but Associates with Clinical Characteristics - A Cross-Sectional Study. J Clin Med. 2021 Feb 4;10(4):590. doi: 10.3390/jcm10040590.
    OpenUrlCrossRef
Back to top
Posted April 22, 2024.
Download PDF
Data/Code
Email
Ribosomal, Satellite III (1q12) and Telomere Repeat Copy Number Variations in Cystic Fibrosis Patients Used as a Model of Permanent Stress and Survivability
E.S. Ershova, E.I. Kondratyeva, L.N. Porokhovnik, A.Y. Voronkova, Yu. L. Melyanovskaya, S.A. Krasovsky, R.V. Veiko, E.K. Zhekaite, M.A. Starinova, N.N. Veiko, S.V. Kostyuk
medRxiv 2024.04.21.24306126; doi: https://doi.org/10.1101/2024.04.21.24306126
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area