Previous studies have shown that mothers of naturally conceived dizygotic (DZ) twins tend to be taller, older, and smoke more than mothers of naturally conceived monozygotic (MZ) twin and mothers of singletons. Here, we investigate whether mothers of naturally conceived DZ twins differ from mothers who conceived their DZ twins after medically assisted reproduction (MAR) in eight maternal traits related to fertility based on observational survey data. We include data from 33,648 mothers from the Netherlands Twin Register (NTR) and 1660 mothers of twins from the Norwegian Mother, Father and Child Cohort Study (MoBA). We contrast mothers of naturally conceived DZ twins with mothers of MAR DZ twins. Next, we further segment the MAR group into mothers who underwent hormonal induction of ovulation but not in vitro fertilization (IVF) and those who IVF twins, comparing them both to each other and against the mothers of naturally conceived DZ twins. Mothers of naturally conceived DZ twins smoke more often, differ in body composition, have a higher maternal age and have more offspring before the twins than mothers of MZ twins. Compared to MAR DZ twin mothers, mothers of naturally conceived DZ twins have fewer miscarriages, lower maternal age and increased height, more offspring and are more often smokers. BMI before the twin pregnancy is similar in both natural and MAR DZ twin mothers. Mothers who received hormonal induction of ovulation (OI) have a lower maternal age, fewer miscarriages, and a higher number of offspring before their twin pregnancy than twin mothers who received IVF and/or intracytoplasmic sperm injection (ICSI) treatments. Our study shows that twin mothers are a heterogenous group and the differences between twin mothers should be taken into account in epidemiological and genetic research that includes twins.

While it is known that vitamin D deficiency is associated with adverse bone outcomes, it remains unclear whether low vitamin D status may increase the risk of a wider range of health outcomes. We had the opportunity to explore the association between common genetic variants associated with both 25 hydroxyvitamin D (25OHD) and the vitamin D binding protein (DBP, encoded by the GC gene) with a comprehensive range of health disorders and laboratory tests in a large academic medical center. We used summary statistics for 25OHD and DBP to generate polygenic scores (PGS) for 66,482 participants with primarily European ancestry and 13,285 participants with primarily African ancestry from the Vanderbilt University Medical Center Biobank (BioVU). We examined the predictive properties of PGS25OHD, and two scores related to DBP concentration with respect to 1322 health-related phenotypes and 315 laboratory-measured phenotypes from electronic health records. In those with European ancestry: (a) the PGS25OHD and PGSDBP scores, and individual SNPs rs4588 and rs7041 were associated with both 25OHD concentration and 1,25 dihydroxyvitamin D concentrations; (b) higher PGS25OHD was associated with decreased concentrations of triglycerides and cholesterol, and reduced risks of vitamin D deficiency, disorders of lipid metabolism, and diabetes. In general, the findings for the African ancestry group were consistent with findings from the European ancestry analyses. Our study confirms the utility of PGS and two key variants within the GC gene (rs4588 and rs7041) to predict the risk of vitamin D deficiency in clinical settings and highlights the shared biology between vitamin D-related genetic pathways a range of health outcomes.

Social support is often considered an environmental factor affecting health, especially in aging populations. However, its genetic underpinnings suggest a more complex origin. This study investigates the heritability of social support through applying a threshold model on data of a large adult sample of twins (N = 8019) from the Netherlands Twin Register, collected between 2009 and 2011. The study employed the Duke – UNC Functional Social Support Questionnaire to assess social support quality. Our analysis revealed genetic contributions to social support, with heritability estimated at 37%, without a contribution of shared environment and no differences between men and women in heritability. The study’s results underscore the complexity of social support as a trait influenced by genetic and environmental factors, challenging the notion that it is solely an environmental construct.

