Over the past decade, epigenetic analyses have made important contributions to our understanding of healthy development and a wide variety of adverse conditions such as cancer and psychopathology. There is increasing evidence that DNA methylation is a mechanism by which environmental factors influence gene transcription and, ultimately, phenotype. However, differentiating the effects of the environment from those of genetics on DNA methylation profiles remains a significant challenge. Monozygotic (MZ) twin study designs are unique in their ability to control for genetic differences because each pair of MZ twins shares essentially the same genetic sequence with the exception of a small number of de novo mutations and copy number variations. Thus, differences within twin pairs in gene expression and phenotype, including behavior, can be attributed in the majority of cases to environmental effects rather than genetic influence. In this article, we review the literature showing how MZ twin designs can be used to study basic epigenetic principles, contributing to understanding the role of early in utero and postnatal environmental factors on the development of psychopathology. We also highlight the importance of initiating longitudinal and experimental studies with MZ twins during pregnancy. This approach is especially important to identify: (1) critical time periods during which the early environment can impact brain and mental health development, and (2) the specific mechanisms through which early environmental effects may be mediated. These studies may inform the optimum timing and design for early preventive interventions aimed at reducing risk for psychopathology.

Low birth weight (LBW) can have an impact on health outcomes in later life, especially in relation to pre-disposition to metabolic disease. Several studies suggest that LBW resulting from restricted intrauterine growth leaves a footprint on DNA methylation in utero, and this influence likely persists into adulthood. To investigate this further, we performed epigenome-wide association analyses of blood DNA methylation using Infinium HumanMethylation450 BeadChip profiles in 71 adult monozygotic (MZ) twin pairs who were extremely discordant for birth weight. A signal mapping to the IGF1R gene (cg12562232, p = 2.62 × 10−8), was significantly associated with birth weight discordance at a genome-wide false-discovery rate (FDR) of 0.05. We pursued replication in three additional independent datasets of birth weight discordant MZ pairs and observed the same direction of association, but the results were not significant. However, a meta-analysis across the four independent samples, in total 216 birth-weight discordant MZ twin pairs, showed a significant positive association between birth weight and DNA methylation differences at IGF1R (random-effects meta-analysis p = .04), and the effect was particularly pronounced in older twins (random-effects meta-analysis p = .008, 98 older birth-weight discordant MZ twin pairs). The results suggest that severe intra-uterine growth differences (birth weight discordance >20%) are associated with methylation changes in the IGF1R gene in adulthood, independent of genetic effects.

The loss of estrogen during menopause causes changes in the female body, with wide-ranging effects on health. Estrogen-containing hormone replacement therapy (HRT) leads to a relief of typical menopausal symptoms, benefits bone and muscle health, and is associated with tissue-specific gene expression profiles. As gene expression is controlled by epigenetic factors (including DNA methylation), many of which are environmentally sensitive, it is plausible that at least part of the HRT-associated gene expression is due to changes in DNA methylation profile. We investigated genome-wide DNA methylation and gene expression patterns of white blood cells (WBCs) and their associations with body composition, including muscle and bone measures of monozygotic (MZ) female twin pairs discordant for HRT. We identified 7,855 nominally significant differentially methylated regions (DMRs) associated with 4,044 genes. Of the genes with DMRs, five (ACBA1, CCL5, FASLG, PPP2R2B, and UHRF1) were also differentially expressed. All have been previously associated with HRT or estrogenic regulation, but not with HRT-associated DNA methylation. All five genes were associated with bone mineral content (BMC), and ABCA1, FASLG, and UHRF1 were also associated with body adiposity. Our study is the first to show that HRT associates with genome-wide DNA methylation alterations in WBCs. Moreover, we show that five differentially expressed genes with DMRs associate with clinical measures, including body fat percentage, lean body mass, bone mass, and blood lipids. Our results indicate that at least part of the known beneficial HRT effects on body composition and bone mass may be regulated by DNA methylation associated alterations in gene expression in circulating WBCs.

