The Scientific World Journal

The Scientific World Journal / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 212183 |

Gordon S. K. Yau, Jacky W. Y. Lee, Victor T. Y. Tam, Stan Yip, Edith Cheng, Catherine C. L. Liu, Benjamin C. Y. Chu, Can Y. F. Yuen, "Differences in Risk Factors for Retinopathy of Prematurity Development in Paired Twins: A Chinese Population Study", The Scientific World Journal, vol. 2014, Article ID 212183, 7 pages, 2014.

Differences in Risk Factors for Retinopathy of Prematurity Development in Paired Twins: A Chinese Population Study

Academic Editor: Mary Elizabeth Hartnett
Received05 Jul 2014
Revised20 Sep 2014
Accepted20 Sep 2014
Published23 Oct 2014


Purpose. To determine the differences in risk factors for retinopathy of prematurity (ROP) in paired twins. Methods. A retrospective medical record review was performed for all paired twins screened for ROP between 2007 and 2012. Screening was offered to very low birth weight (≤1500 grams) and preterm (≤32 weeks) neonates. Twins 1 and 2 were categorized based on the order of delivery. Maternal and neonatal covariates were analyzed using univariate and multivariate regression analyses for both ROP and Type 1 ROP. Results. In 34 pairs of Chinese twins, the mean gestational age (GA) was 30.2 ± 2.0 weeks. In Twin 1, smaller GA (OR = 0.44, ), higher mean oxygen concentration (OR = 1.34, ), presence of thrombocytopenia (OR = 1429.60, ), and intraventricular hemorrhage (OR = 18.67, ) were significant risk factors for ROP. For Twin 2, a smaller GA (OR = 0.45, ) was the only risk factor. There were no significant risk factors for ROP in Twin 1 or Twin 2 on multivariate analysis. Conclusion. In Chinese twin pairs, smaller GA was the only common risk factor for ROP while Twin 1 was more susceptible to the postnatal risks for ROP.

1. Introduction

Retinopathy of prematurity (ROP) is a vasoproliferative retinal disease targeting low birth weight, preterm neonates [1]. ROP is one of the leading causes of childhood blindness worldwide [2]. Several risk factors have been demonstrated to be associated with the development of ROP including small gestational age (GA), low birth weight (BW), anemia, blood transfusion, mechanical ventilation especially with high concentration of supplementary oxygen (FiO2), hypoxemia, perinatal sepsis, use of inotropes, intraventricular hemorrhage (IVH), in vitro fertilization (IVF), and multiple pregnancies [37].

There are only a few publications in the literature investigating the risk factors for ROP development among twin pairs [810]. Paired, preterm twins share the same prenatal factors and have identical gestational ages; they differ only by their birth weight and postnatal exposures. Previous studies have primarily categorized paired twins according to their BW. In a study by Woo et al. [9], it was shown that there was no significant effect of BW on ROP development or on the severity of ROP among twins, while, in Azad et al.’s [8] series, the heavier twin was found to be at greater risk for a more severe stage of ROP. The purpose of this study was to determine the differences in risk factors for ROP development in paired twins that are categorized by their order of delivery.

2. Patients and Methods

The study was approved by the Institutional Review Board of the Hospital Authority of Hong Kong. The study was conducted in accordance with the Declaration of Helsinki and no patient personal data was disclosed in the study. The authors declare no financial or proprietary interests.

This was a retrospective study conducted at Caritas Medical Centre, Hong Kong, which provides ophthalmological service to 2 neonatal intensive care units (NICU) for a population of 1.9 million.

Medical records for consecutive newborns of paired twins that were screened for ROP between the period of January 2007 and December 2012 were retrieved using the Clinical Data Record System of the Hospital Authority of Hong Kong.

2.1. ROP Screening Criterion