In this cohort profile article we describe the lifetime major depressive disorder (MDD) database that has been established as part of the BIObanks Netherlands Internet Collaboration (BIONIC). Across the Netherlands we collected data on Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) lifetime MDD diagnosis in 132,850 Dutch individuals. Currently, N = 66,684 of these also have genomewide single nucleotide polymorphism (SNP) data. We initiated this project because the complex genetic basis of MDD requires large population-wide studies with uniform in-depth phenotyping. For standardized phenotyping we developed the LIDAS (LIfetime Depression Assessment Survey), which then was used to measure MDD in 11 Dutch cohorts. Data from these cohorts were combined with diagnostic interview depression data from 5 clinical cohorts to create a dataset of N = 29,650 lifetime MDD cases (22%) meeting DSM-5 criteria and 94,300 screened controls. In addition, genomewide genotype data from the cohorts were assembled into a genomewide association study (GWAS) dataset of N = 66,684 Dutch individuals (25.3% cases). Phenotype data include DSM-5-based MDD diagnoses, sociodemographic variables, information on lifestyle and BMI, characteristics of depressive symptoms and episodes, and psychiatric diagnosis and treatment history. We describe the establishment and harmonization of the BIONIC phenotype and GWAS datasets and provide an overview of the available information and sample characteristics. Our next step is the GWAS of lifetime MDD in the Netherlands, with future plans including fine-grained genetic analyses of depression characteristics, international collaborations and multi-omics studies.

Gestational diabetes mellitus (GDM) is a frequent complication of pregnancy. The specific mechanisms underlying GDM have not yet been fully elucidated. Contemporary research indicates a potential association between liver enzyme irregularities and an increased risk of metabolic disorders, including diabetes. The alanine aminotransferase (ALT) level is recognized as a sensitive marker of liver injury. An increase in ALT levels is hypothesized to be linked to the pathogenesis of insulin resistance and diabetes. Nonetheless, the definitive causal link between ALT levels and GDM still needs to be determined. This investigation utilized two-sample Mendelian randomization (MR) to examine the genetic causation between alanine aminotransferase (ALT) and GDM. We acquired alanine aminotransferase (ALT)-related GWAS summary data from the UK Biobank, Million Veteran Program, Rotterdam Study, and Lifeline Study. Gestational diabetes data were obtained from the FinnGen Consortium. We employed various MR analysis techniques, including inverse-variance weighted (IVW), MR Egger, weighted median, simple, and weighted weighting. In addition to MR-Egger intercepts, Cochrane’s Q test was also used to assess heterogeneity in the MR data, and the MR-PRESSO test was used to assess horizontal pleiotropy. To assess the association’s sensitivity, a leave-one-out approach was employed. The IVW results confirmed the independent risk factor for GDM development, as indicated by the ALT level (p = .011). As shown by leave-one-out analysis, horizontal pleiotrophy did not significantly skew the causative link (p > .05). Our dual-sample MR analysis provides substantiated evidence of a genetic causal relationship between alanine aminotransferase (ALT) levels and gestational diabetes.

Although the relationship between gaming addiction (GA) and attention deficit hyperactivity disorder (ADHD) is well established, the causal mechanism of this relationship remains ambiguous. We aimed to investigate whether common genetic and/or environmental factors explain the GA-ADHD relationship. We recruited 1413 South Korean adult twins (837 monozygotic [MZ], 326 same-sex dizygotic [DZ], and 250 opposite-sex DZ twins; mean age = 23.1 ± 2.8 years) who completed an online survey on GA and related traits. Correlational analysis and bivariate model-fitting analysis were conducted. Phenotypic correlation between GA and ADHD in the present sample was 0.55 (95% CI [0.51, 0.59]). Bivariate model-fitting analysis revealed that genetic variances were 69% (95% CI [64%, 73%]) and 68% (95% CI [63%, 72%]) for ADHD and GA respectively. The remaining variances (ADHD: 31%; GA: 32%) were associated with nonshared environmental variances, including measurement error. Genetic and nonshared environmental correlations between ADHD and GA were 0.68 (95% CI [0.62, 0.74]) and 0.22 (95% CI [0.13, 0.30]) respectively, which indicates that shared genes can explain 82% of the phenotypic correlation between ADHD and GA. Our study demonstrated that the ADHD-GA association was largely due to shared genetic vulnerability.