Diurnal preference is an individual's preference for daily activities and sleep timing and is strongly correlated with the underlying circadian clock and the sleep-wake cycle validating its use as an indirect circadian measure in humans. Recent research has implicated DNA methylation as a mechanism involved in the regulation of the circadian clock system in humans and other mammals. In order to evaluate the extent of epigenetic differences associated with diurnal preference, we examined genome-wide patterns of DNA methylation in DNA from monozygotic (MZ) twin-pairs discordant for diurnal preference. MZ twins were selected from a longitudinal twin study designed to investigate the interplay of genetic and environmental factors in the development of emotional and behavioral difficulties. Fifteen pairs of MZ twins were identified in which one member scored considerably higher on the Horne–Ostberg Morningness–Eveningness Questionnaire (MEQ) than the other. Genome-wide DNA methylation patterns were assessed in twins’ buccal cell DNA using the Illumina Infinium HumanMethylation450 BeadChips. Quality control and data pre-processing was undertaken using the wateRmelon package. Differentially methylated probes (DMPs) were identified using an analysis strategy taking into account both the significance and the magnitude of DNA methylation differences. Our data indicate that DNA methylation differences are detectable in MZ twins discordant for diurnal preference. Moreover, downstream gene ontology (GO) enrichment analysis on the top-ranked diurnal preference associated DMPs revealed significant enrichment of pathways that have been previously associated with circadian rhythm regulation, including cell adhesion processes and calcium ion binding.

Identification of individuals within pairs of monozygotic (MZ) twins remains unresolved using common forensic DNA typing technology. For some criminal cases involving MZ twins as suspects, the twins had to be released due to inability to identify which of the pair was the perpetrator. In this study, we performed a genome-wide scan on whole blood-derived DNA from four pairs of healthy phenotypically concordant MZ twins using the methylated DNA immunoprecipitation sequencing technology to identify candidate DNA methylation markers with capacity to distinguish MZ twins within a pair. We identified 38 differential methylation regions showing within-pair methylation differences in all four MZ pairs. These are all located in CpG islands, 17 of which are promoter-associated, 17 are intergenic islands, and four are intragenic islands. Genes associated with these markers are related with cell proliferation, differentiation, and growth and development, including zinc finger proteins, PRRX2, RBBP9, or are involved in G-protein signaling, such as the regulator of G-protein signaling 16. Further validation studies on additional MZ twins are now required to evaluate the broader utility of these 38 markers for forensic use.

Monozygotic (MZ) twins provide a natural system for investigating developmental plasticity and the potential epigenetic origins of disease. A major difference in the intrauterine environment between MZ pairs is whether they share a common placenta or have separate placentas. Using DNA methylation measured at >400,000 points in the genome on the Illumina HumanMethylation450 array, we demonstrate that the co-twins of MZ pairs (average age of 14) that shared a common placenta (n = 18 pairs) have more similar DNA methylation levels in blood throughout the genome relative to those with separate placentas (n = 16 pairs). Functional annotation of the genomic regions that show significantly different correlation between monochorionic (MC) and dichorionic (DC) MZ pairs found an over-representation of genes involved in the regulation of transcription, neuronal development, and cellular differentiation. These results support the idea that prenatal environmental exposures may have a lasting effect on an individual's epigenetic landscape, and the potential for these changes to have functional consequences.

Aggressive behavior is highly heritable, while environmental influences, particularly early in life, are also important. Epigenetic mechanisms, such as DNA methylation, regulate gene expression throughout development and adulthood, and may mediate genetic and environmental effects on complex traits. We performed an epigenome-wide association study (EWAS) to identify regions in the genome where DNA methylation level is associated with aggressive behavior. Subjects took part in longitudinal survey studies from the Netherlands Twin Register (NTR) and participated in the NTR biobank project between 2004 and 2011 (N = 2,029, mean age at blood sampling = 36.4 years, SD = 12.4, females = 69.2%). Aggressive behavior was rated with the ASEBA Adult Self-Report (ASR). DNA methylation was measured in whole blood by the Illumina HM450k array. The association between aggressive behavior and DNA methylation level at 411,169 autosomal sites was tested. Association analyses in the entire cohort showed top sites at cg01792876 (chr8; 116,684,801, nearest gene = TRPS1, p = 7.6 × 10−7, False discovery rate (FDR) = 0.18) and cg06092953 (chr18; 77,905,699, nearest gene = PARD6G-AS1, p = 9.0 ×10−7, FDR = 0.18). Next, we compared methylation levels in 20 pairs of monozygotic (MZ) twins highly discordant for aggression. Here the top sites were cg21557159 (chr 11; 107,795,699, nearest gene = RAB39, p = 5.7 × 10−6, FDR = 0.99), cg08648367 (chr 19; 51,925,472, nearest gene = SIGLEC10, p = 7.6 × 10−6, FDR = 0.99), and cg14212412 (chr 6; 105,918,992, nearest gene = PREP, p = 8.0 × 10−6, FDR = 0.99). The two top hits based on the entire cohort showed the same direction of effect in discordant MZ pairs (cg01792876, Pdiscordant twins = 0.09 and cg06092953, Pdiscordant twins = 0.24). The other way around, two of the three most significant sites in discordant MZ pairs showed the same direction of effect in the entire cohort (cg08648367, Pentire EWAS = 0.59 and cg14212412, Pentire EWAS = 3.1 × 10−3). Gene ontology analysis highlighted significant enrichment of various central nervous system categories among higher-ranking methylation sites. Higher-ranking methylation sites also showed enrichment for DNase I hypersensitive sites and promoter regions, showing that DNA methylation in peripheral tissues is likely to be associated with aggressive behavior.