All preterm babies admitted to these 2 NICU with a birth weight ≤ 1500 grams and/or gestational age ≤ 32 weeks were referred to a pediatric ophthalmologist for evaluation. All eligible preterms were examined according to the screening protocol recommended by the Royal College of Ophthalmologists and United Kingdom-ROP (UK-ROP) guidelines [11, 12]. Subjects were first screened at 4 to 8 weeks of postnatal age (≥30 weeks of GA) and were examined weekly to biweekly, until retinal vascularisation reached zone 3 or feature of established ROP regression [11]. Our institution has adopted the UK screening guideline since 1996 since, in Asian countries, severe ROP can be present even in more mature and heavier babies. Treatment was diode laser implemented when the disease progressed to Type 1 ROP as per the “Early treatment for retinopathy of prematurity” (ETROP) study [13]. The staging of ROP was recorded according to the revised International Classification of ROP, including the extent, zone, and presence or absence of “plus” disease [14].

All examinations were performed by 3 experienced pediatric ophthalmologists (Gordon S. K. Yau, Victor T. Y. Tam, and Benjamin C. Y. Chu). Each infant was screened by an indirect ophthalmoscope using a 30-dioptre (D) lens after full pharmacological pupil dilatation with tropicamide 0.5% and phenylephrine 1% eye drops. A lid speculum with scleral indentation after topical anesthesia (amethocaine) was routinely used. All screening was performed under oxygen saturation monitoring and the screening was temporarily withheld in case of desaturations.

The inclusion criteria included all preterm Chinese infants of paired twins irrespective of their ultimate development of ROP or Type 1 ROP. Neonates with incomplete clinical data, non-Chinese ethnicity, and those that deceased before the completion of ROP screening were excluded.

The primary outcome measures included the severity of ROP (the extent, zone, and presence or absence of “plus” disease) as well as the 34 risk factors (both maternal and neonatal) for the development of ROP as follows.

Antenatal maternal risk factors (Tables 1 and 2):(i)maternal diseases: preeclampsia (PET), gestational diabetes mellitus (GDM);(ii)IVF;(iii)use of antenatal steroid (ANS);(iv)preterm premature rupture of membrane.

Univariate logistic analysisMultivariate logistic analysis
Twin 1Twin 2Twin 1Twin 2
Odds ratioCoefficient estimates valueOdds ratioCoefficient estimates valueOdds ratioCoefficient estimates valueOdds ratioCoefficient estimates value

Gender (male versus female)−0.410.700.99−0.010.711.000.000.92
Gestational age0.44−0.820.02*0.45−0.790.03*
Birth weight
Preeclampsia 0.01−−3.931.
Gestational diabetes mellitus2.170.770.541.670.510.681.000.000.971.000.001.00
In vitro fertilization1.090.090.932.000.690.521.000.000.961.000.001.00
Postnatal hypotension3.371.220.360.02−3.931.
Inotropes use 7.001.950.209.672.
Antenatal steroid use 76.104.331.0050.513.921.
Apgar score 1 min0.75−
Apgar score 5 min0.71−0.340.301.160.150.751.000.000.961.000.001.00
Apgar score 10 min0.51−0.670.132.090.740.481.
Respiratory distress syndrome85.994.451.0053.823.991.
Bronchopulmonary dysplasia2.100.740.461.330.290.821.
Surfactant use85.994.451.0053.823.991.
Invasive mechanical ventilation4.711.550.132.330.850.431.000.000.951.000.001.00
Oxygen supplement0.14−1.950.2050.513.921.
Mean oxygen concentration1.340.290.03*0.95−0.050.721.000.000.991.000.000.99
Congenital heart disease1.560.450.722.000.690.521.000.000.991.000.000.99
Patent ductus arteriosus5.751.750.093.921.370.
NSAID use5.781.750.110.92−0.090.941.
Thrombocytopenia 1429.607.27<0.0001*0.02−3.931.
Blood transfusion10.502.350.050.92−0.090.941.000.000.991.000.001.00
Intraventricular hemorrhage 18.672.930.03*0.02−3.990.
Necrotizing colitis 0.01−4.310.099.672.
Neonatal jaundice1.270.240.840.25−1.390.
Total parenteral nutrition2.100.740.461.330.290.821.000.000.991.000.001.00
Hypoglycemia 0.83−0.180.8850.513.921.
Sepsis 0.96−0.040.970.78−0.250.841.
Preterm premature rupture of membrane0.41−0.890.452.000.690.521.
Acidosis 0.71−0.340.731.330.290.821.