We investigated whether an observed reduction in overall childhood cancer risk (<15 years of age) in twins has been sustained, and how this extends into young adulthood. We searched for English language publications reporting childhood cancer risk in twins, obtained unpublished data directly from some authors, and updated a meta-analysis. We used the Swedish Multigeneration Register to investigate the age to which the reduced overall risk of childhood cancer (observed previously using that Swedish dataset and in this and earlier meta-analyses) persisted into the teenage/young adult years, and which specific tumors accounted for the overall risk reduction beyond childhood. Our meta-analysis of studies of aggregate childhood cancer risk in twins confirmed their approximate 15% reduction in cancer mortality and incidence. Further analysis of Swedish Multigeneration Register data for 1958 to 2002 suggested these reduced risks of cancer (particularly leukaemias and renal tumors) extended from childhood to young adult ages. Reduced risks of these and some other specific tumor types occurring across childhood/teenage/young adult years appeared to account for most of the overall risk reduction. Our results suggest a persistent reduction of overall childhood cancer risk in twins and that this extends into young adulthood. Risk reductions for several specific tumors might account for this and, although there are several potential explanations, intrauterine growth patterns of twins might be a major contributor.

Familial twinning and fertility traits were investigated in Nigerian mothers of dizygotic (DZ) twins (MoDZT; N = 972) and controls (N = 525) who responded to our person-to-person interview, which included questions on pregnancy history and family history of DZ twinning. Controls were defined as women who are not twins themselves and do not have twins in their first-degree relatives. Over 95% of the participants were Yoruba. We found that Nigerian MoDZT had an average of 4.0 (±2.6) pairs of twins among their relatives, and of these, the prevalence of DZ twins was significantly higher than that of monozygotic (MZ) twins (45.9% vs. 25.8%). Controls had an average of 0.5 (±0.4) pairs, and over 95% of the controls had no twins in their relatives. These results suggest genetic influences on DZ twinning in Nigerians. MoDZT were significantly younger in their mean age at first child, and had higher parity than controls, suggesting increased fertility in MoDZT. As compared to mothers with a single set of twins, mothers (N = 130) with multiple sets had significantly more twins among their relatives (5.4 pairs vs. 3.7 pairs) and had their first twins at a younger age (28.4 vs. 30.7 years), indicating that mothers with multiple sets of twins might have higher genetic propensity for twinning associated with earlier age at twin pregnancy. Our findings argue for genomewide association studies for DZ twinning in Nigerians, and may help to develop intervention strategies to overcome infertility/subfertility problems.