Tic disorders are moderately heritable common psychiatric disorders that can be highly troubling, both in childhood and in adulthood. In this study, we report results obtained in the first epigenome-wide association study (EWAS) of tic disorders. The subjects are participants in surveys at the Netherlands Twin Register (NTR) and the NTR biobank project. Tic disorders were measured with a self-report version of the Yale Global Tic Severity Scale Abbreviated version (YGTSS-ABBR), included in the 8th wave NTR data collection (2008). DNA methylation data consisted of 411,169 autosomal methylation sites assessed by the Illumina Infinium HumanMethylation450 BeadChip Kit (HM450k array). Phenotype and DNA methylation data were available in 1,678 subjects (mean age = 41.5). No probes reached genome-wide significance (p < 1.2 × 10−7). The strongest associated probe was cg15583738, located in an intergenic region on chromosome 8 (p = 1.98 × 10−6). Several of the top ranking probes (p < 1 × 10−4) were in or nearby genes previously associated with neurological disorders (e.g., GABBRI, BLM, and ADAM10), warranting their further investigation in relation to tic disorders. The top significantly enriched gene ontology (GO) terms among higher ranking methylation sites included anatomical structure morphogenesis (GO:0009653, p = 4.6 × 10−15) developmental process (GO:0032502, p = 2.96 × 10−12), and cellular developmental process (GO:0048869, p = 1.96 × 10−12). Overall, these results provide a first insight into the epigenetic mechanisms of tic disorders. This first study assesses the role of DNA methylation in tic disorders, and it lays the foundations for future work aiming to unravel the biological mechanisms underlying the architecture of this disorder.

Wellbeing (WB) is a major topic of research across several scientific disciplines, partly driven by its strong association with psychological and mental health. Twin-family studies have found that both genotype and environment play an important role in explaining the variance in WB. Epigenetic mechanisms, such as DNA methylation, regulate gene expression, and may mediate genetic and environmental effects on WB. Here, for the first time, we apply an epigenome-wide association study (EWAS) approach to identify differentially methylated sites associated with individual differences in WB. Subjects were part of the longitudinal survey studies of the Netherlands Twin Register (NTR) and participated in the NTR biobank project between 2002 and 2011. WB was assessed by a short inventory that measures satisfaction with life (SAT). DNA methylation was measured in whole blood by the Illumina Infinium HumanMethylation450 BeadChip (HM450k array) and the association between WB and DNA methylation level was tested at 411,169 autosomal sites. Two sites (cg10845147, p = 1.51 * 10−8 and cg01940273, p = 2.34 * 10−8) reached genome-wide significance following Bonferonni correction. Four more sites (cg03329539, p = 2.76* 10−7; cg09716613, p = 3.23 * 10−7; cg04387347, p = 3.95 * 10−7; and cg02290168, p = 5.23 * 10−7) were considered to be genome-wide significant when applying the widely used criterion of a FDR q value < 0.05. Gene ontology (GO) analysis highlighted enrichment of several central nervous system categories among higher-ranking methylation sites. Overall, these results provide a first insight into the epigenetic mechanisms associated with WB and lay the foundations for future work aiming to unravel the biological mechanisms underlying a complex trait like WB.