Statistically significant.

Univariate logistic analysisMultivariate logistic analysis
Twin 1Twin 2Twin 1Twin 2
Odds ratioCoefficient estimates valueOdds ratioCoefficient estimates valueOdds ratioCoefficient estimates valueOdds ratioCoefficient estimates value

Gender (male versus female )31.513.451.000.07−2.641.
Gestational age0.26−1.330.160.70−0.361.−0.011.00
Birth weight−
Preeclampsia 0.04−−
Gestational diabetes mellitus0.05−2.951.0044.083.791.
In vitro fertilization19.392.961.0025.343.
Postnatal hypotension0.04−−
Inotrope use 0.03−3.511.00314.195.751.
Antenatal steroid use 17.292.851.0019.162.951.
Apgar score 1 min0.29−−0.170.711.
Apgar score 5 min0.49−0.710.230.42−0.860.
Apgar score 10 min0.68−0.380.640.36−1.020.441.000.000.971.000.000.94
Respiratory distress syndrome17.042.841.0016.622.811.
Bronchopulmonary dysplasia0.06−2.871.0058.744.
Surfactant use17.042.841.0016.622.811.
Invasive mechanical ventilation45.383.821.0033.643.521.
Oxygen supplement33.563.511.0019.162.951.
Mean oxygen concentration1.390.330.180.86−0.150.621.
Congenital heart disease63.094.141.0025.343.
Patent ductus arteriosus61.684.121.009.512.
NSAID use0.06−2.851.0061.684.
Thrombocytopenia 0.01−4.981.000.04−
Blood transfusion25.343.231.0061.684.
Intraventricular hemorrhage 161.315.081.000.03−3.511.
Necrotizing colitis 0.03−3.511.000.03−3.511.
Neonatal jaundice17.042.841.000.03−3.521.
Total parenteral nutrition0.06−2.871.0058.744.
Hypoglycemia 16.622.811.0019.162.951.
Sepsis 0.06−2.811.0045.383.821.
Preterm premature rupture of membrane0.05−3.091.0025.343.
Acidosis 0.03−3.521.000.06−2.811.

Neonatal risk factors (Tables 1 and 2):(i)demographic information (GA, BW, and gender);(ii)Apgar scores at 1, 5, and 10 minutes;(iii)postnatal interventions: surfactant administration; mechanical ventilation; use of supplementary oxygen; maintenance FiO2; use of nonsteroidal anti-inflammatory agents (NSAID) for patent ductus arteriosus (PDA) closure; blood transfusion; phototherapy; and total parental nutrition (TPN);(iv)neonatal diseases: respiratory distress syndrome (RDS); bronchopulmonary dysplasia; postnatal hypotension; congenital heart disease; PDA; anemia (defined as hemoglobin of <110 grams/litre, hematocrit < 25%); thrombocytopenia; neonatal jaundice (NNJ); IVH; necrotizing enterocolitis (NEC); hypoglycaemia; sepsis (culture positive or use of antibiotics for more than 7 days); meningitis; and acidosis based on blood Astrup results.

2.2. Statistics

Data were analyzed using standard univariate and multiple logistic regression analyses as implemented in the R Programming Language [15]. Most of the covariates were either binary or continuous. Estimates are also calculated independently by pregnancy order (Twin 1 or Twin 2) in order to compare the differences between their risk factors for ROP and Type 1 ROP. The responses and the nonresponses in the data were completely separated by predicators. Owing to the phenomenon of separation as described by Heinze and Schemper [16], the Penalized Maximum Likelihood Estimation method was implemented to reduce bias for related parameters as described by both Firth and Heinze et al. [17, 18]. A stepwise forward selection was used for choosing covariates for the multivariable analysis. Results were presented as odds ratios with the corresponding 95% confidence interval and values. value less than 0.05 was considered as statistically significant.

3. Results