Obesity is an established risk factor for hypertension, but the mechanisms are only partially understood. We examined whether body mass index (BMI)-related DNA methylation (DNAm) variation would mediate the association of BMI with blood pressure (BP). We first conducted a genomewide DNA methylation analysis in monozygotic twin pairs to detect BMI-related DNAm variation and then evaluated the mediating effect of DNAm on the relationship between BMI and BP levels using the causal inference test (CIT) method and mediation analysis. Ontology enrichment analysis was performed for CpGs using the GREAT tool. A total of 60 twin pairs for BMI and systolic blood pressure (SBP) and 58 twin pairs for BMI and diastolic blood pressure (DBP) were included. BMI was positively associated with SBP (β = 1.86, p = .0004). The association between BMI and DNAm of 85 CpGs reached p < 1×10–4 level. Eleven BMI-related differentially methylated regions (DMRs) within LNCPRESS1, OGDHL, RNU1-44P, NPHS1, ECEL1P2, LLGL2, RNY4P15, MOGAT3, PHACTR3, and BAI2 were found. Of the 85 CpGs, 9 mapped to C10orf71-AS1, NDUFB5P1, KRT80, BAI2, ABCA2, PEX11G and FGF4 were significantly associated with SBP levels. Of the 9 CpGs, 2 within ABCA2 negatively mediated the association between BMI and SBP, with a mediating effect of −0.24 (95% CI [−0.65, −0.01]). BMI was also positively associated with DBP (β = 0.60, p = .0495). The association between BMI and DNAm of 193 CpGs reached p < 1×10−4 level. Twenty-five BMI-related DMRs within OGDHL, POU4F2, ECEL1P2, TTC6, SMPD4, EP400, TUBA1C and AGAP2 were found. Of the 193 CpGs, 33 mapped to ABCA2, ADORA2B, CTNNBIP1, KDM4B, NAA60, RSPH6A, SLC25A19 and STIL were significantly associated with DBP levels. Of the 33 CpGs, 12 within ABCA2, SLC25A19, KDM4B, PTPRN2, DNASE1, TFCP2L1, LMNB2 and C10orf71-AS1 negatively mediated the association between BMI and DBP, with a total mediation effect of −0.66 (95% CI [−1.07, −0.30]). Interestingly, BMI might also negatively mediate the association between the DNAm of most CpG mediators mentioned above and BP. The mediating effect of DNAm was also found when stratified by sex. In conclusion, DNAm variation may partially negatively mediate the association of BMI with BP. Our findings may provide new clues to further elucidate the pathogenesis of obesity to hypertension and identify new diagnostic biomarkers and therapeutic targets for hypertension.

Although there is evidence that social status has a genetic basis, it is less known whether the genetic predisposition differs between men and women as well as among different status indicators and whether there are any intercorrelations among predispositions of status indicators. We therefore investigated the genetic predisposition for different indicators of social status separately for men and women, using polygenic scores obtained from the Wisconsin Longitudinal Study. We used multivariate polygenic regression of 7 different social status indicators on a total of 24 different polygenic scores. We find that in both men and women, wages and education show more associations with polygenic scores than the other status indicators. Also, the genetic predispositions for education and wages are correlated in both men and women, whereas in men more than in women, the genetic predispositions seem to cluster into wages and education on the one hand, and status indicators of position in the hierarchy, on the other hand, with being in a management position somewhere in between. These findings are consistent with an assumption of two different forms of selection pressure associated with either cognitive skill or dominance, which holds true particularly in men. We conclude that the genetic predisposition to higher social status may have changed even though the importance of the cultural trait of social status may have been very constant. Social status may thus be an example of a social trait of constant importance, but with a changing genetic predisposition.

As with many other musical traits, the social environment is a key influence on the development of singing ability. While the familial singing environment is likely to be formative, its role relative to other environmental influences such as training is unclear. We used structural equation modeling to test relationships among demographic characteristics, familial environmental variables (early and current singing with family), vocal training, and singing ability in a large, previously documented sample of Australian twins (N = 1163). Notably, early singing with family, and to a lesser extent vocal training, predicted singing ability, whereas current singing with family did not. Early familial singing also mediated the relationship between sex and singing ability, with men who sang less with family during childhood showing poorer ability. Bivariate twin models between early familial singing and singing ability showed the phenotypic correlation was largely explained by shared environmental influences. This raises the possibility of a sensitive period for singing ability, with sociocultural expectations around singing potentially differentiating the developmental trajectories of this skill for men and women.