The disease- and mortality-related difference between biological age based on DNA methylation and chronological age (Δage) has been found to have approximately 40% heritability by assuming that the familial correlation is only explained by additive genetic factors. We calculated two different Δage measures for 132 middle-aged female twin pairs (66 monozygotic and 66 dizygotic twin pairs) and their 215 sisters using DNA methylation data measured by the Infinium HumanMethylation450 BeadChip arrays. For each Δage measure, and their combined measure, we estimated the familial correlation for MZ, DZ and sibling pairs using the multivariate normal model for pedigree analysis. We also pooled our estimates with those from a former study to estimate weighted average correlations. For both Δage measures, there was familial correlation that varied across different types of relatives. No evidence of a difference was found between the MZ and DZ pair correlations, or between the DZ and sibling pair correlations. The only difference was between the MZ and sibling pair correlations (p < .01), and there was marginal evidence that the MZ pair correlation was greater than twice the sibling pair correlation (p < .08). For weighted average correlation, there was evidence that the MZ pair correlation was greater than the DZ pair correlation (p < .03), and marginally greater than twice the sibling pair correlation (p < .08). The varied familial correlation of Δage is not explained by additive genetic factors alone, implying the existence of shared non-genetic factors explaining variation in Δage for middle-aged women.

Neurophysiological measurements of the response to pre-pulse and startle stimuli have been suggested to represent an important endophenotype for both substance dependence and other select psychiatric disorders. We have previously shown, in young adult Mexican Americans (MA), that presentation of a short delay acoustic pre-pulse, prior to the startle stimuli can elicit a late negative component at about 400 msec (N4S), in the event-related potential (ERP), recorded from frontal cortical areas. In the present study, we investigated whether genetic factors associated with this endophenotype could be identified. The study included 420 (age 18–30 years) MA men (n = 170), and women (n = 250). DNA was genotyped using an Affymetrix Axiom Exome1A chip. An association analysis revealed that the CCKAR and CCKBR (cholecystokinin A and B receptor) genes each had a nearby variant that showed suggestive significance with the amplitude of the N4S component to pre-pulse stimuli. The neurotransmitter cholecystokinin (CCK), along with its receptors, CCKAR and CCKBR, have been previously associated with psychiatric disorders, suggesting that variants near these genes may play a role in the pre-pulse/startle response in this cohort.

Structural brain magnetic resonance imaging (MRI) traits share part of their genetic variance with cognitive traits. Here, we use genetic association results from large meta-analytic studies of genome-wide association (GWA) for brain infarcts (BI), white matter hyperintensities, intracranial, hippocampal, and total brain volumes to estimate polygenic scores for these traits in three Scottish samples: Generation Scotland: Scottish Family Health Study (GS:SFHS), and the Lothian Birth Cohorts of 1936 (LBC1936) and 1921 (LBC1921). These five brain MRI trait polygenic scores were then used to: (1) predict corresponding MRI traits in the LBC1936 (numbers ranged 573 to 630 across traits), and (2) predict cognitive traits in all three cohorts (in 8,115–8,250 persons). In the LBC1936, all MRI phenotypic traits were correlated with at least one cognitive measure, and polygenic prediction of MRI traits was observed for intracranial volume. Meta-analysis of the correlations between MRI polygenic scores and cognitive traits revealed a significant negative correlation (maximal r = 0.08) between the HV polygenic score and measures of global cognitive ability collected in childhood and in old age in the Lothian Birth Cohorts. The lack of association to a related general cognitive measure when including the GS:SFHS points to either type 1 error or the importance of using prediction samples that closely match the demographics of the GWA samples from which prediction is based. Ideally, these analyses should be repeated in larger samples with data on both MRI and cognition, and using MRI GWA results from even larger meta-analysis studies.

Longitudinal studies of neuroticism have shown that, on average, neuroticism scores decrease from adolescence to adulthood. The heritability of neuroticism is estimated between 0.30 and 0.60 and does not seem to vary greatly as a function of age. Shared environmental effects are rarely reported. Less is known about the role of genetic and environmental influences on the rank order stability of neuroticism in the period from adolescence to adulthood. We studied the stability of neuroticism in a cohort sequential (classical) twin design, from adolescence (age 14 years) to young adulthood (age 32 years). A genetic simplex model that was fitted to the longitudinal neuroticism data showed that the genetic stability of neuroticism was relatively high (genetic correlations between adjacent age bins >0.9), and increased from adolescence to adulthood. Environmental stability was appreciably lower (environmental correlations between adjacent age bins were between 0.3 and 0.6). This low stability was largely due to age-specific environmental variance, which was dominated by measurement error. This attenuated the age-to-age environmental correlations. We constructed an environmental covariance matrix corrected for this error, under the strong assumption that all age-specific environmental variance is error variance. The environmental (co)variance matrix corrected for attenuation revealed highly stable environmental influences on neuroticism (correlations between adjacent age bins were between 0.7 and 0.9). Our results indicate that both genetic and environmental influences have enduring effects on individual differences in neuroticism.