During the study period, a total of 120 preterm infants of twin pregnancies were screened. Out of the 120 screened infants, 4 (3.3%) were of non-Chinese ethnicity, 72 infants were paired twins (60.0%), 1 infant (0.8%) did not survive before completion of ROP screening, and 1 (0.8%) had insufficient clinical information. The remaining 68 eligible infants of paired twins (34 pairs) were included for analysis.

3.1. Demographics

Of these 34 paired twin infants included in the study, all were of Chinese ethnicity. Eleven pairs were of different genders and 23 pairs were identical twins. Twenty-six infants (13 pairs) were the products of IVF and there was no case of twin-to-twin transfusion. Eighteen twin pairs underwent Caesarean section and 16 twin pairs had normal spontaneous deliver. There was no documented evidence of the use of instrumental delivery from the medical records. There were 39 male (57.4%) and 29 female (42.6%) subjects. The mean GA was 30.2 ± 2.0 weeks (range 25.28 to 35.0 weeks) and the mean BW was 1304.9 ± 278.6 grams (range 610.0 to 1950.0 grams). ROP of any stage developed in 9 (13.2%) infants and Type 1 ROP developed in 2 (2.9%) infants.

The mean GA of the twins was 30.2 ± 2.0 weeks (range 25.28 to 35.0 weeks). Among Twin 1’s, the mean BW was 1308.9 ± 279.6 grams (range 700.0 to 1950.0 grams). ROP of any stage developed in 5 (14.7%) infants and Type 1 ROP developed in 1 (2.9%) infant belonging to Twin 1. Among Twin 2’s, the mean BW was 1301.0 ± 281.8 grams (range 610.0 to 1720.0 grams). ROP of any stage developed in 4 infants (11.8%) and Type 1 ROP developed in 1 (2.9%) infant belonging to Twin 2. There was no significant difference in BW between the 2 twins (, -test).

3.2. Risk Factors for ROP

For ROP development in Twin 1, a younger GA (OR = 0.44, ), a higher FiO2 (OR = 1.34, ), presence of thrombocytopenia (OR = 1429.60, ), and IVH (OR = 18.67, ) were significant independent risk factors.

For Twin 2, a younger GA (OR = 0.45, ) was the only independent risk factor for ROP development and GA was also the only consistent independent risk factor for ROP in both Twin 1 and Twin 2. There were no significant dependent risk factors for ROP development in Twin 1 or Twin 2 on multivariate analysis (Table 1).

For Type 1 ROP, there were no significant independent or dependent risk factors found in Twin 1 or Twin 2 on both univariate and multivariate regression analyses (Table 2).

The relationship between gestational age and birth weight with ROP staging has been summarized in Figures 1 and 2.

4. Discussion

In univariate analysis, smaller GA was the only independent risk factor for ROP development irrespective of the order of delivery while the BW was not a significant risk factor. The significance of GA as a predictor for ROP was consistent with Woo et al.’s [9] study on discordant twin pairs where infants with GA ≤ 28 weeks were associated with significantly () higher risk of ROP development, stage 3 ROP, and threshold ROP requiring laser treatment (76.5%, 52.9%, and 29.4%, resp.) compared to those with GA > 28 weeks (15.8%, 0.0%, and 0.0%, resp.). On the other hand, in a twin pair study by Azad et al. [8], it was shown that there were no significant differences in both GA and BW ( for GA, for BW) among the twins with severe ROP (mean GA of 29.3 ± 1.3 weeks and mean BW of 1239 ± 206 grams) compared to those with less severe ROP (mean GA of 29.3 ± 1.3 weeks and mean BW of 1297 ± 188 grams).

Unlike in previous studies that categorized paired twins based on ROP severity or BW, we categorized the twins based on their order of delivery which has not been reported before in the literature. We feel that this classification is both clinically relevant for clinicians during their management of the twins and important for parents who may have questions about the risk of ROP in their twins. In our study, there were discrepancies in the susceptibility of postnatal risk factors for ROP between the twins despite exposure to similar postnatal conditions. Twin 1 was more likely to develop ROP when exposed to the following factors that were found to be independent ROP risk factors in Twin 1 but not Twin 2: presence of IVH, thrombocytopenia, and higher FiO2. These findings suggest that the first delivered infant of twin pregnancies may be more prone to develop ROP when exposed to the above risk factors as compared to the second infant. This prompts clinicians to be even more vigilant when screening for ROP in the first delivered twin, especially in the presence of the above risk factors.