TwinsMX registry is a national research initiative in Mexico that aims to understand the complex interplay between genetics and environment in shaping physical and mental health traits among the country’s population. With a multidisciplinary approach, TwinsMX aims to advance our knowledge of the genetic and environmental mechanisms underlying ethnic variations in complex traits and diseases, including behavioral, psychometric, anthropometric, metabolic, cardiovascular and mental disorders. With information gathered from over 2800 twins, this article updates the prevalence of several complex traits; and describes the advances and novel ideas we have implemented such as magnetic resonance imaging. The future expansion of the TwinsMX registry will enhance our comprehension of the intricate interplay between genetics and environment in shaping health and disease in the Mexican population. Overall, this report describes the progress in the building of a solid database that will allow the study of complex traits in the Mexican population, valuable not only for our consortium, but also for the worldwide scientific community, by providing new insights of understudied genetically admixed populations.

This study examines the temporal and geographical evolution of polygenic scores (PGSs) across cognitive measures (Educational Attainment [EA], Intelligence Quotient [IQ]), Socioeconomic Status (SES), and psychiatric conditions (Autism Spectrum Disorder [ASD], schizophrenia [SCZ]) in various populations. Our findings indicate positive directional selection for EA, IQ, and SES traits over the past 12,000 years. Schizophrenia and autism, while similar, showed different temporal patterns, aligning with theories suggesting they are psychological opposites. We observed a decline in PGS for neuroticism and depression, likely due to their genetic correlations and pleiotropic effects on intelligence. Significant PGS shifts from the Upper Paleolithic to the Neolithic periods suggest lifestyle and cognitive demand changes, particularly during the Neolithic Revolution. The study supports a mild hypothesis of Gregory Clark’s model, showing a noticeable rise in genetic propensities for intelligence, academic achievement and professional status across Europe from the Middle Ages to the present. While latitude strongly influenced height, its impact on schizophrenia and autism was smaller and varied. Contrary to the cold winters theory, the study found no significant correlation between latitude and intelligence.

Voice production can be influenced by interindividual variations related to genetic, physiological, behavioral, and several environmental factors. Here we examined the effect of zygosity on speaking fundamental frequency (F0) statistical descriptors. Our aims were: (1) to determine whether the genetic similarity between monozygotic (MZ) and dizygotic (DZ) twins affects F0 characteristics, and (2) to quantify the contribution of genetic factors to these characteristics. The study involved 79 same-sex twin pairs of Brazilian Portuguese speakers, comprising 65 MZ and 14 DZ twins, aged 18 to 66 years (31.7 ± 11.6 years), with 21 male and 58 female pairs. Participants were recorded while uttering a greeting phrase and the Brazilian Portuguese version of the ‘Happy Birthday to You’ song. Speech segments were analyzed using Praat free software, and F0 measures were automatically extracted in both Hertz and semitone scales. Statistical descriptors, including centrality, dispersion, and extreme values of F0 were examined, and the ACE model (i.e., total genetic effects, A; shared environmental influences, C; and nonshared environmental influences, E) was employed to estimate the additive effect;ts of monozygosity. As anticipated, we observed a zygosity effect on several F0 parameters, with more similarity between MZ twins compared to DZ twins. We discuss the genetic influences on F0 parameters and the absence of a monozygosity effect in two of them. Additionally, we briefly address potential biases associated with the selected measurement scale for statistical modeling. Finally, we explore the influence of genetic factors on F0 patterns, as well as environmental, life history and linguistic factors, particularly concerning F0 variation in speech.

Researchers have proposed that culture significantly influences perceived stress (PS). To date, however, twin studies on PS have been conducted mostly in western, individualistic cultures, which demonstrate that PS due to controllable (personal) life events is more heritable than PS due to uncontrollable (network) life events. This study aimed to investigate genetic and environmental influences on PS in South Korean twins. South Korea practices a dominant collectivist culture. In total, 1372 twin individuals (mean age = 22.4 ± 2.5 years) completed an online survey on PS, which consisted of the scales, Friendship, Academic Stress, Future Career, Family Conflicts, and Family Financial Difficulties (FFD). Friendship, Academic Stress, and Future Career can be considered PS due to personal life events, and Family Conflict and FFD, PS due to network life events. The general sex-limitation model-fitting analysis revealed the absence of qualitative or quantitative sex differences in genetic and environmental influences. Specifically, additive genetic influences were predominant for Friendship (63%), Academic Stress (67%), and Future Career (57%) for both sexes, with the remaining variance attributable to nonshared environmental influences. In contrast, shared environmental influences were largest for Family Conflict (47% for both genders) and FFD (64% for males, 63% for females) with no significant genetic effects. Despite known cultural differences in the means and variances of PS, South Korean twins exhibited significant genetic effects in PS due to personal life events and large shared environmental effects in PS due to network life events, which is similar to western samples.