Approximately 60% of children with reading difficulties (RD) meet criteria for at least one co-occurring disorder. The most common of these, attention deficit-hyperactivity disorder (ADHD), occurs in 20–40% of individuals with RD. Recent studies have suggested that genetic influences are responsible. To assess the genetic etiologies of RD and the comorbidity of RD and two ADHD symptom dimensions –– inattention (IN) and hyperactivity/impulsivity (H/I) –– we are conducting the first longitudinal twin study of RD and ADHD. Data from twin pairs in which at least one member of the pair met criteria for proband status for RD at initial assessment, and were reassessed 5 years later, were subjected to DeFries-Fulker (DF) analysis. Analyses of reading composite data indicated that over 60% of the proband deficit at initial assessment was due to genetic influences, and that reading deficits at follow-up were due substantially to the same genetic influences. When a bivariate DF model was fitted to reading performance and IN data, genetic influences accounted for 60% of contemporaneous comorbidity and over 60% of the longitudinal relationship. In contrast, analysis of the comorbidity between reading performance and H/I indicated that common genetic influences accounted for only about 20% of the contemporaneous and about 10% of the longitudinal relationships. Results indicate that (1) genetic influences on RD are substantial and highly stable; (2) the comorbidity between RD and IN is due largely to genetic influences, both contemporaneously and longitudinally; and (3) genetic influences contribute significantly less to the comorbidity between RD and H/I.

Twin pairs discordant for disease may help elucidate the epigenetic mechanisms and causal environmental factors in disease development and progression. To obtain the numbers of pairs, especially monozygotic (MZ) twin pairs, necessary for in-depth studies while also allowing for replication, twin studies worldwide need to pool their resources. The Discordant Twin (DISCOTWIN) consortium was established for this goal. Here, we describe the DISCOTWIN Consortium and present an analysis of type 2 diabetes (T2D) data in nearly 35,000 twin pairs. Seven twin cohorts from Europe (Denmark, Finland, Norway, the Netherlands, Spain, Sweden, and the United Kingdom) and one from Australia investigated the rate of discordance for T2D in same-sex twin pairs aged 45 years and older. Data were available for 34,166 same-sex twin pairs, of which 13,970 were MZ, with T2D diagnosis based on self-reported diagnosis and medication use, fasting glucose and insulin measures, or medical records. The prevalence of T2D ranged from 2.6% to 12.3% across the cohorts depending on age, body mass index (BMI), and national diabetes prevalence. T2D discordance rate was lower for MZ (5.1%, range 2.9–11.2%) than for same-sex dizygotic (DZ) (8.0%, range 4.9–13.5%) pairs. Across DISCOTWIN, 720 discordant MZ pairs were identified. Except for the oldest of the Danish cohorts (mean age 79), heritability estimates based on contingency tables were moderate to high (0.47–0.77). From a meta-analysis of all data, the heritability was estimated at 72% (95% confidence interval 61–78%). This study demonstrated high T2D prevalence and high heritability for T2D liability across twin cohorts. Therefore, the number of discordant MZ pairs for T2D is limited. By combining national resources, the DISCOTWIN Consortium maximizes the number of discordant MZ pairs needed for in-depth genotyping, multi-omics, and phenotyping studies, which may provide unique insights into the pathways linking genes to the development of many diseases.

An impressive literature has revealed that variation in virtually every measurable phenotype is the result of a combination of genetic and environmental influences. Based on these findings, studies that fail to use genetically informed modeling strategies risk model misspecification and biased parameter estimates. Twin- and adoption-based research designs have frequently been used to overcome this limitation. Despite the many advantages of such approaches, many available datasets do not contain samples of twins, siblings or adoptees, making it impossible to utilize these modeling strategies. The current study proposes a measurement strategy for estimating the intergenerational transmission of antisocial behavior (ASB) within a nationally representative sample of singletons using an extended pedigree risk approach that relies on information from first- and second-degree relatives. An evaluation of this approach revealed a pattern of findings that directly aligned with studies examining ASB using more traditional twin- and adoption-based research designs. While the proposed pedigree risk approach is not capable of effectively isolating genetic and environmental influences, this overall alignment in results provides tentative evidence suggesting that the proposed pedigree risk measure effectively captures genetic influences. Future replication studies are necessary as this observation remains preliminary. Whenever possible, more traditional quantitative genetic methodologies should be favored, but the presented strategy remains a viable alternative for more limited samples.