A previous study by Prins [19] investigating the postnatal comorbidities among twins found that the second delivered infant was associated with higher risks of postnatal complications including intubation and RDS when compared to their twin. The author postulated that this observation was partly attributed by the smaller birth weight of the second twin rather than due to the order of delivery. Wadhawan et al. [20], on the other hand, found that there was no difference in neurodevelopmental impairment among twin pairs regardless of the order of delivery. In our study, there was no significant difference in BW between Twins 1 and 2 (); thus, the observed differences in ROP risk factors among the twins were most probably not related to BW. Our findings suggest that, in contrast to previous beliefs, Twin 2 was not necessarily the weaker twin and that Twin 1 should receive an equal amount if not a more vigilant monitoring and control for the risk factors leading to ROP. Larger scale prospective cotwin pair studies are warranted to provide an exact pathophysiological explanation to the differences in susceptibility to ROP among the twins.

To the best of our knowledge, this is the first study reporting the incidence and risk factors of ROP and Type 1 ROP in paired twins of Chinese infants based on the order of delivery. This study serves as a platform for future prospective and multicenter studies in this area as we anticipate that the number of preterm twins will likely continue to grow as a result of the growing demands of assisted reproduction.

Our study had its limitations. Firstly, the retrospective nature of this study inevitably generates inconsistencies in data although every effort was made to exclude subjects with incomplete clinical data. Secondly, subjects were screened by 3 pediatric ophthalmologists at different time points and minor interobserver variability can exits, which has not been objectively accounted for in this study but as all were trained to follow a strict ROP screening guideline and given the large population requiring screening, it was the most optimal balance in terms of providing clinical service and standardization for research. Lastly, there were only a small number of infants included in our study due to the rarity of preterm twins; with over 34 risk factors and only 34 pairs of twins, the distribution of risk factors among twin pairs were likely to be uneven and hence negative risk factors should be interpreted with care taking into account this limitation. Nevertheless, this study provides important data on the incidence and risk factors of ROP in the preterm Chinese infants from paired twins using more updated ROP screening guidelines than what currently exists in the literature.

5. Conclusion

In preterm, Chinese paired twins, a smaller GA was the only common independent risk factor for the development of ROP. The twin that was delivered first was more susceptible than the other twin to ROP development if they had IVH, thrombocytopenia, or higher FiO2 exposure.

Conflict of Interests