Women with twin pregnancies experience greater sleep disturbance compared to women with singleton pregnancies. The aims of this study were to explore the sleep quality in women with twin pregnancies and to compare their sleep dimensions with coetaneous single pregnancies. This was an observational study in which women were enrolled at the end of pregnancy in the Obstetric Service of Hospital La Paz (Spain). The women were classified as single (n = 143) or twin pregnancy (n = 62). Pregnant women responded to the Pittsburgh Sleep Quality Index to evaluate sleep quality, latency, duration, efficiency, perturbance, use of medication, and daytime dysfunction. The higher the index, the greater the alteration of sleep quality. Without statistical differences, a poor sleep quality was higher in women with single (66.7%) than women with twin pregnancies (22.8%). The good sleeper slept 6.8 h/day in single pregnancy and 7.3 h/day in twin pregnancy. The sleep perturbation and dysfunctionality were higher in women with twin than single pregnancies. The use of medication to sleep was significantly lower in women with twin than single pregnancies. In women with twin pregnancy, the body weight gain during first trimester had a positive correlation with worse sleep quality and sleep perturbations. Twin pregnancy needed more than 7 h/day to have a high sleep quality, showing greater sleep perturbations and daytime dysfunction than single pregnancies. The control of gestational body weight can improve the sleep quality, disturbances, and duration in twin gestations. Sleep screening during pregnancy would be necessary to handle sleep issues and increase benefits in twin gestational outcomes.

There are sex-dependent differences in hematological and biochemical variables in adulthood attributed to the predominant effects of testosterone in males and estrogen in females. The Twin Testosterone Transfer (TTT) hypothesis proposes that opposite-sex females may develop male-typical traits due to exposure to relatively higher levels of prenatal testosterone than same-sex females. Additionally, prenatal testosterone exposure has been suggested as a correlate of current circulating testosterone levels. Consequently, opposite-sex females might exhibit male-typical patterns in their hematological and biochemical variables. Despite this hypothesis, routine laboratory investigations assign the same reference range to all females. Our cross-sectional study, conducted in Tamale from January to September 2022, included 40 twins, comprising 10 opposite-sex (OS) males (25%), 10 OS females (25%), and 20 same-sex (SS) females (50%), all aged between 18 and 27 years. Fasting venous blood samples were collected and analyzed using automated hematology and biochemistry laboratory analyzers. Results indicated that levels of hemoglobin, serum creatinine, gamma-glutamyl transferase, total protein, globulins, and total testosterone were significantly higher in OS males than OS females. Conversely, total cholesterol and low-density lipoprotein cholesterol were significantly higher in OS females than OS males. Unexpectedly, levels of low-density lipoprotein cholesterol and total testosterone were significantly higher in SS females than OS females. Contrary to expectations, opposite-sex females did not exhibit male-typical patterns in their hematological and biochemical variables. This suggests that the TTT effect may not occur or may not be strong enough to markedly affect hematological and biochemical variables in OS females.