Age difference in the etiology of the relationship between childhood negative emotionality (NE) and conduct problems (CP) has not been previously investigated. Mothers of 662 pairs of twins completed questions on the emotionality (NE) scale of the EAS temperament survey and the CP scale of the Strengths and Difficulties Questionnaires (SDQ) via a telephone interview. Twin data were analyzed separately in younger (ages 3 to 7 years; mostly pre-schoolers) and older children (ages 8 to 13 years; mostly elementary school children). The phenotypic correlation between NE and CP increased from 0.33 among younger twins to 0.43 among older twins. Bivariate model-fitting analysis was performed to determine age difference in the etiology of the relationship between NE and CP. Among younger twins, the correlation between NE and CP was entirely explained by additive genetic factors common to NE and CP. Among older children, however, a small but significant amount of unique environmental correlation emerged to account for about 47% of the phenotypic correlation between NE and CP. The remaining 53% of the phenotypic correlation was due to shared additive genetic factors. We speculate that environmental factors associated with school adjustment may exert influences on the relationship between NE and CP among elementary school children.

Food liking-disliking patterns may strongly influence food choices and health. Here we assess: (1) whether food preference patterns are genetic/environmentally driven; and (2) the relationship between metabolomics profiles and food preference patterns in a large population of twins. 2,107 individuals from TwinsUK completed an online food and lifestyle preference questionnaire. Principle components analysis was undertaken to identify patterns of food liking-disliking. Heritability estimates for each liking pattern were obtained by structural equation modeling. The correlation between blood metabolomics profiles (280 metabolites) and each food liking pattern was assessed in a subset of 1,491 individuals and replicated in an independent subset of monozygotic twin pairs discordant for the liking pattern (65 to 88 pairs). Results from both analyses were meta-analyzed. Four major food-liking patterns were identified (Fruit and Vegetable, Distinctive Tastes, Sweet and High Carbohydrate, and Meat) accounting for 26% of the total variance. All patterns were moderately heritable (Fruit and Vegetable, h2[95% CI]: 0.36 [0.28; 0.44]; Distinctive Tastes: 0.58 [0.52; 0.64]; Sweet and High Carbohydrate: 0.52 [0.45, 0.59] and Meat: 0.44 [0.35; 0.51]), indicating genetic factors influence food liking-disliking. Overall, we identified 14 significant metabolite associations (Bonferroni p < 4.5 × 10−5) with Distinctive Tastes (8 metabolites), Sweet and High Carbohydrate (3 metabolites), and Meat (3 metabolites). Food preferences follow patterns based on similar taste and nutrient characteristics and these groupings are strongly determined by genetics. Food preferences that are strongly genetically determined (h2 ≥ 0.40), such as for meat and distinctive-tasting foods, may influence intakes more substantially, as demonstrated by the metabolomic associations identified here.

Background: The occurrence of multiple births has been recognized as a risk factor for child maltreatment. However, few population-based studies have examined the relationship between multiple births and child maltreatment. This study aimed to evaluate the degree of risk of child maltreatment among singletons and multiple births in Japan and to identify factors associated with increased risk. Methods: Using population-based data, we analyzed the database of records on child maltreatment and medical checkups for infants aged 1.5 years filed at Nishinomiya City Public Health Center between April 2007 and March 2011. To protect personal information, the data were transferred to anonymized electronic files for analysis. Results: After adjusting by logistic regression for each associated factor and gestation number, multiples themselves were not associated with the risk of child maltreatment. However, compared with singletons, multiples had a significantly higher rate of risk factors for child maltreatment, including low birth weight and neural abnormality. Moreover, compared with mothers of singleton, mothers of twins had a significantly higher rate of poor health, which is a risk factor of child maltreatment. Conclusion: Multiples were not associated with the risk of child maltreatment. However, compared with singletons, multiples and their mothers had a significantly higher rate of risk factors of child maltreatment.