The authors declared no potential conflict of interests with respect to the research, authorship, and/or publication of this paper.


  1. T. L. Terry, “Retrolental fibroplasia in premature infants. Further studies on fibrolastic overgrowth of tunica vasculosa lentis,” Archives of Ophthalmology, vol. 33, pp. 203–208, 1945. View at: Google Scholar
  2. P. G. Steinkuller, L. Du, C. Gilbert, A. Foster, M. L. Collins, and D. K. Coats, “Childhood blindness,” Journal of American Association for Pediatric Ophthalmology and Strabismus, vol. 3, no. 1, pp. 26–32, 1999. View at: Publisher Site | Google Scholar
  3. E. M. Lad, T. C. Nguyen, J. M. Morton, and D. M. Moshfeghi, “Retinopathy of prematurity in the United States,” British Journal of Ophthalmology, vol. 92, no. 3, pp. 320–325, 2008. View at: Publisher Site | Google Scholar
  4. A. A. Binkhathlan, L. A. Almahmoud, M. J. Saleh, and S. Srungeri, “Retinopathy of prematurity in Saudi Arabia: incidence, risk factors, and the applicability of current screening criteria,” British Journal of Ophthalmology, vol. 92, no. 2, pp. 167–169, 2008. View at: Publisher Site | Google Scholar
  5. R. Karkhaneh, S. Z. Mousavi, M. Riazi-Esfahani et al., “Incidence and risk factors of retinopathy of prematurity in a tertiary eye hospital in Tehran,” British Journal of Ophthalmology, vol. 92, no. 11, pp. 1446–1449, 2008. View at: Publisher Site | Google Scholar
  6. M. B. Mizoguchi, T. G. Chu, F. M. Murphy, N. Willits, and L. S. Morse, “Dopamine use is an indicator for the development of threshold retinopathy of prematurity,” British Journal of Ophthalmology, vol. 83, no. 4, pp. 425–428, 1999. View at: Publisher Site | Google Scholar
  7. D. R. Brown, A. W. Biglan, and M. M. Stretavsky, “Retinopathy of prematurity: the relationship with intraventricular hemorrhage and bronchopulmonary dysplasia,” Journal of Pediatric Ophthalmology and Strabismus, vol. 27, no. 5, pp. 268–271, 1990. View at: Google Scholar
  8. R. Azad, P. Chandra, S. D. Patwardhan, and A. Gupta, “Profile of asymmetrical retinopathy of prematurity in twins,” Indian Journal of Ophthalmology, vol. 58, no. 3, pp. 209–211, 2010. View at: Publisher Site | Google Scholar
  9. S. J. Woo, K. H. Park, J. Ahn et al., “A co-twin study of the relative effect of birth weight and gestational age on retinopathy of prematurity,” Eye, vol. 25, no. 11, pp. 1478–1483, 2011. View at: Publisher Site | Google Scholar
  10. G. Sanghi, M. R. Dogra, S. Dutta, D. Katoch, and A. Gupta, “Intersibling variability of retinopathy of prematurity in twins and its risk factors,” International Ophthalmology, vol. 32, no. 2, pp. 113–117, 2012. View at: Publisher Site | Google Scholar
  11. Royal College of Paediatrics and Child Health, Guideline for the Screening and Treatment of Retinopathy of Prematurity, 2013,
  12. “Retinopathy of prematurity: guidelines for screening and treatment. The report of a Joint Working Party of the Royal College of Ophthalmologists and the British Association of Perinatal Medicine,” Early Human Development, vol. 46, no. 3, pp. 239–258, 1996. View at: Google Scholar
  13. Early Treatment For Retinopathy Of Prematurity Cooperative Group, “Revised indications for the treatment of retinopathy of prematurity: results of the early treatment for retinopathy of prematurity randomized trial,” Archives of Ophthalmology, vol. 121, no. 12, pp. 1684–1694, 2003. View at: Google Scholar
  14. International Committee for the Classification of Retinopathy of Prematurity, “The international classification of retinopathy of prematurity revisited: an international committee for the classification of retinopathy of prematurity,” Archives of Ophthalmology, vol. 123, no. 7, pp. 991–999, 2005. View at: Publisher Site | Google Scholar
  15. R Core Team, R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, 2014,
  16. G. Heinze and M. Schemper, “A solution to the problem of separation in logistic regression,” Statistics in Medicine, vol. 21, no. 16, pp. 2409–2419, 2002. View at: Publisher Site | Google Scholar
  17. D. Firth, “Bias reduction of maximum likelihood estimates,” Biometrika, vol. 80, no. 1, pp. 27–38, 1993. View at: Publisher Site | Google Scholar | MathSciNet
  18. G. Heinze, P. Meinhard, D. Dunkler, and H. Southworth, Logistf: Firth's bias reduced logistic regression, R package version 1.21, 2013,
  19. R. P. Prins, “The second-born twin: can we improve outcomes?” American Journal of Obstetrics and Gynecology, vol. 170, no. 6, pp. 1649–1657, 1994. View at: Publisher Site | Google Scholar
  20. R. Wadhawan, W. Oh, R. L. Perritt et al., “Twin gestation and neurodevelopmental outcome in extremely low birth weight infants,” Pediatrics, vol. 123, no. 2, pp. e220–e227, 2009. View at: Publisher Site | Google Scholar

Copyright © 2014 Gordon S. K. Yau et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.