The current study was motivated by an interest in deepening understanding of Brazilian twin research, which is underrepresented internationally, in an effort to rectify this situation. Our aim was threefold: (1) to carry out a comprehensive investigation of Brazilian research on twins according to the area of knowledge; (2) to evaluate the representation of research in the field of psychology in comparison with other areas; (3) to evaluate characteristics of the research that may have contributed to its exclusion from the comprehensive meta-analysis of 50 years of twin research. A scoping review was performed according to PRISMA guidelines. Titles and abstracts were searched up to 2022 in six databases: CAPES, BDLTD, PePSIC, PubMed, Google Scholar, and SciELO, using selected keywords both in Portuguese and in English (e.g., ‘twins’ and ‘Brazil’; ‘twinning’ and ‘Brazil’; ‘gemelaridade’ [twinning], and ‘gêmeos’ [twins]). Three hundred and forty publications were included in the review. Approximately half (53.8‰) used the classic twin design to investigate the heritability of several traits, and the other half (46.2%) used other research designs. The scoping review showed that the number of publications doubled approximately every 10 years. Most publications were from the health area, with medicine accounting for approximately half of the studies, followed by psychology, odontology, and biology. We found that the interest in studying twins among Brazilian scientists is increasing over the years and there are reasons to be enthusiastic about the potential impact of this trend in the global scenario.

Immunity activation and inflammation are the main characteristics of rheumatoid arthritis and clonal hematopoiesis. However, it remains unclear whether rheumatoid arthritis increase the risk of clonal hematopoiesis. Here, a Mendelian randomization (MR) analysis was conduct to explore the causal effects of rheumatoid arthritis on clonal hematopoiesis. Summary statistics data of rheumatoid arthritis (13,838 cases and 33,742 controls) and clonal hematopoiesis (10,203 cases and 173,918 controls) derived from a genomewide association study were selected to analyze. We selected inverse-variance weighted, MR-Egger, weighted median, simple mode, and weighted mode to evaluate the causal effect of rheumatoid arthritis on clonal hematopoiesis. The two-sample MR analysis suggested a strong causal relationship between rheumatoid arthritis and clonal hematopoiesis by inverse-variance weighted (OR = 1.002311673, 95% CI [1.000110757, 1.004517433], p = .039706) and weighted median (OR = 1.002311673, 95% CI [1.000110757, 1.004517433], p = .039518447) methods. No significant pleiotropy or heterogeneity was found in the sensitivity analysis. These results supported a potentially causal relationship between rheumatoid arthritis and clonal hematopoiesis, and the exposure of rheumatoid arthritis increased the risks of clonal hematopoiesis. Our findings highlight the importance of how chronic inflammation and immune activation induced rheumatoid arthritis enhances the risks of clonal hematopoiesis, and that early intervention with rheumatoid arthritis patients might reduce the clonal hematopoiesis risks in rheumatoid arthritis patients. Moreover, our study provides clues for prediction of risk factors and potential mechanisms of clonal hematopoiesis.

This present study investigated the parental characteristics of multiple births using national birth data in Japan. This study included birth data from Vital Statistics: Occupational and Industrial Aspects every five fiscal years from 1995 to 2020. The multiple birth rates were defined as the number of live-birth deliveries with multiple fetuses (e.g., twins, triplets) per total live-birth deliveries. Parental ages, nationalities, occupations and household occupation (occupation of the top earner of the household) were considered as parental characteristics. The multiple birth rates were calculated based on parental characteristics for each year, and a log-binomial regression model was used to assess the association between parental characteristics and multiple births. The multiple birth rate for Japanese mothers consistently exceeded that for non-Japanese mothers over the years, and the rate increased progressively from manual workers to lower non-manual workers and then to upper non-manual workers for both maternal and paternal occupations. The regression results indicated that the risk ratio (RR) for multiple births among non-Japanese mothers was significantly lower than that among Japanese mothers. Moreover, concerning household occupation, the RRs of self-employed individuals, full-time employees at smaller companies, others, and the unemployed were significantly lower than those of full-time employees at larger companies. Furthermore, the RRs of lower non-manual and manual workers were significantly lower than those of upper non-manual workers in maternal and paternal occupations. The results suggested an association between multiple births and parental socioeconomic status in Japan.